第七章 扭转-01

1、第七章第七章 柱体扭转柱体扭转圆轴扭转平面假设非圆截面柱体横截面翘曲柱体扭转精确求解是十分困难的 目录目录7.1 7.1 扭转问题的位移解法扭转问题的位移解法7.2 7.2 扭转问题的应力解法扭转问题的应力解法7.4 7.4 椭圆截面杆件的扭转椭圆截面杆件的扭转7.3 7.3 薄膜比拟法薄膜比拟法7.5 7.5 开口薄壁杆件的扭转开口薄壁杆件的扭转7.6 7.6 闭口薄壁杆件的扭转闭口薄壁杆件的扭转7.1 扭转问题的位移解法扭转问题的位移解法柱体扭转柱体扭转自由扭转自由扭转翘曲不受限制翘曲不受限制约束扭转约束扭转翘曲受到限制翘曲受到限制弹性力学讨论自由扭转自由扭转 7.1 位移解法位移解法 2

2、关于圆截面直杆的扭转问题，材料力学的解答，即；剪应力与点到圆心的距离r成正比且与之垂直，变形后横截面仍保持为平面，是正确的。这一结论是库仑于1784年提出的。后来法国人纳维叶（Navier）将这一结果应用于非圆截面直杆的扭转，但这样作的结果是错误的，这可由受扭矩形截面杆的剪应力与边界条件不符合得到证明。 7.1 位移解法位移解法 2对于柱体的自由扭转，假设柱体的位移约束为固定左端面任意一点和相应的两个微分线素，使得柱体不产生刚体位移。柱体右端面作用一力偶T，侧面不受力。设柱体左端面形心为坐标原点，柱体轴线为z 轴建立坐标系。柱体扭转时发生变形，设坐标为 z 的横截面的扭转角为a a ，则柱体单

3、位长的相对扭相对扭转角转角为而横截面的扭转角而横截面的扭转角a a = = j j z zaj= 7.1 位移解法位移解法 3根据观察，对柱体内部位移作以下的假设：1. 刚截面假设: 柱体扭转当横截面翘曲时，它在Oxy平面上的投影形状保持不变，横截面作为整体绕 z 轴转动，如图所示。 当扭转角 a 很小时，设OP=r，则P点的位移为,coscoscoscossinsincossin,sinsinsincoscossincoscoszyu x yv x yzxrarraarrararraajjrra= =2 横截面的翘曲位移与单位长度的相对扭转角j成正比，而且各个截面的翘曲相同,即w=jF (x

4、，y)。 7.1 位移解法位移解法 4 1855年，圣维男（Saint-venant）采用半逆解法，从分析位移入手，提出了等截面直杆扭转的解答。 F(x，y)称为圣维南(Saint Venant)扭转函数，或者称为翘曲函数。柱体自由扭转计算模型柱体自由扭转计算模型l扭转假设：xzvyzujj=),(yxwj=1. 1. 刚截面刚截面假设假设 扭转与圆轴扭转变形相同扭转与圆轴扭转变形相同2. 2. 翘曲翘曲假设假设各个截面的翘曲相同各个截面的翘曲相同 7.1 位移解法位移解法 5根据上面的分析可见，对非圆截面柱体扭转问题的根据上面的分析可见，对非圆截面柱体扭转问题的位移分量位移分量u,v,wu,

5、v,w仅是坐标仅是坐标x,yx,y的二维函数。的二维函数。根据几何方程，应变分量为根据几何方程，应变分量为根据本构方程，应力分量为根据本构方程，应力分量为 ,0 xzxyzyyzxyxxyjjFF=,0 xyyyzzxzxGyGxxyjj=FF=,( , ),yzxzyuvwxjjj= =1.平衡微分方程000=bzzyzzbyzyyxybxzxyxxFzyxFzyxFzyx0yzxzxy=不计体力时， 7.1 位移解法位移解法 7l柱体扭转边界条件对于自由扭转问题，在侧边界没有载荷作用。由于 ，只有 和 不等于零。下面分为柱体侧面和端面两部份面力边界条件讨论；0 xyzxy=xzyz2. 面

6、力边界 nmlFnmlFnmlFzyzxzszzyyxysyxzxyxsx=1)1)端面边界端面边界: :L=m=0L=m=0 ,n=1 ,n=1 00,0sxDsyDsyszxDzyDzyzxxDDsxzxsyzyszF dxdyF dxdyFdxdydxdyxy dxxF y dxddyyTFFF=0 ,0 ,0.0sxsyszzyzxFFFlm=自由边界给定面力为零2)2)横截面侧边界横截面侧边界: :L L、mm ，n=n=0 0 sxxzsyzyszxzyzFnFnFlm= 7.1 位移解法位移解法 5对于位移法求解，需要将平衡微分方程用对于位移法求解，需要将平衡微分方程用位移分量表

7、示位移分量表示位移解法基本方程 将位移表达式代入上式，将位移表达式代入上式，则则上式为上式为Laplace Laplace 方程，方程，它表示位移分量如果满足位移表示的平衡微分方程，它表示位移分量如果满足位移表示的平衡微分方程，即即LamLam方程时，则扭转翘曲函数方程时，则扭转翘曲函数F F (x，y)为调和函数。为调和函数。 022222=F=FFyx调和方程 7.1 位移解法位移解法 6对于平衡微分方程平衡微分方程，在不计体力不计体力的条件下，的条件下，前两个方程自然满足，只有最后一个方程，为前两个方程自然满足，只有最后一个方程，为0yzxzxy=,0 xyyyzzxzxGyGxxyjj

8、=FF=,( , ),yzxzyuvwxjjj= = 7.1 位移解法位移解法 7侧面边界条件侧面边界条件:因此，在横截面周界S上，有 0y lx mxyFF=xmylmylx=FF（在S上） 端面边界条件端面边界条件:考察右端面，考察右端面，l=m=0，而，而n=1。面力的合力为外力矩面力的合力为外力矩T。 7.1 位移解法位移解法 800zxDzyDzyzxDdxdydxdyxy dxdyT= 对于上述边界条件的前两式，对于上述边界条件的前两式， 同理同理 所以，在满足侧边界条件下，所以，在满足侧边界条件下， 端部边界条件的前两式是恒满足的。端部边界条件的前两式是恒满足的。 7.1 位移解

9、法位移解法 90zxDDDsxyy ldxdyGy dxdyxGdxdyxyxGxxxyxsxdmyjjjF=FFF=F0zyDdxdy=0y lx mxyFF=对于端部边界条件的第三式，有对于端部边界条件的第三式，有令令则则： T = j jGKD其中其中K KD D表达了横截面的几何特征表达了横截面的几何特征( (抗扭系数抗扭系数) )，GKGKD D 称为柱体的称为柱体的抗扭刚度抗扭刚度。 7.1 位移解法位移解法 922DDTGx xy y dxdyyxGxyxydxdyyxjjFF=FF=22DDKxyxydxdyyxFF=00zxDzyDzyzxDdxdydxdyxy dxdyT=

10、总之，柱体的自由扭转的总之，柱体的自由扭转的位移解法位移解法，归结为在边界条件，归结为在边界条件相对扭转角相对扭转角j j由公式由公式T =j j GKD确定。确定。 7.1 位移解法位移解法 9022222=F=FFyx求解方程求解方程xmylmylx=FF（在S上） 扭转问题的位移解法方程虽然简单，但是边界条件相对扭转问题的位移解法方程虽然简单，但是边界条件相对比较复杂，比较复杂，因此通常使用应力解法求解柱体的扭转问题。因此通常使用应力解法求解柱体的扭转问题。根据扭转问题的平衡微分方程，可得根据扭转问题的平衡微分方程，可得 因此，必然有一个函数因此，必然有一个函数 (x,y)，使得使得7.

11、2 扭转问题的应力解法扭转问题的应力解法yzxzxy=,xzyzyx= 将上述扭转应力分量将上述扭转应力分量代入变形协调方程，代入变形协调方程，所以，函数所以，函数 (x,y) 满足满足220,0 xzyz=220,0 xy=则前四个方程恒满足，则前四个方程恒满足，而后两个方程要求而后两个方程要求101010101010222222222222222= xxyxzyxzxxx yxz yxz x,22Gj= (),()GyGxyxxyjjFF= (x, y)普朗特（Prandtl）扭转应力函数 将扭转应力函数与翘曲函数公式相比较，将扭转应力函数与翘曲函数公式相比较，则扭转应力函数与翘曲函数的关

12、系为则扭转应力函数与翘曲函数的关系为 因此因此 2, x yC=将上式代入变形协调方程，则 C=2Gj 。 泊松方程扭转问题应力解法的基本方程 7.2 应力解法应力解法 2 7.2 应力解法应力解法 3c=const 边界条件侧面侧面端面端面单连域取为0，即：侧面侧面边界条件：0 xzyzlm=0dydxlmyxy dsx dss=但是应该注意，如果柱体横截面为多连域时，应力函数在每一个边界都是常数，但是各个常数一般并不相同。因此，只能将其中一个边界上的 c取为零。c=0 7.2 应力解法应力解法 4端面端面边界条件：XYZxzyz=,0考察右端面，考察右端面，l=m=0l=m=0，而，而n=

13、1n=1。因此面力边界条件：因此面力边界条件： 面力的合力为外力矩面力的合力为外力矩T T，则则端面面力边界条件为端面面力边界条件为： 00DDDXdxdyYdxdyYxXy dxdyT=可以证明，只要c(x,y)=C（C为常数），前二式是自然满足的，即端部应力之合力为零 =SyxTdd2222yzxzDDDcDccDTxy dxdyxy dxdyxyxydxdyxyx ly m dsdxdyxlym dsdxdy= = = = 7.2 应力解法应力解法 5下面推导应力函数, x y表示的扭矩T关系式 横截面为单连通区域， c=0扭矩方程扭矩方程 Prandtl 应力丘方程应力丘方程2V横截面

14、为多连通区域 C00,0Cx y=CCCCn123,KKKKn123,设应力函数在外边界上的值为零，即在内边界的值分别为 7.2 应力解法应力解法 60,123,122nicC C C CCDniCiDTxlym dsdxdyxlym dsdxdy= = 由Green公式计算得 1122iiiiiCCCAiAxlym dsxdyydxxdyydxdxdydxdyA= = = = 反向 122niiiDTdxdyK AV= 7.2 应力解法应力解法 7扭矩方程的这种直观的几何解释，为下一节将要介绍的借助所谓“薄膜比拟”法求解复杂截面柱杆的扭转问题提供了理论依据。横截面为多连通域的横截面为多连通域

15、的 扭矩方程扭矩方程 Prandtl 应力丘方程应力丘方程 7.2 应力解法应力解法 8横截面为单连通域的扭转问题应力解法横截面为单连通域的扭转问题应力解法在截面内要求满足泊松方程在边界上，要求满足 c(x,y)=022Gj= 22DTdxdyV= Prandtl 应力丘方程单位长度扭转角扭转应力,xzyzyx= 7.4 椭圆截面杆件扭转椭圆截面杆件扭转设有椭圆截面直杆，它的横截面为椭圆边界，椭圆的长短半轴分别为a和b，如图所示 xayb22221=根据自由扭转问题的基本方程，应力函数在横截面的边界上应该等于零，所以假设应力函数为： 22221xymab=x,y 7.4 椭圆截面杆件椭圆截面杆

16、件222Gj= 将上式代入基本方程： 222222222,mma b GGmababjj= = 则扭转基本方程满足 22222222,1a b Gxyx yababj= 将应力函数代入端面边界条件公式，则： 2222222222,1a b GxyTx y dxdydxdyababj= 7.4 椭圆截面杆件椭圆截面杆件3 在椭圆区域，有： 计算可得 回代可得应力函数表达式 2222,1Txyx yab ab= 将其代入应力分量计算公式，可以得到横截面应力分量为xbaTyabTyzxz332,2= 7.4 椭圆截面杆件椭圆截面杆件4 横截面上的任意一点的合成切应力为： 最大切应力发生在椭圆边界上，

17、边界切应力最大值在椭圆短轴处，而最小值在椭圆的长轴处，如图所示 4242222byaxabTyzxz=baTabT2min2max2,2= 7.4 椭圆截面杆件椭圆截面杆件5 椭圆截面翘曲椭圆截面翘曲 w(x,y)由于 将上面两式分别对x和y积分，则 7.4 椭圆截面杆件椭圆截面杆件6 比较上述两式，必然有比较上述两式，必然有f1(x) = f1(x)k(常数常数)，所以，所以 2222( , )ab x yxyab= 常数xybGabaTyxw3322),(=j 上式表达了横截面在变形后并不是保持为平面，而是翘曲成上式表达了横截面在变形后并不是保持为平面，而是翘曲成为曲面，成为双曲抛物面，如

18、图所示。为曲面，成为双曲抛物面，如图所示。 曲面的等高线在曲面的等高线在Oxy面上的投影是双曲线，而且这些双曲线面上的投影是双曲线，而且这些双曲线的渐近线是的渐近线是x轴和轴和y轴。轴。 只有当只有当a=b时，即圆截面杆，才有时，即圆截面杆，才有w=0，横截面保持为平面。横截面保持为平面。7.3 薄膜比拟薄膜比拟德国力学家普朗特（Prandtl） 基本思想：作用均匀压力的薄膜与柱体扭转有着相似的微分方程和边界条件 7.3 薄膜比拟薄膜比拟2 薄膜比拟法通过测试薄膜弯曲的情况，分析柱体扭转时横截面的应力分布 T22222FqZyZxZ=薄膜边界垂度 Z=0 =SyxZVdd22薄膜所围的体积调整薄膜的高度， 使2V=T，则 Z= ， 薄膜垂度Z与扭转应力函数具有相同的函数形式c=0 =SyxTdd2 7.3 薄膜比拟薄膜比拟3 这就是薄膜垂度微分方程