Accuracy of a composite implicit time integration scheme for structural dynamics

A comprehensive study of the two sub‐steps composite implicit time integration scheme for the structural dynamics is presented in this paper. A framework is proposed for the convergence accuracy analysis of the generalized composite scheme. The local truncation errors of the acceleration, velocity, and displacement are evaluated in a rigorous procedure. The presented and proved accuracy condition enables the displacement, velocity, and acceleration achieving second‐order accuracy simultaneously, which avoids the drawback that the acceleration accuracy may not reach second order. The different influences of numerical frequencies and time step on the accuracy of displacement, velocity, and acceleration are clarified. The numerical dissipation and dispersion and the initial magnitude errors are investigated physically, which measure the errors from the algorithmic amplification matrix's eigenvalues and eigenvectors, respectively. The load and physically undamped/damped cases are naturally accounted. An optimal algorithm‐Bathe composite method (Bathe and Baig, 2005; Bathe, 2007; Bathe and Noh, 2012) is revealed with unconditional stability, no overshooting in displacement, velocity, and acceleration, and excellent performance compared with many other algorithms. The proposed framework also can be used for accuracy analysis and design of other multi‐sub‐steps composite schemes and single‐step methods under physical damping and/or loading. Copyright © 2016 John Wiley & Sons, Ltd.     

Geometrically nonlinear analysis of shells by quadrilateral flat shell element with drill,shear, and warping

In this paper, a four‐node quadrilateral flat shell element is proposed for geometrically nonlinear analysis based on updated Lagrangian formulation with the co‐rotational kinematics concept. The flat shell element combines the membrane element with drilling degrees of freedom and the plate element with shear deformation. By means of these linearized elements, a simplified nonlinear analysis procedure allowing for warping of the flat shell element and large rotation is proposed. The tangent stiffness matrix and the internal force recovery are formulated in this paper. Several classic benchmark examples are presented to validate the accuracy and efficiency of the proposed new and more proficient element for practical engineering analysis of shell structures. Copyright © 2016 John Wiley & Sons, Ltd.     

Adding transverse normal stresses to layered shell finite elements for the analysis of composite structures

This work presents a formulation developed to add capabilities for representing the through thickness distribution of the transverse normal stresses, σ_z, in first and higher order shear deformable shell elements within a finite element (FE) scheme. The formulation is developed within a displacement based shear deformation shell theory. Using the differential equilibrium equations for two-dimensional elasticity and the interlayer stress and strain continuity requirements, special treatment is developed for the transverse normal stresses, which are thus represented by a continuous piecewise cubic function. The implementation of this formulation requires only C⁰ continuity of the displacement functions regardless of whether it is added to a first or a higher order shell element. This makes the transverse normal stress treatment applicable to the most popular bilinear isoparametric 4-noded quadrilateral shell elements.

To assess the performance of the present approach it is included in the formulation of a recently developed third order shear deformable shell finite element. The element is added to the element library of the general nonlinear explicit dynamic FE code DYNA3D. Some illustrative problems are solved and results are presented and compared to other theoretical and numerical results.     

Variational sensitivity analysis of a nonlinear solid shell element

The paper is concerned with variational sensitivity analysis of a nonlinear solid shell element, which is based on the Hu–Washizu variational principle. The sensitivity information is derived on the continuous level and discretized to yield the analytical expressions on the computational level. Especially, the pseudo load matrix and the sensitivity matrix, which dominate design sensitivity analysis of shape optimization problems, are derived. Because of the mixed formulation, condensation of the pseudo load matrix on the element level is performed to compute the sensitivity matrix. An illustrative example from the field of geometry‐based shape optimization demonstrates the possible application of the presented formulation. Copyright © 2013 John Wiley & Sons, Ltd.     

复合材料双曲率壳屈曲和后屈曲的非线性有限元研究

基于 ABAQUS软件分析平台 , 采用非线性有限元法研究了横向载荷作用下复合材料双曲率壳的屈曲和后屈曲行为。通过在有限元模型中引入 Tsai2Wu失效准则 , 预测了复合材料双曲率壳的初始失效及渐进破坏过程 , 数值结果和试验数据吻合较好 , 表明了该模型的合理有效性 , 并详细讨论了各种参数对屈曲和后屈曲行为的影响。经分析复合材料双曲率加筋壳在均布压载和剪力联合作用下的屈曲和后屈曲行为 , 得到了屈曲载荷的拟合曲线 , 研究表明顺剪力的存在有利于提高屈曲载荷。     

Optimal solid shell element for large deformable composite structures with piezoelectric layers and active vibration control

In this paper, we present an optimal low‐order accurate piezoelectric solid‐shell element formulation to model active composite shell structures that can undergo large deformation and large overall motion. This element has only displacement and electric degrees of freedom (dofs), with no rotational dofs, and an optimal number of enhancing assumed strain (EAS) parameters to pass the patch tests (both membrane and out‐of‐plane bending). The combination of the present optimal piezoelectric solid‐shell element and the optimal solid‐shell element previously developed allows for efficient and accurate analyses of large deformable composite multilayer shell structures with piezoelectric layers. To make the 3‐D analysis of active composite shells containing discrete piezoelectric sensors and actuators even more efficient, the composite solid‐shell element is further developed here. Based on the mixed Fraeijs de Veubeke–Hu–Washizu (FHW) variational principle, the in‐plane and out‐of‐plane bending behaviours are improved via a new and efficient enhancement of the strain tensor. Shear‐locking and curvature thickness locking are resolved effectively by using the assumed natural strain (ANS) method. We also present an optimal‐control design for vibration suppression of a large deformable structure based on the general finite element approach. The linear‐quadratic regulator control scheme with output feedback is used as a control law on the basis of the state space model of the system. Numerical examples involving static analyses and dynamic analyses of active shell structures having a large range of element aspect ratios are presented. Active vibration control of a composite multilayer shell with distributed piezoelectric sensors and actuators is performed to test the present element and the control design procedure. Copyright © 2005 John Wiley & Sons, Ltd.     

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.     

几何非线性高性能复合材料筋混凝土梁Heterosis组合壳单元

对于高性能碳纤维增强聚合物复合材料(CFRP)筋混凝土梁,研究几何非线性组合壳单元模型,对预应力CFRP筋混凝土梁进行了全过程分析。引入Von Karman理论,推导了局部坐标系下Piola2Kirchhoff 应力矩阵和几何刚度矩阵;分别采用组合壳单元和分层壳单元模拟预应力CFRP 筋和玻璃纤维增强聚合物复合材料(GFRP)筋,并推导了CFRP筋对组合壳单元刚度矩阵的贡献,同时采用Heterosis选择积分技术以避免剪切锁定和零能量模式,研制了相应的非线性计算程序。计算结果与试验数据对比可知,挠度发展规律和预应力CFRP筋应变发展规律均吻合良好,说明了研究单元的有效性及研制程序的正确性;CFRP筋具有高强度性能,梁试件破坏时CFRP筋均未失效;利用预应力CFRP筋应变重分布系数研究了梁的刚度退化规律,表明采用GFRP筋代替普通钢筋在加载后期会使梁的刚度退化减小。      

Efficient implicit finite element analysis of sheet forming processes

The computation time for implicit finite element analyses tends to increase disproportionally with increasing problem size. This is due to the repeated solution of linear sets of equations, if direct solvers are used. By using iterative linear equation solvers the total analysis time can be reduced for large systems. For plate or shell element models, however, the condition of the matrix is so ill that iterative solvers do not reach the huge time‐savings that are realized with solid elements. By introducing inertial effects into the implicit finite element code the condition number can be improved and iterative solvers perform much better. An additional advantage is that the inertial effects stabilize the Newton–Raphson iterations. This also applies to quasi‐static processes, for which the inertial effects finally do not affect the results. The presented method can readily be implemented in existing implicit finite element codes. Industrial size deep drawing simulations are executed to investigate the performance of the recommended strategy. It is concluded that the computation time is decreased by a factor of 5 to 10. Copyright © 2003 John Wiley & Sons, Ltd.     

A higher-order facet quadrilateral composite shell element

A C⁰ finite element formulation of flat faceted element based on a higher-order displacement model is presented for the analysis of general, thin-to-thick, fibre reinforced composite laminated plates and shells. This theory incorporates a realistic non-linear variation of displacements through the shell thickness, and eliminates the use of shear correction coefficients. The discrete element chosen is a nine-noded quadrilateral with five and nine degrees of freedom per node. A comparison of results is also made with the 2-D thin classical and 3-D exact analytical results, and finite element solutions with 9-noded first-order element. © 1997 John Wiley & Sons, Ltd.     

Nonlinear hygro-thermo-elastic vibration analysis of doubly curved composite shell panel using finite element micromechanical model

Nonlinear free vibration behavior of laminated composite curved panel under hygrothermal environment is investigated in this article. The mathematical model of the laminated panel is developed using Green–Lagrange-type geometrical nonlinearity in the framework of higher-order mid-plane kinematics. The corrugated composite properties are evaluated through the micromechanical model and all the nonlinear higher-order terms are included in the present model for the sake of generality. The equation of vibrated panel is obtained using Hamilton's principle and discretized with the help of the finite element steps. The solutions are computed numerically using the direct iterative method. The effect of parameters on the nonlinear vibration responses is examined thoroughly by solving the wide variety of numerical examples.     

Structural analysis of composite laminates using a mixed hybrid shell element

The structural analysis of thin composite structures requires robust and effective shell elements. In this paper the variational formulation is based on a Hu–Washizu functional with independent displacements, stress resultants and shell strains. For the independent shell strains an additional interpolation part is introduced. This yields an improved convergence behaviour especially for laminated shells with coupled membrane and bending stiffness. The developed mixed hybrid shell element possesses the correct rank and fulfills the in–plane and bending patch test. The formulation is tested by several nonlinear examples including bifurcation and post–buckling response. The essential feature of the new element is the robustness in nonlinear computations with large rigid body motions. It allows very large load steps in comparison to standard displacement models.     

An isogeometric solid‐like shell element for nonlinear analysis

An isogeometric solid‐like shell formulation is proposed in which B‐spline basis functions are used to construct the mid‐surface of the shell. In combination with a linear Lagrange shape function in the thickness direction, this yields a complete three‐dimensional representation of the shell. The proposed shell element is implemented in a standard finite element code using Bézier extraction. The formulation is verified using different benchmark tests. Copyright © 2013 John Wiley & Sons, Ltd.     

A quadrilateral shell element incorporating thickness-stretch for nearly incompressible hyperelasticity

The authors proposed a quadrilateral shell element enriched with degrees of freedom to represent thickness-stretch. The quadrilateral shell element can be utilized to consider large deformations for nearly incompressible materials, and its performance is demonstrated in small and large deformation analyses of hyperelastic materials in this study. Formulation of the proposed shell element is based on extension of the MITC4 shell element. A displacement variation in the thickness direction is introduced to evaluate the change in thickness. In the proposed approach, the thickness direction is defined using the director vectors at each midsurface node. The thickness-stretch is approximated by the movements of additional nodes, which are placed along the thickness direction from the bottom to the top surface. The transverse normal strain is calculated using these additional nodes without assuming the plane stress condition; hence, a three-dimensional constitutive equation can be employed without any modification. In this work, the authors apply an assumed strain technique to the special shell element to alleviate volumetric locking for nearly incompressible materials. Several numerical examples are presented to examine the effectiveness of the proposed element.     

A methodology for the formulation of error estimators for time integration in linear solid and structural dynamics

In this article, we present a novel methodology for the formulation of a posteriori error estimators applicable to time‐stepping algorithms of the type commonly employed in solid and structural mechanics. The estimators constructed with the presented methodology are accurate and can be implemented very efficiently. More importantly, they provide reliable error estimations even in non‐smooth problems where many standard estimators fail to capture the order of magnitude of the error. The proposed methodology is applied, as an illustrative example, to construct an error estimator for the Newmark method. Numerical examples of its performance and comparison with existing error estimators are presented. These examples verify the good accuracy and robustness predicted by the analysis. Copyright © 2005 John Wiley & Sons, Ltd.     

复合材料圆柱壳中心受横向集中载荷作用的渐进破坏非线性有限元分析

旨在建立能够正确预计复合材料圆柱壳的屈曲和后屈曲渐进破坏行为的模拟策略。采用有限元方法和Hashin失效准则进行模拟,基于该失效准则编写了用户材料子程序,然后插入到商用有限元软件ABAQUS中。分析了中心受横向集中载荷作用复合材料圆柱壳板,壳板的2条直边弹性支持,2条曲边自由。为了探讨弹性边界条件和集中载荷作用点应力集中的影响,将有限元分析结果与文献中的试验结果进行了比较,提出了一种合理的弹性边界选取依据。研究结果表明,在建模中考虑了弹性边界和集中载荷作用点处存在的应力集中后,本文中模拟的结果与文献中的试验结果比较接近,模拟精度明显高于文献中报道的结果。这也验证了本文中建立的模拟策略的合理性。     

Discrete equivalent time integration methods for transient analysis

It was recently shown in a series of papers that the frequency response of temporally discretized finite element equations and, consequently, the achievable accuracy cannot be manipulated independently in different frequency ranges. In addition, there exist limitations on the achievable accuracy of a time integration method no matter what the order of accuracy of the method is. Motivated by this fact, a family of time integration methods is derived in the time domain based on the principle that the exact solution of the semidiscrete equation of a system and the solution of the time integration method match at discrete time steps. It is necessary to pursue an exact match at discrete time steps, i.e. discrete equivalence, since the solutions of semidiscrete equations are obtained only at the time steps. Two time integration methods, that are exact at the time steps, are obtained based on the impulse and ramp response invariance principles. Numerical examples are presented to show the advantage of the proposed methods and to compare the performance of them with the performance of some popular methods. Copyright © 2003 John Wiley & Sons, Ltd.     

基于ANSYS的复合材料螺旋桨叶片有限元建模与分析

针对船用复合材料螺旋桨叶片,分别采用壳单元赋予不同截面参数和实体单元赋予不同等效材料参数两种建模方法进行了有限元分析,对比了壳单元和实体单元模型的模态、均布压力载荷下的响应以及均匀温度场中的热变形。结果表明:对于从桨根到桨尖各位置铺层不同的叶片来说,两种建模方法在建模时间和铺层调整便捷性方面各有利弊;模态分析中两种模型的固有频率和振型计算结果存在一定差异,它们均可在一定条件下用于复合材料螺旋桨叶片模态计算;两种模型中,模拟均布压力载荷下的响应时应采用壳单元模型,而分析叶片在均匀温度场中结构的热变形应采用实体单元模型。     

A hybrid stress nine-node degenerated shell element for geometric nonlinear analysis

A stabilized nine-node degenerated shell element, previously derived for linear analysis by the admissible matrix formulation, is extended to geometrically nonlinear analysis. The assumed stress modes pertinent to stabilization matrix are contravariant in nature and their energy products with the displacement-derived covariant strain can be programmed without resorting to numerical integration and the Gram-Schmidt orthogonalization. Numerical tests show that its accuracy is virtually identical to the uniformly reduced integration element.     

Sandwich shell finite element for dynamic explicit analysis

The contribution of this paper consists of new development of transverse shear stresses through the thickness and finding an expression for the critical time step for explicit time integration of layered shells. This work presents the finite element (FE) formulation and implementation of a higher‐order shear deformable shell element for dynamic explicit analysis of composite and sandwich shells. The formulation is developed using a displacement‐based third‐order shear deformation shell theory. Using the differential equilibrium equations and the interlayer requirements, special treatment is developed for the transverse shear, resulting in a continuous, piecewise quartic distribution of the transverse shear stresses through the shell thickness. Expressions are developed for the critical time step of the explicit time integration for orthotropic homogeneous and layered shells based on the developed third‐order formulation. To assess the performance of the present shell element, it is implemented in the general non‐linear explicit dynamic FE code DYNA3D. Several problems are solved and results are presented and compared to other theoretical and numerical results. Copyright © 2002 John Wiley & Sons, Ltd.