In this dissertation, the extended Melnikov's method has been applied to several nonlinear ship dynamics models, which are related to the new generation of stability criteria in the International Maritime Organization (IMO). The advantage of this extended Melnikov's method is it overcomes the limitation of small damping that is intrinsic to the implementation of the standard Melnikov's method. The extended Melnikv's method is first applied to two published roll motion models. One is a simple roll model with nonlinear damping and cubic restoring moment. The other is a model with a biased restoring moment. Numerical simulations are investigated for both models. The effectiveness and accuracy of the extended Melnikov's method is demonstrated. Then this method is used to predict more accurately the threshold of global surf-riding for a ship operating in steep following seas. A reference ITTC ship is used here by way of example and the result is compared to that obtained from previously published standard analysis as well as numerical simulations. Because the primary drawback of the extended Melnikov's method is the inability to arrive at a closed form equation, a 'best fit'approximation is given for the extended Melnikov numerically predicted result. The extended Melnikov's method for slowly varying system is applied to a roll-heave-sway coupled ship model. The Melnikov's functions are calculated based on a fishing boat model. And the results are compared with those from standard Melnikov's method. This work is a preliminary research on the application of Melnikov's method to multi-degree-of-freedom ship dynamics. In the last part of the dissertation, the method of manufactured solution is applied to systems with chaotic behavior. The purpose is to identify points with potential numerical discrepancies, and to improve computational efficiency. The numerical discrepancies may be due to the selection of error tolerances, precisions, etc. Two classical chaotic models and two ship capsize models are examined. The current approach overlaps entrainment in chaotic control theory. Here entrainment means two dynamical systems have the same period, phase and amplitude. The convergent region from control theory is used to give a rough guideline on identifying numerical discrepancies for the classical chaotic models. The effectiveness of this method in improving computational efficiency is demonstrated for the ship capsize models. / Ph. D.