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

1、1.输出数据控制。指定要输入到D3PLOT、D3PART、D3THDT文件中的二进制数据。* DATABASE. EXTENT. BI NARYNEIPHNEIPSM双NT STRFBG SIGFLG EPSFLG RLTFLG ENGFLGI "I 0"|1? I I 2 I ICMPFLG IEVERPBEAMIP DCOMP SHGESTSSZ N3THDTQ I I 而 1 I 1 I 1 ININTSLD 【NEIPH】一一写入二进制数据的实体单元额外积分点时间变量的数目 。【NEIPS】一一写入二进制数据的壳单元和厚壳单元每个积分点处额外积分点时间变量的数目。M

2、AXINT】一一写入二进布擞据的壳单元积分点数。如果不是默认值3,则得不到中面的结果。STRFLAG】一一设为1会输出实体单元、壳单元、厚壳单元的应变张量，用于后处理绘 图。对于壳单元和厚壳单元，会输出最外和最内两个积分点处的张量，对于实体单元，只输出一个应变张量。【SIGFLG】一一壳单元数据是否包括应力张量。EQ.1 :包括（默认）EQ.2:不包括EPSFLG 壳单元数据是否包括有效塑性应变。EQ.1 :包括（默认）EQ.2:不包括【RLTFLG】一一壳单元数据是否包括合成应力。EQ.1 :包括（默认）EQ.2:不包括ENGFLG 壳单元数据是否包括内能和厚度。EQ.1 :包括（默认）EQ

3、.2:不包括【CMPFLG】一一实体单元、壳单元和厚壳单元各项异性材料应力应变输出时的局部材料坐标系。EQ.0:全局坐标EQ.1 :局部坐标【IEVERP】 限制数据在1000state之内。EQ.0:每个图形文件可以有不止1个stateEQ.1:每个图形文件只能有 1个state【BEAMIP】一一用于输出的梁单元的积分点数。【DCOMP】一一数据压缩以去除刚体数据。EQ.1 :关闭（默认）。没有刚体数据压缩。EQ.2:开启。激活刚体数据压缩。EQ.3:关闭。没有刚体数据压缩，但节点的速度和加速度被去除。EQ.4:开启。激活刚体数据压缩，同时节点的速度和加速度被去除。【SHGE】一一输出壳单

4、元沙漏能密度。EQ.1 :关闭（默认）。不输出沙漏能。EQ.2:开启。输出沙漏能。【STSSZ】一一输出壳单元时间步、质量和增加的质量。EQ.1 :关闭。（默认）EQ.2:只输出时间步长。EQ.3:输出质量、增加的质量、或时间步长。N3THDT】一一为D3THDT数据设置的能量输出选项。EQ.1 :关闭。能量不写入到 D3THDT数据中。EQ.2:开启（默认）。能量写入到 D3THDT数据中。NINTSLD】一一写入LS-DYNA数据的实体单元积分点数目，默认值为1。对于多个积分点的实体单元，该值可能设为8。如果该值设为1,对于多个积分点的实体单元，将输出一个平均值。2. 接触面二进制数据输出

5、控制*DATABASE.BlNARY_INTFORDTLCDT1 1000.0000【DT】一一输出的时间间隔。【LCDT】指定输出时间间隔的曲线。3. 单元子集的时间历程数据输出控制*DAI8ASE_BI NY_D3THDTDT LCDT |1000 OQO0【DT】一一输出的时间间隔。【LCDT】一一指定输出时间间隔的曲线。4. 完全输出控制*DATA3ASE_BI NARY_D3PLOTDT LCDT | NOBEAM NPLTC PSETID I STATS TSTART IAVG |5.0000|01IOOPT 0【DT】一一输出的时间间隔。【LCDT】一一指定输出时间间隔的曲线。N

6、OBEAM 】EQ.0:离散弹簧单元和阻尼单元写入到D3PLOT或D3PART数据中，并显示为梁单元。单元的整体坐标和合力都被写入到数据中。EQ.1 :离散弹簧单元和阻尼单元不会被写入到D3PLOT或D3PART数据中。EQ.2:与等于0时类似。【NPLTC】以下几项与D3PLOT无关PSETIDISTATSTSTARTIAVGIOOPT5.指定输出文件S (DATABASE. ABSTZ.OPTION -AVSFLTCon t rol C&rd& for ASCII du tpu t BNDOITT DEFGEO DEFORC_ EUOUT GCEOITTGLSTAT0 10

7、0 100.10JNTFORCNIATSUMMOVIEMPGS NCFORC NODFORNODOUTRBDOUToioa f0 .ioc0.10RCFORCRWFORCSBTOUTSECFORC SLEOUT SPCFORC SSSTATSWFORC|0.100 |0.1000.10or | TPRINTTRHISTSPHOLTTBinary Optionsreject 1 ASCIIBINART |ABSTAT】一一气囊统计表。 输出体积、压强、内能、气体质量流入率、 气体质量流出率、 质量、温度、密度。AVSFLT】 AVS数据BNDOUT】一一边界环境的力和能量。输出三个方向的力。【D

8、EFGEO】变形的几何体的文件【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】材料能量。GLSTATJNTFORCMAT

9、SUM动能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二方向合位移

10、二方向合力法向力二方向合速度二方向合力三方向合加速度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二

11、方向中心面积合力【TPRINT】一一结构对的热量输出TRHIST 】追踪质点时间历程信息6. 时间步控制*00N7R0L_TIMESTEPDTINiT m"ISDO TSUMFDT2MSLCTM ERODEMS1ST|0 DOO |0 900|Q|0 00|-7.000e-000I o II【DTNIT】一一初始时间步长。EQ.0.0 :由LS-DYNA决定时间步长。【TSSFAC】一一计算时间步长的缩放因子。默认为0.9。【ISDO】一一4节点壳单元时间步长计算的根据。EQ.0：特征长度=面积/min最长边，最长对角线EQ.1:特征长度=面积/最长对角线EQ.2:时间步长取决于条波

12、速度（bar wave speed）和MAX最短边，面积/min （最长 边，最长对角线）。该选项提供的时间步长相对很大，可能导致计算的不稳定，尤其是在 应用三角形单元时。EQ.3:时间步长取决于最大特征值。该选项适用于材料的声音传播速度渐变的结构。 用于计算最大特征值的计算开销是很有意义的，但时间步长的增长通常考虑不用质量缩放的较短的计算周期。【TSUMIT】一一指定壳单元最小时间步长。当某一单元的时间步长小于给定值时，该单元 的材料属性（弹性模量而不是质量）将被调整，使其时间步长不低于给定值。该选项只适用 于以下材料：MAT_PLASTIC_KINEMA TIC , MAT_POWER_L

13、A W_PLASTICITY , MAT_STRAIN_RA TE_DEPENDENT_PLASTICITY,MAT_PIECEWISE_LINEAR_PLASTICITY。不推荐所谓的刚度缩放选项。下面的 DT2MS选项适用于所有材料和所有单元类型，并且是首选的。如果TSUMIT和DT2MS两个选项都被激活并且TSUMIT值为正，则TSUMIT的值自动置为1E-18,使其功能被屏蔽。如果其 值为负并且其绝对值大于IDT2MS 则TSUMIT |优先应用到质量缩放中，如果其绝对值小于IDT2MS | ,贝U TSUMIT的值自动置为 1E-18。【DT2MS】 控制质量缩放的时间步长。DT2M

14、S 值应用范围大于0用于初始影响无关紧要的准静态分析和时间历程分析W 0默认小于0允许的最小时间步长为TSSFAC* DT2MS ,当且仅当时间步长小于判断标准时，质重缩放才会进仃。该选项可用于质重增加影响不大的瞬态分析。 警告：超单元和 ELMENT_DIRECT_MA TRIX_INPUT 不进行质量缩放，所以DT2MS不影响他们的时间*。这种情况下计算会出错终止，DT2MS应输入一个较小的值。【LCTM】一一限制最大时间步长的曲线。ERODE 到达 TSMIN （见下面卡片 CONTROL_TERMINA TION ）时，实体单元和 t- 壳单元的侵蚀标记。如果此项不设，计算会终止。EQ

15、.0:无侵蚀EQ.1 :有侵蚀【MS1ST】一一限制第一步的质量缩放并且根据之前的时间步确定质量矢量。EQ.0:否EQ.1 :是7. 计算终止时间设置* CONTROL 一 TERMI NAT I ONENDTIM ENDCYC DTMIN ENDENG ENIDMASS| i'oo.ooo| FoqoI FoTI roo"ENDTIM】 强制的计算结束时间。ENDCYC】一一计算循环次数。当达到指定循环次数而没有到达ENDTIM指定的计算结束时间时，同样终止计算。循环次数等于时间步的数目。【DTMIN】初始时间步长的缩放因子，用以决定最小时间步长（TSMIN ）, TSMI

16、N=DTSTART*DTMIN 。式中 DTSTART是由LS-DYNA 决定的初始时间步长。当达到 TSMIN时，LS-DYNA计算终止并输出一个重启动文件。【ENDENG】一一控制计算结束的能量比例变化。ENDMASS】一一控制计算结束的质量比例变化。当且仅当用质量缩放控制最小时间步长时，该选项才起作用。8. 单元控制w OONTROL_SHELLWRRANG ESORT IRNX>< ISTUPD THEORY BWC MITER PROJ|?5 0001 L P |2|1|0【WRPANG】一一壳单元翘曲角度。当某个翘曲角度大于给定值时，会输出警告信息。默 认值为20。【E

17、SORT】一一自动挑选退化的四边形单元，并处理为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个循

18、环重新计算法线方向。（只适用于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和2

19、4号材料）EQ.1: 3 次交叉迭代（iterative plastic ity with 3 iterations ）（默认）EQ.2：完全迭代。EQ.3:不迭代。可能导致错误，慎用。【PROJ】9. 输出控制* CONTROL _dUTHJTNPOPT NEECHOJ NREFUP IACCOP OPIFSJ IPNINT IKEDPT FLUSH| o ° I JI OF o p mlNPOPT输出文件打印抑制选项。EQ.0:无抑制。EQ.1 :节点坐标，单元连通性，刚性墙定义，初始速度都不打印。【NEECHO】一一echo file打印抑制选项。EQ.0:打印所有数据。EQ.

20、1 :节点数据被抑制。EQ.2:单元数据被抑制。EQ.3:节点数据和单元数据都被抑制。【NREFUP】一一梁单元参考点坐标更新选项。该选项要求每个参考点都与梁单元不重合。EQ.0:不更新。EQ.1 :更新。IACCOP】一一从“ nodout”文件和“ d3thdt”文件中的速度得到的平均加速度。EQ.0:不求平均加速度。（默认）EQ.1 :输出间隔求平均值。EQ.2:嵌入的用户自定义的过滤器。【OPIFS】一一接触面文件的输出间隔。【IPNINT】一一打印第一次循环中所有单元的初始时间步长。EQ.0:打印100个最小时间步长单元的步长。EQ.1:打印所有单元的时间步长。【IKEDIT】一一写

21、入D3HSP文件的问题情况报告的间隔步数。如果GLSTAT文件已经输出,则该项忽略不计。【IFLUSH】写入缓存的时间步间隔数。缺省值为5000。10. 沙漏控制* CONTROL_HOURGLASSIHQ QH4 I 0.020IHQ沙漏控制类型：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时，激活翘曲刚度，以得到精 确解。该选项会增加

22、25%的计算开销。在壳单元中，IHQ<4的是基于Belyschko-Tsay公式的粘性沙漏控制模式，IHQ=4 , 5, 6 为刚度控制模式。刚度控制模式在大变形问题中可能使响应变得过于刚硬，使用时要注意。在高速问题中推荐采用粘性模式，在低速问题中推荐采用刚度模式。对于大变形问题，推荐使用选项3或5。【QH】一一沙漏系数。该值如果超过0.15可能引起计算不稳定。缺省值为0.1,可适用于除IHQ=6以外的所有选项。备注：对个别组件的沙漏控制，可通过先建立沙漏属性集合器，再从组件集合器中调用沙漏属性的方法实现。11. 能量耗散控制ACONTROL.ENERGY HGEN RWEN SLNTE

23、N RYLENHGEN沙漏能计算选项。该选项需要大量存储空间，并增加10%的计算开销。EQ.1 :不计算沙漏能。（默认）EQ.2:计算沙漏能并包含在能量平衡中，计算结果写入GLSTAT和MATSUM文件中。【RWEN】一一刚性墙能量耗散选项。EQ.1 :不计算刚性墙能量耗散。EQ.2:计算刚性墙能量耗散并包含在能量平衡中，计算结果写入 认）GLSTAT文件中。（默【SLNTEN】一一滑移面能量耗散控制选项。EQ.1:不计算滑移面能量耗散。（当接触激活时，该选项自动设为EQ.2:计算滑移面能量耗散并包含在能量平衡中，计算结果写入 文件中。2)GLSTAT 和 SLEOUT【RYLEN】一一阻尼衰

24、减能量耗散控制选项EQ.1 :不计算阻尼衰减能量耗散。（默认）EQ.2:计算阻尼衰减能量耗散并包含在能量平衡中，计算结果写入GLSTAT文件中。12.CPU时间控制* CONTROL .CPUCPLTTIM |I oTooo"【CPUTIM】一一用于电流相位分析或重启动。EQ.0:没有CPU时间限制。13.接触控制CONTROLCONTACTSLSFACRWPNALISLCHKSHLTHKJ PENOPT THKCHGORIENENMASS0 id。I 1 ooT 2 I 0 I 1 I 0 I 2USRSTRUSRFRCNSBCSINTERM XPENE SSTHKECDT71ED

25、PRJr |cT |f u | Oo g I o I oSLSFAC 滑移面惩罚因子。EQ.0:缺省值=0.1【RWPNAL】一一刚性墙惩罚因子，用来处理刚体与固定刚性墙的相互作用。 EQ.0.0:不考虑刚体与刚性墙的相互作用。GT.0.0 :刚体与刚性墙进行相互作用，推荐值为1。每个从节点将保存 7个变量。如果所有节点对于刚性墙都是从节点，则会大大增加对内存的要求。【ISLCHK】接触面初始穿透检查。EQ.0:自动设为1EQ.1 :不检查。EQ.2:完全检查初始穿透。SHLTHK】一一在面-面接触和点-面接触类型中考虑壳单元厚度的选项。选项 1和2会激 活新的接触算法。厚度偏置通常包括在单面

26、接触、约束算法、自动面面接触和自动点面接触类型中。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】一一初始化过程中接触

27、面截面自动再定位选项。EQ.0:自动设为1。EQ.1:仅自动（part）输入时激活。接触面由part定义。EQ.2:手动（segment）和自动输入（part）都激活。EQ.3:不激活。【ENMASS】一一对接触过程中销蚀掉的节点的质量的处理。该选项影响所有当周围单元 失效而自动移除相应节点的接触类型。通常，销蚀掉的节点的移除会使计算更稳定，但是质量的减少会导致错误的结果。EQ.0:从计算中移除销蚀的节点。（默认）EQ.1 :保留体单元销蚀的节点并在接触中继续起作用。EQ.2 :保留体单元和壳单元销蚀的节点并在接触中继续起作用。【USRSTR】USRFRC】【NSBCS】INTERM 】【XP

28、ENE】【SSTHK】【ECDT】【TIEDPRJ 】14,体积粘度控制。 体积粘度是为了解决冲击波。* CONTROL.BULK_VISCOSI TYQ1 Q2 IBQ|1 .500 |。、即-1|Q1缺省的二次粘度系数（1.5）。【Q2】一一缺省的线性粘度系数（0.06）。【IB Q】一一体积粘度类型。EQ.-1 :标准。（对于单元类型为 2, 10, 16的壳单元）EQ.+1 :标准。（默认）15.精确性控制* CONTROL.ACCURACYOSU INN FTOSU |口 I 2 I 0OSU第二 目标应力（objective stress）更新的全局标志。一般来讲，只有像轮胎那样有

29、大的转动的显式计算才会需要该选项。EQ.0:关。（默认）EQ.1 :开。【INN】一一壳单元和体单元节点的编号方式。EQ.1 :关。（默认为显式计算）EQ.2:开。仅对壳单元。（默认为隐式计算）EQ.3:开。仅对体单元。EQ.4:开。对壳单元和体单元。【PIDOSU】一一需要目标应力更新的组件集的ID号。如果该选项给定，则只有相应的组件集会应用目标应力更新，OSU选项将设置被忽略。16.全部数据输出*DATABASE. DTBl N1ARY_D3PLOTLCOT NOBfEAM NPLTCPSET1DISTATS (TSTARTJ(IAVG|0 005IOOPT 010|0101【DT】一一输

30、出的时间间隔。【CYCL】一一时间步内的输出间隔（一个时间步长是一个循环）。一般不用。NOBEAM 】关于是 DATABESE-BINARY -D3POLT 或 DATABESE-BINARY -D3PART 的选择标志。EQ.0 :被描述成 beam单元的离散的弹簧和减震器单元添加到D3POLT或D3PART的数据中。单元的球形坐标 xyz还有合力也添加进去。EQ.1 :非不连续的弹簧和减震器单元添加进去。 旧的数据要转换的 KEYWORD 中时这个选择被选择。在旧的数据中没有必须为beam和弹簧单元创建单独的 ID号，然而弹簧和减震器单元也用 beam表示，这样就会出错。EQ.2 :同0一样。在be