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

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

Stress intensity factors of an oblique edge crack subjected to normal and shear tractions

Stress intensity factors (SIFs) were obtained for an oblique crack under normal and shear traction and remote extension loads. The oblique crack was modeled as the pseudodislocation. The stress field due to tractions was solved by the Flamant solution. The SIR of Mode I and Mode II (K_Iand K_II) were then obtained. Finite element analysis was performed with ABAQUS and compared with the analytical solutions. The analytical solutions were in good agreement with the results of FEM. From investigating SIFs and their ranges, the following results were obtained. The growth rate of an oblique edge crack decreased due to the reduction in the SIF ranges. The crack driving force depended on the obliquity, the normal traction and the ratio of crack to traction length. The peak value of shear traction was found as a key parameter to accelerate the crack growth.     

Stress-intensity factors for materials with interface cracks under harmonic loading

The paper studies the distribution of the mode I, II, and III stress-intensity factors (SIFs) for a penny-shaped interface crack in a bimaterial under normal harmonic loading. The dependence of SIFs on the type of a wave (compression–dilatation or shear), loading frequency, and the mechanical properties of the bimaterial is analyzed     

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.     

The mixed-mode fracture mechanics analysis of an embedded arbitrary oriented crack in a two-dimensional functionally graded material plate

Mixed-mode fracture mechanics analysis of an embedded arbitrarily oriented crack in a two-dimensional functionally graded material using plane elasticity theory is considered. The material properties are assumed to vary exponentially in two planar directions. Then, employing Fourier integral transforms with singular integral equation technique, the problem is solved. The stress intensity factors (SIFs) at the crack tips are calculated under in-plane mechanical loads. Finally, the effects of crack orientation, material non-homogeneity, and other parameters are discussed on the value of SIF in mode I and mode II fracture.     

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     

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.     

MODE III 2-D FRACTURE ANALYSIS BY THE SCALED BOUNDARY FINITE ELEMENT METHOD

The scaled boundary finite element method （SBFEM） is a novel semi-analytical technique that combines the advantages of the finite element method and the boundary element method with unique properties of its own. This method has proven very efficient and accurate for determining the stress intensity factors （SIFs） for mode I and mode II two-dimensional crack problems. One main reason is that the SBFEM has a unique capacity of analytically representing the stress singularities at the crack tip. In this paper the SBFEM is developed for mode III （out of plane deformation） two-dimensional fracture anMysis. In addition, cubic B-spline functions are employed in this paper for constructing the shape functions in the circumferential direction so that higher continuity between elements is obtained. Numerical examples are presented at the end to demonstrate the simplicity and accuracy of the present approach for mode Ⅲ two-dimensional fracture analysis.     

A numerical approach for a crack problem by Gauss–Chebyshev quadrature

In this paper, a problem of a crack in an orthotropic strip is studied under plane strain conditions. It is assumed that normal displacements and shear stresses do not act on neither of the boundaries of the strip. Cauchy-type singular integral equation for the crack problem is derived by using the theory of plane elasticity and the Fourier transformation technique. A quadrature collocation approach is adopted for the numerical solutions of the singular integral equation. The effect of relative thickness and mechanical properties of strip on Mode I stress intensity factors (SIFs) are examined under different loading conditions. Some sample results are given for SIFs; also, material orthotropy and geometrical effects are discussed in detail.     

求解变分型积分方程的一种新型数值方法——有限变分法

将作者提出的多虚拟裂纹扩展法(MVCE法)拓展为求解变分型积分方程问题的一种新型数值方法——有限变分法(FVM)。它的基本思想是,给定有限个(N个)局部变分模式,将所求解的未知量用适当的方法离散化,针对这N个局部变分模式列出N个方程,求解N个未知系数,从而求得未知量。单一未知变量FVM的最终方程组的系数矩阵通常是一个对称的窄带矩阵,对角元是大数,有很好的数值计算性能。用FVM求解了三维I型裂纹前缘的应力强度因子(SIF)分布。利用基于FVM的通用权函数法计算程序,可以高精度和高效率地求解表面力、体积力和温度载荷共同作用情况下三维裂纹前缘SIF的分布及其时间历程。FVM可以被推广到更广泛的领域,是一个求解变分型积分方程问题的普遍适用的新型数值方法。     

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

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

Hypersingular integral equation method for a three-dimensional crack in anisotropic electro-magneto-elastic bimaterials

Using Green’s functions, the extended general displacement solutions of a three-dimensional crack problem in anisotropic electro-magneto-elastic (EME) bimaterials under extended loads are analyzed by the boundary element method. Then, the crack problem is reduced to solving a set of hypersingular integral equations (HIE) coupled with boundary integral equations. The singularity of the extended displacement discontinuities around the crack front terminating at the interface is analyzed by the main-part analysis method of HIE, and the exact analytical solutions of the extended singular stresses and extended stress intensity factors (SIFs) near the crack front in anisotropic EME bimaterials are given. Also, the numerical method of the HIE for a rectangular crack subjected to extended loads is put forward with the extended crack opening dislocation approximated by the product of basic density functions and polynomials. At last, numerical solutions of the extended SIFs of some examples are obtained.     

Dislocation emission criterion from a blunt crack tip

Simple closed-form solutions for critical stress intensity factors were developed for edge dislocation emission from a blunt crack under mode I and mode II loading conditions. The solutions are applicable both for a slender elliptical hole and a long slit with blunted head. When the blunt crack considered degenerates to a sharp crack, the present solutions reduce to those previously reported in literature.     

星形裂纹的应力分析

利用复变函数的方法, 通过构造适当的保角映射研究了星形裂纹的平面弹性问题,给出了裂纹尖端I型与II型问题应力强度因子的解析解.并由此模拟出了经典的Griffith裂纹,共点均匀分布三裂纹,十字裂纹,对称八裂纹问题.      

Mode i and mixed mode fields with incomplete crack tip plasticity

Plane strain slip line fields, in which plasticity does not fully surround the crack tip have been developed for mode I and mixed mode I\II cracks under contained yielding. Analytical solutions have been assembled using slip line theory for the plastic sectors and semi-infinite wedge solutions for the elastic sectors. These solutions are compared with finite element solutions based on modified boundary layer formulations. The analytical solutions agree well with numerical solutions, and form a family of fields with incomplete plasticity around the crack tip.     

On partial sliding along a planar crack: the case of a circular sliding zone

Summary A problem of partial sliding along a planar crack with a local drop in frictional resistance is investigated. A sliding zone initiates in the area of reduced friction, and then propagates as the applied shear load is monotonously increased. The problem is formulated in general terms, and then solved for the case when sliding spreads as a penny-shaped zone. Conditions under which the front of the zone stays circular during sliding are analyzed. It is observed that the axisymmetry of the profile of frictional resistance does not necessarily guarantee uniform propagation of sliding in the radial direction. The circular shape becomes the most favorable growth condition only if the shear modes are related in a certain way. The problem is studied based on the criterion of propagation that stress intensity factors(SIFs) for II and III modes vanish on the boundary of the sliding zone. The singular integrals in expressions for the SIFs are reduced to non-singular ones. Analytical solutions are derived for a number of special cases where the radius of the sliding zone is related to the applied shear load.This work is supported by the National Science Foundation through grant DGE-0209543 to the University of New Mexico.     

基于Erdogan基本解边界元法计算应力强度因子

引入含裂纹问题基本解(Erdogan基本解),提出了基于Erdogan基本解的样条虚边界 元法,并阐述了该法在实施过程中的特点与具体做法. 采用该方法详细分析了若干 典型裂纹问题,全面考察了方法的计算精度和收敛情况,以及在求解复杂裂纹问题方面 的能力. 结果显示,该方法具有精度高、收敛快、计算能力强等优点,是裂纹问题分析中 一种具有竞争力的通用计算方法.     

Stress intensity factor measurement of cracks using a piezoelectric element

We present a stress intensity factor (SIF) measurement method of cracks using a piezoelectric element and electrostatic voltmeter. In the method, an isotropic piezoelectric element is first adhered near the crack tip. Then, the surface electrodes are attached to the three different positions on the piezoelectric element. The electric potentials of the surface electrodes, which are proportional to the strain sum (ɛ_x+ɛ_y) on the structural member, are measured by an electrostatic voltmeter during load cycling. Mode I and mode II SIFs of the crack are estimated using the relationship between the SIF and (σ_x+σ_y). The applicability of the proposed method is examined through experiments and numerical analysis.     

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.     

A new finite element method for strain gradient theories and applications to fracture analyses

A new compatible finite element method for strain gradient theories is presented. In the new finite element method, pure displacement derivatives are taken as the fundamental variables. The new numerical method is successfully used to analyze the simple strain gradient problems – the fundamental fracture problems. Through comparing the numerical solutions with the existed exact solutions, the effectiveness of the new finite element method is tested and confirmed. Additionally, an application of the Zienkiewicz–Taylor C1 finite element method to the strain gradient problem is discussed. By using the new finite element method, plane-strain mode I and mode II crack tip fields are calculated based on a constitutive law which is a simple generalization of the conventional J₂ deformation plasticity theory to include strain gradient effects. Three new constitutive parameters enter to characterize the scale over which strain gradient effects become important. During the analysis the general compressible version of Fleck–Hutchinson strain gradient plasticity is adopted. Crack tip solutions, the traction distributions along the plane ahead of the crack tip are calculated. The solutions display the considerable elevation of traction within the zone near the crack tip.