统一物态方程参数计算的研究.pdfVIP

64 第七届全国冲击动力学讨论会文集 统一物态方程参数计算研究 赵 峰 中国工程物理研究院流体物理研究所 冲击波物理与爆轰物理国防科技重点实验室 四川绵阳 621900 依据由 Vinet 、Parsafar 和 Mason （简称P-M ）等人提出的两种统一物态方程， 本文对以其为等温压缩线的完全状态方程进行了研究。根据对热力学状态量的分 析，从热压系数和定容比热入手，对完全物态方程的建立进行了讨论，给出了在 已知物质等温线条件下构成物态方程的方法，并推导了 Vinet 和 P-M 两种表达式 的完全物态方程形式。以几种材料为例，通过实验测量的冲击绝热线和有关参数， 采用最小二乘法，对其物态方程进行了拟合。 1.概述 在涉及动高压研究工作中，使用最多 的当然是以冲击绝热线为参考态的 Guneison 物态方程，该方程将压力和能量分为冷和热两个部分，该方程也有受到 限制的地方，如存在极限压缩度、要求 Guneison 系数只能是常数或比容的函数， [1] 而研究结果 认为 Gruneison 系数与温度和比容都相关，是温度的缓变函数，对比 容的变化依赖更大。 Macdonald[2]在1969年的全面的评论中指出，完全不太可能存在具有简明而对 所有凝聚物质有效的统一状态方程。尽管存在这一似乎暗淡的预言，Vinet 及其合 作者[3-5]仍然在1986年提出了统一固体状态方程（以下简称 Vinet 物态方程）。统一 状态方程的基础是固体材料结合能与相互分子间距离存在统一的关系，根据这个 关系状态方程可以利用如下热力学关系 ∂p ∂E ⎡ ⎤ ⎡ ⎤ p T − （1） ⎢ ⎥ ⎢ ⎥ ∂T ∂v ⎣ ⎦ ⎣ ⎦ v T 计算，其中 p 是压力，T 是温度，v 是摩尔比容，E 是摩尔内能。Vinet 及其合作者 ∂ 略去T ⎛⎜ p ∂v⎞⎟项，即热压项，并利用结合能的体积偏导数近似获得(∂E ∂v) ，即内 ⎝ ⎠T T 压，导出的状态方程如下： 第七届全国冲击动力学讨论会文集 65 [ 2 ] [ ( )] p 3B (1−x) / x expη 1−x （2） 0 1/ 3 ⎛ ⎞ ⎛∂p ⎞ v 其中x ⎜ ⎟ ，等温体积压缩模量B0 −v⎜ ⎟ ，和体积模量对压力偏导数 v0 ⎝ ∂v⎠ ⎝ ⎠ T ⎡⎛ ⎞ ⎤ 3 ∂B η ⎢⎜ ⎟−1 ⎥，由此获