|
发表于 2016-12-7 15:51:11
|
1. 梁理论介绍
应用梁单元之前,首先应当看考察问题是否适合用梁单元建模。梁理论是一维近似三维,近似的前提是细长比假定,即梁的截面尺寸相比于梁轴线方向的典型尺寸足够小,所谓典型尺寸是一种整体尺寸,而非单元长度。如支座间的距离、发生显著改变的梁截面间的距离等。
梁单元是三维或二维空间中的一维线单元,具有一定的抵抗线(梁轴线)变形的刚度。这种变形包含轴向变形、弯曲变形,发生于梁轴线与横截面间的横向剪切变形,空间中还包含扭转变形。abaqus/standard中的部分梁单元还包含翘曲变形(发生在梁横截面上的不均匀的水平面变形)。梁单元的优势在于几何结构简单,并且自由度较少。它以“参考线”(梁轴线)的运动代替实际上的三维实体梁的运动,这种几何上简化的前提是假定梁的全部变形可以仅仅从“参考线”(梁轴线)位置的函数得到。应用梁单元的关键在于判断这种一维建模方法是否合适。
基本假定是垂直于梁轴线的平面,即梁截面不能在其自身平面内变形(除了梁截面面积的不断变化,这可能发生在几何非线性分析中,并且梁截面上各个方向具有相同的应变)。使用任何梁单元前都应仔细考察此假定是否成立,尤其是分析承受大量弯矩或轴向拉、压荷载的非实体截面,比如管道、i形及u形截面。这些截面可能发生崩溃(section collapse),使截面性能变的很差,不能被梁理论预测。类似地,薄壁弯曲管道的抗弯性能更弱,也不能被梁理论预测,因为管壁很容易在自身平面内弯曲,这是由于上述基本假定导致的梁理论所不能考虑的另一种效应。
在动力学分析和特征值分析中,对于细长梁结构,通常梁横截面的转动惯量影响并不显著(绕梁轴线的扭转除外)。因此,abaqus/standard忽略了欧拉-伯努力(euler-bernoulli)梁元弯曲变形时梁截面的转动惯量。对于较厚的梁,转动惯量对动力分析有一定影响,但影响程度要比剪切变形小。
对于铁木辛克(timoshenko)梁单元,程序通过横截面的几何特性计算惯性特性,与扭转模态和弯曲模态相关的转动惯量是不同的。对于非对称横截面,转动惯量在各个弯曲方向也不相同。abaqus允许用户选择铁木辛克梁元的转动惯量公式,如果选择近似的各向同性公式,在abaqus/standard中,所有的转动自由度都被赋予和扭转模态关联的转动惯量,而在abaqus/explicit中,则被赋予放大了的弯曲惯量,通过放大系数的选择使稳定时间增量最大化;梁横截面的质量中心位于节点处。当采用精确的各向异性公式时,弯曲和扭转对应的转动惯量是不同的,而且当梁截面的质量中心不在节点处时,在梁截面定义中包含平动自由度和转动自由度的耦合。使用精确的转动惯量公式(默认)时,可以定义附加质量和附加转动惯量,它们仅影响梁的惯性反应,不增加结构的刚度。
2. 梁单元的选择
abaqus中的梁单元分为欧拉-伯努力梁元和铁木辛克梁元两类,支持实心截面、薄壁闭口截面、薄壁开口截面。
abaqus/standard中的梁单元包括:
l 平面及空间的欧拉-伯努力(细长)梁; ·
l 平面及空间的铁木辛克(剪切变形)梁; ·
l 线性、二次、三次插值公式; ·
l 翘曲(开口截面)梁; ·
l 管单元;
l 杂交梁,通常用于具有明显转动的非常刚硬的梁;
abaqus/explicit中的梁单元包括:
l 平面及空间的铁木辛克(剪切柔性)梁; ·
l 线性及二次插值公式;
(1) 欧拉-伯努利梁
欧拉-伯努力梁元(b23, b23h, b33, b33h)仅在abaqus/standard中提供。不允许横向剪切变形;初始垂直于梁轴线的平截面变形后依然保持平面(如果没有翘曲),并垂直于梁轴线。只能用于模拟细长梁:梁截面尺寸相比于梁轴线方向的典型尺寸,即长细比较小。对于由均一材料构成的梁,只有当典型梁截面尺寸小于梁轴线方向典型尺寸的1/15时,横向剪切变形才可以忽略。单元不包含由压力产生的荷载刚度。
(2) 铁木辛克梁
铁木辛克梁(b21, b22, b31. b32, b31os, b32os, pipe21, pipe22, pipe31, pipe32及对应的杂交单元)允许横向剪切变形。既可分析厚梁,又可细长梁。对于由均一材料构成的梁,对于截面尺寸达到轴向典型尺寸或对于结构反应有显著贡献的最高阶振动波长的1/8的梁,剪切变形梁理论能够提供有用的结果。在这个比例以外,如果还仅仅用梁轴线位置的函数来描述构件的行为,将得不到足够准确的结果。
abaqus假定铁木辛克梁的横向剪切变形是线弹性的,具有固定的模量,因此独立于梁截面的轴向拉、压和弯曲反应。
对于大部分梁截面,abaqus会自动计算横向剪切刚度,用户也可以自定义。如果程序不能从输入部分得到剪切模量值,将无法计算缺省的剪切刚度,比如使用子程序(umat, uhypel, uhyper, vumat)定义材料的情况。对于这些情况,用户必须自己定义剪切刚度值。
铁木辛克梁元可以承受很大的轴向变形,并假定由扭转引起的轴向应变很小。在拉(压)-扭组合荷载下,只有当轴向应变不太大的时候,才能精确计算扭转引起的剪应变。关于剪切刚度的计算可参考帮助文档。 |
|