This paper has shown a numerical motion simulation method which can be employed to study on parametric rolling of ships in a seaway. The method takes account of the main nonlinear terms in the rolling equation which stabilize parametric rolling, including the nonlinear shape of the righting arm curve, nonlinear damping and cross coupling among all 6 degrees of freedom. For the heave, pitch, sway and yaw motions, the method uses response amplitude operators determined by means of the strip method, whereas the roll and surge motions of the ship are simulated, using nonlinear motion equations coupled with the other 4 degrees of freedom. For computing righting arms in seaways, Grim's effective wave concept is used. Using these transfer functions of effective wave together with the heave and pitch transfer functions, the mean ship immersion, its trim and the effective regular wave height are computed for every time step during the simulation. The righting arm is interpolated from tables, computed before starting the simulation, depending on these three quantities and the heel angle. The nonlinear damping moment and the effect of bilge keels are also taken into account. The numerical simulation tool has shown to be able to model the basic mechanism of parametric rolling motions. Some main characteristics of parametric rolling of ships in a seaway can be good reproduced by means of the method. Comprehensive parametric analyses on parametric rolling amplitude in regular waves have been carried out, with that the complicated parametric rolling phenomena can be understood better.