Analytical solutions for plane problem of functionally graded magnetoelectric cantilever beam

In this paper, an exact analytical solution is presented for a transversely isotropic functionally graded magneto-electro-elastic (FGMEE) cantilever beam, which is subjected to a uniform load on its upper surface, as well as the concentrated force and moment at the free end. This solution can be applied for any form of gradient distribution. For the basic equations of plane problem, all the partial differential equations governing the stress field, electric, and magnetic potentials are derived. Then, the expressions of Airy stress, electric, and magnetic potential functions are assumed as quadratic polynomials of the longitudinal coordinate. Based on all the boundary conditions, the exact expressions of the three functions can be determined. As numerical examples, the material parameters are set as exponential and linear distributions in the thickness direction. The effects of the material parameters on the mechanical, electric, and magnetic fields of the cantilever beam are analyzed in detail.     

Static analysis for functionally graded piezoelectric actuators or sensors under a combined electro-thermal load

Based on the theory of piezoelasticity, a functionally graded piezoelectric sandwich cantilever under an applied electric field and/or a heat load is studied. All materials may be arbitrary functional gradients in the thickness direction. The static solution for the mentioned problems is presented by the Airy stress function method. As a special case, assuming that the material composition varies continuously in the direction of the thickness according to a power law distribution, a comprehensive parametric study is conducted to show the influence of electromechanical coupling (EMC), functionally graded index, temperature change and thickness ratio on the bending behavior of actuators or sensors. The distribution of electric field or normal stress in present FGPM actuators is continuous along the thickness, which overcomes the problem in traditional layered actuators. The solution facilitates the design optimization for different piezoelectric actuators and has another potential application for material parameter identification.     

Static analysis of anisotropic functionally graded magneto-electro-elastic beams subjected to arbitrary loading

This paper considers the analytical and semi-analytical solutions for anisotropic functionally graded magneto-electro-elastic beams subjected to an arbitrary load, which can be expanded in terms of sinusoidal series. For the generalized plane stress problem, the stress function, electric displacement function and magnetic induction function are assumed to consist of two parts, respectively. One is a product of a trigonometric function of the longitudinal coordinate (x) and an undetermined function of the thickness coordinate (z), and the other a linear polynomial of x with unknown coefficients depending on z. The governing equations satisfied by these z-dependent functions are derived. The analytical expressions of stresses, electric displacements, magnetic induction, axial force, bending moment, shear force, average electric displacement, average magnetic induction, displacements, electric potential and magnetic potential are then deduced, with integral constants determinable from the boundary conditions. The analytical solution is derived for beam with material coefficients varying exponentially along the thickness, while the semi-analytical solution is sought by making use of the sub-layer approximation for beam with an arbitrary variation of material parameters along the thickness. The present analysis is applicable to beams with various boundary conditions at the two ends. Two numerical examples are presented for validation of the theory and illustration of the effects of certain parameters.     

Magnetic field and magneto elastic stress in an infinite plate containing an elliptical hole with an edge crack under uniform electric current

Two-dimensional solutions of the electric current, magnetic field and magneto elastic stress are presented for a magnetic material of a thin infinite plate containing an elliptical hole with an edge crack under uniform electric current. Using a rational mapping function, the each solution is obtained as a closed form. The linear constitutive equation is used for the magnetic field and the stress analyses. According to the electro-magneto theory, only Maxwell stress is caused as a body force in a plate which raises a plane stress state for a thin plate and the deformation of the plate thickness. Therefore the magneto elastic stress is analyzed using Maxwell stress. No further assumption of the plane stress state that the plate is thin is made for the stress analysis, though Maxwell stress components are expressed by nonlinear terms. The rigorous boundary condition expressed by Maxwell stress components is completely satisfied without any linear assumptions on the boundary. First, electric current, magnetic field and stress analyses for soft ferromagnetic material are carried out and then those analyses for paramagnetic and diamagnetic materials are carried out. It is stated that the stress components are expressed by the same expressions for those materials and the difference is only the magnitude of the permeability, though the magnetic fields H_x, H_y are different each other in the plates. If the analysis of magnetic field of paramagnetic material is easier than that of soft ferromagnetic material, the stress analysis may be carried out using the magnetic field for paramagnetic material to analyze the stress field, and the results may be applied for a soft ferromagnetic material. It is stated that the stress state for the magnetic field H_x, H_y is the same as the pure shear stress state. Solving the present magneto elastic stress problem, dislocation and rotation terms appear, which makes the present problem complicate. Solutions of the magneto elastic stress are nonlinear for the direction of electric current. Stresses in the direction of the plate thickness are caused and the solution is also obtained. Figures of the magnetic field and stress distribution are shown. Stress intensity factors are also derived and investigated for the crack length and the electric current direction.     

Fatigue crack growth induced by domain switching under electromechanical load in ferroelectrics

Reliability calls for a better understanding of the failure of ferroelectric ceramics. The fracture and fatigue of ferroelectric ceramics under an electric field or a combined electric and mechanical loading are investigated. The small-scale domain-switching model is modified to analyze failure due to fracture and fatigue. Effects of anisotropy and electromechanical load coupling are taken into account. Analytical expressions are obtained for domain-switching regions near the crack tip such that of 90° domain switching can be distinguished from 180° domain switching in addition to different initial poling directions. The crack tip stress intensity variation of ferroelectric ceramics due to the domain switching is analyzed. A positive electric field tends to enhance the propagation of an insulating crack perpendicular to the poling direction, while a negative field impedes it. Fatigue crack growth under various coupling loads and effects of the stress field and electric field on near field stress intensity variation are analyzed. Predicted crack growth versus cyclic electric field agrees well with experiment.     

3D elasticity solution for static analysis of functionally graded piezoelectric annular plates on elastic foundations using SSDQM

Three-dimensional elasticity solution for static analysis of functionally graded piezoelectric (FGP) annular plates with and without elastic foundations through using state-space based differential quadrature method (SSDQM) at different boundary conditions is presented in this paper. The material properties are assumed to have an exponent-law variation along the thickness. A semi-analytical approach which makes use of state-space method in thickness direction and one-dimensional differential quadrature method in radial direction is utilized to obtain the mechanical behavior of FGP annular plates. The state variables include a combination of electric potential, electric displacement, three mechanical displacement parameters and three stress parameters. Numerical results are given to demonstrate the convergency and accuracy of the present method. Both closed circuit and open circuit effects are studied and the influences of the Winkler and shearing layer elastic coefficients of the foundations, the material property graded index, radius, thickness, mechanical load and boundary conditions on the deflection response of the FGP annular plates are investigated. The new results can be used as a benchmark solutions for future researches.     

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.     

Elliptical crack normal to functionally graded interface of bonded solids

This paper presents an analysis of an elliptical crack that is perpendicular to a functionally graded interfacial zone between two fully bonded solids. The functionally graded interfacial zone is treated as a non-homogeneous solid layer with its elastic modulus varying in the thickness direction. A generalized Kelvin solution based boundary element method is employed for the calculation of the stress intensity factors associated with the three-dimensional crack problem. The elliptical crack surface is subject to either uniform normal traction or uniform shear traction. The stress intensity factors are examined by taking into account the effects of the non-homogeneity parameter and thickness of the functionally graded interfacial zone, as well as the crack distance to the zone. The SIF values are further incorporated into the S-criterion for prediction of crack growth. The paper presents the most possible direction and location of the elliptical crack growth under an inclined tensile (or compressive) load. The paper further presents results of the critical external loads that would cause the elliptical crack to grow at the most possible location and along the most possible direction. The paper also examines the effects of external load direction and material and geometrical parameters on the critical loads.     

A Simplified Treatise of the Scott Bond Testing Method

The Scott bond test method has been used extensively in the paper industry over the years as a means to assess the bond strength of paper. The method has been a subject of some controversy lately since it does not always correlate to the sensitivity of the material to fracture by delamination. To gain some further insight into which parameters govern the fracture process in a Scott bond test a simplified approach has been chosen in order to formulate an analytical mathematical/mechanical model of the test. The model is dynamic in the sense that inertia effects are included. The material model utilised is a simple cohesive theory that assumes a linear behaviour between stress and crack opening when the material has started to degrade. This choice of material model makes the mathematical model very nonlinear. In fact, a system of three coupled nonlinear second order partial differential equations have to be solved and adjusted to the correct initial conditions. The material parameters needed for the model are the elastic modulus in the thickness direction, the transverse shear (elastic) modulus, the tensile strength (in the thickness direction) and the fracture work (per unit area) for a delamination crack. To investigate the ability of the model, a Scott bond testing apparatus have been equipped with a piezoelectric load sensor. The load cell was mounted on the apparatus’ pendulum so that the load acting on the sample holder could be recorded during the whole impact stage. This was done for a number of different initial velocities of the pendulum and it is found that the model gives a fair prediction of the contact load.     

Propagation behavior of Love waves in a functionally graded half-space with initial stress

The propagation behavior of Love waves in a functionally graded material layered non-piezoelectric half-space with initial stress is taken into account. The Wentzel–Kramers–Brillouin (WKB) technique is adopted for the theoretical derivations. The analytical solutions are obtained for the dispersion relations and the distributions of the mechanical displacement and stress along the thickness direction in the layered structure. First, these solutions are used to study the effects of the initial stress on the dispersion relations and the group and phase velocities, then the influences of the initial stress on the distributions of the mechanical displacement and shear stresses along the thickness direction are discussed in detail. Numerical results obtained indicate that the phase velocity of the Love waves increases with the increase in the magnitude of the initial tensile stress, while decreases with the increase in the magnitude of the initial compression stress. The effects on the dispersion relations of the Love wave propagation are negligible as the magnitudes of the initial stress are less than 100 MPa. Some other results are obtained for the distributions of field quantities along thickness direction. The results obtained are not only meaningful for the design of functionally graded structures with high performance but also effective for the evaluation of residual stress distribution in the layered structures.     

Interplay between polarization rotation and crack propagation in PMN-PT relaxor single crystals

Investigations on the interconnection between the polarization rotation and crack propagation are performed for [110]-oriented 74Pb(Mg_1/3Nb_2/3)O₃-26PbTiO₃ relaxor ferroelectric single crystal under electric loadings along [001] direction. The crystal is of predominantly monoclinic M_A phase with scatter distributed rhombohedral (R) phase under a moderate poling field of 900 V/mm in [001] direction. With magnitude of 800 V/mm, a through thickness crack is initiated near the electrode by electric cycling. Static electric loadings is then imposed to the single crystal. As the applied static electric field increases, domain switching in the monoclinic M_A phase and phase transition from M_A to R phase occur near the crack. The results indicate that the crack features a conducting one. Whether domain switching or phase transition occurs depends on the intensity of the electric field component that is perpendicular to the applied electric field.     

Three dimensional phase field study on the thickness effect of ferroelectric polymer thin film

The electromechanical behavior of poly(vinylidene fluoride-trifluoroethylene) [P(VDF-TrFE)] ferroelectric thin film was investigated using the three dimensional (3D) phase-field method. Various energetic contributions, including elastic, electrostatic, and domain wall energy were taken into account in the variational functional of the phase field model. Evolution of the microscopic domain structures of P (VDF-TrFE) polymer film was simulated. Effects of the in-plane residual stress, the film thickness and externally applied electric bias field on the electromechanical properties of the film were explored. The obtained numerical results showed that the macroscopic responses of the electric hysteresis loops are sensitive to the residual stress and electric bias field. It was also found that thickness has a great effect on the electric hysteresis loops and remanent polarization.     

Magneto elastic stress subjected to uniform magnetic field over a thin infinite plate containing an elliptical hole with an edge crack

Two dimensional solutions of the magnetic field and magneto elastic stress are presented for a magnetic material of a thin infinite plate containing an elliptical hole with an edge crack subjected to uniform magnetic field. Using a rational mapping function, each solution is obtained as a closed form. The linear constitutive equation is used for these analyses. According to the electro-magneto theory, only Maxwell stress is caused as a body force in a plate. In the present paper, it raises a plane stress state for a thin plate, the deformation of the plate thickness and the shear deflection. Therefore the magneto elastic stress is analyzed using Maxwell stress. No further assumption of the plane stress state that the plate is thin is made for the stress analysis, though Maxwell stress components are expressed by nonlinear terms. The rigorous boundary condition expressed by Maxwell stress components is completely satisfied without any linear assumptions on the boundary. First, magnetic field and stress analyses for soft ferromagnetic material are carried out and then those analyses for paramagnetic and diamagnetic materials are carried out. It is stated that those plane stress components are expressed by the same expressions for those materials and the difference is only the magnitude of the permeability, though the magnetic fields H_x, H_y are different each other in the plates. If the analysis of magnetic field of paramagnetic material is easier than that of soft ferromagnetic material, the stress analysis may be carried out using the magnetic field for paramagnetic material to analyze the stress field, and the results may be applied for a soft ferromagnetic material. It is stated that the stress state for the magnetic field H_x, H_y is the same as the pure shear stress state. Solutions of the magneto elastic stress are nonlinear for the direction of uniform magnetic field. Stresses in the direction of the plate thickness and shear deflection are caused and the solutions are also obtained. Figures of the magnetic field and stress distribution are shown. Stress intensity factors are also derived and investigated for the crack length.     

Piezoelectric study on circularly cylindrical layered media

In this paper, an analytical solution in series form for the problem of a circularly cylindrical layered piezoelectric composite consisting of N dissimilar layers is presented within the framework of linear piezoelectricity. Each layer of the composite is assumed to be transversely isotropic with respect to the longitudinal direction (x₃ direction), and the composite is subject to arbitrary electromechanical singularities infinitely extended in a direction perpendicular to the x₁–x₂ plane such that only in-plane electric fields and out-of-plane displacement are produced. The alternating technique in conjunction with the method of analytical continuation is applied to derive the general multilayered media solution in an explicit series form, whose convergence is guaranteed numerically. The distributions of the shear stress and electric field are found to be dependent on the material combinations and the magnitude and position of the electromechanical singularities. An exactly closed form solution is obtained and discussed graphically for a practical example.     

Stress-induced butterfly and square-like magnetostriction loops transition

The stress-induced magnetic domain switching in Fe Ga thin films is studied using phase-field method.In particular,the magnetic field is applied along the [110] direction and biaxial stresses are applied along[100] and [010].A compressive pre-stress corresponds to a smaller coercive magnetic field while a tensile pre-stress corresponded to a larger coercive field.At the same time,it is also found that the transition between butterfly and square-like magnetostriction loops occurs at the critical opposite biaxial stress state.The two different evolutions correspond to two different mechanisms:one is that the single domain swings across a fan area back and forth;the other is that the single domain turns a clockwise circle.The results can be explained by the stress tuned anisotropy energy well.     

铁电陶瓷PZT53复杂力电耦合行为的实验研究

通过实验研究了平行和垂直于极化方向的正应力对铁电陶瓷锆钛酸铅(PZT53)的电滞回线(E3-P3)和电致应变曲线(E-ε)的影响. 实验发现平行于极化方向的压应力对PZT53陶瓷的电滞回线、电致应变曲线形状以及矫顽场大小都有明显的影响,但是垂直于极化方向的拉、压应力只对PZT53陶瓷的电致应变曲线形状有明显的影响,但对电滞回线形状和矫顽场大小都没有显著影响. 采用畴翻转的模型详细解释了观察到的实验现象,所得结果为建立铁电陶瓷的多轴力、电耦合本构模型,提供了物理基础.     

Three dimensional dynamic electromechanical constitutive relations for ferroelectric materials

We propose, in this paper, three dimensional constitutive relations suitable for describing the dynamic electromechanical responses of ferroelectric materials. Our approach is phenomenological and is based on the assumptions that the macroscopic electrical properties can be represented by electric dipoles whose magnitude and direction can be altered by external stimuli. The constitutive relations for the stress and the electric displacement take into account the transient and instantaneous responses of the dipole moments and the transient response of domain switching. A specific respesentation of the theory is also derived; this representation is suitable for studying the responses of ferroelectic materials to dynamic mechanical loading and rapidly applied electric field.     

Dielectric crack problem for a magnetoelectroelastic strip with functionally graded properties

The paper presents a fracture analysis for an electromagnetically dielectric crack in a functionally graded magnetoelectroelastic strip. It is considered that the material properties are varying exponentially along the width direction. Under the assumption of the in-plane magneto-electro-mechanical loadings, the dielectric crack is simulated by using the semi-permeable crack-face boundary conditions. The Fourier transform technique is applied to solve the boundary-value problem and four coupling singular integral equations are determined. A nonlinear system of algebraic equations is further derived and solved numerically to determine the electromagnetic field inside the crack. Then the field intensity factors of stress, electric displacement, and magnetic induction are given. Through the numerical computations, the effects of the material non-homogeneity and the permeability of crack interior on the electric displacement and the magnetic induction at the crack faces are studied. The variations of the intensity factors of stress, electric displacement, and magnetic induction versus the geometry of the crack, the strip width, and the material non-homogeneity are presented in graphics respectively.     

风机叶片涂层雨蚀研究

研究了风电叶片在雨水冲击下的失效过程,利用有限元方法建立雨水冲击的数值模型,分析了雨滴的冲击角度、直径、涂层性能等对涂层受力的影响,并对双雨滴及四雨滴耦合模型进行了研究。结果表明:风雨耦合流场中叶片最大载荷在叶尖处,涂层的冲击受力与雨滴冲击角度和雨滴直径的大小呈正相关。低模量涂层能大幅降低涂层受冲击时的拉伸应力,小幅增加其压缩应力。同一平面内的双雨滴耦合时,拉伸应力受影响的范围较大,压缩应力受影响的范围较小。不同高度差下的双雨滴耦合中,随着高度差的增大,材料拉伸和压缩应力都呈现先增后减的趋势,其中压缩应力随高度差的变化幅度明显大于拉伸应力。四雨滴耦合时压缩应力对间距的变化更加敏感。涂层失效的过程主要为:涂层表面拉应力疲劳产生微裂纹;液压渗透引起裂纹扩展,涂层质量减少,表面粗糙度增加;表面小范围破坏后形成的凹面结构促进涂层的失效。     

Fracture analysis of circular-arc interface cracks in piezoelectric materials

The anti-plane problem of N arc-shaped interfacial cracks between a circular piezoelectric inhomogeneity and an infinite piezoelectric matrix is investigated by means of the complex variable method. Cracks are assumed to be permeable and then explicit expressions are presented, respectively, for the electric field on the crack faces, the complex potentials in media and the intensity factors near the crack-tips. As examples, the corresponding solutions are obtained for a piezoelectric bimaterial system with one or two permeable arc-shaped interfacial cracks, respectively. Additionally, the solutions for the cases of impermeable cracks also are given by treating an impermeable crack as a particular case of a permeable crack. It is shown that for the case of permeable interfacial cracks, the electric field is jumpy ahead of the crack tips, and its intensity factor is always dependent on that of stress. Moreover all the field singularities are dependent not only on the applied mechanical load, but also on the applied electric load. However, for the case of a homogeneous material with permeable cracks, all the singular factors are related only to the applied stresses and material constants.