脆性材料裂纹扩展的分形运动学

大量实验结果表明脆性材料裂纹扩展路径是弯弯曲曲的，不规则的，断裂表面粗糙不平，表现出分形特征。本文应用分形几何推导出脆性材料裂纹扩展的分形运动学公式，得到了动、静态断裂韧性与分形裂纹扩展速度、裂纹长度、微结构特征量δ＿ｃ以及分形维数之间的关系；推导出了裂纹扩展速度的计算公式，得出宏观量测的裂纹扩展速度Ｖ＿０应小于分形裂纹扩展速度Ｖ；并推导出直接根据材料动、静态断裂韧性和裂纹扩展速度Ｖ＿０估算材料裂纹扩展路径的分维计算公式．本文理论分析结果与实验值在定性定量上达到了较好的一致。     

Limiting velocity of crack propagation in dynamically fractured materials

A model of a fractal crack is considered. It is found that the limiting velocity of crack propagation is determined by the fractal dimension of the crack contour. It is shown that for commercial steels, the limiting crack velocity is in the range V _lim = (0.155–0.537)c ₁ (c ₁ is the speed of sound).     

不同加载速率下砂岩弯曲破坏的细观机理

岩石细观破裂形貌是岩石破坏机制的重要反映,为研究不同加载速率对砂岩弯曲破坏的影响,通过三点弯曲实验和扫描电镜方法,对某煤矿关键层砂岩弯曲破断裂纹细观形态以及裂纹的自相似性进行了研究。选取6个不同加载速率对岩样进行三点弯曲实验,观察其宏观断裂情况,并利用扫描电镜对弯曲断裂面表面裂纹细观结构进行观察,并拍摄不同倍数下的扫描电镜图片。对图片进行图像处理后得到砂岩弯曲断裂破坏细观裂纹信息,并计算得到微裂纹的分形盒维数值。结果显示:随着加载速率的提高,砂岩穿晶断裂的比例也随之升高,裂纹分形维数亦随着加载速率的增大而增加,同时,分形维数还与弯曲断裂破坏荷载和抗弯强度成正比。可见,加载速率对断裂方式有一定的影响,且加载速率越大断裂所需的破坏能越大,裂纹分布越广,表明开采速度与岩爆等岩体动力灾变有密切关系。     

Asymptotic analysis of the finite cone-and-plate flow of a non-newtonian fluid

Consider the shearing flow of a viscoelastic fluid trapped by surface tension between a cone and a plate. An asymptotic analysis of this problem in the limit of small gap angle has been done. This limit is realized in many practical situations. It is assumed that the Deborah number De, the Reynolds number Re, and the retardation parameter β are all order unity and that the shape of the free surface is very nearly spherical. Closed form analytic expressions are obtained for the leading terms of the primary and weak secondary motion of the fluid as well as the meniscus shape. It is found that the velocity field is bounded and continuous if and only if . There is a family of curves in the De-β plane on which the velocity field has a removable singularity at the origin. The secondary flow is made up of either one or two toroidal vortices. The meniscus has a bulge near the rotating cone and a trough near the stationary plate.     

Macroscopic and microscopic examination of the relationship between crack velocity and path and Rayleigh surface wave speed in single crystal silicon

The speed of Rayleigh surface waves, denoted C_R, is the accepted upper limit for Mode I crack velocity in monolithic solids. In the current contribution, we discuss several critical issues associated with the velocity of Rayleigh surface waves and crack velocity in single crystal (SC) brittle solids, and the global and local influence of C_R on crack path selection in particular.Recent cleavage experiments in SC silicon showed that crack velocity at certain cleavage planes and crystallographic orientations cannot exceed a small fraction of C_R, and thereafter the crack deflects to other cleavage planes. Indeed, C_R defined by the continuum mechanics ignores atomistic phenomena occurring during rapid crack propagation, and therefore is limited in predicting the crack velocity. Examination of these anomalies shows that this limitation lies in microstructural lattice arrangement and in anisotropic phonon radiation during rapid crack propagation. Globally, C_R has no influence on the crack deflection phenomenon. However, the misfit in C_R between the original plane of propagation and the deflected plane generates local instabilities along the deflection zone.     

Effects of surface diffusion and resolved shear stress intensity factor on environmentally assisted cracking

Hydrogen induced crack-tip plastic deformation has been known as the primary mechanism of hydrogen assisted cracking and stress corrosion cracking. It has been systematically shown that the same mechanism of environmentally assisted crack-tip dislocation emission causes hydrogen assisted cracking, stress corrosion cracking, and liquid metal embrittlement cracking.An embrittling chemical species has to reach a crack tip in order to accelerate crack growth. Very close to a sharp crack tip, surface diffusion is shown to be the dominant transport process of embrittling species for stage-II crack growth. The role of surface diffusion in stage II crack growth is analyzed. The constant cracking velocity is proportional to the surface diffusion coefficient of an embrittling species and inversely proportional to a length parameter, , which is related to the transport process upstream.Dislocation emission at a crack tip is driven by crack-tip resolved shear stress. Crack-tip resolved shear stress field is characterized by resolved shear stress intensity factor, K_RSS·K_RSS is defined, the procedure for its calculation outlined, and its applications to crack-tip dislocation emission and environmentally assisted cracking discussed.     

A new experimental method in determining crack propagation transition temperature in steel

This study concentrates on a new experimental method for determining propagation transition temperature (PTT) in carbon steel (0.45% C). The method is based on the Hopkinson bar concept, by which the nature of pressure-time relationship can be studied when a pulse is applied at the active end of the bar. This concept enables one to determine terminal crack velocity using the round tensile specimen of a proper geometry with the V-notch at the half of its length; the total speimen length is 900 mm. When such a specimen is pulled in a testing machine, the fracture occurs at the notch cross section. The fracture event produces elastic waves which can be analyzed, and thus, the crack velocity can be estimated. This technique has been applied in several temperatures to monitor crack velocity changes. The analysis of oscillograms indicated a sharp change in the average crack velocity as well as in the average acceleration within a very narrow range of temperatures. It has been observed that the propagation transition temperature occurs within the limits ΔT ≈ 10 K, and for steel 0.45%C PPT ≈ 252 K.     

Fatigue damage evolution in silicon films for micromechanical applications

In this paper we examine the conditions for surface topography evolution and crack growth/fracture during the cyclic actuation of polysilicon microelectromechanical systems (MEMS) structures. The surface topography evolution that occurs during cyclic fatigue is shown to be stressassisted and may be predicted by linear perturbation analyses. The conditions for crack growth (due to pre-existing or nucleated cracks) are also examined within the framework of linear elastic fracture mechanics. Within this framework, we consider pre-existing cracks in the topical SiO₂ layer that forms on the Si substrate in the absence of passivation. The thickening of the SiO₂ that is normally observed during cyclic actuation of Si MEMS structures is shown to increase the possibility of stable crack growth by stress corrosion cracking prior to the onset of unstable crack growth in the SiO₂ and Si layers. Finally, the implications of the results are discussed for the prediction of fatigue damage in silicon MEMS structures.     

Near threshold fatigue behaviour of flake graphite cast irons microstructures

This paper describes the influence of material toughness degradation, through reversed temper embrittlement (RTE) and mean stress on the near threshold fatigue crack growth characteristics of a CrMoV turbine bolting steel at ambient and elevated temperatures. It was established at ambient temperatures that strong effects of R-ratio and material condition (toughness) were observed on near threshold fatigue crack growth characteristics. At elevated temperatures it was shown that for the non-embrittled material that only under low R-ratio conditions did increased temperature increase the level of threshold stress intensity ΔK_th, by some 20%. In the case of embrittled material, increasing the temperature increased ΔK_th levels by around 30% and decreased near threshold growth rates by an order of magnitude at low to intermediate R-ratio levels.The effects of R-ratio on ΔK_th for all material and mechanical testing conditions could be simply expressed by the difference between ΔK_th at R = O and a constant B multiplied by R.Quantitative fractographic observations indicated that, generally, the incidence of intergranular failure prevalent in embrittled and non-embrittled steels exhibited a maximum at some specific ΔK level. Also in embrittled steels large effects of environmental assisted crack (EAC) growth were observed at near threshold fatigue crack growth rates. It was suggested that this was the result of the much reduced material cohesive strength which was caused by the presence of both impurity and hydrogen atoms.     

Effects of fractal crack

Experimental results indicate that propagation paths of cracks in geomaterials are often irregular, producing rough fracture surfaces which are fractal. In this paper, crack tip motion along a fractal crack trace is discussed. A fractal kinking model of the crack extension path is established to describe irregular crack growth. The length, velocity and kinking effects of the fractal crack are analysed. A formula is derived to describe the effects of fractal crack propagation on the dynamic stress intensity factor and on crack velocity. Finally, expressions of stress and displacement fields near the fractal crack tip are given.     

Effect of T-stress on crack growth under mixed mode I–III loading

For a crack subjected to combined mode I and III loading the influence of a T-stress is analyzed, with focus on crack growth. The solid is a ductile metal modelled as elastic–plastic, and the fracture process is represented in terms of a cohesive zone model. The analyzes are carried out for conditions of small scale yielding, with the elastic solution applied as boundary conditions on the outer edge of the region analyzed. For several combinations of the stress intensity factors K_I and K_III and the T-stress crack growth resistance curves are calculated numerically in order to determine the fracture toughness. In all situations it is found that a negative T-stress adds to the fracture toughness, whereas a positive T-stress has rather little effect. For given values of K_I and T the minimum fracture toughness corresponds to K_III = 0.     

Fatigue crack growth studies in rail steels and associated weld metal

Fatigue crack growth studies in rail steels and associated weld metal have shown that (a) deformed rail steel exhibited fatigue crack growth rates that are slightly faster than undeformed rail steel and (b) weld metal growth data are appreciably faster than rail steel growth results and exhibit growth rate plateaux that reside above the upper bound reported for rail steel fatigue crack growth.In rail steel microstructures at low ΔK levels fatigue crack extension occurred by a ductile striated growth mechanism. However at K_max values approaching 40 MPa √m transgranular cleavage facets initially formed and their incidence increased with K_max until final fast fracture. The average cleavage facet size agreed well with pearlite nodule dimensions of 60–100 μm.The weld metal microstructure was much coarser than the rail steel and contained highly directional columnar grain growth. At all ΔK levels the dominant fracture mode was transgranular cleavage containing small isolated regions of ductile striated fatigue crack growth. The cleavage facet size varied from 150 to 600 μm; such a large variation was explained by the fact that in general crack extension tended to occur in association with the proeutectoid ferrite phase.     

Numerical Estimation of the Degree of Reservoir Permeability Heterogeneity Using Pressure Drawdown Tests

This study aims to correlate the response of pressure transient test to permeability distribution type. For this purpose, correlated permeability distributions in x–y direction are generated using fractional Brownian motion (fBm) as it has been shown in literature that permeability in carbonate reservoirs exhibits an fBm type distribution horizontally. 2-D fBm permeability distributions created using mid point displacement method are employed as data to a black oil simulator. The intermittence exponent, H or fractal dimension of the distribution, D, as defined by D=2 – H, characterizes the distribution type. All permeability distributions are normalized to represent the same arithmetic mean (20, 100, and 500 mD) and uniform variance so that only their fractal dimension that underlies the smoothness of the distribution distinguishes them. Many different realizations of permeability distributions are generated based on the random number seeds used and pressure transient (drawdown) tests are simulated using a black oil simulator package (ECLIPSE 100). Pressure transient analysis is performed using PanSystem package. As a base case and for the comparison purpose, the same procedure is repeated for the totally homogeneous case (the same permeability for all grids) and a random (normally distributed) permeability distribution with the same mean and uniform variance. The effects of permeability distribution type on the pressure response are clarified. A strong impact of heterogeneity is observed as an increase in skin effect with increasing fractal dimension of permeability distribution. This additional (or pseudo) skin effect due to heterogeneity is correlated to the fractal dimension of the permeability distribution. As a further step, the procedure is repeated for different flow rates applied during the drawdown test. The correlation between the fractal dimension of permeability distribution and additional skin is improved by incorporating the rate into it. The methodology followed can be used in the assessment of reservoir heterogeneity quantitatively using pressure transient response.     

Fatigue crack extension behaviour in porous plasma spray 80%Ni---20%Cr material: The influence of R-ratio

This paper attempts to describe the fatigue crack growth response of a plasma spray 80%---20%Cr material, utilised in the corrosion protection of engineering components, whose microstructure consisted of (i) an austenitic matrix, (ii) a secondary dispersion of a chromite non-metallic inclusion phase and (iii) regions of closed and connected porosity.It was demonstrated that little or no effect of R-ratio was observed on the threshold stress intensity range ΔK_th, which was attributed to both the materials fine to intermediate grain size and probable plain stress testing conditions which significantly decrease crack closure effects. At intermediate fatigue crack growth rates high ratio results were an order of magnitude faster than the low R-ratio data. This was the result of the high R-ratio case seeking out more regions of porosity which then increased the local ΔK on the remaining ligaments leading to accelerated crack growth rates.Porosity was shown to significantly decrease the value of ΔK_th and the extent of porosity observed on fatigue fracture surfaces increased with ΔK level and was well in excess of that of 5% recorded by metallography. Hence the growing fatigue crack preferentially sought out regions of porosity as they represented locations of low fracture energy.     

Limitations of growth regimes for cracks initiated in smooth fatigue specimens

The stages of the growth of small cracks initiating at natural flaws in smooth specimens subjected to fatigue loading are characterized and the dominant propagation mechanisms and corresponding fracture paths are described. Characteristic crack lengths are introduced to assess the transition between the regimes of microstructurally small cracks, physically small cracks and long cracks. A log Δσ-log a-diagram is used to derive estimates of these crack lengths. It is shown that simple formulate can be found which relate these characteristic crack lengths to mechanical and material parameters that can be measured using standard fracture mechanics specimens and fatigue tests.     

Edge fracture in cone-plate and parallel plate flows

Edge fracture is an instability of cone-plate and parallel plate flows of viscoelastic liquids and suspensions, characterised by the formation of a `crack' or indentation at a critical shear rate on the free surface of the liquid. A study is undertaken of the theoretical, experimental and computational aspects of edge fracture. The Tanner-Keentok theory of edge fracture in second-order liquids is re-examined and is approximately extended to cover the Criminale-Ericksen-Filbey (CEF) model. The second-order theory shows that the stress distribution on the semi-circular crack is not constant, requiring an average to be taken of the stress; this affects the proportionality constant, K in the edge fracture equation −N _2c = KΓ/a, where N _2c is the critical second normal stress difference, Γ is the surface tension coefficient and a is the fracture diameter. When the minimum stress is used, K = 2/3 as found by Tanner and Keentok (1983). Consideration is given to the sources of experimental error, including secondary flow and slip (wall effect). The effect of inertia on edge fracture is derived. A video camera was used to record the inception and development of edge fracture in four viscoelastic liquids and two suspensions. The recorded image was then measured to obtain the fracture diameter. The edge fracture phenomenon was examined to find its dependence on the physical dimensions of the flow (i.e. parallel plate gap or cone angle), on the surface tension coefficient, on the critical shear rate and on the critical second normal stress difference. The critical second normal stress difference was found to depend on the surface tension coefficient and the fracture diameter, as shown by the theory of Tanner and Keentok (1983); however, the experimental data were best fitted by the equation −N _2c = 1.095Γ/a. It was found that edge fracture in viscoelastic liquids depends on the Reynolds number, which is in good agreement with the inertial theory of edge fracture. Edge fracture in lubricating grease and toothpaste is broadly consistent with the CEF model of edge fracture. A finite volume method program was used to simulate the flow of a viscoelastic liquid, obeying the modified Phan-Thien-Tanner model, to obtain the velocity and stress distribution in parallel plate flow in three dimensions. Stress concentrations of the second normal stress difference (N ₂) were found in the plane of the crack; the velocity distribution shows a secondary flow tending to aid crack formation if N ₂ is negative, and a secondary flow tending to suppress crack formation if N ₂ is positive. Received: 4 January 1999 Accepted: 19 May 1999     

Study on stress wave propagation in fractured rocks with fractal joint surfaces

This paper presents the experimental and theoretical investigation of property of stress wave propagation in jointed rocks by means of SHPB technique and fractal geometry method. Our aim focuses on the influence of the rough joint surface configuration on stress wave propagation. The comparison of behavior of reflection and transmission waves, deformation and energy dissipation of a rough joint surface characterized by its fractal feature with that of a smooth plane joint has been carried out. It has shown that the rough joint surface distinctly affects the stress wave propagation and energy dissipation in the jointed rocks. The rougher the joint surface was, the more permanent deformation occurred and the more attenuation stress wave took place as well. A nonlinear relationship between the normalized energy dissipation ratio W_J/W_I of the jointed rock and the joint roughness in terms of the fractal dimension has been formulated. It seems that the ratio W_J/W_I, presenting how much energy has been dissipated in the joint, nonlinearly increased with the increment of the fractal dimension D of the jointed surface. The ratio W_J/W_I of a roughly jointed rock, however, tends to be the same as that of a smoothly jointed rock if the fractal dimension is less than a critical value D_c = 2.20. The energy dissipation ratio at the critical point D_c seem to be a constant, not dependent of rock type but fractal joint configuration.     

Geometric Characteristics of Fracture‐Associated Space and Crack Propagation in a Material

It is shown that fracture can be treated as a process occurring in the Finsler space. The use of the Finsler space allows one to construct a delaminated manifold whose characteristics are related to the defect structure of the medium. A method of determining the fractal dimension of fracture is developed using the concept of crack propagation along geodesics.     

Enhancement of fracture toughness of steel with defects by warm-prestressing

The technique of warm-prestressing to improve the resistance of structural steel with defects against low temperature fracture has received considerable attention. It is found that warm-prestressing can improve the fracture toughness and change the COD or δ_c, especially the crack tip plastic opening δ_p.The experimental results obtained from three-point bending tests of 42Mn2 steel specimens at −60°C and −20°C are analyzed. Experiments are also made on the bursting of pressure vessels manufactured from #20 steel. The results indicate that warm-prestressing at room temperature increased the bursting pressure at −40°C for d/t = 0.2 to 0.4, where d is the depth of surface crack and t the vessel thickness.     

An elastic-plastic analysis of an asymmetric cracked specimen subjected to longitudinal shear

Some recent elastic-plastic analyses of cracked specimens subjected to symmetric mode III loading are extended to include asymmetric loading and geometry. Solutions are given for arbitrary work hardening behaviour in any specimen that is amenable to a linear elastic analysis. It is shown that asymmetry has a major influence on the shape of the plastic zone, but does not affect the J-integral unil the loading is well into the large scale yielding range. In particular the “plastic zone corrected” estimate of J, obtained by elastically solving a problem for a crack longer than the actual one, is shown to remain a valid two-term asymptotic expansion in the presence of asymmetry. The general results are applied to a crack at an angle to a uniform stress field in a power law hardening material. The growth of the plastic zone is displayed graphically for various hardening exponents and crack orientations. No other asymmetric solution is available, but values of J are compared with those obtained from a fully plastic analysis in the symmetric case.