Determination of stress intensity factor with direct stress approach using finite element analysis

In this article,a direct stress approach based on finite element analysis to determine the stress intensity fac-tor is improved.Firstly,by comparing the rigorous solution against the asymptotic solution for a problem of an infinite plate embedded a central crack,we found that the stresses in a restrictive interval near the crack tip given by the rigorous solution can be used to determine the stress intensity fac-tor,which is nearly equal to the stress intensity factor given by the asymptotic solution.Secondly,the crack problem is solved numerically by the finite element method.Depending on the modeling capability of the software,we designed an adaptive mesh model to simulate the stress singularity.Thus, the stress result in an appropriate interval near the crack tip is fairly approximated to the rigorous solution of the corre-sponding crack problem.Therefore,the stress intensity factor may be calculated from the stress distribution in the appro-priate interval,with a high accuracy.     

拉伸载荷作用下单边缺陷-边裂纹应力强度因子研究

利用杂交位移不连续法研究拉伸载荷作用下矩形板中单边缺陷-边裂纹(半圆孔裂纹和半方孔裂纹)问题,给出了这三种平面弹性裂纹问题的应力强度因子的详细数值解。通过半圆孔裂纹问题和半方孔裂纹问题与单边裂纹问题的应力强度因子的比较,发现半圆孔和半方孔对单边裂纹有屏蔽影响。此外,本文的研究结果表明,杂交位移不连续法用于分析平面弹性有限体中复杂裂纹问题的应力强度因子简单且又准确。     

焊接节点半椭圆表面裂纹应力强度因子的权函数计算

本文采用一种改进的权函数法来计算焊接节点半椭圆表面裂纹应力强度因子ＫＩ值，并给出了相应的数值处理方法，就Ｔ型板节点进行了数值验算。     

基于子模型法的带有表面裂纹钢丝应力强度因子研究

钢丝裂纹应力强度因子是进行钢丝疲劳断裂寿命评估、疲劳裂纹扩展分析和钢丝断裂强度评估等工作的重要参数。本文首先介绍了裂纹扩展分析软件FRANC3D,然后基于子模型法模拟研究了拉伸荷载作用下带有表面裂纹钢丝的应力强度因子,裂纹种类包括直线形裂纹和半圆形裂纹,最后拟合得到拉伸荷载作用下带表面裂纹钢丝的应力强度因子形状修正系数表达式,分析了利用该公式进行承载力评估时产生误差的原因。研究结果表明,利用子模型模拟分析拉伸荷载作用下带有表面裂纹的钢丝应力强度因子时计算精度高,计算速度快,对计算机硬件要求低;利用该方法得到的钢丝裂纹应力强度因子,在进行索承式桥梁吊索安全性能评估时,评估结果更精确。     

Micromechanics of particle-reinforced elasto-viscoplastic composites: Finite element simulations versus affine homogenization

The micromechanics of elasto-viscoplastic composites made up of a random and homogeneous dispersion of spherical inclusions in a continuous matrix was studied with two methods. The first one is an affine homogenization approach, which transforms the local constitutive laws into fictitious linear thermo-elastic relations in the Laplace–Carson domain so that corresponding homogenization schemes can apply, and the temporal response is computed after a numerical inversion of Laplace transform. The second method is the direct numerical simulation by finite elements of a three-dimensional representative volume element of the composite microstructure. The numerical simulations carried out over different realizations of the composite microstructure showed very little scatter and thus provided – for the first time – “exact” results in the elasto-viscoplastic regime that can be used as benchmarks to check the accuracy of other models. Overall, the predictions of the affine homogenization model were excellent, regardless of the volume fraction of spheres, of the loading paths (shear, uniaxial tension and biaxial tension as well as monotonic and cyclic deformation), particularly at low strain rates. It was found, however, that the accuracy decreased systematically as the strain rate increased. The detailed information of the stress and strain microfields given by the finite element simulations was used to analyze the source of this difference, so that better homogenization methods can be developed.     

Analysis of stress intensity factor in orthotropic bi-material mixed interface crack

Adopting the complex function approach, the paper studies the stress intensity factor in orthotropic bi-material interface cracks under mixed loads. With consideration of the boundary conditions, a new stress function is introduced to transform the problem of bi-material interface crack into a boundary value problem of partial differential equations. Two sets of non-homogeneous linear equations with 16 unknowns are constructed. By solving the equations, the expressions for the real bi-material elastic constant εt and the real stress singularity exponents λt are obtained with the bi-material engineering parameters satisfying certain conditions. By the uniqueness theorem of limit,undetermined coefficients are determined, and thus the bi-material stress intensity factor in mixed cracks is obtained. The bi-material stress intensity factor characterizes features of mixed cracks. When orthotropic bi-materials are of the same material, the degenerate solution to the stress intensity factor in mixed bi-material interface cracks is in complete agreement with the present classic conclusion. The relationship between the bi-material stress intensity factor and the ratio of bi-material shear modulus and the relationship between the bi-material stress intensity factor and the ratio of bi-material Young's modulus are given in the numerical analysis.     

Asymptotic fields in the finite element analysis of electrically permeable interface cracks in piezoelectric bimaterials

Summary The interface crack problem for a piezoelectric bimaterial based on permeable conditions is studied numerically. To find the singular electromechanical field at the crack tip, an asymptotic solution is derived in connection with the conventional finite element method. For mechanical and electrical loads, the complex stress intensity factor for an interface crack is obtained. The influence of the applied loads on the electromechanical fields near the crack tip is also studied. For a particular case of a short crack with respect to the bimaterial size, the numerical results are compared with the exact analytical solutions, obtained for a piezoelectric bimaterial plane with an interface crack.One author (V.G.) gratefully acknowledges the support provided by the Alexander von Humboldt Foundation of Germany.accepted for publication 7 June 2004     

A different approach for calculation of stress intensity factors in continuous fiber reinforced metal matrix composites

In this paper, a numerical model developed for the analysis of a cylindrical element of matrix containing a single fiber is presented. A ring-shaped crack is assumed at interface of fiber and matrix. Both layers in the model are bonded perfectly with the exception of the crack faces. Contact elements, which have bonded feature, are used between fiber and matrix. Displacement correlation method is used to calculate opening-mode and sliding-mode stress intensity factors. These results obtained from the analysis help to understand the debonding phenomenon between fiber and matrix interface. Effects of the mechanical properties of fiber and matrix on direction of crack propagation are also discussed.     

空心车轴不同深度表面裂纹应力强度因子研究

分析了空心车轴的旋转弯曲载荷的特点,建立了空心轴表面周向半椭圆裂纹的模型,给出了半椭圆裂纹的构形参数定义,即形状比、深度比和裂纹前缘相对位置。采用四分之一20节点等参退化奇异单元,通过有限元计算,模拟裂纹前沿的应力奇异性。在此基础上,计算了裂纹前缘表面点和中心点的应力强度因子随着裂纹扩展深度和旋转角度的变化。计算结果表明,对于给定的裂纹构形,在车轴的一个载荷循环中,裂纹前缘同一相对位置的应力强度因子是不断变化的,不同位置的应力强度因子在达到最大值的角度也是不同的,这就导致了裂纹前缘表面点在一些角度下的扩展是不对称的。这些结果为进一步研究空心轴表面裂纹的扩展路径和寿命提供了参考。     

Effect of electric displacement saturation on the stress intensity factor for a crack in a ferroelectric ceramic

A crack in a ferroelectric ceramic with perfect saturation under electric loading is analyzed. The boundary of the electric displacement saturation zone ahead of the crack tip is assumed to be ellipse in shape. The shape and size of ferroelectric domain switching zone near a crack tip is determined based on the nonlinear electric theory. The stress intensity factor induced by ferroelectric domain switching under small-scale conditions is numerically obtained as a function of the electric saturation zone parameter and the ratio of the coercive electric field to the yield electric field. It is found that the stress intensity factor increases as the ratio of the semi-axes of the saturation ellipse increases.     

Anti-plane stress intensity, energy release and energy density at crack tips in a functionally graded strip with linearly varying properties

The fracture problem of a crack in a functionally graded strip with its properties varying in a linear form along the strip thickness under an anti-plane load is considered. The embedded anti-plane crack is located in the middle of strip half way through the thickness. The third mode stress intensity factor is derived using two different methods. In the first method, by employing Fourier integral transforms, the governing equation is converted to a singular integral equation, which is subsequently solved numerically by the collocation method based on Chebyshev polynomials. Then, the problem is solved by means of finite element method in which quadrilateral 8-node singular elements around each crack tip are used. After inspecting the validity of the solution technique, effects of crack geometry and non-homogeneous material parameter on the stress intensity, energy release and energy density are studied and the results of analytical and FEM solutions are compared.     

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.     

求解界面裂纹应力强度因子的围线积分法

本文基于Ｂｅｔｔｉ功互等定理和双材料界面裂纹辅助场，提出了一种求解界面裂纹应力强度因子的方法，即远场围线积分法。此方法与积分径的选择无关，用有元元法计算出远离裂纹尖端的位移场和应力场，应可通过计算绕裂尖围线的积分，精确地给出界面裂纹应力强度因子ＫＩ和ＫⅡ。     

Stress intensity factors for interface crack in finite size specimen using a generalized variational approach

A generalized variational approach together with eigenfunction expansion is applied to determine the stress intensity factors for interface crack in finite size specimen. Application is also made of the complex potentials such that a complex stress intensity factor with components corresponding to the Mode I and II stress intensity factors can be identified with one of the leading coefficients in the eigenfunction expansion. Obtained are the numerical values of the stress intensity factors for an interface edge crack in a bimaterial rectangular specimen. The outside boundary is subjected to uniform stress normal and parallel to the crack. Solutions are also obtained for the same crack aand specimen geoinetry is subjected to a pair of equal and opposite concentrated forces along the open end away from the edge crack. The third example pertains to the case of three-point bending where the centre concentrated load is directed along the interface dividing the two materials. Numerical results are obtained for four different combinations of the bimaterial specimen with an interface edge crack.     

Computing bounds to mixed-mode stress intensity factors in elasticity

A bounding procedure combined with an effective error bound method for linear functionals of the displacements and a simple two points displacement extrapolation method is presented to compute the lower and upper bounds to the stress intensity factors in elastic fracture problems. First, the displacements of two nodes (or node pairs) on the crack edges are used to construct the linear extrapolation to obtain the stress intensity factors at the crack tip, so that stress intensity factors are explicitly expressed as linear functionals of the displacements. Then, a posteriori bounding method is utilized to compute the bounds to the stress intensity factors. Finally, the bounding procedure is verified by a mixed-mode homogenous elastic fracture problem and a bimaterial interface crack problem.     

Variation of the stress intensity factor along the front of a 3-D rectangular crack subjected to mixed-mode load

Summary  The singular integral equation method is applied to the calculation of the stress intensity factor at the front of a rectangular crack subjected to mixed-mode load. The stress field induced by a body force doublet is used as a fundamental solution. The problem is formulated as a system of integral equations with r ⁻³-singularities. In solving the integral equations, unknown functions of body-force densities are approximated by the product of polynomial and fundamental densities. The fundamental densities are chosen to express two-dimensional cracks in an infinite body for the limiting cases of the aspect ratio of the rectangle. The present method yields rapidly converging numerical results and satisfies boundary conditions all over the crack boundary. A smooth distribution of the stress intensity factor along the crack front is presented for various crack shapes and different Poisson's ratio. Received 5 March 2002; accepted for publication 2 July 2002     

残余应力下厚壁筒表面裂纹的应力强度因子计算

本文首先介绍了边界元法计算裂纹尖端应力强度因子的基本理论,接着利用边界元法计算了在残余应力下不同厚壁筒内表面椭圆裂纹的应力强度因子,研究了其大不随椭圆裂纹不同而变化的规律,为厚壁筒结构的设计,制造以及疲劳寿命分析提供了许多有价值的参考资料。     

Determination of stress intensity factors by the finite element discretized symplectic method

When rewriting the governing equations in Hamiltonian form, analytical solutions in the form of symplectic series can be obtained by the method of separation of variable satisfying the crack face conditions. In theory, there exists sufficient number of coefficients of the symplectic series to satisfy any outer boundary conditions. In practice, the matrix relating the coefficients to the outer boundary conditions is ill-conditioned unless the boundary is very simple, e.g., circular. In this paper, a new two-level finite element method using the symplectic series as global functions while using the conventional finite element shape functions as local functions is developed. With the available classical finite elements and symplectic series, the main unknowns are no longer the nodal displacements but are the coefficients of the symplectic series. Since the first few coefficients are the stress intensity factors, post-processing is not required. A number of numerical examples as well as convergence studies are given.     

三维裂纹体应力强度因子的计算

本文采用塌缩三棱柱形奇异单元的位移计算应力强度因子，给出了一个新的全三维外推公式，它是Ｃｈｅｎ和Ｋｕａｎｇ公式〔１３〕的全三维推广，特例证明，它的精度比Ｉｎ－ｇｒａｆｆｅａ和Ｍａｎｕ的公式〔８〕高一阶。数值计算表明，结果稳定和对单元尺寸改变不敏感     

The 3-D stress intensity factor for a half plane crack in a transversely isotropic solid due to the motion of loads on the crack faces

Three-dimensional analysis of a half plane crack in a transversely isotropic solid is performed. The crack is subjected to a pair of normal point loads moving in a direction perpendicular to the crack edge on its faces. Transform methods are used to reduce the boundary value problem to a single integral equation that can be solved by the Wiener-Hopf technique. The Cagniard-de Hoop method is employed to invert the transforms. An exact expression is derived for the mode I stress intensity factor as a function of time and position along the crack edge. Some features of the solution are discussed through numerical results. The project supported by the Guangdong Provincial Natural Science Foundation and the Science Foundation of Shantou University