CTS试样I/II混合型断裂特性数值计算研究

采用CTS试样研究I/II混合型断裂特性计算裂纹前缘应力强度因子时可采用解析公式,一旦裂纹发生扩展,解析公式便不再适用。文中采用有限元法研究紧凑拉伸剪切（CTS）试样在I/II平面混合型加载下的裂纹扩展行为。采用ANSYS建立CTS试样I/II混合型测试系统有限元模型,为模拟真实受力状态,在CTS试样-销-扇型夹具以及扇型夹具-销-U型夹具之间分别建立接触对进行接触力学分析。通过与解析公式结果进行对比验证了该数值方法的可靠性。采用最大环向应力准则（MTS）,模拟了CTS试样不同加载角度下的裂纹扩展路径,获得了裂纹扩展路径中应力强度因子随裂纹长度的变化曲线,解释了裂纹扩展路径不与外载荷方向垂直的原因。结合文中计算结果,在CTS试样I/II混合型裂纹扩展速率实验测得裂纹长度与寿命的关系曲线a-N的基础上,便可得到材料I/II型混合型裂纹扩展速率曲线。     

Investigation of mixed-mode stress intensity factors for nonhomogeneous materials using an interaction integral method

An interaction (energy) integral is derived for the computation of mixed-mode stress intensity factors (SIFs) in nonhomogeneous materials with continuous or discontinuous properties. This method is based on a conservation integral that relies on two admissible mechanical states (actual and auxiliary fields). In general, the interaction energy contour integral is converted into an equivalent domain integral in numerical computations. It can be seen from the equivalent domain integral, the integrand does not involve any derivatives of material properties. Moreover, the formulation can be proved valid even when the integral domain contains material interfaces. Therefore, it is not necessary to limit the material properties to be continuous for the present method. Due to these advantages the application range of the interaction integral method can be greatly enlarged. The numerical implementation of the derived expression is combined with the extended finite element method (XFEM). Using this method, the influences of material properties on the mixed-mode SIFs are investigated for four types of material properties selected in this work. Numerical results show that the mechanical properties and their first-order derivatives can affect mode I and II SIFs greatly, while the higher-order derivatives affect the SIFs very slightly.     

MESHLESS METHOD FOR 2D MIXED-MODE CRACK PROPAGATION BASED ON VORONOI CELL

A meshless method integrated with linear elastic fracture mechanics (LEFM) is presented for 2D mixed-mode crack propagation analysis. The domain is divided automatically into sub-domains based on Voronoi cells, which are used for quadrature for the potential energy. The continuous crack propagation is simulated with an incremental crack-extension method which assumes a piecewise linear discretization of the unknown crack path. For each increment of the crack extension, the meshless method is applied to carry out a stress analysis of the cracked structure. The J-integral, which can be decomposed into mode I and mode II for mixed-mode crack, is used for the evaluation of the stress intensity factors (SIFs). The crack-propagation direction, predicted on an incremental basis, is computed by a criterion defined in terms of the SIFs. The flowchart of the proposed procedure is presented and two numerical problems are analyzed with this method. The meshless results agree well with the experimental ones, which validates the accuracy and efficiency of the method.     

Measurement of stress intensity factors by digital camera

This paper developed a least-squares method to evaluate the mixed-mode stress intensity factors (SIFs) of the isotropic material using the computer images from the digital-camera experiment. This experiment measures the crack opening displacement (COD) and then evaluates SIFs by using the least-squares method. The attractions of this method include: (1) specimen preparation and experiment procedures are not complicated and (2) the isolation of the micro-vibration is not necessary in the experiment. Both finite element simulations and laboratory experiments were applied to validate the current least-squares method with acceptable accuracy, if the even terms of the Irwin’s equation are removed.     

A new boundary element approach of modeling singular stress fields of plane V-notch problems

In this paper, a new boundary element (BE) approach is proposed to determine the singular stress field in plane V-notch structures. The method is based on an asymptotic expansion of the stresses in a small region around a notch tip and application of the conventional BE in the remaining region of the structure. The evaluation of stress singularities at a notch tip is transformed into an eigenvalue problem of ordinary differential equations that is solved by the interpolating matrix method in order to obtain singularity orders (degrees) and associated eigen-functions of the V-notch. The combination of the eigen-analysis for the small region and the conventional BE analysis for the remaining part of the structure results in both the singular stress field near the notch tip and the notch stress intensity factors (SIFs).Examples are given for V-notch plates made of isotropic materials. Comparisons and parametric studies on stresses and notch SIFs are carried out for various V-notch plates. The studies show that the new approach is accurate and effective in simulating singular stress fields in V-notch/crack structures.     

基于p型有限元法和围线积分法计算复合型应力强度因子

传统有限元法由于采用低阶插值计算应力强度因子时,需要划分的网格数较多,收敛速度较慢,得到的应力强度因子精度不足。p型有限元法在网格确定时通过增加插值多项式的阶数来提高计算精度,具有网格划分少、收敛速度快、精度高、自适应能力强等特点。本文采用基于p型有限元法的有限元计算软件StressCheck计算得到应力场和位移场,并由围线积分法导出混合型应力强度因子(SIFs)。通过几个经典算例,分析了围线的选择对计算精度的影响,计算了不同裂纹长度、不同裂纹角度和裂纹在应力集中区域不同位置时的应力强度因子。并将数值结果、理论解与文献中其他数值计算方法所得的部分结果进行了对比分析,结果表明自由度数不大于7000时,导出的应力强度因子相对误差最大不超过1.2%,数值解表现出较高的精度及数值稳定性。     

Computations of modes I and II stress intensity factors of sharp notched plates under in-plane shear and bending loading by the fractal-like finite element method

The fractal-like finite element method (FFEM) is used to compute the stress intensity factors (SIFs) for different configurations of cracked/notched plates subject to in-plane shear and bending loading conditions. In the FFEM, the large number of unknown variables in the singular region around a notch tip is reduced to a small set of generalised co-ordinates by performing a fractal transformation using global interpolation functions. The use of exact analytical solutions of the displacement field around a notch tip as the global interpolation functions reduces the computational cost significantly and neither post-processing technique to extract SIFs nor special singular elements to model the singular region are required. The results of numerical examples of various configurations of cracked/notched plates are presented and validated via published data. Also, new results for cracked/notched plate problems are presented. These results demonstrate the accuracy and efficiency of the FFEM to compute the SIFs for notch problems under in-plane shear and bending loading conditions.     

Three-dimensional mixed-mode stress intensity factors for cracks in functionally graded materials using enriched finite elements

Three-dimensional enriched finite elements are used to compute mixed-mode stress intensity factors (SIFs) for three-dimensional cracks in elastic functionally graded materials (FGMs) that are subject to general mixed-mode loading and constraint conditions. The method, which advantageously does not require special mesh configuration/modifications and post-processing of finite element results, is an enhancement of previous developments applied so far on isotropic homogeneous and isotropic interface cracks. The spatial variation of FGM material properties is taken into account at the level of element integration points. To validate the developed method, two- and three-dimensional mixed-mode fracture problems are selected from the literature for comparison. Two-dimensional cases include: inclined central crack in a large FGM medium under uniform tensile strain and stress loadings, a slanted crack in a finite-size FGM plate under exponentially varying tensile stress loading and an edge crack in a finite-size plate under shear traction load. The three-dimensional example models a deflected surface crack in a finite-size FGM plate under uniform tensile stress loading. Comparisons between current results and those from analytical and other numerical methods yield good agreement. Thus, it is concluded that the developed three-dimensional enriched finite elements are capable of accurately computing mixed-mode fracture parameters for cracks in FGMs.     

孔洞对平面裂纹应力强度因子影响分析的广义参数Williams单元

本文采用圆形奇异区广义参数Williams单元（W单元）建立了中心裂纹与圆孔共存的平面应力模型,奇异区外围利用ABAQUS有限元软件自动网格离散技术与FORTRAN95编程前处理相结合,克服了自主编程中网格离散的局限性．算例分析了圆孔位置和几何参数对I-II混合型裂纹尖端应力强度因子(SIFs)的影响,并与扩展有限元法(XFEM)计算结果进行比较．结果表明：靠近圆孔一侧的裂尖SIFs大于远离圆孔一侧的裂尖SIFs;控制圆孔左边缘到裂纹中心的距离,则两侧裂尖SIFs随圆孔半径的增大而增大;圆孔中心与裂纹中心水平距离越远,圆孔对裂纹扩展的影响越小．同时,基于圆形奇异区的W单元直接计算得到的裂尖SIFs与扩展有限元法得到的解吻合较好,证明了W单元对奇异区离散形状不敏感,且具有高效率和高精度．     

Application of scaled boundary finite element method in static and dynamic fracture problems

The scaled boundary finite element method (SBFEM) is a recently developed numerical method combining advantages of both finite element methods (FEM) and boundary element methods (BEM) and with its own special features as well. One of the most prominent advantages is its capability of calculating stress intensity factors (SIFs) directly from the stress solutions whose singularities at crack tips are analytically represented. This advantage is taken in this study to model static and dynamic fracture problems. For static problems, a remeshing algorithm as simple as used in the BEM is developed while retaining the generality and flexibility of the FEM. Fully-automatic modelling of the mixed-mode crack propagation is then realised by combining the remeshing algorithm with a propagation criterion. For dynamic fracture problems, a newly developed series-increasing solution to the SBFEM governing equations in the frequency domain is applied to calculate dynamic SIFs. Three plane problems are modelled. The numerical results show that the SBFEM can accurately predict static and dynamic SIFs, cracking paths and load-displacement curves, using only a fraction of degrees of freedom generally needed by the traditional finite element methods.The project supported by the National Natural Science Foundation of China (50579081) and the Australian Research Council (DP0452681)The English text was polished by Keren Wang.     

直接计算应力强度因子的扩展有限元法

系统地给出了直接计算应力强度因子的扩展有限元法。该方法以常规有限元法为基础,利用单位分解法思想,通过在近似位移表达式中增加能够反映裂纹面的不连续函数及反映裂尖局部特性的裂尖渐进位移场函数,间接体现裂纹面的存在,从而无需使裂纹面与有限元网格一致,无需在裂尖布置高密度网格,也不需要后处理就可以直接计算出应力强度因子,并且大大简化了前后处理工作。最后通过两个简单算例验证了该方法的精度,分析了影响计算结果的因素,并与采用J积分计算的应力强度因子作了对比,得出了两种方法计算精度相当的结论。     

Determination of mixed-mode stress intensity factors, fracture toughness, and crack turning angle for anisotropic foam material

A numerical and experimental investigation for determining mixed-mode stress intensity factors, fracture toughness, and crack turning angle for BX-265 foam insulation material, used by NASA to insulate the external tank (ET) for the space shuttle, is presented. BX-265 foam is a type of spray-on foam insulation (SOFI), similar to the material used to insulate attics in residential construction. This cellular material is a good insulator and is very lightweight. Breakup of segments of this foam insulation on the shuttle ET impacting the shuttle thermal protection tiles during liftoff is believed to have caused the space shuttle Columbia failure during re-entry. NASA engineers are interested in understanding the processes that govern the breakup/fracture of this material from the shuttle ET. The foam is anisotropic in nature and the required stress and fracture mechanics analysis must include the effects of the direction dependence on material properties. Material testing at NASA Marshall Space Flight Center (MSFC) has indicated that the foam can be modeled as a transversely isotropic material. As a first step toward understanding the fracture mechanics of this material, we present a general theoretical and numerical framework for computing stress intensity factors (SIFs), under mixed-mode loading conditions, taking into account the material anisotropy. We present SIFs for middle tension – M(T) – test specimens, using 3D finite element stress analysis (ANSYS) and FRANC3D fracture analysis software. SIF values are presented for a range of foam material orientations. Mode I fracture toughness of the material is determined based on the SIF value at failure load. We also present crack turning angles for anisotropic foam material under mixed-mode loading. The results represent a quantitative basis for evaluating the strength and fracture properties of anisotropic foam insulation material.     

Numerical analysis of the dual actuator load test applied to fracture characterization of bonded joints

The dual actuator load test was numerically analyzed in order to assess its adequacy for fracture characterization of bonded joints under different mixed-mode loading conditions. This test enables asymmetric loading of double cantilever beam specimens, thus providing a large range of mixed-mode combinations. A new data reduction scheme based on specimen compliance, beam theory and crack equivalent concept was proposed to overcome several difficulties inherent to the test. The method assumes that the dual actuator test can be viewed as a combination of the double cantilever beam and end loaded split tests, which are used for pure modes I and II fracture characterization, respectively. A numerical analysis including a cohesive mixed-mode damage model was performed considering different mixed-mode loading conditions to evaluate the test performance. Some conclusions were drawn about the advantages and drawbacks of the test.     

A domain-independent integral for computation of stress intensity factors along three-dimensional crack fronts and edges by BEM

The present work deals with an evaluation of stress intensity factors (SIFs) along straight crack fronts and edges in three-dimensional isotropic elastic solids. A new numerical approach is developed for extraction, from a solution obtained by the boundary element method (BEM), of those SIFs, which are relevant for a failure assessment of mechanical components. In particular, the generalized SIFs associated to eigensolutions characterized by unbounded stresses at a neighbourhood of the crack front or a reentrant edge and also that associated to T-stress at the crack front can be extracted. The method introduced is based on a conservation integral, called H-integral, which leads to a new domain-independent integral represented by a scalar product of the SIF times some element shape function defined along the crack front or edge. For sufficiently small element lengths these weighted averages of SIFs give reasonable pointwise estimation of the SIFs. A proof of the domain integral independency, based on the bi-orthogonality of the classical two-dimensional eigensolutions associated to a corner problem, is presented. Numerical solutions of two three-dimensional problems, a crack problem and a reentrant edge problem, are presented, the accuracy and convergence of the new approach for SIF extraction being analysed.     

Determination of V-notch SIFs in multi-material anisotropic wedges by digital correlation experiments

In this paper, theoretical formulations based on the Stroh’s complex function approach were used to find the displacement field and H-integral of a sharp V-notch formed from several anisotropic materials. Displacements from the image-correlation experiments are then substituted into the least-squares formulation to find V-notch stress intensity factors (SIFs) in multi-material anisotropic wedges. Validations using the H-integral indicate that the experimental SIFs evaluated from the proposed method of acceptable accuracy. The major advantage is that the proposed method only requires displacements inside the specimen, and displacements near the notch tip, specimen boundaries, or notch surfaces are not necessary.     

The thermoelastic equilibrium of a transversally isotropic medium with an elliptic crack under symmetric loading

The stress intensity factors (SIFs) are evaluated for flat elliptical cracks located in a transversally isotropic material (cracks are assumed perpendicular to the transtropy axis) under an arbitrary load and symmetric temperature. The SIFs for an elliptical crack in a transversally isotropic medium are determined using the formulas (derived by the author in his previous studies) of transition from an isotropic to transversally isotropic material and the relative problem for an isotropic medium. It is proved that these formulas can be employed for an arbitrary homogeneous transversally isotropic material (no matter whether the roots of some characteristic equation of the material are real or complex) with an arbitrary flat crack or a system of coplanar flat cracks, including elliptical ones, under an arbitrary load and symmetric temperature. A transversally isotropic material with two coplanar elliptical cracks is considered as an illustrative example. The dependences of the SIFs on the parameters of cracks and their arrangement at a decreasing temperature are presented. S. P. Timoshenko Institute of Mechanics, National Academy of Sciences of Ukraine Kiev. Translated from Prikladnaya Mekhanika, Vol. 36, No. 4, pp. 96–105, April, 2000.     

Application of hypersingular integral equation method to three-dimensional crack in electromagnetothermoelastic multiphase composites

A three-dimensional crack problem in electromagnetothermoelastic multiphase composites (EMTE-MCs) under extended loads is investigated in this paper. Using Green’s functions, the extended general displacement solutions are obtained by the boundary element method. This crack problem is reduced to solving a set of hypersingular integral equations coupled with boundary integral equations, in which the unknown functions are the extended displacement discontinuities. Then, the behavior of the extended displacement discontinuities around the crack front terminating at the interface is analyzed by the main-part analysis method of hypersingular integral equations. Analytical solutions for the extended singular stresses, the extended stress intensity factors (SIFs) and the extended energy release rate near the crack front in EMTE-MCs are provided. Also, a numerical method of the hypersingular integral equations for a rectangular crack subjected to extended loads is put forward with the extended displacement discontinuities approximated by the product of basic density functions and polynomials. In addition, distributions of extended SIFs varying with the shape of the crack are presented. The results show that the present method accurately yields smooth variations of extended SIFs along the crack front.     

垂直磁压电材料界面三维裂纹的超奇异积分法

应用有限部积分概念和广义位移基本解,垂直于磁压电双材料界面三维复合型裂纹问题被转 化为求解一组以裂纹表面广义位移间断为未知函数的超奇异积分方程问题. 进而,通过主部 分析法精确地求得裂纹尖端光滑点附近的奇性应力场解析表达式. 然后,通过将裂纹表面 位移间断未知函数表达为位移间断基本密度函数与多项式之积,使用有限部积分法对超奇异 积分方程组建立了数值方法. 最后,通过典型算例计算,讨论了广义应力强度因子的变化规 律.     

An interaction integral method for 3D curved cracks in nonhomogeneous materials with complex interfaces

This work derives an interaction integral for the computation of mixed-mode stress intensity factors (SIFs) in three-dimensional (3D) nonhomogeneous materials with continuous or discontinuous properties. The present method is based on a two-state integral by the superposition of actual and auxiliary fields. In 3D domain formulation of the interaction integral derived here, the integrand does not involve any derivatives of material properties. Furthermore, the formulation can be proved to be still valid even when the integral domain contains material interfaces. Therefore, it is not necessary to limit the material properties to be continuous for the present formulation. On account of these advantages, the application range of the interaction integral can be greatly enlarged. This method in conjunction with the finite element method (FEM) is employed to solve several representative fracture problems. According to the comparison between the results and those from the published lectures, good agreement demonstrates the validation of the interaction integral. The results show that the present interaction integral is domain-independent for nonhomogeneous materials with interfaces.     

Contour integrals for mixed-mode crack analysis: effect of nonsingular terms

An integral formulation for computing the nonsingular stresses (NSS) in a cracked body under mixed-mode static and dynamic loads is presented. The reciprocity theorems are applied to find the integral formula. The auxiliary fields are selected to eliminate the singular terms in the asymptotic expansion of the stresses near the crack tip. For elastodynamic crack problems, the integral representation of the NSS is presented in both the time and Laplace transform domain. Required variables along the integration path and region enclosed by the integration contour are obtained from the boundary element analysis. Influence of the NSS on predicting the crack growth direction is investigated for cracks under mixed-mode load conditions.