With the expectation of hull girder asymmetry and corresponding shift in elastic neutral axis resulting from collision damages and other forms of structural deteriorations, the interaction of vertical and horizontal hull girder capacities become quite significant in the assessment of ship structural safety. This paper therefore extends the application of a previously proposed interactive-numerical probabilistic based methodology for structural safety to assess the hull girder ultimate strength reliability of a damaged ship by means of a user-defined numerical framework. Hull girder capacity is calculated using the NS94D ultimate strength code, which is based on the Smith’s progressive collapse method. The resulting deterministic responses have been interactively linked to the NESSUS probabilistic framework so that the reliability of the damaged hull girder is predicted using an implicit limit state function defined based on a transformation of coordinates to appropriately account for any shift in the neutral axis. Random deviations of the constituent variables are directly applied to calculate the ultimate strength deterministic responses, thereby circumventing the need to characterize any correlated strength variable, which is at best subjective. The conventional approach of characterizing ultimate strength by an assumed coefficient of variation and distribution type was found to be conservative in predicting structural safety of ships relative to the proposed method. Application of the interactive-numerical technique for structural reliability is therefore considered significant for problems involving correlated random variables with unknown statistical characteristics. The method is being considered to predict the safety of cracked hull girders by accounting for the residual strength and further load bearing capabilities of deteriorated and adjacent elements.