Numerical results reveal the effects of electric loadings and material parameters of composite on crack propagation and growth. 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. An interaction integral method for computing mixed mode. For example, ansys only provides automatic meshing capability for twodimensional problems. Many experimental investigations demonstrate that profuse microcrack propagations can be observed during the evolution of fracture process in the highstress concentration regions. Thus, the analysis of stresses near the crack tip constitutes an essential part of fracture mechanics. A quasicrystal column containing a pennyshaped crack is investigated to verify the validity of the numerical. Threedimensional static and dynamic stress intensity. Formulation for the mode i stress intensity factor. In this paper, the mintegral is extended for calculating intensity factors for cracked piezoelectric ceramics using the exact boundary conditions on the crack faces. For example, ansys only provides automatic meshing capability for two dimensional problems. 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 poling direction is taken at an angle to the crack faces within the plane. Stress intensity factors for the pennyshaped crack under.
The numerical approach is then applied to study interface crack problems in. 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 obtained results andor conclusions could be of particular interest to the analysis and design of smart sensorsactuators constructed from magnetoelectroelastic composite laminates. The maximum stress intensity factors with respect to a time variable are approximated by a function of biots number and crack length. 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. Fatigue analysis fatigue life of components using ansys high cycle fatigue 6.
A periodic array of cracks in a functionally graded material under transient mechanical loading is investigated. This work shows that it is possible to predict crack growth for a penny shaped surfacecrack subjected to mixedmode loading in general agreement with. Kic is a measure of the fracture toughness of the material. Mintegral for calculating intensity factors of cracked. In 1957, irwin 2 introduced the sifs to describe the stress and displacement. The stress intensity factors for a periodic array of. 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. 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.
Introduction to fracture mechanics david roylance department of materials science and engineering. Effect of fiber volume fraction on fracture mechanics in. The shape of the crack front was also predicted and compared to the crack front. 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.
It illustrates the manner in which the stress intensity factor can be in. Functionally graded materials with periodic cracks. Interacting cracks analysis using finite element method. The problem is solved by the method of boundary integral equations using an iterative algorithm. 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. The pennyshaped crack at a bonded plane with localized. Cracks in transversely isotropic and inhomogeneous elastic. The weight function method is also used to determine closedform stress intensity factors at the deepest point. Analysis of stress intensity factors for threedimensional interface. 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. The stress intensity factors are studied for different values of the wave frequency. Sif plays a dominate role because it indicates the singular intensity of linear elastic crack field stress and strain.
Stress intensity factors in shafts subjected to torsion and axial loading. The shape of crack was mostly defined by local stressdistribution. 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. Development of a numerical procedure for mixed mode k. The rest of the local model is modeled by using 20node three dimensional bri ck elements. It is observed that the normalized sif at the crack tip b increase with the decrease of the distance between cracks.
R f is the final crack length length at failure, r 0 is the initial crack length, m. The stress intensity factor sif at the macrocrack tip and the microcrack tips is obtained. 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. 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. The general values used in the theoretical equation were r f 0. Stress intensity factors for external and penny shaped cracks in transversely isotropic cylinders subjected to thermal shock springerlink. Sif for a pennyshaped crack in a finiteradius cylinder. Stress and reliability analysis of a metalceramic dental. Stress intensity factor for steel hollow pipe with axial crack duration. 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 special elements also defined as singularity function methods where stress singularity at crack tip is modeled. Dimensionless stress intensity factor for a pennyshaped crack in an in.
In this study, a tetrahedral enriched element method and related applications are presented that demonstrate employment of fully unstructured. The effect of crack interactions on stress intensity factors is examined for a periodic array of coplanar penny shaped cracks. 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. Pdf stress intensity factor for a centrally located. Residual fatigue life estimation using a nonlinear. The singular stresses at a point s along the crack front can be written in the local coordinate as 2. The stress intensity factor criterion the singular stress field around a crack tip ki is the mode i stress intensity factor. A theoretical model for hydraulic fracturing through two. Effect of fiber volume fraction on fracture mechanics in continuously reinforced fiber composite materials thomas wasik. Based on linear elastic fracture mechanics, stress intensity factors along the crack front of the pennyshaped crack are evaluated as functions of the crackfront. Thiscritical stress intensity factor isthenameasureof materialtoughness.
Materials free fulltext determination of the enhancement or. Because of the importance of sif, its solutions have b een paid very high attention since the beginning of fracture mechanics. The commonly used methods are based on the stress intensity factor dominated power law. Ansys is used for the initial elastic stress analysis of the uncracked model, and.
The bulk of the experimental work57,9 has dealt with the effects of stress concentrations of fatigue crack. 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. In this study, a technique is described that combines the use of the ansys. Handbook of stressintensity factors handbook of stressintensity factors, george c. In addition, when tetrahedral enriched elements are used, the crack front region. The paper concerns the contact interaction of the opposite faces of a pennyshaped crack under harmonic wave at oblique incidence. Kassir mk, bregman am 1972 the stress intensity factor for a pennyshaped crack between two. Ksolutions and fatigue crack growth in fcc single crystal superalloys by srikant r lii in. Finite element analysis of a subsurface pennyshaped crack with. Threedimensional fracture analysis using tetrahedral. The finite element model is established by ansys 12. Pdf stress intensity factor for a centrally located penny. Kachanovs approximate method for crack interactions kachanov, m.
Threedimensional static and dynamic stress intensity factor. Fracture problems with ansys a very brief introduction. Verify the 12 slope loglog and calculate the sif for a plate with a lateral crack b 100mm. Threedimensional fracture analysis using tetrahedral enriched. Ansys engineering analysis system, theoretical manual 1989.
A singular integral equation is derived, in which the crack surface displacement is the unknown function. In the numerical procedure normalized complex stress intensity. It is the amplitude of stress singularity at the crack tip. Durability evaiuiiti0n 6tceramic components using careslife. These include an inclined through crack, an embedded pennyshaped crack, an. The eightnode plane element, with the software ansys version 18. Analysis of stress intensity factors for threedimensional. 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.
Stress intensity factors numerical calculations for two. Stress intensity factors for external and pennyshaped. 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. 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. The finite element meshes for these models are generated using ansys. 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. The condition of continuity between elements is the most important. Stress intensity factor for a centrally located penny shaped crack. The stress intensity factor and energy release rate are determined.
Approximate determination of stress intensity factor for multiple. International journal of applied electromagnetics and. Od discontinuities in thick cylinders have been the subject of several experimental59 and theoretical120 studies. 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. Previous studies have investigated multiple crack interactions and their effects on the stress intensity factor 5,6,7,8,9,10,11,12,14. Pdf more accurate stress intensity factor derived by.
The stress intensity factor for the crack face pressure case is also shown in figure 5a. The case of a pennyshaped crack with arbitrary normal displacements prescribed at its faces inclusion problem was solved exactly and in closed form by. Lehigh university, institute of fracture and solid mechanics, 1973. Numerical results are obtained to illustrate the variations of the stress intensity factors as a function of the crack.
A penny shaped crack on the 111 plane evolved into a hexagonal crack front, because. The results seem useful for design of the piezoelectric composite structures and devices of high performance. The stress intensity factor at the tip of a pennyshaped crack of radius in an infinite domain under uniaxial tension is k i 2. Multiple pennyshaped cracks interaction in a finite bodyand their.
409 1475 902 1366 563 98 1006 202 536 137 510 200 661 385 846 193 568 756 887 1517 782 76 1544 1662 1277 171 143 1655 1324 1621 1489 327 918 1375 622 936 1158 1101 1304 1010 1476 778 183