笔记LS-DYNA模拟中能量及其显示及单元与网格的区别.docx

《笔记LS-DYNA模拟中能量及其显示及单元与网格的区别.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。

1、GLSTAT(参见*database_glstat)文件中报告的总能量是下面几种能量的和内能internalenergy动能kineticenergy接触(滑移)能contact(sliding)energy沙漏能hourglassenergy系统阻尼能systemdampingenergy刚性墙能量rigidwallenergyGLSTAT中报告的弹簧阻尼能“Springanddamperenergy”是离散单元(discreteelements)、安全带单元(seatbeltelements)内能及和铰链刚度相关的内能(*constrained_joint_stiffn

2、ess…)的总和。而内能“InternalEnergy”包含弹簧阻尼能“Springanddamperenergy”和所有其它单元的内能。因此弹簧阻尼能“Springanddamperenergy”是内能“Internalenergy”的子集。由SMP5434a版输出到glstat文件中的铰链内能“jointinternalenergy”跟*constrained_joing_stiffness不相关。它似乎与*constrained_joint_revolute(_spherical,etc)的罚值刚度相关。这是SMP5434a之前版本都存在的缺失的能量项，对MPP543

3、4a也一样。这种现象在用拉格朗日乘子(LagrangeMultiplier)方程时不会出现。与*constrained_joint_stiffness相关的能量出现在jntforc文件中，也包含在glstat文件中的弹簧和阻尼能和内能中。回想弹簧阻尼能“springanddamperenergy”，不管是从铰链刚度还是从离散单元而来，总是包含在内能里面。在MATSUM文件中能量值是按一个part一个part输出的(参见*database_matsum)。沙漏能Hourglassenergy仅当在卡片*control_energy中设置HGEN项为2时才会计算和输出。同样，刚

4、性墙能和阻尼能仅当上面的卡片中RWEN和RYLEN分别设置为2时才会计算和输出。刚性阻尼能集中到内能里面。质量阻尼能以单独的行“systemdampingenergy”出现。由于壳的体积粘性(bulkviscosity)而产生的能量耗散(energydissipated)在版本970.4748之前是不计算的。在后续子版本中，设置TYPE=-2来在能量平衡中包含它。最理想的情况下能量平衡：总能量totalenergy＝初始总能量＋外力功externalwork。换句话说，如果能量比率energyratio(指的是glstat中的totalenergy/initialenerg

5、y，实际上是totalenergy/(initialenergy+externalwork))等于1.0。注意，质量缩放而增加质量可能会导致能量比率增加。注意在LS-prepost的History>Globalenergies中不包含删掉的单元(erodedelements)的能量贡献，然而GLSTAT文件中的能量包含了它们。注意它们的贡献可以通过ASCII>glstat中的“ErodedKineticEnergy”&“ErodedInternalEnergy”来绘制。侵蚀能量(Erodedenergy)是与删掉的单元相关的内能和删掉的节点相关的动能。典型来说，如果没有单元

6、删掉“energyratiow/oerodedenergy”等于1，如果有单元被删掉则小于1。删掉的单元与“totalenergy/initialenergy”比率没有关系。总能量比率增加要归于其它原因，比如增加质量。重述一下，将一个单元删掉时，文件glstat中的内能和动能不会反映能量的丢失。取而代之的是能量的丢失记录在glstat文件的“erodedinternalenergy”&“erodedkineticenergy”中。如果用内能减去“erodedinternalenergy”将得到分析中还存在的单元的内能。对动能也一样。matsum文件中的内能和动能只包含余下(

7、noneroded)的单元的贡献。注意，如果在*control_contact卡中将ENMASS设置为2，则与删掉的单元的相关的节点不会删掉，“erodedkineticenergy”是0。在LS-prepost中History>Global只是动能和内能的简单相加，因此不包含接触能和沙漏能等的贡献。壳的负内能克服这种不真实效应的措施：a.关掉考虑壳的减薄(ISTUPDin*control_shell)b.调用壳的体积粘性(setTYPE=-2(是-1吧！)在*control_bulk_viscosity卡中)c.