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 (C44) and shear elastic modulus (Cp=(C11-C12)/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 C44 and Cp 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 (utot) and P, B, C44, Cp 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.