A crack perpendicular to the bimaterial interface in finite solid

The dislocation simulation method is used in this paper to derive the basic equations for a crack perpendicular to the bimaterial interface in a finite solid. The complete solutions to the problem, including the T stress and the stress intensity factors are obtained. The stress field characteristics are investigated in detail. It is found that when the crack is within a weaker material, the stress intensity factor is smaller than that in a homogeneous material and it decreases when the distance between the crack tip and interface decreases. When the crack is within a stiffer material, the stress intensity factor is larger than that in a homogeneous material and it increases when the distance between the crack tip and interface decreases. In both cases, the stress intensity factor will increase when the ratio of the size of a sample to the crack length decreases. A comparison of stress intensity factors between a finite problem and an infinite problem has been given also. The stress distribution ahead of the crack tip, which is near the interface, is shown in details and the T stress effect is considered.     

Permeable cracks between two dissimilar piezoelectric materials

An interfacial crack with electrically permeable surfaces between two dissimilar piezoelectric ceramics under electromechanical loading is investigated. An exact expression for singular stress and electric fields near the tip of a permeable crack between two dissimilar anisotropic piezoelectric media are obtained. The interfacial crack-tip fields are shown to consist of both an inverse square root singularity and a pair of oscillatory singularities. It is found that the singular fields near the permeable interfacial crack tip are uniquely characterized by the real valued stress intensity factors proposed in this paper. The energy release rate is obtained in terms of the stress intensity factors. The exact solution of stress and electric fields for a finite interfacial crack problem is also derived.     

Concentrated impact loading on one face of a semi-infinite crack in an anisotropic elastic material

The response of an unbounded anisotropic elastic body containing a semi-infinite crack subjected to a concentrated impact force on one of the crack faces is studied. An exact solution of the dynamic stress intensity factors is obtained from a linear superposition of the solution of Lamb’s problem and a solution of a dislocation emitting from the crack tip. The stress intensity factors exhibit square-root singularity upon the arrival of the Rayleigh wave at the crack tip. As the Rayleigh wave passes through the crack tip, the stress intensity factors either instantaneously assume the static values or gradually approach to zero. Several numerical examples are given for isotropic, cubic and orthotropic materials.     

Observation of Stress Field Around an Oscillating Crack Tip in a Quenched Thin Glass Plate

The variation of stress field around an oscillating crack tip in a quenched thin glass plate is observed using instantaneous phase-stepping photoelasticity. The successive images around the propagating crack are recorded by a CCD camera that is equipped with a pixelated micro-retarder array. Then, the phase maps of the principal stress difference and the principal direction are easily obtained even though the photoelastic fringes cannot be visualized. The path of the crack growth as well as the stress intensity factors and the crack tip constraint are obtained from these phase distributions. Results show that the mode I stress intensity factor and the crack tip constraint vary remarkably with the crack growth. In addition, the results show that the mode-II stress intensity factor exists even though the crack propagates smoothly.     

The weight function for cracks in piezoelectrics

The weight function in fracture mechanics is the stress intensity factor at the tip of a crack in an elastic material due to a point load at an arbitrary location in the body containing the crack. For a piezoelectric material, this definition is extended to include the effect of point charges and the presence of an electric displacement intensity factor at the tip of the crack. Thus, the weight function permits the calculation of the crack tip intensity factors for an arbitrary distribution of applied loads and imposed electric charges. In this paper, the weight function for calculating the stress and electric displacement intensity factors for cracks in piezoelectric materials is formulated from Maxwell relationships among the energy release rate, the physical displacements and the electric potential as dependent variables and the applied loads and electric charges as independent variables. These Maxwell relationships arise as a result of an electric enthalpy for the body that can be formulated in terms of the applied loads and imposed electric charges. An electric enthalpy for a body containing an electrically impermeable crack can then be stated that accounts for the presence of loads and charges for a problem that has been solved previously plus the loads and charges associated with an unsolved problem for which the stress and electric displacement intensity factors are to be found. Differentiation of the electric enthalpy twice with respect to the applied loads (or imposed charges) and with respect to the crack length gives rise to Maxwell relationships for the derivative of the crack tip energy release rate with respect to the applied loads (or imposed charges) of the unsolved problem equal to the derivative of the physical displacements (or the electric potential) of the solved problem with respect to the crack length. The Irwin relationship for the crack tip energy release rate in terms of the crack tip intensity factors then allows the intensity factors for the unsolved problem to be formulated, thereby giving the desired weight function. The results are used to derive the weight function for an electrically impermeable Griffith crack in an infinite piezoelectric body, thereby giving the stress intensity factors and the electric displacement intensity factor due to a point load and a point charge anywhere in an infinite piezoelectric body. The use of the weight function to compute the electric displacement factor for an electrically permeable crack is then presented. Explicit results based on a previous analysis are given for a Griffith crack in an infinite body of PZT-5H poled orthogonally to the crack surfaces.     

一种新的计算各向异性材料裂纹尖端应力强度因子杂交元法

利用有限元特征分析法研究了平面各向异性材料裂纹端部的奇性应力指数以及应力场和位移场的角分布函数,以此构造了一个新的裂纹尖端单元。文中利用该单元建立了研究裂纹尖端奇性场的杂交应力模型,并结合Hellinger-Reissner变分原理导出应力杂交元方程,建立了求解平面各向异性材料裂纹尖端问题的杂交元计算模型。与四节点单元相结合,由此提出了一种新的求解应力强度因子的杂交元法。最后给出了在平面应力和平面应变下求解裂纹尖端奇性场的算例。算例表明,本文所述方法不仅精度高,而且适应性强。     

压电材料反平面界面裂纹的杂交元分析

基于奇异性电弹场数值特征解开发了一种新型反平面界面裂纹尖端单元。将新型单元与四节点压电P-S单元组装,求解从绝缘到导通的任意电边界条件下,压电结构反平面界面裂纹尖端电弹场的数值解。考察了层厚、载荷类型和裂纹面间电边界条件等对反平面界面裂纹尖端断裂参数的影响。算例证明新型单元能使P-S单元数显著降低,计算结果更为精确。     

Drag load modeling for multiple fractures in a borehole

Fracture mechanics modeling of multiple fractures generated by a tailored pressure pulse in a borehole is a complex, coupled gas dynamics/solid mechanics problem. The specific problem addressed in this paper is the determination of crack tip stress intensity factors for arbitrary drag loadings along the crack surfaces. Fracture networks containing an even number of equally spaced radial cracks of alternating lengths are used as models. Finite element calculations and integration schemes similar to those used for arbitrary normal loads along the crack surfaces are performed. Stress intensity factors arising from the drag loadings tend to increase with increasing crack length ratios. Thus, effects of drag loading may be significant in determining crack tip stress intensity factors for a multiple fracture network generated by a controlled gas pulse.     

INTERACTION OF A DISLOCATION DIPOLE WITH A MODE Ⅲ INTERFACE CRACK

Interaction between a screw dislocation dipole and a mode Ⅲ interface crack is investigated. By using the complex variable method, the closed form solutions for complex potentials are obtained when a screw dislocation dipole lies inside a medium. The stress fields and the stress intensity factors at the tip of the interface crack produced by the screw dislocation dipole are given. The influence of the orientation, the dipole arm and the location of the screw dislocation dipole as well as the material mismatch on the stress intensity factors is discussed. zThe image force and the image torque acting on the screw dislocation dipole center are also calculated. The mechanical equilibrium position of the screw dislocation dipole is examined for various material property combinations and crack geometries. The results indicate that the shielding or anti-shielding effect on the stress intensity factor increases abruptly when the dislocation dipole approaches the tip of the crack. Additionally, the disturbation of the interface crack on the motion of the dislocation dipole is also significant.     

On the higher order asymptotic analysis of a non-uniformly propagating dynamic crack along an arbitrary path

Transient mixed-mode elastodynamic crack growth along arbitrary smoothly varying paths is considered. Asymptotically, the crack tip stress field is square root singular with the angular variation of the singular term depending weakly on the instantaneous values of the crack tip speed and on the mode-I and mode-II stress intensity factors. However, for a material particle at a small distance away from the moving crack tip, the local stress field will depend not only on the instantaneous values of the crack tip speed and stress intensity factors, but also on the past history of these time dependent quantities. In addition, for cracks propagating along curved paths the stress field is also expected to depend on the nature of the curved crack path. Here, a representation of the crack tip fields in the form of an expansion about the crack tip is obtained in powers of radial distance from the tip. The higher order coefficients of this expansion are found to depend on the time derivative of crack tip speed, the time derivatives of the two stress intensity factors as well as on the instantaneous value of the local curvature of the crack path. It is also demonstrated that even if cracks follow a curved path dictated by the criterion K ₁₁ ^d =0, the stress field may still retain higher order asymmetric components related to non-zero local curvature of the crack path.     

INTERACTION OF A DISLOCATION DIPOLE WITH A MODE III INTERFACE CRACK

Interaction between a screw dislocation dipole and a mode III interface crack is investigated. By using the complex variable method, the closed form solutions for complex potentials are obtained when a screw dislocation dipole lies inside a medium. The stress fields and the stress intensity factors at the tip of the interface crack produced by the screw dislocation dipole are given. The influence of the orientation, the dipole arm and the location of the screw dislocation dipole as well as the material mismatch on the stress intensity factors is discussed. The image force and the image torque acting on the screw dislocation dipole center are also calculated. The mechanical equilibrium position of the screw dislocation dipole is examined for various material property combinations and crack geometries. The results indicate that the shielding or anti-shielding effect on the stress intensity factor increases abruptly when the dislocation dipole approaches the tip of the crack. Additionally, the disturbation of the interface crack on the motion of the dislocation dipole is also significant.     

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.     

Effect of a flat inclusion on stress intensity factor of a semi-infinite crack

The elastic plane interaction between an arbitrarily located and oriented flat inclusion and a semi-infinite crack subjected to a remote Mode I loading is considered. The method uses distributions of edge dislocations to formulate integral expressions of flat inclusion (including crack) tractions and is shown to be very accurate by a test problem. The stress intensity factors of the main crack tip are presented for a variety of crack inclusion geometries. It is seen that the flat inclusion could either yield a stress enhancement or stress shielding effect to the main crack tip depending upon the location, orientation and thickness of the flat inclusion, and depending upon the modulus ratios of the flat inclusion to matrix.     

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

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

基于数字散斑相关方法测定Ⅰ型裂纹应力强度因子

提出了一种通过数字散斑相关方法测定金属材料Ⅰ型裂纹尖端位置和应力强度因子的实验方法.实验采用疲劳试验机对含Ⅰ型缺口的Cr12MoV钢试件预制裂纹,通过数字散斑相关方法测试试件在三点弯曲加载条件下裂纹的扩展过程及裂尖区域的位移场.将位移场数据代入裂尖位移场方程组,采用牛顿-拉普森方法求解含未知参量的裂尖非线性位移场方程组,计算裂尖位置和应力强度因子.实验结果表明,采用该方法可以准确地测定金属材料Ⅰ型裂纹应力强度因子、裂尖位置及裂纹扩展长度,解决了以往研究中因不能准确测定裂纹尖端位置,而无法准确计算Ⅰ型裂纹裂尖断裂参数的难题,揭示了金属材料裂纹扩展过程中应力强度因子演化特征.     

Interface crack problems for mode Ⅱ of double dissimilar orthotropic composite materials

The fracture problems near the interface crack tip for mode Ⅱ of double dissimilar orthotropic composite materials are studied. The mechanical models of interface crack for mode Ⅱ are given. By translating the governing equations into the generalized bi-harmonic equations,the stress functions containing two stress singularity exponents are derived with the help of a complex function method. Based on the boundary conditions,a system of non-homogeneous linear equations is found. Two real stress singularity exponents are determined be solving this system under appropriate conditions about himaterial engineering parameters. According to the uniqueness theorem of limit,both the formulae of stress intensity factors and theoretical solutions of stress field near the interface crack tip are derived. When the two orthotropic materials are the same,the stress singularity exponents,stress intensity factors and stresses for mode Ⅱ crack of the orthotropic single material are obtained.     

Photoelastic modeling of the fracture of viscoelastic orthotropic plates with a crack

The well-known equations of photoelasticity of linear viscoelastic bodies are used to describe the photoelastic behavior of a viscoelastic orthotropic plate with a crack. Expressions for the stress intensity factors (SIFs) at the crack tip are obtained using photoelastic measurements. The time dependence of the SIFs is analyzed and shown to be determined by the angles between directions of the crack and tension     

双材料界面裂纹应力强度因子的边界元分析

采用双材料基本解建立边界元法基本方程,计算双材料界面裂纹尖端附近的应用力和位移场。不离散界面,并设置面力奇异四分之一点裂尖单元以提高计算精度。数值结果表明,本文的方法具有较高的精度和效率。     

A technique for studying interaction between a screw dislocation dipole or a concentrated load and a mode III crack crossing an interface

A method of potentially wide application is developed for deriving analytical expressions of the elastic interaction between a screw dislocation dipole or a concentrated force and a crack cutting perpendicularly across the interface of a bimaterial. The cross line composed of the interface and the crack is mapped into a line, and then the complex potentials are educed. The Muskhelishvili method is extended by creating a Plemelj function that matches the singularity of the real crack tips, and eliminates the pseudo tips’ singularity induced by the conformal mapping. The stress field is obtained after solving the Riemann–Hilbert boundary value problem. Based on the stress field expressions, crack tip stress intensity factors, dislocation dipole image forces and image torque are formulated. Numerical curves show that both the translation and rotation must be considered in the static equilibrium of the dipole system. The crack tip stress intensity factor induced by the dipole may rise or drop and the crack may attract or reject the dipole. These trends depend not only on the crack length, but also on the dipole location, the length and the angle of the dipole span. Generally, the horizontal image force exerted at the center of the dislocation dipole is much smaller than the vertical one. Whether the dipole subjected to clockwise torque or anticlockwise torque is determined by whether the Burgers vector of the crack-nearby dislocation of the dipole is positive or negative. A concentrated load induces no singularity to crack tip stress fields as the load is located at the crack line. However, as the concentrated force is not located on the crack line but approaches the crack tip, the nearby crack tip stress intensity factor K_IIIu increases steeply to infinity.     

压电材料裂纹顶端条状电饱和区模型的力学分析

在线性压电本构方程框架下,对裂纹顶端条状电饱和区模型进行了严格的数学分析．完整地考虑了各向异性力电耦合效应．建立了电饱和区尺寸与外加电场的依赖关系．证实了当裂纹垂直极化轴时,压电材料的断裂应力随着外加正电场的增加而减小,随着外加负电场的增加而增加．当裂纹平行于极化轴时,与极化轴平行的外加电场对断裂应力无影响