三维裂纹扩展轨迹的边界元数值模拟

提出了一种对三维裂纹扩展轨迹进行数值模拟的新方法。采用一种新的具有C^1连续性、高精度的单节点二次边界单元，使边界元(BEM)的分析效率和裂纹张开位移(COD)、应力强度因子(SIF)的精度大大提高。采用裂纹张开位移全场拟合法(GCDFP)求出裂纹面前缘的SIF，所得到的SIF达到与所用的COD资料同样的精度。使用Paris公式求出裂纹前缘各点的裂纹扩展增量，并用三次B样条函数对这些增量进行拟合，得到新的光滑裂纹前缘。根据以上思想方法，开发了具有较高的计算效率和精度的数值模拟软件。此软件可以自动跟踪裂纹扩展，得到裂纹扩展的轨迹。运用该软件对椭圆和矩形裂纹的扩展轨迹进行了数值模拟。其结果与理论上的预言完全一致，裂纹最后都趋于一个圆裂纹，具有实际指导意义。     

A boundary element application for mixed modeloading idealized sawtooth fracture surface

In this paper, a 2-D elastic-plastic BEM formulation predicting the reduced mode IIand the enhanced mode I stress intensity factors are presented. The dilatant boundary conditions (DBC) are assumed to be idealized uniform sawtooth crack surfaces and an effective Coulombsliding law. Three types of crack face boundary conditions, i.e. (1) BEM sawtooth model-elasticcenter crack tip; (2) BEM sawtooth model-plastic center crack tip; and (3) BEM sawtoothmodel-edge crack with asperity wear are presented. The model is developed to attempt todescribe experimentally observed non-monotonic, non-linear dependence of shear crack behavioron applied shear stress, superimposed tensile stress, and crack length. The asperity sliding isgoverned by Coulombs law of friction applied on the inclined asperity surface which hascoefficient of friction μ. The traction and displacement Greens functions which derive fromNaviers equations are obtained as well as the governing boundary integral equations for an infiniteelastic medium. Accuracy test is performed by comparison stress intensity factors of the BEMmodel with analytical solutions of the elastic center crack tip. The numerical results show thepotential application of the BEM model to investigate the effect of mixed mode loading problemswith various boundary conditions and physical interactions.     

On the half-infinite crack problem in thermo-electro-elasticity

This paper concerns itself with the fundamental solutions of the thermo-electro-elastic field in an infinite medium, weakened by a half-infinite plane crack with two identical point thermal loads applied on the crack surfaces. The corresponding mixed boundary value problem is solved by virtue of the potential theory method conjugated with the general solutions. The boundary governing equations are solved by using the results available in literature. Exact and complete three-dimensional (3D) fundamental solutions are presented in terms of elementary functions. The singularity at the crack tip is analyzed explicitly. The obtained solutions will be of high significance to the related BEM analysis.     

The Analysis of Dynamic Stress Intensity Factor for Semi-Circular Surface Crack Using Time-Domain BEM Formulation

ＩｎｔｒｏｄｕｃｔｉｏｎＴｈｅｔｒａｎｓｉｅｎｔａｎａｌｙｓｉｓｆｏｒ 3＿Ｄｅｌａｓｔｏｄｙｎａｍｉｃｃｒａｃｋｐｒｏｂｌｅｍｓｉｓａｎｉｍｐｏｒｔａｎｔｃｏｎｔｅｎｔｉｎｆｒａｃｔｕｒｅｍｅｃｈａｎｉｃｓ.Ｂｅｃａｕｓｅｏｆｔｈｅｄｉｆｆｉｃｕｌｔｉｅｓｏｆｔｈｅｐｒｏｂｌｅｍｓ,ｔｈｅｃａｌｃｕｌａｔｉｏｎｏｆＤＳＩＦｆｏｒ 3＿Ｄｃｒａｃｋｉｓａｐｒｏｂｌｅｍｕｎｓｏｌｖｅｄｃｏｍｐｌｅｔｅｌｙｕｐｔｏｔｈｅｐｒｅｓｅｎｔ.Ｃｈｅｎ[1]ｆｉｒｓｔｌｙｕｓｅｄｔｈｅｆｉｎｉｔｅｄｉｆｆｅｒｅｎｃｅ…     

Application of a new fast multipole BEM for simulation of 2D elastic solid with large number of inclusions

A fast multipole method (FMM) is applied for BEM to reduce both the operation and memory requirement in dealing with very large scale problems. In this paper, a new version of fast multipole BEM for 2D elastostatics is presented and used for simulation of 2D elastic solid with a large number of randomly distributed inclusions combined with a similar subregion approach. Generalized minimum residual method (GMRES) is used as an iterative solver to solve the equation system formed by BEM iteratively. The numerical results show that the scheme presented is applicable to certain large scale problems. The project supported by the National Nature Science Foundation of China (10172053) and the Ministry of Education     

三维热权函数法和多虚拟裂纹扩展技术

热权函数法直接利用温度场与热权函数的乘积的积分来求应力强度因子(SIF)的过渡过程，它可以免除对每一时刻进行热弹性力学有限元或边界元应力分析，计算效率大大提高。本文给出了三维热权函数法的基本方程，并提出了求解三维热权函数法基本方程的多虚拟裂纹扩展法(MVCE法)。在MVCE法中，可以引入无穷多个虚拟裂纹扩展模式；虚拟裂纹扩展模式与应力强度因子的插值直接相联系，所得到的方程组的系数矩阵是一个三对角矩阵，具有良好的计算性能。它对于裂纹前缘SIF分布急剧变化的情况，有良好的数值模拟能力。实例计算表明，MVCE法具有权高的计算效率，并具有很高的计算精度。     

Dynamic interaction between a sub-interface crack and the interface in a bi-material: time-domain BEM analysis

Summary Transient response of a sub-interface crack in a bi-material is studied with emphasis on the dynamic interaction between the crack and the interface, by combining the traditional time-domain displacement boundary element method (BEM) and the non-hypersingular traction BEM. Computations are performed for an unbounded bi-material with a crack subjected to impact tensile loading on its faces or incident impact waves and a bounded rectangular bi-material plate under remote impact tensile loading. Numerical results of the dynamic stress intensity factors (DSIFs) and dynamic interface tractions are presented for various material combinations and crack locations. It is shown that pronounced increases in DSIFs and the interface tractions may be caused in some cases because of the dynamic interaction between the crack and the interface.This work was initialized during the second author's stay at Institute of Mechanics, TU Darmstadt, Germany under the support of the Alexander von Humboldt Foundation. Discussion on the BEM formulation with Dr. Seelig is gratefully acknowledged. The first two authors are also grateful for the partial support by the China National Natural Science Foundation under Grant No. 10025211 and the NJTU Scientific Paper Fund (PD195).     

Energy domain integral applied to solve center and double-edge crack problems in three-dimensions

A three-dimensional Boundary Element Method (BEM) implementation of the energy domain integral for the numerical computation of the crack energy release rate is presented in this paper. The domain expression of the energy release rate is naturally compatible with the BEM, since stresses, strains and derivatives of displacements at internal points can be evaluated using the appropriate boundary integral equations. The pointwise crack energy release rate is evaluated along the three-dimensional crack front over a cylindrical domains that surround a segment of the crack front. The accuracy of the implementation is demonstrated by solving several problems, which include geometries containing straight as well as curved crack fronts.     

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.     

Mode III fracture analysis by Trefftz boundary element method

This paper presents a hybrid Trefftz (HT) boundary element method (BEM) by using two indirect techniques for mode III fracture problems. Two Trefftz complete functions of Laplace equation for normal elements and a special purpose Trefftz function for crack elements are proposed in deriving the Galerkin and the collocation techniques of HT BEM. Then two auxiliary functions are introduced to improve the accuracy of the displacement field near the crack tips, and stress intensity factor (SIF) is evaluated by local crack elements as well. Furthermore, numerical examples are given, including comparisons of the present results with the analytical solution and the other numerical methods, to demonstrate the efficiency for different boundary conditions and to illustrate the convergence influenced by several parameters. It shows that HT BEM by using the Galerkin and the collocation techniques is effective for mode III fracture problems. The project supported by the National Natural Science Foundation of China(10472082). The English text was polished by Keren Wang.     

BEM for simulation of a 2D elastic body with randomly distributed circular inclusions

Based on our 2D BEM software THBEM2 which can be applied to the simulation of an elastic body with randomly distributed identical circular holes, a scheme of BEM for the simulation of elastic bodies with randomly distributed circular inclusions is proposed. The numerical examples given show that the boundary element method is more accurate and more effective than the finite element method for such a problem. The scheme presented can also be successfully used to estimate the effective elastic properties of composite materials. Project supported by the National Natural Science Foundation of China (No. 19772025).     

Finite element analyses of cracks in piezoelectric structures: a survey

Piezoelectric materials have widespread applications in modern technical areas such as mechatronics, smart structures or microsystem technology, where they serve as sensors or actuators. For the assessment of strength and reliability of piezoelectric structures under combined electrical and mechanical loading, the existence of cracklike defects plays an important role. Meanwhile, piezoelectric fracture mechanics has been established quite well, but its application to realistic crack configurations and loading situations in piezoelectric structures requires the use of numerical techniques as finite element methods (FEM) or boundary element methods (BEM). The aim of this paper is to review the state of the art of FEM to compute the coupled electromechanical boundary value problem of cracks in 2D and 3D piezoelectric structures under static and dynamic loading. In order to calculate the relevant fracture parameters very precisely and efficiently, the numerical treatment must account for the singularity of the mechanical and electrical fields at crack tips. The following specialized techniques are presented in detail (1) special singular crack tip elements, (2) determination of intensity factors K _I –K _IV from near tip fields, (3) modified crack closure integral, (4) computation of the electromechanical J-integral, and (5) exploitation of interaction integrals. Special emphasis is devoted to a realistic modeling of the dielectric medium inside the crack, leading to specific electric crack face boundary conditions. The accuracy, efficiency, and applicability of these techniques are examined by various example problems and discussed with respect to their advantages and drawbacks for practical applications.     

构件三维断裂与疲劳力学及其在航空工程中的应用

本文总结评述了断裂力学由二维理论到三维理论的发展历程。介绍了三维裂纹端部场的K-Tz和J-Tz双参数描述、K-T-Tz和J-QT-Tz三参数描述,以及三维弹塑性断裂准则和三维疲劳裂纹闭合模型等。通过具体的实例介绍了三维断裂理论在航空结构损伤容限分析中的应用。正确考虑离面约束和面内约束对裂纹端部场以及材料断裂和疲劳裂纹扩展性能的影响,能够发展航空结构损伤容限可预测设计能力。     

Three-dimensional cohesive fracture modeling of non-planar crack growth using adaptive FE technique

In this paper, the three-dimensional adaptive finite element modeling is presented for cohesive fracture analysis of non-planer crack growth. The technique is performed based on the Zienkiewicz–Zhu error estimator by employing the modified superconvergent patch recovery procedure for the stress recovery. The Espinosa–Zavattieri bilinear constitutive equation is used to describe the cohesive tractions and displacement jumps. The 3D cohesive fracture element is employed to simulate the crack growth in a non-planar curved pattern. The crack growth criterion is proposed in terms of the principal stress and its direction. Finally, several numerical examples are analyzed to demonstrate the validity and capability of proposed computational algorithm. The predicted crack growth simulation and corresponding load-displacement curves are compared with the experimental and other numerical results reported in literature.     

2D and 3D finite element analysis of crack growth under compressive residual stress field

In the present study, plasticity induced crack closure (PICC) concept and three dimensional (3D) finite element method (FEM) were used to study the effect of compressive residual stress field on the fatigue crack growth from a hole. Furthermore, a new methodology on the basis of a correction factor was presented to increase the PICC precision. The result obtained was compared to two dimensional (2D) FEM, superposition method and Liu’s experimental data. To simulate the elasto-plastic behavior of the material, isotropic hardening was assumed and the Von-Mises yield criterion was implemented. A 3D mesh was built using eight-node hexahedral elements and one half of the specimen was modeled. The simulation results were fairly well correlated with experimental data. Furthermore, the 3D elasto-plastic FEM predicted a slightly smaller fatigue life than a 2D plane stress FEM. Applying the modified PICC method reduces the 3D FEM fatigue life prediction errors.     

Three dimensional microstructural effects on plane strain ductile crack growth

Ductile crack growth under mode I, plane strain, small scale yielding conditions is analyzed. Overall plane strain loading is prescribed, but a full 3D analysis is carried out to model three dimensional microstructural effects. An elastic-viscoplastic constitutive relation for a porous plastic solid is used to model the material. Two populations of second-phase particles are represented, large inclusions with low strength, which result in large voids near the crack tip at an early stage, and small second-phase particles, which require large strains before cavities nucleate. The larger inclusions are represented discretely and the effects of different three dimensional distributions on the crack path and on the overall crack growth rate are analyzed. For comparison purposes, a two dimensional distribution of cylindrical inclusions is analyzed. Crack growth occurs off the initial crack plane in all 3D computations, whereas straight ahead crack growth occurs with the two dimensional cylindrical inclusions. As a consequence, the three dimensional distributions of spherical inclusions exhibit an increased crack growth resistance as compared to the two dimensional distribution of cylindrical inclusions.     

Crack propagation in concrete: Linear elastic fracture mechanics and boundary element method

Crack propagation in concrete is investigated by application of Linear Elastic Fracture Mechanics (LEFM) and the Boundary Element Method (BEM). Acoustic Emission (AE) measurement is made for detecting the load level at which unstable crack initiation occurs in a notched beam under bending. The corresponding critical stress intensity factor K_Ic does not vary appreciably with notch depth. This result is in contrast to that found by the conventional procedure.The critical stress intensity factor K_Ic is employed to analyze the process of crack propagation in an arbitraty direction. The BEM prediction is in good agreement with the experiments.The detection and assessment of crack initiation in concrete structures is of considerable technical interest.     

表面裂纹准静态扩展的有限元模拟计算

由于现有断裂准则未能考虑表面裂纹前沿各点实际断裂阻力的变化，所作的表面裂纹准静态扩展模拟计算不能很好地预测几何形貌的真实变化规律。本文提出了一个考虑应力状态不同对实际断裂阻力影响的变阻力断裂准则；并以该准则作为表面裂纹前沿局部断裂扩展的判据，通过三维有限元法进行了表面裂纹准静态扩展的模拟计算；还将模拟预测值与用多试件法测得的实验值进行比较，验证了变阻力断裂准则作为裂纹稳态扩展判据的有效性     

Fracture Mechanics of Corrosion Cracking in Concrete by Acoustic Emission

Nucleation of microcracks can be detected and analyzed by acoustic emission (AE), by which crack kinematics of locations, types and rientations are quantitatively estimated. The procedure was applied to clarify mechanisms of corrosion cracking. Based on fracture mechanics, numerical analysis was conducted by the boundary element method (BEM). Relations between the stress intensity factors and crack types were investigated by BEM solutions. In experiments, four types of crack patterns were nucleated by employing expansive agent. Following the surface crack, the diagonal crack and/or the horizontal crack propagated. The internal crack extended after the surface crack was terminated. Depending on the crack types, contributions of mode-I and mode-II were varied. According to AE results, four crack patterns observed differently consisted of tensile, mixed-mode and shear cracks. It is demonstrated that mechanisms of corrosion cracking in concrete are dominantly of mode-I failure along with a minor contribution of mixed-mode and mode-II.     

Numerical and experimental determination of three-dimensional multiple crack growth in fatigue

This work is concerned with the assessment of propagation of multiple fatigue cracks in three-dimensions. Computational modelling of fatigue crack propagation is made together with detection and monitoring of the crack shape development. The boundary element method (BEM) is used for automating the modelling of crack propagation in linear elastic as well as elastic–plastic regimes. Strain at several positions on the specimen surface near the crack mouth is measured to monitor crack initiation, shape development and closure levels. Examples are provided to validate the model by comparing the experimental results with those obtained by numerical predictions.