The interaction of doubly periodic cracks

In this paper, an extremely accurate and efficient method for computing the interaction of a set of or multiple sets of general doubly periodic cracks has been presented on the basis of superposition principle, pseudo-traction method, and isolating analysis technique. A great number of typical examples are given in this paper. The stress intensity factors (SIF), the minimum strain energy density factors (SED) of crack tip and the critical stress (CRS) of crack growth are calculated with the accuracy of six significant digits for the rectangularly distributed periodic cracks and five significant digits for the general doubly periodic cracks. The relation of the interaction effect of the double periodic cracks with the periods and the ratio of crack length to crack spacing is analyzed. Also in this paper, the key technique problems for this method are discussed.     

Interactions of Penny-shaped Cracks in Three-dimensional Solids

The interaction of arbitrarily distributed penny-shaped cracks in three-dimensional solids is analyzed in this paper. Using oblate spheroidal coordinates and displacement functions, an analytic method is developed in which the opening and the sliding displacements on each crack surface are taken as the basic unknown functions. The basic unknown functions can be expanded in series of Legendre polynomials with unknown coefficients. Based on superposition technique, a set of governing equations for the unknown coefficients are formulated from the traction free conditions on each crack surface. The boundary collocation procedure and the average method for crack-surface tractions are used for solving the governing equations. The solution can be obtained for quite closely located cracks. Numerical examples are given for several crack problems. By comparing the present results with other existing results, one can conclude that the present method provides a direct and efficient approach to deal with three-dimensional solids containing multiple cracks.The English text was polished by Keren Wang     

Multiple cracks in pressurized hollow cylinder

Presented is a particular solution of the hollow cylinder with one crack; it consists of two parts. The first corresponds to a pair of equal and opposite normal and tangential concentrated forces acting on a crack in an infinite plane region and the second to distributed tractions on both crack surfaces such that the sum of the first and second parts satisfies the prescribed traction boundary conditions on the cracks and cylinder surfaces. The particular solution can be expressed in terms of a density function for each crack giving rise to a system of Fredholm integral equations for the multiple crack system. Several numerical examples will be provided to illustrate the method of solution.     

基于修正权函数的多裂纹无网格模拟

本文采用了一种基于不连续场修正权函数的无网格方法来处理二维平面多裂纹问题。相较于传统的无网格断裂不连续场和奇异场模拟方法,修正权函数法算法简便易实现。采用修正权函数处理多裂纹时,只需要对每一段裂纹周围节点的权函数进行修正,就能同时模拟多裂纹不连续位移场和多裂尖奇异场。本文采用基于不连续场修正权函数的无单元Galerkin方法（EFGM）,对Y型裂纹板、十字型裂纹板和孔边双裂纹板进行了分析。数值结果表明,在不引入扩展基函数情况下,通过修正权函数法能够得到精度较高的应力强度因子解,能较好地拟合多裂纹的裂尖奇异场。     

INTERACTION OF THREE PARALLEL CRACKS IN RECTANGULAR PLATE UNDER CYCLIC LOADS

This paper presents a numerical approach for modeling the interaction between multiple cracks in a rectangular plate under cyclic loads. It involves the formulation of fatigue growth of multiple crack tips under ruixed-mode loading and an extension of a hybrid displacement discontinuity method （a boundary element method） to fatigue crack growth analyses. Because of an intrinsic feature of the boundary element method, a general growth problem of multiple cracks can be solved in a single-region formulation. In the numerical simulation, remeshing of existing boundaries is not necessary for each increment of crack extension. Crack extension is conveniently modeled by adding new boundary elements on the incremental crack extension to the previous crack boundaries. As an example, the numerical approach is used to analyze the fatigue growth of three parallel cracks in a rectangular plate. The numerical results illustrate the validation of the numerical approach and can reveal the effect of the geometry of the cracked plate on the fatigue growth.     

Basic solutions of multiple parallel symmetric mode-III cracks in functionally graded piezoelectric/piezomagnetic material plane

The Schmidt method is adopted to investigate the fracture problem of multiple parallel symmetric and permeable finite length mode-III cracks in a functionally graded piezoelectric/piezomagnetic material plane. This problem is formulated into dual integral equations, in which the unknown variables are the displacement jumps across the crack surfaces. In order to obtain the dual integral equations, the displacement jumps across the crack surfaces are directly expanded as a series of Jacobi polynomials. The results show that the stress, the electric displacement, and the magnetic flux intensity factors of cracks depend on the crack length, the functionally graded parameter, and the distance among the multiple parallel cracks. The crack shielding effect is also obviously presented in a functionally graded piezoelectric/piezomagnetic material plane with mul- tiple parallel symmetric mode-III cracks.     

Automated numerical simulation of the propagation of multiple cracks in a finite plane using the distributed dislocation method

In this paper, an automated numerical simulation of the propagation of multiple cracks in a finite elastic plane by the distributed dislocation method is developed. Firstly, a solution to the problem of a two-dimensional finite elastic plane containing multiple straight cracks and kinked cracks is presented. A serial of distributed dislocations in an infinite plane are used to model all the cracks and the boundary of the finite plane. The mixed-mode stress intensity factors of all the cracks can be calculated by solving a system of singular integral equations with the Gauss–Chebyshev quadrature method. Based on the solution, the propagation of multiple cracks is modeled according to the maximum circumferential stress criterion and Paris' law. Several numerical examples are presented to show the accuracy and efficiency of this method for the simulation of multiple cracks in a 2D finite plane.     

THE CONSTITUTIVE RELATION OF CRACK-WEAKENED ROCK MASSES UNDER AXIAL-DIMENSIONAL UNLOADING

An accurate and efficient numerical method for solving the crack-crack interaction problem is presented. The method is mainly by means of the dislocation model, stress superposition principle and Chebyshev polynomial expansion of the pseudo-traction. This method can be applied to compute the stress intensity factors of multiple kinked cracks and multiple rows of periodic cracks as well as the overall strains of rock masses containing multiple kinked cracks under complex loads. Many complex computational examples are given. The dependence of the crack-crack interaction on the crack configuration, the geometrical and physical parameters, and loads pattern, is investigated. By comparison with numerical results under confining pressure unloading, it is shown that the crack-crack interaction under axial-dimensional unloading is weaker than those under confining pressure unloading. Numerical results for single faults and crossed faults show that the single faults are more unstable than the crossed faults. It is found from numerical results for different crack lengths and different crack spacing that the interaction among kinked cracks decreases with an increase in length of the kinked cracks and the crack spacing under axial-dimensional unloading.     

Extended meshless method based on partition of unity for solving multiple crack problems

An extended meshless method based on partition of unity was used in this study to simulate multiple cracks. The cracks are implicitly denoted by a jump in the displacement field function, which has nodes that have domains of influence completely segmented by cracks. Nodes whose domains of influence are partially segmented by cracks are extended by the crack tip singularity function. The influence domain of a node is independent of cracks so that the sparsity of the system equations should not be affected by cracks and the computing time should not increase with the effect of the cracks. Additionally, r ^?1/2 singularity can be accurately reproduced at the crack tip. Compared with the modified intrinsic enriched meshless method, our method has a higher computational efficiency and precision. Several numerical examples show that the extended meshless method based on partition of unity is feasible and effective in simulating multiple cracks.     

近置多裂纹相互作用的渐近分析方法

考虑到近置裂纹的强相互作用,提出了一种多裂纹相互作用的渐近分析方法. 经典Kachanov 方法将裂纹表面伪面力分解为两部分：均匀分布部分和非均匀部分,并假设裂纹的相互作用 仅由均匀部分引起,而忽略非均匀部分的影响. 该假设大大简化了分析过程,而且当裂纹间 距不是很小时,有很好的精度. 但当裂纹非常靠近或者沿主荷载方向重叠时,由于裂纹尖端 进入了其它裂纹的应力强化区或者应力屏蔽区,强相互作用使得该假设不再合理,从而带来 较大的误差. 为了提高分析近置多裂纹问题的精度,将裂纹表面伪面力分解为抛物线型分布 部分及高阶部分,考虑抛物线型分布张力对其他裂纹的影响,同样忽略高阶部分的影响. 通 过对三共线裂纹及两平行偏置裂纹两个实例的分析,验证了对于近置裂纹,新渐近方法具有 良好的精度.     

双周期裂纹压电材料反平面剪切问题断裂力学分析

研究压电材料双周期裂纹反平面剪切与平面电场作用的问题．运用复变函数方法，获得了该问题严格的闭合解，并由此给出了裂纹尖端应力强度因子和电位移强度因子的精确公式．数值算例显示了裂纹分布特征对材料断裂行为的重要影响．叠间小裂纹能够对主裂纹的应力和电位移场起着屏蔽作用，相反行间小裂纹却起着放大作用，至于钻石形分布裂纹的影响规律则更为复杂．对于某些特殊情形给予了解答并导出一系列有意义的结果。     

基于均匀荷载面曲率的简支裂纹梁损伤识别

为了采用模态参数对结构裂纹进行定位与定量,基于集中柔度模型,采用无质量的扭转弹簧模拟裂纹,建立简支裂纹梁的振动微分方程。针对现有柔度曲率指标仅能判断裂纹的大致范围,基于线性插值理论,建立裂纹位置与相邻测点均匀荷载面曲率差的关系,提出裂纹进一步定位公式,实现裂纹位置的精确定位。针对现有大多数损伤识别方法无法实现裂纹的损伤定量,基于位移曲率与结构刚度和弯矩的关系,理论推导了均匀荷载面曲率的结构刚度损伤程度识别方法,基于弹簧串联原理和线刚度思想,首次提出串联等效线刚度模型,建立裂纹深度与均匀荷载面曲率的关系,实现裂纹深度的定量。通过简支裂纹梁数值算例,考虑多裂纹的损伤情况,验证了新方法对裂纹定位与定量的有效性。     

基于局域分析的疲劳短裂纹群体演化随机模型

采用局域裂纹数密度描述金属材料中不同局部区域的疲劳短裂纹群体损伤的发展情况通过考虑在不同局域存在的材料性质的随机涨落及局部损伤对损伤总量发展的影响，建立了局域裂纹数密度演化随机方程对方程数值求解从而模拟了材料的疲劳短裂纹损伤过程结果显示出主裂纹出现的随机性，并讨论了裂纹总数与最大裂纹尺度在统计意义上的演化特征     

Environmentally assisted initiation and growth of multiple surface cracks

The initial stages of stress corrosion on an amorphous polymer is investigated. This is done by exposing stressed specimens of polycarbonate to an acetone and water solution. The surface develops two distinct features of degradation that appear on different length scales when subjected to tensile stress. Small pits form on the surface and make it rough. These pits are in the order of micrometers, and are found to be randomly distributed. They occur even without load and seem to slightly increase in number with increasing stress. In the millimeter domain, visible to the bare eye, surface cracks are formed transverse to the direction of loading. The occurrence of cracks is seen to have a positive stress threshold value, exceeding which, a linear increase of number of cracks with stress is found. The manners in which the cracks grow and coalesce on the surface are examined. It is seen that they do not meet crack tip to crack tip. Instead, they avoid each other initially and coalesce crack tip to crack side. The results are discussed in the light of mechanical considerations. A stress analysis for a few configurations of meeting cracks supports the experimental observations. With assumptions of stress corrosion crack growth and coalescence, a simulation of cracks growing from randomly distributed initiation sites is performed. Similar crack patterns as obtained in the experiments are found.     

DYNAMIC BEHAVIOR OF SEVERAL CRACKS IN FUNCTIONALLY GRADED STRIP SUBJECTED TO ANTI-PLANE TIME-HARMONIC CONCENTRATED LOADS

This paper contains a theoretical formulations and solutions of multiple cracks sub- jected to an anti-plane time-harmonic point load in a functionally graded strip. The distributed dislocation technique is used to construct integral equations for a functionally graded material strip weakened by several cracks under anti-plane time-harmonic load. These equations are of Cauchy singular type at the location of dislocation, which are solved numerically to obtain the dislocation density on the faces of the cracks. The dislocation densities are employed to evaluate the stress intensity factor and strain energy density factors （SEDFs） for multiple cracks with differ- ent configurations. Numerical calculations are presented to show the effects of material properties and the crack configuration on the dynamic stress intensity factors and SEDFs of the functionally graded strip with multiple curved cracks.     

Thermoelastic analysis of multiple defects with the extended finite element method

In this paper, the extended finite element method (XFEM) is adopted to analyze the interaction between a sin-gle macroscopic inclusion and a single macroscopic crack as well as that between multiple macroscopic or micro-scopic defects under thermal/mechanical load. The effects of different shapes of multiple inclusions on the material thermomechanical response are investigated, and the level set method is coupled with XFEM to analyze the interaction of multiple defects. Further, the discretized extended finite element approximations in relation to thermoelastic prob-lems of multiple defects under displacement or temperature field are given. Also, the interfaces of cracks or materials are represented by level set functions, which allow the mesh assignment not to conform to crack or material interfaces. Moreover, stress intensity factors of cracks are obtained by the interaction integral method or the M-integral method, and the stress/strain/stiffness fields are simulated in the case of multiple cracks or multiple inclusions. Finally, some numer-ical examples are provided to demonstrate the accuracy of our proposed method.     

A coupled IEM/FEM approach for solving elastic problems with multiple cracks

In this paper, the detailed two-dimensional infinite element method (IEM) formulation with infinite element (IE)–finite element (FE) coupling scheme for investigating mode I stress intensity factor in elastic problems with imbedded geometric singularities (e.g. cracks) is presented. The IE–FE coupling algorithm is also successfully extended to solve multiple crack problems. In this method, the domain of the primary problem is subdivided into two sub-domains modeled separately using the IEM for the multiple crack sub-domain, and the FEM for the uncracked sub-domain. In the IE sub-domain, the similarity partition concept together with the IEM formulation are employed to automatically generate a large number of infinitesimal elements, layer by layer, around the tip of each crack. All degrees of freedom related to the IE sub-domain, except for those associated with the coupling interface, are condensed and transformed to form a finite master IE for each crack with master nodes on sub-domain boundary only. All of the stiffness matrices constructed in the IE sub-domains are assembled into the system stiffness matrix for the FE sub-domain. The resultant FE solution with a symmetrical stiffness matrix, having the singularity effect of imbedded cracks in IEs, is required only for solving multiple crack problems.Using these efficient numerical techniques a very fine mesh pattern can be established around each crack tip without increasing the degree of freedom of the global FEM solution. One is easily allowed to conduct parametric analyses for various crack sizes without changing the FE mesh. Numerical examples are presented to show the performance of the proposed method and compared with the corresponding known results where available.     

Numerical computation and analytical estimation of stress-intensity factors for strips with multiple edge cracks

In this paper, stress-intensity factors for a two-dimensional problem are determined. Strips with multiple symmetrical edge cracks in tension are investigated. A simple analytical estimation is compared to numerical results. The influence of penetration of the crack faces and mixed-mode loading on the numerical results is investigated. A simple method to estimate stress-intensity factors for strips with multiple edge cracks is proposed.     

Evolving crack patterns in thin films with the extended finite element method

This paper develops the extended finite element method (XFEM) to evolve patterns of multiple cracks, in a brittle thin film bonded to an elastic substrate, with a relatively coarse mesh, and without remeshing during evolution. A shear lag model describes the deformation in three dimensions with approximate field equations in two-dimensions. The film is susceptible to subcritical cracking, obeying a kinetic law that relates the velocity of each crack to its energy release rate. At a given time, the XFEM solves the field equations and calculates the energy release rate of every crack. For a small time step, each crack is extended in the direction of maximal hoop stress, and by a length set by the kinetic law. To confirm the accuracy of the XFEM, we compare our simulation to the exiting solutions for several simple crack patterns, such as a single crack and a set of parallel cracks. We then simulate the evolution of multiple cracks, initially in a small region of the film but of different lengths, showing curved crack propagation and crack tip shielding. Starting with multiple small cracks throughout the film, the XFEM can generate the well-known mud crack pattern.     

弹性半平面中边界裂纹研究的新方法

讨论了载荷作用在裂纹面上的弹性半平面边界裂纹问题.研究以线弹性断裂力学为基础,采用复变函数方法以及Riemann-Hilbert(R-H)边值问题的一般理论,将问题分拆为含有限裂纹的全平面问题与无裂纹的半平面问题的叠加,计算得到裂纹尖端的应力强度因子.与文献结果比较,该方法具有精度高的优点.