碰撞模拟中的控制卡片设置及意义

1.指定输出文件【ABSTAT】——气囊统计表。输出体积、压强、内能、气体质量流入率、气体质量流出率、质量、温度、密度。【AVSFLT】——AVS数据【BNDOUT】——边界环境的力和能量。输出三个方向的力。【DEFGEO】——变形的几何体的文件【DEFORC】——离散单元。输出三个方向的力。【ELOUT】——单元数据。（见DATABASE_HISTORY_OPTION）梁单元平面应力块平面应变轴向合力xx,yy,zz应力xx,yy,zz应力xx,yy,zz应变S方向剪切合力xy,yz,zx应力xy,yz,zx应力xy,yz,zx应变T方向剪切合力塑性应变有效应力下表面应变S方向合力矩屈服函数上表面应变T方向合力矩扭力合力【GCEOUT】——几何接触实体。包含三个方向力和力矩。【GLSTAT】——总体数据。【JNTFORC】——运动副力文件【NATSUM】——材料能量。GLSTATJNTFORCMATSUM动能x,y,z三方向的力动能内能x,y,z三方向的力矩内能总能量沙漏能比率x,y,z三方向的动量刚性墙能量x,y,z三方向的刚体速度弹簧和阻尼能量总动能沙漏能总内能阻尼能总沙漏能滑移面能量外功x,y,z三方向速度时间步单元ID号控制的时间步【MOVIE】——【MPGS】——【NCFORC】——接触面节点力【NODFOR】——节点力组【NODOUT】——节点数据NCFORCNODOUTNODFORx方向力位移x,y,z三方向力y方向力速度z方向力加速度转动量角速度角加速度【RBDOUT】——刚体数据【RCFORC】——接触面合成力【RWFORC】——刚性墙所受的力RBDOUTRCFORCRWFORC三方向合位移三方向合力法向力三方向合速度三方向合力三方向合加速度【SBTOUT】——安全带输出文件【SECFORC】——横截面通过的力（见DATABASE_CROSS_SECTION_OPTION）【SLEOUT】——滑移面的能量【SPCFORC】——单点约束的反作用力【SPHOUT】——SPH数据（见DATABASE_HISTORY_OPTION）【SSSTAT】——子系统数据【SWFORC】——节点约束反力（焊点和铆钉）SECFORCSLEOUTSPCFORCSWFORCx,y,z三方向力Slave能量x,y,z三方向力轴向力x,y,z三方向力矩Master能量x,y,z三方向力矩剪切力x,y,z三方向中心面积合力【TPRINT】——结构对的热量输出【TRHIST】——追踪质点时间历程信息2.时间步控制【DTNIT】——初始时间步长。EQ.0.0：由LS-DYNA决定时间步长。【TSSFAC】——计算时间步长的缩放因子。默认为0.9。【ISDO】——4节点壳单元时间步长计算的根据。EQ.0：特征长度=面积/min{最长边，最长对角线}EQ.1：特征长度=面积/最长对角线EQ.2：时间步长取决于条波速度（barwavespeed）和MAX{最短边，面积/min（最长边，最长对角线）}。该选项提供的时间步长相对很大，可能导致计算的不稳定，尤其是在应用三角形单元时。EQ.3：时间步长取决于最大特征值。该选项适用于材料的声音传播速度渐变的结构。用于计算最大特征值的计算开销是很有意义的，但时间步长的增长通常考虑不用质量缩放的较短的计算周期。【TSUMIT】——指定壳单元最小时间步长。当某一单元的时间步长小于给定值时，该单元的材料属性（弹性模量而不是质量）将被调整，使其时间步长不低于给定值。该选项只适用于以下材料：MAT_PLASTIC_KINEMATIC，MAT_POWER_LAW_PLASTICITY，MAT_STRAIN_RATE_DEPENDENT_PLASTICITY，MAT_PIECEWISE_LINEAR_PLASTICITY。不推荐所谓的刚度缩放选项。下面的DT2MS选项适用于所有材料和所有单元类型，并且是首选的。如果TSUMIT和DT2MS两个选项都被激活并且TSUMIT值为正，则TSUMIT的值自动置为1E-18，使其功能被屏蔽。如果其值为负并且其绝对值大于│DT2MS│，则│TSUMIT│优先应用到质量缩放中，如果其绝对值小于│DT2MS│，则TSUMIT的值自动置为1E-18。【DT2MS】——控制质量缩放的时间步长。DT2MS值应用范围大于0用于初始影响无关紧要的准静态分析和时间历程分析等于0默认小于0允许的最小时间步长为TSSFAC*│DT2MS│，当且仅当时间步长小于判断标准时，质量缩放才会进行。该选项可用于质量增加影响不大的瞬态分析。警告：超单元和ELMENT_DIRECT_MATRIX_INPUT不进行质量缩放，所以DT2MS不影响他们的时间步长。这种情况下计算会出错终止，DT2MS应输入一个较小的值。【LCTM】——限制最大时间步长的曲线。【ERODE】——到达TSMIN（见下面卡片CONTROL_TERMINATION）时，实体单元和t-壳单元的侵蚀标记。如果此项不设，计算会终止。EQ.0：无侵蚀EQ.1：有侵蚀【MS1ST】——限制第一步的质量缩放并且根据之前的时间步确定质量矢量。·····EQ.0：否EQ.1：是3.计算终止时间设置【ENDTIM】——强制的计算结束时间。【ENDCYC】——计算循环次数。当达到指定循环次数而没有到达ENDTIM指定的计算结束时间时，同样终止计算。循环次数等于时间步的数目。【DTMIN】——初始时间步长的缩放因子，用以决定最小时间步长（TSMIN），TSMIN=DTSTART*DTMIN。式中DTSTART是由LS-DYNA决定的初始时间步长。当达到TSMIN时，LS-DYNA计算终止并输出一个重启动文件。【ENDENG】——控制计算结束的能量比例变化。【ENDMASS】——控制计算结束的质量比例变化。当且仅当用质量缩放控制最小时间步长时，该选项才起作用。4.单元控制【WRPANG】——壳单元翘曲角度。当某个翘曲角度大于给定值时，会输出警告信息。默认值为20。【ESORT】——自动挑选退化的四边形单元，并处理为CO三角形单元公式，以保证求解稳定。EQ.0：不挑选（默认）EQ.1：完全挑选并处理【IRNXX】——单元法线更新选项。该选项影响Hughes-Liu，Belytschko-Wong-Chiang，和Belytschko-Tsay单元公式。当且仅当翘曲刚度选项被激活时，即BCW=1时，以上单元公式才受影响。对于Hughes-Liu壳单元类型1，6和7，IRNXX必须设为-2以调用上表面或下表面作为参考面。EQ.-2：··········EQ.-1：每个循环都重新计算法线方向。EQ.0：自动设为-1。EQ.1：重启动时计算。EQ.n：每n个循环重新计算法线方向。（只适用于Hughes-Liu壳单元类型）【ISTUPD】——单元厚度变化选项。该选项影响所有单元公式EQ.0：不变化。EQ.1：膜变形引起厚度改变。该选项对金属板料成型或所有膜片拉伸作用很大的情况都很重要。【THEORY】——缺省的壳单元理论。默认值为2，Belytschko-Tsay单元公式。【BWC】——Belytschko-Tsay单元公式的翘曲刚度选项。EQ.1：增加Belytschko-Wong-Chiang公式的翘曲刚度。EQ.2：Belytschko-Tsay单元公式。不增加翘曲刚度。（默认）【MITER】——平面应力塑性选项。（适用于3，18，19和24号材料）EQ.1：3次交叉迭代（iterativeplasticitywith3iterations）（默认）EQ.2：完全迭代。EQ.3：不迭代。可能导致错误，慎用。【PROJ】——········5.沙漏控制【IHQ】——沙漏控制类型：EQ.1：标准LS-DYNA类型。（默认）EQ.2：Flanagan-Belyschko积分类型。EQ.3：有精确体积的Flanagan-Belyschko积分类型。EQ.4：类型2的刚度形式。EQ.5：类型3的刚度形式。EQ.6：········EQ.8：适用于单元类型为16的全积分壳单元。当IHQ=8时，激活翘曲刚度，以得到精确解。该选项会增加25%的计算开销。在壳单元中，IHQ4的是基于Belyschko-Tsay公式的粘性沙漏控制模式，IHQ=4，5，6为刚度控制模式。刚度控制模式在大变形问题中可能使响应变得过于刚硬，使用时要注意。在高速问题中推荐采用粘性模式，在低速问题中推荐采用刚度模式。对于大变形问题，推荐使用选项3或5。【QH】——沙漏系数。该值如果超过0.15可能引起计算不稳定。缺省值为0.1，可适用于除IHQ=6以外的所有选项。备注：对个别组件的沙漏控制，可通过先建立沙漏属性集合器，再从组件集合器中调用沙漏属性的方法实现。6.能量耗散控制【HGEN】——沙漏能计算选项。该选项需要大量存储空间，并增加10%的计算开销。EQ.1：不计算沙漏能。（默认）EQ.2：计算沙漏能并包含在能量平衡中，计算结果写入GLSTAT和MATSUM文件中。【RWEN】——刚性墙能量耗散选项。EQ.1：不计算刚性墙能量耗散。EQ.2：计算刚性墙能量耗散并包含在能量平衡中，计算结果写入GLSTAT文件中。（默认）【SLNTEN】——滑移面能量耗散控制选项。EQ.1：不计算滑移面能量耗散。（当接触激活时，该选项自动设为2）EQ.2：计算滑移面能量耗散并包含在能量平衡中，计算结果写入GLSTAT和SLEOUT文件中。【RYLEN】——阻尼衰减能量耗散控制选项EQ.1：不计算阻尼衰减能量耗散。（默认）EQ.2：计算阻尼衰减能量耗散并包含在能量平衡中，计算结果写入GLSTAT文件中。7.CPU时间控制【CPUTIM】——用于电流相位分析或重启动。EQ.0：没有CPU时间限制。8.接触控制【SLSFAC】——滑移面惩罚因子。EQ.0：缺省值=0.1【RWPNAL】——刚性墙惩罚因子，用来处理刚体与固定刚性墙的相互作用。·······EQ.0.0：不考虑刚体与刚性墙的相互作用。GT.0.0：刚体与刚性墙进行相互作用，推荐值为1。每个从节点将保存7个变量。如果所有节点对于刚性墙都是从节点，则会大大增加对内存的要求。【ISLCHK】——接触面初始穿透检查。EQ.0：自动设为1EQ.1：不检查。EQ.2：完全检查初始穿透。【SHLTHK】——在面-面接触和点-面接触类型中考虑壳单元厚度的选项。选项1和2会激活新的接触算法。厚度偏置通常包括在单面接触、约束算法、自动面面接触和自动点面接触类型中。EQ.0：不考虑厚度偏置。EQ.1：考虑厚度偏置但刚体除外。EQ.2：考虑厚度偏置，包括刚体。【PENOPT】——惩罚刚度值选项。EQ.0：自动设为1EQ.1：主动截面和从节点刚度的最小值。（默认）EQ.2：用主动截面的刚度值。（过去的方法）EQ.3：用从节点的刚度值。EQ.4：用从节点的刚度值，面积或质量加权。EQ.5：与4相同，但是厚度加权。通常不推荐使用。选项4和5推荐在金属成型计算中使用。【THKCHG】——单面接触中考虑壳单元厚度变化的选项。EQ.0：不考虑。（默认）EQ.1：考虑壳单元厚度变化。【ORIEN】——初始化过程中接触面截面自动再定位选项。EQ.0：自动设为1。EQ.1：仅自动（part）输入时激活。接触面由part定义。EQ.2：手动（segment）和自动输入（part）都激活。EQ.3：不激活。【ENMASS】——对接触过程中销蚀掉的节点的质量的处理。该选项影响所有当周围单元失效而自动移除相应节点的接触类型。通常，销蚀掉的节点的移除会使计算更稳定，但是质量的减少会导致错误的结果。EQ.0：从计算中移除销蚀的节点。（默认）EQ.1：保留体单元销蚀的节点并在接触中继续起作用。EQ.2：保留体单元和壳单元销蚀的节点并在接触中继续起作用。【USRSTR】——·········【USRFRC】——·········【NSBCS】——【INTERM】【XPENE】【SSTHK】【ECDT】【TIEDPRJ】9