Model order reduction is an effective method to generate macromodels for system-level simulation. But it is difficult to deal with the electro-mechanical-damping coupling problems. So we presents a new approach to model the capacitive microaccelerometers with squeeze film damping and electrostatic effects using model order reduction (MOR) method. In this approach, the mechanical, squeeze film damping and electrostatic domains of the devices are modeled separately and then coupled at system-level. The macromodel for squeezed film damping effects could account for slip condition of the flow at low pressure and edge effects. In addition, some important parameters are preserved as symbol. The extracted macromodels are translated into the hardware description language and imported into a circuit simulator. An accelerometer is used to demonstrate the feasibility and efficiency of the proposed approach. Numerical simulation results show that the extracted macromodel can dramatically reduce the computation cost while capturing the device behavior accurately.