受冲板壳结构的弹塑性力学特性分析

采用基于随动坐标系的 4节点假设应变场壳单元及显式有限元方法分析受冲板壳结构的弹塑性力学特性。材料模型采用弹塑性等向强化模型 ,接触搜寻采用一体化接触搜寻方法 ,接触力由罚参数法计算 ,算例表明 :该方法简明、直观、快捷、方便     

A general high‐order finite element formulation for shells at large strains and finite rotations

For hyperelastic shells with finite rotations and large strains a p‐finite element formulation is presented accommodating general kinematic assumptions, interpolation polynomials and particularly general three‐dimensional hyperelastic constitutive laws. This goal is achieved by hierarchical, high‐order shell models. The tangent stiffness matrices for the hierarchical shell models are derived by computer algebra. Both non‐hierarchical, nodal as well as hierarchical element shape functions are admissible. Numerical experiments show the high‐order formulation to be less prone to locking effects. Copyright © 2003 John Wiley & Sons, Ltd.     

A robust non‐linear mixed hybrid quadrilateral shell element

In the paper a non‐linear quadrilateral shell element for the analysis of thin structures is presented. The variational formulation is based on a Hu–Washizu functional with independent displacement, stress and strain fields. The interpolation matrices for the mid‐surface displacements and rotations as well as for the stress resultants and strains are specified. Restrictions on the interpolation functions concerning fulfillment of the patch test and stability are derived. The developed mixed hybrid shell element possesses the correct rank and fulfills the in‐plane and bending patch test. Using Newton's method the finite element approximation of the stationary condition is iteratively solved. Our formulation can accommodate arbitrary non‐linear material models for finite deformations. In the examples we present results for isotropic plasticity at finite rotations and small strains as well as bifurcation problems and post‐buckling response. The essential feature of the new element is the robustness in the equilibrium iterations. It allows very large load steps in comparison to other element formulations. Copyright © 2005 John Wiley & Sons, Ltd.     

An Assumed Strain Triangular Solid Element for Efficient Analysis of Plates and Shells with Finite Rotation

A simple triangular solid shell element formulation is developed for efficient analysis of plates and shells undergoing finite rotations. The kinematics of the present solid shell element formulation is purely vectorial with only three translational degrees of freedom per node. Accordingly, the kinematics of deformation is free of the limitation of small angle increments, and thus the formulation allows large load increments in the analysis of finite rotation. An assumed strain field is carefully selected to alleviate the locking effect without triggering undesirable spurious kinematic modes. In addition, the curved surface of shell structures is modeled with flat facet elements to obviate the membrane locking effect. Various numerical examples demonstrate the efficiency and accuracy of the present element formulation for the analysis of plates and shells undergoing finite rotation. The present formulation is attractive in that only three points are needed for numerical integration over an element.     

A finite element formulation for piezoelectric shell structures considering geometrical and material non‐linearities

In this paper, we present a non‐linear finite element formulation for piezoelectric shell structures. Based on a mixed multi‐field variational formulation, an electro‐mechanical coupled shell element is developed considering geometrically and materially non‐linear behavior of ferroelectric ceramics. The mixed formulation includes the independent fields of displacements, electric potential, strains, electric field, stresses, and dielectric displacements. Besides the mechanical degrees of freedom, the shell counts only one electrical degree of freedom. This is the difference in the electric potential in the thickness direction of the shell. Incorporating non‐linear kinematic assumptions, structures with large deformations and stability problems can be analyzed. According to a Reissner–Mindlin theory, the shell element accounts for constant transversal shear strains. The formulation incorporates a three‐dimensional transversal isotropic material law, thus the kinematic in the thickness direction of the shell is considered. The normal zero stress condition and the normal zero dielectric displacement condition of shells are enforced by the independent resultant stress and the resultant dielectric displacement fields. Accounting for material non‐linearities, the ferroelectric hysteresis phenomena are considered using the Preisach model. As a special aspect, the formulation includes temperature‐dependent effects and thus the change of the piezoelectric material parameters due to the temperature. This enables the element to describe temperature‐dependent hysteresis curves. Copyright © 2011 John Wiley & Sons, Ltd.     

An assumed strain triangular curved solid shell element formulation for analysis of plates and shells undergoing finite rotations

A formulation for 36‐DOF assumed strain triangular solid shell element is developed for efficient analysis of plates and shells undergoing finite rotations. Higher order deformation modes described by the bubble function displacements are added to the assumed displacement field. The assumed strain field is carefully selected to alleviate locking effect. The resulting element shows little effect of membrane locking as well as shear locking, hence, it allows modelling of curved shell structures with curved elements. The kinematics of the present formulation is purely vectorial with only three translational degrees of freedom per node. Accordingly, the present element is free of small angle assumptions, and thus it allows large load increments in the geometrically non‐linear analysis. Various numerical examples demonstrate the validity and effectiveness of the present formulation. Copyright © 2001 John Wiley & Sons, Ltd.     

A reduced integration solid‐shell finite element based on the EAS and the ANS concept—Large deformation problems

In this paper we address the extension of a recently proposed reduced integration eight‐node solid‐shell finite element to large deformations. The element requires only one integration point within the shell plane and at least two integration points over the thickness. The possibility to choose arbitrarily many Gauss points over the shell thickness enables a realistic and efficient modeling of the non‐linear material behavior. Only one enhanced degree‐of‐freedom is needed to avoid volumetric and Poisson thickness locking. One key point of the formulation is the Taylor expansion of the inverse Jacobian matrix with respect to the element center leading to a very accurate modeling of arbitrary element shapes. The transverse shear and curvature thickness locking are cured by means of the assumed natural strain concept. Further crucial points are the Taylor expansion of the compatible cartesian strain with respect to the center of the element as well as the Taylor expansion of the second Piola–Kirchhoff stress tensor with respect to the normal through the center of the element. Copyright © 2010 John Wiley & Sons, Ltd.     

A triangular finite shell element based on a fully nonlinear shell formulation

This work presents a fully nonlinear six-parameter (3 displacements and 3 rotations) shell model for finite deformations together with a triangular shell finite element for the solution of the resulting static boundary value problem. Our approach defines energetically conjugated generalized cross-sectional stresses and strains, incorporating first-order shear deformations for an inextensible shell director (no thickness change). Finite rotations are treated by the Euler–Rodrigues formula in a very convenient way, and alternative parameterizations are also discussed herein. Condensation of the three-dimensional finite strain constitutive equations is performed by applying a mathematically consistent plane stress condition, which does not destroy the symmetry of the linearized weak form. The results are general and can be easily extended to inelastic shells once a stress integration scheme within a time step is at hand. A special displacement-based triangular shell element with 6 nodes is furthermore introduced. The element has a nonconforming linear rotation field and a compatible quadratic interpolation scheme for the displacements. Locking is not observed as the performance of the element is assessed by several numerical examples, which also illustrate the robustness of our formulation. We believe that the combination of reliable triangular shell elements with powerful mesh generators is an excellent tool for nonlinear finite element analysis.Fellowship funding from FAPESP (Fundação de Amparo à Pesquisa do Estado de São Paulo) and CNPq (Conselho Nacional de Pesquisa), together with the material support and stimulating discussions in IBNM (Institut für Baumechanik und Numerische Mechanik), are gratefully acknowledged in this work.     

A triangular finite element for thin plates and shells

A finite element modelling technique which utilizes a triangular element with 45 degrees-of-freedom and seven-point integration has been tested for analysis of thin plate and shell structures. The element is based on the degenerate solid shell concept and the mixed formulation with assumed independent inplane and transverse shear strains. Numerical result indicates effectiveness of the present modelling technique which features combined use of elements with kinematic modes and those without kinematic modes in an attempt to eliminate both locking and spurious kinematic modes at global structural level.     

Ductile damage analysis of elasto‐plastic shells at large inelastic strains

The objective of this contribution is to model ductile damage phenomena under consideration of large inelastic strains, to couple the corresponding constitutive law with a multi‐layer shell kinematics and to give finally an adequate finite element formulation. An elastic–plastic constitutive law is formulated by using a spatial hyperelasto‐plastic formulation based on the multiplicative decomposition of the deformation gradient. To include isotropic ductile damage the continuum damage model of Rousselier is modified so as to consider large strains and additionally extended by various void nucleation and macro‐crack criteria. In order to achieve numerical efficiency, elastic strains are supposed to be sufficiently small providing a numerical effective integration based on the backward Euler rule. Finite element formulation is enriched by means of the enhanced strain concept. Thus the well‐known deficiencies due to incompressible deformations and the inclusion of transverse strains are avoided. Several examples are given to demonstrate the performance of the algorithms developed concerning large inelastic strains of shells and ductile damage phenomena. Copyright © 2000 John Wiley & Sons, Ltd.     

On the use of an enhanced <Emphasis Type="Italic">transverse</Emphasis> shear strain shell element for problems involving large rotations

 This paper presents the extension of a previously developed formulation for shell elements in order to account for non-linear geometric effects, particularly in the presence of large rotations. To eliminate transverse shear locking, the developed shell formulation provides an enlargement of the transverse shear strain field coming from the usual degenerated concept. Doing so, additional transverse shear strain terms are included into the original displacement-based functional, following the enhanced strain approach. To reproduce the behavior of shell structures under large rotations and displacements, a rotation-free configuration is considered, where constitutive relations are stated. Dealing with finite incremental rotations, a singularity-free procedure is employed, characterizing the evolution of normal vectors to shell's mid-surface. Representative non-linear examples are considered, providing the validation of the enhanced shell element as well as the algorithmic procedures implemented, when compared to other formulations in the literature.     

Cyclic plastic behavior analysis based on the micromorphic mixed hardening plasticity model

In this paper, a small strain micromorphic elasto-plastic model with isotropic/kinematic hardening is presented for modeling the size effect and Bauschinger effect in material with microstructure. A nonlinear kinematic hardening model is embedded into the micromorphic framework by employing a backstress, a micro-backstress and a micro-couple-backstress in a physical way. The material intrinsic length scale is introduced in the constitutive law, leading to the presence of higher order stress. The present model is further implemented into a 2D plane strain finite element frame with a fully implicit stress integration scheme. The generalized consistent tangent modulus is derived to achieve the parabolic convergence of the global nodal force equilibrium equation. Two numerical examples, including a thin film and a plate with underlying structures subjected to cyclic loading, are analyzed to verify the theoretical developments and numerical formulations. Plastic behaviors in micromorphic continuum, such as size effect, Bauschinger effect, ratcheting effect and plastic shakedown phenomenon, are investigated.     

Large inelastic strain analysis by multilayer shell elements

Summary The objective of this contribution is to model large inelastic strains of ductile metals, to couple this material model with a multilayer shell kinematics and finally to achieve a finite element formulation applicable in general form to shell analysis. Elasto-plastic constitutive law is formulated by using the multiplicative decomposition of the deformation gradient and Neo-Hookean model for elastic strains assuming an overall isotropic material behavior. These 3D-material model is then enforced directly into a multilayer shell kinematics which provides a very accurate consideration of local effects, particularly stresses across the thickness. Finite element formulation is accomplished by means of the enhanced strain concept. Thus the well known deficiencies due to incompressible deformations and the inclusion of transverse strains are avoided. Several examples are given to demonstrate the performance of the algorithms developed concerning various aspects.     

Finite element formulation of the Ogden material model with application to rubber-like shells

The present contribution proposes a variational procedure for the numerical implementation of the Ogden material model. For this purpose the strain energy density originally formulated in terms of the principal stretches is transformed as variational quantities into the invariants of the right Cauchy–Green tensor. This formulation holds for arbitrary three-dimensional deformations and requires neither solving eigenvalue problems nor co-ordinate system transformations. Particular attention is given to the consideration of special cases with coinciding eigenvalues. For the analysis of rubber-like shells this material model is then coupled with a six parametric shells kinematics able to deal with large strains and finite rotations. The incompressibility condition is considered in the strain energy, but it is additionally used as 2-D constraint for the elimination of the stretching parameter at the element level. A four node isoparametric finite element is developed by interpolating the transverse shear strains according to assumed strain concept. Finally, examples are given permitting to discuss the capability of the finite element model developed concerning various aspects. © 1998 John Wiley & Sons, Ltd.     

Thin shells with finite rotations formulated in biot stresses: Theory and finite element formulation

A bending theory for thin shells undergoing finite rotations is presented, and its associated finite element model is described. The kinematic assumption is based on a shear elastic Reissner-Mindlin theory. The starting point for the derivation of the strain measures are the resultant equilibrium equations and the associated principle of virtual work. Within this formulation the polar decomposition of the shell material deformation gradient leads to symmetric strain measures. The associated work-conjugate stress resultants and stress couples are integrals of the Biot stress tensor. This tensor is invariant with respect to rigid body motions and, therefore, appropriate for the formulation of constitutive equations. Finite rotations are introduced via Eulerian angles. The finite element discretization of arbitrary shells is based on the isoparametric concept formulated with respect to a plane reference configuration. The numerical model is applied to different non-linear plate and shell problems and compared with existing formulations. Due to a consistent linearization, the step size of a load increment is only limited by the local convergence behaviour of Newton's method.     

A nine node finite element for analysis of geometrically non-linear shells

A nine node finite element is presented for the analysis of thin shell structures undergoing large deflection. The finite element formulation is based on the concept of degenerate solid shell element and the Hellinger-Reissner principle with independent strain. Three versions of assumed independent strain are selected to suppress spurious kinematic modes. One version leads to a finite element model which is kinematically stable at element level while the other two give globally stable models. Numerical tests indicate that the finite element model which is stable at element level may reveal the locking effect in certain cases. However, the other two models are free of locking.     

THE FOUR-NODE C0 SHELL ELEMENT REFORMULATED

A reformulated four-node shell element, based on the analysis of the moment redistribution mechanism development by C⁰ plate bending and shell elements, is presented. The moment redistribution mechanism of a finite shell element model is shown to be predominantly activated by the membrane flexural action of the shell. This action is triggered through the membrane strain components which participate in the moment equilibrium equations of the finite element assembly system. An equivalent elastic foundation action, along with the activation of the in-plane twisting stiffness of the shell, may also contribute to the moment redistribution mechanism of the finite shell element model. The proposed shell element formulation aims at retaining the non-spurious contribution of the transverse shear/membrane strain energy to the flexural behaviour of the shell, through the activation of the moment redistribution mechanism. Yet, any potentially spurious, whether locking or kinematic, mechanism is rejected. In warped configurations, the element activates appropriate coupling mechanisms of the bending terms to nodal translations. The so-obtained reformulated four-node shell element exhibits an excellent behaviour without experiencing any locking phenomena or zero-energy modes, while its formulation is kept simple, based on physical considerations. The proposed formulation performs equally well in flat as well as in warped shell element applications.     

Standard formulation of elastic-plastic deformation laws

Summary The present article introduces the standard formulation of the deformation law of inelastic bodies applied to large deformations, taking into account the thermomechanical problem.Detailed explanations concerning elastic-plastic material behaviour are given. The effects of different tensor rates in modelling kinematic hardening are analysed numerically under the conditions of an initially isotropic state. Final studies deal with simulating yield behaviour of sheet metals in plane principal stress.