THE PROBLEMS OF LIMIT LOAD ANALYSIS AND OPTIMIZATION USING EQUILIBRIUM FINITE ELEMENTS
Authors: Kalanta Stanislovas 
Published: March 2019

The equilibrium dual discrete mathematical models of the problems of limit load analysis and optimization are investigated in the article. These models are presented in terms of static and kinematic formulation using equilibrium finite elements. In these mathematical models the possible discontinuities of displacement velocities are evaluated and the velocity of energy dissipation is estimated not only within the volume of finite elements, but at the plastic surfaces between elements.

At first, on the basis of the energy principle of the maximum external power [1,2] the general mathematical models (3) and (7) of static formulation of limit load analysis and optimization problems are created. In these models the yield conditions are controlled not only within the volume, but also at the surfaces of finite elements. The equilibrium finite elements and interpolation functions of strains (9) are used for discretization of these models. The constancy of external power is taken as the optimum criterion. The discrete expressions of fundamental relationships—equilibrium and geometric equations, yield conditions (10)-(12) for finite element and (14)-(17) for the discrete model of a body are developed. The discrete expressions of yield conditions are given using the classic collocation methods: collocation at the point, collocation at the sphere (element) and Bubnov-Galiorkin's collocation method [11,12]. The equilibrium equations of discrete structure are developed on the basis of virtual displacement principle while geometrical equations are derived using virtual force principle.

In contrast to the approach of other authors, yield conditions and geometrical equations are described not only within finite elements, but also at the surfaces between elements. That helps to design the dual discrete mathematical models of the problems (21)-(26), in which the discontinuities of displacement velocities and the velocity of energy dissipation in the place of those discontinuities are estimated. The mathematical models (22), (24), (26), (29) and (32) of kinematic problem formulation are developed from sensible static formulations by Lagrange's multiplier method. The modified mathematical models (27)-(29) are presented. In these models the equilibrium equations are eliminated or the geometrical equations are transformed into compatible equations of plastic stress velocities, in this way decreasing the number of equations and unknown values.

The dependence of the numerical results (limit load) of the frame on the approximation degree of bending moments, as well as on the discretization method of yield conditions are illustrated. In table 1 the values of limit loading parameter F 0 and their error of calculation ΔF 0 (per cent, in comparison with the analytic solution F 0 = 0,4662M 0) are presented. They are given for the first and second order finite elements with linear and parabolic distribution of bending moments using different discrete yield conditions and a different number of finite elements. The numerical result shows, that the discretization of yield conditions by Bubnov-Galiorkin's method gives the best accuracy and stable solutions. By discretizing yield conditions using the point's collocation and collocation at the element, the accurancy of numerical results depends not only on the number of elements, but also on a more or less successful choice of finite elements net.