《弹性力学》第五章平面问题的复变函数法严选课资课件

1第五章 平面问题的复变函数法1第五章第五章 平面问题的复变函数法平面问题的复变函数法 直角坐标及极坐标求解平面问题，所涉直角坐标及极坐标求解平面问题，所涉及的物体边界是直线或圆弧形。对于其他一及的物体边界是直线或圆弧形。对于其他一些边界，例如椭圆形、双曲形、非同心圆等些边界，例如椭圆形、双曲形、非同心圆等就要用不同的曲线坐标。应用复变函数可使就要用不同的曲线坐标。应用复变函数可使该类问题得以简化。本章只限于介绍复变函该类问题得以简化。本章只限于介绍复变函数方法在弹性力学中的简单应用。数方法在弹性力学中的简单应用。2平面问题的复变函数法第五章 平面问题的复变函数法 直角 5-4 5-4 多连通域内应力与位移的单值条件多连通域内应力与位移的单值条件5-3 5-3 边界条件的复变函数表示边界条件的复变函数表示5-2 5-2 应力和位移的复变函数表示应力和位移的复变函数表示5-1 5-1 应力函数的复变函数表示应力函数的复变函数表示5-6 5-6 含孔口的无限大板问题含孔口的无限大板问题5-5 5-5 无限大多连体的情形无限大多连体的情形第五章第五章 平面问题的复变函数法平面问题的复变函数法3 5-4 多连通域内应力与位移的单值条件5 5-1 5-1 应力函数的复变函数表示应力函数的复变函数表示 在第二章中已经证明，在平面问题里，如果体力是常量，就一定存在一个应力函数，它是位置坐标的重调和函数，即现在，引入复变数z=xiy和zxiy以代替实变数x 和y。注意4 5-1 应力函数的复变函数表示 在第二章中已 可以得到变换式进而5 可以得到变换式进而平面问题的复变函数法5令于是可将方程式变换成为由6令于是可将方程式变换成为由平面问题的复变函数法6 可知，P是调和函数可由解析函数的实部得到。设f(z)为解析函数，可令由令得则7 可知，P是调和函数可由解析函数的实部得到。设f(z)为 将上式对 积分，得到再对z积分，得到令即则8 将上式对 积分，得到再对z积分，得到令即则平面问题的 注意上式左边的重调和函数是实函数，可见该式右边的四项一定是两两共轭，前两项已经是共轭的，后两项也应是共轭的：令即得有名的古萨公式也可以写成9 注意上式左边的重调和函数是实函数，可见该式右边的四项 于是可见，在常量体力的平面问题中，应力函数总可以用复变数z的两个解析函(z)和(z)来表示，称为K-M 函数。而求解各个具体的平面问题，可归结为适当地选择这两个解析函数，并根据边界条件决定其中的任意常数。10 于是可见，在常量体力的平面问题中，应力函数总 5-2 5-2 应力和位移的复变函数表示应力和位移的复变函数表示根据应力分量和应力函数的关系一一 应力分量的复变函数表示应力分量的复变函数表示11 5-2 应力和位移的复变函数表示根据应力分量和应力 可得到应力分量的复变函数表示由可得而由12 可得到应力分量的复变函数表示由可得而由平面问题的复变函 可得或13 可得或平面问题的复变函数法13 只要已知(z)及(z)，就可以把上述公式右边的虚部和实部分开，由虚部得出xy，由实部得出y-x。和就是应力分量的复变函数表示。当然也可以建立公式，把x、y、xy三者分开用(z)和(z)来表示，但那些公式将比较冗长，用起来很不方便。14 只要已知(z)及(z)，就可以把 二二 位移分量的复变函数表示位移分量的复变函数表示 假定为平面应力问题。由几何方程及物理方程可得15 二 位移分量的复变函数表示 假定为平面应力问题 由于并注意到同理可得16 由于并注意到同理可得平面问题的复变函数法16将上两式分别对x及y积分，得其中的f1及f2为任意函数。将上式代入式17将上两式分别对x及y积分，得其中的f1及f2为任意函数。将上 由于18 由于平面问题的复变函数法18 从而得到于是得到刚体位移 f1(y)u0y，f2(x)v 0 x故有19 从而得到于是得到刚体位移 故有平面问题的复变函数法19 若不计刚体位移，则有由式得到20 若不计刚体位移，则有由式得到平面问题的复变函数法20 这就是位移分量的复变函数表示。若已知(z)及(z)，就可以将该式右边的实部和虚部分开，从而得出u和v。将结果回代,并两边除以 得 上述公式是针对平面应力情况导出的。对于平面应变情况，须将式中的E改换为 ，改换为 。21 这就是位移分量的复变函数表示。若已知(z)及(z5-3 5-3 边界条件的复变函数表示边界条件的复变函数表示 为了求得边界上各结点处的值，须要应用应力边界条件，即：而代入上式，即得：225-3 边界条件的复变函数表示 为了求得边 由图可见，l=cos(N,x)=dy/ds,m=cos(N,y)=-dx/ds,于是，前式可改写为：由此得：23 由图可见，l=cos(N,x)=dy/ds,m=c 设A是边界上的固定点,B为任意一点，则从A到B边界上的合力，可用上式从A点到B点对s积分得到：将式24 设A是边界上的固定点,B为任意一点，则从A到 代入，整理得：把应力函数加上一个复常数，并不影响应力。因此，可把应力函数A处的值设为零，于是对于边界上的有或这就是应力边界条件。25 代入，整理得：把应力函数加上一个复常数，并不 对于位移边界条件将其代入下式即得平面应力情况下位移边界条件的复变函数表示 对于平面应变，须将式中的E改换为 ，改换为 。26 对于位移边界条件将其代入下式即得平面应力情况 5-4 5-4 多连通域内应力与位移的单值条件多连通域内应力与位移的单值条件 应力确定后，应力函数仍可差一个任意的线性函数，这时K-M函数并未完全确定。对于单连通区域，可以通过选取适当坐标系等办法，使得K-M函数完全确定；但对于多连通区域仍不能完全确定。本节讨论K-M函数在多连通区域内满足单值的条件。设有多连通区域，有一内边界C，设在边界C上的外力矢量已给定。通常的多值函数是对数函数，我们设27 5-4 多连通域内应力与位移 DC这里zk为内部边界内的任意一点，f和f为单值的解析函数（全纯函数），而Ak，Bk为常数：28 DC这里zk为内部边界内的任意一点，f和f 前面的函数的导数是单值的，但他们本身是多值的，当z绕周边一周时，函数值ln(zk)产生一个增量2i,于是(z)和(z)的增量分别是2i Ak和2iBk,这时应力主矢量按照公式左边将得到应力主矢量(沿整个边界），右边得到一增量：29 前面的函数的导数是单值的，但他们本身是多值的，这时位移按照公式也将得到增量，根据单值性这个增量应为零：结合可得到30 这时位移按照公式也将得到增量，根据单值性这个增 于是当有m个内边界时，取31 于是当有m个内边界时，取平面问题的复变函数法31 5-5 5-5 无限大多连体的情形无限大多连体的情形 当多连体的外边界趋于无限远时，该多连体成为无限大的多连体，除上述条件外，还需考虑无限远的极限情况。以坐标原点为圆心，作充分大的圆周sR，将所有的内边界包围在其内，对于sR之外，弹性体之内的任意一点，可得到 在在s sR R之外的解析函数之外的解析函数32 5-5 无限大多连体的情形 当多连体的外边 于是可写为其中Px,Py为m个边界上沿x,y方向的面力之和。33 于是可写为其中Px,Py为m个边界上沿x,y方向的面 于是 由于在无穷远处的应力分量应该是有限的，级数中n2的系数应为零。将多连通区域内的全纯函数 和 展开为罗郎级数：34 于是 由于在无穷远处的应力分量应该是有限的，级数 同样从中，由于在无穷远处的应力分量应该是有限的，故有其中略去了和应力无关的常数项。35 同样从中，由于在无穷远处的应力分量应该是有限的 于是其中与应力计算无关，可取为零，而36 于是其中与应力计算无关，可取为零，而平面问题 这时当z时，可得同样当z时，由可得从中可求得相应的系数，并可以看到在无限远处，应力的分布是均匀的。37 这时当z时，可得同样当z时，由可得从中 系数则38 系数则平面问题的复变函数法385-6 5-6 含孔口的无限大板问题含孔口的无限大板问题 以坐标原点为圆心，作充分大的圆周sR，将所有的内边界包围在其内，对于sR之外，弹性体之内的任意一点，可得到395-6 含孔口的无限大板问题 以坐标原点为圆40平面问题的复变函数法40 改写为其中41 改写为其中平面问题的复变函数法41 对于孔边上的点42 对于孔边上的点平面问题的复变函数法42 将上列各式代入就得到极坐标下圆周边界上的级数形式的应力边界条件。设周边上的外力为已知，并将其展开为傅氏级数43 将上列各式代入就得到极坐标下圆周边界上的级数形式的应 比较两边eik和e-ik的系数，可得44 比较两边eik和e-ik的系数，可得平面问题的复 由无限远处的应力条件，可得45 由无限远处的应力条件，可得45由位移的单值条件有及可求得再由46由位移的单值条件有及可求得再由平面问题的复变函数法46 可求得至此，全部系数均已求出。例例 设孔周边为均匀压力p，无限远处的应力为零。47 可求得至此，全部系数均已求出。例 设孔周边为均匀压 则有于是可求得48 则有于是可求得平面问题的复变函数法48 最后得到根据上述方法，圆孔口无限大板的一般问题都可以得到解决。49 最后得到根据上述方法，圆孔口无限大板的一般问题都可以练习5.1 试考察下列复变函数所解决的问题(1)(2)解:基本公式为(1)将分别代入(a)、(b)式50平面问题的复变函数法练习5.1 试考察下列复变函数所解决的得联立求解以上两式,得 所给的函数可以解决矩形薄板在x方向受均布拉力q的问题.如图5.1(a)所示(2)将代入(a),(b)两式,得xyqq图5.1(a)51平面问题的复变函数法得联立求解以上两式,得 所给的函数联立求解以上两式,得 所给的函数可以解决矩形薄板受纯剪切问题.如图5.1(b)示.qqxy图5.1(b)练习5.2 如图所示.试证矩形截面梁的纯弯曲问题可用如下的复变函数求解.其中I为梁截面的惯矩,M为作用的弯矩.Myxzy解:基本公式为52平面问题的复变函数法联立求解以上两式,得 所给的函数可将代入(1)、(2)式由(1)式得即53平面问题的复变函数法将代入(1)、(2)式由(1)式得即53或由(2)式得即将(4)、(5)式联立求得54平面问题的复变函数法或由(2)式得即将(4)、(5)式联立求验证边界条件(3)在侧面:所以由得55平面问题的复变函数法验证边界条件(3)在侧面:所以由得55由得故即(3)式恒成立.由解答 所表示的是一个纯弯时,梁横截面上的应力状态.56平面问题的复变函数法由得故即(3)式恒成立.由解答 练习5.3 试导出用复变函数 及 表示极坐标中应力分量的公式解:因为在平面问题中所以又因为在平面问题中,有57平面问题的复变函数法练习5.3 试导出用复变函数 则58平面问题的复变函数法则58因为所以练习5.4 试用公式由 导出半平面体在边界上受集中力作用时的应力分量公式.59平面问题的复变函数法因为所以练习5.4 试用公式由 ryroP解：由得因为60平面问题的复变函数法ryroP解：由得因为60而所以61平面问题的复变函数法而所以61即由（1）、（2）、（3）式得62平面问题的复变函数法即由（1）、（2）、（3）式得6263结 束 平面问题的复变函数法63