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.     

<Emphasis Type="Italic">T</Emphasis>-stress analysis for a Griffith crack in a magnetoelectroelastic solid

T-stress as an important parameter characterizing the stress field around a cracked tip has attracted much attention. This paper concerns the T-stress near a cracked tip in a magnetoelectroelastic solid. By applying the Fourier transform, we solve the associated mixed boundary-value problem. Adopting crack-faces electromagnetic boundary conditions nonlinearly dependent on the crack opening displacement, coupled dual integral equations are derived. Then, the closed-form solution for the T-stress is obtained. A comparison of the T stresses for a cracked magnetoelectroelastic solid and for a cracked purely elastic material is made. Obtained results reveal that in addition to applied mechanical loading, the T-stress is dependent on electric and magnetic loadings for a vacuum crack.     

The elastic T-stress for slightly curved or kinked cracks

This work presents a solution for the elastic T-stress at the tip of a slightly curved or kinked crack based on a perturbation approach. Compared to other exact or numerical solutions the present solution is accurate for considerable deviations from straightness. The T-stress variation as crack extends along a curved trajectory is subsequently examined. It is predicted that T-stress always keeps negative during crack extension when the crack has an initial negative T-stress. In the case of a positive T-stress and non-zero first and second stress intensity factors initially accompanying the crack, the T-stress is not positive with increasing the extension length until a threshold is exceeded. Based on directional stability criterion with respect to the sign of the T-stress, this result implies that for a straight crack with a positive T-stress, the crack extension path will not turn immediately and instead keep a stable growth until a critical length is reached. This prediction is consistent with experimental observations.     

Elastic stresses in bodies with elliptic,parabolic, and tunnel cracks

The paper gives explicit expressions of the elastic T-stress components T _I, T _II, and T _III for an elliptic crack in an unbounded body under uniform pressure and bending and expressions of all the T-stress components for parabolic and tunnel cracks under uniform loading. These formulas are derived by analyzing the asymptotic behavior of the stress components near the crack front using special harmonic functions. The dependence of the T-stresses on Poisson’s ratio, semiaxes and parametric angle of the elliptic crack is studied. The expressions of T _I, T _II, and T _III for a penny-shaped crack under arbitrary uniform pressure and bending follow as a special case from the respective expressions for an elliptic crack __________ Translated from Prikladnaya Mekhanika, Vol. 43, No. 8, pp. 57–70, August 2007.     

A method for measuring mode I crack tip constraint under static and dynamic loading conditions

A novel experimental technique for measuring crack tipT-stress, and hence in-plane crack tip constraint, in elastic materials has been developed. The method exploits optimal positioning of stacked strain gage rosette near a mode I crack tip such that the influence of dominant singular strains is negated in order to determineT-stress accurately. The method is demonstrated for quasi-static and low-velocity impact loading conditions and two values of crack length to plate width ratios (a/W). By coupling this new method with the Dally-Sanford single strain gage method for measuring the mode I stress intensity factorK _I , the crack tip biaxiality parameter is also measured experimentally. Complementary small strain, static and dynamic finite element simulations are carried out under plane stress conditions. Time histories ofK _I andT-stress are computed by regression analysis of the displacement and stress fields, respectively. The experimental results are in good agreement with those obtained from numerical simulations. Preliminary data for critical values ofK _I and β for dynamic experiments involving epoxy specimens are reported. Dynamic crack initiation toughness shows an increasing trend as β becomes more negative at higher impact velocities.     

The lowest order solution of the crack–inhomogeneity interaction for mode I crack

The lowest order solution to crack–inhomogeneity interaction is derived for mode I crack. The basic solution evaluates the variation of near-tip stress intensity factors induced by an inhomogeneity of arbitrary shape. A set of simplified forms of the basic solution is also obtained for several special inhomogeneity shapes. As validated by numerical examples, the approximate solutions have good accuracy.     

An Experimental Study of Mixed Mode Crack Initiation and Growth in Functionally Graded Materials

Quasi-static mixed mode crack initiation and growth in functionally graded materials (FGMs) was studied through fracture experiments on polymer-based FGMs manufactured by selective ultraviolet irradiation poly(ethylene carbon monoxide)—a photo-sensitive copolymer that becomes more brittle and stiffer under ultraviolet irradiation. The objective of the study was to determine whether crack kinking criteria for homogeneous materials, e.g., maximum hoop stress criterion, also hold for FGMs. Single edge notched tension specimens with different spatial variations of Young's modulus, failure stress and failure strain, were tested. Near tip mode mixity was introduced either by inclining the crack to the remote loading direction, as in the case of homogeneous materials, or to the direction of material gradient, or both. A full-field digital image correlation technique was used to measure in real-time the displacement field around the crack tip while it propagated through the graded material, and to extract the fracture parameters of stress intensity factor K _I and K _II , and the T-stress. It was found that the nonsingular T-stress term in the asymptotic expansion for stresses plays a very important role in accurately measuring fracture parameters. It was also found that the maximum tangential stress criterion can be applied to the case of FGMs to predict crack kinking provided that the effect of the T-stress is accounted for and the process zone size is small compared to the intrinsic material gradient length scale. However, for accurate crack path prediction at a length scale comparable to the material gradient, detailed material property information is required. In general, the crack will propagate towards a region that exhibits less fracture toughness, but, unlike the case of homogeneous materials, along a path where K _II is not necessarily equal to zero.     

Effects of partial crack–face contact for the bending of thin shell structures

The physical occurrence that crack surfaces are in contact at the compressive edges when a flat or a shell is subjected to a bending load has been recognized. This article presents a theoretical analysis of crack–face contact effect on the stress intensity factor of various shell structures such as spherical shell, cylindrical shell containing an axial crack, cylindrical shell containing a circumferential crack and shell with two non-zero curvatures, under a bending load. The formulation of the problem is based on the shear deformation theory, incorporating crack–face contact by introducing distributed force at the compressive edge. Material orthotropy is concerned in this analysis. Three-dimensional finite element analysis (FEA) is conduced to compare with the theoretical solution. It is found that due to curvature effect crack–face contact behavior in shells differs from that in flat plates, in that partial contact of crack surfaces may occur in shells, depending on the shell curvature and the nature of the bending load. Crack–face contact has significant influence on the stress intensity factor and it increases the membrane component but decreases the bending component.     

Green’s function for T-stress of semi-infinite plane crack

Green’s function for the T-stress near a crack tip is addressed with an analytic function method for a semi-infinite crack lying in an elastical, isotropic, and infinite plate. The cracked plate is loaded by a single inclined concentrated force at an interior point. The complex potentials are obtained based on a superposition principle, which provide the solutions to the plane problems of elasticity. The regular parts of the potentials are extracted in an asymptotic analysis. Based on the regular parts, Gre...     

Modeling of crack growth in ductile solids: a three-dimensional analysis

A population of several spherical voids is included in a three-dimensional, small scale yielding model. Two distinct void growth mechanisms, put forth by [Int. J. Solids Struct. 39 (2002) 3581] for the case of a two-dimensional model containing cylindrical voids, are well contained in the model developed in this study for spherical voids. A material failure criterion, based on the occurrence of void coalescence in the unit cell model, is established. The critical ligament reduction ratio, which varies with stress triaxiality and initial porosity, is used to determine ligament failure between the crack tip and the nearest void. A comparison of crack initiation toughness of the model containing cylindrical voids with the model containing spherical voids reveals that the material having a sizeable fraction of spherical voids is tougher than the material having cylindrical voids. The proposed material failure determination method is then used to establish the fracture resistance curve (J–R curve) of the material. For a ductile material containing a small volume fraction of microscopic voids initially, the void by void growth mechanism prevails, which results in a J–R curve having steep slope. On the other hand, for a ductile material containing a large volume fraction of initial voids, the multiple voids interaction mechanism prevails, which results in a flat J–R curve. Next, the effect of T-stress on fracture resistance is examined. Finally, nucleation and growth of secondary microvoids and their effects on void coalescence are briefly discussed.     

T-stress evaluation for slightly curved crack using perturbation method

A general formulation for evaluating the T-stress at crack tips in a curved crack is introduced. In the formulation, a singular integral equation with the distribution of dislocation along the curve is suggested. For a slightly curved crack, a small parameter is generally assumed for the crack configuration. By using the assumption for the small parameter, the perturbation method is suggested and it reduces the singular integral equation into many successive singular integral equations. If the cracked plate has a remote loading and the curve configuration is a quadratic function, the mentioned successive singular integral equations can be solved in a closed form. Therefore, the solution for the T-stress in a closed form is obtained. The obtained results for T-stress are shown by figures. It is found that if the involved parameter is not too small, the influence of the curve configuration is significant. Comparison for T-stresses obtained from a quadratic-shaped curved crack and an arc crack is presented.     

Some plane adiabatic ideal gas flows containing shocks

We obtain the solution describing adiabatic flows of an ideal gas characterized by the two parameters a and b such that [a]=L ^m+1 T ^–1, [b]=ML ^–2–2m where m is arbitrary (m > 0).h This solution permits the construction of flows containing shocks.Translated from Zhurnal Prikladnoi Mekhaniki i Tekhnicheskoi Fiziki, Vol. 10, No. 3, pp. 71–73, May–June, 1969.     

Symmetric-Galerkin boundary element analysis of the dynamic T-stress for the interaction of a crack with an auxetic inclusion

In this work, the symmetric-Galerkin boundary element method (SGBEM) for 2-D elastodynamics in the Laplace-space frequency domain (Laplace domain) is employed to study the dynamic stress intensity factors (DSIFs) and the dynamic T-stress (DTS) during the interaction between a crack and an auxetic inclusion under impact loading conditions. It is found that, while the auxeticity has virtually no effect on the DSIFs, its influence on the DTS is noticeable. This finding is particularly important as it implies the imperative need of fracture criteria based on both the DSIFs and DTS for predicting crack growth in composite materials with auxetic phases.     

Numerical investigation of the effect of rotational relaxation rate on the nitrogen shock wave structure

The nitrogen rotational relaxation process is numerically investigated on the basis of the Monte-Carlo method and the calculation of the classic three-dimensional equations of motion of rigid-rotator molecules. The dependence of the relaxation parameterZ _R on the initial values of the translational and rotational temperaturesT _t andT _R and the asymmetry parameter in the molecular interaction potential is obtained. The problem of the shock wave structure is solved using the Pullin phenomenological model whose parameters are determined on the basis of the solution of the relaxation problem with exact calculation of the molecular interaction. The dependence of the flow fields and the relative thickness of the density profile on the quantityZ _R is obtained and the calculated data are compared with the results of experiments. It is shown that the data are in better agreement for larger values of the asymmetry parameter.Moscow. Translated from Izvestiya Rossiiskoi Akademii Nauk, Mekhanika Zhidkosti i Gaza, No. 4, pp. 159–167, July–August, 1995.     

On the interaction between a dislocation and a circular inhomogeneity with imperfect interface in antiplane shear

The solution of appropriate elasticity problems involving the interaction between inclusions and dislocations plays a fundamental role in many practical and theoretical applications, namely, it increases the understanding of material defects thereby providing valuable insight into the mechanical behavior of composite materials.Although the problem of a three-phase circular inclusion interacting with a dislocation in antiplane shear has been presented [Xiao and Chen, Mech. Mater. 32 (2000) 485], the analysis is limited to the classical perfect bonding condition. The current paper considers the solution for a homogeneous circular inclusion interacting with a dislocation under thermal loadings in antiplane shear. The bonding along the inhomogeneity–matrix interface is considered to be imperfect with the assumption that the interface imperfections are constant. It is found that when the inhomogeneity is soft, regardless of the level of interface imperfection, the inhomogeneity will always attract the dislocation. As a result, no equilibrium positions are available. Alternatively, when the inhomogeneity is hard, an unstable equilibrium position is found which depends on the imperfect interface condition and the shear moduli ratio μ₂/μ₁.     

A Study of Crack–inclusion Interactions and Matrix–inclusion Debonding Using Moiré Interferometry and Finite Element Method

Failure behavior of composite materials in general and particulate composites in particular is intimately linked to interactions between a matrix crack and a second phase inclusion. In this work, surface deformations are optically mapped in the vicinity of a crack–inclusion pair using moiré interferometry. Edge cracked epoxy beams, each with a symmetrically positioned cylindrical glass inclusion ahead of the tip, are used to simulate a compliant matrix crack interacting with a stiff inclusion. Processes involving microelectronic fabrication techniques are developed for creating linear gratings in the crack–inclusion vicinity. The debond evolution between the inclusion–matrix pair is successfully mapped by recording crack opening displacements under quasi-static loading conditions. The surface deformations are analyzed to study evolution of strain fields due to crack–inclusion interactions. A numerical model based on experimental observations is also developed to simulate debonding of the inclusion from the matrix. An element stiffness deactivation method in conjunction with critical radial stress criterion is successfully demonstrated using finite element method. The proposed methodology is shown to capture the experimentally observed debonding process well.

Effect of T-stress on the mode-I fracture toughness of concrete

T-stress expressions are provided for three-point bending (TPB) beams and compact tension (CT) specimens and then its influence on mode I fracture toughness of concrete is investigated. The study shows that T-stress is dependent on the specimen's geometry and the material's property as well, and for TPB and CT specimens of regular size, T-stress is so small that its consequences can be neglected. The study also indicates that concrete specimen size should be carefully chosen to make sure the existence of K-dominance ahead of the crack tip, thus fracture toughness extracted from these specimen configurations can be reliable.     

<Emphasis Type="Italic">T</Emphasis>-stress estimation by the two-parameter approach for a specimen with a V-shaped notch

In the present research, T-stress solutions are provided for a V-shaped notch in the case of surface defects in a pressurised pipeline. The V-shaped notch is analyzed with the use of the finite element method by the Castem2000 commercial software to determine the stress distribution ahead of the notch tip. The notch aspect ratio is varied. In contrast to a crack, it is found that the T-stress is not constant and depends on the distance from the notch tip. To estimate the T-stress in the case of a notch, a novel method is developed, inspired by the volumetric method approach proposed by Pluvinage. The method is based on averaging the T-stress over the effective distance ahead of the notch tip. The effective distance is determined by the point with the minimum stress gradient in the fracture process zone. This approach is successfully used to quantify the constraints of the notch-tip fields for various geometries and loading conditions. Moreover, the proposed T-stress estimation creates a basis for analyzing the crack path under mixed-mode loading from the viewpoint of the two-parameter fracture mechanics.     

Dissociation of a heavy quarkonium in quark–gluon plasma

In this work we calculate the dissociation of heavy mesons such as the upsilon due to absorption of a thermal gluon in a quark–gluon plasma. We also calculate the dissociation of heavy mesons due to the effect of color charge screening. The lifetime of quarkonium moving with velocity v through a quark–gluon plasma at temperature T is computed. An explicit, configuration–space potential is found for the screened interaction between the quarks constituting the meson. This potential is non-spherical, but axially symmetric about the direction of v. We solve the Schrödinger equation for the relative motion of the quarks in this potential, and use the bound-state wavefunction as the initial state for the dissociation of the meson. The meson lifetime is thus determined as a function of v and T, and conclusions are drawn concerning the possibility of detection of the meson in a high-energy heavy-nucleus collision.