ls-dyna计算稳定性的讨论

lzkhnu

作为新人，没事看看帖子，浏览浏览博文，看到有关dyna计算不稳定的介绍。在这里将所学知识总结一下。
我们知道dyna分单精度版和双精度版。
单精度的是S  g! |6 T1 B5 Y7 G$ A8 ?  S+ L# h
ls971_s_R4.2_win32_p2 W+ Z& ]* w. M5 S! b  ^  J* Y& R) s
双精度的是D+ T: f/ B2 }) w1 B
ls971_d_R4.2_win32_p
从字面上理解双精度版精度会高些。
然后就是一些特定的问题需要双精度结果才会准确，如隐式求解。具体见
http://forum.simwe.com/viewthread.php?tid=875713&highlight=%CB%AB%BE%AB%B6%C8
关于dyna（即便是双精度）运算常出现的不稳定。
一些计算不稳定的提示消息：
“out-of-range velocities” 速度超出范围
“negative volume in brick element” 体单元负体积
“termination due to mass increase” 因质量增加而终止
用来克服显式求解中的不稳定的方法如下：
增加d3plot的输出频率到可以显示出不稳定的出现过程。这可以提供导致不稳定性发生的线索。
其它的一些些解决数值不稳定性的技巧：
* 试着减小时间步(timestep)缩放系数(即使使用了质量缩放mass-scaling)
* 单元类型和/或沙漏(hourglass)控制。对出现不稳定的减缩体和壳单元，试着用沙漏控制type 4 和沙漏系数0.05
。或者试着用类型16的壳单元，沙漏控制type 8。如果壳响应主要是弹性，设置BWC=1 和 PROJ=1 (仅对B-T壳)。避免使用type=2体单元。对体单元部件，在厚度方向最少用两个体单元。
* 接触。设置接触的bucket sorts之间周期数为0，这样会使用缺省的分类间隔。如果参与接触的两个部件的相对速度异常的大，可能需要减小bucket sort的间隔(比如减小到5,2甚至1)。
如果仿真过程中有明显的接触穿透出现，转换到使用*contact_automatic_surface_to_surface或者*contact_automatic_single_surface，并设置SOFT=1。确保几何考虑了壳单元的厚度。如果壳非常薄，比如小于1mm，放大或者设置接触厚度到一个更加合理的值。
* 避免冗余的接触定义，也就是说不要对同样的两个部件定义多于一个的接触对。
* 查找出现不稳定的部件的材料定义中的错误(比如误输入，不一致的单位系统等)
* 关掉所有的*damping
这些技巧是一些通用的方法，可能并不适合于所有的情况。
具体见链接
http://blog.sina.com.cn/s/blog_6817db3a0100llgu.html
关于双精度的不稳定讨论还有
http://forum.simwe.com/viewthread.php?tid=857455&highlight=%CB%AB%BE%AB%B6%C8
看着可能有点用，就整理了下。若有重复或嫌水平低也莫口水啊。

lojade

我的课题里用单精度版本的算就是不对，很容易出现负体积，只能用双精度版的计算。
我做的是动力学方面的，莫非我遇到的就是“一些特定的问题”之一？

lzkhnu

把你的问题贴出来让各位大牛们帮解决下啊。顺便我也好取取经

lojade

3# lzkhnu
我的问题用双精度版本算就解决了呀。

qgchen

偶是新手，问下楼主单双精度是如何转换的？

qgchen

刚才自己找到了，在求解界面里找到设置项了，刚学DYNA不久，让大家见笑了