对两端受扭矩T作用矩形截面杆探讨.docVIP

对两端受扭矩T作用矩形截面杆探讨 　　摘要: 对两端受扭矩T作用的矩形截面杆进行了ANSYS有限元分析,对中间截面上的最大剪应力和扭矩之间的关系以及该关系式随高度h的变化规律进行了研究分析,并验证平截面假设在此情形下不成立. 　　关键词：矩形截面杆;扭转;ANSYS;有限元分析;扭矩;最大剪应力;平截面假设 　　Abstract:The rectangular section bar ends by the torque of T function by ANSYS finite element analysis, the middle section of the maximum shear stress and torque as well as the relationship between the variation with height of H are analyzed, and verified the assumption of plane section in this case is not established. 　　Keywords:Rectangular section bar; torsion; ANSYS; finite element analysis; torque; the maximum shear stress plane section assumption; 　　中图分类号：O343.2文献标识码:A 　　 　　非圆截面杆受扭转时横截面将成为曲面,产生所谓翘曲现象,如图1所示矩形截面杆的扭转.所以对于非圆截面杆,平面假设不再成立,根据平面假设所建立的扭转应力公式显然不再适用.非圆截面杆的扭转问题只能用弹性力学的方法去研究. 　　 　　图1 矩形截面杆的扭转 　　为了更好地反映矩形截面杆扭转的规律,本文通过对两端受扭矩T作用,高度为h, 宽度为b,长度为a的矩形截面杆进行ANSYS有限元分析,得出了中间截面上的最大剪应力和扭矩之间的关系以及该关系式随高度h的变化规律,并验证平截面假设在此情形下不成立. 　　1 弹性力学方法分析 　　 　　图2 最大剪应力分布 　　根据弹性力学方法分析的结论, 矩形截面杆在扭矩T的作用下,切应力的方向为:周边处的切应力与周边平行;对称轴处的切应力与对称轴垂直.最大剪应力发生在截面长边的中点处(如图2),其精确计算公式为: 　　(1) 　　其中,T为扭矩,h为截面高度,b为截面宽度. 令 　　(2) 　　则式（1）可写成 　　 (3) 　　表1由电算程序计算得到的不同的h/b对应的γ值 　　 　　 　　因子γ只与比值h/b有关.通过自编的C++电算程序计算,因子γ与h/b的关系如表1. 　　该表数据与《弹性力学》教材所提供的相关数据符合. 为便于有限元分析建模,这里选用h/b值分别为1.0、1.2、1.4、1.6、2.0、3.0时的模型. 　　2 有限元分析 　　2.1 有限元模型 　　在有限元模型分析中,我们选用了kN—m的单位制.模型长度a为0.5 m,截面宽度b为0.1 m,截面高度h随h/b值的变化分别取为0.1 m,0.12 m,0.14 m,0.16 m,0.2 m,0.3 m. 　　2.2 网格划分 　　选用Structural solid 8 nodes 45单元,单元长度取为0.01 m,采用映射式划分,划分后的网格由规则的正方体组成. 　　2.3 材料属性 　　选用Q235钢,其弹性模量E为210 GPa,泊松比取为0.3. 　　2.4 边界条件 　　 在求解这个问题过程中,对矩形截面杆的两端截面中心点处施加了位移=0的约束,使其不发生刚体位移而绕中轴扭转.一方面,只固定了两个点(理论上该两点位移为零),不影响结果分析.另一方面, 对实心的矩形等截面直杆,约束扭转引起的正应力通常很小,可忽略不计,更何况其约束远未及约束扭转严格.故仍可作为自由扭转问题看待. 　　2.5 载荷施加 　　 　　图3 等效扭矩施加 　　两端的扭矩通过在两端截面边缘添加等效节点力施加.例如,在截面尺寸为0.1 m×0.1 m的情况下,划分网格后每边有11个节点,在每个节点上施加1 kN的节点力,则等效扭矩为2.2 KNm.根据圣维南原理,在离两端截面足够远处的中截面,分析结果将与直接加扭矩的情形相同. 为研究中间截面上的最大剪应力和扭矩之间的关系,以便于通过线性拟合确定两者间关系,等效扭矩取5级,分别为2.2 KNm,3.3 KNm,4.4 KNm,5.5 KNm,6.6 KNm.图3为截面尺寸为0.1 m×0.2 m模型等效扭矩的施加示意图. 2.6 计算结果 　　 　　图4