Anti-plane fracture analysis for the weak-discontinuous interface in a non-homogeneous piezoelectric bi-material structure

The concept of weak discontinuity is extended to functionally graded piezoelectric bi-material interface, and fracture analysis for the weak discontinuous interface is performed by the methods of Fourier integral transform and Cauchy singular integral equation. Numerical results of the total energy release rate (TERR) and the mechanical strain energy release rate (MSERR) are obtained to show the effects of non-homogeneity parameters, geometrical parameters and loads. Parametric studies yield three conclusions: (1) To reduce the weak-discontinuity of the interface is beneficial to resisting interfacial fracture. The effect of the weak-discontinuity of the interface on TERR and MSERR still depends on the strip width. The wider the strip, the more sensitive the TERR and MSERR will be to the weak-discontinuity of the interface. (2) To predict the effect of electric load on crack propagation, MSERR is more appropriate than TERR to be used as a fracture parameter. To predict the effect of mechanical load on crack propagation, both of them could be used as fracture parameters, and MSERR is more conservative. (3) Mechanical load and negative electric displacement load would promote crack propagation, but positive electric displacement load would retard it. For the structure applied by combined mechanical and positive electric displacement loads, crack propagation may be impeded by appropriately selecting the strip width and the ratio of non-homogeneity parameters.     

Perturbation analysis of Mode III interfacial cracks advancing in a dilute heterogeneous material

The paper addresses the problem of a Mode III interfacial crack advancing quasi-statically in a heterogeneous composite material, that is a two-phase material containing elastic inclusions, both soft and stiff, and defects, such as microcracks, rigid line inclusions and voids. It is assumed that the bonding between dissimilar elastic materials is weak so that the interface is a preferential path for the crack. The perturbation analysis is made possible by means of the fundamental solutions (symmetric and skew-symmetric weight functions) derived in Piccolroaz et al. (2009). We derive the dipole matrices of the defects in question and use the corresponding dipole fields to evaluate “effective” tractions along the crack faces and interface to describe the interaction between the main interfacial crack and the defects. For a stable propagation of the crack, the perturbation of the stress intensity factor induced by the defects is then balanced by the elongation of the crack along the interface, thus giving an explicit asymptotic formula for the calculation of the crack advance. The method is general and applicable to interfacial cracks with general distributed loading on the crack faces, taking into account possible asymmetry in the boundary conditions.The analytical results are used to analyse the shielding and amplification effects of various types of defects in different configurations. Numerical computations based on the explicit analytical formulae allows for the analysis of crack propagation and arrest.     

Experimental and numerical analysis of fretting crack formation based on 3D X-FEM frictional contact fatigue crack model

Nowadays, numerical simulation of 3D fatigue crack growth is easily handled using the eXtended Finite Element Method coupled with level set techniques. The finite element mesh does not need to conform to the crack geometry. Most difficulties associated to complex mesh generation around the crack and the re-meshing steps during the possible propagation are hence avoided. A 3D two-scale frictional contact fatigue crack model developed within the X-FEM framework is presented in this article. It allows the use of a refined discretization of the crack interface independent from the underlying finite element mesh and adapted to the frictional contact crack scale. A stabilized three-field weak formulation is also proposed to avoid possible oscillations in the local solution linked to the LBB condition when tangential slip is occurring. Two basic three-dimensional numerical examples are presented. They aim at illustrating the capacities and the high level of accuracy of the proposed X-FEM model. Stress intensity factors are computed along the crack front. Finally an experimental 3D ball/plate fretting fatigue test with running conditions inducing crack nucleation and propagation is modeled. 3D crack shapes defined from actual experimental ones and fretting loading cycle are considered. This latter numerical simulation demonstrates the model ability to deal with challenging actual complex problems and the possibility to achieve tribological fatigue prediction at a design stage based on the fatigue crack modeling.     

A generalized screw dislocation interacting with interfacial cracks along a circular inhomogeneity in magnetoelectroelastic solids

The interaction of a generalized screw dislocation with circular arc interfacial cracks under remote antiplane shear stresses, in-plane electric and magnetic loads in transversely isotropic magnetoelectroelastic solids is dealt with. By using the complex variable method, the general solutions to the problem are presented. The closed-form expressions of complex potentials in both the inhomogeneity and the matrix are derived for a single circular-arc interfacial crack. The intensity factors of stress, electric displacement and magnetic induction are provided explicitly. The image forces acting on the dislocation are also calculated by using the generalized Peach–Koehler formula. For the case of piezoelectric matrix and piezomagnetic inclusion, the shielding and anti-shielding effect of the dislocation upon the stress intensity factors is evaluated in detail. The results indicate that if the distance between the dislocation and the crack tip remains constant, the dislocation in the interface will have a largest shielding effect which retards the crack propagation. In addition, the influence of the interfacial crack geometry and materials magnetoelectroelastic mismatch upon the image force is discussed. Numerical computations show that the perturbation effect of the above parameters upon the image force is significant. The main result shows that a stable or unstable equilibrium point may be found when a screw dislocation approaches the surface of the crack from infinity which differs from the perfect bonded case under the same conditions. The present solutions contain a number of previously known results which can be shown to be special cases.     

Dynamic crack deflection and penetration at interfaces in homogeneous materials: experimental studies and model predictions

We examine the deflection/penetration behavior of dynamic mode-I cracks propagating at various speeds towards inclined weak planes/interfaces of various strengths in otherwise homogeneous isotropic plates. A dynamic wedge-loading mechanism is used to control the incoming crack speeds, and high-speed photography and dynamic photoelasticity are used to observe, in real-time, the failure mode transition mechanism at the interfaces. Simple dynamic fracture mechanics concepts used in conjunction with a postulated energy criterion are applied to examine the crack deflection/penetration behavior and, for the case of interfacial deflection, to predict the crack tip speed of the deflected crack. It is found that if the interfacial angle and strength are such as to trap an incident dynamic mode-I crack within the interface, a failure mode transition occurs. This transition is characterized by a distinct, observable and predicted speed jump as well as a dramatic crack speed increase as the crack transitions from a purely mode-I crack to an unstable mixed-mode interfacial crack.     

Stress intensity factor analysis of a three-dimensional interfacial corner between anisotropic bimaterials under thermal stress

A numerical method using a path-independent H-integral based on the conservation integral was developed to analyze the singular stress field of a three-dimensional interfacial corner between anisotropic bimaterials under thermal stress. In the present method, the shape of the corner front is smooth. According to the theory of linear elasticity, asymptotic stress near the tip of a sharp interfacial corner is generally singular as a result of a mismatch of the materials’ elastic constants. The eigenvalues and the eigenfunctions are obtained using the Williams eigenfunction method, which depends on the anisotropic materials’ properties and the geometry of an interfacial corner. The order of the singularity related to the eigenvalue is real, complex or power-logarithmic. The amplitudes of the singular stress terms can be calculated using the H-integral. The stress and displacement around an interfacial corner for the H-integral are obtained using finite element analysis. In this study, a proposed definition of the stress intensity factors of an interfacial corner, which includes those of an interfacial crack and a homogeneous crack, is used to evaluate the singular stress fields. Asymptotic solutions of stress and displacement around an interfacial corner front are uniquely obtained using these stress intensity factors. To prove the accuracy of the present method, several different kinds of examples are shown such as interfacial corners or cracks in three-dimensional structures.     

Semi-permeable Yoffe-type interfacial crack analysis in MEE composites based on the strip electro-magnetic polarization saturation model

Fracture analysis of a semi-permeable Yoffe-type interfacial crack propagating subsonically in magneto-electro-elastic(MEE) composites is presented based on the strip electromagnetic polarization saturation(SEMPS) model. The electro-magnetic fields inside the crack are considered under the semi-permeable boundary condition. Nonlinear effects near the interfacial crack tip are represented by different electro-magnetic saturation zones. Utilizing the extended Stroh's method, we derive the moving dislocation densities as well as intensity factor and energy release rate for Yoffe-type MEE interfacial crack. Numerical results through an iterative approach are presented to show the characteristics of fracturedominant parameters with respect to propagation velocity and boundary condition category. The fracture-dominant parameters under the semi-permeable boundary condition are lower than those under the impermeable one, which implies that the electro-magnetic fields in the crack gap can retard the propagation of MEE interfacial crack.     

The deformation of the front of a 3D interface crack propagating quasistatically in a medium with random fracture properties

One studies the evolution in time of the deformation of the front of a semi-infinite 3D interface crack propagating quasistatically in an infinite heterogeneous elastic body. The fracture properties are assumed to be lower on the interface than in the materials so that crack propagation is channelled along the interface, and to vary randomly within the crack plane. The work is based on earlier formulae which provide the first-order change of the stress intensity factors along the front of a semi-infinite interface crack arising from some small but otherwise arbitrary in-plane perturbation of this front. The main object of study is the long-time behavior of various statistical measures of the deformation of the crack front. Special attention is paid to the influences of the mismatch of elastic properties, the type of propagation law (fatigue or brittle fracture) and the stable or unstable character of 2D crack propagation (depending on the loading) upon the development of this deformation.     

ANALYSIS OF A SCREW DISLOCATION INSIDE A CIRCULAR INCLUSION WITH INTERFACIAL CRACKS IN PIEZOELECTRIC COMPOSITES

IntroductionPiezoelectric materials have potentials for use in many modern devices and compositestructures. The presence of various defects, such as inclusions, holes, dislocations andcracks, can greatly influence their characteristics and coupled behavio…     

Theoretical analysis of crack front instability in mode I+III

This paper focusses on the theoretical prediction of the widely observed crack front instability in mode I+III, that causes both the crack surface and crack front to deviate from planar and straight shapes, respectively. This problem is addressed within the classical framework of fracture mechanics, where the crack front evolution is governed by conditions of constant energy-release-rate (Griffith criterion) and vanishing stress intensity factor of mode II (principle of local symmetry) along the front. The formulation of the linear stability problem for the evolution of small perturbations of the crack front exploits previous results of Movchan et al. (1998) (suitably extended) and Gao and Rice (1986), which are used to derive expressions for the variations of the stress intensity factors along the front resulting from both in-plane and out-of-plane perturbations. We find exact eigenmode solutions to this problem, which correspond to perturbations of the crack front that are shaped as elliptic helices with their axis coinciding with the unperturbed straight front and an amplitude exponentially growing or decaying along the propagation direction. Exponential growth corresponding to unstable propagation occurs when the ratio of the unperturbed mode III to mode I stress intensity factors exceeds some “threshold” depending on Poisson's ratio. Moreover, the growth rate of helical perturbations is inversely proportional to their wavelength along the front. This growth rate therefore diverges when this wavelength goes to zero, which emphasizes the need for some “regularization” of crack propagation laws at very short scales. This divergence also reveals an interesting similarity between crack front instability in mode I+III and well-known growth front instabilities of interfaces governed by a Laplacian or diffusion field.     

岩石偏心圆孔单裂纹平台圆盘的动态裂纹扩展与止裂

针对平台圆环构型的优点, 提出偏心圆孔单裂纹平台圆盘(cracked eccentrically holed flattened disc, CEHFD), 该试样具有更长的断裂路径。利用霍普金森压杆加载系统, 径向冲击CEHFD试样, 完成Ⅰ型动态断裂实验。砂岩试样表面粘贴应变片和裂纹扩展计, 用于监测裂纹动态起裂、扩展和止裂的全过程。实验表明, 在整个断裂过程中, 裂纹非匀速扩展, 裂纹扩展速度在裂纹起裂后加速上升, 在裂纹止裂前有明显的减速, 与地震时断层的动态破裂全过程完全吻合。采用实验-数值-解析法得到动态应力强度因子, 其时间历程呈现先增大后减小的趋势。根据断裂过程不同时刻, 得到相应的动态起裂韧度、扩展韧度及止裂韧度。在动态断裂全过程中, 动态扩展韧度为速度的函数, 变化趋势与速度一致, 随着时间先增大后减小; 动态起裂韧度大于动态止裂韧度, 止裂韧度随着裂纹最大扩展速度的增大而降低, 并且有较大的离散性。     

层状粘接复合结构界面力学参数声导波定征方法

在推导层状粘接复合结构良好粘接及存在弱界面、滑移界面和脱层等几种不同界面条件下声导波的广义频散方程的基础上,分析了界面径向与轴向力学参数对声导波传播特性的影响,进一步提出以频散特性为基础的超声导波定征方法和在最小二乘意义下的反向算法对粘接复合结构层间界面力学参数进行了估计,分析了影响估计准确性的各种因素,研究了超声导波定征方法对粘接复合结构层间力学参数的灵敏度及其在误差传递中的意义。     

Bridging micro to macroscale fracture properties in highly heterogeneous brittle solids: weak pinning versus fingering

The effect of strong toughness heterogeneities on the macroscopic failure properties of brittle solids is investigated in the context of planar crack propagation. The basic mechanism at play is that the crack is locally slowed down or even trapped when encountering tougher material. The induced front deformation results in a selection of local toughness values that reflect at larger scale on the material resistance. To unravel this complexity and bridge micro to macroscale in failure of strongly heterogeneous media, we propose a homogenization procedure based on the introduction of two complementary macroscopic properties: An apparent toughness defined from the loading required to make the crack propagate and an effective fracture energy defined from the rate of energy released by unit area of crack advance. The relationship between these homogenized properties and the features of the local toughness map is computed using an iterative perturbation method. This approach is applied to a circular crack pinned by a periodic array of obstacles invariant in the radial direction, which gives rise to two distinct propagation regimes: A weak pinning regime where the crack maintains a stationary shape after reaching an equilibrium position and a fingering regime characterized by the continuous growth of localized regions of the fronts while the other parts remain trapped. Our approach successfully bridges micro to macroscopic failure properties in both cases and illustrates how small scale heterogeneities can drastically affect the overall failure response of brittle solids. On a broader perspective, we believe that our approach can be used as a powerful tool for the rational design of heterogeneous brittle solids and interfaces with tailored failure properties.     

热辐射应力图像分析确定裂纹尖端塑性区

本文借助于热辐应力图像分析技术记录了铝质紧养伤拉伸试件的热辐射应力图像。根据裂尖周围温度变化曲线的分布规律确定其特征区的,从而得以裂尖塑性区的开头曲线,并与理论值比较,符合良好。对裂纹扩展面进行ＳＥＭ分析,得出该试件的断裂特性,对于进行细观损伤研究有重要的意义。     

Interaction between rigid-disc inclusion and penny-shaped crack under elastic time-harmonic wave incidence

Influence of a rigid-disc massive inclusion on a neighboring penny-shaped crack induced by the time-harmonic wave propagation in an infinite elastic matrix is investigated by the numerical solution of associated 3D elastodynamic problem. No restrictions on the mutual orientation of interacting objects and direction of wave incidence are assumed. The inclusion is perfectly bonded with a matrix and supposes the translations and rotations, the crack faces are load-free. Frequency-domain problem is reduced to a system of boundary integral equations (BIEs) relative to the interfacial stress jumps (ISJs) on the inclusion and the crack opening displacements (CODs). The subtraction technique in conjunction with mapping technique, under taking into account the structure of solution at the fronts of inclusion and crack, is applied for regularization of BIEs obtained. A discrete analogue of equations is constructed by using the collocation scheme. Numerical calculations are carried out for the grazing incidence of a plane P-wave on the crack, where the interacting inclusion is coplanar and perpendicular to the crack, and has the same radius. The shielding and amplification effects of inclusion are assessed by the analysis of mode-I stress intensity factor (SIF) in the crack vicinity depending on the wave number, incident wave direction, position of the crack front point, inclusion mass, crack-inclusion orientation and distance.     

表面裂纹疲劳扩展的数值模拟(Ⅱ)

根据Paris疲劳裂纹扩展规律,对拉伸和纯弯曲疲劳载荷下表面裂纹扩展进行了数值模拟。数值模型中,用三次样条函数曲线拟合裂纹尖端,在裂纹扩展增量计算中考虑了裂纹闭合影响。裂纹形状演化的模拟结果与Newman和Raju经验公式预测结果进行了比较,表明了所采用的数值模拟方法的实用性。研究发现,裂纹闭合对疲劳裂纹扩展过程中的裂纹形状演化以及裂纹尖端的应力强度因子(SIF)分布都有明显影响。同裂纹形状演化一样,疲劳裂纹扩展过程中裂纹尖端的SIF分布表现出明显的特征。最后,建议了一个简单函数来统一描述表面裂纹尖端的SIF分布。     

Mode-3 spontaneous crack propagation along functionally graded bimaterial interfaces

The effects of combining functionally graded materials (FGMs) of different inhomogeneous property gradients on the mode-3 propagation characteristics of an interfacial crack are numerically investigated. Spontaneous interfacial crack propagation simulations were performed using the newly developed spectral scheme. The numerical scheme derived and implemented in the present work can efficiently simulate planar crack propagation along functionally graded bimaterial interfaces. The material property inhomogeneity was assumed to be in the direction normal to the interface. Various bimaterial combinations were simulated by varying the material property inhomogeneity length scale. Our parametric study showed that the inclusion of a softening type FGM in the bimaterial system leads to a reduction in the fracture resistance indicated by the increase in crack propagation velocity and power absorbed. An opposite trend of increased fracture resistance was predicted when a hardening material was included in the bimaterial system. The cohesive tractions and crack opening displacements were altered due to the material property inhomogeneity, but the stresses ahead of the cohesive zone remained unaffected.     

类岩石材料中表面裂纹扩展模式的实验研究

采用粘贴式轴拉方法实现了类岩石材料的直接拉伸实验,通过在类岩石材料试件中预制表面裂纹研究了类岩石材料中表面裂纹的扩展模式.研究发现,粘贴式轴拉方法可以满足一般类岩石材料的直接拉伸实验;类岩石材料中表面裂纹首先从试件正面开始扩展;表面裂纹的存在极大地影响着材料的破坏模式,类岩石材料在试件的正面上会表现出明显的二维穿透性裂纹的扩展形态,但是受其影响裂纹在厚度方向扩展会发生偏转,从而使得试件背面的最终贯通方向没有与加载方向垂直,而是呈一倾角,并且和表面裂纹的倾角和深度有关.     

Three-dimensional effects in interfacial crack propagation

The paper describes the use of crack-opening interferometry for examining the variation in normal crack-opening displacements (NCOD) along the front of an interfacial crack in an edge-cracked bimaterial strip under biaxial loading. For the glass/epoxy combination considered here, the crack front was concave in the direction of crack growth, in contrast to previous observations with a glass/polyurethane/glass sandwich specimen and cracks in homogeneous materials. The NCOD were greatest in the interior of the specimen for all mode-mixes considered and the exponents in a power-law fit of NCOD versus distance from the crack front decreased towards the free surface. The exponents varied with mode-mix, suggesting that interfacial crack-front geometries could be similarly affected.