LSDYNA刚体转动问题.doc

几乎看了本版块所有关于转动的讨论，现在终于对刚体的转动彻底的了解了，现在把前一段时间的心得拿出来，希望能对和我一样正在为这个问题迷茫的人有所帮助，以下两个文件是用ansys命令流写的，因为关于转动的k文件在“关于刚体转动的问题”这个主题里面已经讨论的很清除了，相信有很多人都是用ansys做前处理和后处理的，所以即使看懂了k文件，也未必就知道用ansys的前处理怎么来实现。关于刚体转动，主要要注意以下几点： ??????? 1. 如果是绕自己的质心转动，那很简单，不用设转动惯量和转动轴，因为ls-dyna默认就是绕着质心转动的。??????? 2.如果是绕任意轴转动，那么就必须设新的质心center、惯性张量和转动轴。这主要是用edlcs和edipart两个命令来实现，前一个设转动轴， 后一个设center和惯性张量，关于惯性张量可以用solidworks来计算（wswu在那个主题曾提到这个方法），当然惯性张量也可以自己估算,如 果质量体相对坐标轴比较远就可以把质量体当作集中质量，比如计算IXY=m*x*y(x和y是把质量中心投影到xy平面后的那个点距x和y的距离）。有人 曾提到惯性张量可以随便给，当然随便给可能能算出结果来，但是那个结果可能并不是你所想要的结果，所以，如果是针对课题，这些东西不妨花点时间，弄出一个 确定的值。??????? 其实也没有什么太多要注意的，以上两点都做好了，你肯定就ok了。??????? 我这两个文件借助了iamaxelong提出的用刚体带动弹性体的方法，因为iamaxelong老大那个ansys命令流贴在网页上被弄了一些乱码在里 面，复制过来之后要做大量的修改，而且分析所需要的时间太长，所以我这里整了一个特别简单的例子，我的电脑分析只需要1分钟，可以帮助你快速看到分析的结 果，相信即使比我还菜的人也能很快了解刚体转动的问题啦。 /PREP7 !定义单元和材料?ET,1,SOLID164??MP,DENS,1,7800?MP,EX,1,2.1e11?MP,NUXY,1,0.3??EDMP,RIGI,2,7,4MP,DENS,2,7800?MP,EX,2,2.1e11?MP,NUXY,2,0.3? !画刚体和弹性体的实体模型BLOCK,-0.1,0.1,0.4,0.6,0,0.05,??BLOCK,-0.05,0.05,0.45,0.55,0.05,0.1, KWPAVE,?????? 9wpro,,90.000000,VSBW,?????? 1?? KWPAVE,????? 10VSBW,?????? 3?? ?VGLUE,all? !确定实体模型上线条的网格密度LESIZE,all,0.025, , , , , , ,1?lesize,11,,,2lesize,12,,,2lesize,31,,,2lesize,32,,,2lesize,9,,,2lesize,10,,,2lesize,39,,,2lesize,40,,,2lesize,21,,,2lesize,22,,,2lesize,23,,,2lesize,24,,,2 numcmp,volu !划分单元vsweep,2vsweep,3vsweep,4 mat,2vsweep,1 !建立partedpart,create !定义载荷数组*DIM,time,ARRAY,2,1,1, , ,??*SET,TIME(2,1,1) , 10??*DIM,wz,ARRAY,2,1,1, , ,*SET,WZ(1,1,1) , 5?*SET,WZ(2,1,1) , 5 !加载EDLOAD,ADD,RBOZ,0,?????? 2,TIME,WZ, 0, , , , , !定义分析选项FINISH?/SOLTIME,0.5,??EDRST,30,??EDHTIME,50,EDDUMP,1,????EDENERGY,1,1,1,1/STATUS,SOLUSOLVE?? !输出动画FINISH?/POST1?SET,LAST/DSCALE,1,1.0??/REPLOT?PLDI, ,ANTIME,50,0.1, ,1,2,0,0?FINISH? ? ? /PREP7 !定义单元和材料?ET,1,SOLID164??MP,DENS,1,7800?MP,EX,1,2.1e11?MP,NUXY,1,0.3??EDMP,RIGI,2,7,4MP,DENS,2,7800?MP,EX,2,2.1e11?MP,NUXY,2,0.3? !画刚体和弹性体的实体模型BLOCK,-0.1,0.1,0.4,0.6,0,0.05,??BLOCK,