有初始几何缺陷的壳体弹塑性大变形分析

本文在等温小变形弹塑性内时本构方程偏量形式的基础上，导出了适用于大位移、大转动、小应变分析的弹塑性内时本构方程，进一步推导出了带有初始几何缺陷的非线性弹塑性问题的有限元方程，可用于分析缺陷对结构非线性弹塑性反应的影响，也可用于带缺陷的非线性问题求解及稳定性分析．     

Elastoplastic Analysis with Limited Plastic Deformations and Displacements*

ABSTRACT

This paper presents solution methods for elastoplastic and shakedown analysis of linearly elastic, perfectly plastic bodies for which the conventional classical formulations of these problems are completed by constraints on overall plastic deformation and elastoplastic displacement. The methods are described in terms of nonlinear mathematical programming and provide solutions when the plastic reserves of the body are not fully exhausted, and the plastic performance and the plastic deformations are controlled. Application of the method is illustrated by an example.     

The General Mathematical Theory of Plasticity and the Il’yushin Postulates of Macroscopic Definability and Isotropy

The physical laws characterizing the relation between stresses and strains are considered and analyzed in the general modern theory of elastoplastic deformations and in its postulates of macroscopic definability and isotropy for initially isotropic continuous media. The fundamentals of this theory in continuum mechanics were developed by A.A. Il’yushin in the mid-twentieth century. His theory of small elastoplastic deformations under simple loading became a generalization of Hencky’s deformation theory of flow, whereas his theory of elastoplastic processes which are close to simple loading became a generalization of the Saint-Venant–Mises flow theory to the case of hardening media. In these theories, the concepts of simple arid complex loading processes arid the concept of directing form change tensors are introduced; the Bridgman law of volume elastic change and the universal Roche–Eichinger laws of a single hardening curve under simple loading are adopted; and the Odquist hardening for plastic deformations is generalized to the case of elastoplastic hardening media for the processes of almost simple loading without consideration of a specific history of deformations for the trajectories with small arid mean curvatures. In this paper we discuss the possibility of using the isotropy postulate to estimate the effect of forming parameters in the stress-strain state appeared due to the strain-induced anisotropy during the change of the internal structures of materials. We also discuss the possibility of representing the second-rank symmetric stress and strain tensors in the form of vectors in the linear coordinate six-dimensional Euclidean space. An identity principle is proposed for tensors and vectors.     

Thermoelastoplastic Deformation of a Multilayer Ball

The problem of centrally symmetric deformation of a multilayer elastoplastic ball in the process of successive accretion of preheated layers to its outer surface is considered in the framework of small elastoplastic deformations. The problems of residual stress formation in the elastoplastic ball with an inclusion and a cavity are solved under various mechanical boundary conditions on the inner surface and for prescribed thermal compression distributions. The graphs of residual stress and displacement fields are constructed.     

Finite elastoplastic deformation of an internally pressurized hollow sphere

This paper considers large elastoplastic deformations of an internally pressurized hollow sphere of dilatant soil. A complete analytical solution for the expansion of a hollow sphere is developed. The soil is modelled as an elastic-perfectly plastic material obeying the Mohr-Coulomb yield criterion. A non-associated plastic flow rule is used and therefore the dilation of the material is fully taken into account. Closed form solutions are obtained for the stresses and the elastic-plastic deformations of arbitrary magnitude when a hollow sphere of soil is subjected to constant external pressure and monotonically increasing internal pressure. A selection of numerical results is presented to indicate the effects of various key parameters     

Stress resultants plasticity with general closest point projection

In this paper, general closest point projection algorithm is derived for the elastoplastic behavior of a cross-section of a beam finite element. For given section deformations, the section forces (stress resultants) and the section tangent stiffness matrix are obtained as the response for the cross-section. Backward Euler time integration rule is used for the solution of the nonlinear evolution equations. The solution yields the general closest projection algorithm for stress resultants plasticity model. Algorithmic consistent tangent stiffness matrix for the section is derived. Numerical verification of the algorithms in a mixed formulation beam finite element proves the accuracy and robustness of the approach in simulating nonlinear behavior.     

A FINITE ELEMENT—MATHEMATICAL PROGRAMMING METHOD FOR ELASTOPLASTIC PROBLEMS BASED ON THE PRINCIPLE OF VIRTUAL WORK

By expanding the yielding function according to Taylor series and neglecting the high order terms, the elastoplastic constitutive equation is written in a linear complementary form. Based on this linear complementary form and the principle of virtual work, a finite element-complementary method is derived for elastoplastic problem. This method is available for materials which satisfy either associated or nonassociated flow rule. In addition, the existence and uniqueness of solution for the method are also discussed and some useful conclusions are given. The project is supported by the National Natural Science Foundation of China     

Formation of the residual stress field under local thermal actions

The one-dimensional process of material deformation due to local heating and subsequent cooling is analyzed in the framework of the classical theory of elastoplastic deformations. The problem of formation of residual stresses in a thin plate made of an elastoplastic material under a given thermal action is solved. The graphs of fields of residual stresses and displacements are constructed.     

An endochronic plasticity formulation for filled rubber

The nature of elastomeric material demands the consideration of finite deformations, nonlinear elasticity including damage as well as rate-dependent and rate-independent dissipative properties. While many models accounting for these effects have been refined over time to do better justice to the real behavior of rubber-like materials, the realistic simulation of the elastoplastic characteristics for filled rubber remains challenging.The classical elastic-ideal-plastic formulation exhibits a distinct yield-surface, whereas the elastoplastic material behavior of filled rubber components shows a yield-surface free plasticity. In order to describe this elastoplastic deformation of a material point adequately, a physically based endochronic plasticity model was developed and implemented into a Finite Element code. The formulation of the ground state elastic characteristics is based on Arruda and Boyce (1993) eight-chain model. The evolution of the constitutive equations for the nonlinear endochronic elastoplastic response are derived in analogy to the Bergström–Boyce finite viscoelasticity model discussed by Dal and Kaliske (2009).     

Deformation and Heating of an Elastoviscoplastic Cylindrical Layer Moving Owing to a Varying Pressure Drop

Within the framework of the theory of large elastoplastic deformations generalized to the case of viscous and thermophysical properties of materials, we give a solution of a sequence of coupled problems on the onset and development of a flow in a material layer filling the gap between two rigid coaxial cylindrical surfaces under increasing pressure drop and on the subsequent flow deceleration under decreasing pressure gradient. Here the thermophysical and deformation processes are coupled, and the yield stress depends on temperature. Heat production due to the layer material friction against the rough cylindrical boundary surfaces is taken for an additional heat source.     

有初始几何缺陷的壳体弹塑性大变形分析

在等温小变形弹塑性内时本构方程偏量形式的基础上，导出了适用于大位移、小应变分析的弹塑性内时本构方程。并导出了带有初始几何缺陷的非线性弹塑性问题的有限元方程。文中给出的算例表明本方法是可行有效的。     

A FINITE ELEMENT—MATHEMATICAL PROGRAMMING METHOD FOR ELASTOPLASTIC PROBLEMS BASED ON THE PRINCIPLE OF VIRTUAL WORK

By expanding the yielding function according to Tay lor series and neglecting the high order terms,the elastoplastic constitutive equation is written in a linear complementary form.Based on this linear complementary form and the principle of virtual work,a finite element-complementary method is derived for elastoplastic problem.This method is available for materials which satisfy either associated or nonassociated flow rule.In addition,the existence and uniqueness of solution for the method are also discussed and some useful conclusions are given.     

Elastoplastic Deformation of Flexible Cylindrical Shells With Two Circular Holes under Axial Tension

The elastoplastic state of thin cylindrical shells with two circular holes under axial tension is analyzed considering finite deflections. The distributions of stresses along the contours of the holes and in the zone of their concentration are studied by solving doubly nonlinear boundary-value problems. The solution obtained is compared with the solutions that account for either physical nonlinearity (plastic deformations) and geometrical nonlinearity (finite deflections) alone and with a numerical solution of the linearly elastic problem. The stress-strain state near the two holes is analyzed depending on the distance between the holes and the nonlinear factors accounted for__________Translated from Prikladnaya Mekhanika, Vol. 41, No. 5, pp. 52–57, May 2005.     

A convergent bounding principle for a class of elastoplastic strain-hardening solids

For a class of elastoplastic solids exhibiting a mixed kinematic-isotropic behavior, a convergent bounding principle is presented. This principle, formulated by means of a suitable perturbation method, is able to provide upper bounds on different kinds of the actual plastic deformation induced by a specified load history. These bounds, expressed in terms of a “fictitious solution,” can be rendered as close to the actual deformations as desired, but at the cost of increasing computational efforts. A bounding principle for repeated loads and a shakedown theorem are also presented as special cases of the convergent bounding principle.     

Physically and Geometrically Nonlinear Deformation of Spherical Shells with an Elliptic Hole

The elastoplastic state of thin spherical shells with an elliptic hole is analyzed considering that deflections are finite. The shells are made of an isotropic homogeneous material and subjected to internal pressure of given intensity. Problems are formulated and a numerical method for their solution with regard for physical and geometrical nonlinearities is proposed. The distribution of stresses (strains or displacements) along the hole boundary and in the zone of their concentration is studied. The results obtained are compared with the solutions of problems where only physical nonlinearity (plastic deformations) or geometrical nonlinearity (finite deflections) is taken into account and with the numerical solution of the linearly elastic problem. The stress—strain state in the neighborhood of an elliptic hole in a shell is analyzed with allowance for nonlinear factors __________ Translated from Prikladnaya Mekhanika, Vol. 41, No. 6, pp. 95–104, June 2005.     

Thermal Stresses in an Elastoplastic Tube Depending on the Choice of Yield Conditions

We use the solution of a one-dimensional problem of the theory of thermal stresses in an elastoplastic tube heated on its interior surface and maintained at a constant temperature on the exterior surface as an example to make a comparison of both the results and solution methods depending on the choice of each of three conventional yield criteria: piecewise linear criteria of maximum shear and maximum reduced stresses and a smooth criterion of maximum octahedral stresses. It is established that while the transition of stresses from the face of the Tresca prism to its edge (change in the flow regime) in the first of the piecewise linear yield criteria takes place at the plastic flow onset, in the second one, this transition occurs on the elastoplastic boundary. The yield stress is assumed to be temperature dependent.     

Use of the Boussinesq solution in geotechnical and road engineering: influence of plasticity

The Boussinesq solution for the distribution of stresses in a half-space resulting from surface loads is largely used in geotechnical and road engineering. It is based on the assumption of a linear–elastic homogeneous isotropic half-space for the soil media. Since the soil exhibits nonlinear and irreversible behavior, it is of major interest to study the validity of this solution for elastoplastic soils. This paper includes an investigation of this issue using finite element modeling. The study is conducted by comparing the elastic stress distribution to that obtained using elastoplastic finite element analyses. Results show that the plasticity reduces the attenuation of the vertical stresses in the soil mass, which means that the Boussinesq solution underestimates the stresses in an area which contributes to the soil settlement. To cite this article: M. Sadek, I. Shahrour, C. R. Mecanique 335 (2007).     

A Simplified Model for Calculating Hydrodynamic Lubrication Film Thickness in Elastoplastic Line Contacts

The present paper proposes a simplified model for calculating hydrodynamic lubrication film thickness in elastoplastic line contacts. According to the Saint-Venant’s principle, the pressure in the contact is taken as uniformly distributed, this gives the contact surface elastic deformations in the inlet zone far away from the contact center close to real ones while gives those close to the contact center greater than real ones. This treatment is validated for hydrodynamic lubricated elastic contacts for relatively light loads and high rolling speeds. It gives the film thickness at the contact center a little higher than that calculated based on the real elastic model. The treatment is extended to a hydrodynamic lubricated elastoplastic line contact. The contact surfaces in the inlet zone are assumed as elastic and their deformations are calculated based on the uniform pressure distribution in the elastoplastic contact area. An inlet zone analysis is taken for obtaining the calculating equation of the hydrodynamic film thickness at the contact center. The equation overestimates the central film thickness but gives a satisfactory film thickness prediction for the heavy load which gives significant plastic deformations in the elastoplastic contact. It is found that when the load is lighter than 0.6 w _pc , the contact can be taken as elastic when calculating the central film thickness, while when the load is heavier than 0.6 w _pc , the contact can be taken as fully plastic; Here w _pc is the critical load for the contact fully plastic deformation. The plastic deformation in an elastoplastic line contact is found to reduce the hydrodynamic lubrication film thickness in the contact. This reduction is greater for higher rolling speeds and heavier loads. However, it is significantly dropped with increasing surface hardness.     

Analysis of Stretching of Elastoplastic Samples and Necking with Edge Effects

Large deformations of elastoplastic cylindrical rods and shells of different thicknesses under tension are studied. The effect of sample geometry and the dependence of stresses on deformations under uniaxial tension on edge effects and necking are estimated. The applicability of the Considere criterion for determining the instance of stability loss of plastic tensile strain is analyzed. The computational experiments confirm that the role of similarity in the processes of nonuniform tensile strain of samples can be played by the ratio of the tangent of the slope of the true strain diagram to the stress intensity.     

Bending of an elastoplastic Hencky bar-chain: from discrete to nonlocal continuous beam models

The static behavior of an elastoplastic one-dimensional lattice system in bending, also called a microstructured elastoplastic beam or elastoplastic Hencky bar-chain (HBC) system, is investigated. The lattice beam is loaded by concentrated or distributed transverse monotonic forces up to the complete collapse. The phenomenon of softening localization is also included. The lattice system is composed of piecewise linear hardening–softening elastoplastic hinges connected via rigid elements. This physical system can be viewed as the generalization of the elastic HBC model to the nonlinear elastoplasticity range. This lattice problem is demonstrated to be equivalent to the finite difference formulation of a continuous elastoplastic beam in bending. Solutions to the lattice problem may be obtained from the resolution of piecewise linear difference equations. A continuous nonlocal elastoplastic theory is then built from the lattice difference equations using a continualization process. The new nonlocal elastoplastic theory associated with both a distributed nonlocal elastoplastic law coupled to a cohesive elastoplastic model depends on length scales calibrated from the spacing of the lattice model. Differential equations of the nonlocal engineering model are solved for the structural configurations investigated in the lattice problem. It is shown that the new micromechanics-based nonlocal elastoplastic beam model efficiently captures the scale effects of the elastoplastic lattice model, used as the reference. The hardening–softening localization process of the nonlocal continuous model strongly depends on the lattice spacing which controls the size of the nonlocal length scales.