Stress intensity factors calculation in anti-plane fracture problem by orthogonal integral extraction method based on femol

For an anti-plane problem, the differential operator is self-adjoint and the corresponding eigenfunctions belong to the Hilbert space. The orthogonal property between eigenfunctions (or between the derivatives of eigenfunctions) of anti-plane problem is exploited. We developed for the first time two sets of radius-independent orthogonal integrals for extraction of stress intensity factors (SIFs), so any order SIF can be extracted based on a certain known solution of displacement (an analytic result or a numerical result). Many numerical examples based on the finite element method of lines (FEMOL) show that the present method is very powerful and efficient.

Unloading characteristics of anti-plane shear crack in softening materials

The local characteristics of the anti-plane shear stress and strain field are determined for a material where the stress increases linearly with strain up to a limit and then softens nonlinearly. Two unloading models are considered such that the unloading path always returns to the origin while the other assumes the unloading modulus to be that of the initial shear modulus. As the applied shear increases, an unloading zone is found to prevail between a zone in which the material softens and another zone in which the material is linear-elastic even though the crack does not propagate. The divisions of these zones are displayed graphically.

The relationship of SIF between plate and plane fracture problems and the effect of the plate thickness on SIF

In the present paper, it is shown that the zero series eigenfunctions of Reissner plate cracks/notches fracture problems are analogous to the eigenfunctions of anti-plane and in-plane. The singularity in the double series expression of plate problems only arises in zero series parts. In view of the relationship with eigen-values of anti-plane and in-plane problem, the solution of eigen-values for Reissner plates consists of two parts: anti-plane problem and in-plane problem. As a result the corresponding eigen-values or the corresponding eigen-value solving programs with respect to the anti-plane and in-plane problems can be employed and many aggressive SIF computed methods of plane problems can be employed in the plate. Based on those, the approximate relationship of SIFs between the plate and the plane fracture problems is figured out, and the effect relationship of the plate thickness on SIF is given.

Computation of stress intensity factors by the sub-region mixed finite element method of lines

Based on the sub-region generalized variational principle, a sub-region mixed version of the newly-developed semi-analytical 'finite element method of lines' (FEMOL) is proposed in this paper for accurate and efficient computation of stress intensity factors (SIFs) of two-dimensional notches/cracks. The circular regions surrounding notch/crack tips are taken as the complementary energy region in which a number of leading terms of singular solutions for stresses are used, with the sought SIFs being among the unknown coefficients. The rest of the arbitrary domain is taken as the potential energy region in which FEMOL is applied to obtain approximate displacements. A mixed system of ordinary differential equations (ODEs) and algebraic equations is derived via the sub-region generalized variational principle. A singularity removal technique that eliminates the stress parameters from the mixed equation system eventually yields a standard FEMOL ODE system, the solution of which is no longer singular and is simply and efficiently obtained using a standard general-purpose ODE solver. A number of numerical examples, including bi-material notches/cracks in anti-plane and plane elasticity, are given to show the generally excellent performance of the proposed method.

Dynamic SIF of interface crack between two dissimilar viscoelastic bodies under impact loading

In this paper, the transient dynamic stress intensity factor (SIF) is determined for an interface crack between two dissimilar half-infinite isotropic viscoelastic bodies under impact loading. An anti-plane step loading is assumed to act suddenly on the surface of interface crack of finite length. The stress field incurred near the crack tip is analyzed. The integral transformation method and singular integral equation approach are used to get the solution. By virtue of the integral transformation method, the viscoelastic mixed boundary problem is reduced to a set of dual integral equations of crack open displacement function in the transformation domain. The dual integral equations can be further transformed into the first kind of Cauchy-type singular integral equation (SIE) by introduction of crack dislocation density function. A piecewise continuous function approach is adopted to get the numerical solution of SIE. Finally, numerical inverse integral transformation is performed and the dynamic SIF in transformation domain is recovered to that in time domain. The dynamic SIF during a small time-interval is evaluated, and the effects of the viscoelastic material parameters on dynamic SIF are analyzed.

导电薄板内电流密度分布与反平面剪切的比拟

定量分析电流密度在含裂纹载流薄板内的分布是当前利用电流热效应止裂技术中一个首先要解决的问题.由于裂纹的存在,电流密度在裂尖形成带奇异性分布的高度密集.现有的分析方法往往比较复杂或局限于特殊布置形式的裂纹.通过电流密度分布与弹性力学里反平面剪切问题的比拟,把分析含裂纹载流薄板内电流密度的分布等效于考虑相应的Ⅲ型裂纹问题,并比照Ⅲ型裂纹的应力强度因子来定义电流密度因子.而对于裂纹问题的处理可采用分布位错法这一断裂力学里便利有效的分析手段.由给出的算例可见,所提出的比拟解法可以方便精确地求解电流密度在裂尖附近的奇异分布,并有助于对这一奇异性在概念上的直观理解.

Asymptotic solution of mode-iii crack in damaged softening materials

In this paper, a constitutive model of elasticity coupled with damage suggested by Lemaitre et al, [1] is used. The macroscopic stress-strain response of the model includes two stages: strain hardening and strain softening. The basic equation is derived for the anti-plane shear problem. Several lowest order asymptotic solutions are obtained, and assembled for the crack-tip fields.

Inhomogeneous behavior of antiplane shear crack in softening material with local unloading

Examined in this work is the anti-plane stress and strain near a crack in a material that softens beyond the elastic peak and unloads on a linear path through the initial state. The discontinuity in the constitutive relation is carried into the analysis such that one portion of the local solution is elliptic in character and the other hyperbolic. Material elements in one region may cross over to another as the loading is increased. Local unloading can thus prevail. Presented are the inhomogeneous character of the asymptotic stress and strain in the elliptic and hyperbolic region, in addition to the region in which the material elements had experienced unloading. No one single stress or strain coefficient would be adequate for describing crack instability.

Plane strain asymptotic fields for cracks terminating at the interface between elastic and pressure-sensitive dilatant materials

The plane strain asymptotic fields for cracks terminating at the interface between elastic and pressure-sensitive dilatant material are investigated in this paper. Applying the stress-strain relation for the pressure-sensitive dilatant material, we have obtained an exact asymptotic solution for the plane strain tip fields for two types of cracks, one of which lies in the pressure-sensitive dilatant material and the other in the elastic material and their tips touch both the bimaterial interface. In cases, numerical results show that the singularity and the angular variations of the fields obtained depend on the material hardening exponent n, the pressure sensitivity parameter mu and geometrical parameter lambda.

Interacting multiple cracks in piezoelectric materials

In this paper, a method is presented to calculate the plane electro-elastic fields in piezoelectric materials with multiple cracks. The cracks may be distributed randomly in locations, orientations and sizes. In the method, each crack is treated as a continuous distributed dislocations with the density function to be determined according to the conditions of external loads and crack surfaces. Some numerical examples are given to show the interacting effect among multiple cracks.

On The Interaction Of Channeling/Segmentation Cracks Of Coating: Study On The Critical Spacing

Channeling/segmentation cracks may arise in the coating subjected to in-plane tensile stress. The interaction between these multiple cracks, say the effect of the spacing between two adjacent cracks oil the behaviors of channels themselves and the interface around the interface corners, attracts wide interest. However, if the spacing is greater than a specific magniture,, namely the Critical Spacing (CS), there should be no interaction between such channeling/segmentation cracks. In this study, file mechanism of the effect of the crack spacing oil the interfacial stress around the interface corner will be Interpreted firstly. Then the existence of the CS will be verified and the relationship between the CS and the so-called stress transfer length Ill coating will be established for plane strain condition. Finally, the dependence of the stress transfer length, simultaneously of the CS, on the sensitive parameters will be investigated with finite element method and expressed with a simple empirical formula. (C) 2007 Elsevier Ltd. All rights reserved.

Asymptotic analysis of plane-strain mode I steady-state crack growth in transformation toughening ceramics (II)

Based on a constitutive law which includes the shear components of transformation plasticity, the asymptotic solutions to near-tip fields of plane-strain mode I steadity propagating cracks in transformed ceramics are obtained for the case of linear isotropic hardening. The stress singularity, the distributions of stresses and velocities at the crack tip are determined for various material parameters. The factors influencing the near-tip fields are discussed in detail.

Micromechanics of elastic solids weakened by a large number of cracks

This paper presents a micromechanics analysis of the elastic solids weakened by a large number of microcracks in a plane problem. A new cell model is proposed. Each cell is an ellipse subregion and contains a microcrack. The effective moduli and the stress intensity factors for an ellipse cell are obtained. The analytic closed formulas of concentration factor tensor for an isotropic matrix containing an anisotropic inclusion are derived. Based on a self-consistent method, the effective elastic moduli of the solids weakened by randomly oriented microcracks are obtained.

Complex stress functions and plastic zone sizes for notch cracks subjected to various loading conditions

In this paper, the conformal mapping method is used to solve the plane problem of an infinite plate containing a central lip-shaped notch subjected to biaxial loading at a remote boundary or a surface uniform pressure on the notch. The stress intensity factors KI and KII are obtained by the derived complex stress functions. The simple analytical expressions can be applied to the situation of cracks originating from a circular or an elliptical notch. The plastic zone sizes for such notch cracks are subsequently evaluated in light of the Dugdale strip yield concept. The results are consistent with available numerical data.

Experimental Study on Vortex-Induced Vibrations of a Pipeline Near a Plane Boundary in Steady Flow

对单向水流作用下近壁管道横向涡激振动进行了实验模拟,重点探讨了管道与壁面间隙比(e/D)对管道涡激振动幅值和涡激振动频率响应特性的影响规律.实验结果表明,管道与壁面间隙宽度对管道涡激振动特性有较明显影响.在较大间隙比(e/D＞0.66)下,管道振幅随着Vr数的增大先快速增长到最大值,然后平缓下降;在振动初期(即Vr数较小时),管道振动频率变化基本符合Strouhal规律;在振动中后期(即Vr数较大时),管道振动频率变化不符合Strouhal规律,而在管道固有频率附近缓慢增长.在较小间隙比(e/D＜0.30)下,管道振幅随Vr数的增大先平缓上升到最大值,随后较快速下降;在振动初期,管道振动频率变化不遵循Strouhal规律;在整个振动范围内,与较大间隙比情况相比,随着Vr数增加,管道振动频率增长幅度明显较大.