In vehicle multi-body dynamics (MBD) modeling, the stiffness parameterization of leaf spring (LS) is an unavoidable challenge regardless of the selection of modeling methods. On the contrary, the parameterization of Coil Spring (CS) stiffness is easy to achieve by adjusting the scale factor. Therefore, a novel LS stiffness parameterization method by treating the suspension stiffness as an intermediate variable through a CS stiffness is proposed based on the virtual displacement theory. The proposed method is then implemented in the vehicle-level modeling of a commercial Van with front transverse leaf spring suspensions and rear longitudinal parabolic leaf spring suspensions. The MBD model is validated by natural frequency tests and suspension stiffness simulations. Furthermore, the vertical acceleration of the car body is also verified. Results show that the root mean square (RMS) values of body vertical acceleration in the equivalent CS model are just slightly lower than that in the LS suspensions. The applicability and capability of the proposed method are proven to address the limitation of LS stiffness parameterization in MBD modeling. It lays the groundwork for efficiently simulating the LS suspension in vehicle ride and handling design and optimization.