The aim of this paper is to develop a simulation model of large deformation problems following a semi-analytical method, incorporating the complications of geometric and material non-linearity in the formulation. The solution algorithm is based on the method of energy principle in structural mechanics, as applicable for conservative systems. A one-dimensional solid circular bar problem has been solved in post-elastic range assuming linear elastic, linear strain hardening material behaviour. Type of loading includes uniform uniaxial loading and gravity loading due to body force, whereas the geometry of the bar is considered to be non-uniformly taper. Results are validated successfully with benchmark solution and some new results have also been reported. The location of initiation of elasto-plastic front and its growth are found to be functions of geometry of the bar and loading conditions. Some indicative results have been presented for static and dynamic problems and the solution methodology developed for one-dimension has been extended to the elasto-plastic analysis of two-dimensional strain field problems of a rotating disk.