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

1年前浏览3516

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

在Workbench中，无法直接生成Link单元（图惜注：2020之后的版本已改进），需要在模型中插入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中不需要再创建点焊连接。