Magneto-electro-elastic antiplane analysis of a bimaterial BaTiO3–CoFe2O4 composite wedge with an interface crack

The antiplane analysis is made for a bimaterial BaTiO₃–CoFe₂O₄ composite wedge containing an interface crack. The coupled magneto-electro-elastic field is induced by the piezoelectric/piezomagnetic BaTiO₃–CoFe₂O₄ composite materials. For the crack problems, the intensity factors of stress, strain, electric displacement, electric field, magnetic induction and magnetic field at crack tips are derived analytically. Also, the energy density criterion is applied to predict the fracture behavior of the interface crack. The numerical results also show that the energy release rate for a crack in a single wedge is negative.     

Analysis method of planar cracks of arbitrary shape in the isotropic plane of a three-dimensional transversely isotropic magnetoelectroelastic medium

The hyper-singular boundary integral equation method of crack analysis in three-dimensional transversely isotropic magnetoelectroelastic media is proposed. Based on the fundamental solutions or Green’s functions of three-dimensional transversely isotropic magnetoelectroelastic media and the corresponding Somigliana identity, the boundary integral equations for a planar crack of arbitrary shape in the plane of isotropy are obtained in terms of the extended displacement discontinuities across crack faces. The extended displacement discontinuities include the displacement discontinuities, the electric potential discontinuity and the magnetic potential discontinuity, and correspondingly the extended tractions on crack face represent the conventional tractions, the electric displacement and the magnetic induction boundary values. The near crack tip fields and the intensity factors in terms of the extended displacement discontinuities are derived by boundary integral equation approach. A solution method is proposed by use of the analogy between the boundary integral equations of the magnetoelectroelastic media and the purely elastic materials. The influence of different electric and magnetic boundary conditions, i.e., electrically and magnetically impermeable and permeable conditions, electrically impermeable and magnetically permeable condition, and electrically permeable and magnetically impermeable condition, on the solutions is studied. The crack opening model is proposed to consider the real crack opening and the electric and magnetic fields in the crack cavity under combined mechanical-electric-magnetic loadings. An iteration approach is presented for the solution of the non-linear model. The exact solution is obtained for the case of uniformly applied loadings on the crack faces. Numerical results for a square crack under different electric and magnetic boundary conditions are displayed to demonstrate the proposed method.     

An interface crack with contact zones in a piezoelectric/piezomagnetic bimaterial

An interface crack with a frictionless contact zone at the right crack tip between two semi-infinite piezoelectric/piezomagnetic spaces under the action of a remote mechanical loading, magnetic and electric fluxes as well as concentrated forces at the crack faces is considered. Assuming that all fields are independent on the coordinate x ₂ co-directed with the crack front, the stresses, the electrical and the magnetic fluxes as well as the derivatives of the jumps of the displacements, the electrical and magnetic potentials are presented via a set of analytic functions in the (x ₁, x ₃)-plane with a cut along the crack region. Two cases of magneto-electric conditions at the crack faces are considered. The first case assumes that the crack is electrically and magnetically permeable, and in the second case the crack is assumed electrically permeable while the open part of the crack is magnetically impermeable. For both these cases due to the above-mentioned representation the combined Dirichlet–Riemann boundary value problems have been formulated and solved exactly. Stress, electric and magnetic induction intensity factors are found in a simple analytical form. Transcendental equations and a closed form analytical formula for the determination of the real contact zone length have been derived for both cases of magnetic conditions in the crack region. For a numerical illustration of the obtained results a bimaterial BaTiO₃–CoFe₂O₄ with different volume fractions of BaTiO₃ has been used, and the influence of the mechanical loading and the intensity of the magnetic flux upon the contact zone length and the associated intensity factors as well as the energy release rate has been shown.     

Dynamic analysis of a penny-shaped dielectric crack in a magnetoelectroelastic solid under impacts

The transient response of a magneto-electro-elastic material with a penny-shaped dielectric crack subjected to in-plane magneto-electro-mechanical impacts is made. To simulate an opening crack with a dielectric interior, the crack-face electromagnetic boundary conditions are supposed to depend on the crack opening displacement and the jumps of electric and magnetic potentials across the crack. Four ideal crack-face electromagnetic boundary conditions involving a combination of electrically permeable or impermeable and magnetically permeable or impermeable assumptions can be reduced. The Laplace and Hankel transform techniques are further utilized to solve the mixed initial-boundary-value problem. Three coupling Fredholm integral equations are obtained and solved by the composite Simpson's rule. Dynamic field intensity factors of stress, electric displacement, magnetic induction, crack opening displacement (COD), electric potential and magnetic potential are given in the Laplace transform domain. By means of a numerical inversion of the Laplace transform, numerical results are calculated to show the variations of the physical parameters of concern versus the normalized time in graphics. The effects of applied electric and magnetic loads on the dynamic intensity factors of stress and COD, and the dynamic energy release rate for a BaTiO₃-CoFe₂O₄ composite with a penny-shaped vacuum crack are discussed in detail.     

Transient response of an interfacial crack between dissimilar magnetoelectroelastic layers under magnetoelectromechanical impact loadings: Mode-I problem

The dynamic response of an interfacial crack between two dissimilar magnetoelectroelastic layers is investigated under magnetic, electrical and mechanical impact loadings. Four kinds of ideal crack-face assumptions, i.e., magnetoelectrically impermeable (Case 1), magnetically impermeable and electrically permeable (Case 2), magnetically permeable and electrically impermeable (Case 3) and magnetoelectrically permeable (Case 4), are adopted separately. The dynamic field intensity factors and energy release rates are derived. The effects of loading combinations and crack configurations especially for the former on the dynamic response are examined according to energy release rate criterion. The numerical results show that, among others, a negative magnetic (or electrical) loading is generally prone to inhibit the crack extension rather than a positive one for a magnetically (or electrically) impermeable interfacial crack. Results presented in this paper should have potential applications to the design of multilayered magnetoelectroelastic structures.     

Transient response of a magnetoelectroelastic solid with two collinear dielectric cracks under impacts

Dynamic analysis of two collinear electro-magnetically dielectric cracks in a piezoelectromagnetic material is made under in-plane magneto-electro-mechanical impacts. Generalized semi-permeable crack-face boundary conditions are proposed to simulate realistic opening cracks with dielectric. Ideal boundary conditions of a combination of electrically permeable or impermeable and magnetically permeable or impermeable assumptions are several limiting cases of the semi-permeable dielectric crack. Utilizing the Laplace and Fourier transforms, the mixed initial-boundary-value problem is reduced to solving singular integral equations with Cauchy kernel. Dynamic intensity factors of stress, electric displacement, magnetic induction and crack opening displacement (COD) near the inner and outer crack tips are determined in the Laplace transform domain. Numerical results for a special magnetoelectroelastic solid are calculated to show the influences of the dielectric permittivity and magnetic permeability inside the cracks on the crack-face electric displacement and magnetic induction. By means of a numerical inversion of the Laplace transform, the variations of the normalized intensity factors of stress and COD are discussed against applied magnetoelectric impact loadings and the geometry of the cracks for fully impermeable, vacuum, fully permeable cracks and shown in graphics.     

Magnetoelectroelastic composite with poling parallel to plane of line crack under out-of-plane deformation

The three-dimensional field equations can in general be regarded as the sum of in-plane and out-of-plane deformation. The method for the general solution is the same for both although the boundary conditions could make a difference. If a particular solution in exact form may be found for the out-of-plane case, the same may not hold for the in-plane case. Hence, there may be a good reason for discussing the out-of-plane crack problem in certain situations that should be emphasized. Otherwise, the reason may lie in the exploration of possible application to the in-plane problem, a direct solution of which would have required a considerable effort. The contribution of this work rests on the new findings for the case of poling parallel to the crack in a magnetoelectroelastic composite made of BaTiO₃–CoFe₂O₄. The inclusions are BaTiO₃ and the matrix is CoFe₂O₄. Several new features of the solution were not expected before hand.Unlike in-plane deformation with poling normal to the crack plane, maximum crack growth enhancement is found to occur in the BaTiO₃–CoFe₂O₄ composite for a volume fraction of about 50%. Crack retardation increases as the volume fraction of the inclusions either increase or decrease. The occurrence of this same phenomenon in Mode I and II remain to be investigated. Poling direction of magnetic and electric field for line defects can have a significant effect on crack growth for magnetoelectroelastic materials. The foregoing conclusions are based on predictions made from the strain energy density criterion.     

Exact solution for elliptical inclusion in magnetoelectroelastic materials

This paper considers the magneto-electro-mechanical coupling between an inclusion and matrix, which are both of magnetoelectroelastic materials. The general cases including the mode I, mode II and mode III are studied. Analytical solutions for an elliptical cylinder inclusion inside an infinite magnetoelectroelastic medium under combined mechanical–electrical–magnetic loadings are formulated via the Stroh formalism. Crack problem is also investigated and the stress, electric and magnetic fields in the vicinity of the crack tip are determined by a complex vector of intensity factors. Various special cases, including an impermeable inclusion, a permeable crack, a rigid and permeable inclusion, a rigid and permeable line inclusion, a permeable cavity, and an impermeable cavity are discussed.     

Magnetic and electric poling effects associated with crack growth in BaTiO3–CoFe2O4 composite

Magnetoelectroelastic composite possesses the dual feature that the application of magnetic field induces electric polarization and electric field induces magnetization. The poling directions introduced magnetically and electrically can be different in addition to those for the applied magnetic and electric field. Their choices can influence the character of crack growth which could be enhanced or retarded. The details of how the directions of poling and applied field would affect crack initiation and growth are discussed in relation to the volume fraction of inclusions for a BaTiO₃–CoFe₂O₄ two phase composite. The multi-functional aspects of magnetoelectroelastic materials are involved since they entail multi-scaling features. Failure criteria that applies to isotropic elastic materials may not hold for composites exhibiting piezomagnetic and piezoelectric properties. For instance, a negative energy release rate has been obtained for cracks in piezoelectric materials.In view of what has been said with reference to the energy release rate approach, it is desirable to use the strain energy density function as a failure criterion, even if it is only for its positive definiteness character. Physically speaking, it is attractive to have a function that could rank the proportion of energy related to volume and shape change. They determine the proportion of the hard and soft phase of the composite and hence the volume fraction of the constituent. Strength and toughness parameters used for ranking isotropic and homogeneous materials will not apply for anisotropic and/or nonhomogeneous materials if these microstructure effects could not be suppressed to a lower scale and represented as an average at the macroscopic scale. Too much emphases cannot be placed on the need to clarify the multi-scaling aspects of piezoelectric and piezomagnetic materials. Their behavior as affected by the presence of crack-like defects should be understood prior to deciding whether the material characterization approach would be suitable. That is whether simplicity could justify at the expense of conceptual rigor. Much of this would depend on scaling the time and size related to loading and material structure interaction. The magnetoelectroelastic crack model selected in the work to follow perhaps will provide an insight into the complexicity of the state of affairs for treating the finer details of material behavior with rigor.The proposed test model shows that crack growth in the magnetoelectroelastic materials can be suppressed by increasing the magnitude of the piezomagnetic constants in relation to those for piezoelectricity. A more rational means of evaluating the resistance of materials against fracture is thus proposed, particularly when anisotropy and inhomogeneity might be present.     

Strip electric-magnetic breakdown model in a magnetoelectroelastic medium

Extending the polarization saturation model [Gao et al., 1997. Local and global energy release rates for an electrically yielded crack in a piezoelectric ceramic. J. Mech. Phys. Solids 45, 491-510] and the dielectric breakdown (DB) model [Zhang et al., 2005. The strip dielectric breakdown model. Int. J. Fract. 132, 311-327] in piezoelectric materials, the Strip Electric-Magnetic Breakdown (SEMB) model is proposed for electrically and magnetically impermeable crack in a magnetoelectroelastic medium to study the effect of the nonlinear character of electric field and magnetic field on fracture of magnetoelectroelastic materials. In the SEMB model, the electric field in the strip of the electric breakdown zone ahead of the crack tip is equal to the electric breakdown strength, while the magnetic filed in the strip of the magnetic breakdown zone is equal to the magnetic breakdown strength. By using the extended Stroh formalism and the extended dislocation modeling of a crack, the Griffith crack problem under the electrically and magnetically elastic-plastic condition in a magnetoelectroelastic medium is reduced to a set of dual integral equations. The sizes of the electric breakdown zone and the magnetic breakdown zone, the extended intensity factors and the local J-integral are obtained. The effect of the combined mechanical-electric-magnetic loadings on the local J-integral is studied.     

Discussion on electromagnetic crack face boundary conditions for the fracture mechanics of magneto-electro-elastic materials

This paper discusses electromagnetic boundary conditions on crack faces in magneto- electroelastic materials, where piezoelectric, piezomagnetic and magnetoelectric effects are coupled. A notch of finite thickness in these materials is also addressed. Four idealized electromagnetic boundary conditions assumed for the crack-faces are separately investigated, i.e. (a) electrically and magnetically impermeable (crack-face), (b) electrically impermeable and magnetically permeable, (c) electrically permeable and magnetically impermeable, and (d) electrically and magnetically permeable. The influence of the notch thickness on important parameters, such as the field intensity factors, the energy release rate at the notch tips and the electromagnetic fields inside the notch, are studied and the results are obtained in closed-form. Results under different idealized electromagnetic boundary conditions on the crack-face are compared, and the applicability of these idealized assumptions is discussed.The project supported by the National Natural Science Foundation of China (10102004) The English text was polished by Yunming Chen.     

Piezomagnetic and piezoelectric poling effects on mode I and II crack initiation behavior of magnetoelectroelastic materials

The ferrite and ferroelectric phase of magnetoelectroelastic (MEE) material can be selected and processed to control the macroscopic behavior of electron devices using continuum mechanics models. Once macro- and/or microdefects appear, the highly intensified magnetic and electric energy localization could alter the response significantly to change the design performance. Alignment of poling directions of piezomagnetic and piezoelectric materials can add to the complexity of the MEE material behavior to which this study will be concerned with.Appropriate balance of distortional and dilatational energy density is no longer obvious when a material possesses anisotropy and/or nonhomogeneity. An excess of the former could result in unwanted geometric change while the latter may lead to unexpected fracture initiation. Such information can be evaluated quantitatively from the stationary values of the energy density function dW/dV. The maxima and minima have been known to coincide, respectively, with possible locations of permanent shape change and crack initiation regardless of material and loading type. The direction of poling with respect to a line crack and the material microstructure described by the constitutive coefficients will be specified explicitly with reference to the applied magnetic field, electric field and mechanical stress, both normal and shear. The crack initiation load and direction could be predicted by finding the direction for which the volume change is the largest. In contrast to intuition, change in poling directions can influence the cracking behavior of MEE dramatically. This will be demonstrated by the numerical results for the BaTiO₃–CoFe₂O₄ composite having different volume fractions where BaTiO₃ and CoFe₂O₄ are, respectively, the inclusion and matrix.To be emphasized is that mode I and II crack behavior will not have the same definition as that in classical fracture mechanics where load and crack extension symmetry would coincide. A striking result is found for a mode II crack. By keeping the magnetic poling fixed, a reversal of electric poling changed the crack initiation angle from θ₀=+80° to θ₀=−80° using the line extending ahead of the crack as the reference. This effect is also sensitive to the distance from the crack tip. Displayed and discussed are results for r/a=10⁻⁴ and 10⁻¹. Because the theory of magnetoelectroelasticity used in the analysis is based on the assumption of equilibrium where the influence of material microstructure is homogenized, the local space and temporal effects must be interpreted accordingly. Among them are the maximum values of (dW/dV)_max and (dW/dV)_min which refer to as possible sites of yielding and fracture. Since time and size are homogenized, it is implicitly understood that there is more time for yielding as compared to fracture being a more sudden process. This renders a higher dW/dV in contrast to that for fracture. Put it differently, a lower dW/dV with a shorter time for release could be more detrimental.     

Fracture mechanics for a mode III crack in a magnetoelectroelastic solid

In this paper, we developed a Stroh-type formalism for anti-plane deformation and then investigated the fracture mechanics for an elliptical cavity in a magnetoelectroelastic solid under remotely uniform in-plane electromagnetic and/or anti-plane mechanical loading, which allowed us to take the electromagnetic field inside the cavity into account. Reducing the cavity into a crack, we had explicit solutions in closed forms for a mode III crack, which included the extreme cases for an impermeable crack and a permeable crack. The results were illustrated with plots, showing that in the absence of mechanical loads, an applied electric or magnetic field, positive or negative, always tended to close the crack. On the other hand, in the presence of a mechanical load, a negative electric or magnetic field retarded crack growth, while a positive field could either enhance or retard crack propagation, depending on the strengths of the applied electric/magnetic fields and the level of the mechanical load as well. In other words, the effect of electric/magnetic fields on the fracture behavior is mechanical load-dependent.     

Stress intensity factor range and propagation mode of surface cracks under rolling–sliding contact

Finite element analyses were conducted in order to evaluate the mode I and mode II stress intensity factors for inclined edge cracks under cyclic contact load under rolling and rolling–sliding condition. The SIF range depends on crack orientation, crack length to Hertzian contact zone half-width ratio, friction between the crack faces and friction on the contact surface. The results were combined in two compact functions that determine the ΔK_I and ΔK_II values. The crack propagation mode and direction were investigated using both the maximum stress criterion and the minimum strain energy density criterion. The results are displayed in graph form, which allows a fast evaluation of the crack growth condition.     

ANTI-PLANE ANALYSIS FOR ELLIPTICAL INCLUSION IN MAGNETOELECTROELASTIC MATERIALS

This paper considers the multi-field coupling in magneroelectroelastic composite materials consisting of the inclusion and the matrix are magnetoelectroelastic materials. The mechanical,electric and magnetic fields around an elliptical cylinder inclusion are formulated by complex potentials. Inside the inclusion,the strain,electric and magnetic fields are found to be uniform and vary with the shape of the ellipse. When the inclusion is reduced to a crack,along the interface,the strain,electric field strength and magnetic field strength equal the corresponding remote ones,which can be used as the boundary condition. Special cases,such as a rigid and permeable inclusion,a soft and impermeable inclusion,a line inclusion and a crack problem are discussed in detail.     

The Stress State of a Ferromagnetic with an Elliptic Crack in a Homogeneous Magnetic Field

The problem on the stress–strain state of an infinite isotropic body made of a magnetically soft material and containing an elliptic crack is considered. It is assumed that the body is under an external magnetic field perpendicular to the crack plane. The basic characteristics of the stress–strain state and the magnetic field induced are determined and their singularities near the elliptic crack are studied. Formulas are given for the stress intensity factors for the force and magnetic fields near the crack tip     

Angle cracks in two bonded functionally graded piezoelectric material under anti-plane shear

Consider two bonded functionally graded piezoelectric material (FGPM) with finite height. Each material contains an arbitrary oriented crack. The material properties are assumed in exponential forms in the direction normal to the interface. The crack surface condition is assumed to be electrically impermeable or permeable. Using the Fourier transform technique, the problem can be reduced to a system of singular integral equations, which are then solved numerically by applying the Gauss-Chebyshev integration formula to obtain the stress intensity factors at the crack tips. Numerical calculations are carried out to obtain the energy density factor S and the energy release rate G. In impermeable case, the energy release rate has been shown to be negative as the electric loads are applied. The positive definite characteristic of the energy density factor makes it possible for predicting the fracture behavior of the cracked structure. The influences of the non-homogeneous parameters and crack orientation on the energy density factors at the crack tips are discussed in detail. The results show that the energy density factor at the crack tip will be increased when the crack tip is located within the softer material.     

Dynamic fracture analysis of a penny-shaped crack in a magnetoelectroelastic layer

This paper analyzes the dynamic magnetoelectroelastic behavior induced by a penny-shaped crack in a magnetoelectroelastic layer subjected to prescribed stress or prescribed displacement at the layer surfaces. Two kinds of crack surface conditions, i.e., magnetoelectrically impermeable and permeable cracks, are adopted. The Laplace and Hankel transform techniques are employed to reduce the problem to Fredholm integral equations. Field intensity factors are obtained and discussed. Numerical results of the crack opening displacement (COD) intensity factors are presented and the effects of magnetoelectromechanical loadings, crack surface conditions and crack configuration on crack propagation and growth are examined. The results indicate that among others, the fracture behaviors of magnetoelectroelastic materials are affected by the sizes and directions of the prescribed magnetic and/or electric fields, and the effects are strongly dependent on the elastic boundary conditions.     

Wave band gaps in two-dimensional piezoelectric/piezomagnetic phononic crystals

In this paper, the elastic wave propagation in phononic crystals with piezoelectric and piezomagnetic inclusions is investigated taking the magneto-electro-elastic coupling into account. The electric and magnetic fields are approximated as quasi-static. The band structures of three kinds of piezoelectric/piezomagnetic phononic crystals—CoFe₂O₄/quartz, BaTiO₃/CoFe₂O₄ and BaTiO₃–CoFe₂O₄/polymer periodic composites are calculated using the plane-wave expansion method. The piezoelectric and piezomagnetic effects on the band structures are analyzed. The numerical results show that in CoFe₂O₄/quartz structures, only one narrow band gap exists along the Γ–X direction for the coupling of xy-mode and z-mode for the filling fraction f being 0.4; while in BaTiO₃/CoFe₂O₄ composites, only one narrow band gap exists along the Γ–X direction forxy-mode and no band gap exists for z-mode as the filling friction f is 0.5. Moreover, for the new type of magneto-electro-elastic phononic crystal—BaTiO₃–CoFe₂O₄/polymer periodic composite, the band gap characteristics are more superior in the whole considered frequency regions due to the big contrast of the material properties in the two constituents and the effects of the piezoelectricity and piezomagneticity on the band gap structures are remarkable.