In order to compute the force transmitted through rubber bushings with a commercial multibody dynamic analysis program, a rubber bushing model is needed. The default rubber bushing model of multibody dynamic analysis program such as ADAMS or DADS is the Voigt model which is a parallel spring-damper system. The damping characteristic of the Voigt model is proportional to the frequency. However, the frequency-dependent characteristic of rubber bushing is not well represented using the Voigt model. Alternatively, the viscoelastic characteristic of rubber bushings can be denoted by complex stiffness. It is the purpose of this paper to develop a bushing model which can represent more realistic bushing behavior for multibody dynamic analysis program, especially ADAMS. The bushing is implemented using the ADAMS Transfer Function, Single Input, Single-Output feature. The proposed method is applied to the rubber bushing of the lower control arm of a rear suspension of 1994 Ford Taurus. and the transmitted force through rubber bushing model are analyzed and compared with the Voigt model