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. The stress intensity factor for the crack face pressure case is also shown in figure 5a. The finite element meshes for these models are generated using ansys. 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. An interaction integral method for computing mixed mode. It illustrates the manner in which the stress intensity factor can be in.
In 1957, irwin 2 introduced the sifs to describe the stress and displacement. Mode i stress intensity factor estimate for a halfpenny shaped crack in a trilayer structure by reduction to a bilayer model engineering fracture mechanics january 1, 2009. Because of the importance of sif, its solutions have b een paid very high attention since the beginning of fracture mechanics. Ansys is used for the initial elastic stress analysis of the uncracked model, and.
The problem is solved by the method of boundary integral equations using an iterative algorithm. Lehigh university, institute of fracture and solid mechanics, 1973. The effect of crack interactions on stress intensity factors is examined for a periodic array of coplanar penny shaped cracks. Durability evaiuiiti0n 6tceramic components using careslife. Crack propagation composite interface ansys interface failure modes dissertations, academic mechanical engineering masters usf title effect of fiber volume fraction on fracture mechanics in continuously reinforced fiber composite materials aggregation usf. In this study, a tetrahedral enriched element method and related applications are presented that demonstrate employment of fully unstructured. The stress intensity factor at the tip of a pennyshaped crack of radius in an infinite domain under uniaxial tension is k i 2. The paper concerns the contact interaction of the opposite faces of a pennyshaped crack under harmonic wave at oblique incidence.
The elastic stress field around a crack tip 3 brittle fracture in a solid in the form of crack growth is governed by the stress. In this paper, the mintegral is extended for calculating intensity factors for cracked piezoelectric ceramics using the exact boundary conditions on the crack faces. A physical modeling of threedimensional solid media by an idealized mathematical domain that occupies the full space is standard and widely used when inputs and responses of inte. The obtained results andor conclusions could be of particular interest to the analysis and design of smart sensorsactuators constructed from magnetoelectroelastic composite laminates. Interacting cracks analysis using finite element method. Materials free fulltext determination of the enhancement or. For example, ansys only provides automatic meshing capability for two dimensional problems. The finite element model is established by ansys 12.
Pdf stress intensity factor for a centrally located. Based on linear elastic fracture mechanics, stress intensity factors along the crack front of the pennyshaped crack are evaluated as functions of the crackfront. These include an inclined through crack, an embedded pennyshaped crack, an. Based on linear elastic fracture mechanics, stress intensity factors along the crack front of the pennyshaped crack are evaluated as functions of the crack front. The commonly used methods are based on the stress intensity factor dominated power law. The poling direction is taken at an angle to the crack faces within the plane. Thus, the analysis of stresses near the crack tip constitutes an essential part of fracture mechanics. The condition of continuity between elements is the most important. This paper presents the results of investigations conducted to analyze the singular stress field and the associated stress intensity factor sif for a centrally located penny shaped crack around a. Calibration of threedimensional sif using ansys sif for a pennyshaped crack in a finiteradius cylinder submodel method this is a simple threedimensional crack problem in finite domain, a pennyshaped crack in a finiteradius cylinder subjected to remote uniform tension. Kamaya 5,6 performed the linearelastic and elasticplastic analysis by the finite element method for the interactions between the semicircular and semielliptical surface cracks under a tensile or bending load, and obtained the relationship between.
Formulation for the mode i stress intensity factor. Ksolutions and fatigue crack growth in fcc single crystal superalloys by srikant r lii in. Fracture problems with ansys a very brief introduction. For example, ansys only provides automatic meshing capability for twodimensional problems. In this section, we apply the approximate method of kachanov and laures 1989 to the analysis of interacting cracks to obtain the mode i sif of an array of periodic coplanar pennyshaped cracks in an infinite medium subjected to a remote normal traction at infinity, fig. The bulk of the experimental work57,9 has dealt with the effects of stress concentrations of fatigue crack. Kassir mk, bregman am 1972 the stress intensity factor for a pennyshaped crack between two. A penny shaped crack on the 111 plane evolved into a hexagonal crack front, because. Stress intensity factor for steel hollow pipe with axial crack duration. A pennyshaped crack is introduced on the optical fiber. Introduction to fracture mechanics david roylance department of materials science and engineering.
The rest of the local model is modeled by using 20node three dimensional bri ck elements. Stress intensity factor for a centrally located penny shaped crack. In the numerical procedure normalized complex stress intensity. Kic is a measure of the fracture toughness of the material. In this study the fracture mechanics parameters, including the strain energy release rate, the stress intensity factors and phase angles, along the curvilinear front of a threedimensional bimaterial interface crack in electronic packages are considered by using finite element method with the virtual crack closure technique vcct. In the absence of automated and customized methods and tools, some of todays existing methods for solving threedimensional fracture problems require comprehensive finite element meshing, laborintensive analysis and postprocessing efforts. Fatigue analysis fatigue life of components using ansys high cycle fatigue 6. Analysis of stress intensity factors for threedimensional. This paper presents a threedimensional viscoelastic model to study the interactions of a pennyshaped 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.
The general values used in the theoretical equation were r f 0. In this section, we apply the approximate method of kachanov and laures 1989 to the analysis of interacting cracks to obtain the mode i sif of an array of periodic coplanar pennyshaped cracks in an infinite medium subjected to a remote normal traction at infinity, figure 1 a. Previous studies have investigated multiple crack interactions and their effects on the stress intensity factor 5,6,7,8,9,10,11,12,14. The shape of the crack front was also predicted and compared to the crack front. Stress intensity factors for the pennyshaped crack under.
Development of a numerical procedure for mixed mode k. A quasicrystal column containing a pennyshaped crack is investigated to verify the validity of the numerical. The stress intensity factor and energy release rate are determined. Verify the 12 slope loglog and calculate the sif for a plate with a lateral crack b 100mm. The numerical approach is then applied to study interface crack problems in. Stress intensity factors numerical calculations for two. The stress intensity factors are studied for different values of the wave frequency. Crack propagation calculations for optical fibers under. This work shows that it is possible to predict crack growth for a penny shaped surfacecrack subjected to mixedmode loading in general agreement with. It was also found that unbalanced shaft stiffness has a maximum value with a crack at 0. Threedimensional static and dynamic stress intensity factor. Analysis of stress intensity factors for threedimensional interface.
Residual fatigue life estimation using a nonlinear. The effects of both the ratio of crack radius to cylinder radius and the applied magnetoelectrical loads at infinity on the stress intensity factor are evaluated. Stress intensity factors for external and pennyshaped. Numerical results reveal the effects of electric loadings and material parameters of composite on crack propagation and growth. A periodic array of cracks in a functionally graded material under transient mechanical loading is investigated. Od discontinuities in thick cylinders have been the subject of several experimental59 and theoretical120 studies. The pennyshaped crack at a bonded plane with localized.
Handbook of stressintensity factors handbook of stressintensity factors, george c. Threedimensional fracture analysis using tetrahedral. Multiple pennyshaped cracks interaction in a finite bodyand their. Mintegral for calculating intensity factors of cracked. Stress intensity factors in shafts subjected to torsion and axial loading.
Many experimental investigations demonstrate that profuse microcrack propagations can be observed during the evolution of fracture process in the highstress concentration regions. Sif for a pennyshaped crack in a finiteradius cylinder. Finite element analysis of a subsurface pennyshaped crack with. It is the amplitude of stress singularity at the crack tip.
In addition, when tetrahedral enriched elements are used, the crack front region. The stress intensity factor criterion the singular stress field around a crack tip ki is the mode i stress intensity factor. International journal of applied electromagnetics and. R f is the final crack length length at failure, r 0 is the initial crack length, m. A singular integral equation is derived, in which the crack surface displacement is the unknown function. The main objective of this work is to calculate the mixed mode energy release rates around the crack front for a submodel of a global shell m odel by. Thiscritical stress intensity factor isthenameasureof materialtoughness. Problems of fracture mechanics and fatigue a solution.
Threedimensional fracture analysis using tetrahedral enriched. The stress intensity factors for a periodic array of. The stress intensity factor sif at the macrocrack tip and the microcrack tips is obtained. Crack overview fracture mechanics is essential to the mechanical safety of structures, in which cracks and the corresponding stress intensity factors sifs near their tips fronts are important 1. Kachanovs approximate method for crack interactions kachanov, m. A theoretical model for hydraulic fracturing through two. The singular stresses at a point s along the crack front can be written in the local coordinate as 2. Pdf stress intensity factor for a centrally located penny. Effect of fiber volume fraction on fracture mechanics in continuously reinforced fiber composite materials thomas wasik. Dimensionless stress intensity factor for a pennyshaped crack in an in. Threedimensional static and dynamic stress intensity. Functionally graded materials with periodic cracks. Numerical results are obtained to illustrate the variations of the stress intensity factors as a function of the crack. The results seem useful for design of the piezoelectric composite structures and devices of high performance.
The maximum stress intensity factors with respect to a time variable are approximated by a function of biots number and crack length. The special elements also defined as singularity function methods where stress singularity at crack tip is modeled. The case of a pennyshaped crack with arbitrary normal displacements prescribed at its faces inclusion problem was solved exactly and in closed form by. Approximate determination of stress intensity factor for multiple. Stress and reliability analysis of a metalceramic dental. Since an analytical solution exists, the problem of a finite length crack in an infinite body subjected to crack face traction and electric flux. The eightnode plane element, with the software ansys version 18. Sif plays a dominate role because it indicates the singular intensity of linear elastic crack field stress and strain. Ansys engineering analysis system, theoretical manual 1989. According to figure 2, the modeii stress intensity factor at the tip of the fracture shown in figure 1c can be fitted by the following equations. The weight function method is also used to determine closedform stress intensity factors at the deepest point.
Effect of fiber volume fraction on fracture mechanics in. It is observed that the normalized sif at the crack tip b increase with the decrease of the distance between cracks. Alavistress intensity factors for a penny shaped crack in a half space. In this study, a technique is described that combines the use of the ansys. Cracks in transversely isotropic and inhomogeneous elastic.
877 437 35 869 152 863 988 1427 1126 625 632 1527 1280 793 1527 139 71 1033 208 757 779 719 257 1097 453 956 1014 664 1616 847 202 1077 1267 1423 1628 427 795 188 562 478 463 150 283 1431 621 188 355