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.     

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.     

Mixed-mode fracture analysis of orthotropic functionally graded materials under mechanical and thermal loads

Mixed-mode fracture problems of orthotropic functionally graded materials (FGMs) are examined under mechanical and thermal loading conditions. In the case of mechanical loading, an embedded crack in an orthotropic FGM layer is considered. The crack is assumed to be loaded by arbitrary normal and shear tractions that are applied to its surfaces. An analytical solution based on the singular integral equations and a numerical approach based on the enriched finite elements are developed to evaluate the mixed-mode stress intensity factors and the energy release rate under the given mechanical loading conditions. The use of this dual approach methodology allowed the verifications of both methods leading to a highly accurate numerical predictive capability to assess the effects of material orthotropy and nonhomogeneity constants on the crack tip parameters. In the case of thermal loading, the response of periodic cracks in an orthotropic FGM layer subjected to transient thermal stresses is examined by means of the developed enriched finite element method. The results presented for the thermally loaded layer illustrate the influences of the material property gradation profiles and crack periodicity on the transient fracture mechanics parameters.     

Transient thermal shock fracture analysis of functionally graded piezoelectric materials by the extended finite element method

Transient thermal dynamic analysis of stationary cracks in functionally graded piezoelectric materials (FGPMs) based on the extended finite element method (X-FEM) is presented. Both heating and cooling shocks are considered. The material properties are supposed to vary exponentially along specific direction while the crack-faces are assumed to be adiabatic and electrically impermeable. A dynamic X-FEM model is developed in which both Crank–Nicolson and Newmark time integration methods are used for calculating transient responses of thermal and electromechanical fields respectively. The generalized dynamic intensity factors for the thermal stresses and electrical displacements are extracted by using the interaction integral. The accuracy of the developed approach is verified numerically by comparing the calculated results with reference solutions. Numerical examples with mixed-mode crack problems are analyzed. The effects of the crack-length, poling direction, material gradation, etc. on the dynamic intensity factors are investigated. It shows that the transient dynamic crack behaviors under the cooling shock differ from those under the heating shock. The influence of the thermal shock loading on the dynamic intensity factors is significant.     

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

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

A new variable-order singular boundary element for calculating stress intensity factors in three-dimensional elasticity problems

A new three-dimensional variable-order singular boundary element has been constructed for stress analysis of three-dimensional interface cracks and internal material junctions. The singular fields in the vicinity of crack front or junction have been accurately represented by the singular elements by taking account the variable order of singularities and the angular profiles of field variables. Both the singular stress fields and displacement fields are independently formulated by the element’s shape functions. Different kinds of displacement formulations are investigated. The formulation combining singular and linear terms is found to be the most accurate one. The mixed-mode stress intensity factors are treated as nodal unknowns. The variation of stress intensity factors along the line of singularity can be obtained directly from the final system of equations and thus no post processing, such as three-dimensional J-integral or domain integral, is necessary. Numerical examples involving stress singularity, such as penny-shaped cracks in homogeneous and dissimilar material interface, plates with through-thickness cracks, and a dissimilar inclusion, are investigated. The analysis results are in good agreement with those reported in the literature.     

基于扩展有限元的应力强度因子的位移外推法

针对平面裂纹问题,阐述了扩展有限元法的单元位移模式、推导了扩展有限元法的控制方程、介绍了特殊单元的数值积分技术.基于最小二乘法,建立了应力强度因子位移外推法的计算公式.利用MATLAB编写计算程序,对平面裂纹问题用扩展有限元法进行了计算.基于扩展有限元法的计算结果,分别利用位移外推法和相互作用积分法,对平面裂纹的应力强度因子进行了计算.计算结果表明,位移外推法比相互作用积分法能更方便和准确地计算平面裂纹的应力强度因子.     

Computation of thermal fracture parameters for orthotropic functionally graded materials using Jk-integral

This article introduces a computational method based on the J_k-integral for mixed-mode fracture analysis of orthotropic functionally graded materials (FGMs) that are subjected to thermal stresses. The generalized definition of the J_k-integral is recast into a domain independent form composed of line and area integrals by utilizing the constitutive relations of plane orthotropic thermoelasticity. Implementation of the domain independent J_k-integral is realized through a numerical procedure developed by means of the finite element method. The outlined computational approach enables the evaluation of the modes I and II stress intensity factors, the energy release rate, and the T-stress. The developed technique is validated numerically by considering two different problems, the first of which is the problem of an embedded crack in an orthotropic FGM layer subjected to steady-state thermal stresses; and the second one is that of periodic cracks under transient thermal loading. Comparisons of the mixed-mode stress intensity factors evaluated by the J_k-integral based method to those calculated through the displacement correlation technique (DCT) and to those available in the literature point out that, the proposed form of the J_k-integral possesses the required domain independence and leads to numerical results of high accuracy. Further results are presented to illustrate the influences of the geometric and material constants on the thermal fracture parameters.     

In-plane stress analysis of an orthotropic plane containing multiple defects

The solution of Volterra type climb and glide edge dislocations is utilized to formulate integral equations for an orthotropic homogeneous infinite plane weakened by multiple smooth cracks and/or cavities. Cavities are considered as closed curved cracks without singularity. The integral equations are of Cauchy singular type which are converted to hypersingular integral equations. These equations are then solved numerically to determine stress intensity factors for cracks and hoop stress on the cavities. The results for isotropic and orthotropic planes are compared with available solutions in literature and excellent agreement is observed. The formulation allows stress analysis of orthotropic planes with several arbitrarily oriented cracks and cavities.     

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

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

An embedded crack in a functionally graded orthotropic coating bonded to a homogeneous substrate under a frictional Hertzian contact

In this paper, we consider the elasto-static problem of an embedded crack in a graded orthotropic coating bonded to a homogeneous substrate subject to statically applied normal and tangential surface loading. The crack direction is parallel to the free surface. The coating is graded in the thickness direction and is orthogonal to the crack direction. This coating is modelled as a non-homogeneous medium with an orthotropic stress–strain law. The equivalent crack surface stresses are first obtained and substituted in the plane elasticity equations. Using integral transforms, the governing equations are converted into singular integral equations which are solved numerically to yield the displacement field as well as the crack-tip stress intensity factors. This study presents a complete theoretical formulation for the problem in the static case. A numerical predictive capability for solving the singular integral equations and computing the crack-tip stress intensity factors is proposed. Since the loading is compressive, a previously developed crack-closure algorithm is applied to avoid interpenetration of the crack faces. The main objective of the paper is to investigate the effects of the material orthotropy and non-homogeneity of the graded coating on the crack-tip stress intensity factors, with and without using the crack-closure algorithm, for the purpose of gaining better understanding on the behavior and design of graded coatings.     

A domain-independent interaction integral for magneto-electro-elastic materials

Magneto-electro-elastic (MEE) materials usually consist of piezoelectric (PE) and piezomagnetic (PM) phases. Between different constituent phases, there exist lots of interfaces with discontinuous MEE properties. Complex interface distribution brings a great difficulty to the fracture analysis of MEE materials since the present fracture mechanics methods can hardly solve the fracture parameters efficiently of a crack surrounded by complex interfaces. This paper develops a new domain formulation of the interaction integral for the computation of the fracture parameters including stress intensity factors (SIFs), electric displacement intensity factor (EDIF) and magnetic induction intensity factor (MIIF) for linear MEE materials. The formulation derived here does not involve any derivatives of material properties and moreover, it can be proved that an arbitrary interface in the integral domain does not affect the validity and the value of the interaction integral. Namely, the interaction integral is domain-independent for material interfaces and thus, its application does not require material parameters to be continuous. Due to this advantage, the interaction integral becomes an effective approach for extracting the fracture parameters of MEE materials with complex interfaces. Combined with the extended finite element method (XFEM), the interaction integral is employed to solve several representative problems to verify its accuracy and domain-independence. Good results show the effectiveness of the present method in the fracture analysis of MEE materials with continuous and discontinuous properties. Finally, the particulate MEE composites composed of PE and PM phases are considered and four schemes of different property-homogenization level are proposed for comparing their effectiveness.     

求解混合型裂纹应力强度因子的围线积分法

本文用复变函数理论推导出裂纹的辅助场，并用Betti功互等定理给出求解混合型裂纹应力强度因子的远场围绕积分法．此方法与积分路径的选择无关，用有限元法计算出远离裂纹尖端的位移场和应力场，就可通过计算绕裂端的围线积分，精确地给出混合型裂纹的应力强度因子KⅠ和KⅡ的数值解．     

Boundary element analysis of anisotropic Kirchhoff plates

In this paper, the radial integration method is used to obtain a boundary element formulation without any domain integral for general anisotropic plate bending problems. Two integral equations are used and the unknown variables are assumed to be constant along each boundary element. The domain integral which arises from a transversely applied load is exactly transformed into a boundary integral by a radial integration technique. Uniformly and linearly distributed loads are considered. Several computational examples concerning orthotropic and general anisotropic plate bending problems are presented. The results show good agreement with analytical and finite element results available in the literature.     

A domain-independent interaction integral for fracture analysis of nonhomogeneous piezoelectric materials

Piezoelectric materials and structures contain more or less electromechanical interfaces in engineering applications. It is difficult to obtain the fracture parameters efficiently of the piezoelectric materials with complex interfaces. This paper presents a domain-independent interaction integral for material nonhomogeneity and discontinuity which can be used for solving the stress intensity factors (SIFs) and the electric displacement intensity factor (EDIF) of piezoelectric materials with complex interfaces efficiently. The interaction integral is based on the J-integral by superimposition of two admissible states and the present formulation does not involve any derivatives of mechanical and electric properties. Moreover, it is proved that the interface in the integral domain does not affect the value of the interaction integral and thus, the present method does not require electromechanical parameters of piezoelectric materials to be continuous. The interaction integral method combined with the extended finite element method (XFEM) is used to investigate the influences of material continuity on the SIF and the EDIF and the results show that the material parameters and their first-order derivatives affect both the SIF and the EDIF greatly, while the higher-order derivatives affect both of them slightly.     

SH波在正交各向异性功能梯度无限长条中心裂缝处的散射

研究了正交各向异性功能梯度材料无限长条中心裂缝对SH波的散射问题,为方便起见,材料两个方向的剪切模量和密度假定为指数模型.通过Fourier积分变换,将问题转化为对偶积分方程的求解.然后,用Cop-son方法求解对偶积分方程,定义了标准动应力强度因子,通过数值算例,讨论了在SH波作用下,裂缝尖端的标准动应力强度因子与入射波的频率、材料参数之间的关系.     

Dependence of stress intensity factors on elastic constants for cracks in an orthotropic bimaterial with a thin film

An interface crack and a subinterface crack in an orthotropic bimaterial structure consisting of a thin film and a half plane substrate are analyzed. The orthotropic bimaterial structure is subjected to compressive load and bending moment per unit thickness. Complete expressions of stress intensity factors for the two cracks are obtained based on the path independence of the J integral, apart from one dimensionless parameter undetermined each. The dependence of the dimensionless parameters on material constants is examined. A reduction of the number of necessary material parameters for the parameters is made based upon the modified Stroh formalism. The explicit dependence of the dimensionless parameters on one orthotropic parameter for the film is determined by using the orthotropy rescaling technique. Variations of the dimensionless parameters with the other material parameters are also obtained through numerical computations.     

On interface crack in laminated anisotropic medium

In this paper, the interface crack problems of a multilayered anisotropic medium under a state of generalized plane deformation are considered within the framework of anisotropic theory. A general solution procedure is introduced such that it can be uniformly applied to media with transversely isotropic, orthotropic, monoclinic, etc. layers. The problem is reduced to the solution of a system of singular integral equations by means of Fourier transform method and the stiffness matrix formulation. A Jacobi polynomial technique is then used to solve the integral equations numerically. The stress intensity factors are provided. The stress intensity factors have been calculated numerically and displayed graphically.