Green's functions for the stress intensity factor evolution in finite cracks in orthotropic materials

The elastodynamic response of an infinite orthotropic material with finite crack under concentrated loads is examined. Solution for the stress intensity factor history around the crack tips is found. Laplace and Fourier transforms are employed to solve the equations of motion leading to a Fredholm integral equation on the Laplace transform domain. The dynamic stress intensity factor history can be computed by numerical Laplace transform inversion of the solution of the Fredholm equation. Numerical values of the dynamic stress intensity factor history for some example materials are obtained. This solution can be used as a Green's function to solve dynamic problems involving fini te cracks.     

Thermal fracture analysis of a functionally graded orthotropic strip with a crack

This paper is concerned with the thermal fracture problem of a functionally graded orthotropic strip, where the crack is situated parallel to the free edges. All the material properties are assumed to be dependent only on the coordinate y (perpendicular to the crack surfaces). By using Fourier transform, the thermoelastic problem is reduced to those that involve a system of singular integral equations. Numerical results are presented to show the effects of the crack position and the material distribution on the thermal stress intensity factors.     

Dynamic stress intensity factor due to concentrated loads on a propagating semi-infinite crack in orthotropic materials

The elastodynamic response of an infinite orthotropic material with a semi-infinite crack propagating at constant speed under the action of concentrated loads on the crack faces is examined. Solution for the stress intensity factor history around the crack tip is found for the loading modes I and II. Laplace and Fourier transforms along with the Wiener-Hopf technique are employed to solve the equations of motion. The asymptotic expression for the stress near the crack tip is analyzed which lead to a closed-form solution of the dynamic stress intensity factor. It is found that the stress intensity factor for the propagating crack is proportional to the stress intensity factor for a stationary crack by a factor similar to the universal function k(v) from the isotropic case. Results are presented for orthotropic materials as well as for the isotropic case.     

Mixed mode asymptotic crack tip fields in orthotropic materials: Derivation and range of dominance

In the first part of this work the asymptotic stress field distribution surrounding a crack in a generally orthotropic solid (i.e. one in which the material and loading axes do not coincide) is derived. The resulting crack tip stress intensity factors are also related to the energy release rate of the cracked solid. In any physical situation, however, the range of dominance of these asymptotic fields will be limited, and will depend on specific geometry and loading parameters. Thus, the second part of this work deals with the range of dominance of the derived stress fields in edge cracked bending loaded fiber reinforced composite plates. The lower and upper limits of the range of dominance of the solution are respectively determined from three-dimensional and two-dimensional full field finite element solutions. The comparison of full field solutions with the asymptotic result provides information on the latter's range of dominance. In all cases the effects of mixed mode loading are considered.     

Developing new enrichment functions for crack simulation in orthotropic media by the extended finite element method

New enrichment functions are proposed for crack modelling in orthotropic media using the extended finite element method (XFEM). In this method, Heaviside and near‐tip functions are utilized in the framework of the partition of unity method for modelling discontinuities in the classical finite element method. In this procedure, by using meshless based ideas, elements containing a crack are not required to conform to crack edges. Therefore, mesh generation is directly performed ignoring the existence of any crack while the method remains capable of extending the crack without any remeshing requirement. Furthermore, the type of elements around the crack‐tip remains the same as other parts of the finite element model and the number of nodes and consequently degrees of freedom are reduced considerably in comparison to the classical finite element method. Mixed‐mode stress intensity factors (SIFs) are evaluated to determine the fracture properties of domain and to compare the proposed approach with other available methods. In this paper, the interaction integral (M‐integral) is adopted, which is considered as one of the most accurate numerical methods for calculating stress intensity factors. Copyright © 2006 John Wiley & Sons, Ltd.     

Kinking of an interface crack in an orthotropic bimaterial

We investigated the asymptotic problem of a kinked interface crack in an orthotropic bimaterial under in‐plane loading conditions. The stress intensity factors at the tip of the kinked interface crack are described in terms of the stress intensity factors of the interface crack prior to the kink combined with a dimensionless matrix function. Using a modified Stroh formalism and an orthotropy rescaling technique, the matrix function was obtained from the solutions of the corresponding problem in transformed bimaterial. The effects of orthotropic and bimaterial parameters on the matrix function were examined. A reduction in the number of dependent material parameters on the matrix function was made using the modified Stroh formalism. Moreover, the explicit dependence of one orthotropic parameter on the matrix function was determined using an orthotropic rescaling technique. The effects of the other material parameters on the matrix function were numerically examined. The energy release rate was obtained for a kinked interface crack in an orthotropic bimaterial.     

Thermal stress intensity factors of an infinite orthotropic layer with a crack

Stress intensity factors are determined for a crack in an infinite orthotropic layer. The crack is situated parallel to the plane surfaces of the layer. Stresses are solved for two kinds of the boundary conditions with respect to temperature field. In the first problem, the upper surface of the layer is heated to maintain a constant temperature T ₀, while the lower surface is cooled to maintain a constant temperature –T ₀. In the other problem, uniform heat flows perpendicular to the crack. The surfaces of the crack are assumed to be insulated. The boundary conditions are reduced to dual integral equations using the Fourier transform technique. To satisfy the boundary conditions outside the crack, the difference in temperature at the crack surfaces and differences in displacements are expanded in a series of functions that vanish outside the crack. The unknown coefficients in each series are evaluated using the Schmidt method. Stress intensity factors are then calculated numerically for a steel layer that behaves as an isotropic material and for a tyrannohex layer that behaves as an orthotropic material.     

Fracture analysis of interfacial crack by global-local finite element

An eigenfunction expansion is used to formulate the global element on the crack tip. The global-local finite element method employs both conventional finite element and classical Rayleigh-Ritz kinematic approach. The hybrid Ritz method not only preserves the finite element modelling capability but adds the advantage of using prior information regarding the anticipate behaviour of the particular problem. Thus, it is able to achieve better accuracy with fewer elements in comparison with conventional finite element. Several examples relative to crack problems are presented to demonstrate the global-local finite element method.     

Dynamic Stress Intensity Factor of Interfacial Finite Cracks in Orthotropic Materials

The elastodynamic response of an infinite non-homogeneous orthotropic material with an interfacial finite crack under distributed normal and shear impact loads is examined. Solution for the stress intensity factor history around the crack tips is found. Laplace and Fourier transforms are employed to solve the equations of motion leading to a Fredholm integral equation on the Laplace transform domain. The dynamic stress intensity factor history can be computed by numerical Laplace transform inversion of the solution of the Fredholm equation. Numerical values of the dynamic stress intensity factor history for some materials are obtained. Interfacial cracks between two different materials and between two pieces of the same material but different fiber orientation are considered. Bimaterial formulation of a crack problem is shown to converge to the mono-material formulation, derived independently, in the limiting case when both materials are the same.     

Finite element analysis of fretting crack propagation

In this work, the finite elements method (FEM) is used to analyse the growth of fretting cracks. FEM can be favourably used to extract the stress intensity factors in mixed mode, a typical situation for cracks growing in the vicinity of a fretting contact. The present study is limited to straight cracks which is a simple system chosen to develop and validate the FEM analysis. The FEM model is tested and validated against popular weight functions for straight cracks perpendicular to the surface. The model is then used to study fretting crack growth and understand the effect of key parameters such as the crack angle and the friction between crack faces. Predictions achieved by this analysis match the essential features of former experimental fretting results, in particular the average crack arrest length can be predicted accurately.     

An analytical expression for the J2‐integral of an interfacial crack in orthotropic bimaterials

This paper presents a new analytical expression relating the J₂‐integral and stress intensity factors (SIF) in an in‐plane traction‐free crack between two orthotropic elastic solids using the complex function method. The singular oscillatory near tip field of a bimaterial interfacial crack is usually characterized by a pair of SIFs. In linear elastic interfacial fracture mechanics, the majority of numerical and experimental methods rely on the analytical equations relating J_k‐integrals and SIFs. Although an analytical equation relating J₁‐integral or strain energy release rate and SIFs is available, a similar relation for J₂‐integral in debonded anisotropic solids is non‐existent. Using this new analytical expression, in conjunction with the values of J_k, the SIFs can be computed without the need for an auxiliary relation. An example with known analytical solutions for SIFs is presented to show the variation of the J₂‐integral near the crack tip of a bimaterial orthotropic plate. Different bimaterial combinations are considered, and the effect of material mismatch on J_k is demonstrated.     

Elastodynamic stress intensity factors for a semi-infinite crack due to three-dimensional concentrated shear loadings on the crack faces

The dynamic stress intensity factors for a semi-infinite crack in an otherwise unbounded elastic body is investigated. The crack is subjected to a pair of suddenly-applied shear point loads on its faces at a distance l away from the crack tip. This problem is treated as the superposition of two problems. The first problem considers the disturbance by a concentrated shear force acting on the surface of an elastic half space, while the second problem discusses a half space with its surface subjected to the negative of the tangential surface displacements induced by the first problem in the front of the crack edge. A fundamental problem is proposed and solved by means of integral transforms together with the application of the Wiener–Hopf technique and Cagniard–de Hoop method. Exact expressions are then derived for the mode II and III dynamic stress intensity factors by taking integration over the fundamental solution. Some features of the solutions are discussed. This revised version was published online in July 2006 with corrections to the Cover Date.     

Estimation of high‐temperature fracture parameters for small punch specimen with a surface crack

An outline of a newly proposed methodology for evaluating creep crack growth (CCG) parameters using cracked small‐punch (SP) specimens is explained. Three‐dimensional finite element analyses were performed to calculate the stress intensity factor along the crack front for a surface crack formed at the centre of a SP specimen. Effects of crack ratio, (a/t); crack aspect ratio, (a/c); and thickness of the specimen, (t), on the fracture parameters were studied. It was observed that the minimum variation of K‐value along the crack front can be achieved when a/c was 0.50 except the location very near the intersection of the crack and free surface. This condition is similar to the case of constant K‐values along the crack front of the conventional compact tension specimen. Thus, it can be argued that the SP specimen with a surface crack is a suitable specimen geometry for CCG testing. The proposed CCG test method was found to be practically applicable for the crack geometry of 0.10 to 0.30 of a/t with constant aspect ratio of 0.50. An estimation of the K and C_t‐parameter under the small scale creep condition was derived. Future work for further development of the suggested CCG testing is discussed.     

Fractal finite element method based shape sensitivity analysis of multiple crack system

This paper presents fractal finite element based continuum shape sensitivity analysis for a multiple crack system in a homogeneous, isotropic, and two dimensional linear-elastic body subjected to mixed-mode (modes I and II) loading conditions. The salient feature of this method is that the stress intensity factors and their derivatives for the multiple crack system can be obtained efficiently since it only requires an evaluation of the same set of fractal finite element matrix equations with a different fictitious load. Three numerical examples are presented to calculate the first-order derivative of the stress intensity factors or energy release rates.     

非均匀材料界面裂纹的Cell-Based光滑有限元法

为提高非均匀材料界面裂纹尖端断裂参数的求解精度,基于非均匀材料界面断裂力学、Cell-Based光滑有限元(Cell-SFEM)和非均匀材料的互交作用积分法,提出了求解非均匀材料界面裂纹尖端断裂参数的CellBased光滑有限元法,推导了基于Cell-Based光滑有限元法的非均匀材料的互交作用积分法,对非均匀材料间的界面裂纹尖端处正则应力强度因子进行了求解,并与参考解进行了比较,讨论了互交积分区域大小和光滑子元个数与正则应力强度因子的关系。数值算例结果表明:本方法具有很高的计算精度,对积分区域大小不敏感,可为设计、制造抗破坏非均匀材料提供依据。     

Transient analysis of dynamic crack propagation in piezoelectric materials

Abstract

In this paper, the transient analysis of semi‐infinite propagating cracks in piezoelectric materials subjected to dynamic anti‐plane concentrated body force is investigated. The crack surface is assumed to be covered with an infinitesimally thin, perfectly conducting electrode that is grounded. In analyzing this problem, it has characteristic lengths and a direct attempt towards solving this problem by transform and Wiener‐Hopf techniques (Noble, 1958) is not applicable. In order to solve this problem, a new fundamental solution for propagating cracks in piezoelectric materials is first established and the transient response of the propagating crack is obtained by superposition of the fundamental solution in the Laplace transform domain. The fundamental solution to be used is the responses of applying exponentially distributed traction in the Laplace transform domain on the propagating crack surface. Taking into account the quasi‐static approximation, exact analytical transient solutions for the dynamic stress intensity factor and the dynamic electric displacement intensity factor are obtained by using the Cagniard‐de Hoop method (Cagnard, 1939; de Hoop, 1960) of Laplace inversion and are expressed in explicit forms. Numerical calculations of dynamic intensity factors are evaluated and the results are discussed in detail. The transient solutions for stationary cracks have been shown to approach the corresponding static values after the shear wave of the piezoelectric material has passed the crack tip.     

On the intersonic crack propagation in an orthotropic medium

Motivaded by recent theoretical studies the elastodynamic response of an orthotropic material with a semi-infinite line crack, which propagates intersonically. is revisited through an approach which differs from those used in previous studies. The near tip stress and displacement fields are obtained for Mode I and Mode II of steady state crack propagation. The strain energy release rate analysis confirms that the Mode I is physically impossible due to the order of stress singularity, which is larger then one half. For Model II the order of stress is less than one half and it is shown that a steady state intersonic propagation is allowed only for a particular crack tip velocity which is a function of the material orthotropy.     

Shear loaded interface crack under the influence of friction: a finite difference solution

Formulation of the elastic two‐dimensional problem of contact with friction is presented. Two‐dimensional equilibrium equations and boundary conditions in an orthogonal curvilinear co‐ordinate system are written explicitly. The above formulation is solved with the aid of the finite difference technique. An iterative algorithm which does not require load increments is employed for solving interface fracture problems with contact and friction subjected to a monotonically increasing load. The J‐integral is extended for problems in which there is friction along the crack faces. Stress intensity factors are calculated by means of the J‐integral, as well as an asymptotic expansion of the tangential shift. Two problems are analysed: (1) a crack in homogeneous material in the presence of friction involving stationary contact; and (2) an interface crack in the presence of friction involving receding contact. Results are compared to those found by analytical and semi‐analytical methods which are presented in the literature, as well as to those obtained by means of the finite element method. The accuracy of the results establishes the reliability of the finite difference analysis, as well as the post‐processors. In addition, a problem involving stick conditions is considered. It is observed that with increasing friction, the normal gaps and tangential shifts decrease. The size of the contact zone increases and values of the stress intensity factor decrease. Copyright © 2004 John Wiley & Sons, Ltd.     

Stress intensity factor analysis for an interface crack between dissimilar isotropic materials under thermal stress

Thermal stresses, one of the main causes of interfacial failure between dissimilar materials, arise from different coefficients of linear thermal expansion. Two efficient numerical procedures in conjunction with the finite element method (FEM) for the stress intensity factor (SIF) analysis of interface cracks under thermal stresses are presented. The virtual crack extension method and the crack closure integral method are modified using the superposition method. The SIF analyses of some interface crack problems under mechanical and thermal loads are demonstrated. Very accurate mode separated SIFs are obtained using these methods.     

弯扭组合载荷下圆管半椭圆表面裂纹应力强度因子的有限元分析

对表面裂纹复合型应力强度因子的研究一直是线弹性断裂力学中的重要课题,例如弯扭组合载荷下圆管半椭圆表面裂纹应力强度因子的计算,到现在也没有一个正确的分析解。考虑到裂尖的应力奇异性,在裂纹前沿手动设置三维奇异单元,用三维有限元法中的1/4点位移法计算弯扭组合载荷下圆管表面椭圆裂纹前沿的Ⅰ型、Ⅱ型和Ⅲ型应力强度因子,并分析其随裂纹深度增加时的变化规律。运用该方法计算了有关模型的应力强度因子,并与该模型的实验值进行了比较,计算结果和实验结果吻合良好。