Ansys Workbench工程实例之——Link杆单元静力学分析

二力杆是工程中常见的基本构件,它只能承受轴向拉伸和压缩,不能承受剪力和弯矩,杆单元可用于模拟桁架、拉索、链杆、弹簧、弹性支座、橡胶垫等。比如在现代桥梁中广泛使用了钢索、拉杆等结构,通过梁单元不能很真实第模拟它们,此时就需要通过Link(杆)单元模拟这些构件的特性。

在Workbench中,无法直接生成Link单元,需要在模型中插入Command命令(即APDL命令)来生成Link单元,常用LINK 180单元,具有塑性、蠕变、旋转、大变形和大应变功能,支持弹性、各向同性强化塑性、随动强化塑性等材料定义。Link180单元只有X、Y、Z三个方向的平移自由度,没有旋转自由度,一般都基于大变形状态。以下通过实例详解Link单元的使用。

实例1,如下模型,横杆长1m,斜杆与之夹角30°,两杆的截面积为100mm²,在两杆交汇处吊有1000N的重物,求两杆的变形与应力(不考虑压杆稳定性)。

Step1 建模,本例模型较简单,直接使用DM建模。

在DM的XY平面创建草图,绘制一条折线,水平线长1000mm,斜线与水平线夹角30°。然后通过Concept——Lines From Sketches创建线体。

通过Concept——Cross Section——Rectangular创建正方形截面,边长为10mm。给梁线赋予截面。

注意:也可创建其他形状的截面,因为杆单元只有轴向变形、轴向应力和轴力,它们与截面积有关系,而与截面形状无关系。此处创建的截面积也不是必须满足100mm²,因为Link单元的截面积是通过Command命令定义的。此处创建截面只是为了导入Mechanical后不报错且在后处理云图中更直观。

Step2 分析前处理。

(1)模型网格类型定义。梁的默认网格是Beam188,此处需要将拉索改为Link180,并且设置其他属性。

右击特征树中拉索的线体——Insert——Commands,在文本最后加入以下信息:

et,matid,link180

r,matid,100

sstif,on

其中:

et,matid,link180 表示定义当前几何体为Link180单元。

r,matid,100是r,matid,100,0,0的省略写法。表示定义当前几何体的截面积为100mm²,注意校核下图中的单位。省略的第一个0表示keyopt(2)参数,一般保持默认0或不填,表示保持体积不变(拉长时截面收缩、压缩时截面扩张),若改为1则表示截面积不变。第二个0表示keyopt(3)参数,默认0或不填则表示能承受拉伸和压缩,即二力杆状态,如果为1则表示只能承受拉伸,即拉索状态,如果为2则表示只能承受压缩。

sstif,on表示应力刚化,对于Link单元,如果是线性分析(未打开大变形开关),那么需要添加此行,如果已打开大变形,那么系统默认此选项打开。所以建议link单元都添加此句。

(2)网格划分。

Link单元通常只划分一份网格,如果涉及找形分析的话,可以划分为多个单元,但常常需要调试才能收敛。在Mesh下添加Sizing控制,在Geometry中选择2条杆线,Type为Number of Divisions,数量为1,Behavior为Hard(强制,试情况可不设置)。

(3)设置边界条件。

固定左端两端点,对右端端点施加向下1000N的力。

Step3 结果后处理。

计算后可直接添加变形结果,总变形云图如下。

但是无法通过Srees、Beam Tool添加Link应力结果,也无法通过Beam Results添加轴力结果。通过Woorksheet也无法添加以上结果。要添加轴向应力与轴力结果,需要导入到Ansys经典截面,操作如下。

设置Analysis Settings属性,Analysis Data Management下的Save MAPDL db修改为yes,表示保存为DB文件用于与APDL交互。

右击Solution——Clear Generated Data,清除计算数据,然后再右击Solution——Solve重新计算。

退出Mechanical,回到WB主界面,将Toolbox工具箱中Component Systems下的Mechanical APDL拖入到本项目的Solution栏中。

随后出现APDL项目,点击Update Project更新项目,使原项目所有栏出现√。右击Analysis栏——Edit in Mechanical APDL...,进入APDL界面。

进入APDL界面后,点击Toolbar中的RESUM_DB刷新DB数据,将在图形区出现模型。

(1)应力添加方法如下:

显示截面:显示应力前需要先将截面形状显示出来,点击菜单栏PlotCtrls——Style——Size and Shape,出现对话框,点击Display of element后方的勾选框,使其变为On,点击对话框下方ok按钮退出。

读取结果文件:点击General Postproc——Read Results——Last Set,即可读取结果文件。

显示应力云图:点击General Postproc——Plot Results——Contour Plot——Element Solu,跳出对话框。在对话框中选择Stress——X-Component of stressl或von Mises stress,再点击ok按钮,即可添加轴向(X向)应力结果,由于杆只能受轴向力,所以等效应力数值与轴向应力相同,但是X-Component of stressl可显示正负,正为拉应力,负为压应力。

轴向应力结果如下图,斜杆拉应力为20Mpa,平杆压应力为17.32Mpa。

(2)轴力添加方法如下:

求轴力前需要先定义单元表:点击General Postproc——Plot Results—Element Table——Define Table,在跳出的Element Table Data窗口中选择ADD...,跳出Define Addition Element Table Items窗口,在User Label for item后输入forcee,在Rrsults data item的左窗口中选择最下端By sequence num选项,右方窗口默认,在右下窗口的SMISC,后输入1。点击ok按钮。出现警告窗口,点击Close按钮,点击Element Table Data窗口Close按钮。

显示轴力云图:点击General Postproc——Plot Results——Contour Plot——Elem Table,跳出对话框。在对话框中选择FORCEE(默认),再点击ok按钮,即可添加轴力结果。正值表示受拉,负值表示受压。

结果显示如下,斜拉杆受拉力2000N,水平杆受压力1732N。

此处使用力学公式可验证上面的计算结果:

斜拉杆受到的拉力=1000/sin30°=2000N

水平杆受到的压力=1000/tan30°=1732N

杆的应力=受力/截面积

斜拉杆拉伸应力=2000/100=20Mpa

水平杆压缩应力=1732/100=13.32Mpa

可见与软件计算一致。

实例2,某信号塔主体由直径为Φ50mm的圆钢构成,安装后高度为10m,在距离地面8m处均布架设有4根Φ10mm,相互对称夹角为20°的钢绳拉紧在地面,每根钢绳初始预紧量为全长的1%(81.2mm)。型号塔顶端安装的设备最大受到1000N的横向风力,求信号塔主梁与拉绳的变形与应力。

Step1 建模,本例模型较简单,直接使用DM建模。

在DM的XY平面创建草图,绘制一条竖直的10000mm的线条。然后通过Concept——Lines From Sketches创建线体。

通过Create——Point在刚才的线体上创建一个加载点(Point Load),属性中Type选择Point Load,Definition选择Single(单个),Guide Edges选择线体,FD1输入8000。边在线体上创建了加载点。

再在XY面创建草图,绘制一条折现,对称角为20°,折点高8000,通过Concept——Lines From Sketches创建线体,创建时注意Operation选择Add Frozen,否则新创建的线体将和之前的合并为一个零件。

另外一条折现通过在YZ平面绘制再提取完成。也可通过Create——Body Transformation——Rotate旋转复制刚才的折现完成。

通过Concept——Cross Section——Circular分别创建两个圆形截面,Circular1半径为25mm,Circular2半径为5mm。给竖直梁线赋予截面Circular1,给拉索赋予截面Circular2。

Step2 分析前处理。

(1)模型网格类型定义。梁的默认网格是Beam188,此处需要将拉索改为Link180,并且设置单元其他属性。

右击特征树中拉索的线体——Insert——Commands,在文本最后加入以下信息:

et,matid,link180 !定义单元类型为Link180。

r,matid,78.5,,1 !截面积为78.5mm²,仅能承受拉伸。

sstif,on !刚度硬化打开。

inis,set,dtyp,epel !初始状态说明

inis,defi,,,,,0.01 !初始应变为0.01(即1%)

其中:

r,matid,78.5,,1 表示定义当前几何体的截面积为78.5mm²,保持体积不变,仅能承受拉伸,即拉索状态。

对于inis,set,dtyp,epel与inis,defi,,,,,0.01也可不定义,而使用两步载荷,第一步使用强制位移预紧,第二步加载边界条件,此处不作详解。注意在实际工程中,拉索初始应变远远达不到0.01,此例只为说明操作方法。

(2)连接的定义。

由于梁单元有6各自由度,而Link单元只有3个自由度,一般不通过Form New Part共享节点,因此本例在竖梁上建立了Point Load用于连接梁与Link。梁与Link之间通常采用Spot Weld(点焊)连接。

右击特征树Connections——Insert——Spot Weld,分别选择两条Link上的折点与梁上的点,如果不好选择可先将其余线体隐藏(F9)。最终使两条Link的拐点都与梁上的Point Load连接。

(3)网格划分。

每条Link线体仅划分一个单元,在Mesh下添加Sizing控制,在Geometry中选择四条拉索,Type为Number of Divisions,数量为1,Behavior为Hard(强制)。

同样,给竖梁划分500份网格。

(4)设置边界条件。

在有预紧的静力学分析中,一般需要将计算分为两步,即对Analysis Settings设置2步,一步预紧,一步加载工况。比如螺栓连接、拉索连接等。

但是在本例中,由于在Link中已经设置了初始预紧应变,不需要再对4条拉索设置力或强制位移,所以也不需要分步。

但是在Analysis Settings属性中需要打开大变形选项。

固定竖梁和4条拉索下顶点,对上顶点施加1000N沿X方向的力。

Step3 结果后处理。

总变形如下,梁与Link单元的变形云图都可以通过Mechanical显示。

添加梁应力结果:选择特征树Solution,点击工具栏Tool——Beam Tool,可在结果中添加梁结果,最大组合应力云图如下左。

添加等效应力结果:选择特征树Solution,在属性中将Beam Section Results改为Yes,即可添加工具栏Stress——Equivalent (von-Mises),云图如下右。

对于Link单元的应力结果需要到APDL界面查看,操作同实例1,此例为了使结果显示清晰,我们需要隐藏网格轮廓,点击标题栏PlotCtrls——Style——Edge Options——Element outline style选择none。等效应力结果如下图。

本例扩展:

本例拉索与竖直梁连接处采用了点焊连接,也可以采用Form New Part共享节点,此时需要在竖直梁的草绘图中作两段线,下段长度8m,上端长度2m。提取草图之后不需要在竖直梁上添加Point Load点。作折现草图时注意折点处需要出现P标志才表示与竖直梁的中间端点重合了。最后在通过Form New Part组合三条线体即可,其余操作同上文。在Mechanical中不需要再创建点焊连接。

以上两例中Link单元均只划分了1个单元,杆索中间不受垂直力和弯矩,那么Link180是否可以模拟由于受重力而产生向下的挠度的高空电缆呢?

实例3,某电缆长10m,截面直径20mm,固定两端,初始应变0.0001,求最大挠度。

Step1 建模,本例模型较简单,直接使用DM建模。

在DM的XY平面创建草图,绘制一条水平直线长10000mm,然后通过Concept——Lines From Sketches创建线体。

通过Concept——Cross Section——Rectangular创建正方形截面,边长为10mm。给梁线赋予截面。

通过Concept——Cross Section——Circular分别创建半径为10mm的圆形截面,给电缆赋予截面Circular1。

Step2 分析前处理。

(1)模型网格类型定义。右击特征树中拉索的线体——Insert——Commands,在文本最后加入以下信息:

et,matid,link180 !定义单元类型为Link180。

r,matid,1256,,1 !截面积为1256mm²,仅能承受拉伸。

sstif,on !刚度硬化打开。

inis,set,dtyp,epel !初始状态说明。

inis,defi,,,,,0.0001 !初始应变为0.0001。

(2)网格划分。

Link线体划分100个单元:在Mesh下添加Sizing控制,在Geometry中选择四条拉索,Type为Number of Divisions,数量为100,Behavior为Hard(强制)。

(3)设置边界条件。

Analysis Settings属性中需要打开大变形选项。

固定电缆两端点,对电缆施加-Y方向的247N的均布力(模拟重力作用)。

Step3 结果后处理。

总变形如下。

其实线缆弯曲并不是因为每个Link单元产生弯曲,而是单元之间的节点产生转动,如果将单元划分为2~4个,则弯曲云图如下。由此也可知,Link单元之间相当于是通过铰点连接在一起,铰点不限制单元的转动。所以有时也形象第称Link单元为链单元。当单元数量足够多时,就可以通过多段直线模拟电缆的真实弯曲。

扩展:

除了通过r,定义截面积为定常数外,也可以让软件自动识别我们模型的截面积,命令为“*get,用户定义的名称,secp,matid,pro,area”,然后通过SECTYPE与SECDATA调用用户定义的名称。比如实例2的指令可按下图修改:

写在最后:本文案例仅为说明操作方法,不能直接作为工程计算的结论。关于Link单元笔者平时很少使用,也有很多问题没有搞懂,本文难免有很多纰漏,敬请批评指正。

1.《[ansys如何中单元]ansys如何选择单元类型》援引自互联网,旨在传递更多网络信息知识,仅代表作者本人观点,与本网站无关,侵删请联系页脚下方联系方式。

2.《[ansys如何中单元]ansys如何选择单元类型》仅供读者参考,本网站未对该内容进行证实,对其原创性、真实性、完整性、及时性不作任何保证。

3.文章转载时请保留本站内容来源地址,https://www.lu-xu.com/keji/3218667.html