LS-DYNA常见问题汇总2.0 yuminhust2005 2008-10-12 3. 质量缩放
质量缩放指的是通过增加非物理的质量到结构上从而获得大的显式时间步的技术。
在一个动态分析中，任何时候增加非物理的质量来增大时间步将会影响计算结果（因为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）。
3. 如果DT2MS＝0且DT≠0，初始时间步将不考虑增加点焊的质量，但是之后每一个周期时间步都会增加10%，直到时间步达到正确的值（考虑点焊质量增加）。glstat & matsum不包含‖added mass‖的行。
注意质量增加会引起能量比率增长。
第 6 页
LS-DYNA常见问题汇总2.0 yuminhust2005 2008-10-12 4. 长分析时间
当用显式时间积分时，对于仿真非常小的部件而分析时间又要相当长时没有好的方法。质量缩放(mass-scaling)增加了需要确认非物理质量的增加不会显著影响计算结果的负担。当使用时间缩放(time-scaling)时也有同样的问题。时间缩放(time-scaling)是指为了减小需要的时间步数，通过增加加载速率而缩短仿真时间。
要确认时间步不是仅由很少的小单元或者刚度大单元控制，可以通过在d3hsp文件中搜索‖smallest‖来显示100个最小的时间步单元。如果只有很少的几个单元控制时间步，可以把那些单元及邻近区域重新remesh或者把它们变成刚体。
可是仅运行必要长的时间是很明显的。这意味着在一个跌落分析的情况时，给跌落物体一个初速度，把它放在离地面一个非常小的距离。冲击之后，仅运行足够得到需要的结果的时间。
值得注意的是对于一个长时间的仿真，如果时间步数超过了50万步，最好使用双精度版本的LS-DYNA求解器，使截断误差的影响最小化。运行双精度版本要增加30%的时间。
对于长时间的分析，自动显式/隐式转换可能是一个选择。使用这种方法，用户可以指定在一个时间段内使用隐式积分。隐式积分的优点是时间步不由单元尺寸控制，所以可以得到大的时间步。当然，隐式计算也非常点用cpu时间。而且，目前并不是所有的LS-DYNA的功能和材料都在隐式分析中实现(大部分已经实现)。下面的FEA information newsletter里讨论了显式/隐式转换(http://www.feapublications.com/pages/pdfnews/3feadec.pdf)。
See also: mass_scaling, quasistatic
第 7 页
LS-DYNA常见问题汇总2.0 yuminhust2005 2008-10-12 5. 准静态
动态松驰(Dynamic relaxation)并不是有意为一般的准静态(quasi-static)分析设置的。它适合于当预载只产生小的弹情况应变的施加预载，或者初始化系统到一个预定义的几何形状[1]。但对其它更多情况并不适合。
你可以通过做一个常规的显示仿真来模拟准静态分析，通过按需要调用时间/质量缩放(time-scaling,mass-scaling)来在可接受的时间内得到结果，但这种方法是需要技巧地。你必须监测系统动能按希望的使惯性效应最小化。基本上动能相对内能应该保持在一个较小的值。时间缩放是指加载比在准静态实验里更快，以减少总的仿真时间。关于质量缩放更多内容可以看‖mass_scaling‖一节。或者你可以尝试用LS-DYNA运行一个隐式静力分析。可以看用户手册里的卡片*control_implicit_…‖和Appendix M。
See also: gravity.txt, readme.preload, mass_scaling, long_run_times, implicit.general, quick_initialization. Note[1]：初始化到预定义的几何
1. 从第一次分析的最终状态输出一个节点位移文件。（这一部分未按原文翻译）注意d3plot文件里不包含节点转动信息，因此转动输出为0。这对初始化壳和梁单元会是个问题。LS-Prepost有一个选项是输出节点位移，在Output->Nodal Displacements里。但是这个输出是i8,3e16格式的，但需要的是i8,3e15，所以要注意修改一下。
如果你做了一个正常的动态松驰分析来得到初始状态，一个预定义位移和转动的drdisp.sif文件在DR阶段结束时会自动创建。
2. 在第二次分析时，快速的初始化到第一步输出的预定义的几何。你需要设置卡片*control_dynamic_relaxation里的参数IDRFLG=2，而且在命令行里指定‖m=filename‖(其中filename指第一步创建的文件)。这样在瞬态分析之前，LS-DYNA会自动做一个100步的预分析来使节点根据文件filename指定的数据移动到指定值。
第 8 页
LS-DYNA常见问题汇总2.0 yuminhust2005 2008-10-12 6. 计算不稳定
一些表示计算不稳定的消息如：
―out-of-range velocities‖ 速度超出范围
―negative volume in brick element‖ 体单元负体积 ―termination due to mass increase‖ 因质量增加而终止 用来克服显式求解中的不稳定的方法如下：
首先(也是最重要的)是使用可获得的最新的LS-DYNA版本。最新的执行块可以从ftp://user@ftp.lstc.com上下载(注：前提是你有访问权限)。联系LSTC获得user帐号的密码。最新的BETA版执行块可以在ftp://ftp.lstc.com/outgoing/ls971上找到(不需要密码，但lstc公司对ftp访问有IP限制)。
其次是增加d3plot的输出频率到可以显示出不稳定的出现过程。这可以提供导致不稳定性发生的线索。 其它的不些解决数值不稳定性的技巧： * 试着用双精度LS-DYNA版本运行一次
* 试着减小时间步(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
这些技巧是一些通用的方法，可能并不适合于所有的情况。
第 9 页