A finite element investigation of stable crack growth in anti-plane shear

Although the anti-plane strain case is of minor practical value in engineering applications, such an idealization facilitates mathematical investigations of strain and displacement fields accompanying extending cracks. This paper presents finite element solutions to anti-plane strain crack propagation problems and contrasts the numerical results with available analytic solutions in an effort to assess the accuracy of the numerical procedures. The nature of dominant strain singularities for stationary and moving cracks, the question of stableversus unstable or catastrophic crack growth and the implications of various proposed fracture criteria are discussed.

---

Résumé Bien que le case des dilatations antiplanaires soit d'une importance pratique mineure dans les applications de la construction, une telle idéalisation facilite les investigations mathématiques sur les champs de dilatation et de déplacement qui accompagnent des fissures en cours d'extension. Le mémoire présente des solutions par élément fini aux problèmes de la propagation de fissure sous des dilatations antiplanaires et fait apparaître le contraste entre les résultats numériques et les solutions analytiques disponibles, dans une tentative de faire valoir ou de constater la précision des procédures numériques. La nature des singularités déterminantes de dilatation dans le cas de fissures stationnaires et de fissures en mouvement, la question de la croissance catastrophique instable ou stable d'une fissure et les implications que les divers critères proposés de rupture contiennent sont discutées.

     

A special bending crack tip finite element

The general expression of the Williams expansion around a bending crack tip is presented. Based on this a special circular bending crack tip finite element of an arbitrary high order is developed, and the element stiffness matrix is explicitly given, which is especially convenient for engineering analyses. Numerical examples are presented and compared with previous results to demonstrate the efficiency and accuracy of the special element.     

Three-dimensional finite element simulations of mixed-mode stable tearing crack growth experiments

Mixed-mode stable tearing crack growth events in Arcan plate specimens made of aluminum alloy 2024-T3 are simulated using three-dimensional (3D) finite element methods. A modeling/simulation procedure utilizing a mixed-mode CTOD fracture criterion and the custom 3D crack growth simulation software, CRACK3D, with an automatic local re-meshing option is demonstrated. Simulation predictions of the load-crack extension curve and the in-plane curvilinear crack growth path are compared with experimental measurements for various mixed-mode loading cases. Issues such as the effects of near-tip finite element size and crack extension increment size on simulation predictions are investigated.     

A finite element analysis of stable crack growth in an aluminium alloy

Extensive stable cracking has been observed in large test pieces of 25 mm thick weldable AlMgZn alloy which is used in the construction of a portable bridge. Standard fracture specimens produced valid K_IC values, with short cracks exhibiting unstable fracture. Finite element analysis of the large specimens determined a valid J_R-curve that can increase the effective K_C by several times the K_IC value. The R-curve has an unusual ‘concave’ shape that is associated with the change from initially flat fracture to fully slant fracture. The early stages of the R-curve are affected by in-plane constraint that can be indexed by the T-stress. The R-curve can be used to explain the stability of long cracks in full-scale tests on a bridge prototype, compared with the instability of short cracks in small, standard test pieces.     

A finite element study of crack growth in WC-Co

A finite element analysis of crack growth in tungsten carbide cobalt has been used to study the plastic deformation of binder ligaments bridging the crack faces in the wake of a matrix crack. A multiligament zone observed on the intersection of an arrested crack with a free surface is used to model the crack tip region. The plane stress results demonstrate that plasticity is confined to a band linking the ligament tips well in accordance with experimentally found deformation patterns. The plane strain calculations for the same microstructure supply information about hole nucleotion and growth which are known to control the failure process of the ligaments. It is concluded from a recent analysis of void growth in homogeneous materials, that plastic deformation in the binder of WC-Co is concentrated in the neck between a blunting crack tip and a void growing ahead of it. Thus in both cases, plane stress and plane strain, non bridging binder regions deform purely elastically in contrast to the results of recent finite element calculations. It is seen that the previously used concept of a plastic zone size in the binder of cemented carbides equal to or larger than the mean intercept length of the binder, 305-1, must be modified. 305-2 constitutes only an upper limit for the mean size of the plastic zone while the actual extension of plasticity is smaller.

---

Résumé On utilise une analyse par éléments finis de la croissance d'une fissure dans un alliage Cobalt-Carbure de Tungstène pour l'étude de la déformation plastique des ligaments de liaison entre les faces d'une fissure prenant naissance dans une matrice. Pour représenter la région à l'extrémité de la fissure, on utilise la zone multiligamentaire observée à l'intersection d'une fissure ouverte et d'une surface libre. Les résultats correspondent à un état plan de tension démontrant que la plasticité est confinée à une bande reliant les extrémités des ligaments, ce qui est bien en accord avec les aspects de déformation trouvés expérimentalement. Les calculs en état plan de déformation pour la même microstructure fournissent une information sur le processus de nucléation et de croissance des lacunes, qui est connu pour contrôler le processus de rupture d'un ligament. On conclut de l'analyse pour des matériaux homogènes que la déformation plastique dans une liaison de WC-Co est concentrée dans la position rétrécie entre l'extrémité arrondie d'une fissure et une lacune en croissance en amont de celle-ci. Dès lors, dans les deux cas de tension plane ou de déformation plane, les régions de liaison qui ne relient pas les faces de la rupture se déforment de manière purement élastique, ceci en contraste avec les résultats de calculs récents par éléments finis. On constate que le concept de taille de zone plastique égale ou supérieure à la longueur moyenne 317-3 de la liaison doit être modifié. 317-4 ne représente qu'une limite supérieure pour la taille moyenne de la zone plastique, tandis que l'étendue réelle de la plasticité est plus petite.

     

Numerical study of dynamic crack growth by the finite element method

Recent developments in numerical techniques for dynamic transient stress analysis have ensured that realistic models can now be employed in crack propagation studies. In this paper transient dynamic finite element solutions are undertaken for both double cantilever beam (DCB) and pipeline problems with propagation of the crack being permitted. Standard parabolic isoparametric elements are employed for spatial discretization with an explicit (central difference) scheme being employed for time integration. Both critical stress and energy balance crack propagation criteria are considered.The pressurised pipeline problem is solved for as a fully three-dimensional solid. Firstly, a stationary crack is considered and both large deformations and plasticity effects are accounted for. The transient case of a dynamically propagating crack is then modelled, employing both a stress and energy criterion. Elastic large deformation behaviour is permitted for this case.

---

Résumé Des développements récents dans les techniques numériques pour l'analyse des contraintes dynamiques transitoires ont permis d'utiliser à présent des modèles réalistes dans les études de propagation des fissures. Dans ce mémoire, on envisage des solutions par éléments finis pour les transitoires dynamiques dans les cas de la poutre double cantilever et de problèmes de pipelines où l'on autorise la propagation d'une fissure. On recourt aux éléments paramétriques paraboliques standards pour réaliser une division discrète de l'espace, et l'on utilise pour l'intégration dans le temps un schéma explicite à différence centrale. On considére à la fois les critères de contraintes critiques et d'équilibre d'énergie lors de la propagation de la fissure. Le problème du pipeline pressurisé est solutionné en considérant ce dernier comme un solide tridimensionnel. En premier lieu, on considère une fissure stationnaire et l'on tient compte des effets des grandes déformations et de la plasticité. On met ensuite en équation le cas transitoire d'une fissure en propagation dynamique, en utilisant un critère de contrainte et un critère d'énergie. Ce cas permit d'envisager le comportement sous des déformations élastiques importantes.

     

A finite element algorithm for creep crack growth

A finite element algorithm involving the “Breakable element” concept is proposed for the prediction of the growth of a crack in a solid subject to combined thennoelastic-plastio-creep load. The unique advantage of this algorithm is its ability to provide detail stress and strain distributions, the kinematics of the inelastic zones, as well as the profiles of the growing crack. A numerical example with three assigned effective rupture strains as fracture criteria, is included to illustrate these special features.     

Extended finite element method for cohesive crack growth

The extended finite element method allows one to model displacement discontinuities which do not conform to interelement surfaces. This method is applied to modeling growth of arbitrary cohesive cracks. The growth of the cohesive zone is governed by requiring the stress intensity factors at the tip of the cohesive zone to vanish. This energetic approach avoids the evaluation of stresses at the mathematical tip of the crack. The effectiveness of the proposed approach is demonstrated by simulations of cohesive crack growth in concrete.     

Experimental and finite element studies on mode I and mixed mode (I and II) stable crack growth—II. finite element analysis

Finite element studies are presented on both mode I and mixed mode stable crack growth under static loadings through an aluminium (D16AT) alloy. A COD based criterion has been used to predict the load-displacement diagram from initiation to instability. The theoretical predictions are compared with experimental results presented in Part I. Results on computed crack profiles, stress-strain distribution ahead of the crack tip, J integrals, J resistance curves, plastic zones, etc., are included. The study indicates that the load-displacement diagram associated with a mixed mode stable crack growth in a compact tension type of specimen geometry can be predicted reasonably accurately using the criterion of a fixed crack opening displacement at a finite distance behind the crack tip provided the crack is allowed to grow in the direction of initial growth in the finite element analysis. The crack assumes a more blunted profile in a mixed mode than in the mode I at all the stages of stable extension. The distributions of normal stress and strain in the direction perpendicular to the crack extension line, ahead of the current crack tip, have similarities between the mode I and mixed mode, irrespective of loading angle. Both the stress and strain levels increase as the crack extension proceeds. In a mixed mode, the J integral at the onset of crack extension is the lowest compared with the values at the later stages of the extension. Further, the tearing modulus associated with initial kinking is very small; it becomes close to the mode I values at the later stages. The tearing modulus remained approximately constant during the whole mode I stable growth and it had a similar trend subsequent to kinking in a mixed mode. The specific work of crack extension is zero as Δa → 0 and it increases gradually with Δa irrespective of the mode of loading; the actual variation depends on the loading angle. The plastic zone size grows as the stable extension progresses; the growth is approximately the maximum along the crack extension line.     

A study of stable crack growth using experimental methods,finite elements and fractography

Systematic analysis of the in-plane constraint influence on J-resistance curves is presented. J_R curves were also recorded and analyzed beyond the limits of crack extension inside which the stress field can be assumed to be dominated by J-integral. Three steels and four types of specimen: SEN(B), SEN(T), CCT and DENT have been tested. Along with the J_R curves the fracture mechanisms have been analyzed with the help of scanning microscopy. The numerical, finite element analysis has been adopted to compute the Q-stresses, as a measure of the in-plane constraint prior to the onset of crack growth. The analysis of the stress field in front of the crack has been performed to check whether the state of stress prior to the crack growth can predetermine the way the crack will grow. It turns out that characteristic features in the J_R curves runs can be predicted qualitatively from the Q(a/W) and Q(J) curves. However, there is a good correlation between Q-stress and voids diameters on fractured surfaces. Several patterns in J_R curves runs have been observed for tested specimens; e.g. no influence of specimen thickness on J_R curves runs was observed for side-grooved specimens. Strong influence of specimen thickness on J_R curve shape was observed for non-side-grooved specimens. J_R curve run higher for thinner specimens unless they are dominated by plane stress. For bent specimens J_R curves run higher for shorter cracks but they run lower for specimens in tension.     

Experimental and finite element analyses on stress intensity factors of an elliptical surface crack in a circular shaft under tension and bending

Experimental backtracking technique and finite element analysis have been employed to evaluate the stress intensities along the front of an elliptical surface crack in a cylindrical rod. The finite element solution covers a wide range of crack shapes loaded under end-free and end-constrained axial tension and pure bending. Convenient closed form stress intensity expressions along the whole crack front for each of the loading cases have been given in terms of the crack aspect ratio, crack depth ratio and place ratio.The closed form solutions have been compared against a number of representative solutions collected from the literature. It has been found that different finite element results for the interior points are generally in good mutual agreement, while solutions derived from other methods may sometimes indicate different trends. At the surface interception point agreement is less good because of a complication in the interpretation of stress intensity there.Experimental backtracking results on the end-constrained axial tension case corroborate well with the closed form solution presented. It suggests that the current closed form solution is adequate in describing the stress intensities along the whole crack front of real surface cracks in cylindrical rods.     

A finite element method for crack growth without remeshing

An improvement of a new technique for modelling cracks in the finite element framework is presented. A standard displacement‐based approximation is enriched near a crack by incorporating both discontinuous fields and the near tip asymptotic fields through a partition of unity method. A methodology that constructs the enriched approximation from the interaction of the crack geometry with the mesh is developed. This technique allows the entire crack to be represented independently of the mesh, and so remeshing is not necessary to model crack growth. Numerical experiments are provided to demonstrate the utility and robustness of the proposed technique. Copyright © 1999 John Wiley & Sons, Ltd.     

A new enriched finite element for fatigue crack growth

This paper presents a numerical method for fatigue crack growth in the framework of finite element method, i.e. a new enriched element is presented in which only the analytical solutions around crack tips are used to describe the displacements and stresses fields. A special variational principle is introduced to simplify the mathematical derivations for discrete equation, and the stiffness matrix of the new enriched element is given in a compact form. Moreover, the stiffness matrix is found to be independent on the element size. Numerical examples on fatigue crack growth are given to illustrate the validity of the present method.     

Prediction of fatigue crack growth rates using crack closure finite element analysis

A three-dimensional finite element fatigue crack closure model of a corner crack and of a through thickness crack has been developed to evaluate the range of effective stress intensity factor from the distribution of the range of stress ahead of the crack tip. The corresponding fatigue crack growth rate was evaluated from a Paris law fit to experimental data from high stress ratio tests. The point of origin for the range of stress distribution was adjusted in accordance with Irwin’s plastic zone correction. Encouraging comparisons of finite element predictions of fatigue crack growth rate incorporating closure effects with experimental measurements were obtained.     

Experimental and finite element studies on mode I and mixed mode (I and II) stable crack growth—I. Experimental

Experimental results on mode I and mixed mode stable crack growth under static loadings through an aluminium alloy (D16AT) are presented. The compact tension type of geometry was employed for both the sets of tests. Data pertaining to load-deflection diagrams, crack opening displacements, crack front geometry, etc., are included. There is a greater spurt of crack growth at the initiation stage in a mixed mode than in mode I. The crack opening angle (COA) remained nearly constant during the whole stable growth. There is a substantial tunneling, the extent of which increases as the extension progresses in both mode I and mixed mode. The tunneling reduces as the ratio a₀/W increases. Because of this tunneling, the COD at a point finite distance behind the crack tip and on the specimen surface is much more than expected. At the maximum load the tunneling is 2 to 3.5 mm in the case of mode I. The crack extends initially almost along a straight line at an angle with the initial crack in a mixed mode. The maximum to initiation load ratio varied in the range 1.50 to 1.75 for the whole range of tests.     

Crack growth in a mixed-mode loading on marble beams under three point bending

A combined theoretical and experimental study of a crack growth in a mixed-mode I–II loading is presented. A 160×40×20 mm marble beam, with an artificial crack 8 mm and 10 mm long each, was subjected to three point bending. The crack was located vertically to the beam's lower longitudinal fiber, through the whole width of the beam. The position of the crack was displaced from the center of the beam to one of the supporting points. The vertical force P, placed on the middle of the upper fiber of the beam, imposed the combination of the opening (mode I) and the sliding (mode II) modes on the crack mouth, creating the mixed-mode I–II loading case. The stress intensity factors K _I and K _II, which describe the local stress and strain field around the crack tip, were determined by a suitable finite element program. The crack growth was defined by two classical fracture criteria of LEFM; the minimum strain energy density and the maximum circumferential stress criteria. The initial crack growth angle () was calculated from both criteria and the critical load (P _c) from the minimum strain energy density (SED) criterion. These theoretical predictions were compared with some experimental results found from three marbles with different elastic constants; the Krystallina of Kavala, the Snow-white of Thassos and the White of Piges Drama. The theoretical results showed the same trend of and P _c as the experimental ones and they are in good agreement.Presented at Fourth Greek National Congress on Mechanics, 26–29 June 1995, held at Xanthi, Greece.     

Micromechanical multiplane finite element modeling of crack growth in fiber reinforced materials

The purpose of this study was to devise and verify a scheme of analysis which can be used to investigate the micromechanical failure mechanisms and determine an effective fracture toughness for a class of fiber reinforced materials. The material of primary interest in this study consists of a linearly elastic matrix material reinforced with rows of parallel, linearly elastic and straight fibers. Micromechanical multiplane finite element and experimental studies of the stress conditions near a crack front in a side cracked fiber reinforced epoxy tensile specimen were conducted. The 2-D multiplane method of analysis, recently developed at Syracuse University for approximate analysis of a class of 3-D problems, was the basis of the micromechanical finite element analytical technique developed in this study. Since failure of a member fabricated from a fiber reinforced material is generally proceeded by local failures, sequential finite element analyses were performed to model the progressive failure mechanism. Local failure modes considered in the analysis are yield in either the matrix material or fibers, crack extension in the matrix material, and failure of the matrix to fiber bond. The agreement between the multiplane analytical and laboratory test results show that the multiplane method provides a useful tool for micromechanical study of fiber reinforced composite materials.     

Modelling crack growth by level sets in the extended finite element method

An algorithm which couples the level set method (LSM) with the extended finite element method (X‐FEM) to model crack growth is described. The level set method is used to represent the crack location, including the location of crack tips. The extended finite element method is used to compute the stress and displacement fields necessary for determining the rate of crack growth. This combined method requires no remeshing as the crack progresses, making the algorithm very efficient. The combination of these methods has a tremendous potential for a wide range of applications. Numerical examples are presented to demonstrate the accuracy of the combined methods. Copyright © 2001 John Wiley & Sons, Ltd.