Postbuckling of FGM plates with piezoelectric actuators under thermo-electro-mechanical loadings

A postbuckling analysis is presented for a simply supported, shear deformable functionally graded plate with piezoelectric actuators subjected to the combined action of mechanical, electrical and thermal loads. The temperature field considered is assumed to be of uniform distribution over the plate surface and through the plate thickness and the electric field considered only has non-zero-valued component E_Z. The material properties of functionally graded materials (FGMs) are assumed to be graded in the thickness direction according to a simple power law distribution in terms of the volume fractions of the constituents, and the material properties of both FGM and piezoelectric layers are assumed to be temperature-dependent. The governing equations are based on a higher order shear deformation plate theory that includes thermo-piezoelectric effects. The initial geometric imperfection of the plate is taken into account. Two cases of the in-plane boundary conditions are considered. A two step perturbation technique is employed to determine buckling loads and postbuckling equilibrium paths. The numerical illustrations concern the postbuckling behavior of perfect and imperfect, geometrically mid-plane symmetric FGM plates with fully covered or embedded piezoelectric actuators under different sets of thermal and electric loading conditions. The effects played by temperature rise, volume fraction distribution, applied voltage, the character of in-plane boundary conditions, as well as initial geometric imperfections are studied.     

Parameters controlling fracture resistance in functionally graded materials under mode I loading

In this investigation the fracture behavior of functionally graded materials (FGMs) was studied by means of experiments carried out on model polymer-based FGMs. Model graded materials were manufactured by selective ultraviolet irradiation of ECO [poly(ethylene carbon monoxide)], a photo-sensitive ductile copolymer that becomes more brittle and stiffer under exposure to ultraviolet light. The mechanical response of the graded material was characterized using uniaxial tensile tests. Single edge notched tension graded ECO specimens possessing different spatial variations of Young’s modulus, failure stress and failure strain were tested under remote opening loading. A full-field digital image correlation technique was used to measure in real-time the displacement field around the crack tip while it propagated through the graded material. The measured displacement field was then used to extract fracture parameters such as stress intensity factor and T-stress, and thus construct resistance curves for crack growth in the FGMs. For this loading configuration it was found that the nonsingular T-stress term in the asymptotic expansion for stresses needs to be accounted for in order to accurately measure the fracture resistance in FGMs. In addition, the influence of local failure properties (i.e., failure stress and failure strain) on crack growth resistance was investigated in detail. It was found that depending on the combined effects of the spatial variation of these two failure parameters, regardless of the spatial variation of the Young’s modulus, the FGM fracture resistance can either increase, decrease or remain constant with continued crack growth.     

An improved numerical method for computation of stress intensity factors along 3D curved non-planar cracks in FGMs

A simplified strategy based on the interaction energy integral is implemented in the finite element framework to evaluate mixed mode Stress Intensity Factors (SIFs) in 3D non-planar cracks. The proposed approach does not require any a priori information about crack front and crack surface curvatures, therefore different arbitrary non-planar cracks can be easily investigated. In particular, both conical and lens-shaped cracks in homogeneous materials are considered as case studies in order to demonstrate the accuracy of the present approach. Finally, the computational strategy is extended to Functionally Graded Materials (FGMs) and the effect of graded material properties (Young’s modulus and Poisson’s ratio) on the SIFs is studied in detail.     

SH波在压电材料条中垂直界面裂纹处的散射

研究了SH波在压电材料条中裂纹处的散射.压电材料条两侧涂有相同梯度参数的两个半无限大功能梯度材料,裂纹垂直于界面.通过Fourier变换,利用边界条件把问题转化为柯西核奇异积分方程,然后利用Chebyshev多项式对奇异积分方程进行数值求解.通过数值计算,分析讨论了压电条的几何参数和SH波频率对标准动应力强度因子的影响.     

Time dependent behavior of ferroelectric materials undergoing changes in their material properties with electric field and temperature

This paper presents a time dependent polarization constitutive model suitable for predicting nonlinear polarization and electro-mechanical strain responses of ferroelectric materials subject to various histories of electric fields. The constitutive model is derived based on a single integral form with nonlinear (electric field and temperature dependent) integrand. The total polarization consists of the time-dependent and residual components. The residual component of the polarization is due to polarization switching in the ferroelectric materials. We use an ‘internal clock’ concept to incorporate the effect of electric field on the rate of polarization. The corresponding strain response is determined through the use of third order piezoelectric constant and/or fourth order electrostrictive constant that vary with polarization stage. It is assumed that in absence of polarization, both piezoelectric and electrostrictive constants are zero. To incorporate the effect of temperature on the overall polarization behavior all material parameters in the constitutive model are allowed to change with the ambient temperature. We present numerical studies on the effect of time, temperature, and electric field on the response of ferroelectric material followed by verification of the constitutive model. Experimental data on lead zirconate titanate (PZT) materials available in the literature are used to verify the model.     

Mode III crack problems for two bonded functionally graded piezoelectric materials

This paper investigates the singular electromechanical field near the crack tips of an internal crack. The crack is perpendicular to the interface formed by bonding two half planes of different functionally graded piezoelectric material. The properties of two materials, such as elastic modulus, piezoelectric constant and dielectric constant, are assumed in exponential forms and vary along the crack direction. The singular integral equations for impermeable and permeable cracks are derived and solved by using the Gauss–Chebyshev integration technique. It shows that the stresses and electrical displacements around the crack tips have the conventional square root singularity. The stress intensity and electric displacement intensity factors are highly affected by the material nonhomogeneity parameters β and γ. The solutions for some degenerated problems can also be obtained.     

Postbuckling analysis of axially-loaded laminated cylindrical shells with piezoelectric actuators

A compressive postbuckling analysis is presented for a laminated cylindrical shell with piezoelectric actuators subjected to the combined action of mechanical, electric and thermal loads. The temperature field considered is assumed to be a uniform distribution over the shell surface and through the shell thickness, and the electric field is assumed to be the transverse component E_Z only. The material properties are assumed to be independent of the temperature and the electric field. The governing equations are based on the classical shell theory with von Kármán–Donnell-type kinematic nonlinearity. The nonlinear prebuckling deformations and initial geometric imperfections of the shell are both taken into account. A boundary layer theory of shell buckling, which includes the effects of nonlinear prebuckling deformations, large deflections in the postbuckling range, and initial geometric imperfections of the shell, is extended to the case of hybrid laminated cylindrical shells. A singular perturbation technique is employed to determine the buckling loads and postbuckling equilibrium paths. The numerical illustrations concern the compressive postbuckling behavior of perfect and imperfect, cross-ply laminated cylindrical thin shells with fully covered or embedded piezoelectric actuators under different sets of thermal and electric loading conditions. The effects played by temperature rise, applied voltage, shell geometric parameter, stacking sequence, as well as initial geometric imperfections are studied.     

Two-dimensional sliding frictional contact of functionally graded materials

A multi-layered model for sliding frictional contact analysis of functionally graded materials (FGMs) with arbitrarily varying shear modulus under plane strain-state deformation has been developed. Based on the fact that an arbitrary curve can be approached by a series of continuous but piecewise linear curves, the FGM is divided into several sub-layers and in each sub-layers the shear modulus is assumed to be a linear function while the Poisson's ratio is assumed to be a constant. In the contact area, it is assumed that the friction is one of Coulomb type. With this model the fundamental solutions for concentrated forces acting perpendicular and parallel to the FGMs layer surface are obtained. Then the sliding frictional contact problem of a functionally graded coated half-space is investigated. The transfer matrix method and Fourier integral transform technique are employed to cast the problem to a Cauchy singular integral equation. The contact stresses and contact area are calculated for various moving stamps by solving the equations numerically. The results show that appropriate gradual variation of the shear modulus can significantly alter the stresses in the contact zone.     

Torsional buckling and postbuckling of FGM cylindrical shells in thermal environments

A postbuckling analysis is presented for a functionally graded cylindrical shell subjected to torsion in thermal environments. Heat conduction and temperature-dependent material properties are both taken into account. The temperature field considered is assumed to be a uniform distribution over the shell surface and varied in the thickness direction. The material properties of functionally graded materials (FGMs) are assumed to be graded in the thickness direction according to a simple power law distribution in terms of the volume fractions of the constituents, and are assumed to be temperature-dependent. The governing equations are based on a higher order shear deformation theory with a von Kármán–Donnell-type of kinematic non-linearity. The non-linear prebuckling deformations and initial geometric imperfections of the shell are both taken into account. A singular perturbation technique is employed to determine the buckling load and postbuckling equilibrium paths. The numerical illustrations concern the postbuckling behavior of twist, perfect and imperfect, FGM cylindrical shells under different sets of thermal fields. The results reveal that the volume fraction distribution of FGMs has a significant effect on the buckling load and postbuckling behavior of FGM cylindrical shells subjected to torsion. They also confirm that the torsional postbuckling equilibrium path is weakly unstable and the shell structure is virtually imperfection–insensitive.     

Shear wave velocity in sand considering the effects of frequency based on particle contact theory

Since the shear waves involved in in-situ and laboratory measurement methods vary significantly in terms of the frequency range, it is necessary to consider the effects of frequency on the shear wave velocity. In this study, sand particles are assumed to be spherical solid particles with an equal radius and identical material properties, and sand skeletons are regarded as granular aggregations generated through the random packing of sand particles. It is also assumed that the sand particles only undergo elastic deformation during shear wave propagation. Based on a spherical particle model, a formula is obtained for calculating the shear wave velocity in sand, with the shear wave frequency as an extra influencing parameter. The quantitative calculations demonstrate that the shear wave velocity decreases with an increase of sand porosity, and accelerates with increases of vertical effective stress and elastic modulus of the sand particles. It is also indicated that both the particle density and Poisson’s ratio of the sand particles have negligible effects on the shear wave propagation. The frequency dispersion characteristics of shear wave propagating in sand are also discussed. Moreover, the critical frequency is defined and its analytical expression is derived. The calculation results obtained using the proposed equations agree well with the in-situ measurement results and bender element test data.     

Anti-plane shear crack normal to and terminating at the interface of two bonded piezoelectric ceramics

The electroelastic analysis of two bonded dissimilar piezoelectric ceramics with a crack perpendicular to and terminating at the interface is made. By using Fourier integral transform, the associated boundary value problem is reduced to a singular integral equation with generalized Cauchy kernel, the solution of which is given in closed form. Results are presented for a permeable crack under anti-plane shear loading and in-plane electric loading. Obtained results indicate that the electroelastic field near the crack tip in the homogeneous piezoelectric ceramic is dominated by a traditional inverse square-root singularity, while the electroelastic field near the crack tip at the interface exhibits the singularity of power law r^−α, r being distance from the interface crack tip and α depending on the material constants of a bi-piezoceramic. In particular, electric field has no singularity at the crack tip in a homogeneous solid, whereas it is singular around the interface crack tip. Numerical results are given graphically to show the effects of the material properties on the singularity order and field intensity factors.     

An analytical study on the nonlinear vibration of functionally graded beams

Nonlinear vibration of beams made of functionally graded materials (FGMs) is studied in this paper based on Euler-Bernoulli beam theory and von Kármán geometric nonlinearity. It is assumed that material properties follow either exponential or power law distributions through thickness direction. Galerkin procedure is used to obtain a second order nonlinear ordinary equation with quadratic and cubic nonlinear terms. The direct numerical integration method and Runge-Kutta method are employed to find the nonlinear vibration response of FGM beams with different end supports. The effects of material property distribution and end supports on the nonlinear dynamic behavior of FGM beams are discussed. It is found that unlike homogeneous beams, FGM beams show different vibration behavior at positive and negative amplitudes due to the presence of quadratic nonlinear term arising from bending-stretching coupling effect.     

Nonlinear axisymmetric thermomechanical response of piezo-FGM shallow spherical shells

Thermomechanical instability of shallow spherical shells made of functionally graded material (FGM) and surface-bonded piezoelectric actuators is studied in this paper. The governing equations are based on the classical shell theory of shells and the Sanders nonlinear kinematics equations. It is assumed that the property of the FGMs varies continuously through the thickness of the shell according to a power law distribution of the volume fraction of the constituent materials. The constituent materials of the functionally graded shell are assumed to be mixture of ceramic and metal. The analytical solutions are obtained for uniform external pressure, thermal loading, and constant applied actuator voltage.     

Stress-intensity factors for surface cracks in functionally graded materials under mode-I thermomechanical loading

This paper describes the development and application of a general domain integral method to obtain J-values along crack fronts in three-dimensional configurations of isotropic, functionally graded materials (FGMs). The present work considers mode-I, linear-elastic response of cracked specimens subjected to thermomechanical loading, although the domain integral formulation accommodates elastic–plastic behavior in FGMs. Finite element solutions and domain integral J-values for a two-dimensional edge crack show good agreement with available analytical solutions for both tension loading and temperature gradients. A displacement correlation technique provides pointwise stress-intensity values along semi-elliptical surface cracks in FGMs for comparison with values derived from the proposed domain integral. Numerical implementation and mesh refinement issues to maintain path independent J-values are explored. The paper concludes with a parametric study that provides a set of stress-intensity factors for semi-elliptical surface cracks covering a practical range of crack sizes, aspect ratios and material property gradations under tension, bending and spatially-varying temperature loads.     

Stochastic nonlinear bending response of laminated composite plates with system randomness under lateral pressure and thermal loading

In this paper, the effect of sensitivity of randomness in system parameters on the nonlinear transverse central deflection response of laminated composite plates subjected to transverse uniform lateral pressure and thermal loading is examined. System parameters such as the lamina material properties, expansion of thermal coefficients, lamina plate thickness and lateral load are modelled as basic random variables. A higher order shear deformation theory in the von-Karman sense is used to model the system behavior of the laminated plate. A direct iterative-based C ⁰ nonlinear finite element method in conjunction with the first-order perturbation technique developed by the authors is extended for thermal problem to obtain the second-order response statistics, i.e., mean and variance of the nonlinear transverse deflection of the plate. Typical numerical results of composite plates with temperature independent and dependent material properties subjected to uniform temperature and combination of uniform and transverse temperature are obtained for various combinations of geometric parameters, uniform lateral pressures, staking sequences and boundary conditions. The results have been compared with those available in the literature and an independent Monte Carlo simulation.     

Vibration analysis of FG annular sector in moderately thick plates with two piezoelectric layers

Free vibration of functionally graded(FG) annular sector plates embedded with two piezoelectric layers is studied with a generalized differential quadrature(GDQ)method. Based on the first-order shear deformation(FSD) plate theory and Hamilton's principle with parameters satisfying Maxwell's electrostatics equation in the piezoelectric layers, governing equations of motion are developed. Both open and closed circuit(shortly connected) boundary conditions on the piezoelectric surfaces, which are respective conditions for sensors and actuators, are accounted for. It is observed that the open circuit condition gives higher natural frequencies than a shortly connected condition. For the simulation of the potential electric function in piezoelectric layers, a sinusoidal function in the transverse direction is considered. It is assumed that properties of the FG material(FGM) change continuously through the thickness according to a power distribution law.The fast rate convergence and accuracy of the GDQ method with a small number of grid points are demonstrated through some numerical examples. With various combinations of free, clamped, and simply supported boundary conditions, the effects of the thicknesses of piezoelectric layers and host plate, power law index of FGMs, and plate geometrical parameters(e.g., angle and radii of annular sector) on the in-plane and out-of-plane natural frequencies for different FG and piezoelectric materials are also studied. Results can be used to predict the behaviors of FG and piezoelectric materials in mechanical systems.     

Two-dimensional contact mechanics of functionally graded materials with arbitrary spatial variations of material properties

A multi-layered model for frictionless contact analysis of functionally graded materials (FGMs) with arbitrarily varying elastic modulus under plane strain-state deformation has been developed. Based on the fact that an arbitrary curve can be approached by a series of continuous but piecewise linear curves, the FGM is divided into several sub-layers and in each sub-layers the shear modulus is assumed to be linear function while the Poisson’s ratio is assumed to be a constant. With the model, the frictionless contact problem of a functionally graded coated half-space is investigated. By using the transfer matrix method and Fourier integral transform technique, the problem is reduced to a Cauchy singular integral equation. The contact pressure, contact region and indentation are calculated for various indenters by solving the equations numerically.     

粘贴压电层功能梯度材料Timoshenko梁的热过屈曲分析

研究了上下表面粘贴压电层的功能梯度材料Timoshenko梁在升温及电场作用下的过屈曲行为。在精确考虑轴线伸长和一阶横向剪切变形的基础上,建立了压电功能梯度Timoshenko层合梁在热-电-机械载荷作用下的几何非线性控制方程。其中,假设功能梯度的材料性质沿厚度方向按照幂函数连续变化,压电层为各向同性均匀材料。采用打靶法数值求解所得强非线性边值问题,获得了在均匀电场和横向非均匀升温场内两端固定Timoshenko梁的静态非线性屈曲和过屈曲数值解。并给出了梁的变形随热、电载荷及材料梯度参数变化的特性曲线。结果表明,通过施加电压在压电层产生拉应力可以有效地提高梁的热屈曲临界载荷,延缓热过屈曲发生。由于材料在横向的非均匀性,即使在均匀升温和均匀电场作用下,也会产生拉-弯耦合效应。但是对于两端固定的压电-功能梯度材料梁,在横向非均匀升温下过屈曲变形仍然是分叉形的。     

Stress analysis of functionally graded rotating discs: analytical and numerical solutions

This study deals with stress analysis of annular rotating discs made of functionally graded materials(FGMs).Elasticity modulus and density of the discs are assumed to vary radially according to a power law function,but the material is of constant Poisson’s ratio.A gradient parameter n is chosen between 0 and 1.0.When n = 0,the disc becomes a homogeneous isotropic material.Tangential and radial stress distributions and displacements on the disc are investigated for various gradient parameters n by means of the diverse elasticity modulus and density by using analytical and numerical solutions.Finally,a homogenous tangential stress distribution and the lowest radial stresses along the radius of a rotating disc are approximately obtained for the gradient parameter n = 1.0 compared with the homogeneous,isotropic case n = 0.This means that a disc made of FGMs has the capability of higher angular rotations compared with the homogeneous isotropic disc.     

Stochastic Response of Laminated Composite Shell Panel in Hygrothermal Environment

In this article, the effect of random system properties is accounted to estimate the free vibration of multilayered composite shell panel in hygrothermal environment. The majority of previous investigations assumed that the material properties are independent of temperature and moisture. Establishing margin on material properties for design, analysis, and conceptualizing the material is highly difficult and uncertain. To perturb that, a first-order perturbation technique is adopted to obtain the response statistics of the structure by obtaining the mean and variance of random natural frequencies. The higher-order shear deformation theory with 9 degrees of freedom per node with the von-Karman sense of nonlinear kinematics is employed for generating basic formulation. The analysis is carried out by using quadratic C⁰ eight-noded isoparametric element. The governing equation for free vibration of laminated composite panel is derived using variational principle, which is a generalization of the principle of virtual displacement. The solution methodology is validated with published results. Mean and variance are obtained for different cross-ply of spherical and cylindrical laminates with different stacking sequence and boundary conditions.