Analysis of dynamic stress intensity factors of three-point bend specimen containing crack

A new formula is obtained to calculate dynamic stress intensity factors of the three-point bending specimen containing a single edge crack in this study. Firstly, the weight function for three-point bending specimen containing a single edge crack is derived from a general weight function form and two reference stress intensity factors, the coefficients of the weight function are given. Secondly, the history and distribution of dynamic stresses in uncracked three-point bending specimen are derived based on the vibration theory. Finally, the dynamic stress intensity factors equations for three-pointing specimen with a single edge crack subjected to impact loadings are obtained by the weight function method. The obtained formula is verified by the comparison with the numerical results of the finite element method (FEM). Good agreements have been achieved. The law of dynamic stress intensity factors of the three-point bending specimen under impact loadings varing with crack depths and loading rates is studied.     

Analysis of dynamic stress intensity factors of three-point bend specimen containing crack

A new formula is obtained to calculate dynamic stress intensity factors of the three-point bending specimen containing a single edge crack in this study. Firstly, the weight function for three-point bending specimen containing a single edge crack is derived from a general weight function form and two reference stress intensity factors, the coefficients of the weight function are given. Secondly, the history and distribution of dynamic stresses in uncracked three-point bending specimen are derived based on the vibration theory. Finally, the dynamic stress intensity factors equations for three-pointing specimen with a single edge crack subjected to impact loadings are obtained by the weight function method. The obtained formula is verified by the comparison with the numerical results of the finite element method (FEM). Good agreements have been achieved. The law of dynamic stress intensity factors of the three-point bending specimen under impact loadings varing with crack depths and loading rates is studied.     

The influence of elastic waves on dynamic stress intensity factors (two-dimensional problems)

Summary In this paper, an indirect boundary integral equation method for the solution of dynamic crack problems is presented. The Laplace transform method is used to derive the fundamental solutions for the opening mode (mode I) and the sliding mode (mode II) displacement discontinuity. Accurate dynamic stress intensity factorsK _N (t) (N=I,II) resulting from different time-dependent loads on the crack surface are obtained. The specific influences of the various elastic waves on the stress intensity factors can be clearly seen from the results.On leave Central-South University of Technology Changsha, P.R. China     

Analysis of dynamic stress intensity factors of thick-walled cylinder under internal impulsive pressure

Dynamic stress intensity factors are evaluated for thick-walled cylinder with a radial edge crack under internal impulsive pressure. Firstly, the equation for stress intensity factors under static uniform pressure is used as the reference case, and then the weight function for a thick-walled cylinder containing a radial edge crack can be worked out. Secondly, the dynamic stresses in uncracked thick-walled cylinders are solved under internal impulsive pressure by using mode shape function method. The solution consists of a quasi-static solution satisfying inhomogeneous boundary conditions and a dynamic solution satisfying homogeneous boundary condi- tions, and the history and distribution of dynamic stresses in thick-walled cylinders are derived in terms of Fourier-Bessel series. Finally, the dynamic stress intensity factor equations for thick-walled cylinder containing a radial edge crack sub- jected to internal impulsive pressure are given by dynamic weight function method. The finite element method is utilized to verify the results of numerical examples, showing the validity and feasibility of the proposed method.     

Convolution of the impact three-dimensional elasto-dynamics and dynamic stress intensity factor for an elliptic cracks

This paper presents a formulation for three-dimensional elastodynamics with an elliptic crack based on the Laplace and Fourier transforms and the convolution theorem. The dynamic stress intensity factor for the crack is determined by solving a Fredholm integral equation of the first kind. The results of this paper are very close to those given by the two-dimensional dual integral equation method. The project supported by the National Natural Science Foundation of China (K19672007)     

Dynamic stress intensity factor K III and dynamic crack propagation characteristics of anisotropic materials

Based on the mechanics of anisotropic materials, the dynamic propagation problem of a mode Ⅲ crack in an infinite anisotropic body is investigated. Stress, strain and displacement around the crack tip are expressed as an analytical complex function, which can be represented in power series. Constant coefficients of series are determined by boundary conditions. Expressions of dynamic stress intensity factors for a mode Ⅲ crack are obtained. Components of dynamic stress, dynamic strain and dynamic displacement around the crack tip are derived. Crack propagation characteristics are represented by the mechanical properties of the anisotropic materials, i.e., crack propagation velocity M and the parameter ~. The faster the crack velocity is, the greater the maximums of stress components and dynamic displacement components around the crack tip are. In particular, the parameter α affects stress and dynamic displacement around the crack tip.     

压电界面裂纹裂尖强度因子的显式计算公式

对常见横观各向同性压电材料(TIP)中界面裂纹的裂纹面与压电材料的极化方向成任意夹角的一般情况进行了研究,通过推导得到了计算裂尖强度因子的显式外推公式,同时给出了裂纹面与极化方向垂直的典型情况下的外推公式.这些显式计算公式为常见数值方法如有限元法及边界元法在压电材料断裂力学中的应用提供了便利.     

热、机械载荷作用下夹杂对应力强度因子的影响

推导了远场应力、热应力耦合作用下含夹杂裂纹体的应力强度因子求解公式,改进了体积力法中的裂纹面合力平衡条件,将应力强度因子的求解归结为解一组积分方程,再将积分方程转化为线性方程组进行数值求解。算例分析结果表明方法正确、有效。将该算法应用于含Al2O3夹杂的FGH95材料应力强度因子分析中,计算结果表明热应力对应力强度因子影响很小。     

Determination of the dynamic stress intensity factorsK I d andK II d,for a mixed-mode propagating crack

In this paper, the dynamic propagation problem of a mixed-mode crack was studied by means of the experimental method of caustics. The initial curve and caustic equations were derived under the mixed-mode dynamic condition. A multi-point measurement method for determining the dynamic stress intensity factors,K _I ^d , andK _II ^d , and the position of the crack tip was developed. Several other methods were adopted to check this method, and showed that it has a good precision. Finally, the dynamic propagating process of a mixed-mode crack in the three-point bending beam specimen was investigated with our method.     

动、静裂纹作用偏置效应的动焦散冲击实验

为了研究冲击荷载作用下脆性材料中运动裂纹与静止裂纹的相互作用,选取动态载荷下断裂行为与岩石材料类似且本身光学特性较好的有机玻璃（PMMA）作为实验材料,试件尺寸为220 mm×50 mm×5 mm,采用激光切割制作长度5 mm的预制裂纹和长度10 mm的静止裂纹,预制裂纹位于试件的底部边缘中心,静止裂纹的中心位于试件水平轴线。将静止裂纹偏置距离作为单一变量,采用数字激光动态焦散实验系统对含不同缺陷的PMMA进行三点弯曲实验,并结合几何分形理论研究不同偏置距离下运动裂纹的分形规律。实验结果表明：存在预制裂纹与静止裂纹的临界偏置距离（6 mm）,该条件下裂纹轨迹对应的分形维数值最大,裂纹轨迹的规则程度最低,裂纹破坏形态最复杂。当预制裂纹与静止裂纹的偏置距离在0～6 mm时,裂纹Ⅰ起裂后垂直向上扩展一段距离,与静止裂纹交汇,并停滞一段时间后发生二次起裂,直至贯穿试件,偏置距离和交汇点竖向坐标值呈近似线性函数关系。偏置距离的存在不会影响裂纹Ⅰ的起裂时间和应力强度因子,但会显著减小裂纹Ⅱ的动态应力强度因子,且停滞时长随偏置距离的增大而逐渐缩短。当偏置距离大于临界偏置距离时,运动裂纹不再与静止裂纹交汇而是呈拱状向试件上边缘扩展直至贯穿,裂纹的起偏时间、起偏位置也会出现明显的滞后现象。

     

Determining stress intensity factors in orthotropic composites from far-field measured temperatures

We demonstrate the ability to determine stress intensity factors in orthotropic materials directly from measured temperatures away from the crack and using far-field expressions for the stresses. This is advantageous, recognizing that recorded thermoelastic data can be very unreliable near the tip of a crack. In addition to singular terms that govern in the immediate vicinity of the crack tip, the present series expressions for the stresses contain higher-order finite terms. Little measured input information is needed and data acquisition positions can be selected largely at the user's discretion.     

一种确定结合材料界面端应力强度因子的数值方法

基于弹性力平面问题的基本方程,给出了结合材料界面端的应力奇异性特征方程以及位移场和奇异应力场。提出了一种确定结合材料界面端应力强度因子的数值外插方法。对界面端区域进行了有限元网格单元划分。经过具体实例检验进一步确定了求解应力强度因子的最佳方向,该数值外插法的计算结果精度符合工程应用的要求,为工程材料强度的评价提供了有效的计算途径。     

Numerical simulation of stress intensity factors via the thermoelastic technique

The purpose of this study is to investigate the accuracy of the least squares method for finding the in-plane stress intensity factorsK _I andK _II using thermoelastic data from isotropic materials. To fully understand the idealized condition ofK _I andK _II calculated from thermoelastic experiments, the total stress field calculated from finite element analysis is used to take the place of data obtained from real thermoelastic experiments. In the finite element analysis, theJ-integral is also calculated to compare with (K _I ² +K _II ² )/E evaluated by the least squares method. The stress fields near the crack tip are dominated by the two stress intensity factors; however, the edge effect will cause inaccuracy of the thermoelastic data near the crack tip. Furthermore, the scan area of thermoelastic experiments cannot be too small. Therefore, we suggest that three or four terms of stress function be included in the least squares method for evaluating stress intensity factors via the thermoelastic technique. In the idealized condition, the error can be smaller than 3 percent from our numerical simulations. If only ther ^–1/2 term (K _I andK _II ) is included in the least squares method, even in the idealized case the error can be up to 20 percent.     

Dynamic stress intensity factors of multiple cracks in a functionally graded orthotropic half-plane

In the present paper dynamic stress intensity factor and strain energy density factor of multiple cracks in the functionally graded orthotropic half-plane under time-harmonic loading are investigated. By utilizing the Fourier transformation technique the stress fields are obtained for a functionally graded orthotropic half-plane containing a Volterra screw dislocation. The variations of the material properties are assumed to be exponential forms which the equilibrium has an analytical solution. The dislocation solution is utilized to formulate integral equation for the half-plane weakened by multiple smooth cracks under anti-plane deformation. The integral equations are of Cauchy singular type at the location of dislocation which are solved numerically to obtain the dislocation density on the faces of the cracks. The dislocation densities are employed to determined stress intensity factor and strain energy density factors (SEDFs) for multiple smooth cracks under anti-plane deformation. Numerical examples are provided to show the effects of material properties and the crack configuration on the dynamic stress intensity factors and SEDFs of the functionally graded orthotropic half-plane with multiple curved cracks.     

Computing lower and upper bounds on stress intensity factors in bimaterials

This paper presents a finite element approach for finding complementary bounds of stress intensity factors (SIFs) in bimaterials. The SIF is formulated as an explicit computable linear function of displacements by means of the two-point extrapolation method. An appropriate and computable form of the SIF plays a crucial role in the dual problem involved in the computing procedure of both lower and upper bounds. In our discussions, computable forms of stress intensity factors, K₀ and K_r, are derived, which are related to the energy release rate, and the ratio of the open mode and shear mode SIFs, respectively. Based on a posteriori finite element error estimation, a bounding procedure is used to compute the bounds on the two stress intensity factors. Finally, bounds on the SIFs in a bimaterial interface crack problem are provided to verify the procedure.     

Mode III stress intensity factors in an interfacial crack in dissimilar bonded materials

The antiplane shear deformation problem of two edge-bonded dissimilar isotropic wedges is considered. In the case when the sum of the two apex angles is equal to 2π, the problem reduces to that of two edge-bonded dissimilar materials with an interfacial crack subjected to concentrated antiplane shear tractions on the crack faces. An explicit expression is extracted for the stress intensity factor at the crack tip. In the special cases of different combinations of the apex angles, the obtained expression for the stress intensity factor may be simplified and relations of a simpler form are given for the stress intensity factor. It is shown that the stress intensity factor is dependent on the material properties as well as the geometry and loading. However, in special cases of equal apex angles as well as the case of similar materials the dependency of the stress intensity factor on the material properties disappears.     

Residual stress field and reduction of stress intensity factors in cold-worked holes

Closed-form and semi-analytical solutions are obtained for the residual stress distributions in a plate caused by pressure acting on a central circular hole, representing the cold-work process. The material is elastic–perfectly plastic. Both Tresca and von Mises yield criteria are used and the corresponding residual stress distributions are compared. The relation between the dimension of the plastic zone and the value of internal pressure is presented. The relation between the magnitude of the residual stresses and the remote uniform tensile stress required to open symmetrical radial cracks is also presented. The reduction of the stress intensity factors of cracked open and riveted holes as a function of the internal pressure applied (or mandrel radial displacement) is investigated using numerical models for both an elastic–perfectly plastic material and for an Al 2024-T3 Alclad aluminum alloy.     

3D dynamic stress intensity factors at three rectangular cracks in an infinite elastic medium subjected to a time-harmonic stress wave

Summary Dynamic stresses around three coplanar cracks in an infinite elastic medium are determined in the paper. Two of the cracks are equal, rectangular and symmetrically situated on either side of the centrally located rectangular crack. Time-harmonic normal traction acts on each surface of the three cracks. To solve the problem, two kind of solutions are superposed: one is a solution for a rectangular crack in an infinite elastic medium, and the other one is that for two rectangular cracks in an infinite elastic medium. The unknown coefficients in the combined solution are determined by applying the boundary conditions at the surfaces of the cracks. Finally, stress intensity factors are calculated numerically for several crack configurations. Received 14 July 1998; accepted for publication 2 December 1998     

数值积分法求解厚壁筒的动态应力强度因子

用数值积分法求解了厚壁筒表面裂纹的动态应力强度因子，其结果与有限元的计算结果作了比较，表明该方法简单有效，对工程应用极有参考价值     

A fractional differential constitutive model for dynamic stress intensity factors of an anti-plane crack in viscoelastic materials

Fractional differential constitutive relationships are introduced to depict the history of dynamic stress inten- sity factors （DSIFs） for a semi-infinite crack in infinite viscoelastic material subjected to anti-plane shear impact load. The basic equations which govern the anti-plane deformation behavior are converted to a fractional wave-like equation. By utilizing Laplace and Fourier integral transforms, the fractional wave-like equation is cast into an ordinary differential equation （ODE）. The unknown function in the solution of ODE is obtained by applying Fourier transform directly to the boundary conditions of fractional wave-like equation in Laplace domain instead of solving dual integral equations. Analytical solutions of DSIFs in Laplace domain are derived by Wiener-Hopf technique and the numerical solutions of DSIFs in time domain are obtained by Talbot algorithm. The effects of four parameters α, β, b1, b2 of the fractional dif- ferential constitutive model on DSIFs are discussed. The numerical results show that the present fractional differential constitutive model can well describe the behavior of DSIFs of anti-plane fracture in viscoelastic materials, and the model is also compatible with solutions of DSIFs of anti-plane fracture in elastic materials.