文档介绍：用 ABAQUS 分析扭转问题发表时间: 2005-6-15 来源:有限元在线关键字: CAE ABAQUS 信息化应用调查在线投稿加入收藏发表评论好文推荐打印文本工程中承受扭转载荷的构件有很多, 如车床的光杆、搅拌机和主轴、汽车传动轴、钻机和钻杆等。本文以 ABAQUS 软件为平台, 旨在解决两个问题, 1 、如何分析橡胶材料的性能, 2 、如何正确的完成扭转分析。 1 、引言在处理扭转问题时, 常规的计算方法, 往往会伴随一些假设, 这会降低了结果的准确程度。根据有限元理论, 使用有限元软件求解扭转问题会大大提高求解的精确度,特别是对复杂的结构,效果更为明显。本文以橡胶产品为例,讨论的在 ABAQU S 软件中,如何正确完成扭转分析,并提取需要的分析结果。 2 、问题描述受扭转件结构由钢筒和橡胶筒组成,产品尺寸如图 1 所示。图1 产品结构简图 3 、有限元建模加载时内芯固定,在外圈施加扭转位移。根据产品的 CAD 结构建立有限元模型如图 2 所示: 图2 产品有限元模型图 4 、材料性质定义钢:弹性模量 EX=2 × 105MPa ,泊松比μ= 橡胶: 橡胶是一种超弹性材料, 对于超弹性材料, 不用杨氏模量和泊松比, 而用应变势能(U) 来表达应力—应变关系。 ABAQU S 软件中有两种应变势能可利用,分别是多项式模型和奥根( Ogden )模型,本例中使用多项式模型,表达式如下: 式中: U—应变势能, Jel —弹性体积比; I1、 I2—应变不变量; Di—定义材料的压缩性; Cij — Rinvlin 系数。本例中取 N=1 , 以橡胶材料的单轴拉伸,单轴压缩和平面剪切实验数据为依据,并考虑到橡胶的不可压缩性,输入方程系数值: C01= , C10= , D1=0 注意事项: 橡胶的特性错综复杂, 材料特性和几何特性均呈非线性变化的。如果要准确预测模型中发生变形或应变部分的行为,那么提供的试验数据的范围要涵盖计算模型中可能会出现的变形状态和应变范围。 5 、加载求解加载时,内钢筒的内套固定,即 UX=UY=UZ=0, 将外钢筒的最外层结点的坐标系定义为柱坐标系。在此柱坐标系中施加扭转载荷。注意事项: 在分析过程中最外层结点的柱坐标系的方向始终保持与结点原始柱坐标系的方向相同, 并不随着结点的移动而改变方向。因此,施加的位移量并不是结点始末位置间的弧线长度,而应是结点始末位置间的直线距离。所以此例中所加的位移载荷应为: UR=60*sin30 ° =30(mm) UF=60*cos30 ° -60=-(mm) 6 、结果分析将模型信息提交给求解器,经迭代求解出结果:橡胶部分的 Von — Mises 应力云图如图 3 所示: 图3 橡胶的 Von Mises 应力云图由上图看出, 橡胶上的应力呈层状分布, 且随着半径的增加承递减趋势。其递减情况如图 4 所示。应力的这种随半径递减的分布趋势,与实际使用中,产品率先在中间发生破坏的情况是相吻合的。由分析求得的反作用力结果如下表所示计算扭矩需要沿圆周切线方向的分量,而这里的 F 合是沿弦线的合力,其沿切线方向的分量应为,F切=F合*cos15 °。计算的扭矩结果与实验结果的对照情况如下表所示产生误差的原因大致如下: 实测操作的本身;计算过程