Abstract
This paper establishes a three-dimensional mathematical model for a tubing/riser system subjected to external loadings. Effects such as mud weights, internal pressure in the tubing and annulus, and cable forces are all included in the model. The derivations of the mathematical model are based on principals of beam-column stability theory, which lead to a set of ordinary differential equations with the loading and displacements as functions of the axial location of the riser. With the system equations established, the rest of the paper is devoted to the devise of a numerical algorithm to solve the equations with boundary conditions imposed. The problem is re-casted as two points boundary value problem (TPVBP). Modified Quasilinearization Algorithms (MQA) together with Method of Particular Solution (MPS) are employed in the search of numerical solution for the TPBVP. Results from these algorithms are compared to those obtained from commercial finite element program to verify the validity of the mathematical model and accuracy of the numerical algorithm.