上海交通大学：《弹塑性力学》课程教学资源（讲义）弹性与塑性力学解题方法（2/3）

二、平面问题的基本方程 平面应力问题： 00x + 应力边界条件 >平衡方程 0c匹+f:=0 ay 0t对 2+j,=0 + lox+mt=Fx ay Ou Itx+moy =Fy >几何方程： v sy-by 1g= ax ay E o:-uo3] 平面应变问题： >物理方程 6y= El,-uo.] 1 ,E→1- μ→1-4 E Yo-GTo 国上大 MB6011弹性塑性力学 25 三、相容方程和AirY应力函数 du 00x+ s+fx=0 6x= 8-Eox-o,] ax ay av OTx ,= 6,= + Ox a0=0 ay El,-uos] ay = + 8x ay Yo G a28x+ ap2 相容方程（应变协调方程） 8x2 axay +6.+o,-+( +r） ay ■应力分量表示的相容方程 82 62 6.+o,）-1+( ax+ 体力为常量：又2口x+o,)=0 圈上清大生 ME6011弹性塑性力学 26 1

1 ME6011 弹性塑性力学 平衡方程   0      x x yx f x y     0      y xy y f x y   二、平面问题的基本方程                      y u x v y v x u xy y    x 几何方程：                   xy xy y y x x y G E E         1 1 1 x 物理方程 x m yx Fx l    xy m y Fy l    应力边界条件 平面应变问题： 2 1 , 1         E E 25 平面应力问题： ME6011 弹性塑性力学 三、相容方程和Airy应力函数 相容方程（应变协调方程） 0 0                      y xy y x x xy f x y f x y                          y u x v y v x u xy y    x                   xy xy y y x x y G E E         1 1 1 x y x x y x y xy             2 2 2 2 2                                 y f x f x y x y x y (1 ) 2 2 2 2    应力分量表示的相容方程                     y f x f x y x y (1 ) 2    2 2 2 2 2 x  y       体力为常量：   0 2   x   y  26

00x Ox ots+fx=0 +不含材料常数，只要两物 ay 体具有相同的形状，受相 Ofx Ox 三+f,=0 同面力，则不论何种材料 ay 何种平面问题，其应力分 82 Or?+ 82 布是相同的，数值亦相同。 x+o,)=0 解答普遍性。 卫两末有裤系数非齐次线性衡撒分方程，其通解等于对应的齐次通解加 00x+ OTx =0 特解： Ox ay 0y-0 Ox=-fx·x Ox=0 Ox ay o,=-f,.y 0y=0 =0 g =-yfx-xfy 国上唐成大生 ME6011弹性塑性力学 27 00x+ 0rg=0 00:0(-1g) 0A 0A Ox ay 0x= ay Tx = Ox OTx 80y=0 a(-tg)_0o上 aB aB Ox dy y= Ox Ox ay 存在一个函数仰(K》 (A= bp by 8'o 0o 8p Ox= ay2 0y= ar? Tx=- B= ap axoy 0x 4齐次通解 4全解： a2中 86 0= -xf. dy2 0,= -y, To=- OxCy 8o 8o 0= 02 6, dx? To=- OxOy -才-可 。p化，)：平面问愿的应力函数：不论伞（化，》取何函数，上式得到的应力分量恒满 足平衡徽分方程。(Airy stress function). :上潘久毛大警 ME6011弹性塑性力学 2

2 ME6011 弹性塑性力学 0 0                    x y x y xy y x xy       0 2 2 2 2                x y x y   0 0                      y xy y x x xy f x y f x y     不含材料常数，只要两物 体具有相同的形状，受相 同面力，则不论何种材料、 何种平面问题，其应力分 布是相同的，数值亦相同。 解答普遍性。 前两式为常系数非齐次线性偏微分方程，其通解等于对应的齐次通解加 上一个特解。 特解：  0            xy y y x x f y f x    0 0           xy x y y x  yf xf   27 ME6011 弹性塑性力学 0 0                    x y x y xy y x xy     ( ) x y x xy        ( ) x y xy y        xy x A y A x          xy y B x B y          B x y A         y x x y x y xy                  2 2 2 2 2 x,y) :平面问题的应力函数：不论  x,y) 取何函数，上式得到的应力分量恒满 足平衡微分方程。 （Airy stress function） 齐次通解 全解： 22 2 2 2 x x y y xy xf yf y x xy                22 2 2 2 x y xy x y yf xf y x xy              28 存在一个函数x,y)

72a+c,)=0 + + a+ 82)2 =0 7272p=0 ÷应力函数表示 的协调条件 讨论： lox+mty=Fx lts+moy=Fy 。应力函数（化y为双调和函数。满足P2P2p=0相当于满足平衡微分方程 和变形协调条件，即满足了平面问题的八个基本方程。 求出P(化，后，可根据通解求出应力分量，如果在边界上满足应力边界条件， 则得到的就是正确解答。 。求出应力分量后，由物理方程求应变分量，再由几何方程求位移分量。 国上本大学 ME6011弹性塑性力学 29 。当体力为零时： Ox= o'g 2 0,= 0'g dr? 8'p g= axay Q平面应变问题： →， 1-4 ,E→1-μ 07272p=0 为四阶偏微分方程，直接求解比较困难，故常用逆解法和 半逆解法。 逆解法： 设定印(y) 中 边界条件>解决的问题 满足口4p=0 半逆解法： 解决的问题 边界形状→设定 →p(x,y 叉4p=0 x→正确解答 受力情况 边界条件 :上潘久毛大学 ME6011弹性塑性力学 30 3

3 ME6011 弹性塑性力学 0 2 2 2 2 2 2 2 2                             x y x y 0 2 2     应力函数表示 的协调条件   0 2   x   y  0 2 2 2 2 2                x y 应力函数x,y)为双调和函数。满足 2 2 = 0 相当于满足平衡微分方程 和变形协调条件，即满足了平面问题的八个基本方程。 求出x,y)后，可根据通解求出应力分量，如果在边界上满足应力边界条件， 则得到的就是正确解答。 求出应力分量后，由物理方程求应变分量，再由几何方程求位移分量 。 讨论： x m yx Fx l    xy m y Fy l    29 ME6011 弹性塑性力学                     x y x y xy y x       2 2 2 2 2 当体力为零时： 平面应变问题： 2 1 , 1         E E  2 2 = 0 为四阶偏微分方程，直接求解比较困难，故常用逆解法和 半逆解法。 逆解法： 设定x,y) 满足 4 = 0 半逆解法：      xy y x    解决的问题 解决的问题 边界形状 受力情况      xy y x    x,y) 4 = 0 边界条件 设定 正确解答 边界条件 30

≤例题 图示很长的矩形柱体，材料的比重为Y,将其放入 形状相同的刚性槽内若不考虑摩擦力，设应力函数 的形式为0=Ar2y+By2+Cy2+Dx2 试求各应力分量、应变分量以及位移分量。 因上海大华 ME6011弹性塑性力学 31 解： 根据Airy应力函数可得通解 0o 0x= =6By+2C oy2 8o 0,= dr? =2Ay+2D 8o =-2Ax axdy f=0,f=-) 特解为：O=0,O,=0,Ty=yx 【w=-f-对， 所以全解为： 0,=6B+2C,0,=2Ay+2D, Ty=-2Ax+yx 圈上本大学 ME6011弹性塑性力学 32 4

4 ME6011 弹性塑性力学  例题 图示很长的矩形柱体，材料的比重为γ，将其放入 形状相同的刚性槽内若不考虑摩擦力，设应力函数 的形式为 2 3 2 2   Ax y  By Cy  Dx 试求各应力分量、应变分量以及位移分量。 h a a x y o 31 ME6011 弹性塑性力学 解： 根据Airy应力函数可得通解 32                             Ax x y Ay D x By C y xy y x 2 2 2 6 2 2 2 2 2 2       特解为： 0 0 x y xy   ， ，   x ( 0, ) x y f f    xy x y   yf  xf 所以全解为： 6By 2C,  x   2Ay 2D,  y   2  xy    Ax x 

应力边界条件y=h处，tm=0→A=y/2 0,=0→D=-h/2 刚性槽的条件 6,k=0三和变量x无关6，=0 ￡=1+40-4)a,-uo,]=0→o,=1- o:=r(y-h) a0-创一B片名c=n变 1-46 1-42 1+0-2p0y-m) 8,=1+0-0o,-uo,=E1-4 u=J8,= y=0时，v=0 v=∫e,dy 11+I-2四-例+K ◆K=0 E 1-4 国上活庆大峰 ME6011弹性塑性力学 33 §4一5逆解法和半逆解法 求解弹性力学问题在数学上是比较复杂的，因此不 得不采用逆解法和半逆解法。 优点：这样求解，在数学上比较容易； 缺点：带有一定的偶然性和相当大的局限性，有时 需要进行多次反复试算。 圈上清大峰 ME6011弹性塑性力学 34 5

5 ME6011 弹性塑性力学 应力边界条件 y  h 处， 0 0   y xy   刚性槽的条件  0  a a x  dx  x 和变量 x 无关  0 x  A   2 D  h 2 x  y      1 [(1 ) ] 0 1     x  x y E      ( y h)  y    ( ) 1 y h x        1 6      B  1 2 h C        ( ) 1 1 (1 )(1 2 ) [(1 ) ] 1 y h E E y y x                     0  u dx x  hy K y E v dy y         ) 2 ( 1 1 (1 )(1 2 ) 2      y  0时，v  0 K  0 33 ME6011 弹性塑性力学 34 §4－5 逆解法和半逆解法 求解弹性力学问题在数学上是比较复杂的，因此不 得不采用逆解法和半逆解法。 优点：这样求解，在数学上比较容易； 缺点：带有一定的偶然性和相当大的局限性，有时 需要进行多次反复试算