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.     

Constant moving crack in a magnetoelectroelastic material under anti-plane shear loading

Analytical solutions for an anti-plane Griffith moving crack inside an infinite magnetoelectroelastic medium under the conditions of permeable crack faces are formulated using integral transform method. The far-field anti-plane mechanical shear and in-plane electrical and magnetic loadings are applied to the magnetoelectroelastic material. Expressions for stresses, electric displacements and magnetic inductions in the vicinity of the crack tip are derived. Field intensity factors for magnetoelectroelastic material are obtained. The stresses, electric displacements and magnetic inductions at the crack tip show inverse square root singularities. The moving speed of the crack have influence on the dynamic electric displacement intensity factor (DEDIF) and the dynamic magnetic induction intensity factor (DMIIF), while the dynamic stress intensity factor (DSIF) does not depend on the velocity of the moving crack. When the crack is moving at very lower or very higher speeds, the crack will propagate along its original plane; while in the range of M_c1 < M < M_c2, the propagation of the crack possibly brings about the branch phenomena in magnetoelectroelastic media.     

Fracture of magnetoelectroelastic composite materials using boundary element method (BEM)

The behavior of cracked linear magnetoelectroelastic solids is analysed by means of the dual Boundary Element Method (BEM) approach. Media possessing fully coupled piezoelectric, piezomagnetic and magnetoelectric effects are considered. An explicit 2-D Green’s function in terms of the extended Stroh formalism for magnetoelectroelastic full-plane under static loading is implemented. Hypersingular integrals arising in the traction boundary integral equations are computed through a regularization technique. Evaluation of fracture parameters directly from computed nodal values is discussed. The stress intensity factors (SIF), the electric displacement intensity factor (EDIF), the magnetic induction intensity factor (MIIF) as well as the mechanical strain energy release rate (MSERR) are evaluated for different crack configurations in both finite and infinite solids subjected to in-plane combined magnetic–electric–mechanical loading conditions. The accuracy of the boundary element solution is confirmed by comparison with selected analytical solutions in the literature. The new results that can be of interest in the design and maintenance of novel magnetoelectroelastic devices are also discussed.     

Analysis of cracked magnetoelectroelastic composites under time-harmonic loading

This paper presents a numerical model for the analysis of cracked magnetoelectroelastic materials subjected to in-plane mechanical, electric and magnetic dynamic time-harmonic loading. A traction boundary integral equation formulation is applied to solve the problem in combination with recently obtained time-harmonic Green’s functions (Rojas-Diaz et al., 2008). The hypersingular boundary integral equations appearing in the formulation are first regularized via a simple change of variables that permits to isolate the singularities. Relevant fracture parameters, namely stress intensity factors, electric displacement intensity factor and magnetic induction intensity factor are directly evaluated as functions of the computed nodal opening displacements and the electric and magnetic potentials jumps across the crack faces. The method is checked by comparing numerical results against existing solutions for piezoelectric solids. Finally, numerical results for scattering of plane waves in a magnetoelectroelastic material by different crack configurations are presented for the first time. The obtained results are analyzed to evaluate the dependence of the fracture parameters on the coupled magnetoelectromechanical load, the crack geometry and the characteristics of the incident wave motion.     

Torsional impact response of a penny-shaped crack in a functional graded strip

The torsional impact response of a penny-shaped crack in a nonhomogeneous strip is considered. The shear modulus is assumed to be functionally graded such that the mathematics is tractable. Laplace and Hankel transforms were used to reduce the problem to solving a Fredholm integral equation. The crack tip stress field is obtained by considering the asymptotic behavior of Bessel function. Explicit expressions of both the dynamic stress intensity factor and the energy density factor were derived. And it is shown that, as crack driving force, they are equivalent for the present crack problem. Investigated are the effects of material nonhomogeneity and strip‘s highness on the dynamic fracture behavior.Numerical results reveal that the peak of the dynamic stress intensity factor can be suppressed by increasing the nonhomogeneity parameter of the shear modulus, and that the dynamic behavior varies little with the adjusting of the strip‘ s highness.     

Stress state of a transversely isotropic magnetoelectroelastic body with a plane crack under antisymmetric loads

The static equilibrium of a transversely isotropic magnetoelectroelastic body with a plane crack of arbitrary shape in the isotropy plane under antisymmetric mechanical loading is studied. The relationships between the stress intensity factors (SIFs) for an infinite magnetoelectroelastic body and the SIFs for a purely elastic body with the same crack and under the same antisymmetric loading are established. This enables the SIFs for a magnetoelectroelastic body to be found directly from the analogous problem of elasticity. As an example of using this result, the SIFs for penny-shaped, elliptic, and parabolic cracks in a magnetoelectroelastic body under antisymmetric mechanical loading are found Translated from Prikladnaya Mekhanika, Vol. 44, No. 10, pp. 37–51, October 2008.     

A contour integral method for dynamic stress intensity factors

The contour integral method previously used to determine static stress intensity factors is applied to dynamic crack problems. The required derivatives of the traction in the reference problem are obtained numerically by the displacement discontinuity method. Stress intensity factors are determined by an integral around a contour which contains a crack tip. If the contour is chosen as the outer boundary of the body, the stress intensity factor is obtained from the boundary values of traction and displacement. The advantage of this path-independent integral is that it yields directly both the opening-mode and sliding-mode stress intensity factors for a straight crack. For dynamic problems, Laplace transforms are used and the dynamic stress intensity factors in the time domain are determined by Durbin's inversion method. An indirect boundary element method, incorporating both displacement discontinuity and fictitious load techniques, is used to determine the boundary or contour values of traction and displacement numerically.     

Transient response of a cracked magnetoelectroelastic strip under anti-plane impact

The transient anti-plane problem of a magnetoelectroelastic strip containing a crack vertical to the boundary is considered. Singular integral equations for the impermeable crack are obtained by using Fourier and Laplace transforms. Numerical results show the effects of the relative loading parameters κ_D and κ_B, and the crack configuration on the dynamic fracture behavior. The results obtained indicate that for the impermeable crack, the electric and magnetic impacts have significant influences on the dynamic stress intensity factor and the dynamic energy density factor.     

Pre-kinking of a moving crack in a magnetoelectroelastic material under in-plane loading

A constant moving crack in a magnetoelectroelastic material under in-plane mechanical, electric and magnetic loading is studied for impermeable crack surface boundary conditions. Fourier transform is employed to reduce the mixed boundary value problem of the crack to dual integral equations, which are solved exactly. Steady-state asymptotic fields near the crack tip are obtained in closed form and the corresponding field intensity factors are expressed explicitly. The crack speed influences the singular field distribution around the crack tip and the effects of electric and magnetic loading on the crack tip fields are discussed. The crack kinking phenomena is investigated using the maximum hoop stress intensity factor criterion. The magnitude of the maximum hoop stress intensity factor tends to increase as the crack speed increases.     

STUDY ON DYNAMIC STRESS INTENSITY FACTORS OF DISK WITH A RADIAL EDGE CRACK SUBJECTED TO EXTERNAL IMPULSIVE PRESSURE

A dynamic weight function method is presented for dynamic stress intensity factors of circular disk with a radial edge crack under external impulsive pressure. The dynamic stresses in a circular disk are solved under abrupt step external pressure using the eigenfunction method. The solution consists of a quasi-static solution satisfying inhomogeneous boundary conditions and a dynamic solution satisfying homogeneous boundary conditions. By making use of Fourier-Bessel series expansion, the history and distribution of dynamic stresses in the circular disk are derived. Furthermore, the equation for stress intensity factors under uniform pressure is used as the reference case, the weight function equation for the circular disk containing an edge crack is worked out, and the dynamic stress intensity factor equation for the circular disk containing a radial edge crack can be given. The results indicate that the stress intensity factors under sudden step external pressure vary periodically with time, and the ratio of the maximum value of dynamic stress intensity factors to the corresponding static value is about 2.0.     

Experimental method of dynamic caustics and its application in fracture mechanics

The equation of static and dynamic caustics, and the formulae determining the position of crack tip and stress intensity factor are given. It is proven that for the case of low speed of crack propagation the static formula is applicable in calculation. A simple method to measure the static stress-optical constants is proposed. An Optical system which is suitable for the experiments of dynamic caustics was set-up and used to study the fracture in beam and rings with initial crack under impact loading. A series of dynamic caustics' photographs and curves showing the variations of corresponding crack lengths and dynamic stress intensity factors with time, are presented.     

The dynamic behaviour of slender structures with cross-sectional cracks

A dynamic model for beams with cross-sectional cracks is discussed. It is shown that a crack can be represented by a consistent, static flexibility matrix. Two different methods for the determination of the flexibility matrix are discussed. If the static stress intensity factors are known, the flexibility matrix can be determined from an integration of these stress intensity factors. Alternatively, static finite element calculations can be used for the determination of the flexibility matrix. Both methods are demonstrated in the present paper. The mathematical model was applied to an edge-cracked cantilevered beam and the eigenfrequencies were determined for different crack lengths and crack positions. These results were compared to experimentally obtained eigenfrequencies. In the experiments, the cracks were modelled by sawing cuts. The theoretical results were, for all crack lengths, in excellent agreement with the experimental data. The dynamic stress intensity factor for a longitudinally vibrating, centrally cracked bar was determined as well. The results compared very well with dynamic finite element calculations. The crack closure effect was experimentally investigated for an edge-cracked beam with a fatigue crack. It was found that the eigenfrequencies decreased, as functions of crack length, at a much slower rate than in the case of an open crack.     

动态焦散线实验方法及其在断裂力学中的初步应用

本文给出了静、动态的焦散线方程,和确定裂纹尖端位置、动应力强度因子K_1~d的计算公式,并证明了在裂纹扩张速度较低时,可用静态公式计算,还研究了带预裂纹的三点弯曲梁和圆环在冲击载荷下的断裂问题,得到了系列动态焦散线照片、裂尖位置和动应力强度因子随时间的变化曲线。     

Transient response of the crack-tip field in a magnetoelectroelastic half-space with a functionally graded coating under impacts

This paper has twofold aims. One is to study the dynamic response of a magnetoelectroelastic half-space with functionally graded coating containing crack at the interface when subjected to sudden impacts. Two different loading positions, where the material and crack surfaces are loaded respectively, are considered. By using the integral transform method, the problem is reduced to solving singular integral equations. Obtained numerical results show that the overshoots of dynamic fracture parameters are strongly amplified or reduced depending on negative or positive gradient, respectively for the case of the material surface being loaded suddenly. This implies that a functionally graded coating with a positive gradient index is preferable in engineering design due to its capability of preventing the structure from cracking. The second objective is to give a comparison of relevant dynamic parameters such as the intensity factors of stress and strain, energy release rate, and energy density factor, and their features are elucidated under dynamic combined loadings. It indicates that the strain intensity factor can overcome the drawbacks of the rest parameters, and may be chosen as an effective fracture parameter, while three others cannot be adopted as fracture criteria to describe the feature of onset of crack growth.     

功能梯度材料裂纹尖端动态应力场

研究受反平面剪切作用的功能梯度材料动态裂纹问题,通过积分变换-对偶积分方程方法推出了裂纹尖端动态应力场,时间域内的动态应力强度因子由Laplace数值反演获得,研究结果表明功能梯度材料的梯度越大,相应的裂纹问题的动态应力强度因子值越低。     

A mode III crack crossing the magnetoelectroelastic bimaterial interface under concentrated magnetoelectromechanical loads

A mode III crack cutting perpendicularly across the interface between two dissimilar semi-infinite magnetoelectroelastic solid is studied under the combined loads of a line force, a line electric charge and a line magnetic charge at an arbitrary location. The impermeable conditions are implied on the crack faces. The technique developed in literature for the elastic bimaterial with a crack cutting interface is exploited to treat the magnetoelectroelastic bimaterial. The Riemann-Hilbert problem can be formulated and solved based on complex variable method. Analytical solutions can be obtained for the entire plane. The intensity factors around crack tips can be defined for the elastic, electric and magnetic fields. It shows that, no matter where the load position is, the electric displacement intensity factors (EDIFs), as well as the magnetic induction intensity factors (MIIFs), are identical in magnitude but opposite in sign for both crack tips, on condition that a line force is solely applied. Alternatively, if only a line electric charge is considered, then the stress intensity factors (SIFs) and the MIIFs exhibit the behavior. Likewise, if only a line magnetic charge is applied, it turns to the SIFs and the EDIFs instead. In addition, the dependence of the intensity factors is graphically shown with respect to the location of a line force. It is found that the SIF for a crack tip tends to be infinite if the applied force is approaching the tip itself, but the EDIF, with the complete opposite trend, tends to be vanishing. Finally, focusing on the more practical case of piezoelectric/piezomagnetic bimaterial, variation of the SIF along with the moduli as well as the piezo constitutive coefficients is explored. These analyses may provide some guidance for material selection by minimizing the SIF. It is also believed that the results obtained in this paper can serve as the Green’s function for the dissimilar magnetoelectroelastic semi-infinite bimaterial with a crack cutting the interface under general magnetoelectromechanical loads.     

Dynamic analysis of a cracked magnetoelectroelastic medium under antiplane mechanical and inplane electric and magnetic impacts

The transient analysis of a magnetoelectroelastic medium containing a crack is made under antiplane mechanical and inplane electric and magnetic impacts. The crack is assumed to penetrate through the solid along the poling direction. By using the Fourier and Laplace transforms, the associated mixed boundary value problem is reduced to a Fredholm integral equation of the second kind, which is solved numerically. By means of a numerical inversion of the Laplace transform, dynamic field intensity factors are obtained in the time domain. Numerical results are presented graphically to show the effects of the material properties and applied electric and magnetic impacts on the dynamic intensity factors of COD and stress, and dynamic energy density factors. The results indicate that except for the intensity factors of electric displacement and magnetic induction, other field intensity factors exhibit apparent transient feature. Moreover, they depend strongly on mechanical input as well as electric and magnetic impacts.     

爆炸应力波作用下动、静裂纹相互作用的实验研究

采用数字激光动态焦散线测试系统,研究爆炸应力波作用下动裂纹与预制静裂纹(水平夹角为90°、150°)相互作用机理,以及裂纹扩展的动态行为。结果表明:(1)在动、静裂纹贯通之前,静裂纹两端便出现焦散斑,动、静裂纹贯通以后,静裂纹沿爆炸应力波传播方向扩展,并且扩展速度小于动裂纹扩展速度,也小于无静裂纹时动裂纹扩展速度; (2)静裂纹存在时,动裂纹扩展的总体长度减小。动裂纹起裂时间缩短,扩展速度基本不受静裂纹的影响,裂纹应力强度因子值大于静裂纹两端值; (3)随着静裂纹水平夹角的增大,动、静裂纹贯通时动裂纹沿水平方向偏转距离增大,静裂纹B端反向扩展与动裂纹相互“咬合”,C端裂纹扩展位移和速度增大。     

An electrically impermeable and magnetically permeable interface crack with a contact zone in magnetoelectroelastic bimaterials under a thermal flux and magnetoelectromechanical loads

An interface crack with a frictionless contact zone at the right crack-tip between two dissimilar magnetoelectroelastic materials under the action of a thermal flux and remote magnetoelectromechanical loads is considered. The open part of the crack is assumed to be electrically impermeable and magnetically permeable, and the crack faces are assumed to be heat insulted. The inhomogeneous combined Dirichlet–Riemann and Hilbert boundary value problems are, respectively, formulated and solved analytically. Stress, electrical displacement intensity factors as well as energy release rate are found in analytical forms, and analytical expressions for the contact zone length have been obtained for both the general case and the case of small contact zone length. Some numerical results are presented, which show clearly the effects of thermal and magnetoelectromechanical loads on the contact zone length, stress intensity factor and energy release rate. Results presented in this paper should have potential applications to the design of multilayered magnetoelectroelastic structures and devices.     

三点弯曲试样动态应力强度因子计算研究

利用Hopkinson压杆对三点弯曲试样进行冲击加载,采集了垂直裂纹面距裂尖2mm和与裂纹面成60°距裂尖5mm处的应变信号。根据裂尖附近测试的应变信号计算试样的动态应力强度因子,并与有限元计算结果进行比较,结果表明由于裂尖有一段疲劳裂纹区,通过裂尖附近应变信号来计算动态应力强度因子时,如果裂尖位置确定不准及粘贴应变片位置不够准确对计算结果将带来很大影响。因此利用应变片法计算动态应力强度因子时,为了获得更准确的计算结果,在实验后应对试件裂纹面进行分析测量,重新确定裂尖位置,必要时需对应变片至裂尖距离进行修正后再计算动态应力强度因子值。