Numerical analysis of highly deformable elastoplastic beams

Abstract

The objective of the present study is to develop a numerical formula-tion to predict the behavior of highly deformable elastoplastic thin beams. Following the Euler-Bernoulli bending, the axial and shear effects are neglected, and the nonlinear second-order differential equation regarding the angle of rotation is defined based on the spe-cific moment-curvature relationship. Although the formulation can be used for general materials, three constitutive models are employed: linear-elastic, bilinear elastoplastic, and linear-elastic with Swift isotropic hardening. The resultant boundary value problem is solved by means of the fourth-order Runge-Kutta integration procedure and the one-parameter nonlinear shooting method. The performance of the present formulation is investigated via three numerical problems involving finite bending of slender beams composed of elastoplastic materials. For these problems, numerical solutions regarding rota-tions, displacements and strains for the loading, unloading and re-loading phases are provided. Finally, it is shown that the present methodology can also be used to determine the post-buckling behav-ior of elastoplastic thin beams.