A 3D cylindrical finite element model for thick curved beam stress analysis

A finite element model is proposed to perform stress analysis for thick curved beams and panels subjected to various types of loadings. The model has 18 nodes in a three‐dimensional cylindrical co‐ordinates system. Three stress components on radial surface (σ_rr, τ_rθ, and τ_rz) and three displacement components (u_r, u_θ, and u_z) are used as nodal degrees of freedom. Therefore, the continuity condition for both stresses and displacements is achieved in the radial direction. Formulation of nodal shape functions and equilibrium equations are based on three‐dimensional elasticity theory and a minimum potential energy method. The accuracy of the method is verified with the standard test problems and exact solutions from the theory of elasticity. The model shows no locking phenomena. Convergence is investigated and the application to layered composite panel is illustrated. Copyright © 2003 John Wiley & Sons, Ltd.     

A C1 finite element for flexural and torsional analysis of rectangular piezoelectric laminated/sandwich composite beams

This work deals with the development of a new C¹ finite element for analysing the bending and torsional behaviour of rectangular piezoelectric laminated/sandwich composite beams. The formulation includes transverse shear, warping due to torsion, and elastic–electric coupling effects. It also accounts for the inter-layer continuity condition at the interfaces between layers, and the boundary conditions at the upper and lower surfaces of the beam. The shear strain is represented by a cosine function of a higher order in nature and thus avoiding shear correction factors. The warping function obtained from a three-dimensional elasticity solution is incorporated in the present model. An exact integration is employed in evaluating various energy terms due to the application of field consistency approach while interpolating the transverse shear and torsional strains. The variation of the electric potential through the thickness is taken care of in the formulation based on the observation of three-dimensional solution. The performance of the laminated piezoelectric element is tested comparing with analytical results as well as with the reference solutions evaluated using three-dimensional finite element procedure. A detailed study is conducted to highlight the influence of length-to-thickness ratio on the displacements, stresses and electric potential field of piezoelectric laminated beam structures subjected to flexural and torsional loadings. Copyright © 2004 John Wiley & Sons, Ltd.     

基于有限元法的周期压电智能梁滤波特性分析

周期结构具有通频和禁频特性,使其在动态载荷的滤波器、具有主动控制功能的结构研究中得到了重要应用。基于Timoshenko梁理论,考虑基梁和压电片的转动惯量和剪切效应,采用有限元法和传递矩阵法推导了波在周期性地粘贴压电片的Timoshenko梁中的传播模型,分析了几何尺寸和材料特性对其频带性质的影响,并与Bernoulli-Euler梁理论得到的结果进行了对比。研究表明,当基梁与压电层厚度比达到40时,禁带带宽减小了54%,因此对于周期结构中的深梁,应舍弃Bernoulli-Euler梁理论而采用Timoshenko梁理论建立的模型;对于不同尺寸和材料特性的压电周期结构,频带性质会有很大不同,可以通过调整结构的参数来改变其频带性质,从而改变波动在结构中的传播特性。     

The closed form general solutions of 2-D curved laminated beams of variable curvatures

The closed form general solutions of in-plane static laminated curved beams of variable curvatures are derived. All the quantities of axial force, shear force, radial and tangential displacements are decoupled and expressed as harmonic functions in terms of angle of tangent slope. The first and second moments of laminated curved beams with respect to x- and y-axes are defined as fundamental geometric properties. The curved beams can be any curve. Some closed form solutions of non-circular laminated beams, rings are presented. Applications of elliptic, parabola, catenary, cycloid, and exponential spiral laminated curved beams are shown. The laminated circular ring under point load and distributed load are studied as well. Two kinds of laminate (0/θ) and (θ/−θ) are investigated in this paper. Effects of aspect ratio, thickness ratios, orthotropy ratios and stacking sequence on the laminated ring are presented. Results of isotropic curved beams are compared with existed solutions. The results are consistent.     

A semi-analytical finite element for laminated composite plates

This paper presents a semi-analytical finite element solution for laminated composite plates. The method is based on a mixed variational principle that involves both displacements and stresses. Finite element meshes are only used in the plane of plate, while the through thickness distributions of displacements and stresses are obtained using the method of state equations. Numerical results show that the rate of convergence of the new method is fast and the solutions can be very close to corresponding exact three-dimensional ones. The use of a recursive formulation of the state equations leads to an algebra equation system, from which solution are sought, whose dimension is independent of the numbers of layers of the plate considered.     

An accurate two-node finite element for shear deformable curved beams

An accurate two-node (three degrees of freedom per node) finite element is developed for curved shear deformable beams. The element formulation is based on shape functions that satisfy the homogeneous form of the partial differential equations of motion which renders it free of shear and membrane locking. The element is demonstrated to converge to the results obtained from a shear deformable straight beam when the beam becomes shallower. Numerical examples were performed to demonstrate the accuracy and efficiency with respect to previously published formulations. © 1998 John Wiley & Sons, Ltd.     

一种基于Layerwise理论的压电复合材料层合圆柱壳机电耦合有限单元

基于Reddy的Layerwise理论,对含压电铺层的复合材料层合壳的静力响应特性进行了理论研究。基于Layerwise理论,推导了含压电层的复合材料层合壳的应变分量与电场强度表达式。利用Hamilton原理和变分法,推导了压电智能层合壳的欧拉-拉格朗日方程,并采用有限元解法,建立了相应的有限元控制方程及其机电耦合刚度矩阵。通过算例结果与文献中的精确解和试验值进行了对比,表明相较于传统的经典层合板壳理论,本文理论方法的有效性和优势性;并分析了径厚比等参量对两端简支压电智能层合壳静力响应值的影响规律。      

A new higher-order hybrid-mixed curved beam element

The purpose of this work is to show the successful use of nodeless degrees of freedom in developing a highly accurate, locking free hybrid-mixed C⁰ curved beam element. In the performance evaluation process of the present field-consistent higher-order element, the effect of field consistency and the role of higher-order interpolation on both displacement-type and hybrid-mixed-type elements are carefully examined. Several benchmark tests confirm the superior behaviour of the present element. © 1998 John Wiley & Sons, Ltd.     

具任意脱层压电弹性层合梁的谐响应分析

基于Timoshenko梁理论,考虑剪切变形和力电耦合效应,对压电弹性层合脱层梁进行了分区处理,可方便地描述脱层长度和脱层位置,建立了弯曲波在梁中传播的基本方程式;通过长波情况下的电势分布简化,考虑无外加电压时,对该类悬臂梁进行了求解分析。给出了梁的频散关系,得到了各位移分量和电势的解析解,分析了脱层参数对位移和电势的影响规律。     

A semi‐analytical curved element for linear elasticity based on the scaled boundary finite element method

This work introduces a semi‐analytical formulation for the simulation and modeling of curved structures based on the scaled boundary finite element method (SBFEM). This approach adapts the fundamental idea of the SBFEM concept to scale a boundary to describe a geometry. Until now, scaling in SBFEM has exclusively been performed along a straight coordinate that enlarges, shrinks, or shifts a given boundary. In this novel approach, scaling is based on a polar or cylindrical coordinate system such that a boundary is shifted along a curved scaling direction. The derived formulations are used to compute the static and dynamic stiffness matrices of homogeneous curved structures. The resulting elements can be coupled to general SBFEM or FEM domains. For elastodynamic problems, computations are performed in the frequency domain. Results of this work are validated using the global matrix method and standard finite element analysis. Copyright © 2016 John Wiley & Sons, Ltd.     

Time and frequency domain analysis of piezoelectric energy harvesters by monolithic finite element modeling

The successful design of piezoelectric energy harvesting devices relies upon the identification of optimal geometrical and material configurations to maximize the power output for a specific band of excitation frequencies. Extendable predictive models and associated approximate solution methods are essential for analysis of a wide variety of future advanced energy harvesting devices involving more complex geometries and material distributions. Based on a holistic continuum mechanics modeling approach to the multi‐physics energy harvesting problem, this article proposes a monolithic numerical solution scheme using a mixed‐hybrid 3‐dimensional finite element formulation of the coupled governing equations for analysis in time and frequency domain. The weak form of the electromechanical/circuit system uses velocities and potential rate within the piezoelectric structure, free boundary charge on the electrodes, and potential at the level of the generic electric circuit as global degrees of freedom. The approximation of stress and dielectric displacement follows the work by Pian, Sze, and Pan. Results obtained with the proposed model are compared with analytical results for the reduced‐order model of a cantilevered bimorph harvester with tip mass reported in the literature. The flexibility of the method is demonstrated by studying the influence of partial electrode coverage on the generated power output.     

含裂纹压电材料的Cell-Based光滑扩展有限元法

为精确而有效地求解机电耦合作用下含裂纹压电材料的断裂参数,首先,通过将复势函数法、扩展有限元法和光滑梯度技术引入到含裂纹压电材料的断裂机理问题中,提出了含裂纹压电材料的Cell-Based光滑扩展有限元法;然后,对含中心裂纹的压电材料强度因子进行了模拟,并将模拟结果与扩展有限元法和有限元法的计算结果进行了对比。数值算例结果表明:Cell-Based光滑扩展有限元法兼具扩展有限元法和光滑有限元法的特点,不仅单元网格与裂纹面相互独立,且裂尖处单元不需精密划分,与此同时,Cell-Based光滑扩展有限元法还具有形函数简单且不需求导、对网格质量要求低且求解精度高等优点。所得结论表明Cell-Based光滑扩展有限元法是压电材料断裂分析的有效数值方法。      

A three‐noded shear‐flexible curved beam element based on coupled displacement field interpolations

An efficient shear‐flexible three‐noded curved beam element is proposed herein. The shear flexibility is based on Timoshenko beam theory and the element has three degrees of freedom, viz., tangential displacement (u), radial displacement (w) and the section‐rotation (θ). A quartic polynomial interpolation for flexural rotation ψ is assumed a priori. Making use of the physical composition of θ in terms of ψ and u, a novel way of deriving the polynomial interpolations for u and w is presented, by solving force‐moment and moment‐shear equilibrium equations simultaneously. The field interpolation for θ is then constructed from that of ψ and u. The procedure leads to high‐order polynomial field interpolations which share some of the generalized degrees of freedom, by means of coefficients involving material and geometric properties of the element. When applied to a straight Euler–Bernoulli beam, all the coupled coefficients vanish and the formulation reduces to classical quintic‐in‐w and quadratic‐in‐u element, with u, w, and ?w/?x as degrees of freedom. The element is totally devoid of membrane and shear locking phenomena. The formulation presents an efficient utilization of the nine generalized degrees of freedom available for the polynomial interpolation of field variables for a three‐noded curved beam element. Numerical examples on static and free vibration analyses demonstrate the efficacy and locking‐free property of the element. Copyright © 2001 John Wiley & Sons, Ltd.     

压电层合板的B样条小波有限元半解析法

利用小波有限元法的优越性可方便地求解压电材料与复合材料混合层合板的某些静力学问题。根据层合结构的特点,将区间B样条尺度函数作为插值函数离散结构的平面域,应用压电材料修正后的H-R(Hellinger-Reissner)变分原理推导了压电材料的Hamilton正则方程的区间B样条小波(BSWI)元列式。该BSWI元的主要特点之一是厚度方向是解析解形式的。针对具体问题的求解,为了保证各层之间力学量和电学量的连续性,进一步应用了状态转移矩阵技术。数值算例表明所提出的区间B样条小波单元是成功的。采用推导压电材料BSWI元的方法可建立磁电弹性材料类似的BSWI元。     

A nonlocal higher-order curved beam finite model including thickness stretching effect for bending analysis of curved nanobeams

In the present work, a finite element approach is developed for the static analysis of curved nanobeams using nonlocal elasticity beam theory based on Eringen formulation coupled with a higher-order shear deformation accounting for through-thickness stretching. The formulation is general in the sense that it can be used to compare the influence of different structural theories, through static and dynamic analyses of curved nanobeams. The governing equations derived here are solved introducing a 3-nodes beam element. The formulation is validated considering problems for which solutions are available. A comparative study is done here by different theories obtained through the formulation. The effects of various structural parameters such as thickness ratio, beam length, rise of the curved beam, loadings, boundary conditions, and nonlocal scale parameter are brought out on the static bending behaviors of curved nanobeams.     

Free vibration analysis of arches using curved beam elements

The natural frequencies and mode shapes for the radial (in‐plane) bending vibrations of the uniform circular arches were investigated by means of the finite arch (curved beam) elements. Instead of the complicated explicit shape functions of the arch element given by the existing literature, the simple implicit shape functions associated with the tangential, radial (or normal) and rotational displacements of the arch element were derived and presented in matrix form. Based on the relationship between the nodal forces and the nodal displacements of a two‐node six‐degree‐of‐freedom arch element, the elemental stiffness matrix was derived, and based on the equation of kinetic energy and the implicit shape functions of an arch element the elemental consistent mass matrix with rotary inertia effect considered was obtained. Assembly of the foregoing elemental property matrices yields the overall stiffness and mass matrices of the complete curved beam. The standard techniques were used to determine the natural frequencies and mode shapes for the curved beam with various boundary conditions and subtended angles. In addition to the typical circular arches with constant curvatures, a hybrid beam constructed by using an arch segment connected with a straight beam segment at each of its two ends was also studied. For simplicity, a lumped mass model for the arch element was also presented. All numerical results were compared with the existing literature or those obtained from the finite element method based on the conventional straight beam element and good agreements were achieved. Copyright © 2003 John Wiley & Sons, Ltd.     

Non‐linear spatial Timoshenko beam element with curvature interpolation

The paper presents a spatial Timoshenko beam element with a total Lagrangian formulation. The element is based on curvature interpolation that is independent of the rigid‐body motion of the beam element and simplifies the formulation. The section response is derived from plane section kinematics. A two‐node beam element with constant curvature is relatively simple to formulate and exhibits excellent numerical convergence. The formulation is extended to N‐node elements with polynomial curvature interpolation. Models with moderate discretization yield results of sufficient accuracy with a small number of iterations at each load step. Generalized second‐order stress resultants are identified and the section response takes into account non‐linear material behaviour. Green–Lagrange strains are expressed in terms of section curvature and shear distortion, whose first and second variations are functions of node displacements and rotations. A symmetric tangent stiffness matrix is derived by consistent linearization and an iterative acceleration method is used to improve numerical convergence for hyperelastic materials. The comparison of analytical results with numerical simulations in the literature demonstrates the consistency, accuracy and superior numerical performance of the proposed element. Copyright © 2001 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.     

Coupled mechanics and finite element for non‐linear laminated piezoelectric shallow shells undergoing large displacements and rotations

A theoretical framework is presented for analysing the coupled non‐linear response of shallow doubly curved adaptive laminated piezoelectric shells undergoing large displacements and rotations. The formulated mechanics incorporate coupling between in‐plane and flexural stiffness terms due to geometric curvature, coupling between mechanical and electric fields, and encompass geometric non‐linearity effects due to large displacements and rotations. The governing equations are formulated explicitly in orthogonal curvilinear co‐ordinates and are combined with the kinematic assumptions of a mixed‐field shear‐layerwise shell laminate theory. Based on the above formulation, a finite element methodology together with an incremental‐iterative technique, based on Newton–Raphson method is formulated. An eight‐node coupled non‐linear shell element is also developed. Various evaluation cases on laminated curved beams and cylindrical panels illustrate the capability of the shell finite element to predict the complex non‐linear behaviour of active shell structures including buckling, which is not captured by linear shell models. The numerical results also show the inherent capability of piezoelectric shell structures to actively induce large displacements through piezoelectric actuators, by jumping between multiple equilibrium states. Copyright © 2005 John Wiley & Sons, Ltd.