LS-DYNA常见问题汇总1.0
资料来源:网络和自己的总结 yuminhust2005
Copyright of original English version owned by relative author. Chinese version owned by Aifea.com/Kevin 目录
1.Consistent system of units 单位制度 .............................................................................. 2 2.Mass Scaling 质量缩放.................................................................................................. 2 3.Long run times 长分析时间........................................................................................... 6 4.Quasi-static 准静态 ....................................................................................................... 8 5.Instability 计算不稳定 ................................................................................................ 10 6.Negative Volume 负体积.............................................................................................. 12 7.Energy balance 能量平衡 ............................................................................................ 14 8.Hourglass control 沙漏控制 ........................................................................................ 19 9.Damping 阻尼 ............................................................................................................ 23 10.ASCII output for MPP via binout ............................................................................... 27 11.Contact Overview 接触概述 ...................................................................................... 30 12.Contact Soft 1 接触Soft=1 ........................................................................................ 33 13.LS-DYNA中夹层板(sandwich)的模拟 ........................................................................ 35 14. 怎样进行二次开发 ................................................................................................... 37
1.Consistent system of units 单位制度
相信做仿真分析的人第一个需要明确的就是一致单位系统(Consistent Units)。计算机只认识0&1、只懂得玩数字,它才不管你用的数字的物理意义。而工程师自己负责单位制的统一,否则计算出来的结果没有意义,不幸的是大多数老师在教有限元数值计算时似乎没有提到这一点。见下面LS-DYNA FAQ中的定义: Definition of a consistent system of units (required for LS-DYNA):
1 force unit = 1 mass unit * 1 acceleration unit 1 力单位 = 1 质量单位 × 1 加速度单位 1 acceleration unit = 1 length unit / (1 time unit)^2 1 加速度单位 = 1 长度单位/1 时间单位的平方
The following table provides examples of consistent systems of units.
As points of reference, the mass density and Young‘s Modulus of steel are provided in each system of units. ―GRAVITY‖ is gravitational acceleration.
MASS LENGTH TIME FORCE STRESS ENERGY DENSITY YOUNG’s kg kg kg kg kg gm gm gm gm ton lbfs2/in slug kgfs2/mm kg gm m cm cm cm mm cm cm mm mm mm in ft mm mm cm s s N 1.e-02N Pa GPa dy/cm2 Mbar Pa MPa MPa psi psf 1000Pa 100000Pa Joule erg 7.83E+03 7.83E-03 7.83E-03 7.83E-03 7.83E+00 2.07E+11 2.07E+09 2.07E+03 2.07E-03 2.07E+02 2.07E+12 2.07E+00 2.07E+11 2.07E+05 2.07E+05 3.00E+07 4.32E+09 7.00E+02 7.83E-06 7.83E+00 Velocity GRAVITY (56.3KMPH) 15.65 1.56E+03 1.56 1.56E-03 15.65 1.56E+03 1.56E-03 1.56E+04 15.65 1.56E+04 6.16E+02 51.33 1.56E+04 2.07E+08 2.07E+06 9.806 9.81E+02 9.81E-04 9.81E-10 9.81E-03 9.81E+02 9.81E-10 9.81E+03 9.81E-03 9.81E+03 386 32.17 (Japan) 9.81E+02 ms 1.e+04N us 1.e+10N ms s s ms s s s s s ms KN dyne 1.e-06N N N lbf lbf kgf mN KN-mm 7.83E-06 1.e7Ncm 7.83E+00 7.83E-03 N-mm N-mm lbf-in lbf-ft 7.83E-03 7.83E-09 7.33E-04 15.2 us 1.e+07N kgf/mm2 kgf-mm 8.02E-10 2.Mass Scaling 质量缩放
质量缩放指的是通过增加非物理的质量到结构上从而获得大的显式时间步的技术。
在一个动态分析中,任何时候增加非物理的质量来增大时间步将会影响计算结果(因为F=ma)。有时候这种影响不明显,在这种情况下增加非物理的质量是无可非议的。比如额外的质量只增加到不是关键区域的很少的小单元上或者准静态的分析(速度很小,动能相对峰值内能非常小)。总的来说,是由分析者来判断质量缩放的影响。你可能有必要做另一个减小或消除了质量缩放的分析来估计质量增加对结果的灵敏度。
你可以通过人工有选择的增加一个部件的材料密度来实现质量缩放。这种手动质量缩放的方法是独立于通过设置*Control_timestep卡DT2MS项来实现的自动质量缩放。
当DT2MS设置为一个负值时,质量只是增加到时间步小于TSSFAC*|DT2MS|的单元上。通过增加这些单元的质量,它们的时间达到TSSFAC*|DT2MS|。有无数种TSSFAC和DT2MS的组合可以得到同样的乘积,因而有相同的时间步,但是对于每一种组合增加的质量将是不一样的。一般的趋势是TSSFAC越小,增加的质量越多。 作为回报,当TSSFAC减小时计算稳定性增加(就像在没有做质量缩放的求解中一样)。 如果TSSFAC缺省的值0.9会导致稳定性问题,可以试试0.8或者0.7。 如果你减小TSSFAC,你可以相应增加|DT2MS|,这样还是可以保证时间步乘积不变。
为了确定什么时候和位置质量自动增加了,可以输出GLSTAT和MATSUM文件。这些文件允许你绘出完整的模型或者单独部件所增加的质量对时间的曲线。为了得到由壳单元组成的部件增加的质量云图,将*database_extent_binary卡的STSSZ项设置为3。 这样你可以用ls-prepost绘出每个单元的质量增加量的云图,具体方法是通过选择Fcomp>Misc>time step size。
在*control_timestep中设置DT2MS正值和负值的不同之处如下:
负值:初始时间步将不会小于TSSFAC*-DT2MS。质量只是增加到时间步小于TSSFAC*|DT2MS|的单元上。当质量缩放可接受时,推荐用这种方法。用这种方法时质量增量是有限的。过多的增加质量会导致计算任务终止。
正值:初始时间将不会小于DT2MS。 单元质量会增加或者减小以保证每一个单元的时间步都一样。这种方法尽管不会因为过多增加质量而导致计算终止,但更难以作出合理的解释。
*control_timestep卡中的参数MS1ST控制是否只是在初始化时增加一次质量(MS1ST=1)还是任何需要维持由DT2MS所指定的时间步时都增加质量(MS1ST=0)。
你可以通过在*control_termination卡片中设置参数ENDMAS来控制当质量增加到初始质量一定比率时终止计算(只对自动质量缩放有效)
-------------------------- 可变形点焊梁的质量缩放
*mat_spotweld卡的质量缩放参数DT只影响点焊单元。如果*control_timestep卡中没有指定质量缩放(DT2MS=0),而且时间由可变形点焊控制,可以用参数DT来在初始化时增加惯量到点焊单元上来提高时间步达到DT指定的值。当DT不为0时,增加到可变形点焊梁元上的质量会输出到d3hsp文件里。MATSUM 中动量和动能不受增加到可变形点焊上的质量的影响。GSLTAT中DOES和总的KE受增加的质量的影响。
考虑三种调用可变形点焊的质量缩放的情况:
1.当DT2MS为负值*mat_spotweld卡DT=0时,尽管在d3hsp文件中可变形点焊质量增量百分比不真实。下面几个值是正确的:d3hsp中‖added spotweld mass‖; 第一个时间步之后的‖added mass‖ & ―percentage increase‖; glstat和matsum中的‖added mass‖。
2. 当DT2MS为负值且*mat_spotweld卡DT≠0时,可变形点焊质量增加不会包含在d3hsp、glstat、matsum文件中的‖added mass‖里。这非常容易令人误解。用户必须检查d3hsp文件的‖added spotweld mass‖。建议不要同时使用两种质量缩放标准,推荐使用第一种方法(即负的DT2MS&DT=0)。