Large viscoplastic deformations of shells. Theory and finite element formulation

The paper is concerned with large viscoplastic deformations of shells when the constitutive model is based on the concept of unified evolution equations. Specifically the model due to Bodner and Partom is modified so as to fit in the frame of multiplicative viscoplasticity. Although the decomposition of the deformation gradient in elastic and inelastic parts is employed, no use is made of the concept of the intermediate configuration. A logarithmic elastic strain measure is used. An algorithm for the evaluation of the exponential map for nonsymmetric arguments as well as a closed form of the tangent operator are given. On the side of the shell theory itself, the shell model is chosen so as to allow for the application of a three-dimensional constitutive law. The shell theory, accordingly, allows for thickness change and is characterized by seven parameters. The constitutive law is evaluated pointwise over the shell thickness to allow for general cyclic loading. An enhanced strain finite element method is given and various examples of large shell deformations including loading-unloading cycles are presented.     

An inelastic constitutive model of blood vessels

Summary The present paper is concerned with the formulation of a constitutive model describing the hysteresis loops of stress-strain relations of blood vessels under cyclic loading conditions. It is assumed that the hysteresis loop is composed of elastic deformation and viscoplastic deformation. Hence the total strain is expressed as the sum of the elastic part and the inelastic part on the basis of a finite deformation theory. Then the elastic part is established by postulating a strain energy function of an exponential type, while the inelastic part is formulated by modifying the nonlinear kinematic hardening rule in the viscoplastic model proposed by Chaboche et al. A comparison of the numerical result with the literature shows that the present model can describe the hysteresis loop qualitatively.     

板壳非线性分析的新理论新方法

提出了板壳非线性分析的新理论新方法。首先建立了下列几个新的本构关系:塑性应变向量增量与总应变向量增量的新关系,热塑性应变向量增量与总应变向量增量及温度应变向量增量的新关系,粘塑性应变向量增量与总应变向量增量的新关系,热粘塑性应变向量增量与总应变向量增量及温度应变向量增量的新关系。这些关系分别称为弹塑性应变增量理论、热弹塑性应变增量理论、弹粘塑性应变增量理论及热弹粘塑性应变增量理论,避开了屈服曲面、加载曲面、流动法则及复杂的非线性应力应变关系。其次建立了非线性样条无网格法,这种方法是以新的本构关系、几何非线性理论、变分原理、广义变分原理、加权残数法及样条离散化为基础建立的,避免了经典本构关系及有限元法带来的巨大困难及缺陷,不仅计算简便,而且精度高,收敛速度很快。建立了板壳非线性分析的统一格式,对板壳的几何非线性分析、材料非线性分析及双重非线性分析都适用。     

A facet-like shell theory

A linear theory for facet-like thin elastic shells is derived where strain/displacement, curvature change/displacement and constitutive relations appear the same as for flat plates. Application of Koiter's arguments shows that the theory is a valid first approximation. The theory is of interest for limiting cases of faceted finite element analysis of smooth shells. Although the final equations of facet-like shell theory do not have quite as simple a form as more conventional equations it is possible that their derivation from equations for flat plates may appeal to engineers. A specialization of the equations is given to circular cylindrical shells where four simple illustrative examples show no essential differences with results from more conventional theory.     

A FEM model for the active control of curved FGM shells using piezoelectric sensor/actuator layers

In this paper, a generic finite element formulation is developed for the static and dynamic control of FGM (functionally graded material) shells with piezoelectric sensor and actuator layers. The properties of the FGM shell are graded in the thickness direction according to a volume fraction power‐law distribution. The proposed finite element model is based on variational principle and linear piezoelectricity theory. A constant displacement and velocity feedback control algorithm coupling the direct and inverse piezoelectric effects is applied in a closed‐loop system to provide feedback control of the integrated FGM shell structure. Both static and dynamic control of FGM shells are simulated to demonstrate the effectiveness of the proposed active control scheme within a framework of finite element discretization and piezoelectric integration. Copyright © 2002 John Wiley & Sons, Ltd.     

Non‐linear analysis of shells with arbitrary evolving cracks using XFEM

A new formulation and numerical procedures are developed for the analysis of arbitrary crack propagation in shells using the extended finite element method. The method is valid for completely non‐linear problems. Through‐the‐thickness cracks in sandwich shells are considered. An exact shell kinematics is presented, and a new enrichment of the rotation field is proposed which satisfies the director inextensibility condition. To avoid locking, an enhanced strain formulation is proposed for the 4‐node cracked shell element. A finite strain plane stress constitutive model based on the logarithmic corotational rate is employed. A cohesive zone model is introduced which embodies the special characteristics of the shell kinematics. Stress intensity factors are calculated for selected problems and crack propagation problems are solved. Copyright © 2004 John Wiley & Sons, Ltd.     

Development of geometrically exact new shell elements based on general curvilinear co‐ordinates

In the present study first‐order shear deformable shell finite elements based on general curvilinear co‐ordinates are proposed. For the development of the present shell elements, a partial mixed variational functional with independently assumed strains is provided in order to avoid the severe locking troubles known as transverse shear and membrane lockings. Bubble functions are included in the shape function of displacement to improve the performance of the developed element. The proposed assumed strain four‐ and nine‐node elements based on the general tensor shell theory provide an efficient linkage framework for shell surface modelling and finite element analysis. In the several benchmark problems, the present shell elements with exact geometric representations demonstrate their performance compared to previously reported results. Copyright © 2002 John Wiley & Sons, Ltd.     

Large strain analysis of rubber-like materials based on a perturbed Lagrangian variational principle

A mixed finite element method is presented for the large strain analysis of rubber-like materials, which are considered to be nearly incompressible. Two types of constitutive relations are included: generalized Rivlin and Ogden's models. The finite element equations are derived on the basis of a perturbed Lagrangian variational principle from which both the displacement and pressure fields are independently approximated by appropriate shape functions. A physically meaningful pressure parameter is introduced in the expression of complementary energy. In the paper, a special effort is made to split the deformation energy into two distinct parts: isochoric and hydrostatic parts. By doing this, a quadratic convergence rate of nonlinear iterative solution is achieved, particularly for problems deformed in the large strain range. The finite element equations are specialized for a two-dimensional 9-node Lagrange element with three-term pressure parameters. Five examples are given to demonstrate the application of the proposed numerical algorithm.Research work supported by National Science Foundation under the grant number EET-8714628     

Inelastic material behavior of polymers – Experimental characterization,formulation and implementation of a material model

In this contribution an experimental procedure based on displacement controlled tensile tests at different rates of loading, relaxation experiments and deformation controlled loading and unloading processes with intermediate relaxations to experimentally characterize and classify the nonlinear, inelastic mechanical behavior of polymers is presented. These experiments provide data for a structured approach to parameter identification. In line with the experiments, a small strain uniaxial viscoplastic material model is derived, subsequently generalized to multiaxial loadings and implemented into a finite element program. The combination of the experimental procedure and the proposed material model is then used to characterize and model the mechanical behavior of the thermoplastic polypropylene. After the identification of the necessary material parameters, stress–strain curves have been computed for different uni- and multiaxial loadings and are compared with experimental results.     

A mixed finite element for the nonlinear analysis of in-plane loaded masonry walls

A mixed membrane eight-node quadrilateral finite element for the analysis of masonry walls is presented. Assuming that a nonlinear and history-dependent 2D stress-strain constitutive law is used to model masonry material, the element derivation is based on a Hu-Washizu variational statement, involving displacement, strain, and stress fields as primary variables. As the behavior of masonry structures is often characterized by strain localization phenomena, due to strain softening at material level, a discontinuous, piecewise constant interpolation of the strain field is considered at element level, to capture highly nonlinear strain spatial distributions also within finite elements. Newton's method of solution is adopted for the element state determination problem. For avoiding pathological sensitivity to the finite element mesh, a novel algorithm is proposed to perform an integral-type nonlocal regularization of the constitutive equations in the present mixed formulation. By the comparison with competing serendipity displacement-based formulation, numerical simulations prove high performances of the proposed finite element, especially when coarse meshes are adopted.     

An integral equation analysis of inelastic shells

This paper presents an integral equation approach for the analysis of deformation and stresses in inelastic shells of arbitrary shape subjected to arbitrary loading. The proposed mathematical model is completely consistent and is derived by transforming the three-dimensional equations from the Cartesian to the appropriate curvilinear coordinates of the shell. Appropriate kinematic assumptions for the dependence of the displacements on the thickness coordinate of the shell and assumptions regarding the loads at the ends are introduced consistently in the model to take advantage of the thinness of the shell. Numerical implementation and numerical results are presented for elastic and inelastic deformation of axisymmetric shells subjected to axisymmetric loading. These results are compared against exact elasticity, Love-Kirchoff model analysis of inelastic cylindrical shells and finite element solutions.     

A numerical study of hypoelastic and hyperelastic large strain viscoplastic Perzyna type models

For the case of metals with large viscoplastic strains, it is necessary to define appropriate constitutive models in order to obtain reliable results from the simulations. In this paper, two large strain viscoplastic Perzyna type models are considered. The first constitutive model has been proposed by Ponthot, and the elastic response is based on hypoelasticity. In this case, the kinematics of the constitutive model is based on the additive decomposition of the rate deformation tensor. The second constitutive model has been proposed by García Garino et al., and the elasticresponse is based on hyperelasticity. In this case, the kinematics of the constitutive model is based on the multiplicativedecomposition of the deformation gradient tensor. In both cases, the resultant numerical models have been implemented in updated Lagrangian formulation. In this work, global and local numerical results of the mechanical response of both constitutive models are analyzed and discussed. To this end, numerical experiments are performed and different parameters of the constitutive models are tested in order to study the sensitivity of the resultantalgorithms. In particular, the evolution of the reaction forces, the effective plastic strain, the deformed shapes and the sensitivity of the numerical results to the finite element mesh discretization have been compared and analyzed. The obtained results show that both models have a very good agreement and represent very well the characteristic of the viscoplastic constitutive model.     

带有沙漏控制的相对自由度壳元

针对一点积分的八节点相对自由度壳单元存在的沙漏现象，提出采用拟应变法解决该问题的方法，并对锁死问题进行研究。给出了带有沙漏控制的八节点相对自由度壳元内的坐标、位移插值公式，推导了拟应变的表达式，通过Hu-Washizu变分原理，建立了有限元求解方程。利用Wilson非协调位移模式，单元的计算精度得到了明显改善。算例表明：基于八节点相对自由度壳单元，本文给出的沙漏控制算法能够有效的解决线性静力问题，并且具有较高的计算精度。     

A quasi-conforming triangular laminated composite shell element based on a refined first-order theory

A quasi-conforming triangular laminated shell element based on a refined first-order shear deformation theory is presented. The Hu-Washizu variational principle, involving strain and displacement fields as variables, with stresses being considered as Lagrange multipliers, is used to develop the laminate composite shell element. Both strains and displacements are discretized in the element, while displacements alone are discretized at the boundary. The inter-element C ¹ continuity is satisfied a posteriori in a weak form. Due to the importance of rotations and shear deformation in the geometrically non-linear analyses of shells, 7 degrees of freedom per node are chosen, viz. three displacements, two first-derivatives in the in-plane directions of the out-of-plane displacement, and two transverse shear strains at each node. To consider the effect of transverse shear deformation on the global behavior of the laminated composite shell, the Reissner-Mindlin first-order theory, with shear correction factors of Chow and Whitney, is adopted. The transverse shear stresses are obtained through the integration of the 3-D equilibrium equations; and the warping induced by transverse shear is considered in the calculation of the in-plane stresses to improve their accuracy. Numerical examples show that the element has good convergence properties and leads to highly accurate stresses.     

Scaling of cylindrical shells under axial impact

A technique of correcting scaled strain rate sensitive structures subject to dynamic loads is described and applied to a shell under axial impact. The prototype and models responses were simulated using finite element method. Final deformed shape, collapse mode, displacement at the top of the shell and maximum force of scaled models are compared to the respective prototype. It is shown that the models response are quite different from the prototype if no correction is performed. By correcting the initial impact velocity, the models scale quite well and it is even possible to obtain the transition phenomenon between global and progressive buckling of shells subject to the axial impact of a mass.     

基于几何非线性不协调元变分原理的圆柱壳非线性初始稳定性分析

本文根据几何非线性不协调元增量变分原理,按严格的壳体方程,建立了高精度的圆柱壳几何非线性20参数矩形精化不协调元RCSR4,并用于圆柱壳非线性初始稳定性分析。计算结果表明,该方法收敛性良好。     

复合材料四结点四边形多层退化壳单元

本文提出了一个复合材料四结点四边形多层退化壳单元。单元从修正的Hu-Washizu变分原理导出,独立假设位移场、每层内部应变场和应力场。文中给出了几个复合材料板和壳静力线性问题实例。计算结果表明本文的单元是准确和有效的。