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

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

A hybrid multiscale finite element/peridynamics method

A hybrid method is presented that uses a representative volume element-based multiscale finite element technique combined with a peridynamics method for modeling fracture surfaces. The hybrid method dynamically switches from finite element computations to peridynamics based on a damage criterion defined on the peridynamics grid, which is coincident with the nodes of the finite element mesh. Nodal forces are either computed by the finite element method or peridynamics, as appropriate. The multiscale finite element method used here is a representative volume element-based approach so that inhomogeneous local scale material properties can be derived using homogenization. In addition, automatic cohesive zone insertion is used at the local scale to model fracture initiation. Results demonstrate that local scale flaw distributions can alter fracture patterns and initiation times, and the use of cohesive zone insertion can improve accuracy of crack paths.     

采用等效刚度有限元模型的复合材料机翼颤振分析

采用等效刚度方法,研究了一种适用于机翼初步设计阶段的动力学和颤振分析的结构有限元模型。该方法首先计算不同布局形式的加筋壁板的刚度矩阵,然后将其赋予与加筋壁板平面形状相同的光板（等效板）上,使加筋壁板和等效板具相同的力学性能。该方法的优点是避免了加强筋的有限元建模,从而使有限元模型的复杂程度大大降低,但同时等效刚度结构有限元模型仍能反映机翼加筋壁板的结构特性。以某客机概念方案的机翼为例,建立了反映实际结构详细有限元及其等效刚度有限元模型。计算结果和对比分析表明,两种模型的固有频率、振动模态和颤振分析结果吻合得很好,从而验证了等效刚度方法在机翼结构动力学和颤振分析方面的准确性。由于该方法具有简单快速和准确的优点,可用于机翼初步设计阶段对颤振特性的评估。     

Simulations of piezoelectric Lamb wave delay lines using a finite element method

The method of piezoelectric finite elements was applied to the simulation of piezoelectric Lamb-wave delay lines with and without acoustical absorbers. In this finite-element analysis, free as well as electrically driven vibrations were computed. The shapes of the symmetric Lamb modes were determined by the solution of eigenproblems, and the transient mechanical build-up process was studied for a switched electrical sine wave excitation. The transfer function, the group delay time, and the impedance matrix of devices of different designs are calculated. The good agreement between simulation and experimental results indicates the suitability of the finite-element method for optimizing acoustic delay-line devices.     

Modeling of composite piezoelectric structures with the finite volume method

Piezoelectric devices, such as piezoelectric traveling- wave rotary ultrasonic motors, have composite piezoelectric structures. A composite piezoelectric structure consists of a combination of two or more bonded materials, at least one of which is a piezoelectric transducer. Piezoelectric structures have mainly been numerically modeled using the finite element method. An alternative approach based on the finite volume method offers the following advantages: 1) the ordinary differential equations resulting from the discretization process can be interpreted directly as corresponding circuits; and 2) phenomena occurring at boundaries can be treated exactly. This paper presents a method for implementing the boundary conditions between the bonded materials in composite piezoelectric structures modeled with the finite volume method. The paper concludes with a modeling example of a unimorph structure.     

Analysis of cracked anisotropic plates using the hybrid finite element method

This paper presents a convenient and efficient method to obtain accurate stress intensity factors for cracked anisotropic plates. In this method, a complex variable formulation in conjunction with a hybrid displacement finite element scheme is used to carry out the stiffness and stress calculations of finite cracked plates subjected to general boundary and loading conditions. Unlike other numerical methods used for local analysis such as the boundary element method, the present method results in a symmetric stiffness matrix, which can be directly incorporated into the stiffness matrix representing other structural parts modeled by conventional finite elements. Therefore, the present method is ideally suited for modeling cracked plates in a large complex structure.     

Mechanical properties of hybrid composites using finite element method based micromechanics

A micromechanical analysis of the representative volume element of a unidirectional hybrid composite is performed using finite element method. The fibers are assumed to be circular and packed in a hexagonal array. The effects of volume fractions of the two different fibers used and also their relative locations within the unit cell are studied. Analytical results are obtained for all the elastic constants. Modified Halpin–Tsai equations are proposed for predicting the transverse and shear moduli of hybrid composites. Variability in mechanical properties due to different locations of the two fibers for the same volume fractions was studied. It is found that the variability in elastic constants and longitudinal strength properties was negligible. However, there was significant variability in the transverse strength properties. The results for hybrid composites are compared with single fiber composites.     

Solving the nonlinear Poisson-type problems with F-Trefftz hybrid finite element model

A hybrid finite element model based on F-Trefftz kernels (fundamental solutions) is formulated for analyzing Dirichlet problems associated with two-dimensional nonlinear Poisson-type equations including nonlinear Poisson-Boltzmann equation and diffusion-reaction equation. The nonlinear force term in the Poisson-type equation is frozen by introducing the imaginary terms at each Picard iteration step, and then the induced Poisson problem is solved by the present hybrid finite element model involving element boundary integrals only, coupling with the particular solution method with radial basis function interpolation. The numerical accuracy of the present method is investigated by numerical experiments for problems with complex geometry and various nonlinear force functions.     

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.     

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

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

Axisymmetric dynamics of composite spherical shells with active piezoelectric/composite stiffeners

Summary The paper presents a theoretical formulation for spherical shells reinforced by meridional and circumferential stiffeners. Active damping of the shell is introduced through control action of piezoelectric coupled pairs bonded to the meridional stiffeners. The induced loads can include radial pressure and a thermal field that are independent of the circumferential coordinate. Neglecting local deformations between adjacent meridional stiffeners, the response of the shell will be axisymmetric. The analysis employs the Donnell-Mushtari-Vlasov version of Love's theory of shells together with a smeared stiffeners technique. The paper also considers a particular case of shell mounted piezoelectic coupled pairs without conventional stiffeners. A closed form solution is derived for spherical panels without conventional stiffeners within the range of the meridional coordinate between 75° and 90° using a version of the Geckeler approximation.     

A study of finite size effects on cracked 2-D piezoelectric media using extended finite element method

The extended finite element method is applied on the two-dimensional (2-D) finite piezoelectric media weakened by a crack. The fourfold standard enrichment functions are taken in conjugation with the interaction integral to evaluate the intensity factors. Four sequence of analysis, namely crack–mesh alignment, aspect ratio, mesh with local refinement and domain independency is done on the center and edge crack problems. These four analyses when combined together give an optimum result to study the finite specimen. It is observed that for smaller values of strip width to crack length ratio the finiteness of the specimen size affects the intensity factors.     

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.     

二维压电材料动强度因子的扩展有限元计算

采用扩展有限元求解二维弹性压电材料动断裂问题。扩展有限元的网格独立于裂纹,因此网格生成可大大地简化,且裂纹扩展时不需重构网格。采用相互作用积分技术计算动强度因子。比较了标准的力裂尖加强函数和力-电裂尖加强函数对动强度因子的影响,结果表明标准的力裂尖加强函数能有效地分析压电材料动断裂问题。分析了极化方向对动强度因子的影响。数值分析表明采用扩展有限元获得的动强度因子与其他数值方法解吻合得很好。     

Multiscale design of three-dimensional nonlinear composites using an interface-enriched generalized finite element method

A computational framework is developed to model and optimize the nonlinear multiscale response of three-dimensional particulate composites using an interface-enriched generalized finite element method. The material nonlinearities are associated with interfacial debonding of inclusions from a surrounding matrix which is modeled using C⁻¹ continuous enrichment functions and a cohesive failure model. Analytic material and shape sensitivities of the homogenized constitutive response are derived and used to drive a nonlinear inverse homogenization problem using gradient-based optimization methods. Spherical and ellipsoidal particulate microstructures are designed to match a component of the homogenized stress-strain response to a desired constructed macroscopic stress-strain behavior.     

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.     

The finite element method in problems of nonlinear optics

The Possibility of the application of the finite element method to some problems of nonlinear optics is investigated in this paper. The self-action of a light beam in a nonlinear medium is considered. The general approach to the cretion of conservative computation schemes is presented, based on varitional principles. Definite schemes, which are applicable for the problem of thermal self-action, are described in detail both in the case of cylindrical and or rectangular co-ordinates. The accuracy and convergence of the models are analysed. The results of computation of the self-action problems in motionless and moving media are presented.     

An efficient finite element with layerwise mechanics for smart piezoelectric composite and sandwich shallow shells

In this work, we present a new efficient four-node finite element for shallow multilayered piezoelectric shells, considering layerwise mechanics and electromechanical coupling. The laminate mechanics is based on the zigzag theory that has only seven kinematic degrees of freedom per node. The normal deformation of the piezoelectric layers under the electric field is accounted for without introducing any additional deflection variables. A consistent quadratic variation of the electric potential across the piezoelectric layers with the provision of satisfying the equipotential condition of electroded surfaces is adopted. The performance of the new element is demonstrated for the static response under mechanical and electric potential loads, and for free vibration response of smart shells under different boundary conditions. The predictions are found to be very close to the three dimensional piezoelasticity solutions for hybrid shells made of not only single-material composite substrates, but also sandwich substrates with a soft core for which the equivalent single layer (ESL) theories perform very badly.