The solution of some typical examples is given, and the significant effect of the stress intensity factor of the rock on the crack propagation is shown. The influence of layer thickness on the stress intensity. The accuracy of the stress intensity factor calculation is satisfactorily examined for rectangular, pennyshaped and elliptical planar cracks. The energy release rate is directly related to the stress intensity factor associated with a given twodimensional loading mode modei, modeii, or modeiii when the crack grows straight ahead. If the crack is located centrally in a finite plate of width and height, an approximate relation for the stress intensity factor is. Today, it is the critical stress intensity factor k ic, found in the plane strain condition, which is accepted as the defining property in linear elastic fracture mechanics. Therefore there is no need to derive ready made stress intensity factor for each load system and resulting internal stress distribution. Does anyone know how to calculate the stress intensity factor of 3d. However, in most practical cases the deepest point a fig.
This paper presents a threedimensional viscoelastic model to study the interactions of a penny shaped interfacial crack and a center of dilatation in the infinite viscoelastic bimaterial, which can model the rock fracture subjected to stress and thermal dilatation during some engineering process. Fracture mechanics of throughcracked cylindrical pressure. In order to validate the boundary collocation method and the calculation procedure, the mode stress intensity factors for a square plate with a single crack and an airfoil shaped plate with single crack under mixed boundary condition as shown in figures 1a and 2a are calculated. Crack tip stress dependence on the stress intensity factor k s y x w xy v yy v xx m u y r a s s s 3 yy 2 k x v s 2 yy 2 k x v s k 1 stress components, v ij, at the crack tip depend on the stress intensity factors k i which is influenced by. The mode i stress intensity factor along the crack fronts of a rectangular discontinuity in an infinite body is independent of youngs modulus 45.
Stress intensity factors calculation by using abaqus. Stress intensity factor determination plays a central role in linearly elastic fracture mechanics problems. The sifs which are exact for a penny shaped crack are based on the well known solution for a point load acting normally to such a crack. Quantification of fracture interaction using stress intensity factor. It is found that the boundary element scheme can evaluate the stress intensity factor at the crack tip to within an accuracy of 0. Calculate the stress intensity factors by using abaqus. Mechanics of a rock anchor with a penny shaped basal crack. Does anyone know how to calculate the stress intensity factor of 3d penny shaped internal crack under uniaxial compression. The interface pennyshaped crack reconsidered 169 where 41 42 it can be shown that kw,g has the following property we note that in the present form the size of the contact region, l. In order to test the accuracy, we also list the results in tables 3 5. In the former paper the analysis given applies to an axisymmetric distribution of pressure p. The jintegral, stress intensity factors, and t stress should be path independent, and abaqus provides for its evaluation on as many contours as you request. Stress intensity factor determination for threedimensional crack. Stress intensity factor and effective stiffness of a solid.
Stress intensity factors for external and penny shaped cracks in transversely isotropic cylinders subjected to thermal shock. In theory the stress at the crack tip where the radius is nearly zero, would tend to infinity. The stress intensity factor for a penny shaped crack between two dissimilar materials m. Evolution and interaction of surface cracks under thermal. If possible, please recommend some references or classical books which. The pennyshaped crack at a bonded plane with localized. The case of a pennyshaped crack with arbitrary normal displacements prescribed at its faces inclusion problem was solved exactly and in closed form by. Oct 01, 2003 read multiple penny shaped cracks interaction in a finite body and their effect on stress intensity factor, engineering fracture mechanics on deepdyve, the largest online rental service for scholarly research with thousands of academic publications available at your fingertips. The influence of layer thickness on the stress intensity factor of a penny shaped crack in a sandwiched viscoelastic bimaterial. Stress intensity factor of a pennyshaped crack with small.
Fracture propagation is controlled by the stress field near the crack tip. Please help me how can i start and provide me any tutorialexample if you have. Numerical results for the stress intensity factors are. Calculation of stress intensity factors and crack opening a. Introduction the assessment of surface flaws is of great importance in structural integrity analyses.
Numerical calculations of the stress intensity factors are carried out for a beryllium oxide cylinder, which possesses transverse isotropy. Crack tip stress distribution and stress intensity factor. An expression is developed for the stress intensity factor of a pennyshaped crack in an infinite elastic solid subjected to nonaxisymmetric normal loading. The stress intensity is given as stress kiy sqrtc where c is flaw depth and y 1. Stress intensity factors are given for various penny. If the weight function is known for a crack in a component, the stress intensity. Usually these cracks have a shape which is close to that of a semiellipse. Introduction to fracture mechanics david roylance department of materials science and engineering. Kachanovs approximate method for crack interactions int. Although the energy release rate is calculated directly in abaqusstandard, it is usually not straightforward to compute stress intensity factors from a known jintegral for mixedmode problems. Unlike the energy release rates computed by the vcct, the stress intensity factors for interfacial cracks are convergent and not dependent on the virtual crack extension length mesh size a. Sometimes a crack line may form a closed loop for example, when modeling a full penny shaped crack without invoking symmetry conditions. The first contour is normally at the crack tip, and subsequent contours are generated automatically as contours passing through the nearest neighboring elements, moving out from the crack tip. The weight function procedure developed by buckner 8 simplifies the determination of stress intensity factors.
The effects of the material anisotropy on the stress intensity factor and on the crack shape are investigated for a penny shaped crack in a transversely isotropic plate of finite thickness. There are various stress intensity factor solutions, particularly for flat plates and pressure vessels with various cracked geometries. The stress intensity factor at the tip of a pennyshaped crack of radius a \displaystyle a in an infinite domain under uniaxial tension. A mathematical formulation is presented for the dynamic stress intensity factor mode i of a permeable pennyshaped crack subjected to a timeharmonic propagating longitudinal wave in an infinite poroelastic solid. Because the stress field near the crack tip is asymptotic dominant or singular, it is characterized by the stress intensity factor. Normalized stress intensity factors along crack front for penny shaped crack as a function ofthe angle from the symmetry plane e 0 towards the other symmetry plane e 90 and normalized stress intensity factors assuming. This is applicable to cracks under plane stress, plane strain, and antiplane shear. Pdf stress intensity factor for a centrally located. The stress intensity factors along the front of the pennyshaped crack can be found in figure 19.
This often proves difficult to achieve in practice, due to the unbounded nature of the stress field along the crack front. Most of the numerical methods require separate calculation of the stress intensity factor for each given stress distribution and each crack length. Numerical results were obtained to examine the effects of the elastic. Calculation of stress intensity factors by the westergaard method. The stress intensity factor calculation is valid for cracks in homogeneous, linear elastic materials. Pennyshaped crack in an infinite domain under uniaxial tension. The interface penny shaped crack reconsidered 169 where 41 42 it can be shown that kw,g has the following property we note that in the present form the size of the contact region, l a, is not yet determined.
Kachanovs approximate method for crack interactions kachanov, m. This study describes the evolution of the stress intensity factor k i for fissures in brittle elastic materials subjected to transient thermal actions. When the crack is expanding under a constant flow rate, the classical solution is found to be approximately valid for very large cracks, and nevertheless the crack. Stress intensity factors for external and pennyshaped. For a penny shaped crack embedded within an infinite solid and loaded remotely perpendicular to the crack at a stress. Please help me how can i start and please provide me any tutorialexampleguidance if you have. Plastic zone at the tip of a penny shaped crack according to the dugdale model. He compared the ansys results with the crack growth program franc3d. Tweed international journal of engineering science 1973 173 stress intensity factors for some throughcracked fastener holes a. A formula is derived for the stress intensity factor at the rim of a penny shaped crack in an infinite solid in which there is an axisymmetric distributing of body forces acting in a direction normal to the original crack surfaces. The stress intensity factor is generally expressed in the form k y. Selected papers on crack tip stress fields 1997 sanford.
Well known opening mode stress intensity factor variations for a pennyshaped crack situated close and orthogonaly to the surface e. This assignment was conducted through numerical computations of stress intensity factors sifs for two penny shaped cracks in elastic solid estimated by the. The dynamic stress field under the shwaves at the moving crack tip of functionally graded materials is analyzed, and the influences of parameters such as graded parameter, crack velocity, the angle of the incidence and the number of the waves on dynamic stress intensity factor are also studied. Stress intensity factors for cracked finite plates with. The paper presents the integral equations governing the pennyshaped crack problem and the numerical results for the stress intensity factor at the tip of the crack. Stress intensity factors for cracked finite plates with mixed. This page provides stress intensity factor solutions for common cases. The singular state associated with the embedded crack with finite, nonaxisymmetric normal loading is that of plane strain.
Chell international journal of fracture 1976 114 crack tip, stress intensity factors for plane extension and plate bending problems g. The stress intensity factor sif of a half penny shaped crack normal to the interface in the top layer of a threelayer bonded structure is obtained by the finite element method for a wide range of parameters. Tsang et al 6 investigated the sif of parallel penny shaped cracks under itype load. If the crack is not located centrally along the width, i.
In this paper, the stress intensity factor k for the crack front line. Traditionally, direct methods of calculating the stress intensity factor, such as from the crack opening displacement, rely heavily on accurate representation of the stress field near the crack front. The stressintensity factor for a pennyshaped crack between two dissimilar materials. The stress intensity factor can then be determined for penny shaped cracks in infinite or finite solids subjected to symmetric loading about the plane containing the crack. This critical stress intensity factor is then a measure of. An analytical solution for the axisymmetric problem of a. The coefficients of the polynomial p x1, x2 are determined using the boundary conditions on the surfaces of the pennyshaped cracks where the stresses on the crack surfaces are zero.
This integral equation is solved numerically by using gaussian. Calculation of stress intensity factors for interfacial. Stress intensity factors for penny and halfpenny shaped. Kassir department of civil engineering, the city college of the city university of new york, new york, n. Stress intensity factors numerical calculations for two penny shaped. The stress intensity, k i, represents the level of stress at the tip of the crack and the fracture toughness, k ic, is the highest value of stress intensity that a material under very specific planestrain conditions that a material can withstand without fracture. An expression for the stress intensity factor has been found by cherepanov 1974, a result which was generalized for a nonhomogeneous body by sankar and fabrikant 1983 a. Threedimensional viscoelastic interactions of a center of.
The stress intensity factor for a pennyshaped crack in an elastic. Stress intensity factors for cracks in anisotropic. Numerical determination of the stress intensity factor ki in. If the weight function is known for a given cracked body, the stress intensity factor due to any load system applied to the body can be determined by using the same weight function. A generation of special triangular boundary element shape. Application of displacement and traction boundary integral.
Some of these solutions are based on the use of thinshell theory18, which does not take into account the three dimensional nature of the highly localised stresses in the vicinity of the crack front. Stress intensity factors for cracklineloaded edgecrack. In particular, the effect of the waveinduced fluid flow. A simple empirical equation for evaluating the elasticplastic jintegral for surface flaws is proposed. A pathindependent integral for the calculation of stress. The case studied is a penny shaped crack located in a rock mass and through which runs a fluid whose temperature falls in time. Abstract a pathindependent surface integral is derived for the calculation of the stress intensity factor in three dimensions. Stress intensity factor sifbased fracture modeling has received great. A formula is derived for the stress intensity factor at the rim of a pennyshaped crack in an infinite solid in which there is an axisymmetric distributing. The y term allows the effects of finite width, crack shape, position along crack front, bulging, stress concentration factors, local stress concentration. Jan 01, 20 the effect of crack interactions on stress intensity factors is examined for a periodic array of coplanar pennyshaped cracks. The stress intensity factors along the front of the penny shaped crack can be found in figure 19. Figure 2 stresses around half penny shaped crack in bore of rivet hole in the lap joint assembly stress intensity factors caused by the three different contact conditions the sifs are computed directly from the stresses and displacements in the bem.
Stress intensity factors numerical calculations for two penny shaped cracks in the elastic solid author links open overlay panel d. Stress intensity factors for cracks in anisotropic materials. Fracture materials database stress intensity factor solutions. Calculate the stress under point load with no crack use a book like roark for this.
He simulated crack growth and computed the stresses and the stress intensity factor. The stress intensity factors for a periodic array of. Stress intensity factor calculation point in 3d crack problems a cen ter point of the crack fiont element, b side node of the crack fiont edge 105 figure 7. They are simply the real and imaginary parts of a complex stress intensity factor. Some of these solutions are based on the use of thinshell theory18, which does not take into account the three dimensional nature of the highly localised stresses in the vicinity of the crack.
Stress intensity factor determination for threedimensional. Stress intensity factor variation of a pennyshaped crack. The stressintensity factor for a pennyshaped crack. Kachanov 1987 presented an approximate method for calculating.
The stressintensity factor for a pennyshaped crack between two. Stress intensity factors numerical calculations for two. The stress intensity factor for a pennyshaped crack in an. Dynamic stress intensity factor mode i of a permeable.
Does anyone know how to calculate the stress intensity factor. It is also valid for an interfacial crack between two different isotropic linear elastic materials. Therefore, the weight functions for the point a and b have been derived analogously to. I am now following abaqus example problems manual 6. Stress intensity factors for cracklineloaded edge crack specimens by john e. It illustrates the manner in which the stress intensity factor can be in. If stress intensity is less than critical stress intensity factor, kic. A comparison between six models, calculating stress intensity factor sif mode i for central cracked plate with uniform tensile stress, was made by alansari 17 in order to select.
An expression for the surface displacement of the crack is also given. The maximum stress intensity factors with respect to a time variable are approximated by a function of biots number and crack length. Mode i stress intensity factor estimate for a halfpenny. Stress intensity factor for a centrally located penny shaped crack surrounding a spherical cavity in a cylinder subjected to axial load. Calculating stress intensity factor mode i for composite. The stress intensity factor can then be determined for pennyshaped cracks in infinite or finite solids subjected to symmetric loading about the plane containing the crack. Stress intensity factors at any point on the crack front of penny and half penny shaped cracks subjected to stress gradients are presented. Thiscritical stress intensity factor isthenameasureof materialtoughness. A pathindependent integral for the calculation of stress intensity factors in threedimensional. Mar 01, 2016 as a common parameter to character the ability of crack propagation, the stress intensity factor sif is widely used in investigating the interactions among multicracks. Numerical determination of the stress intensity factor ki. Stress intensity factor for a centrally located penny shaped crack. General stress intensity factor expressions for cracks in mode i the stress intensity factor is defined as. Pennyshaped crack embedded in an infinite solid medium 106 figure 7.
A slightly increase of the critical stress intensity factor has been observed with increasing the testing temperature, and this result is in accordance with the model calculations, which indicate a toughening effect, i. Pdf stress intensity factor for a centrally located penny. Does anyone know how to calculate the stress intensity. Dynamic analysis of a piezoelectric material with an.
701 1512 953 141 1307 896 29 926 207 1262 1689 1008 842 581 360 470 1612 1042 535 1670 1229 1011 2 1163 600 1336 1126 302 429 466 313 1157 1466 481 52 1433 1261 87 1266 941 557 492 1420 144 1306