几何非线性复合材料层合固体壳单元

通过定义广义应力,提出了一个改进的刚度矩阵,以克服固体壳元的厚度自锁问题,并能保证沿复合材料层合结构厚度方向上的连续应力分布;将应力插值函数分为低阶和高阶两部分,建议了一个新的非线性变分泛函,推导了一个用于几何非线性分析的九节点固体壳单元,该单元的计算精度和效率基本上与九节点减缩积分单元相当,与同类型其他单元相比,该单元显著提高了计算效率。     

Nonlinear buckling analysis of hygrothermoelastic composite shell panels using finite element method

Hygrothermal stresses due to the change in environmental condition may induce buckling and dynamic instability in the composite shell structures. In the present investigation, the hygrothermoelastic buckling behavior of laminated composite shells are numerically simulated using geometrically nonlinear finite element method. The orthogonal curvilinear coordinate is used for modeling a general doubly curved deep or shallow shell surface. The geometrically nonlinear finite element formulation is based on general nonlinear strain–displacement relations in the orthogonal curvilinear coordinate system. The present theory can be applicable to thin and moderately thick shells. The mechanical linear and nonlinear stiffnesses, and the nonmechanical nonlinear geometric stiffness matrices and the hygrothermal load vector are presented. It is also observed that during the present numerical solution of nonlinear equilibrium equation, in order to construct the nonlinear stiffness matrices for the first load step, the initial deformation can be assumed as zero or any computer generated small random number or the properly scaled fundamental buckling mode shape. To verify the present formulations and finite element code, the present results are compared well with those available in the open literature. Parametric studies such as thickness ratio and shallowness ratio on buckling are performed for spherical, truncated conical and cylindrical composite shell panels. The buckling behavior and deflection shapes are characterized by multiple wrinkles along unreinforced direction at higher moisture concentrations or temperature rise.     

基于共旋三角形厚薄通用壳元的几何非线性分析

基于共旋列式方法发展了一种用于复合材料层合板结构几何非线性分析的简单高效的三结点三角形平板壳元。该壳元由具有面内转动自由度的广义协调膜元GT9与假设剪切应变场和假设单元转角场的广义协调厚薄通用板元TMT组合而成。为避免薄膜闭锁而采用单点积分计算与薄膜应变有关的项, 同时增加一个稳定化矩阵以消除单点积分导致的零能模式。基于层合板一阶剪切变形理论, 给出了考虑层合板具体铺层顺序的修正的横向剪切刚度, 使该壳元可用于中厚层合板结构的分析。由于共旋列式大转动小应变的假设, 共旋列式内核的几何线性的单元刚阵可仅计算一次而保存下来用于整个几何非线性求解的过程以提高计算效率。数值算例表明提出的壳元进行包括复合材料层合板结构的厚薄壳结构的几何非线性分析的精度高且效率高。     

Dynamic behavior of delaminated plates considering progressive failure process

The paper is to study on dynamic response behavior of the delaminated composite plates considering progressive failure process. A formula of element stiffness and mass matrices for the composite laminates is deduced by using the first-order shear deformation theory combined with the selecting numerical integration scheme. A damping model is constituted by a generalized orthogonal damping model on basis of Adams' strain energy method. A virtual interface linear spring element is also employing for avoiding the overlap and penetration phenomenon between the upper and lower sublaminates at the delamination region. The failure analysis method for the delaminated plates under dynamic loading is established by a modified Newmark direct integral method in conjunction with Tsai's failure criterion and corresponding stiffness degradation scheme. By some numerical examples, the effects of frequency of dynamic load, delamination length and location, and reduction of structure stiffness during the progressive failure process upon dynamic behavior of the delaminated composite plates are discussed. The method and conclusions would be useful for composite structures designers.     

An efficient through-thickness integration scheme in an unlimited layer doubly curved isoparametric composite shell element

This paper presents an efficient numerical integration scheme for evaluating the matrices (stiffness, mass, stress-stiffness and thermal load) for a doubly curved, multilayered, composite, quadrilateral shell finite element. The element formulation is based on three-dimensional continuum mechanics theory and it is applicable to the analysis of thin and moderately thick composite shells. The conventional formulation requires a 2 × 2 × 2 or 2 × 2 × 1 Gauss integration per layer for the calculation of element matrices. This method becomes uneconomical when a large number of layers is used owing to an excessive amount of computations. The present formulation is based on explicit separation of the thickness variable from the shell surface parallel variables. With the through-thickness variables separated, they are combined with the thickness dependent material properties and integrated separately. The element matrices are computed using the integrated material matrices and only a 2 × 2 spatial Gauss integration scheme. The response results using the present formulation are identical to those obtained using the conventional formulation. For a small number of layers, the present method requires slightly more CPU time. However, for a larger number of layers, numerical data are presented to demonstrate that the present formulation is an order-of-magnitude economical compared to the conventional scheme.     

基于改进弧长法的层压复合壳后屈曲反应分析

研究了层压复合壳在横向均匀外压作用下的后屈曲反应。应用全拉格朗日公式描述和9结点退化三维壳单元。提出了一种渐进破坏的模式,并且引入改进的弧长法用于非线性有限元分析。重点研究了不同铺设顺序及方向的层压复合壳后屈曲变形形态和破坏过程。数值计算结果的精度和稳定性得以证实。     

Global format for energy–momentum based time integration in nonlinear dynamics

A global format is developed for momentum and energy consistent time integration of second‐order dynamic systems with general nonlinear stiffness. The algorithm is formulated by integrating the state‐space equations of motion over the time increment. The internal force is first represented in fourth‐order form consisting of the end‐point mean value plus a term containing the stiffness matrix increment. This form gives energy conservation for systems with internal energy as a quartic function of the displacement components. This representation is then extended to general energy conservation via a discrete gradient representation. The present procedure works directly with the internal force and the stiffness matrix at the time integration interval end‐points, and in contrast to previous energy‐conserving algorithms, it does not require any special form of the energy function nor use of mean value products at the element level or explicit use of a geometric stiffness matrix. An optional monotonic algorithmic damping, increasing with response frequency, is developed in terms of a single damping parameter. In the solution procedure, the velocity is eliminated and the nonlinear iterations are based on the displacement components alone. The procedure represents an energy consistent alternative to available collocation methods, with an equally simple implementation. Copyright © 2014 John Wiley & Sons, Ltd.     

复合材料加筋圆拱壳的失稳特性分析

本文应用增量形式的拉格朗日列式法对其有纵横加筋的迭层圆拱壳在均布载荷作用下的稳定性进行了非线性有限元分析。文中应用Sander 壳体理论及横向剪切的影响, 推导了矩形壳元及与该壳元变形相协调的直梁元和曲梁元的切线刚度矩阵。编制了FORTRAN 计算程序。计算并分析了加筋拱壳的局部及整体失稳过程。      

动荷载作用下含损伤复合材料层合板承载能力

研究了含分层损伤层合板的动力响应和承载能力。基于层合板的一阶剪切理论,采用分项等参插值方法推导了复合材料层合板刚度阵、质量阵列式,在瑞利阻尼的基础上构造了相应的阻尼阵列式;建立了用于含分层损伤复合材料层合板动力分析的分层模型和虚拟界面联接单元,以防止低阶模态中在分层处出现的上、下子板不合理的脱离和嵌入现象;同时又采用Tsai提出的0.44刚度退化准则和动力分析的Newmark法,对含分层损伤复合材料层合板结构进行了在动荷载作用下的破坏和承载能力分析;通过典型算例,分别讨论了外载频率,分层位置,以及刚度退化对含损伤复合材料动力响应特征和承载能力的影响。本文中提出的方法和得到的结论对复合材料工程设计具有参考价值。      

复合材料加筋层合板的极限强度分析

基于更新拉格朗日格式,应用非线性层合三维退化壳元,结合有效的复合材料失效准则、刚度退化模型以及提出的刚度矩阵奇异判断准则,对复合材料加筋层合板在轴向载荷作用下的压缩极限强度问题进行深入研究。讨论了铺层方式、板厚等对极限强度的影响。通过与试验结果进行比较,表明基于提出的刚度矩阵奇异判断准则,结合增量更新拉格朗日格式下非线性层合三维退化有限元的计算方法,能有效计算复合材料加筋层合板的轴向压缩极限强度,并具有很高的精度。     

A computational model for the finite element analysis of thermoplasticity coupled with ductile damage at finite strains

An updated Lagrangian implicit FEM model for the analysis of large thermo‐mechanically coupled hyperelastic‐viscoplastic deformations of isotropic porous materials is considered. An appropriate framework for constitutive modelling is introduced that includes a stress‐free thermally expanded configuration and a plastically deformed unstressed damaged configuration. A two‐level iterative scheme is employed at each time increment to solve the field equations governing the conservation of momentum (mechanical step) and the conservation of energy (thermal step) for the coupled thermo‐mechanical problem. Exact linearizations for the calculation of the tangent stiffness are performed in each of these solution steps. A fully implicit, thermo‐mechanically coupled and incrementally objective Euler‐backward radial return based map is developed for the time integration of the constitutive equations. The present model is used to analyse a number of benchmark examples including metal forming processes wherein temperature and the accumulated damage play an important role in influencing the deformation mechanism and the nature of the deformed workpiece. Copyright © 1999 John Wiley & Sons, Ltd.     

An exact conserving algorithm for nonlinear dynamics with rotational DOFs and general hyperelasticity. Part 2: shells

Following the approach developed for rods in Part 1 of this paper (Pimenta et al. in Comput. Mech. 42:715–732, 2008), this work presents a fully conserving algorithm for the integration of the equations of motion in nonlinear shell dynamics. We begin with a re-parameterization of the rotation field in terms of the so-called Rodrigues rotation vector, allowing for an extremely simple update of the rotational variables within the scheme. The weak form is constructed via non-orthogonal projection, the time-collocation of which ensures exact conservation of momentum and total energy in the absence of external forces. Appealing is the fact that general hyperelastic materials (and not only materials with quadratic potentials) are permitted in a totally consistent way. Spatial discretization is performed using the finite element method and the robust performance of the scheme is demonstrated by means of numerical examples.     

Formulation of implicit finite element methods for multiplicative finite deformation plasticity

Some constitutive and computational aspects of finite deformation plasticity are discussed. Attention is restricted to multiplicative theories of plasticity, in which the deformation gradients are assumed to be decomposable into elastic and plastic terms. It is shown by way of consistent linearization of momentum balance that geometric terms arise which are associated with the motion of the intermediate configuration and which in general render the tangent operator non-symmetric even for associated plastic flow. Both explicit (i.e. no equilibrium iteration) and implicit finite element formulations are considered. An assumed strain formulation is used to accommodate the near-incompressibility associated with fully developed isochoric plastic flow. As an example of explicit integration, the rate tangent modulus method is reviewed in some detail. An implicit scheme is derived for which the consistent tangents, resulting in quadratic convergence of the equilibrium iterations, can be written out in closed form for arbitrary material models. All the geometrical terms associated with the motion of the intermediate configuration and the treatment of incompressibility are given explicitly. Examples of application to void growth and coalescence and to crack tip blunting are developed which illustrate the performance of the implicit method.     

Substructuring in the implicit simulation of single point incremental sheet forming

This paper presents a direct substructuring method to reduce the computing time of implicit simulations of single point incremental forming (SPIF). Substructuring is used to divide the finite element (FE) mesh into several non-overlapping parts. Based on the hypothesis that plastic deformation is localized, the substructures are categorized into two groups: the plastic—nonlinear—substructures and the elastic—pseudo-linear—substructures. The plastic substructures assemble a part of the FE mesh that is in contact with the forming tool; they are iteratively updated respecting all nonlinearities. The elastic substructures model the elastic deformation of the rest of the FE mesh. For these substructures, the geometrical and the material behaviour are assumed linear within the increment. The stiffness matrices and the internal force vectors are calculated at the beginning of each increment then they are statically condensed to eliminate the internal degrees of freedom (DOF). In the iteration process the condensed stiffness matrices for the elastic substructures are kept constant. The condensed internal force vectors are updated by the multiplication of the condensed stiffness matrices and the displacement increments. After convergence, any geometrical and material nonlinearity for the elastic substructures are nonlinearly updated. The categorization of substructures in plastic and elastic domains is adapted during the simulation to capture the tool motion. The resulting, plastic and condensed elastic, set of equations is solved on a single processor. In an example with 1600 shell elements, the presented substructuring of the SPIF implicit simulation is 2.4 times faster than the classical implicit simulation.     

A finite element post-processed Galerkin method for dimensional reduction in the non-linear dynamics of solids. Applications to shells

In this paper we introduce the finite element version of the so-called post-processed Galerkin method into the field of solid mechanics and apply the new technique to the dynamics of shells. The proposed post-processed method provides low-cost means to lift low-dimensional solutions to high-dimensional solutions. It is the very fact that the kinematical fields are improved to higher orders which makes the method of great interest. Our shell theory is geometrically exact in the sense that all non-linearities are included in the formulation. For time integration an energy/momentum scheme is used to enhance integration stability. Two hierarchical enhanced finite elements are formulated, on the basis of which a specific post-processed method is then developed. With the help of some examples of non-linear shell vibrations, a critical examination and validation of the post-processed method is carried out.     

弹性膜结构几何非线性分析的刚体准则法

提出了一种新型弹性空间膜结构几何非线性分析方法。根据刚体准则的思想，初始受力平衡的单元在经历刚体位移后，其单元结点力方向随单元发生转动而大小不变，单元仍保持平衡。建立了新型三角形空间膜单元，该膜单元由三根空间杆件组成铰接三角形，并在中间张拉薄膜而成，杆件的材料与薄膜相同。所建立的空间膜单元的整体位形由杆单元空间铰接三角形确定，而膜单元的有限弹性变形由内部张拉的薄膜变形确定。由满足刚体准则的杆单元几何刚度矩阵推导了空间膜单元的几何刚度矩阵。根据刚体准则思想，认为膜单元在变形过程中，其刚体位移对其整体变形的贡献较大，而单元的弹性变形贡献较小。采用更新的拉格朗日格式的增量迭代法，在分析的每个阶段充分考虑刚体转动效应，利用小变形线性化理论处理自然变形的剩余效应。该方法几何刚度矩阵推导简单，无需引入对单元大变形的人为假定，可容易地退化为平面膜单元，增量迭代计算过程充分考虑刚体准则，对若干典型空间膜结构算例的分析及与已有方法的比较，验证了所建单元与方法的准确性以及计算效率。     

Nonlinear viscoelastic analysis of orthotropic beams using a general third-order theory

The displacement based finite element model of a general third-order beam theory is developed to study the quasi-static behavior of viscoelastic rectangular orthotropic beams. The mechanical properties are considered to be linear viscoelastic in nature with a scope to undergo von Kármán nonlinear geometric deformations. A differential constitutive law is developed for an orthotropic linear viscoelastic beam under the assumptions of plane-stress. The fully discretized finite element equations are obtained by approximating the convolution integrals using a trapezoidal rule. A two-point recurrence scheme is developed that necessitates storage of data from the previous time step only, and not from the entire deformation history. Full integration is used to evaluate all the stiffness terms using spectral/hp lagrange polynomials. The Newton iterative scheme is employed to enhance the rate of convergence of the nonlinear finite element equations. Numerical examples are presented to study the viscoelastic phenomena like creep, cyclic creep and recovery for thick and thin beams using classical mechanical analogues like generalized n-parameter Kelvin-Voigt solids and Maxwell solids.     

Explicit thickness integration for three-dimensional shell elements applied to non-linear analysis

The problem of multilayered degenerated 3-D shell elements for which the numerical integration is performed for each ply is that of the high generation time in non-linear analysis when the number of plies is important. But these elements give accurate results for thin and moderately thick shells, so in order to reduce the generation time explicit thickness integration is investigated. We first write an expansion of the strain-displacement matrix in power series of the thickness variable in order to obtain explicit expressions of the tangent stiffness matrix and internal force vector, appearing in the non-linear formulation. Explicit expressions of non-linear stiffness matrices are presented, using the explicit integration-first approximation. Simple expressions of several matrices, sub-matrices and vectors appearing in the formulation are given here in order to obtain an important computing-time gain. Next, some numerical validation tests comparing the classical element with numerical thickness integration and this one are discussed to prove validity of this formulation.