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.     

混合边界多裂隙体在简谐变温场作用下的断裂问题（张开型）

给出了无限平面作用有简谐变化的点热源时的位移场、应力场基本解、用间接法构造出混合边值多裂隙体在简谐变温场作用下的热断裂问题的边界积分方程,并离散求解.数值结果表明,该方法求解多裂隙体的简谐热断裂问题精度好,计算工作量少.文中计算了含边界裂缝的平板、含三条平行裂缝的平板在简谐变温场作用下缝端应力强度因子的变化过程,并与实验结果进行了比较,两者吻合良好.     

BIEM analysis of dynamically loaded anti-plane cracks in graded piezoelectric finite solids

Anti-plane cracks in finite functionally graded piezoelectric solids under time-harmonic loading are studied via a non-hypersingular traction based boundary integral equation method (BIEM). The formulation allows for a quadratic variation of the material properties in two directions. The boundary integral equation (BIE) system is treated by using the frequency dependent fundamental solution based on Radon transforms. Its numerical solution provides the displacements and tractions on the external boundary as well as the crack opening displacements from which the mechanical stress intensity factor (SIF) and the electrical displacement intensity factor (EDIF) are determined. Several examples for single and multiple straight and curved cracks demonstrate the applicability of the method and show the influence of the different system parameters.     

Interaction of microcracks with a macrocrack yielded in a narrow strip

Plastic zone growth of collinear cracks has had a longstanding interest in ductile fracture. This work further considers yield zone growth in an isotropic, homogeneous elastic–perfectly plastic infinite plate containing a macrocrack with several neighboring microcracks. Normal loading is considered at distances far away from the cracks. The strip yield is adopted where the plastic zone is assumed to be confined to two narrow strips extending from the ends of a finite length crack while the microcracks are assumed to be elastic. The plastic zone length and crack opening displacement are found from asymptotic solution and compared with finite element solution.     

Effective Diffusivity of Porous Materials with Microcracks: Self-Similar Mean-Field Homogenization and Pixel Finite Element Simulations

We investigate the influence of distributed microcracks on the overall diffusion properties of a porous material using the self-similar cascade continuum micromechanics model within the framework of mean-field homogenization and computational homogenization of diffusion simulations using a high-resolution pixel finite element method. In addition to isotropic, also anisotropic crack distributions are considered. The comparison of the results from the cascade continuum micromechanics model and the numerical simulations provides a deeper insight into the qualitative transport characteristics such as the influence of the crack density on the complexity and connectivity of crack networks. The analysis shows that the effective diffusivity for a disordered microcrack distribution is independent of the absolute length scale of the cracks. It is observed that the overall effective diffusivity of a microcracked material with the microcracks oriented in the direction of transport is not necessarily higher than that of a material with a random orientation of microcracks, independent of the microcrack density.     

Thermal fracture of macrocrack with closure as influenced by microcracks

The fracture stability of macrocracks under uniform heat flux is analyzed to include the effect of a system of microcracks. The interaction of cracks leads to full or partial closure of crack surfaces. The boundary problem is stated and a solution is obtained. The domains where microcracks are closed and/or affect partial closing of the macrocrack are found. Evaluated is the macrocrack tip stress intensity factor accounting for closure.     

Interactive effect of multiple crack growth on fatigue

Interaction of multiple cracks are found in mud when drying and crazing, thin film varnishes or coatings of aeroengine turbine blades. A two-dimensional multiple crack interaction model is developed to simulate the growth of interacting parallel surface cracks. Density and the initial distribution of the microcracks are accounted for in analyzing the growth of a crack to a pre-determined length. Analytical predictions are discussed with reference to experimental observations of fatigue cracks on coated turbine blades. Introduction of a large density of similar cracks can enhance the fatigue life of structural components.     

局部有限元模型边界条件的直接传递方法

提出一种将整体分析得到的节点力或节点位移直接传递到精细化局部有限元模型的方法,即部分混合单元法。沿精细化局部有限元模型周边建立一组过渡单元,该组过渡单元采用与整体模型一致的单元类型和模拟方式,其外侧边界上的节点与整体模型节点的相对坐标对应,内侧边界与精细化局部有限元模型采用基于面约束的方式连接。在外侧边界上根据节点坐标对应施加整体分析获得的节点力或节点位移,过渡单元就可直接将边界条件传递到精细化局部有限元模型。通过贵州红水河特大桥钢-混结合段的精细化有限元分析,验证了本文方法的实用性和有效性。     

The generalised weight function method for crack problems with mixed boundary conditions

The generalised weight function method for determination of stress intensity factors for crack problems with mixed boundary conditions is presented. The method is based on the application of Betti's reciprocal theorem to equivalent crack problems with special boundary conditions, which are converted from the original problems by the principle of superposition. Expressions of stress intensity factors for center cracks subjected to arbitrary stresses in finite plates with various boundary conditions are derived. Examples of practical interest are given. The results reveal the important roles of the boundary displacement constraint and finite dimensions in the crack parameter evaluation.     

Transient dynamic responses of a cracked solid subjected to in-plane loadings

Transient dynamic responses of an elastic cracked solid subjected to in-plane surface loadings are investigated in this study. Two vertical cracks, a surface-breaking crack and a sub-surface crack, are considered. The frequency responses of the plane strain problem are calculated by the computational mechanics combining the finite element method with the boundary integral equation. The finite element method is used for the near-field enclosing the crack, while the boundary integral equation is applied for the far-field to satisfy the Sommerfeld radiation condition. The transient responses are then obtained using fast Fourier transform. Surface displacements, crack opening displacements, and dynamic stress intensity factors are presented to show the significant effects of the cracks. The interaction between the elastic waves and the cracks as well as the mode conversion phenomena can be observed and understood through the numerical simulations.     

主裂纹与微裂纹相互作用问题的有效数值计算方法

提出了平面弹性介质中主裂纹与微裂纹相互作用问题的有效数值计算 方法. 通过把适于单一裂纹的Bueckner原理扩充到含有多裂纹的一般体系，将原问题分解 为承受远处载荷不含裂纹的均匀问题，和在远处不承受载荷但在裂纹面上承受面力的多裂纹 问题. 于是，以应力强度因子作为参量的问题可以通过考虑后者(多裂纹问题)来解决，而 利用提出的杂交位移不连续法，这种多裂纹问题是容易数值求解的. 列举 Cai和 Faber为评价主裂纹与微裂纹相互作用问题的近似方法而列举的算例，说明 该数值方法对分析平面弹性介质中主裂纹与微裂纹相互作用问题既简单又非常有效.     

基于混凝土裂纹数字图像的有限元网格生成

通过结合数字图像处理技术、边界矢量化技术与有限元网格自动生成技术,提出了基于数字图像的混凝土构件裂纹有限元分析方法。本文以混凝土裂纹为研究对象,先采用Canny算子提取裂纹边缘,对中断的边缘和伪边缘利用顺序边缘连接算法进行连接和剔除;然后,通过矢量化技术将裂纹边界拟合成小线段;最后在矢量化边界的基础上,利用自动生成技术生成裂纹结构的有限元网格,实现对裂纹区域受力情况的分析。实验结果表明:本文提出的方法可实现混凝土结构裂纹的图像监测与力学分析的一体化,能够大幅提高现役结构的评估速度和效率。     

纳米晶体铝在单轴拉伸下的位错堆积对微裂纹扩展的影响

摘要：针对纳米晶体材料,研究了单轴拉伸载荷作用下纳米晶体铝中的裂纹与裂纹尖端发射的位错所形成的滑移面之间的相互作用。通过分布位错法,将裂纹和滑移面等效为均匀分布的连续位错,获得了裂纹面上应力场。并引入裂纹尖端的无位错区,研究了裂纹尖端无位错区对微裂纹的萌生和主裂扩展之间的影响。结果表明,不考虑裂纹尖端无位错区时,裂纹长度较短,会先在晶界处形成微裂纹,主裂纹较长时,主裂纹会直接穿晶扩展。滑移面与裂纹尖端夹角较大时,会增加裂纹尖端发射的位错个数,从而抑制主裂纹的扩展。考虑裂纹尖端无位错区时,无位错区先于晶界处出现微裂纹,通过主裂纹与微裂纹之间位错的相互发射,导致裂纹与尖端处微裂纹汇合,有效加速了主裂纹的扩展。     

Elastic interaction of interfacial spherical-cap cracks in hollow particle filled composites

This work analyzes the elastic interaction between two spherical-cap cracks present along the outer surface of a hollow particle embedded in a dissimilar medium under remote uniaxial tensile loading. A semi-analytical approach based on an enriched Galerkin method is adopted to determine stress and deformation fields as functions of particle wall thickness and cracks’ configuration. The present analysis is limited to multiple interfacial spherical-cap cracks; that is, crack propagation is restrained to the particle-matrix interface and possibility of crack kinking in the matrix is not considered. Interfacial crack growth characteristics, conditions for stable crack propagation, equal crack growth, and shielding are established through energy release rate analysis. The study is relevant to the analysis of tensile and flexural failure of syntactic foams used in marine and aerospace applications. Results specialized to glass-vinyl ester syntactic foams demonstrate that particle wall thickness can be used to control crack stability and growth characteristics as well as tailoring the magnitude of the shielding phenomenon. Predictions are compared to finite element findings for validation and to results for penny-shaped cracks to elucidate the role of crack curvature.     

Méthode couplée CVFE/BE pour les calculs acoustiques en domaine non borné

A new method based on the flux conservation for acoustic computations in unbounded domain is presented. The domain is divided in outer and inner regions by a fictitious boundary. The inner problem is solved by a control volume finite element method while the outer flow is treated with a boundary element method. The coupling of the two sub-domains is performed using an appropriate expression for the outgoing fluxes satisfying the Sommerfeld condition. A two-dimensional configuration with a circular interface is considered to show that the method is effective. Results of this method are compared to an analytical solution and to results from the literature.     

任意多孔多裂纹板的裂纹闭合接触分析

基于摩擦接触问题的数学规划解法，采用各向异性体平面弹性理论中的复势方法，建立了含多椭圆孔及裂纹群有限大各向异性板，在任意载荷作用下裂纹闭合或局部闭合问题的有效分析方法。通过在可能闭合的裂纹边界引入互补变量函数并将其展成Fourier级数形式，以Faber级数为工具，应用保角映射技术和最小二乘边界配点法，导出无卸载情况下裂纹面摩擦接触的线性互补模型，并通过算例验证了方法的有效性。数值结果表明，由于采用级数解描述板应力场和位移场，该方法具有较高的计算精度和效率，便于研究裂纹闭合对应力强度因子等断裂参数的影响。     

任意多孔多裂纹有限大板的应力强度因子分析

采用各向异性体平面弹性理论中的复势方法,以Faber级数为工具,应用保角映射技术和最小二乘边界配点法,导出内边界条件精确满足,外边界条件近似满足的含多椭圆孔及裂纹群有限大板在任意载荷作用下的应力场、位移场的级数解,建立了任意多椭圆孔及裂纹群有限大板应力强度因子的有效分析方法,讨论了各参数对裂尖应力强度因子及孔边应力集中的影响．数值结果表明,该方法具有计算精度高、收敛速度快、方便快捷等优点,有利于全面系统地研究各参数对结构断裂性能的影响．     

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.     

Nonlinear frequency mixing of Lamb wave for detecting randomly distributed microcracks in thin plates

In engineering practice, the microcracks usually generate during the initial period of fatigue, and randomly distribute in the metallic structure. It is great challenge for linear non-destructive testing methods to detect those microscopic damages. In this work, the numerical simulations on nonlinear interaction between Lamb wave and microcracks are performed to investigate the behavior of frequency mixing in aluminum plates. From the simulation results, the nonlinear parameter based upon sum frequency shows a monotonically increase with the length of microcracks, and a monotonically decrease with the width of microcracks. The distribution area and density of the cracks also affect the nonlinear interactions between Lamb wave and microcracks. This work theoretically demonstrates that the nonlinear frequency mixing method is a robust strategy to detect microcracks at the initial fatigue stage in metal plates.     

Interaction of three interfacial Griffith cracks between bonded dissimilar orthotropic half planes

The paper deals with the interaction of a pair of outer cracks on a central crack situated at the interface of two dissimilar orthotropic half-planes. The mixed boundary value problem is reduced to solving a pair of simultaneous singular integral equations which have finally been solved numerically by using Jacobi polynomials. The analytical expressions for stress intensity factors at the central crack tip and the expression of the strain energy release rate have been derived for general loading. Numerical values of the interaction effects of the outer cracks on the central crack have been calculated through stress magnification factors. It is seen that the interaction effects are either shielding or amplification depending on the size of the outer cracks and their spacing from the central crack.