剑桥大学工程系《ABAQUS常见问题解答》翻译连载(3):单元

点亮“在看”,优质内容不错过

注:本文为系列翻译文章第3篇,如有不当之处,欢迎交流指正!尤其是复核材料部分确实不熟悉,请各位审核留言提意见。

Chapter-3:Elements

第三章 单元
Q3.1:如何查看壳单元的正向法向方向?

对于壳单元,法向正方向遵循右手定则。手指弯曲方向为单元节点按顺序编号,大拇指方向为法向,如图所示。节点顺序编号可以根据下面的关键字命令行来获得:

*ELEMENT, TYPE=S4R530, 101, 102, 202, 201

如果采用ABAQUS/CAE则直接通过单击Tools-->Query...选项,在窗口中选择Shell/Membrane NormalShell element normal,abaqus6.14版本),然后选中对应元素完成。

Colour coding :   Brown - Positive   Magenta - NegativeShell/Membrane faces are colored by normal:Brown=positive,Purple=negative(abaqus6.14版本)
如果是使用ABAQUS/Viewer对*.odb文件操作,则在Visualization模块选择Plot Material Orientation显示法向,红线表明壳单元的法向法向,其中n1法向为蓝绿色,n2法向为黄色。
Q3.2:壳单元和2D实体单元区别?

2D平面单元(如CPE4、CPE8——用于平面应变分析,CP4、CPS8——用于平面应力分析)仅适用于平面内受到载荷作用的情况,这些单元仅有平面内的自由ux、uy,这就意味着在Z方向没有变化。壳单元则可以适用于受到平面外部载荷作用的工况,这种情况下平面应力或者应变分析类型的单元则不适用。

壳单元还可以用来分析复合材料受到平面载荷作用的情况,可以通过材料属性设置壳单元厚度进行分析。
Q3.3:如何定义空间梁的n1方向?

可以直接查阅帮助文档中的Beam elements章节中的Beam element cross-section orientation获取相关内容。

当然有很多种方法定义n1方向,如对于一个矩形截面的梁单元,n1方向和宽度方向平行。

如果梁全部处于X-Y平面内,则n1方向默认为-Z方向。下面是其它一些定义n1方向的方法:

  • 1)在定义单元的数据行中指定附加的节点,假设为第三个点,那么从节点1指向节点3的方向称为v,作为n1的近似方向。abaqus则基于txv定义n2方向,有了n2就可以定义n1(根据n2xt)。总结一下:只要v处在t和n1相同的平面内,则矢量方向就不会出错。

注:t方向定义为单元的一个节点指向下一个节点。

  • 2)单元截面指定近似的n1方向,与上述方法一样,abaqus将会首先确定n2方向然后在计算n1方向。

如果既定义了附加节点又指定了近似的n1方向,则abaqus优先采用附加节点。这里有两个方法可以覆盖前文提到的由t和近似n1定义的n2方向:

  • 给出n2的矢量分量作为节点坐标的第4、第5和第6个数据(inp文件中的*NODE关键字块下);

  • 使用关键字*NORMAL直接指定。

如果上述两者包含,则优先采用关键字*NORMAL。当有了n2方向,局部梁单元的切线t将由n1xn2重新定义。

完成分析之后可以利用下面的命令来可视化t\n1\n2方向:

*SET,FILL=ON*VIEW,VIEW=(X1,Y1,Z1),UP=(X2,Y2,Z2)Example:*VIEW,VIEW=(1,1.5,2),UP=(0,1,0)*DRAW,NORMAL
n1显示为蓝色,n2为红色,t为白色。
Q3.4:对于空间梁单元为什么需要定义n1方向?
考虑非圆截面如长方形截面,梁的弯曲刚度受宽度方向的影响很大。n1方向的确定能够避免产生模棱两可的影响。
Q3.5:如何解读壳单元的应力输出值?

使用壳单元分析时,我们应该心里清楚应力是与每个单元的局部坐标相对于整体坐标方向决定的。

对于有关约定方向方面的内容,可以参考Autodes提供的一篇文章,链接地址:http://help.autodesk.com/view/RSAPRO/2020/ENU/?guid=GUID-AE34635A-F2B5-4A62-A6B0-95066594D554
Q3.6:如何定义平面复合材料壳单元的铺层方向?

每一层的局部坐标系采用关键字*ORIENTATION指定,Z方向垂直于矩形板。3层铺层方向都相对于垂向旋转来指定,如下图所示:

铺层方向为-30°的情形如图所示(基于法向的逆时针为正)

Q3.7:如何定义圆柱复材壳单元的铺层方向?

首先定义圆柱坐标系。如图所示定义了全局坐标系,Z'与全局坐标系的X方向一致,关键字命令如下:

*ORIENTATION, SYSTEM=CYLINDRICAL, NAME=CYLIND < xa >, < ya > , < za > , < xb >, < yb > , < zb > < normal-axis > , < rotation angle >

则有

*ORIENTATION, SYSTEM=CYLINDRICAL, NAME=CYLIND 0, 0, 0, 1, 0, 0 1, 0

第二行表明3个轴与圆柱表明垂直。第一个数为轴1也就是X',第二个数是相对于法向轴转动的角度,如果需要对材料的方向进行重新定义,任何转动都是按照顺时针进行的。

不作任何旋转,另外两个轴(Y'和Z')建立了局部方向。Y'和Z'各自是材料1和材料2的方向。如图所示,蓝色为材料1的方向,蓝绿色为材料2的方向。

假设复合材料在圆柱方向有三个铺层组成,方向分别为+10°,0°,-10°,可以采用下面的命令指定:

***ORIENT***SHELL SECTION, ELSET=CYLIND, COMPOSITE, ORIENTATION=ORI0.05, ,WOVEN,  100.05, , PLYTRON, 00.05, ,PLYTRON, 0
*ORIENTATION, SYSTEM=CYLINDRICAL, NAME=ORI0,0,0,0,0,11,0*****SHELL SECTION, ELSET=CYLIND, COMPOSITE, ORIENTATION=ORI0.05, ,WOVEN, -450.05, , PLYTRON, 900.05, ,PLYTRON, 0

再看一个例子,假设0°的铺层圆柱轴向重合,另外两个铺层方向为30°分布在轴向两侧。这种情况下最好使材料1的方向和轴向一致,可以通过相对法向额外旋转90°实现。

***ORIENTATION, SYSTEM=CYLINDRICAL, NAME=ORI0,0,0,0,0,11,90*****SHELL SECTION, ELSET=CYLIND, COMPOSITE, ORIENTATION=ORI0.05, ,WOVEN,  300.05, , PLYTRON, 00.05, ,PLYTRON, -30
Q3.8:2D实体分析(平面应力/平面应变)是否可以利用壳单元?
不行。不可能用一个三维壳单元到二维实体分析中。尝试在Wire草图中绘制一根线从理论上代替壳单元,然后也可以创建一个壳单元截面属性。然而这样操作无法将截面属性赋予模型,截面属性下拉选项中没有选项。可以考虑使用Beam单元创建面。
Q3.9:我想创建一个从压缩到拉伸的变刚度的弹簧单元,如何处理?

可以按照下面的命令创建弹簧单元:

*Element, type=Spring2, elset=Springs/Dashpots-1-spring1, Bay1GroundFloor.5, Part-2-1.12, Bay1GroundFloor.7, Part-2-1.43, Bay1GroundFloor.8, Part-2-1.64, Bay1GroundFloor.3, Part-2-1.2

设置弹簧刚度:

*Spring, elset=Springs/Dashpots-1-spring2, 2**** spring stiffness**7.8E+06,

通过新增关键字nonlinear可以实现弹簧从压缩到拉伸的变刚度定义:

*Spring, elset=Springs/Dashpots-1-spring, nonlinear2, 2**** Force1, Rel. Displacement-1** Force2, Rel. Displacemnet-2**10.E8, -0.110.E6, -0.001 0, 07.8E+06, 1.0
Q3.10:当我们在进行网格划分的时候如何改变单元的自由度?

在Mesh模块单击Element Type,在弹出的窗口中可以选择单元的阶次(线性或者二次)。 当使用二阶二维单元和六面体单元,如果采用低阶单元则选择Incompatible Mode

后台窗口回复关键字“单元章节”获取英文原文PDF文件!

//////////

培训、业务咨询、职位招聘等请移步公众号窗口菜单栏查看,欢迎交流合作!

(0)

相关推荐