在Abaqus中使用梁单元进行计算.docx
- 文档编号:23354940
- 上传时间:2023-05-16
- 格式:DOCX
- 页数:34
- 大小:42.98KB
在Abaqus中使用梁单元进行计算.docx
《在Abaqus中使用梁单元进行计算.docx》由会员分享,可在线阅读,更多相关《在Abaqus中使用梁单元进行计算.docx(34页珍藏版)》请在冰豆网上搜索。
在Abaqus中使用梁单元进行计算
在Abaqus中使用梁单元进行计算
在Abaqus中使用梁单元进行计算
(2012-03-2611:
28:
00)
转载▼
标签:
abaqus
梁
杂谈
分类:
ABAQUS
xiaozity助理工程师:
在练习老庄的Crane例题时,欲提取梁元的截面应力。
反复折腾后,小小体会,总结如下:
(1)书中讲到:
“线性梁元B21、B31及二次梁元B22、B32是考虑剪切变形的Timoshenko梁单元;而三次梁元B23、B33不能模拟剪切变形,属Euler梁单元”。
(2)众所周知,当要考虑剪切变形时,例如深梁,采用Timoshenko梁单元比较合适。
三次梁元由于可模拟轴线方向的三阶变量,因而对static问题,一个构件常常用一个三次单元就足够,特别对于分布载荷的梁,三次梁元的精度相当高。
(3)Abaqus会默认在积分点处的若干截面点输入应力值;但用户可自定义应力输出的截面点位置,这通过property-section-manage-edit-outputpoints来定义输出应力值的截面点;
(4)特别要指出的是,无论B22还是B33还是其它梁元,其输出的应力分量只有S11,如图所示;那么,现在的问题是:
1:
S11代表什么应力,根据经验,大家会认为11是1方向的正应力或主应力等等
2:
为什么没有S22、S33、S12......
下面分别说明:
1:
S11表达的是梁元的弯曲应力,即局部坐标系下截面上的正应力
2:
只输出S11,而无其它应力,这是因为梁元之所以成为梁元,有一基本前提就是用梁元来模拟的构件,其正应力是最主要的,而剪应力是可忽略的;一个基本的佐证就是:
众所周知,在建立梁的总势能方程时,总是讲剪切应变能是小量,因而它总是被忽略掉的;忽略剪应力的一个结果是:
mises应力将与S11在数值上完全相同,不仅Abaqus如此,Ansys也是如此,这也难怪有人讲:
“Timoshenko梁单元是骗人的,它根本没有考虑剪应力”;对这件事情,我想作如下评价:
(A)不仅Timoshenko梁单元,其它梁元(不考虑剪切变形)确实在应力的层面没有考虑剪应力的影响,这可从mises应力与S11的比较看出来;而为什么这样处理,理由如上所述,剪应力是高阶量,可忽略,否则就认为不能用梁元来模拟。
(B)尽管应力层面没有考虑剪应力,但在变形层面Timoshenko梁单元还是考虑了剪切的影响;所以说Timoshenko梁单元是骗人的说法,本人不太同意。
怪兽:
助理工程师
楼上的问题有些是重复的,但为了清晰,还是独立回答
(1)变形层面考虑剪切具体有两个方面:
1:
在单元刚度中考虑了剪切刚度;
2:
从所周知,梁元的变形是由结点位移插值获得;考虑剪切变形的位移的影响就体现在插值过程中考虑剪切变形;比如Timoshenke梁元采用弯曲和剪切独立插值的形状函数;
所谓应力层面不考虑剪切,是因为:
1:
Abaqus/Ansys在计算Mises应力时只考虑了正应力,而没有把剪应力拿进去,所以它的Mises应力并非真正的等效应力,只是正应力取绝对值而已;
2:
另外,说程序在应力层面不考虑剪切,也在于它根本就没有去计算剪应力,也因此而没有剪应力的输出选项;
为了更清楚的回答你的问题,有必要提一下有限元的过程:
结点位移?
插值获得单元位形?
几何方程导出应变?
本构方程获得应力;应力是在结果的最未端,只是程序认为剪应力太不重要,在与应力有关的结果中,没有去理会它而已。
(2)如果第一个问题我说清楚了,我想第二个就不是问题了
(3)所谓考虑剪切变形,具体的就在于:
在梁元的单刚方程和插值函数中考虑了剪切变形的影响;影响谁呢,影响着原本不考虑剪切刚度的弯曲刚度和弯曲变形;也就是说,考虑剪切的影响,其实质在于:
在弯曲刚度中把剪切刚度加进去,以使梁更刚硬一些,以模拟剪切变形不能忽略的“深梁”。
(4)4\5是一个问题,一并回答;具体如何证明,我考虑单独在Abaqus中用算例证明好象行不通;
要想了解特定的梁元是否考虑剪切变形,根据我现在的经验,只能看它的理论手册,一看单元刚度方程,二看插值函数。
如果假定程序给出了插值阶数相同的梁元A、B,其中:
A?
考虑剪切
B?
不考虑剪切
那么,在其它条件相同的情况下,A算出来的位移应该有小一些,因为它更刚硬;
如果一定要通过上述思路进行算例证明,考虑用ansys的梁元与Abaqus梁元进行上述比较,前提是选取的梁元插值阶数要相同。
不知道我说清楚没有~欢迎讨论
megvin助理工程师:
我猜,wong兄所提弯曲刚度提高了,可能是基于插值阶数降低之后,弯曲变得刚硬.这没问题;
我强调的是,比较的前提是:
选取相同的插值阶段;更具体一点,是选取相同的弯曲插值阶次;只有在这个条件下,才能比较两者的差别.
选取相同的弯曲插值阶次的目的在于:
固定其弯曲刚度
chsxiaolin助理工程师:
按照上面的说法,与其说“应力层面不考虑剪应力,变形层面考虑剪切变形”,还不如说“只考虑弹性剪切,不考虑塑性剪切及其影响”。
单说不考虑剪应力又哪里来的剪切刚度呢?
针对楼上提法,稍稍再补充一下:
1:
Abaqus/Ansys(其它程序不清楚)中,考虑剪切的影响,一律认为剪切永远是弹性的;如果说程序一定要计算剪应力,其本构也一定是基于弹性本构确定。
为什么认为剪切不会进入塑性,可能是因为:
如前所述,还是基于梁元的特点,认为剪切一般只会处于弹性;而且对梁元而言,剪切变形总是弯曲的小量,能考虑一下剪切影响已经很不错了~
2:
剪切刚度是构件的抗剪切变形的能力,是刚度层次;剪应力则属于是截面层次;也就是说:
在变形的时候考虑剪切刚度的影响,但在计算应力的时候,并没有去理会到底产生了多大的剪应力。
From:
http:
//58.213.153.47/viewthread.php?
tid=1150914
合理选用梁单元——BeamElementsinABAQUS
作者:
麦田
本文内容主要译自ABAQUSDocuments
1.梁理论
1.1应用梁单元之前,首先应当考察问题是否适合用梁单元建模。
梁理论是用一维近似三维,近似的前提是长细比假定,即梁的截面尺寸相比于梁轴线方向的典型尺寸足够小。
所谓典型尺寸是一种整体尺寸,而非单元长度,比如:
∙支座间的距离;
∙发生显著改变的梁截面间的距离;
∙感兴趣的最高阶振动模态的波长;
梁单元是三维或二维空间中的一维线单元,具有一定的抵抗线(梁轴线)变形的刚度。
这种变形包含轴向变形、弯曲变形,发生于梁轴线与横截面间的横向剪切变形,空间中还包含扭转变形。
ABAQUS/Standard中的部分梁单元还包含翘曲变形(发生在梁横截面上的不均匀的出平面变形)。
梁单元的主要优势在于几何构造简单,并且自由度较少。
它以“参考线”(梁轴线)的运动代替实际上的三维实体梁的运动,这种几何上简化的前提是假定梁的全部变形可以仅仅从“参考线”(梁轴线)位置的函数得到。
应用梁单元的关键在于判断这种一维建模方法是否合适。
基本假定是垂直于梁轴线的平面,即梁截面不能在其在身平面内变形(除了梁截面面积的不断变化,这可能发生在几何非线性分析中,并且梁截面上各个方向具有相同的应变)。
使用任何梁单元前都应仔细考察此假定是否成立,尤其是分析承受大量弯矩或轴向拉、压荷载的非实体截面,比如管道、I形及U形截面。
这些截面可能发生崩溃(sectioncollapse),使截面性能变的很差,不能被梁理论预测。
类似地,薄壁弯曲管道的抗弯性能更弱,也不能被梁理论预测,因为管壁很容易在自身平面内弯曲,这是由于上述基本假定导致的梁理论所不能考虑的另一种效应。
1.2在动力分析和特征值分析中应用梁单元
对于细长梁结构,通常梁横截面的转动惯量影响并不显著(绕梁轴线的扭转除外)。
因此,ABAQUS/Standard忽略了欧拉-伯努力(Euler-Bernoulli)梁元弯曲变形时梁截面的转动惯量。
对于较厚的梁,转动惯量对动力分析有一定影响,但影响程度要比剪切变形小。
对于铁木辛克(Timoshenko)梁元,程序通过横截面的几何特性计算惯性特性,与扭转模态和弯曲模态相关的转动惯量是不同的。
对于非对称横截面,转动惯量在各个弯曲方向也不相同。
ABAQUS允许用户选择铁木辛克梁元的转动惯量公式,如果选择近似的各向同性公式,在ABAQUS/Standard中,所有的转动自由度都被赋予和扭转模态关联的转动惯量,而在ABAQUS/Explicit中,则被赋予放大了的弯曲惯量,通过放大系数的选择使稳定时间增量最大化;梁横截面的质量中心位于节点处。
当采用精确的各向异性公式时,弯曲和扭转对应的转动惯量是不同的,而且当梁截面的质量中心不在节点处时,在梁截面定义中包含平动自由度和转动自由度的耦合。
使用精确的转动惯量公式(默认)时,可以定义附加质量和附加转动惯量,它们仅影响梁的惯性反应,不增加结构的刚度。
2.梁单元的选择
ABAQUS中的梁单元分为欧拉-伯努力梁元和铁木辛克梁元两类,支持实心截面、薄壁闭口截面、薄壁开口截面。
ABAQUS/Standard中的梁单元包括:
∙平面及空间的欧拉-伯努力(细长)梁;
∙平面及空间的铁木辛克(剪切变形)梁;
∙线性、二次、三次插值公式;
∙翘曲(开口截面)梁;
∙管单元
∙杂交梁,通常用于具有明显转动的非常刚硬的梁;
ABAQUS/Explicit中的梁单元包括:
∙平面及空间的铁木辛克(剪切柔性)梁;
∙线性及二次插值公式;
2.1欧拉-伯努力梁
欧拉-伯努力梁元(B23,B23H,B33,B33H)仅在ABAQUS/Standard中提供。
不允许横向剪切变形;初始垂直于梁轴线的平截面变形后依然保持平面(如果没有翘曲),并垂直于梁轴线。
只能用于模拟细长梁:
梁截面尺寸相比于梁轴线方向的典型尺寸,即长细比较小。
对于由均一材料构成的梁,只有当典型梁截面尺寸小于梁轴线方向典型尺寸的1/15时,横向剪切变形才可以忽略。
单元不包含由压力产生的荷载刚度。
插值:
欧拉-伯努力梁元采用三次插值公式,这对于梁上分布荷载是精确的。
因此,适合于动力振动分析——因为d’Alembert荷载(达朗贝尔荷载,即惯性力)也是分布的。
三次梁单元适合于分析小应变、大转动问题。
由于采用近似的方程,不适合分析扭转稳定问题,也不适合分析有特别大转动的问题(比如约180度),此时,应该采用一次或二次梁单元。
质量公式:
欧拉-伯努力梁元采用一致质量公式。
绕梁轴线扭转时的转动惯量和铁木辛克梁元一样。
不能定义附加惯量。
2.2铁木辛克梁
铁木辛克梁(B21,B22,B31.B32,B31OS,B32OS,PIPE21,PIPE22,PIPE31,PIPE32及对应的杂交单元)允许横向剪切变形。
既可分析厚梁,又可细长梁。
对于由均一材料构成的梁,对于截面尺寸达到轴向典型尺寸或对于结构反应有显著贡献的最高阶振动波长的1/8的梁,剪切变形梁理论能够提供有用的结果。
在这个比例以外,如果还仅仅用梁轴线位置的函数来描述构件的行为,将得不到足够准确的结果。
ABAQUS假定铁木辛克梁的横向剪切变形是线弹性的,具有固定的模量,因此独立于梁截面的轴向拉、压和弯曲反应。
对于大部分梁截面,ABAQUS会自动计算横向剪切刚度,用户也可以自定义。
如果程序不能从输入部分得到剪切模量值,将无法计算缺省的剪切刚度,比如使用子程序(UMAT,UHYPEL,UHYPER,VUMAT)定义材料的情况。
对于这些情况,用户必须自己定义剪切刚度值。
铁木辛克梁元可以承受很大的轴向变形,并假定由扭转引起的轴向应变很小。
在拉(压)-扭组合荷载下,只有当轴向应变不太大的时候,才能精确计算扭转引起的剪应变。
关于剪切刚度的计算可参考帮助文档。
插值:
程序为有限轴向应变、剪切变形梁单元提供一次和二次插值。
单元B21,B31,B31OS,PIPE21,PIPE31及对应的杂交单元采用线性插值。
特别适合于分析接触问题,比如沟渠或海底管线的铺设,钻头和井孔的接触,及它们对应的动力学问题。
单元B22,B32,B32OS,PIPE22,PIPE32及对应的杂交单元采用二次插值。
质量公式:
线性铁木辛克梁元采用集中质量公式。
ABAQUS/Standard中,二次铁木辛克梁元采用一致质量公式,但动力分析中将按照1:
4:
1的分布,采用集中质量公式。
ABAQUS/Explicit中的二次铁木辛克梁元也采用上述形式的集中质量公式。
转动惯量和附加惯量:
缺省情况下,铁木辛克梁采用精确的(各向异性,位移-转角耦合)转动惯量公式,也可以采用各向同性的、解耦的近似方法。
例外是,对于静力分析中采用自动稳定控制(automaticstabilization)的情况,程序在计算铁木辛克梁的质量矩阵时,直接假定转动惯量为各向同性,忽略指定的转动惯量类型。
在一些结构应用中,对于具有复杂几何形状和质量分布的截面,梁单元是一种一维的近似。
沿梁长度方向,梁截面上可能存在着对结构刚度没有贡献,但却对惯性有贡献的分布质量,如机器、船舱中的货物、装满流体的容器等。
这种情况下,可以在梁截面定义中定义附加质量和附加转动惯量。
还可以定义与附加惯量相关的质量比例阻尼。
程序根据质量比例,通过对材料阻尼和附加惯性阻尼进行加权平均,确定单元的质量比例阻尼。
翘曲(开口)梁:
在三维空间中使用梁单元时,必须注意梁截面上可能存在的由于扭矩引起的翘曲变形,除了圆形截面外,所有的截面在承受扭矩时都会发生出平面变形。
翘曲将改变截面上的剪应变分布。
如果翘曲没有得到有效约束,开口截面很容易发生扭曲,尤其是梁截面壁厚较薄的情况。
单元B31OS,B32OS以及对应的杂交单元通过在每个节点处增加一个自由度来考虑翘曲,并假定翘曲变形是截面上位置的函数,仅翘曲幅度随截面在梁轴线上的位置而改变,并且上述单元只能在ABAQUS/Standard中使用。
它们可用于分析薄壁、开口截面,对于这些截面,翘曲约束起一定作用,由翘曲引起的轴向变形不可忽略,例如I形截面和任意开口截面。
在其他的梁单元中,翘曲被认为是自由的,翘曲引起的任何轴向应力都被忽略。
用这些单元分析薄壁、开口截面梁时,扭转不能被很好的反映。
通常,只有当通过节点的梁轴线是连续的,并且节点两侧的梁横截面相同时,翘曲幅度才是连续的。
因此,如果开口梁在节点处相交,那么就需要为这个具有不同轴线方向的交叉梁定义单独的节点,并且需要为两个相交构件在这个节点处翘曲幅度定义合理的约束。
例如,如果接点被加强,翘曲变形就可能被消除;那么,在接点连接处,作为边界条件,翘曲自由度7应该在合适的构件上被完全约束。
参考资料:
ABAQUS,Inc.2008.AbaqusTheoryManual,Version6.8.Providence,RhodeIsland:
DassaultSystèmes.
From:
在进行梁单元截面内力的后处理时,会遇到SF1、SF2、SF3、SM1、SM2、SM3,其具体含义是什么呢?
根据23.3.8Beamelementlibrary的解释
Sectionforces,moments,andtransverseshearforces
SF1Axialforce.——梁单元轴力
SF2Transverseshearforceinthelocal2-direction(notavailableforB23,B23H,B33,B33H).——n2方向的横向剪力
SF3Transverseshearforceinthelocal1-direction(availableonlyforbeamsinspace,notavailableforB33,B33H).——n1方向的横向剪力
SM1Bendingmomentaboutthelocal1-axis.——绕n1轴的弯矩
SM2Bendingmomentaboutthelocal2-axis(availableonlyforbeamsinspace).——绕n2轴的弯矩
SM3Twistingmomentaboutthebeamaxis(availableonlyforbeamsinspace).——绕梁单元轴线的扭矩
ABAQUS的任务管理命令可以暂停、恢复、和终止一个正在背景运行的任务,方法如下(在命令行输入并运行):
1、暂停一个正在运行的任务:
abaqussuspendjob=job-name
2、恢复一个暂停的任务:
abaqusresumejob=job-name
3、终止一个正在运行的任务:
abaqusterminatejob=job-name
其中任务暂停(suspend)的时候,windows任务管理栏中仍会保留standard/explicit的计算线程,只是不再使用CPU资源,当任务恢复(resume)的时候继续工作。
任务终止则就像CAE中提交的任务的KILL功能类似,直接cut掉正在运行的任务,不可恢复。
其实也就和在windows任务管理栏中强行终止差不多,但属于合法操作。
ABAQUS中USERTIME/SYSTEMTIME/TOTALCPUTIME/WALLCLOCK的意义
2011-02-2421:
56:
29| 分类:
ABAQUS| 标签:
time cpu total clock wall |字号 订阅
USERTIMEreferstotheCPUtimespentexecutingABAQUS.
SYSTEMTIMEreferstotheamountofOSkernelCPUtimespentbytheoperatingsystemdoingworkonbehalfoftheABAQUSprocess.
TOTALCPUTIMEisthesumofthesetwonumbers.
WALLCLOCKtimereferstotheactualphysicaltimespentfortheanalysisprocesstocomplete.IftheanalysisjobisrunningonasingleCPU,andthejobhasexclusiveaccesstothatCPU,thedifferencebetweenTOTALCPUTIMEandWALLCLOCKTIMEislargelythetimetakentoperformallI/Orequests.
求解器算法比较应当用TotalCPU,对用户感到的时间用Wallclock
在进行壳单元截面内力的后处理时,会遇到SF1、SF2、SF3、SF4、SF5、SF6、SM1、SM2、SM3,其具体含义是什么呢?
根据23.6.8Continuumshellelementlibrary的解释
Sectionforces,moments,andtransverseshearforces
SF1Directmembraneforceperunitwidthinlocal1-direction.——n1方向单位宽度内的正拉(压)力
SF2Directmembraneforceperunitwidthinlocal2-direction.——n2方向单位宽度内的正拉(压)力
SF3Shearmembraneforceperunitwidthinlocal1–2plane.——n1-n2平面内单位宽度内的剪力
SF4Transverseshearforceperunitwidthinlocal1-direction.——n1方向单位宽度内的横向剪力
SF5Transverseshearforceperunitwidthinlocal2-direction.——n2方向单位宽度内的横向剪力
SF6Thicknessstressintegratedovertheelementthickness.——壳单元厚度方向的集中力
SM1Bendingmomentforceperunitwidthaboutlocal2-axis.——绕n2轴的弯矩
SM2Bendingmomentforceperunitwidthaboutlocal1-axis.——绕n1轴的弯矩
SM3Twistingmomentforceperunitwidthinlocal1–2plane.——n1-n2平面内的扭矩
Fileextensions
abq
Statefile,onlyusedbyAbaqus/Explicit.Itiswrittenbytheanalysis,continue,andrecoveroptions.Itisreadbytheconvertandrecoveroptions.Thisfileisrequiredforrestart.
axi
Symmetricmodeldatafile,onlyusedbyAbaqus/Standard.Itiswrittenduringsymmetricmodelgenerationbythedatacheckandanalysisoptions.
bsp
Textfilecontainingbeamcross-sectionpropertiesformeshedsectionprofiles.ItiswrittenbyAbaqus/Standardduringmeshedbeamcross-sectiongeneration.
cid
Auto-releasefile,whichcontainsinformationneededforlicenserecoveryandsuspension.
com
Commandfile,createdbytheAbaqusexecutionprocedure.
dat
Printedoutputfile.Itiswrittenbytheanalysis,datacheck,parametercheck,andcontinueoptions.Abaqus/Explicitdoesnotwriteanalysisresultstothisfile.
eig
Lanczoseigenvectorfile.ThisisatemporaryscratchfilethatisusedtostoretheeigenvectorscalculatedbytheLanczoseigensolverduringthesolutionprocedure.
f
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- Abaqus 使用 单元 进行 计算
![提示](https://static.bdocx.com/images/bang_tan.gif)