Stress singularities for crack normal to and terminating at bimaterial interface on orthotropic half-plates

The problem of a crack normal to and terminating at an interface in two joined orthotropic plates is considered and the eigenequation for the asymptotic behavior of stresses at the crack tip on the interface is given in an explicit form. It is found that the singular stress field around the crack tip can be separated into two independent fields, respectively of the mode I and II. Also it is found that for both the mode I and II deformations the effects of elastic constants on the stress singularity order can be respectively expressed by three material parameters, two of which are the same for both the mode I and mode II deformations.     

宏微观双尺度运动裂纹模型面内拉伸下的解析解*

研究裂纹动态扩展中宏微观因素相互作用机制与微观裂尖区的钝化效应。平面拉伸状态下,宏观主裂纹以恒定速度运动。通过一个介观约束应力过渡区,将宏观主裂纹与微观裂尖区相连接,由此建立了一个宏微观双尺度运动裂纹模型。应用弹性动力学与复变函数理论,分别在宏观与微观尺度下对该模型进行解析求解,获得了解析解。通过裂纹张开位移从宏观到微观的连续性条件与宏微观应力场协调条件,将两个不同尺度下的解相耦合,获得了计算宏微观损伤区特征长度的显式表达式。研究表明,运动裂纹的宏观应力场仍具有通常的r&;#61485;1/2的奇异性。由于微观裂尖的钝化,微观应力场奇异性的阶次有所降低,与宏观应力场相比具有弱奇异性。双尺度运动裂纹模型中,可允许裂纹运动速度达到剪切波速,解除了经典运动裂纹理论中裂纹速度不能超过Rayleigh波速的限制。数值结果表明,介观损伤过渡区与裂尖微观损伤区尺寸,及裂纹张开位移等,与裂纹运动速度、材料性质、约束应力比、裂尖钝化角度等因素有关。     

On the stress intensity factors associated with cracks interacting with an interface between two elastic media

Closed form expressions are obtained for the stresses at a crack tip when a crack is approaching a welded boundary (or a free surface) and when it has just passed through the interface. The solutions which are obtained in terms of a small parameter, the distance from or through the interface, are given in explicit form for the mode 3 situation and for some mode 1 and 2 cases. The importance of the change of stress singularity when the crack meets the interface is demonstrated.     

Linear elastic crack tip modelling by the displacement discontinuity method

A special crack tip element is developed to model the feature of a stress singularity at a crack tip, based on an elastic solution for an arbitrary displacement discontinuity. By incorporating this crack tip element into the conventional displacement discontinity method, the propagation of a crack with an arbitrary configuration can be modelled using the Griffith-Irwin energy criterion in a linear elastic medium. The accuracy and consistency of the method is illustrated by various examples; the results are seen to be in excellent agreement with the exact analytical solutions and the known experimental evidence. The study shows that the displacement discontinuity method is efficient and has a potential for modelling the behaviour of a complex crack.     

压电材料中渗透性周期共线反平面裂纹问题

本文利用复变函数方法,借助于Riemann-Schwarz延拓技术和保形映照方法,研究了渗透性边界条件下周期共线反平面裂纹问题,获得了解的表达式,得到了力学和电学强度因子。结果表明在裂纹尖端应力和电位移的奇异性都与远场作用的应力载荷和裂纹长度有关,其中应力的奇异性与材料无关,电位移的奇异性则与材料有关,电载荷对裂尖的奇异性没有影响。最后,运用数值算例,给出周期裂纹间的干涉效应和裂纹的尺度效应。     

Micro-crack initiation at tip of a rigid line inhomogeneity in piezoelectric materials

For the square-root singularity shear stress found at the tip of a rigid line inhomogeneity (an anti-crack) in piezoelectric media, one possible way of releasing high strain energy is to initiate a micro-crack at the inhomogeneity tip. In our current study, a dislocation pileup model for micro-crack initiation at the inhomogeneity tip is proposed based on Zener-Stroh crack initiation mechanism. An interesting and important physical result that emerges from the analysis is that the critical stress intensity factor for the anti-crack (the line inhomogeneity) can be related to the fracture toughness of a conventional Griffith crack in the same material. Analytical results further show that under mechanical loading, the critical stress and electric displacement intensity factors of an anti-crack are only related to the corresponding intensity factors of stress and electric displacement of the crack, respectively. While if the anti-crack is under displacement loading (with net dislocation pile-up at the inhomogeneity tip), the critical stress and electric displacement intensity factors of an anti-crack depend on both of the total mechanical dislocations b_T and electricity dislocations b_D.     

Investigation of the interaction of two collinear cracks in anisotropic elasticity materials by means of the nonlocal theory

In this paper, the interaction of two collinear cracks in anisotropic elasticity materials subjected to an anti-plane shear loading is investigated by means of the nonlocal theory. By use of the Fourier transform, the problem can be solved with the help of a pair of triple integral equations, in which the unknown variable is the displacement on the crack surface. To solve the triple integral equations, the displacement on the crack surface is expanded in a series of Jacobi polynomials. Unlike the classical elasticity solutions, it is found that no stress singularity is present at the crack tip. The nonlocal elastic solutions yield a finite hoop stress at the crack tip, thus allowing us to use the maximum stress hypothesis as a fracture criterion. The magnitude of the finite stress field depends on the crack length, the distance between two cracks and the lattice parameter of materials.     

Orthogonal meshless finite volume method applied to elastodynamic crack problems

An orthogonal meshless finite volume method has been presented to solve some elastodynamic crack problems. An orthogonal weighted basis function is used to construct shape function so there is no problem of singularity in this new form. In this work, for three-dimensional dynamic fracture problems, a new displacement function is used at the tip of the crack to give a new OMFVM. When the new OMFVM is used, the singularity of the stresses at the tip of the crack can be shown to be better than that in the primal OMFVM. High computational efficiency and precision are other benefits of the method. Solving some sample crack problems of thin-walled structures show a good performance of this method.     

Analytical solution for the singular stress distribution due to V-notch in an orthotropic plate material

On the basis of general solutions of two-dimensional linear elasticity, displacement and singular stress fields near the singular point in orthotropic materials are derived in closed form expressions. According to the presented expressions, analysis formulas of displacement and singular stress fields near the tip of a V-notch under the symmetric and the anti-symmetric modes are obtained subsequently. The open literatures devoted to developing stress singularity near the tip of the V-notch in anisotropic or orthotropic materials. In this study, however, not only direct eigenequations were derived, but also the explicit solutions of displacement and singular stress fields were obtained. At the end, the correctness of the formulas of the singular stress field near the tip of the V-notch has been verified by FEM analysis.     

Strip yield analysis of two collinear unequal cracks in an infinite sheet

Plastic zone sizes and crack tip opening displacements (CTODs) are obtained for two collinear cracks in an infinite sheet subjected to known remote stress. Analysis is conducted by assuming the crack accompanying plastic zones as a fictitious crack and formulating integral equations based upon traction free and no stress singularity conditions. In addition, critical remote stress, plastic zone sizes, and CTODs when the adjacent plastic zones touched are obtained by assuming the coalesced fictitious cracks as a single fictitious crack and formulating integral equations based upon no stress singularity and zero coalesced point displacement conditions. Extensive numerical results are presented.     

The scattering of harmonic elastic anti-plane shear waves by a Griffith crack in a piezoelectric material plane by using the non-local theory

In this paper, the dynamic behavior of a Griffith crack in a piezoelectric material plane under anti-plane shear waves is investigated by using the non-local theory for impermeable crack face conditions. For overcoming the mathematical difficulties, a one-dimensional non-local kernel is used instead of a two-dimensional one for the anti-plane dynamic problem to obtain the stress and the electric displacement near the crack tips. By using the Fourier transform, the problem can be solved with the help of two pairs of dual integral equations. These equations are solved using the Schmidt method. Contrary to the classical elasticity solution, it is found that no stress and electric displacement singularity is present near the crack tip. The non-local dynamic elastic solutions yield a finite hoop stress near the crack tip, thus allowing for a fracture criterion based on the maximum dynamic stress hypothesis. The finite hoop stress at the crack tip depends on the crack length, the circular frequency of incident wave and the lattice parameter. For comparison results between the non-local theory and the local theory for this problem, the same problem in the piezoelectric materials is also solved by using local theory.     

Partition of unity enrichment for bimaterial interface cracks

Partition of unity enrichment techniques are developed for bimaterial interface cracks. A discontinuous function and the two‐dimensional near‐tip asymptotic displacement functions are added to the finite element approximation using the framework of partition of unity. This enables the domain to be modelled by finite elements without explicitly meshing the crack surfaces. The crack‐tip enrichment functions are chosen as those that span the asymptotic displacement fields for an interfacial crack. The concept of partition of unity facilitates the incorporation of the oscillatory nature of the singularity within a conforming finite element approximation. The mixed‐mode (complex) stress intensity factors for bimaterial interfacial cracks are numerically evaluated using the domain form of the interaction integral. Good agreement between the numerical results and the reference solutions for benchmark interfacial crack problems is realized. Copyright © 2004 John Wiley & Sons, Ltd.     

Evaluation of crack tip stress and strain fields under nonproportional loading in a viscoelastic material

The crack tip strain and stress fields in a viscoelastic material under nonproportional loading conditions are evaluated. In order to evaluate the strain field, the crack tip displacement field is measured by using the digital image correlation (DIC) technique. This displacement field is then approximated by using the theoretically obtained crack tip displacement field in viscoelastic materials. The result shows that the approximation method can smoothly reconstruct the experimentally obtained displacement field. From the approximated displacement field, the crack tip strain field can be precisely obtained by using the differential form of the theoretical displacement. On the other hand, the crack tip stress field is analyzed by using the stress function. This suggests that the strain and stress fields can be independently evaluated. In addition, different time dependencies between stress and strain fields near the crack tip are observed. Based on this experiment, we can discuss the several criteria for the crack propagation directions in viscoelastic materials.     

A finite element stress analysis of a crack in a bi-material plate

A stress analysis is presented for the problem of a crack in one material of a bi-material plate located perpendicular to the material interface. A numerical solution using the finite element techniques to determine the force displacement relationships is used. Knowing this, a work integral method is used to determine the stress intensity factors for the crack. Since the work integral is independent of path, the path of integration can be chosen far enough away from the crack tip to avoid the complications of the crack tip singularity. The problem is studied for a number of cases where the crack length to plate width ratio, distance from crack tip to material interface, and the ratio of material constants were varied as parameters.     

Displacements, strains and stresses associated with propagating cracks in materials with continuously varying properties

Summary. Crack tip stress, strain and displacement fields for a propagating crack along the direction of property gradation in functionally graded materials (FGMs) are obtained through an asymptotic analysis coupled with a displacement potential approach. The analysis for the opening mode is performed assuming two types of property variation: (a) linearly varying shear modulus with constant density, and (b) exponentially varying shear modulus and density. The first six terms in the series expansion of the stress, strain and displacement fields for the dynamic crack are derived to explicitly bring out the influence of nonhomogeneity on the structure of the displacement, strain and stress fields. The analysis revealed that crack tip stress fields retain the inverse square root singularity and only the higher order terms in the expansion are influenced by material inhomogeneity. Using these stress, strain and displacement fields, contours of constant maximum shear stress, constant first stress invariant and constant in-plane displacements are generated and the effect of the nonhomogeneity parameter on these contours is discussed.     

Investigation of anti-plane shear behavior of a Griffith crack in a piezoelectric material by using the non-local theory

In this paper, the behavior of a Griffith crack in a piezoelectric material under anti-plane shear loading is investigated by using the non-local theory for impermeable crack surface conditions. By using the Fourier transform, the problem can be solved with two pairs of dual integral equations. These equations are solved using Schmidt method. Numerical examples are provided. Contrary to the previous results, it is found that no stress and electric displacement singularity is present at the crack tip.     

FEM for evaluation of weight functions for SIF, COD and higher-order coefficients with application to a typical wedge splitting specimen

In the evaluation of accurate weight functions for the coefficients of first few terms of the linear elastic crack tip fields and the crack opening displacement (COD) using the finite element method (FEM), singularities at the crack tip and the loading point need to be properly considered. The crack tip singularity can be well captured by a hybrid crack element (HCE), which directly predicts accurate coefficients of first few terms of the linear elastic crack tip fields. A penalty function technique is introduced to handle the point load. With the use of these methods numerical results of a typical wedge splitting (WS) specimen subjected to wedge forces at arbitrary locations on the crack faces are obtained. With the help of appropriate interpolation techniques, these results can be used as weight functions. The range of validity of the so-called Paris equation, which is widely used in the evaluation of the COD from the stress intensity factors (SIFs), is established.     

The collapsed cubic isoparametric element as a ingular element for crack probblems

For the 12?node quadrilateral isoparametric elements, it is shown that the inverse square root singularity of the strain field at the crack tip can be obtained by the simple technique of collapsing the quadrilateral elements into triangular elements around the crack tip and placing the two side nodes of each side of the triangles at 1/9 and 4/9 of the length of the side from the tip. This is analogous to placing the mid?side nodes at quarter points in the vicinity of the crack tip for the quadratic isoparametric elements. The advantages of this method are that the displacement compatibility is satisfied throughout the region and that there is no need of special crack tip elements. The stress intensity factors can be accurately obtained by using general purpose programs having isoparametric elements such as NASTRAN.     

Numerical analysis of plane cracks in strain-gradient elastic materials

The classical linear elastic fracture mechanics is not valid near the crack tip because of the unrealistic singular stress at the tip. The study of the physical nature of the deformation around the crack tip reveals the dominance of long-range atomic interactive forces. Unlike the classical theory which incorporates only short range forces, a higher-order continuum theory which could predict the effect of long range interactions at a macro scale would be appropriate to understand the deformation around the crack tip. A simplified theory of gradient elasticity proposed by Aifantis is one such grade-2 theory. This theory is used in the present work to numerically analyze plane cracks in strain-gradient elastic materials. Towards this end, a 36 DOF C¹ finite element is used to discretize the displacement field. The results show that the crack tip singularity still persists but with a different nature which is physically more reasonable. A smooth closure of the structure of the crack tip is also achieved.