无初位移的ansys实现方法.docVIP

无初位移的ansys实现方法 姚再兴 （1北京力算科技有限公司 北京 100011） 摘要：以一个算例说明了ansys实现初始作用下无初位移、无初应变的方法和过程。对比了重力作用、无初位移、只有初应力三种工况的计算结果。计算结果从竖向位移、竖向应变、竖向应力三个方面表示。并含有全部命令流代码。 关键词：无初位移；无初应变；初应力；ansys 引言 很多计算对象在初始状态下是有力的作用的,从而在其内部的应力并不为零.但这一状态下的位移,认为是处处为零,应变也是处处为零.其后的受力和变形都是相对于这一状态而言的.例如,要计算桥梁在某重型卡车作用下的变形.在卡车作用在桥上之前,桥体已经受到重力的作用,内部有应力存在,相对于不受力的状态,存在位移和应变.但我们通常要保留桥体的应力状态,而要去掉重力引起的位移和变形.ansys的无初位移功能可以很好地实现这类计算。本文通过一个算例说明ansys的这一功能和需要注意的事项。 本文中用到的算例是在ansys14.0下通过的。 实现思路 对这一类问题，ansys首先在初始作用力下进行计算，然后把初始作用力产生的应力导出到文件里备用，其后在所有载荷（包括初始作用力）和初始应力（从先前产生的文件中导入）共同作用下进行计算。这样，初始作用力下的位移和应变消除了，初始作用下的应力得到保留。 算例 问题描述 边长为1m的立方体平放在水平面上，底部完全约束，初始作用力是重力，求在没有其它作用力情况下的位移、应变和应力，如图 1。材料参数如下： 弹性模量300E9Pa 泊松比0.3 密度2000kg/m3 重力加速度10m/s2 图 1计算无初位移的立方体 形成初应力文件 形成初应力文件的命令流如下： 当程序运行到solve命令时，产生扩展名为ist的初应力文件。 在重力作用下，竖向位移、应变和应力分别如下三图。 图 2重力作用下的立方体竖向位移云图 图 3重力作用下的立方体竖向应变云图 图 4重力作用下的立方体竖向应力云图 无初位移计算 在形成扩展名为ist的初应力文件和扩展名为db的数据文件后，进行无初位移的命令流如下。 无初位移计算所得的竖向位移、应变和竖向应力如下两图。 图 5重力无初位移的立方体竖向位移云图 从图 5与图 2比较可看出，位移几乎为零，达到了无初位移。 图 6重力无初位移的立方体竖向应变云图 从图 6与图 3相比较可以看出，应变也几乎为零，达到了无初应变。 图 7重力无初位移的立方体竖向应力云图 从图 7与图 4比较可以看出，应力没什么区别。 初应力单独作用计算 为了对比，将3.3中命令流“acel,,10,”改为“acel,,0,”或直接去掉，重新计算。这样相当于载荷只是重力作用下的初应力，然后去掉重力，其结果如下三图。 图 8初应力单独作用下立方体竖向位移云图 图 8说明，无初位移状态下，去掉重力后达到平衡，立方体有与重力方向相反的位移。 图 9初应力单独作用下立方体竖向应变云图 图 9说明，在无初位移状态下去掉重力，竖向有伸长的现象，原来压缩越大，伸长也越大。 图 10初应力单独作用下立方体竖向应力云图 图 10说明，在无初位移状态下，去掉重力后，应力变为零。 结论 初应力可以在初始作用下求得，然后以载荷的形式作用在物体上。这样，初始作用下的位移和应变就变成了零，而初始应力的作用却保留了下来。 投稿日期：2016-03-24 作者：姚再兴（1974-），男，四子王旗人，副教授，博士后。主要从事工程力学、岩土工程及其科学计算。yaozaixing@163.com 1 /PREP7 ET,1,SOLID185 !选185号单元 !以下定义材料参数 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,EX,1,,3e11 MPDATA,PRXY,1,,0.3 MPTEMP,,,,,,,, MPTEMP,1,0 MPDATA,DENS,1,,2000 !生成立方体 BLOCK,,1,,1,,1, !单元划分 MSHAPE,0,3D MSHKEY,1 CM,_Y,VOLU VSEL, , , , 1 CM,_Y1,VOLU CHKMSH,VOLU CMSEL,S,_Y VMESH,_Y1 CMDELE,_Y CMDELE,_Y1 CMDELE,_Y2 !加约束 FLST,2,1,5,ORDE,1 FITEM,2,3 /GO DA,P51X,ALL, !底部约束 FINISH save !形成扩展名为db的文件，以备后续计算时使用 /SOL acel,,10, !加重力 ISWRITE,1 !形成初应力文件 solve finish /clear resume