LYDYNA能量平衡.doc

LYDYNA能量平衡 GLSTAT(参见*database_glstat)文件中报告的总能量是下面几种能量的和： 内能 internal energy 动能 kinetic energy 接触(滑移)能 contact(sliding) energy 沙漏能 houglass energy 系统阻尼能 system damping energy 刚性墙能量 rigidwall energy GLSTAT 中报告的弹簧阻尼能”Spring and damper energy”是离散单元(discrete elements)、安全带单元 (seatbelt elements)内能及和铰链刚度相关的内能(*constrained_joint_stiffness…)之和。而内能”Internal Energy”包含弹簧阻尼能”Spring and damper energy”和所有其它单元的内能。 因此弹簧阻尼能”Spring and damper energy”是内能”Internal energy”的子集。 由SMP 5434a 版输出到glstat 文件中的铰链内能”joint internal energy”跟*constrained_joing_stiffness 不相 关。它似乎与*constrained_joint_revolute(_spherical,etc)的罚值刚度相关连。这是SMP 5434a 之前版本都存在的 缺失的能量项，对MPP 5434a 也一样。这种现象在用拉格朗日乘子(Lagrange Multiplier)方程时不会出现。 与*constrained_joint_stiffness 相关的能量出现在jntforc 文件中，也包含在glstat 文件中的弹簧和阻尼能和 内能中。回想弹簧阻尼能”spring and damper energy”，不管是从铰链刚度还是从离散单元而来，总是包含在内 能里面。 在MATSUM 文件中能量值是按一个part 一个part 的输出的(参见*database_matsum)。 沙漏能Hourglass energy 仅当在卡片*control_energy 中设置HGEN 项为2 时才计算和输出。同样，刚性墙 能和阻尼能仅当上面的卡片中RWEN 和RYLEN 分别设置为2 时才会计算和输出。刚性阻尼能集中到内能里 面。质量阻尼能以单独的行”system damping energy”出现。由于壳的体积粘性(bulk viscosity)而产生的能量耗散 (energy dissipated)在版本970.4748 之前是不计算的。在后续子版本中，设置TYPE=-2 来在能量平衡中包含它。 最理想的情况下能量平衡： 总能量total energy ＝ 初始总能量 ＋ 外力功external work 换句话说，如果能量比率energy ratio(指的是glstat 中的total energy/initial energy，实际上是total energy/(initial energy + external work)) 等于1.0。注意，质量缩放而增加质量可能会导致能量比率增加。 注意在LSprepost 的HistoryGlobal energies 中不包含删掉的单元(eroded elements)的能量贡献，然而 GLSTAT 文件中的能量包含了它们。注意它们的贡献可以通过ASCIIglstat 中的”Eroded Kinetic Energy” “Eroded Internal Energy”来绘制。侵蚀能量(Eroded energy)是与删掉的单元相关的内能和删掉的节点相关的动 能。 典型来说，如果没有单元删掉”energy ratio w/o eroded energy”等于1，如果有单元被删掉则小于1。删掉 的单元与”total energy/initial energy”比率没有关系。总能量比率增加要归于其它原因，比如增加质量。重述一 下，将一个单元删掉时，文件glstat 中的内能和动能不会反映能量的丢失。取而代之的是能量的丢失记录在 glstat 文件的”eroded internal energy” “eroded kinetic energy”中。 如果用内能减去”eroded internal energy”将得 到分析中还存在的单元的内能。对动能也一样。 matsum 文件中的内能和动能只包含余下(noneroded)的单元的贡献。 注意，如果在*control_contact 卡中将ENMASS 设置为2，则与删掉的单元的相关的节点不会删掉，”eroded kinetic energy”是0。 在LSprepost 中His