《2022年第六章弹性力学平面问题的直角坐标系解答》由会员分享，可在线阅读，更多相关《2022年第六章弹性力学平面问题的直角坐标系解答（18页珍藏版）》请在金锄头文库上搜索。
1、优质资料欢迎下载第六章弹性力学平面问题的直角坐标系解答在第五章讨论了弹性力学问题的基本解法：位移法和应力法， 并结合简单的三维问题，根据问题的特点，猜想问题的应力解或位移解，并验证猜想的解是否满足应力法或位移法的基本方程和边界条件，满足则为问题真解。弹性体都是三维的，而受力（外力）一般也是空间力系，但如果所研究弹性体具有某种特殊形状。例如：一个方向的尺寸远比其他两个方向的尺寸大的多或小的多，并且承受某种特殊规定的外力和约束，则可以把空间问题（三维）作为近似的平面问题（二维）处理，这将使分析和计算大大简化，而所得结果也能满足工程上对精度的要求。第 1 节平面问题的分类平面问题在工程中极为常见，而
2、且平面问题的解析解在整个弹性力学解析解中占有较大比重。因此必须给予足够的重视。平面问题分为平面应力问题和平面应变问题两类。下面将它们分类简要说明一下。1.1平面应力问题固体的形状特点：物体一个方向尺寸比其它两个方向尺寸小的多（等厚度薄板） 。x2x1x3o x2t 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 1 页，共 18 页优质资料欢迎下载受力和约束特点：沿厚度（x3方向）均匀分布，体力f3 = fz= 0, 面力03ZX，在薄板表面无面力，坐标系（ x1 , x2 , x3）放在板厚中间平面中平面，以 z(或 x3）轴垂直板面。满足上述条
3、件的问题称为平面应力问题。由物体几何特点和受力特点知：在2tz处，0ZYXz=zx=zy=0。由于薄板很薄，表面三个应力分量为零，则近似认为在 V 内认为z=zx=zy=0。平面应力问题：应力分量仅存三个x=x(x,y), y=y(x,y),xy=xy(x,y),均为 x,y的函数，待求。将应力分量代入各向同性材料的本构关系存在四个应变分量 (待求量)：x , y , xy ,z（其中z不独立）位移分量待求量：u(x,y) , v(x,y) （考虑平面内位移） . 平面应力问题待求未知函数一共八个：3 个应力 3 个应变 2 个位移1.2平面应变问题形状特点：物体一个方向尺寸（z 或 x3）
4、比其它两个方向（x,y或 x1 ,x2）大的多，如水坝、涵洞。x1(x) x2(y) x3(z) 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 2 页，共 18 页优质资料欢迎下载受力和约束情况：沿z(或x3）轴方向无变化，体力f3 = fz = 0，面力03ZX，这样x3 = z = const面均可看成对称面，对称结构受对称荷载和约束，则此对称面处的位移和变形为零，即w=0(z=0),zx=zy=0所以平面应变问题：应变分量仅有三个x,y,xy=yx，位移分量两个：u(x,y) , v(x,y)，应力分量：x, y,xy,z（其中z 不独立）
5、 。平面应变问题待求未知函数仍然八个：3 应力 3 应变2 位移。第 2 节平面问题的基本方程和边界条件2.1 平衡微分方程（ 2 个）两个平面问题一致：, +f =0, , =1,20Xyxyxx，0Yyxyxy2.2几何方程（ 3 个）两平面问题一致：)(21,uuxux, yvy，xvyuxy2.3相容方程（ 1 个）精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 3 页，共 18 页优质资料欢迎下载两平面问题一致：yxxyxyyx22222对于平面应力问题还应有022xz，022yz，02yxz，但对于薄板厚度尺寸远此三个方程可以不考虑。2
6、.4本构方程（ 3 个）平面应力问题)(1yxxE，)(1xyyE，xyxyE)1(2平面应变问题)1()1(2yxxE，)1()1(2xyyE，xyxyE)1(2两个平面问题的基本方程仅物理方程有所不同，将平面应力物理方程中弹性系数21EE，1，则平面应力问题的物理方程变为平面应变问题的物理方程。所以按平面应力问题求解的结果中弹性系数也如此替换，则可得到平面应变问题解。2.5边界条件位移边界条件：uu( =1,2）uu，vv在Su上力的边界条件：nXyxxmlX，yxymlY在S上精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 4 页，共 18
7、页优质资料欢迎下载第3节平面问题的基本解法3.1 位移法基本未知函数：u(x,y) , v(x,y)基本方程两个：用u , v表示的平衡微分方程。平面应力问题：011,2fuGuG其中22222yx平面应变问题：0211,2fuGuG边界条件：位移边界uu，vv在Su上力的边界yxxmlX，yxymlY在S上（应力需要用位移微分表示）3.2应力法基本未知函数（ 3 个） ：x, y,xy=yx,基本方程（ 3 个） ：2 个平衡微分方程, +f =01 个相容方程：平面应力问题时)(1()(2yfxfyxyx平面应变问题时)(11)(2yfxfyxyx力边界条件：精选学习资料 - - - -
8、- - - - - 名师归纳总结 - - - - - - -第 5 页，共 18 页优质资料欢迎下载nXyxxmlX，yxymlY在S =S上当体力为常数或体力为零时，两个平面问题的相容方程一致2(x+y )=0 (x+y )为调合函数，与弹性系数无关，不管是平面应力（应变）问题，也不管材料如何，只要方程一致，应力解一致，有利实验。3.3 应力函数解法当体力为常量或为零时，按应力法解的基本方程（共三个）为, +f =0 2=0应力法基本方程的前两个为非齐次方程，所以根据微分方程理论， 非齐次微分方程的通解等于其特解加上齐次微分方程的通解。所以,非齐次方程特解可以选x= -fxx, y= -fy
9、y ,xy= 0;（特解还可以选其它形式）下面工作求齐次微分方程, =0的通解，或者0yxyxx，0yxyxy的通解同时通解还需要满足相容方程：2(x+y )=0 对于上面三个齐次微分方程要求出其通解，仍是一个较复杂、困难的问题。1862 年Airy提出将满足三个齐次微分方程的3 个应力分量的齐次解精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 6 页，共 18 页优质资料欢迎下载由一个函数（应力函数）的二阶微分来表示，使之自然满足齐次平衡微分方程, =0 这样应力法的齐次基本方程仅为用应力函数表示的相容方程，使未知函数和基本方程数均减为一个。Ai
10、ry提出应力函数(x,y) 与齐次微分方程中待求应力分量之间满足如下微分关系：22yx，22xy，yxxy2（a）应力函数(x,y) 与待求应力分量齐次解之间的微分关系是由两个齐次平衡微分方程导出的：xxyxyAyx，xyxAyxyyxBxy，xyyB得yAyBxA，xB从而导出 (a)式。则 (a) 式使得齐次的平衡微分方程自然满足，将(a) 式代入相容方程，得0)(42222222xy精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 7 页，共 18 页优质资料欢迎下载应力函数解法的基本方程（一个）基本方程为由应力函数满足的双调合方程。最后应力分
11、量解为其特解加通解：xfyxx2，yfxyy2，yxxy2在边界上应力分量满足力的边界条件（在S上） ，用应力函数表示：)()(222yxmxfylXx，)()(222yfxmyxlYy对于单连域，应力函数(x,y)满足双调和方程4=0，且在S上满足用应力函数二阶偏微分表示的边界条件，则由(x,y)导出应力分量为真解，对于复连域，还要考虑位移的单值条件。3.4 应力函数的特性1.应力函数加上一个线性函数a+bx+cy，并不影响应力，换句话说，某问题的应力函数为，则1=+a+bx+cy也是问题的应力函数。应力函数可确定到只差一个线性函数。2.无体力作用时， 应力函数及其一阶偏导数的边界值可分别由
12、边界的面力的主矩和主矢量来确定yBAABRdSYxx)()(, xBAABRdSXyy)()(BBABBABABMdSXyydSYxx)()(对 B 点取矩)逆时针为正。y y x o ds dy n e1e2-dx x o A B F 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 8 页，共 18 页优质资料欢迎下载下面推导一下 . 对于无体力时fx=fy= 0；力的边界条件为Xyxmyl222, Yxmyxl222而dSdyenl),cos(1, dSdxenm),cos(2代入上式得XdSdxyxdSdyy222XydSd)(YdSdxxd
13、Sdyyx)(222YxdSd)(积分得xBABAABRdsXydSdyy)()()(yBAABRdsYxx)()(根据函数的求导公式dSdyydSdxxdSddSXdSdydSYdSdxdSdCACA)()(而 C 为边界上动点上式对 s 积分得dSdSXdSdydSYdSdxBACACAAB)(采用分部积分精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 9 页，共 18 页优质资料欢迎下载dSXyydSYxxdSXydSYxdSXydSYxdSXydSYxdSXydSYxBABBABBABABABABBBABABACACA)()()()(边界力
14、对 B 点之矩第4节多项式应力函数运用举例例题 1 矩形域无体力作用时应力函数分别为二次项和三次项的结果（而一次项无须考虑） ，采用逆解法。1取为二次项：2322122),(ycxycxcyx，代入4=0， 满足。将代入应力分量与应力函数的关系式，得322cyx，122cxy，22cyxxy可见，矩形域各点应力状态一样，为常量。设c1,c2,c3均为正值。矩形域边界面力如图所示。3.取为三次项：342322316226),(ydxydyxdxdyx，c1x c3y c21 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 10 页，共 18 页优质资
15、料欢迎下载代入4=0， 满足。将代入应力分量与应力函数的关系式，得ydxdyx4322ydxdxy2122应力为 x、y的线性式。ydxdyxxy322仅取一项346ydx= d4y, y=xy= 0。在边界上面力分布与坐标系位置有关。坐标系如下图所示面力分布为纯弯问题，在两端面的面力将产生一个M 。12334223422hdydydxMhhhhxzIMhMd1234，yIMzx（材料力学解）由应力分量求应变分量：yEIMExx1，EIMyExy，0xy；通过几何方程积分及约束条件可以求出位移。本题讨论坐标位置选取不同将导致边界上面力分布不同，从而对应不同y x h/2 h/2 1 d4 h/
16、2 M M 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 11 页，共 18 页优质资料欢迎下载的问题。因此，本题在 边界上面力分布与坐标系位置有关。如果仍取一项346ydx= d4y, y=xy= 0；但坐标位置变了 , 边界上面力分布如下图。例题 2 无体力作用的悬臂梁，在端部受集中力P 作用。本题采用应力函数的半逆解法。半逆解法思路：1.根据受力情况和求解经验， 包括材料力学的解， 定性估计应力分量的变化， 并根据应力分量与应力函数关系，反推出函数的主要项。2.将所设代入4=0和力的边界条件进行检验，如果不满足则进行修正（适当增加项），再代
17、入4=0和力的边界条件进行检验，直至满足所有方程为止。本题求解的基本情况：基本方程4=0，边界条件为混合边界条件：主要边界上，在y= h：0X，0Y（无面力）次要边界上，h y x d4h d4h/2 d4h/2 x 1 y P M P l x 2h 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 12 页，共 18 页优质资料欢迎下载在x=l：0X，hhPdyY在x=0： 严格要求u=0，v=0PdyYPlMdyXdyXhh0解题：1根据受力特点知在x处弯矩：M=P(l-x)，材料力学应力解：yIxlPyIMx)(所以x包含y和xy项，又因为22
18、yx，可设313166ycxya，代入4=0， 满足。将代入应力分量与应力函数的关系式，得2121120yayxycxyaxyyx在y= h时，212haxy为均匀剪力。将应力分量代入边界条件：代 入 主 要 边 界 ：y= h， l=0,m=1精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 13 页，共 18 页优质资料欢迎下载0000yxyYX，如果满足，则a1=0 。代回应力分量表达式001xyyxyc本应力解对应纯弯问题，不是所要求的。2对要进行修正，消去y= h面上均匀剪力设313166ycxya+b1xy，代入4=0，满足。将代入应力分
19、量与应力函数的关系式，得1211120byaycxyaxyyx代入主要边界：y= hy= 0满足；xy= 0 02121bha，或2112hab；代回应力分量表达式精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 14 页，共 18 页优质资料欢迎下载)(2022111hyaycxyaxyyx将应力分量代入次要边界，代入x=l边界：l=1 , m=0,则0xX011yclya， 或lac11而PdyYPdyhyahh)(22213123hPa；2123hPlc，hPb431将惯性矩323hI代入a、b、c表达式，则IPa1，IPhb221，IPlc1
20、代回应力分量表达式yxlIPx)(，0y，)(222yhIPxy与材料力学解相同。注意本题应力解在梁两端不能用。因为用到了圣维南原理。有了应力解后，依次求应变和位移。在位移的确定中，当x=0 , u=v=0不能处处满足，而用到0000000yxyxyxdxdvvu， 将刚体位移去掉， 放松了位移边界处理。例题 3 简支梁（不计体力）上面受均载作用，仍采用应力函数解的半精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 15 页，共 18 页优质资料欢迎下载逆解法。考虑应力特点：y与x无关，y由q引起，且在2hy处y为常数。设)(22yfxy)()(1y
21、fyxfx)()()(2212yfyxfyfx代入基本方程4=002ffxffx21 (4)2(4)1(4)2微分方程对全梁满足。因此。要求044dxfd，0414dxfd，0222424dxfddxfd由前两个常微分方程积分得到f(y)和f1 (y)的表达式，代回第三个常微分方程积分，可得到f2 (y)的表达式。所有待定系数由边界条件定。例题 4 楔形体受重力和液体压力作用，楔形体下端无限长。x 1 y ql ql l h q l x n g gy 精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 16 页，共 18 页优质资料欢迎下载楔形体的体积
22、力fx= X = 0，fy= Y = g；边界条件 : 在x=0处, gyX，0Y，则边界处的应力为x= - gy，xy =0在x = ytg处,0X，0Y从楔形体的受力情况分析,可以认为在楔形体y=c截面上内力为受压力和弯曲组合 , 应力分量y ,xy为x,y的线性项，x为 y 的线性式，由应力分量与应力函数的关系式,可设应力函数为x,y的纯三次式取=ax3+bx2y+cxy2+ey3将代入双调合方程4=0，满足。将代入应力分量与应力函数的关系式，得cybxyxgybyaxgyxeycxyxyyx22266222222a,b,c,d四个系数由四个力边条件确定。由x=0时：l= -1, m=0 ,gyX，0Y得02602cygyeyc解得e= - g/6，c =0，精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 17 页，共 18 页优质资料欢迎下载在x = ytg处：l = cos, m = cos(900+)= -sin, 0X，0Y得22ctggb，336ctggctgga将待定系数的解代回应力表达式,可得到楔形体的应力解。精选学习资料 - - - - - - - - - 名师归纳总结 - - - - - - -第 18 页，共 18 页
