Stability analysis of FGM sandwich plates by using variable-kinematics Ritz models

Abstract

A linearized buckling analysis of functionally graded material (FGM) isotropic and sandwich plates is carried out by virtue of the Hierarchical Trigonometric Ritz Formulation (HTRF). Quasi-3D Ritz models based on equivalent single layer (ESL) and zig zag (ZZ) plate theories are developed within the framework of the Carrera Unified Formulation (CUF). Several in-plane loading conditions accounting for axial, biaxial, and shear loadings are taken into account. Parametric studies are carried out in order to evaluate the effects of significant parameters, such as volume fraction index, length-to-thickness ratio, sandwich plate type, and loading type, on the critical buckling loads.     

Free vibration and buckling analysis of composite laminated plates using layerwise models based on isogeometric approach and Carrera unified formulation

In this paper, layerwise (LW) theory has been utilized along with that of equivalent single layer (ESL) for free vibration and linearized buckling analysis of composite laminated plates. To this end, the isogeometric approach and Carrera unified formulation (CUF) have been combined. In particular, Taylor-like and Legendre-like polynomial expansions have been utilized in the framework of CUF to approximate the solution field in ESL and LW models, respectively. Indeed, CUF provides an ideal tool that facilitates the implementation of higher orders of the solution field expansion. As ESL model cannot inherently provide interlaminar continuity, they are not suitable for analyzing thick laminated plates. However, the LW model not only presents a three-dimensional (3D)-like accurate mathematical model using two-dimensionalplate theories but also considers interlaminar continuity requirements and obviates the need for the use of shear correction factor. In addition, the nonuniform rational B-spline basis functions have been employed to approximate the solution field, due to their interesting attributes in the analysis. These functions are able to describe the exact geometry of the structure and make it technically feasible to provide refinement process during analysis. The presented numerical results confirm the validity of the proposed methodology.     

Extension of MITC to higher‐order beam models and shear locking analysis for compact,thin‐walled,and composite structures

A class of mixed interpolated beam elements is introduced in this paper under the framework of the Carrera Unified Formulation to eliminate the detrimental effects due to shear locking. The Mixed Interpolation of Tensorial Components (MITC) method is adopted to generate locking‐free displacement‐based beam models using general 1D finite elements. An assumed distribution of the transverse shear strains is used for the derivation of the virtual work, and the full Gauss‐Legendre quadrature is used for the numerical computation of all the components of the stiffness matrix. Linear, quadratic, and cubic beam elements are developed using the unified formulation and applied to linear static problems including compact, laminated, and thin‐walled structures. A comprehensive study of how shear locking affects general beam elements when different classical integration schemes are used is presented, evidencing the outstanding capabilities of the MITC method to overcome this numerical issue. Refined beam theories based on the expansion of pure and generalized displacement variables are implemented making use of Lagrange and Legendre polynomials over the cross‐sectional domain, allowing one to capture complex states of stress with a 3D‐like accuracy. The numerical examples are compared to analytic, numerical solutions from the literature, and commercial software solutions, whenever it is possible. The efficiency and robustness of the proposed method demonstrated throughout all the assessments, illustrating that MITC elements are the natural choice to avoid shear locking and showing an unprecedent accuracy in the computation of transverse shear stresses for beam formulations.     

Aerodynamic and mechanical hierarchical aeroelastic analysis of composite wings

Abstract

The coupled bending-torsion flutter is here investigated through Carrera Unified Formulation (CUF). The hierarchical capabilities of CUF offer a procedure to obtain refined one-dimensional models that, by going beyond the assumptions of classical theories, accurately describe the kinematics of structures. Aerodynamic loadings have been determined according to Theodorsen theory, from which the steady formulation can be easily obtained. The displacement variables over the cross section (x-z plane) are approximated by x,z polynomials of any order, N. The finite element method is used to solve the governing equations, which are derived in a weak form through the principle of virtual displacements. The equations are written in terms of “fundamental nuclei,” which do not vary with the theory order, N. Several wing configurations have been studied, giving great attention to thin-walled box beams made of orthotropic material. The effects of sweep angle and lamination scheme on flutter conditions have been investigated, and the results have been compared with solutions obtained from two-dimensional theories, experimental tests, and aeroelastic analyses carried out with the doublet lattice method (DLM). The unsteady theory, combined with advanced beam theories, represents a computationally cheap tool for preliminary aeroelastic studies of complex wing structures.     

An analytical approach for local buckling analysis of initially delaminated composite thin-walled columns with open and closed sections

Local buckling of intact thin-walled columns is generally performed by modeling the wall segments as long plates and by assuming that edges common to two or more plates remain straight. Thus, the buckling load can be determined by considering the wall segments as individual plates rotationally restrained by the adjacent wall segments. This technique is combined with plate theories as a new analytical method to predict the buckling load of an initially delaminated column with any arbitrary sections (open or closed). First, moments at the rotationally restrained edges of delaminated segment (web or flange) are obtained from the curvature and stiffness of the adjacent laminates. Then, the strain energy of this delaminated segment with distributed moment at edges is calculated based on the first-order shear deformation theory. Using the principal of minimum potential energy, the governing equations are obtained and solved by the Rayleigh–Ritz approximation technique. Results of the present approach are compared with three-dimensional finite-element results obtained from eigenvalue buckling analysis in ANSYS software for both box- and channel-section columns with cross-ply and angle-ply stacking sequences. Finally, the effects of delamination size and location are investigated on the buckling loads.     

开口及边界条件对复合材料三分之一柱面壳压缩屈曲性能的影响

研究了完整、开口周边加强及开口加口盖3种型式的复合材料三分之一柱面壳的压缩屈曲性能,考查了3种典型复合材料柱面壳的轴压屈曲强度,分析了开口及口盖对柱面壳压缩稳定性的影响.结果表明:开口大大降低了柱面壳的轴压屈曲强度;口盖可以部分恢复其强度,但很难达到开口之前的水平.进行了开口加口盖经编织物铺层三分之一柱面壳轴向压缩试验,其轴压屈曲强度比用平面织物制造的相同结构的降低很多.为了探究其轴压屈曲强度比同类结构偏低很多的原因,进行了非均匀加载复合材料柱面壳模型有限元分析.结果表明:柱面壳边界不均匀加载会降低其承载能力,根据柱面壳刚度分布制定边界载荷可以提高其承载能力.     

基于密度分布曲线法的复合材料变刚度铺层优化

基于梯度的优化方法对复合材料层合板进行了变刚度铺层优化设计。在优化过程中需确定铺层中各单元的密度以及角度。为了使优化结果具有可制造性,优化结果需满足制造工艺约束并且铺层角度需从预定角度中选取。为了避免在优化问题中引入过多的约束并减少设计变量的数目,提出密度分布曲线法(DDCM)对层合板中各单元的密度进行参数化。根据各单元的密度以及角度设计变量并基于Bi-value Coding Parameterization(BCP)方法中的插值公式确定各单元的弹性矩阵。优化过程中以结构柔顺度作为优化目标,结构体积作为约束,优化算法采用凸规划对偶算法。对碳纤维复合材料的算例结果表明:采用DDCM可得到较理想的优化结果,并且收敛速率较快。     

A variable kinematic doubly-curved MITC9 shell element for the analysis of laminated composites

The present article considers the linear static analysis of composite shell structures with double-curvature geometry by means of a shell finite element with variable through-the-thickness kinematic. The refined models used are grouped in the Unified Formulation by Carrera (CUF) and they permit the distribution of displacements and stresses along the thickness of the multilayered shell to be accurately described. The shell element has nine nodes and the mixed interpolation of tensorial components (MITC) method is used to contrast the membrane and shear locking phenomenon. The governing equations are derived from the principle of virtual displacement (PVD) and the finite element method (FEM) is employed to solve them. Cross-ply spherical shells with simply-supported edges and subjected to bi-sinusoidal pressure are analyzed. Various laminations, thickness ratios, and curvature ratios are considered. The results, obtained with different theories contained in the CUF, are compared with both the elasticity solutions given in the literature and the analytical solutions obtained using the CUF and the Navier’s method. From the analysis, one can conclude that the shell element based on the CUF is very efficient and its use is mandatory with respect to the classical models in the study of composite structures. Finally, shells with different lamination, boundary conditions, and loads are also analyzed using high-order layer-wise theories in order to provide FEM benchmark solutions.     

A direct stiffness analysis of a composite beam with partial interaction

The use of the conventional semi-analytical stiffness method in finite element analysis, in which interpolation polynomials are used to develop the stiffness relationships, leads to problems of curvature locking when beam-type elements are developed for composite members with partial interaction between the materials of which it is comprised. The curvature locking phenomenon that occurs for composite steel–concrete members is quite well reported, and the general approach to minimizing the undesirable ramifications of curvature locking has been to use higher-order polynomials with increasing numbers of internal nodes. This paper presents an alternate formulation based on a direct stiffness approach rather than starting from pre-defined interpolation polynomials, and which does not possess the undesirable locking characteristics. The formulation is based on a more general approach for a bi-material composite flexural member, whose constituent materials are joined by elastic shear connection so as to provide partial interaction. The stiffness relationships are derived, and these are applied to a simply supported and a continuous steel–concrete composite beam to demonstrate the efficacy of the method, and in particular its ability to model accurately both very flexible and very stiff shear connection that causes difficulties when implemented in competitive semi-analytical algorithms. Copyright © 2004 John Wiley & Sons, Ltd.     

中厚椭球壳自由振动动力刚度法分析

介绍精确动力刚度法分析中厚椭球壳自由振动具体实施方法,据环向波数不同将中厚椭球壳自由振动分解为一系列确定环向波数的一维振动;利用控制方程Hamilton形式建立动力刚度关系,用常微分方程求解器COLSYS求解控制方程获得单元动力刚度,用Wittrick-Williams算法求得该环向波数下椭球壳自振频率。数值算例给出中厚圆球壳及椭球壳不同边界条件的自振频率,验证动力刚度法高效、可靠、精确。     

汽车排气系统吊耳动刚度优化方法的研究

基于汽车排气系统吊耳传递的动态载荷最小、吊耳耐疲劳性最好,建立了考虑动力总成在内的排气系统振动分析模型。进行了排气系统的自由模态和约束模态的测试,并和计算值进行了对比分析,证明了所建立的排气系统振动模型的正确性。以吊耳的垂向动态载荷最小和其静变形量在一定范围内为优化目标,建立了排气系统吊耳动刚度优化模型。优化后,在怠速工况和2档全负荷加速工况下对车身底板驾驶员位置进行了振动响应测试,测试结果表明,利用优化后的吊耳刚度,能够有效降低车身底板的振动加速度,表明了阐述的排气系统建模和吊耳动刚度优化方法的有效性。文中建模与优化方法,对排气系统的吊耳动刚度计算与优化具有指导意义。     

Free vibration analysis of sandwich beams using improved dynamic stiffness method

In this paper, free vibration of three-layered symmetric sandwich beam is investigated using dynamic stiffness and finite element methods. To determine the governing equations of motion by the present theory, the core density has been taken into consideration. The governing partial differential equations of motion for one element contained three layers are derived using Hamilton’s principle. This formulation leads to two partial differential equations which are coupled in axial and bending deformations. For the harmonic motion, these equations are combined to form one ordinary differential equation. Closed form analytical solution for this equation is determined. By applying the boundary conditions, the element dynamic stiffness matrix is developed. They are assembled and the boundary conditions of the beam are applied, so that the dynamic stiffness matrix of the beam is derived. Natural frequencies and mode shapes are computed by the use of numerical techniques and the known Wittrick–Williams algorithm. After validation of the present model, the effect of various parameters such as density, thickness and shear modulus of the core for various boundary conditions on the first natural frequency is studied.     

Thermal buckling and postbuckling analysis of a laminated composite beam with embedded SMA actuators

In this paper, the thermal buckling and postbuckling behaviours of a composite beam with embedded shape memory alloy (SMA) wires are investigated analytically. For the purpose of enhancing the critical buckling temperature and reducing the lateral deflection for the thermal buckling, the characteristics of thermal buckling are investigated through the use of the shape recovery force associated with SMA wire actuators. The results of both thermal buckling and postbuckling behaviours present quantitatively how the shape recovery force affects the thermal buckling behaviour. The analytical results show that the shape recovery force reduces the thermal expansion of the composite laminated beam, which results in both an increment of the critical buckling temperature and also a reduction of the lateral deflection of postbuckling behaviours. A new formula is also proposed to describe the critical buckling temperature of the laminated composite beam with embedded SMA wire actuators.     

基于动刚度法的热环境下纤维增强聚合物基复合薄板抗振性能的退化

从动力学角度,提出了动刚度法来研究热环境下纤维增强聚合物基复合薄板的抗振性能退化特性。首先,考虑脉冲激励载荷的影响,利用能量法、板壳理论和振型叠加法等,实现了热环境下动刚度的理论求解。同时,总结和归纳了热环境下复合薄板抗振性能退化的分析流程,并以TC500碳纤维/树脂基复合薄板为例,通过对比理论与测试获得的动刚度、固有频率、阻尼和振型结果,证明了该动态指标的有效性和方法的可行性,可以利用该指标来量化评价复合薄板在热环境下的动态性能退化问题。      

Nonlinear static and dynamic buckling analysis of functionally graded shallow spherical shells including temperature effects

This paper presents an analytical approach to investigate the nonlinear static and dynamic unsymmetrical responses of functionally graded shallow spherical shells under external pressure incorporating the effects of temperature. Governing equations for thin FGM spherical shells are derived by using the classical shell theory taking into account von Karman–Donnell geometrical nonlinearity. Approximate solutions are assumed and Galerkin procedure is applied to determine explicit expressions of static critical buckling loads of the shells. For the dynamical response, motion equations are numerically solved by using Runge–Kutta method and the criterion suggested by Budiansky–Roth. A detailed analysis is carried out to show the effects of material and geometrical parameters, boundary conditions and temperature on the stability and dynamical characteristics of FGM shallow spherical shells.     

Post buckling analysis of shear deformable shallow shells by the boundary element method

This paper presents four boundary element formulations for post buckling analysis of shear deformable shallow shells. The main differences between the formulations rely on the way non‐linear terms are treated and on the number of degrees of freedom in the domain. Boundary integral equations are obtained by coupling boundary element formulation of shear deformable plate and two‐dimensional plane stress elasticity. Four different sets of non‐linear integral equations are presented. Some domain integrals are treated directly with domain discretization whereas others are dealt indirectly with the dual reciprocity method. Each set of non‐linear boundary integral equations are solved using an incremental approach, where loads and prescribed boundary conditions are applied in small but finite increments. The resulting systems of equations are solved using a purely incremental technique and the Newton–Raphson technique with the Arc length method. Finally, the effect of imperfections (obtained from a linear buckling analysis) on the post‐buckling behaviour of axially compressed shallow shells is investigated. Results of several benchmark examples are compared with the published work and good agreement is obtained. Copyright © 2010 John Wiley & Sons, Ltd.     

High-order free vibration analysis of sandwich beams with a flexible core using dynamic stiffness method

In this paper, high-order free vibration of three-layered symmetric sandwich beam is investigated using dynamic stiffness method. The governing partial differential equations of motion for one element are derived using Hamilton’s principle. This formulation leads to seven partial differential equations which are coupled in axial and bending deformations. For the harmonic motion, these equations are divided into two ordinary differential equations by considering the symmetrical sandwich beam. Closed form analytical solutions of these equations are determined. By applying the boundary conditions, the element dynamic stiffness matrix is developed. The element dynamic stiffness matrices are assembled and the boundary conditions of the beam are applied, so that the dynamic stiffness matrix of the beam is derived. Natural frequencies and mode shapes are computed by use of numerical techniques and the known Wittrick–Williams algorithm. Finally, some numerical examples are discussed using dynamic stiffness method.     

A refined dynamic stiffness element for free vibration analysis of cross-ply laminated composite cylindrical and spherical shallow shells

An exact free vibration analysis of doubly-curved laminated composite shallow shells has been carried out by combining the dynamic stiffness method (DSM) and a higher order shear deformation theory (HSDT). In essence, the HSDT has been exploited to develop first the dynamic stiffness (DS) element matrix and then the global DS matrix of composite cylindrical and spherical shallow shell structures by assembling the individual DS elements. As an essential prerequisite, Hamilton’s principle is used to derive the governing differential equations and the related natural boundary conditions. The equations are solved symbolically in an exact sense and the DS matrix is formulated by imposing the natural boundary conditions in algebraic form. The Wittrick–Williams algorithm is used as a solution technique to compute the eigenvalues of the overall DS matrix. The effect of several parameters such as boundary conditions, orthotropic ratio, length-to-thickness ratio, radius-to-length ratio and stacking sequence on the natural frequencies and mode shapes is investigated in details. Results are compared with those available in the literature. Finally some concluding remarks are drawn.     

Recursive second order convergence method for natural frequencies and modes when using dynamic stiffness matrices

When exact dynamic stiffness matrices are used to compute natural frequencies and vibration modes for skeletal and certain other structures, a challenging transcendental eigenvalue problem results. The present paper presents a newly developed, mathematically elegant and computationally efficient method for accurate and reliable computation of both natural frequencies and vibration modes. The method can also be applied to buckling problems. The transcendental eigenvalue problem is first reduced to a generalized linear eigenvalue problem by using Newton's method in the vicinity of an exact natural frequency identified by the Wittrick–Williams algorithm. Then the generalized linear eigenvalue problem is effectively solved by using a standard inverse iteration or subspace iteration method. The recursive use of the Newton method employing the Wittrick–Williams algorithm to guide and guard each Newton correction gives secure second order convergence on both natural frequencies and mode vectors. The second order mode accuracy is a major advantage over earlier transcendental eigenvalue solution methods, which typically give modes of much lower accuracy than that of the natural frequencies. The excellent performance of the method is demonstrated by numerical examples, including some demanding problems, e.g. with coincident natural frequencies, with rigid body motions and large‐scale structures. Copyright © 2003 John Wiley & Sons, Ltd.     

Transient dynamic analysis of tapered FRP composite transmission poles using finite element method

Power transmission poles are subjected to dynamic cantilever bending due to wind gusts and cable unilateral failure, or may also be subjected to vehicle impacts. In this paper, transient dynamic analysis of tapered fiber-reinforced polymer (FRP) composite transmission poles with circular thin-walled cross-section subjected to dynamic cable tension and vehicle impacts is investigated by combination of tapered beam finite element and precise time integration method. It is assumed that the material behavior is linearly elastic and the laminate of the cross-section of the wall is symmetric or antisymmetric angle-ply. The effect of fiber type and orientation, the pole geometry, and the concentrated mass at the pole tip are evaluated by performing the dynamic analysis of FRP poles under step, triangular and sine pulses. There is a good agreement between the results of the present method and those obtained from the poles modeled by ANSYS commercial finite element software and existed literatures. Also, there is a significant shorter run-time in the present method. It is concluded that beyond 10 layers for the laminate with constant thickness of the wall, the pole tip deflection does not change.