Numerical integration of linear boundary value problems in solid mechanics by segmentation method

Abstract

Numerical integration of the system of governing equations which define a boundary value problem written down in the form of a coupled system of first-order ordinary differential equations is shown to be a powerful technique. After presenting the basic approach the paper critically examines the numerical schemes available for situations when the boundary value problem so defined has boundary layer characteristics. One such method which is originally due to Goldberg, Setlur and Alspaugh³ is described in detail, with documentation in the form of a flow diagram and a FORTRAN listing of a working subroutine. The method is shown to be computationally efficient and reliable for the solution of a class of problems in the field of solid mechanics. Potential use of the method for the solution of magnetostatic problems is indicated.