为了确定bcc和fcc晶体的EAM/FS势函数参数，研究了EAM/FS势函数与弹性常数的关系。对bcc和fcc结构晶体，分别导出压强(P)和体积弹性模量(B)、弹性常数(C₄₄)和剪切弹性模量(C_p=(C₁₁-C₁₂)/2)的利用嵌入函数、对势函数和电子密度分布函数的表达式。发现C₄₄和C_p的大小不仅取决于所考虑的原子和周围原子的距离，还受近邻原子排列状况的影响。将关于结合能(u_tot)和P、B、C₄₄、C_p的5个拟合方程转化为求最小值的优化模型，给出了5种典型bcc结构晶体(V、Mo、Nb、Ta、W)和3种典型fcc结构晶体(Cu、γ-Fe、Ni)的结合能中待定参数的值。采用上述参数计算得到的最小结合能与实验结果一致，此时所对应的原子间距与晶格常数相同，表明了方法的有效性。

Potential functions are extensively applied in molecular dynamics (MD) simulation of metals. Selection of them is a very important step in MD simulations due to its effects of the precision and reliability of the simulations. They are one of the most important reference data during the process of calculation. In order to cover the shortage of pairwise potentials for modelling transition metals, EAM/FS many-body potentials have been introduced since 80's of last century. For the sake of determining parameters in the EAM/FS potential functions of bcc and fcc crystals through macro mechanical properties, relations between the EAM/FS potential functions and elastic constants were investigated in this work. Expressions of the pressure (P) and the bulk modulus (B), elastic constant (C₄₄) and shear elastic modulus (C_p=(C₁₁-C₁₂)/2) in terms of the embedding function, pair potential function and the electron density distribution function were deduced for bcc and fcc structures, respectively. It was found that the magnitude of the C₄₄ and C_p depends on the distances between the considered atom and surrounding atoms, but also the configuration of surrounding atoms. Finally, by converting five fitting equations about the cohesive energy (u_tot) and P, B, C₄₄, C_p into an optimization model of finding minimum value, the values of the six undetermined parameters in the cohesive energy were given for five typical bcc crystals (V, Mo, Nb, Ta and W) and three typical fcc crystals (Cu, γ-Fe, Ni), respectively. For each crystal, calculation errors show accuracy of parameter values. The obtained calculation results, for the minimum cohesive energy and the corresponding atomic distance, fit well with the reported experimental data, by adopting the above values of the parameters, which indicates the effectiveness for our method.