Effect of boundary conditions on shakedown analysis of heterogeneous materials

The determination of the ultimate load-bearing capacity of structures made of elastoplastic heterogeneous materials under varying loads is of great importance for engineering analysis and design. Therefore, it is necessary to accurately predict the shakedown domains of these materials. The static shakedown theorem, also known as Melan's theorem, is a fundamental method used to predict the shakedown domains of structures and materials. Within this method, a key aspect lies in the construction and application of an appropriate self-equilibrium stress field(SSF). In the structural shakedown analysis, the SSF is typically constructed by governing equations that satisfy no external force(NEF) boundary conditions. However, we discover that directly applying these governing equations is not suitable for the shakedown analysis of heterogeneous materials. Researchers must consider the requirements imposed by the Hill-Mandel condition for boundary conditions and the physical significance of representative volume elements(RVEs). This paper addresses this issue and demonstrates that the sizes of SSFs vary under different boundary conditions, such as uniform displacement boundary conditions(DBCs), uniform traction boundary conditions(TBCs), and periodic boundary conditions(PBCs). As a result, significant discrepancies arise in the predicted shakedown domain sizes of heterogeneous materials. Built on the demonstrated relationship between SSFs under different boundary conditions, this study explores the conservative relationships among different shakedown domains, and provides proof of the relationship between the elastic limit(EL) factors and the shakedown loading factors under the loading domain of two load vertices. By utilizing numerical examples, we highlight the conservatism present in certain results reported in the existing literature. Among the investigated boundary conditions, the obtained shakedown domain is the most conservative under TBCs.Conversely, utilizing PBCs to construct an SSF for the shakedown analysis leads to less conservative lower bounds, indicating that PBCs should be employed as the preferred boundary conditions for the shakedown analysis of heterogeneous materials.

A SOLUTION PROCEDURE FOR LOWER BOUND LIMIT AND SHAKEDOWN ANALYSIS BY SGBEM

The symmetric Galerkin boundary element method (SGBEM) instead ofthe finite element method is used t perform lower bound limit andshakedown analysis of structures. The self-equilibrium stress fieldsare constructed by a linear combination of several basicself-equilibrium stress fields with parameters to be determined.These basic self-equilibrium stress fields are expressed as elasticresponses of the body to im- posed permanent strains and obtainedthrough elastic-plastic incremental analysis.

SHAKEDOWN ANALYSIS OF SHELL STRUCTURES OF KINEMATIC HARDENING MATERIALS

It is of great practical importance to analyze the shakedown of shell structures under cyclic loading, especially of those made of strain hardening materials.In this paper, same further understanding of the shakedown theorem for kinematic hardening materials has been made, and it is applied to analyze the shakedown of shell structures Though the residual stress of a real stale is related to plastic strain, the time-independent residual stress field as we will show in the theorem may be unrelated to the time-independent kinematically admissible plastic strain field For the engineering application, it will lie much more convenient to point this out clearly and definitely, otherwise it will be very difficult. Also, we have proposed a new method of proving this theorem.The above theorem is applied to the shakedown analysis of a cylindrical shell with hemispherical ends. According to the elastic solution, various possible residual sfcss and plastic strain Jlelds, the shakedown analysis of the structure can be reduced to a mathematical programming problem.The results of calculation show that the shakedown load oj strain hardening materials is about 30-40% higher than that of ideal plastic materials. So it is very important to consider the hardening of materials in the shakedown analysis,for it can greatly increase the structure design capacity, and meanwhile provide ascicntific basis to improve the design of shell structures.

Kinematic Shakedown Analysis for Strain-Hardening Plates with the C1 Nodal Natural Element Method

This paper proposes a novel numerical solution approach for the kinematic shakedown analysis of strain-hardening thin plates using the C^(1)nodal natural element method(C^(1)nodal NEM).Based on Koiter’s theorem and the von Mises and two-surface yield criteria,a nonlinear mathematical programming formulation is constructed for the kinematic shakedown analysis of strain-hardening thin plates,and the C^(1)nodal NEM is adopted for discretization.Additionally,König’s theory is used to deal with time integration by treating the generalized plastic strain increment at each load vertex.A direct iterative method is developed to linearize and solve this formulation by modifying the relevant objective function and equality constraints at each iteration.Kinematic shakedown load factors are directly calculated in a monotonically converging manner.Numerical examples validate the accuracy and convergence of the developed method and illustrate the influences of limited and unlimited strain-hardening models on the kinematic shakedown load factors of thin square and circular plates.

Upper Bound Shakedown Analysis of Plates Utilizing the C^(1) Natural Element Method

This paper proposes a numerical solution method for upper bound shakedown analysis of perfectly elasto-plastic thin plates by employing the C^(1) natural element method.Based on the Koiter’s theorem and von Mises yield criterion,the nonlinear mathematical programming formulation for upper bound shakedown analysis of thin plates is established.In this formulation,the trail function of residual displacement increment is approximated by using the C^(1) shape functions,the plastic incompressibility condition is satisfied by introducing a constant matrix in the objective function,and the time integration is resolved by using the Konig’s technique.Meanwhile,the objective function is linearized by distinguishing the non-plastic integral points from the plastic integral points and revising the objective function and associated equality constraints at each iteration.Finally,the upper bound shakedown load multipliers of thin plates are obtained by direct iterative and monotone convergence processes.Several benchmark examples verify the good precision and fast convergence of this proposed method.

Shakedown Analysis of 3-D Structures Using the Boundary Element Method Based on the Static Theorem

The static shakedown theorem was reformulated for the boundary element method (BEM) rather than the finite element method with Melan's theorem, then used to develop a numerical solution procedure for shakedown analysis. The self-equilibrium stress field was constructed by a linear combination of several basis self-equilibrium stress fields with undetermined parameters. These basis self-equilibrium stress fields were expressed as elastic responses of the body to imposed permanent strains obtained using a 3-D BEM elastic-plastic incremental analysis. The lower bound for the shakedown load was obtained from a series of nonlinear mathematical programming problems solved using the Complex method. Numerical examples verified the precision of the present method.

A Shakedown Strength Based Parametric Optimization Technique and Its Application on an Airtight Module

In aerospace engineering,design and optimization of mechanical structures are usually performed with respect to elastic limit.Besides causing insufcient use of the material,such design concept fails to meet the ever growing needs of the light weight design.To remedy this problem,in the present study,a shakedown theory based numerical approach for performing parametric optimization is presented.Within this approach,strength of the structure is measured by its shakedown limit calculated from the direct method.The numerical method developed for the structural optimization consists of nested loops:the inner loop adopts the interior point method to solve shakedown problems pertained to fxed design parameters,while the outer loop employs the genetic algorithm to fnd optimal design parameters leading to the greatest shakedown limit.The method established is frst verifed by the classic plate-with-a-circular-hole example,and after that it is applied to an airtight module for determining few key design parameters.By carefully analyzing results generated during the optimization process,it is convinced that the approach can become a viable means for designing similar aerospace structures.

Shakedown and Limit Analysis of Defective Pressure Vessels

This paper presents a failure analysis of pressure vessels with defects by a direct method of limit and shakedown analysis. The defects considered are part through slots with various geometric configurations. The engineering situation considered here has practical importance in the pressure vessel industry. The results are compared with those obtained by a step by step procedure using the professional code ABAQUS and where possible, with those provided by semiempirical formulae used by industry. The limit and shakedown analysis methods are found to be more economical and more reliable than marching solutions achieved by step by step elastic plastic analysis. The effects of various part through slots on the load carrying capacities of pressure vessels are investigated.