An analysis of collective damage for short fatigue cracks based on equilibrium of crack numerical density

Collective damage of short fatigue cracks was analyzed in the light of equilibrium of crack numerical density. With the estimation of crack growth rate and crack nucleation rate, the solution of the equilibrium equation was studied to reveal the distinct feature of saturation distribution for crack numerical density. The critical time that characterized the transition of short and long-crack regimes was estimated, in which the influences of grain size and grain-boundary obstacle effect were investigated. Furthermore, the total number of cracks and the first order of damage moment were discussed.

Experimental and theoretical study on numerical density evolution of short fatigue cracks

Fatigue testing was performed using a kind of triangular shaped specimen to obtain the characteristics of numerical density evolution for short cracks at the primary stage of fatigue damage. The material concerned is a structural alloy steel. The experimental results show that the numerical density of short cracks reaches the maximum value when crack length is slightly less than the average grain diameter, indicating grain boundary is the main barrier for short crack extension. Based on the experimental observations and related theory, the expressions for growth velocity and nucleation rate of short cracks have been proposed. With the solution to phase space conservation equation, the theoretical results of numerical density evolution for short cracks were obtained, which were in agreement with our experimental measurements.

An analysis on overall crack-number-density of short-fatigue-cracks

The evolution of dispersed short-fatigue-cracks is analysed based on the equilibrium of crack-number-density (CND). By separating the mean value and the stochastic fluctuation of local CND, the equilibrium equation of overall CND is derived. Comparing with the mean-field equilibrium equation, the equilibrium equation of overall CND has different forms in the expression of crack-nucleation-rate or crack-growth-rate. The simulation results are compared with experimental measurements showing the stochastic analyses provide consistent tendency with experiments. The discrepancy in simulation results between overall CND and mean-field CND is discussed.

A theoretical analysis on stochastic behaviour of short cracks and fatigue damage of metals

In this paper, the closed form of solution to the stochastic differential equation for a fatigue crack evolution system is derived. and the relationship between metal fatigue damage and crack stochastic behaviour is investigated. It is found that the damage extent of metals is independent of crack stochastic behaviour ii the stochastic deviation of the crack growth rate is directly proportional to its mean value. The evolution of stochastic deviation of metal fatigue damage in the stage close to the transition point between short and long crack regimes is also discussed.

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.

The scattering of general SH plane wave by interface crack between two dissimilar viscoelastic bodies

The scattering of general SH plane wave by an interface crack between two dissimilar viscoelastic bodies is studied and the dynamic stress,intensity factor at the crack-tip is computed. The scattering problem can be decomposed into two problems: one is the reflection and refraction problem of general SH plane waves at perfect interface (with no crack); another is the scattering problem due to the existence of crack. For the first problem, the viscoelastic wave equation, displacement and stress continuity conditions across the interface are used to obtain the shear stress distribution at the interface. For the second problem, the integral transformation method is used to reduce the scattering problem into dual integral equations. Then, the dual integral equations are transformed into the Cauchy singular integral equation of first kind by introduction of the crack dislocation density function. Finally, the singular integral equation is solved by Kurtz's piecewise continuous function method. As a consequence, the crack opening displacement and dynamic stress intensity factor are obtained. At the end of the paper, a numerical example is given. The effects of incident angle, incident frequency and viscoelastic material parameters are analyzed. It is found that there is a frequency region for viscoelastic material within which the viscoelastic effects cannot be ignored.

Simulations and Experiments of Stochastic Characteristics for Collective Short Fatigue Cracks in Steels

Stochastic characteristics prevail in the process of short fatigue crack progression. This paper presents a method taking into account the balance of crack number density to describe the stochastic behaviour of short crack collective evolution. The results from the simulation illustrate the stochastic development of short cracks. The experiments on two types of steels show the random distribution for collective short cracks with the number of cracks and the maximum crack length as a function of different locations on specimen surface. The experiments also give the variation of total number of short cracks with fatigue cycles. The test results are consistent with numerical simulations.

粘弹性材料界面裂纹的瞬态响应和散射研究

本文研究粘弹性材料界面裂纹对冲击载荷的瞬态响应和对广义平面波的稳态散射。相对于已有广泛研究的弹性材料裂纹瞬态响应和稳态散射问题，本文的研究有三个突出特点：1）粘弹性材料；2）界面裂纹；3）广义平面波入射。粘弹性材料界面裂纹对冲击载荷的瞬态响应和对广义平面波的散射尚无开展研究，本文在弹性材料相应问题的研究基础上，首先开展了这一问题的研究。对于冲击载荷下粘弹性界面裂纹的瞬态响应问题，利用Laplace积分变换方法，将粘弹性材料卷积型本构方程转化为Laplace变换域内的代数型本构方程，从而可以在Laplace变换域内象处理弹性材料的冲击响应一样，将相应的混合边值问题归结为关于裂纹张开位移COD的对偶积分方程，并进一步引入裂纹位错密度函数CDD （Crack Dislocation Density），将对偶积分方程化成关于CDD的奇异积分方程（SIE）。用数值方法求解奇异积分方程得到变换域内的动应力强度因子数值解，最后利用Laplace积分逆变换数值方法得到时间域内的动应力强度因子的时间响应。理论分析考虑了两种裂纹模型，即Griffith界面裂纹和柱面圆弧型界面裂纹。考虑的载荷包括反平面冲击载荷和平面冲击载荷。对于平面冲击载荷，通过对裂尖应力场的奇性分析，首次发现粘弹性界面裂纹裂尖动应力场奇性指数不是常数0.5，而是与震荡指数一样依赖材料参数。针对反平面冲击载荷给出了一个算例，计算了裂尖动应力强度因子的时间响应，并与弹性材料的结果作了比较，发现粘弹性效应的影响不仅使过冲击峰值降低，而且使峰值点后移。粘性效应较大时，过冲击现象甚至不出现。关于粘弹性界面裂纹对广东省义平面波的散射问题，首先研究广义平面波在无裂纹存在的理想界面的反射和透射，再研究由于界面裂纹的存在而产生的附加散射场。利用粘弹性材料的复模量理论，可将粘弹性材料的卷积型相构方程化成频率域内的代数型本构方程。类似弹性平面波的处理，在频率域内将问题最终归结为关于裂纹位错密度CDD的奇异积分方程。数值方法求解奇异积分方程即可得到频率域内的散射场，并进而得到裂尖动应力强度因子和远场位移型函数和散射截面。理论分析考虑了两种裂纹模型：Griffith界面裂纹和柱面圆弧型界面裂纹。研究的入射波有广义的SH波和P波。对于广义平面P波入射的情况，通过对裂尖应力场的奇性分析，同样发现粘弹性界面裂纹裂尖动应力场奇性指数不地常数0.5，而是与震荡指数一样依赖于材料参数。对柱面裂纹散射远场的渐近分析，发现远场位移和应力除含有几何衰减因子外，都含有一个材料衰减速因子。散射截面由于材料衰减因子的存在也成为依赖散射半径的量。为了使散射截面仍有意义，文中提出一种修正办法。对Griffith界面裂纹，给出了一个广义平面SH波入射的算例；对柱面界面裂纹，给出了一个广义平面P波入射的算例。计算了不同入射角和入射频率下裂纹的张开位移和动就应力强度因子，并分析了其依赖关系。求解奇异积分方程的数值方法和Laplace积分逆变换数值方法是本文的基本数值方法。本文对这两种方法作了大量的调研和系统的研究。在对比分析的基础上，对现有的各种方法从原理，适用范围，计数效率，优势及特点进行了归纳总结。并尝试了奇异积分方程的最新数值方法--分片连续函数法，证实了其适用性和方便性．

Dynamic behavior of a cylindrical crack in a functionally graded interlayer under torsional loading

In this paper the problem of a cylindrical crack located in a functionally graded material (FGM) interlayer between two coaxial elastic dissimilar homogeneous cylinders and subjected to a torsional impact loading is considered. The shear modulus and the mass density of the FGM interlayer are assumed to vary continuously between those of the two coaxial cylinders. This mixed boundary value problem is first reduced to a singular integral equation with a Cauchy type kernel in the Laplace domain by applying Laplace and Fourier integral transforms. The singular integral equation is then solved numerically and the dynamic stress intensity factor (DSIF) is also obtained by a numerical Laplace inversion technique. The DSIF is found to rise rapidly to a peak and then reduce and tend to the static value almost without oscillation. The influences of the crack location, the FGM interlayer thickness and the relative magnitudes of the adjoining material properties are examined. It is found among others that, by increasing the FGM gradient, the DSIF can be greatly reduced.

Electromechanical influence of crack velocity at bifurcation for poled ferroelectric materials

The piezoelastodynamic field equations are solved to determine the crack velocity at bifurcation for poled ferroelectric materials where the applied electrical field and mechanical stress can be varied. The underlying physical mechanism, however, may not correspond to that assumed in the analytical model. Bifurcation has been related to the occurrence of a pair of maximum circumferential stress oriented symmetrically about the moving crack path. The velocity at which this behavior prevails has been referred to as the limiting crack speed. Unlike the classical approach, bifurcation will be identified with finite distances ahead of a moving crack. Nucleation of microcracks can thus be modelled in a single formulation. This can be accomplished by using the energy density function where fracture initiation is identified with dominance of dilatation in relation to distortion. Poled ferroelectric materials are selected for this study because the microstructure effects for this class of materials can be readily reflected by the elastic, piezoelectic and dielectric permittivity constants at the macroscopic scale. Existing test data could also shed light on the trend of the analytical predictions. Numerical results are thus computed for PZT-4 and compared with those for PZT-6B in an effort to show whether the branching behavior would be affected by the difference in the material microstructures. A range of crack bifurcation speed upsilon(b) is found for different r/a and E/sigma ratios. Here, r and a stand for the radial distance and half crack length, respectively, while E and a for the electric field and mechanical stress. For PZT-6B with upsilon(b) in the range 100-1700 m/s, the bifurcation angles varied from +/-6degrees to +/-39degrees. This corresponds to E/sigma of -0.072 to 0.024 V m/N. At the same distance r/a = 0.1, PZT-4 gives upsilon(b) values of 1100-2100 m/s; bifurcation angles of +/-15degrees to +/-49degrees; and E/sigma of -0.056 to 0.059 V m/N. In general, the bifurcation angles +/-theta(0) are found to decrease with decreasing crack velocity as the distance r/a is increased. Relatively speaking, the speed upsilon(b) and angles +/-theta(0) for PZT-4 are much greater than those for PZT-6B. This may be attributed to the high electromechanical coupling effect of PZT-4. Using upsilon(b)(0) as a base reference, an equality relation upsilon(b)(-) < upsilon(b)(0) < upsilon(b)(+) can be established. The superscripts -, 0 and + refer, respectively, to negative, zero and positive electric field. This is reminiscent of the enhancement and retardation of crack growth behavior due to change in poling direction. Bifurcation characteristics are found to be somewhat erratic when r/a approaches the range 10(-2)-10(-1) where the kinetic energy densities would fluctuate and then rise as the distance from the moving crack is increased. This is an artifact introduced by the far away condition of non-vanishing particle velocity. A finite kinetic energy density prevails at infinity unless it is made to vanish in the boundary value problem. Future works are recommended to further clarify the physical mechanism(s) associated with bifurcation by means of analysis and experiment. Damage at the microscopic level needs to be addressed since it has been known to affect the macrocrack speeds and bifurcation characteristics. (C) 2002 Published by Elsevier Science Ltd.

Torsional impact response of a penny-shaped interface crack in bonded materials with a graded material interlayer

In this paper, the dynamic response of a penny-shaped interface crack in bonded dissimilar homogeneous half-spaces is studied. It is assumed that the two materials are bonded together with such a inhomogeneous interlayer that makes the elastic modulus in the direction perpendicular to the crack surface is continuous throughout the space. The crack surfaces art assumed to be subjected to torsional impact loading. Laplace and Hankel integral transforms are applied combining with a dislocation density,function to reduce the mixed boundary value problem into a singular integral equation with a generalized Cauchy kernel in Laplace domain. By solving the singular integral equation numerically, and using a numerical Laplace inversion technique, the dynamic stress intensity factors art obtained. The influences of material properties and interlayer thickness on the dynamic stress intensity factor are investigated.

Ductile crack propagation with the strain-energy density criterion

The strain energy density criterion is used to characterize subcritical crack growth in a thin aluminum alloy sheet undergoing general yielding. A finite element analysis which incorporates both material and geometrical nonlinear behaviors of the cracked sheets is developed to predict fracture loads at varying crack growth increments. The predicted results are in excellent agreement with those measured experimentally, thus confirming the validity of the strain energy density criterion for characterizing ductile crack propagation.

The permanence of density waves of galaxy and the numerical investigation of waser mechanism

This paper deals with in detail the permanence of the spiral structure of galaxies andthe characters of waser mechanism. A simplified model of galaxy is adopted. Variousdynamical characters of density waves are studied using numerical calculation method. Theresults verify very well the switch character f waser and the tunnel effect of density wavesat the potential barrier of corotation circle as is shown in a previous work of the author.

Numerical investigation of crack growth in concrete subjected to compression by the generalized beam lattice model

The beam lattice-type models, such as the Euler-Bernoulli (or Timoshenko) beam lattice and the generalized beam (GB) lattice, have been proved very effective in simulating failure processes in concrete and rock due to its simplicity and easy implementation. However, these existing lattice models only take into account tensile failures, so it may be not applicable to simulation of failure behaviors under compressive states. The main aim in this paper is to incorporate Mohr-Coulomb failure criterion, which is widely used in many kinds of materials, into the GB lattice procedure. The improved GB lattice procedure has the capability of modeling both element failures and contact/separation of cracked elements. The numerical examples show its effectiveness in simulating compressive failures. Furthermore, the influences of lateral confinement, friction angle, stiffness of loading platen, inclusion of aggregates on failure processes are respectively analyzed in detail.

Complex potentials and singular integral equation for curve crack problem in antiplane elasticity

Four types of the fundamental complex potential in antiplane elasticity are introduced: (a) a point dislocation, (b) a concentrated force, (c) a dislocation doublet and (d) a concentrated force doublet. It is proven that if the axis of the concentrated force doublet is perpendicular to the direction of the dislocation doublet, the relevant complex potentials are equivalent. Using the obtained complex potentials, a singular integral equation for the curve crack problem is introduced. Some particular features of the obtained singular integral equation are discussed, and numerical solutions and examples are given.