LS-DYNA系列_Mie-Gruneisen状态方程

引言

计算机仿真计算中，对于真实世界中的力学过程的模拟，一直是研究热点之一。在科研、工程应用方面都有着重要的意义。但计算机专业毕业的学生，由于缺乏力学专业知识，因此在制作仿真程序时对基础概念不清楚。本系列内容旨在帮助专业计算机人员迅速了解力学知识，为国产商业力学仿真软件提供一份力量。
本文主要讲述针对固体热动力学行为提出的Mie-Gruneisen状态方程。

历史发展

Mie-Gruneisen状态方程由Grunesien提出，可以准确地描述金属材料在高温、高压、高应变率条件下表现出来的材料动态行为。该方程早期是大量固体状态金属的一个精确热动力学描述方程，之后的改进形式则可分别描述气体和固体爆炸物、真实气体和高压固体等的本构关系。

状态方程描述

金属材料Grunesien状态方程的形式为：
$$p-p_H=\rho \gamma (E-E_H)$$
式中，$p_H$和$E_H$分别为沿固体Hugoniot曲线的压力和内能，$p$和$E$为受压缩固体材料的压力和内能，$\rho$为材料密度，$\gamma$为Gruneisen系数，可通过以下公式来估计：
$$\rho \gamma ={\rho }_0{\gamma }_0$$
式中，${\rho }_0$、${\gamma }_0$分别为绝对零度时的密度和Gruneisen常数。基于实验数据，冲击波速度$U_s$和粒子速度$U_p$之间的关系为：
$$U_s=c_0+s_1{U}_p$$
式中，$c$和$s$为与材料有关的系数。

冲击波前后的质量方程为：
$$\rho (U_s-U_p)={\rho }_0{U}_s$$
由该式可得$U_s-U_p=\frac{{\rho }_0}{\rho }{U}_s$，即$U_p=(1-\frac{{\rho }_0}{\rho })U_s$。令$\eta =1-\frac{{\rho }_0}{\rho }$，有$U_p=\eta U_s$。将其代入$U_s-U_p$关系式（即$U_s=c_0+s_1{U}_p$）中，可得$U_s=\frac{c_0}{1-s_1\eta }$和$U_p=\frac{\eta c_0}{1-s_1\eta }$。
冲击波前后的动量方程为：
$$p_H-p_0={\rho }_0{U}_s{U}_p$$
式中，$p_0$为冲击波前压强，此处认为$p_0=0$，将$U_s=\frac{c_0}{1-s_1\eta }$和$U_p=\frac{\eta c_0}{1-s_1\eta }$代入该式，可得出Hugoniot压力为：
$$p_H=\frac{{\rho }_0{c}_0^2\eta }{(1-s_1\eta {)}^2}$$
式中，$\eta =1-\frac{{\rho }_0}{\rho }$。
而且，从Hugoniot关系出发，可以得到内能$E_H$：
$$E_H=\frac{1}{2}\frac{p_H\mu }{\rho }$$
式中，$\mu =\frac{{\rho }_0}{\rho }-1$。