平面问题的直角坐标解答

《平面问题的直角坐标解答》由会员分享，可在线阅读，更多相关《平面问题的直角坐标解答（87页珍藏版）》请在装配图网上搜索。

1、要点要点 用逆解法、用逆解法、半逆解法半逆解法求解平面弹性求解平面弹性力学问题。力学问题。3-1 3-1 多项式解答多项式解答3-2 3-2 位移分量的求出位移分量的求出3-3 3-3 简支梁受均布载荷简支梁受均布载荷3-4 3-4 楔形体受重力和液体压力楔形体受重力和液体压力3-5 3-5 级数式解答级数式解答3-6 3-6 简支梁受任意横向载荷简支梁受任意横向载荷主主 要要 内内 容容3-1 3-1 多项式解答多项式解答适用性：适用性：由一些直线边界构成的弹性体。由一些直线边界构成的弹性体。目的：目的：考察一些简单多项式函数作为应力函数考察一些简单多项式函数作为应力函数(x,y)，能解决什

2、么样的，能解决什么样的力学问题。力学问题。逆解法逆解法cbyaxyx),(其中：其中：a、b、c 为待定系数。为待定系数。检验检验(x,y)是否满足双调和方程：是否满足双调和方程：0244224444yyxx显然显然(x,y)满足双调和方程，因而可作为应力函数。满足双调和方程，因而可作为应力函数。（1）1.一次多项式一次多项式（2）（3）对应的应力分量：对应的应力分量：02yxxyXxyx22XxXx 0YyYy 0Yyxy22若体力：若体力：X=Y=0，则有：，则有：0 xyyx结论结论1：（1）（2）一次多项式对应于一次多项式对应于无体力和无应力状态；无体力和无应力状态；在该函数在该函数(

3、x,y)上加上或减去一个一次多项式，对应力无影响。上加上或减去一个一次多项式，对应力无影响。2.二次多项式二次多项式（1）22cybxyax其中：其中：a、b、c 为待定系数。为待定系数。(假定：假定：X=Y =0;a 0,b 0,c 0)检验检验(x,y)是否满足双调和方程，显然有是否满足双调和方程，显然有（2）0,0,02244444yxyx04(可作为应力函数可作为应力函数)（3）由式（由式（2-26）计算应力分量：）计算应力分量：byxxy2cyx222axy222xy2c2c2a2abxy结论结论2：二次多项式对应于二次多项式对应于均匀应力分布。均匀应力分布。xy0202yxy0试求

4、图示板的应力函数。试求图示板的应力函数。例：例：xy00 xyyx0),(202),(yyx3.三次多项式三次多项式（1）3223dycxyybxax其中其中:a、b、c、d 为待定系数。为待定系数。检验检验(x,y)是否满足双调和方程，显然有是否满足双调和方程，显然有（2）0,0,02244444yxyx04(可作为应力函数可作为应力函数)(假定：假定：X=Y =0)（3）由式（由式（2-26）计算应力分量：）计算应力分量：cybxyxxy222dycxyx6222axbyxy6222结论结论3：三次多项式对应于三次多项式对应于线性应力分布。线性应力分布。讨论：讨论：,3dy取)0(YX可算

5、得：可算得：0 xydyx60yxy12h2hll图示梁对应的边界条件：图示梁对应的边界条件：:2hy0,0 xyy:lx0,6xyxdydh3mindh3maxMM3dy可见：可见：对应于矩形截面梁的对应于矩形截面梁的纯弯曲问题纯弯曲问题应力分布。应力分布。常数常数 d 与弯矩与弯矩 M 的关系：的关系：220hhxdy(1)由梁端部的边界条件：由梁端部的边界条件：0622hhdydy(2)Mdyyhhx222226hhMdydyMhd32)2(3hMd 或yIMxyhMx312yhMx)12/(3可见：可见：此结果与材力中结果相同，此结果与材力中结果相同，说明材力中纯弯曲梁的应力结果是正确

6、的。说明材力中纯弯曲梁的应力结果是正确的。0 xydyx60yxy12h2hllMMyIMxdh3mindh3max说明：说明：(1)组成梁端力偶组成梁端力偶 M 的面力的面力须线性须线性分布分布，且中心处为零，结果才，且中心处为零，结果才是是精确的精确的。(2)若按其它形式分布，如：若按其它形式分布，如：则此结果不精确，有误差；则此结果不精确，有误差；但按圣维南原理，仅在两端误差较但按圣维南原理，仅在两端误差较大，离端部较远处误差较小。大，离端部较远处误差较小。(3)当当 l 远大于远大于 h 时，误差较小；反之误差较大。时，误差较小；反之误差较大。4.四次多项式四次多项式（1）432234

7、eydxyycxybxax检验检验(x,y)是否满足双调和方程是否满足双调和方程（2）cyx8244ax2444ey2444代入：代入：04得得033eca024824eca432234eydxyycxybxax可见，对于函数：可见，对于函数：其待定系数，须满足下述关系才能作为应函数：其待定系数，须满足下述关系才能作为应函数：033eca（3）应力分量：应力分量：yxxy222343dycxybx22yx22xy221262eydxycx221262axbxycy 应力分量为应力分量为 x、y 的二次函数。的二次函数。（4）特例：特例：44eyax 212eyx0 xy212axy（须满足：（

8、须满足：a+e=0）总结：总结：（多项式应力函数（多项式应力函数 的性质）的性质）（1）多项式次数多项式次数 n 4 时，则系数可以任意选取，总可满足时，则系数可以任意选取，总可满足 。04多项式次数多项式次数 n 4 时，则系数时，则系数须满足须满足一定条件，才能满足一定条件，才能满足 。04多项式次数多项式次数 n 越高，则系数间需满足的条件越多。越高，则系数间需满足的条件越多。（2）一次多项式，对应于一次多项式，对应于无体力和无应力状态；无体力和无应力状态；任意应力函数任意应力函数(x,y)上加上加上或减去一个上或减去一个一次多项式一次多项式，对应力无影响。，对应力无影响。二次多项式二次

9、多项式，对应，对应均匀应力均匀应力状态，即全部应力为常量；状态，即全部应力为常量；三次多项式三次多项式，对应于对应于线性分布应力线性分布应力。（3）（4）用多项式构造应力函数用多项式构造应力函数(x,y)的方法的方法 逆解法（只能解决简单逆解法（只能解决简单直直线应力边界线应力边界问题）。问题）。按应力求解平面问题，其基本未知量为：按应力求解平面问题，其基本未知量为：，本节说明，本节说明如何由如何由 求出形变分量、位移分量？求出形变分量、位移分量？xyyx,xyyx,问题：问题：3-2 3-2 位移分量的求出位移分量的求出以纯弯曲梁为例，说明如何由以纯弯曲梁为例，说明如何由 求出形变分量、位移

10、分量？求出形变分量、位移分量？xyyx,xyl1hMM1.形变分量与位移分量形变分量与位移分量由前节可知，其应力分量为：由前节可知，其应力分量为：12/3hMyyIMx0 xy0y平面应力情况下的物理方程：平面应力情况下的物理方程：（1）形变分量）形变分量)(1xyyE)(1yxxEGxyxy(a)将式（将式（a）代入得：）代入得：IMyEyIMyEx10 xy(b)（2）位移分量）位移分量将式（将式（b）代入几何方程得：）代入几何方程得：0 xvyuxyIMyExux1IMyEyvy(c)（2）位移分量）位移分量0 xvyuxyIMyExux1IMyEyvy(c)将式（将式（c）前两式积分，

11、得：）前两式积分，得：)(222xfyEIMv)(1yfxyEIMu(d)将式将式(d)代入代入(c)中第三式，得：中第三式，得：)(),(21xfyf式中：式中：为待定函数。为待定函数。)()(12yfxfxEIM整理得：整理得：0)()(21xfyfxEIM（仅为（仅为 x 的函数）的函数）（仅为（仅为 y 的函数）的函数）要使上式成立，须有要使上式成立，须有)(2xfxEIM)(1yf(e)式中：式中：为常数。为常数。积分上式，得积分上式，得01)(uyyf022)(vxxEIMxf将上式代入式（将上式代入式（d），得），得0uyxyEIMu02222vxxEIMyEIMv(f)（1）(

12、f)讨论：讨论：式中：式中：u0、v0、由位移边界条件确定。由位移边界条件确定。常数00 xEIMyuxx当当 x=x0=常数常数xEIMyu（2）位移分量）位移分量0uyxyEIMu02222vxxEIMyEIMvxyl1hMM u 关于铅垂方向的变化率，即铅垂方向线段的转角。关于铅垂方向的变化率，即铅垂方向线段的转角。常数00 xEIMyuxxyu0|xx说明：说明：同一截面上的各铅垂同一截面上的各铅垂线段转角相同线段转角相同。横截面保持平面横截面保持平面 材力中材力中“平面保持平面平面保持平面”的假设成的假设成立立。（2）常数EIMxv22102222vxxEIMyEIMv将下式中的第二

13、式对将下式中的第二式对 x 求二阶导数：求二阶导数：0uyxyEIMu说明：说明：在微小位移下，梁纵向纤维的曲在微小位移下，梁纵向纤维的曲率相同。即率相同。即EIMxv221 材料力学中挠曲线微分方程材料力学中挠曲线微分方程2.位移边界条件的利用位移边界条件的利用（1）两端简支）两端简支02222vxxEIMyEIMv0uyxyEIMu(f)其边界条件：其边界条件：000yxu000yxv将其代入将其代入(f)式，有式，有0202vlEIMl00u00vEIMl2将其代回将其代回(f)式，有式，有ylxEIMu)2(22)(2yEIMxxlEIMv(3-3)梁的挠曲线方程：梁的挠曲线方程：xx

14、lEIMvy)(20 与材力中结果相同与材力中结果相同00ylxv（2）悬臂梁）悬臂梁02222vxxEIMyEIMv0uyxyEIMu(f)边界条件边界条件0lxv0lxu22hyhh/2h/2由式（由式（f）可知，此边界条件无法满足。）可知，此边界条件无法满足。边界条件改写为：边界条件改写为：0,000ylxylxvu（中点不动）（中点不动）00ylxxv（轴线在端部不转动）（轴线在端部不转动）代入式（代入式（f），有），有00u0202vllEIM0lEIM可求得：可求得：00uEIMlv220EIMlyxlEIMu)(222)(2yEIMxlEIMvyxlEIMu)(222)(2yEI

15、MxlEIMv(3-4)h/2h/2挠曲线方程：挠曲线方程：20)(2|xlEIMvy与材料力学中结果相同与材料力学中结果相同说明：说明：（1）求位移的过程：求位移的过程：（a）将应力分量代入物理方程）将应力分量代入物理方程)(1xyyE)(1yxxEGxyxy（b）再将应变分量代入几何方程）再将应变分量代入几何方程xvyuxyxuxyvy（c）再利用位移边界条件，确定常数。）再利用位移边界条件，确定常数。（2）若为平面应变问题，则将材料常数若为平面应变问题，则将材料常数E、作相应替换。作相应替换。（3）若取固定端边界条件为：若取固定端边界条件为：h/2h/20,000ylxylxvu（中点不

16、动）（中点不动）00ylxyu（中点处竖向线段转角为零）（中点处竖向线段转角为零）00u得到：得到：0202vlEIMl0EIMl02222vxxEIMyEIMv0uyxyEIMu求得：求得：00uEIMlv220EIMl此结果与前面情形相同。此结果与前面情形相同。（为什么？为什么？）（1）024422444yyxx(2-27)（2）xyyx,然后将然后将 代入式（代入式（2-26）求出应力分量：）求出应力分量：),(yx先由方程（先由方程（2-27）求出应力函数：）求出应力函数：),(yxYyxy22Xxyx22yxxy2(2-26)（3）再让再让 满足应力边界条件和位移单值条件（多连体问题

17、）。满足应力边界条件和位移单值条件（多连体问题）。xyyx,04按应力求解平面问题的基本步骤：按应力求解平面问题的基本步骤：按应力求解平面问题的方法按应力求解平面问题的方法：逆逆解解法法（1）根据问题的条件根据问题的条件（几何形状、受力特点、边界条件等），（几何形状、受力特点、边界条件等），假设各种满足相容方程（假设各种满足相容方程（2-27）的）的(x,y)的形式；的形式；（2）然后利用应力分量计算式（然后利用应力分量计算式（2-26），求出），求出 （具有待（具有待定系数）；定系数）；xyyx,（3）再利用应力边界条件式（再利用应力边界条件式（2-18），来考察这些应力函数），来考察这些应

18、力函数(x,y)对对应什么样的边界面力问题，从而得知所设应力函数应什么样的边界面力问题，从而得知所设应力函数(x,y)可以求可以求解什么问题。解什么问题。（1）根据问题的条件根据问题的条件（几何形状、受力特点、边界条件等），（几何形状、受力特点、边界条件等），假设部分应力分量假设部分应力分量 的某种函数形式的某种函数形式；xyyx,（2）根据根据 与应力函数与应力函数(x,y)的关系及的关系及 ，求，求出出(x,y)的形式；的形式；xyyx,04（3）最后利用式（最后利用式（2-26）计算出）计算出 并让其满足边界条件和并让其满足边界条件和位移单值条件。位移单值条件。xyyx,半逆解法的数学基

19、础：半逆解法的数学基础：数理方程中分离变量法数理方程中分离变量法。半逆解法半逆解法位移分量求解位移分量求解：（1）将已求得的应力分量将已求得的应力分量（2）（3）xyyx,代入物理方程，求得应变分量代入物理方程，求得应变分量xyyx,将应变分量将应变分量xyyx,代入几何方程，并积分求得位移分量代入几何方程，并积分求得位移分量表达式；表达式；由位移边界条件确定表达式中常数，得最终结果。由位移边界条件确定表达式中常数，得最终结果。3-3 3-3 简支梁受均布载荷简支梁受均布载荷要点要点 用用半逆解法半逆解法求解梁、长板类平面问题。求解梁、长板类平面问题。xyllqlql1yzh/2h/2q1.应

20、力函数的确定应力函数的确定(1)分析：分析：y 主要由弯矩引起；主要由弯矩引起；x 主要由剪力引起；主要由剪力引起；xy由由 q 引起（挤压应力）。引起（挤压应力）。又又 q=常数，图示坐标系和几何对称，常数，图示坐标系和几何对称，不随不随 x 变化。变化。y推得：推得：)(yfy(2)由应力分量表达式确定应力函数由应力分量表达式确定应力函数 的形式：的形式：),(yx)(22yfxy积分得：积分得：)()(1yfyxfx)()()(2212yfyxfyfx(a)(b)(),(),(21yfyfyf 任意的待定函数任意的待定函数xyllqlql1yzh/2h/2q)()(1yfyxfx)()(

21、)(2212yfyxfyfx(a)(b)(),(),(21yfyfyf 任意的待定函数任意的待定函数(3)由由 确定：确定：04)(),(),(21yfyfyf)(22)2(224yfyx044x)()()(2)4(2)4(1)4(244yfyxfyfxy代入相容方程：代入相容方程：444224442yyxx0)(2)()()(2)2()4(2)4(1)4(2yfyfyxfyfx0)(2)()()(2)2()4(2)4(1)4(2yfyfyxfyfxxyllqlql1yzh/2h/2q方程的特点：方程的特点：关于关于 x 的二次方程，且要求的二次方程，且要求 l x l 内方程均成立。内方程均

22、成立。由由“高等代数高等代数”理论，须有理论，须有x 的一、二次的系数、自由项同时为零。即：的一、二次的系数、自由项同时为零。即：0)()4(yf0)(2)()2()4(2yfyf0)()4(1yf对前两个方程积分：对前两个方程积分：GyFyEyyf231)(DCyByAyyf23)(c)此处略去了此处略去了f1(y)中的常数项中的常数项对第三个方程得：对第三个方程得：)(2)()2()4(2yfyfBAy412积分得：积分得：23452610)(KyHyyByAyf(d)GyFyEyyf231)(DCyByAyyf23)(c)23452610)(KyHyyByAyf(d)xyllqlql1y

23、zh/2h/2q)()(1yfyxfx)()()(2212yfyxfyfx(a)(b)将将(c)(d)代入代入(b)，有，有)()(223232GyFyEyxDCyByAyx)610(2345KyHyyByA(e)此处略去了此处略去了f2(y)中的一次项和常数项中的一次项和常数项式中含有式中含有9个待定常数。个待定常数。)()(223232GyFyEyxDCyByAyx)610(2345KyHyyByA(e)2.应力分量的确定应力分量的确定22yxKHyByAyFEyxBAyx2622)26()26(223222xyDCyByAy23yxxy2)23()23(22GFyEyCByAyx(f)(

24、g)(h)3.对称条件与边界条件的应用对称条件与边界条件的应用22yxKHyByAyFEyxBAyx2622)26()26(223222xyDCyByAy23yxxy2)23()23(22GFyEyCByAyx(f)(g)(h)3.对称条件与边界条件的应用对称条件与边界条件的应用（1）对称条件的应用：）对称条件的应用：xyllqlql1yzh/2h/2q由由 q 对称、几何对称：对称、几何对称：yx,x 的偶函数的偶函数xy x 的奇函数的奇函数由此得：由此得：026 FEy0232GFyEy要使上式对任意的要使上式对任意的 y 成立，须有：成立，须有：0GFExyllqlql1yzh/2h/

25、2qKHyByAyBAyxx2622)26(2232DCyByAyy23)23(2CByAyxxy（2）边界条件的应用：）边界条件的应用：(a)上下边界（主要边界）：上下边界（主要边界）：;0,2xyhy;,2qhyy;0,2yhy024823DChBhAhqDChBhAh248230432CBhhA0432CBhhA由此解得：由此解得：,23hqA,0B2qDhqC23代入应力公式代入应力公式xyllqlql1yzh/2h/2qKHyyhqyxhqx26463323223233qyhqyhqyxhqxyhqxy23623(i)(j)(k)(b)左右边界（次要边界）：左右边界（次要边界）：（由

26、于对称，只考虑右边界即可。）（由于对称，只考虑右边界即可。）,lx 未知22hyhlxxy022hyhlxx 难以满足，需借助于圣维南原理。难以满足，需借助于圣维南原理。静力等效条件：静力等效条件：轴力轴力 N=0；弯矩弯矩 M=0；剪力剪力 Q=ql；qldyQhhlxxy22022dyNhhlxx022dyyMhhlxxKHyyhqyxhqx26463323223233qyhqyhqyxhqxyhqxy23623(i)(j)(k)0K02Kh0)2646(223323hhdyKHyyhqylhqqllyhqyhqlhy2232323620)646(24322232dyHyyhqyhqlhh

27、hqhqlH1032qldylhqyhqlhh)236(2223可见，这一条件自动满足。可见，这一条件自动满足。qldyQhhlxxy22022dyNhhlxx022dyyMhhlxxxyllqlql1yzh/2h/2q)534()(622223hyhyqyxlhqx(p)截面上的应力分布：截面上的应力分布：xyxy)()(103)(3222qxlhq三次抛物线三次抛物线q22112hyhyqy22346yhxhqxy4.与材料力学结果比较与材料力学结果比较xyllqlql1yzh/2h/2q)534()(622223hyhyqyxlhqx(p)22112hyhyqy22346yhxhqxy4

28、.与材料力学结果比较与材料力学结果比较材力中几个参数材力中几个参数：截面宽：截面宽：b=1,3121hI截面惯矩：截面惯矩：静矩：静矩：2822yhS弯矩：弯矩：)(222xlqM剪力：剪力：qxQ将其代入式将其代入式(p)，有，有53422hyhyqyIMx22112hyhyqybIQSxy（3-6）xyllqlql1yzh/2h/2q53422hyhyqyIMx22112hyhyqybIQSxy（3-6）比较，得：比较，得：（1）xxy第一项与材力结果相同，为主要项。第一项与材力结果相同，为主要项。第二项为修正项。当第二项为修正项。当 h/l1，该，该项误差很小，可略；当项误差很小，可略；

29、当 h/l较大时，较大时，须修正。须修正。（2）y为梁各层纤维间的挤压应力，材力为梁各层纤维间的挤压应力，材力中不考虑。中不考虑。（3）与材力中相同。与材力中相同。注意：注意：按式（按式（3-6），梁的左右），梁的左右边界存在水平面力：边界存在水平面力：lxxX53422hyhyq说明式（说明式（3-6）在两端不）在两端不适用。适用。解题步骤小结：解题步骤小结：（1）（2）（3）根据问题的条件：几何特点、受力特点、约束特点（面力分布规根据问题的条件：几何特点、受力特点、约束特点（面力分布规律、对称性等），估计律、对称性等），估计某个应力分量某个应力分量（）的变化形）的变化形式。式。xyyx,由

30、由 与应力函数与应力函数 的关系式（的关系式（2-26），求得应），求得应力函数力函数 的具体形式（具有待定函数）。的具体形式（具有待定函数）。xyyx,),(yx),(yx（4）（5）将具有待定函数的应力函数将具有待定函数的应力函数 代入相容方程：代入相容方程：确确定定 中的待定函数形式。中的待定函数形式。),(yx04),(yx由由 与应力函数与应力函数 的关系式（的关系式（2-26），求得应），求得应力分量力分量 。xyyx,),(yxxyyx,由边界条件确定由边界条件确定 中的待定常数。中的待定常数。xyyx,用半逆解法求解用半逆解法求解梁、矩形长板梁、矩形长板类弹性力学平面问题的类弹

31、性力学平面问题的基本步骤基本步骤：应力函数法求解平面问题的基本步骤：应力函数法求解平面问题的基本步骤：（1）024422444yyxx(2-27)（2）xyyx,然后将然后将 代入式（代入式（2-26）求出应力分量：）求出应力分量：),(yx先由方程（先由方程（2-27）求出应力函数：）求出应力函数：),(yxYyxy22Xxyx22yxxy2(2-26)（3）再让再让 满足应力边界条件和位移单值条件（多连体问题）。满足应力边界条件和位移单值条件（多连体问题）。xyyx,04求解方法求解方法：逆逆解解法法（1）根据问题的条件根据问题的条件（几何形状、受力特点、边界条件等），（几何形状、受力特点

32、、边界条件等），假设各种满足相容方程（假设各种满足相容方程（2-27）的）的(x,y)的形式；的形式；（2）然后利用应力分量计算式（然后利用应力分量计算式（2-26），求出），求出 （具有待（具有待定系数）；定系数）；xyyx,（3）再利用应力边界条件式（再利用应力边界条件式（2-18），来考察这些应力函数），来考察这些应力函数(x,y)对对应什么样的边界面力问题，从而得知所设应力函数应什么样的边界面力问题，从而得知所设应力函数(x,y)可以求可以求解什么问题。解什么问题。半逆解法的数学基础：半逆解法的数学基础：数理方程中分离变量法数理方程中分离变量法。（1）根据问题的条件根据问题的条件（几何

33、形状、受力特点、边界条件等），（几何形状、受力特点、边界条件等），假设部分应力分量假设部分应力分量 的某种函数形式的某种函数形式；xyyx,（2）根据根据 与应力函数与应力函数(x,y)的关系及的关系及 ，求，求出出(x,y)的形式；的形式；xyyx,04（3）最后利用式（最后利用式（2-26）计算出）计算出 并让其满足边界条件和并让其满足边界条件和位移单值条件。位移单值条件。xyyx,半逆解法半逆解法位移分量求解位移分量求解：（1）将已求得的应力分量将已求得的应力分量（2）（3）xyyx,代入物理方程，求得应变分量代入物理方程，求得应变分量xyyx,将应变分量将应变分量xyyx,代入几何方程

34、，并积分求得位移分量代入几何方程，并积分求得位移分量表达式；表达式；由位移边界条件确定表达式中常数，得最终结果。由位移边界条件确定表达式中常数，得最终结果。1.应力函数的确定应力函数的确定(1)分析：分析：y 主要由弯矩引起；主要由弯矩引起；x 主要由剪力引起；主要由剪力引起；xy由由 q 引起（挤压应力）。引起（挤压应力）。又又 q=常数，图示坐标系和几何对称，常数，图示坐标系和几何对称，不随不随 x 变化。变化。y推得：推得：)(yfy(2)由应力分量表达式确定应力函数由应力分量表达式确定应力函数 的形式：的形式：),(yx)(22yfxy积分得：积分得：)()(1yfyxfx)()()(

35、2212yfyxfyfx(a)(b)(),(),(21yfyfyf 任意的待定函数任意的待定函数简支梁受均布载荷简支梁受均布载荷xyllqlql1yzh/2h/2q)()(223232GyFyEyxDCyByAyx)610(2345KyHyyByA(e)442244442yyxx0)(2)()()(2)2()4(2)4(1)4(2yfyfyxfyfx0)()4(yf0)(2)()2()4(2yfyf0)()4(1yfxyllqlql1yzh/2h/2q2.应力分量的确定应力分量的确定22yxKHyByAyFEyxBAyx2622)26()26(223222xyyxxy2)23()23(22GF

36、yEyCByAyx(f)(g)(h)DCyByAy233.由边界条件确定待定常数由边界条件确定待定常数xyllqlql1yzh/2h/2q附：附：应力函数确定的应力函数确定的“材料力学方法材料力学方法”要点：要点：利用材料力学中应力与梁内力的关系，假设某个利用材料力学中应力与梁内力的关系，假设某个应力分量的函数形式。应力分量的函数形式。适用性：适用性：直梁、长板条等受连续分布面力、杆端集中直梁、长板条等受连续分布面力、杆端集中力、杆端集中力偶等。力、杆端集中力偶等。应力函数常可表示为：应力函数常可表示为：)()(),(ygxfyx设法由边界面力先确定设法由边界面力先确定 其中之一，然后将其其中

37、之一，然后将其代入代入 确定另外一个函数。确定另外一个函数。)()(ygxf或04材力中，应力分量与梁内力的关系为：材力中，应力分量与梁内力的关系为：)()(2yfxQxy)()(1yfxMx式中：式中：M(x)弯矩方程；弯矩方程；Q(x)剪力方程。剪力方程。当有横向分布力当有横向分布力q(x)作用时，纵向纤维间存在挤压应力作用时，纵向纤维间存在挤压应力 ，y同时，横向分布力同时，横向分布力q(x)的挤压作用时，对轴向应力的挤压作用时，对轴向应力 也也产生影响。产生影响。x应力分量与梁内力的关系可表示为：应力分量与梁内力的关系可表示为：)()()()(21yfxqyfxMx)()(3yfxqy

38、)()(4yfxQxy考虑挤压应力影响导致考虑挤压应力影响导致然后由：然后由：xyxy222xy22yx04确定应力函数确定应力函数 的具体形式。的具体形式。例：例：悬臂梁，厚度为单位悬臂梁，厚度为单位1，=常数。求：常数。求：应力函数应力函数 及梁内应力。及梁内应力。xyObl解：解：(1)应力函数的确定应力函数的确定xQM取任意截面，其内力如图：取任意截面，其内力如图：bxQ)(0)()()(xlbbxlxM取取 作为分析对象，可假设：作为分析对象，可假设：xy)()()(ybfyfxQxy（a）f(y)为待定函数为待定函数xy由由 与应力函数与应力函数 的关系，有：的关系，有：)(2yb

39、fyx（b）对对 x 积分一次，有：积分一次，有：对对 y 再积分一次，有：再积分一次，有：)()()(321xfyfybxf)()(0yfybxfy其中：其中：dyyfyf)()(02dyyfyf)()(1（c）xyOblxQM)()()(321xfyfybxf（c）04由由 确定待定函数：确定待定函数：024422444yyxx0)()()()4(3)4(2)4(1xfyfybxf（d）要使上式对任意的要使上式对任意的x，y成立，有成立，有0)()()4(3)4(2xfyf0)()4(1yf（e）（f）由式（由式（e）求得）求得CyByAyyf231)(（g）由式（由式（f）得）得)()4

40、(3xf)()4(2yf（h）（i）积分式（积分式（h）和（）和（i）得）得2232423)(xCxBxAxf2131412)(yCyByAyf（j）（k）xyOblxQM)(223242xCxBxA)(23CyByAybx)(213141yCyByA(l)包含包含9个待定常数，由边界条件确定。个待定常数，由边界条件确定。(2)应力分量的确定应力分量的确定1121222612)26(CyByABAybxyx)23(22CByAybyxxy2222222612CxBxAxy(m)(3)利用边界条件确定常数利用边界条件确定常数xyOblxQM1121222612)26(CyByABAybxyx)2

41、3(22CByAybyxxy2222222612CxBxAxy(3)利用边界条件确定常数利用边界条件确定常数22,0byxybyylxxylxx,022,0byxybyy(o)代入可确定常数为：代入可确定常数为：0222CBA0111CBABAbC1代入式（代入式（m）得）得xyOblxQMxy0 x0yxy注：注：也可利用也可利用 M（x）=0，考虑，考虑0)()(yfxMx进行分析。此时有：进行分析。此时有：022yx)(1xfy)()(21xfxyf)(),(21xfxf为待定函数，由相容方程确定。为待定函数，由相容方程确定。llqlql1yzh/2h/2qqxxQ)(剪力：剪力：可假设

42、剪应力：可假设剪应力：)(yqxfxy)(yxfy0y)(yfy3-4 3-4 楔形体受重力和液体压力楔形体受重力和液体压力要点要点半逆解法半逆解法（因次或量纲分析法）（因次或量纲分析法）ggxyO问题的提法：问题的提法：楔形体，下部可无限延伸。楔形体，下部可无限延伸。侧面受水压作用：侧面受水压作用：g)m/N(3（水的容重）；（水的容重）；自重作用：自重作用：g)m/N(3（楔形体的容重）；（楔形体的容重）；求：楔形体应力分布规律求：楔形体应力分布规律 。xyyx,1.应力函数及应力分量应力函数及应力分量(1)分析分析：(a),gg,x 的量纲为：的量纲为：,gg)m/N(3 的形式应为：的

43、形式应为：xgygxgygx,的线性组合。的线性组合。x 的量纲为：的量纲为：2N/m(b)由由 推理得：推理得：22yx应为应为 x、y 的三次函数。的三次函数。应力函数可假设为：应力函数可假设为：3223eycxyybxaxgggyxyO(2)应力分量应力分量3223eycxyybxax考虑到：考虑到：X=0，Y=（常体力）（常体力）gcybx22Xxyx22eycx62Yyxy22yxxy2gybyax26(a)显然，上述应力函数满足相容方程。显然，上述应力函数满足相容方程。2.边界条件的利用边界条件的利用(1)x=0 （应力边界）：（应力边界）：gyxx000 xxygyey602 c

44、y0c6ge代入式（代入式（a），则应力分量为：），则应力分量为：gggyxyON2bxxy2gyxgybyaxy26(b)(2)（应力边界）：（应力边界）：tanyx 0YX0tantanxyxxyml0tantanxxyxxmlcosl其中：其中：sin将将(b)代入，有代入，有0)tan2()(bymgyl0)2()(bxmgyl0)26()2(gybyaxmbxl0)2tan6()tan2(gybyaymbyl0)2tan6(tan2gbambl0tan2bmgl)2cos(m代入，可求得：代入，可求得：gggyxyObxxy2gyxgybyaxy26(b)3cot3cot6ggaco

45、t2gb 代入式（代入式（b），有：），有：gyxyggxggy)cot()cot2cot(232cotgxyxxy(3-7)xyx)(y)(李维（李维（Levy）解答）解答沿水平方向的应力分布沿水平方向的应力分布与材力结果比较：与材力结果比较：xyxy 沿水平方向不变，在材力中无法求得。沿水平方向不变，在材力中无法求得。沿水平方向线性分布，与材力中沿水平方向线性分布，与材力中偏心受压偏心受压公式公式算得结果相同。算得结果相同。沿水平方向沿水平方向线性分布线性分布，材力中为，材力中为抛物线分布抛物线分布。gyxyggxggy)cot()cot2cot(232cotgxyxxy(3-7)李维（李

46、维（Levy）解答）解答gggyxyOxyx)(y)(沿水平方向的应力分布沿水平方向的应力分布结果的适用性：结果的适用性：（1）当坝的横截面变化时，不再当坝的横截面变化时，不再为为平面应变问题平面应变问题，其结果误，其结果误差较大。差较大。（2）假定坝下端无限延伸，可自由假定坝下端无限延伸，可自由变形。而实际变形。而实际坝高有限坝高有限，底部，底部与基础相连，有与基础相连，有地基约束地基约束，故，故底部处结果误差较大底部处结果误差较大。（3）实际坝顶实际坝顶非尖顶非尖顶，坝顶处有其它，坝顶处有其它载荷，故载荷，故坝顶处结果误差较大坝顶处结果误差较大。三角形重力坝的精确分析，常借助于三角形重力坝

47、的精确分析，常借助于有限元数值方法有限元数值方法求解。求解。工程应用：工程应用：求使坝稳定时的角度求使坝稳定时的角度 ，称为，称为安息角安息角。因次分析法（量纲分析法）：因次分析法（量纲分析法）：ggxyO楔形体，下部可无限延伸。楔形体，下部可无限延伸。侧面受水压作用：侧面受水压作用：g)m/N(3（水的溶重）；（水的溶重）；自重作用：自重作用：g)m/N(3（楔形体的溶重）；（楔形体的溶重）；求：楔形体应力分布规律求：楔形体应力分布规律 。xyyx,分析思路：分析思路：(a),gg,x 的量纲为：的量纲为：,gg)m/N(3 的形式应为：的形式应为：xgygxgygx,的线性组合。的线性组合

48、。x 的量纲为：的量纲为：2N/m(b)由由 推理得：推理得：22yx应为应为 x、y 的三次函数。的三次函数。应力函数可假设为：应力函数可假设为：3223eycxyybxax平面问题的直角坐标解答平面问题的直角坐标解答一、多项式解答一、多项式解答逆解法逆解法二、梁、长板类弹性体应力函数方法二、梁、长板类弹性体应力函数方法应力分量与梁内力的关系可表示为：应力分量与梁内力的关系可表示为：)()()()(21yfxqyfxMx)()(3yfxqy)()(4yfxQxy考虑挤压应力影响导致考虑挤压应力影响导致然后由：然后由：xyxy222xy22yx04确定应力函数确定应力函数 的具体形式。的具体形

49、式。三、三角形板、楔形体的求解方法三、三角形板、楔形体的求解方法因次分析法（量纲分析法）：因次分析法（量纲分析法）：ggxyO楔形体，下部可无限延伸。楔形体，下部可无限延伸。侧面受水压作用：侧面受水压作用：g)m/N(3（水的溶重）；（水的溶重）；自重作用：自重作用：g)m/N(3（楔形体的溶重）；（楔形体的溶重）；分析思路：分析思路：(a),gg,x 的量纲为：的量纲为：,gg)m/N(3 的形式应为：的形式应为：xgygxgygx,的线性组合。的线性组合。x 的量纲为：的量纲为：2N/m(b)由由 推理得：推理得：22yx应为应为 x、y 的三次函数。的三次函数。应力函数可假设为：应力函数

50、可假设为：3223eycxyybxax例：例：图示矩形板，长为图示矩形板，长为 l，高为，高为 h，体力不计，体力不计，试证以下函数是应力函数，并指出能解决什试证以下函数是应力函数，并指出能解决什么问题。式中么问题。式中k、q为常数。为常数。xyOlhhkxyhkxy23233解：解：(1)应力分量：应力分量：32212hkxyx022xyhkhkyyxxy236322边界条件：边界条件：02326322hkhhkhyxy02hyy显然，上下边界无面力作用。显然，上下边界无面力作用。上下边界上下边界(2)xyOlh:0 x01222322hhhhxdyhkxydy012223222hhhhxd

51、yhkxyydy左边界左边界223222236hhhhxydyhkhkydykhkyhkyhh2233232k右边界右边界:lx 01222322hhhhxdyhklydy22322212hhhhxdyhklyydy2233312hhhklykl223222236hhhhxydyhkhkydykhkyhkyhh2233232kkl结论：结论：可解决悬臂梁左端可解决悬臂梁左端受集中力问题。受集中力问题。例：例：图示矩形截面简支梁，长为图示矩形截面简支梁，长为 l，高为，高为 h，受，受有三角形分布载荷作用，体力不计。试求其有三角形分布载荷作用，体力不计。试求其应力分布。应力分布。解：解：（1）应

52、力函数形式的确定）应力函数形式的确定梁截面上弯矩和剪力为：梁截面上弯矩和剪力为：)()()()(21yfxqyfxMx)()(3yfxqy)()(4yfxQxylxqxM3)(30lxqxQ2)(20lxqxq0)(由材料力学方法可确定应力分量的由材料力学方法可确定应力分量的分离变量分离变量形式：形式：)()(320130yflxqyflxq)(30yflxq)(2420yflxq取应力分量取应力分量 分析，分析，y取应力分量取应力分量 与应力函数的关系：与应力函数的关系：y22xy)(10yflxq对此式积分：对此式积分：22xy)(10yflxq对此式积分：对此式积分：)()(210yfd

53、xyflxqx)()(22120yfyflxq)()()(3210yfdxyfdxyflxq)()()(232120yfdxyfyflxq)()()(632130yfyxfyflxq为待定函数为待定函数)(),(),(321yfyfyf（2）由相容方程确定待定函数）由相容方程确定待定函数044x)()()(6)4(3)4(2)4(13044yfyxfyflxqy)(22)2(10224yflxqyx代入代入024422444yyxx0)()(2)()(6)4(310)4(2)4(130yfyflxqyxfyflxq0)()(2)()(6)4(3)2(10)4(2)4(130yfyflqyfxy

54、flxq要使上述方程对任意的要使上述方程对任意的 x 成立，有成立，有0)()4(1yf0)(2)(210)4(2yflqyf0)()4(3yf)()()(632130yfyxfyflxq(a)(b)(c)积分式（积分式（a），得），得DCyByAyyf231)(将上式代入（将上式代入（b）积分，得）积分，得GyFyyEyByAlqyf2345023610)(积分式（积分式（c），得），得233)(KyHyyf(d)(e)(f)将求得的将求得的)(),(),(321yfyfyf代入应力函数，有代入应力函数，有（3）计算应力分量）计算应力分量22yxBAylxq26630FEyByAylxq22

55、22230KHy2622xyDCyByAylxq230yxxy2CByAylxq232220GFyEyByAylq23222340GyFyyEyByAlxq234503610DCyByAylxq2330623KyHy(g)(h)（3）利用边界条件确定待定常数）利用边界条件确定待定常数上边界上边界：2hy xlqy00 xyDhChBhAlxq222230CBhhAlxq2202320412322340GFhEhBhAhlq124823DhChBhA0832CBhhA041232234GFhEhBhAh(i)(j)(k)下边界下边界：2hy 0y0 xy024823DhChBhA0832CBhh

56、A041232234GFhEhBhAh(l)(m)(n)左边界左边界：0 x0220hhxxdy0220hhxxydy0220hhxxydy0)26(22hhdyKHy0)26(222hhdyKyHy02322222340dyGFyEyByAylqhh左边界左边界：lx 022hhlxxdyMydyhhlxx22620lqQdyhhlxxy2220lq(o)(p)(q)(r)(s)(t)联立求解式（联立求解式（i）（t）,可得具体的可得具体的应力分量应力分量。注：注：位移边界条件转化为应力边界条件。位移边界条件转化为应力边界条件。（1）（2）试按材料力学中确定应力的方法，写出试按材料力学中确定

57、应力的方法，写出图示两梁所有应力分量形式。（含有待定函图示两梁所有应力分量形式。（含有待定函数）数）课堂练习：课堂练习：3-5 3-5 级数式解答级数式解答问题的提出问题的提出多项式解答：多项式解答：只能求解载荷只能求解载荷简单简单，且，且连续分布连续分布的问题。的问题。不能求解载荷不能求解载荷复杂复杂，且，且间断分布间断分布的问题。的问题。级数式解答：级数式解答：),(yx其基本思路是将应力函数其基本思路是将应力函数 分解成关于分解成关于 xy 的的两个单变量函数的乘积两个单变量函数的乘积。分离变量法。分离变量法。（属逆解法）（属逆解法）1.级数形式的应力函数级数形式的应力函数假设：假设：)

58、(sin),(yfxyx(a)式中：式中：为任意常数，其量纲为为任意常数，其量纲为 ,1长度)(yf为为 y 的任意（待定）函数。的任意（待定）函数。)(sin)2(2224yfxyx)(sin444yfxx)(sin)4(44yfxy将其代入将其代入 ：04载荷载荷复杂复杂，且，且间断分布间断分布的问题，可由的问题，可由级数式级数式解答解决。解答解决。有：有：)(cos),(1yfxyx442244442yyxx)(sin)(sin2)(sin)4()2(24yfxyfxyfx0)()(2)(sin4)2(2)4(yfyfyfx(b)解上述方程，得解上述方程，得0)()(2)(4)2(2)4

59、(yfyfyf其中：其中：A、B、C、D 都是任意常数，都是任意常数，将其代入应力函数将其代入应力函数 ，得，得ychyshyBchysh)(DyCyAyf(c)再取如下应力函数：再取如下应力函数：y)chyshyBchysh(sin),(DyCyAxyx式中：式中：也为任意常数也为任意常数,为为 y 的任意（待定）函数。的任意（待定）函数。)(1yf类似于上面的运算，可得应力函数的另一解类似于上面的运算，可得应力函数的另一解：(d)显然，将式显然，将式(c)与与(d)相加，仍为可作为应力函数：相加，仍为可作为应力函数：)chshchsh(cos),(yyDyyCyByAxyx1)chshch

60、sh(cosmmmmmmmmmmyyDyyCyByAx)chshchsh(sin),(yDyyCyyByAxyx(e)取取 和和 的一系列值，即取：的一系列值，即取：m,m)(m将由此构成的将由此构成的 加起来，有加起来，有),(yx)chshchsh(cosyyDyyCyByAx1)chshchsh(sinmmmmmmmmmmyyDyyCyByA(3-8)显然，式显然，式(3-8)满足相容方程，可作为应力函数。且在其上满足相容方程，可作为应力函数。且在其上再加若干个满足相容方程的应力函数，仍可作为应力函数。再加若干个满足相容方程的应力函数，仍可作为应力函数。2.级数形式的应力分量级数形式的应

61、力分量将上述应力函数将上述应力函数 代入应力分量表达式（代入应力分量表达式（2-26），有），有),(yx12ch)2(sh)2(cosmmmmmmmmmmmyCByDAx1222ch)2(sh)2(sinmmmmmmmmmmmxyCByDAxyyyDyyCmmmmchshyyDyyCmmmmchsh122myx12mxyyx(3-9)式（式（3-9）满足相容方程、平衡方程，）满足相容方程、平衡方程，只要适当选取：只要适当选取：mmmmmmmmmmDCBADCBA,;,使其满足边界条件，即为某问题的解。使其满足边界条件，即为某问题的解。3-6 3-6 简支梁受任意横向载荷简支梁受任意横向载荷边

62、界条件边界条件1.边界条件的级数表示边界条件的级数表示上下边界：上下边界：0Hyxy)(0 xqyy)(1xqHyyRdyHxxy00左右边界：左右边界：00yxy0lxx00 xx10RdyHlxxy（a）（b）（c）（d）由边界条件（由边界条件（c），得），得;0mmmmDCBA),3,2,1(mlmm此时应力分量式（此时应力分量式（3-9）简化为）简化为1222ch)2(sh)2(sinmmmmmxlymmlCBlymmlDAlxmlmlymyDlymyCmmchshlymyClymBlymAlxmlmmmmmyshchshsin1222lymyDmch1222ch)(sh)(cosmm

63、mmmxylymmlDAlymmlCBlxmlmlymyClymyDmmchsh（3-10）将此应力分量式（将此应力分量式（3-10）代入边界条件（）代入边界条件（b），有），有0cos12mmmDmlAlxmm（e）12ch)(sh)(cosmmmmmlHmmlDAlHmmlCBlxmm0chshlHmHClHmHDmm（f）0Hyxy00yxy（b）0mmDmlA（i）lHmHlHmmlClHmBlHmAmmmchshshch0shchlHmHlHmmlDm（j）)(sin1222xqBlxmmlmm（g）)(0 xqyy)(1xqHyy（a）1m222chshsinmmlHmBlHmAl

64、xmml)(chsh1xqlHmHDlHmHCmm（h）将此应力分量式（将此应力分量式（3-10）代入边界条件（）代入边界条件（a），有），有将将)(),(1xqxq在区间（在区间（0，l）上展为和等式左边相同的级数，即）上展为和等式左边相同的级数，即lxmsin的级数，由的级数，由Fourier级数的展开法则，有级数的展开法则，有lxmdxlxmxqlxqmlsinsin)(2)(10lxmdxlxmxqlxqmlsinsin)(2)(1011（3-11）比较式（比较式（3-11）与式（）与式（g）和（）和（h）两边的系数，有）两边的系数，有lmdxlxmxqlBml0222sin)(2lm

65、dxlxmxqmlB022sin)(2（k）lHmBlHmAmchshmlHmHDlHmHCmmchshldxlxmxqml0122sin)(2（l）由式由式 (i)、(j)、(k)、(l)可求得全部和系数：可求得全部和系数：，代入式，代入式（3-10）求得应力分量。）求得应力分量。mmmmDCBA,说明：说明：（1）边界条件（边界条件（d）在求解中没有用到，但可以证明是自动满足的。）在求解中没有用到，但可以证明是自动满足的。（2）级数求解计算工作量很大，通常由有关计算软件求解，如：级数求解计算工作量很大，通常由有关计算软件求解，如：MathCAD、Matlab、Mathematica等。等。

66、（3）结果在梁的端部误差较大；另外，当梁的跨度与高度相当时结结果在梁的端部误差较大；另外，当梁的跨度与高度相当时结果误差也较大。果误差也较大。弹性力学平面问题的基本理论弹性力学平面问题的基本理论小结小结一、两类平面问题及其特征一、两类平面问题及其特征名名 称称平面应力问题平面应力问题平面应变问题平面应变问题未知量未知量已知量已知量未知量未知量已知量已知量位位 移移应应 变变应应 力力外外 力力几何形状几何形状xyyx,0zxyz)(yxzExyyx,0zxyz0zxyyx,0zxyz0zxyyx,0zxyz)(yxzvu,0wvu,0w体力、面力的作用面都平行于体力、面力的作用面都平行于xoy平面，且沿板厚不变化。平面，且沿板厚不变化。体力、面力的作用面都平行于体力、面力的作用面都平行于xoy平面，且沿平面，且沿 z 向不变化。向不变化。z 方向的尺寸远方向的尺寸远小小于板面内的于板面内的尺寸（等厚度薄平板）尺寸（等厚度薄平板）z 方向的尺寸远方向的尺寸远大大于于xoy平面平面内的尺寸（等截面长柱体）内的尺寸（等截面长柱体）二、平面问题的基本方程二、平面问题的基本方程（1）平衡微分方程