Strip model plasticity analysis of mixed mode crack opening displacement in aluminum alloy LY12

Under mixed mode loading, the crack tip blunts and undergoes displacements in two directions, the normal and shear component corresponding to Mode I and Mode II loading, respectively. These local displacements are determined by the duplicated film method and used to analyze the behavior of mixed mode fracture in aluminum alloy LY12. The mixed mode resultant crack opening displacement (COD) at fracture initiation tends to increase more rapidly with increase of the Mode II shear component. The fracture initiation value of COD for pure Mode II loading is six (6) times greater than that for Mode 1 loading. The same applies to the maximum effective plastic strain crack growth near the crack front. Observed are two typical morphologies, the equal-axes dimples and the parabolic dimples with evidence of slippage as dominated, respectively, by Mode I and Mode II loading.     

A centrally cracked thin circular disk, Part I: 3D Elastic–plastic finite element analysis

A full field solution, based on small deformation, three-dimensional elastic–plastic finite element analysis of the centrally cracked thin disk under mode I loading has been performed. The solution for the stresses under small-scale yielding and lo!cally fully plastic state has been compared with the HRR plane stress solution. At the outside of the 3D zone, within a distance of rσ_o/J=18, HRR dominance is maintained in the presence of a significant amount of compressive stress along the crack flanks. Ahead of this region, the HRR field overestimate the stresses. These results demonstrate a completely reversed state of stress in the near crack front compared to that in the plane strain case. The combined effect of geometry and finite thickness of the specimen on elastic–plastic crack tip stress field has been explored. To the best of our knowledge, such an attempt in the published literature has not been made yet. For the qualitative assessment of the results some of the field parameters have been compared to the available experimental results of K, gives a fair estimate of the crack opening stress near the crack front at a distance of order 10⁻² in. On the basis of this analysis, the Linear Elastic Fracture Mechanics approach has been adopted in analyzing the fatigue crack extension experiments performed in the disk (Part II).     

Mixed mode (I and II) crack tip fields in bulk metallic glasses

Stationary crack tip fields in bulk metallic glasses under mixed mode (I and II) loading are studied through detailed finite element simulations assuming plane strain, small scale yielding conditions. The influence of internal friction or pressure sensitivity on the plastic zones, notch deformation, stress and plastic strain fields is examined for different mode mixities. Under mixed mode loading, the notch deforms into a shape such that one part of its surface sharpens while the other part blunts. Increase in mode II component of loading dramatically enhances the normalized plastic zone size, lowers the stresses but significantly elevates the plastic strain levels near the notch tip. Higher internal friction reduces the peak tangential stress but increases the plastic strain and stretching near the blunted part of the notch. The simulated shear bands are straight and extend over a long distance ahead of the notch tip under mode II dominant loading. The possible variations of fracture toughness with mode mixity corresponding to failure by brittle micro-cracking and ductile shear banding are predicted employing two simple fracture criteria. The salient results from finite element simulations are validated by comparison with those from mixed mode (I and II) fracture experiments on a Zr-based bulk metallic glass.     

Elastic–plastic near field solution of an eccentric crack under shear in a finite width plate

The near crack line analysis method is used to investigate an eccentric crack loaded by shear forces in a finite width plate, and the analytical solution is obtained in this paper. The solution includes: the unit normal vector of the elastic–plastic boundary near the crack line, the elastic–plastic stress fields near crack line, variations of the length of the plastic zone along the crack line with an external loads, and the bearing capacity of a finite plate with a centric crack loaded by shear stress in the far field. The results obtained in this paper are sufficiently precise near the crack line because the assumptions of small scale yielding theory have not been made and no other assumptions have been taken. Subsequently, the present results are compared with the traditional line elastic fracture mechanical solutions and elastoplastic near field solutions under small scale yielding condition. On the basis of the minimum strain energy density (SED) theory, the minimum values of SED in the vicinity of the crack tip are determined, the initial growth orientation of crack are determined. It is found that the normalized load under large scale yielding condition is higher than those under small scale yielding condition when the length of the plastic zone is the same.     

Perfectly plastic stress field at a stationary crack tip

Under the hypothesis that all the perfectly plastic stress components at a orach tip are the functions of θ only, making use of yield conditions and equilibrium equations. we derive the generally analytical expressions of the perfectly plastic stress field at a crack tip. Applying these generally analytical expressions to the concrete cracks, the analytical expressions of perfectly plastic stress fields at the tips of Mode Ⅰ Mode Ⅱ, Mode Ⅲ and Mixed Mode Ⅰ-Ⅱ cracks are obtained.     

Experimental analysis of crack tip fields in rubber materials under large deformation

A three-nested-deformation model is proposed to describe crack-tip fields in rubber-like materials with large deformation.The model is inspired by the distribution of the measured in-plane and out-of-plane deformation.The inplane displacement of crack-tip fields under both Mode I and mixed-mode(Mode I-II) fracture conditions is measured by using the digital Moire’ method.The deformation characteristics and experimental sector division mode are investigated by comparing the measured displacement fields under different fracture modes.The out-of-plane displacement field near the crack tip is measured using the three-dimensional digital speckle correlation method.     

An inclined surface crack subject to biaxial loading

The elastic–plastic stress fields and mode mixity parameters for semi-elliptical surface cracks on biaxial loaded plates have been investigated using detailed three-dimensional finite element calculations. Different degrees of mode mixity are given by combinations of the far-field stress level, biaxial stress ratio and inclined crack angle. These analyses were performed for different surface flaw geometries to study the combined load biaxiality and mode mixity effects on the crack-front stress fields and the size and shape of the plastic zones. It is clear from considering the local stress distributions along the crack front that the elastic crack tip singularities have been derived for several particular cases of mixed mode biaxial loading. By theoretical analysis, the new formulae have been introduced for both the elastic and plastic mode-mixity parameters, accounting for ratios between the I/II, II/III and III/I modes. Particular attention was paid to the strong variations of the mode-mixity parameters along the semi-elliptical surface crack front. The mixed-mode behavior of the crack growth direction angle along the semi-elliptical crack front for different combinations of biaxial loading and inclination crack angles was also determined. It was done using methods based on the maximum tangential stress and the strain energy density criteria.     

Mixed-mode interface toughness of wafer-level Cu-Cu bonds using asymmetric chevron test

Characterization of interfacial adhesion is critical for the development of wafer bonding processes to manufacture microsystems with high yield and reliability. It is imperative that the test method used in such adhesion studies corresponds to the loading conditions present during processing and operation of the devices. In most applications in which wafers and die are bonded, the interface experiences a combination of shear and normal loading (i.e. mixed-mode loading) with the relative magnitude of the Mode I and II components varying in different scenarios. In the current work, the toughness of Cu-Cu thermocompression bonds, which are of interest for the fabrication of three-dimensional integrated circuits, is analyzed using a bonded chevron specimen with layers of different thickness that allows for the application of interfacial loading with variable mode mixity. The phase angle (a function of the degree of mode mixity at the interface) is varied from 0° to 24° by changing the layer thickness ratio from 1 to 0.48. The Cu-Cu bond toughness increases from 2.68 to 10.1 J/m², as the loading is changed from Mode I (pure tension) to a loading with a phase angle of 24°. The energy of plastic dissipation increases with increasing mode mixity, resulting in the enhanced interface toughness. The Mode I toughness of Cu-Cu bonds is minimally affected by plasticity, and therefore, provides the closest estimate of the interfacial work of fracture under the bonding conditions employed.     

层状磁电复合材料界面共线裂纹平面问题分析

本文研究了面内电磁势载荷作用下双层压电压磁复合材料中共线界面裂纹问题．考虑了压电材料的导磁性质和压磁材料的介电性质,引入了界面电位移和磁感强度的连续性条件．利用Fourier 变换得到一组第二类Cauchy 型奇异积分方程．进一步导出了相应问题的应力强度因子、电位移强度因子和磁感强度强度因子的表达式,给出了应力强度因子的数值结果．结果表明电磁载荷会导致界面裂纹尖端I、II 混合型应力奇异性,同时还伴随着电位移和磁感强度的奇异性．比较了双裂纹左右端的应力强度因子,发现在面内极化方向上施加面内磁势载荷时共线裂纹内侧尖端区域的两个法向应力场发生互相干涉增强．     

Modeling of fatigue crack growth from a notch

When a fatigue crack is nucleated and propagates into the vicinity of the notch, the crack growth rate is generally higher than that can be expected by using the stress intensity factor concept. The current study attempted to describe the crack growth at notches quantitatively with a detailed consideration of the cyclic plasticity of the material. An elastic–plastic finite element analysis was conducted to obtain the stress and strain histories of the notched component. A single multiaxial fatigue criterion was used to determine the crack initiation from the notch and the subsequent crack growth. Round compact specimens made of 1070 steel were subjected to Mode I cyclic loading with different R-ratios at room temperature. The approach developed was able to quantitatively capture the crack growth behavior near the notch. When the R-ratio was positive, the crack growth near a notch was mainly influenced by the plasticity created by the notch and the resulted fatigue damage during crack initiation. When the R-ratio was negative, the contact of the cracked surfaces during a part of a loading cycle reduced the cyclic plasticity of the material near the crack tip. The combined effect of notch plasticity and possible contact of cracked surface were responsible for the observed crack growth phenomenon near a notch.     

A finite element study of the near crack tip deformation of a ductile material under mixed mode loading

In ductile fracture, voids near a crack tip play an important role. From this point of view, a large deformation finite element analysis has been made to study the deformation, stress and strain, and void ratio near the crack tip under mixed mode plane strain loading conditions, employing Gurson's constitutive equation which has taken into account the effects of void nucleation and growth. The results show that: (i) one corner of the crack tip sharpens while the other corner blunts, (ii) the stress and strain distributions except for the near crack tip region, can be superimposed by normalizing distance from the crack tip by a crack tip deformation length, i.e., a steady-state solution under a mixed mode condition has been obtained, (iii) the field near a crack tip can be divided into four characteristic fields (K field, HRR field, blunted crack tip field, and damaged region), and (iv) the strain and void volume fraction become concentrated in the sharpened part of a crack tip with increasing Mode II component.     

Determination of mixed-mode stress intensity factors, fracture toughness, and crack turning angle for anisotropic foam material

A numerical and experimental investigation for determining mixed-mode stress intensity factors, fracture toughness, and crack turning angle for BX-265 foam insulation material, used by NASA to insulate the external tank (ET) for the space shuttle, is presented. BX-265 foam is a type of spray-on foam insulation (SOFI), similar to the material used to insulate attics in residential construction. This cellular material is a good insulator and is very lightweight. Breakup of segments of this foam insulation on the shuttle ET impacting the shuttle thermal protection tiles during liftoff is believed to have caused the space shuttle Columbia failure during re-entry. NASA engineers are interested in understanding the processes that govern the breakup/fracture of this material from the shuttle ET. The foam is anisotropic in nature and the required stress and fracture mechanics analysis must include the effects of the direction dependence on material properties. Material testing at NASA Marshall Space Flight Center (MSFC) has indicated that the foam can be modeled as a transversely isotropic material. As a first step toward understanding the fracture mechanics of this material, we present a general theoretical and numerical framework for computing stress intensity factors (SIFs), under mixed-mode loading conditions, taking into account the material anisotropy. We present SIFs for middle tension – M(T) – test specimens, using 3D finite element stress analysis (ANSYS) and FRANC3D fracture analysis software. SIF values are presented for a range of foam material orientations. Mode I fracture toughness of the material is determined based on the SIF value at failure load. We also present crack turning angles for anisotropic foam material under mixed-mode loading. The results represent a quantitative basis for evaluating the strength and fracture properties of anisotropic foam insulation material.     

Effect of temperature on the fracture modes of BHW-35 steel under mode-II loading

Finite element method (FEM) has been used to analyze the stress and strain fields and the stress tri-axial levels around the tip of the crack under mode- II loading. The results show that: under mode- II loading, the direction of the maximum tensile stress and that of the maximum tri-axial levels (R _o ) exist at an angle of –75. 3° from the original crack plane; the maximum shear stress andR _o = 0 exist along the original crack plane.Mode- II loading experiment using BHW-35 steel at different temperatures show that there are two kinds of fracture mode, opening mode (or tensile mode) and sliding mode (or shear mode). A decrease in temperature causes the fracture mode to change from shear mode to tensile mode. For BHW-35 steel, this critical temperature is about –90 C. Actually, under any kind of loading mode (mode I . mode II , mode III or mixed mode), there always exist several kinds of potenital fracture modes (for example, opening mode, sliding mode, tearing mode or mixed mode). The effect of temperature under mode- II loading is actually related to the change of the elastic-plastic properties of the material.     

The stress field near a blunting crack tip under mixed modes 1 and 2

The distribution of stress near a crack tip deforming in plane strain under mixed Modes 1 and 2 loads is calculated using slip-line field theory. Large geometry change effects are taken into account. The results are compared and contrasted both with known Mode 1 large-deformation solutions and with previous mixed-mode solutions for a sharp crack. The results are then combined with stress-based fracture criteria to discuss mixed-mode fracture on the lower shelf.     

Analytic solutions to crack tip plastic zone under various loading conditions

Based on stress field equations and Hill yield criterion, the crack tip plastic zone is determined for orthotropic materials and isotropic materials under small-scale yielding condition. An analytical solution to calculating the crack tip plastic zone in plane stress states is presented. The shape and size of the plastic zone are analyzed under different loading conditions. The obtained results show that the crack tip plastic zones present “butterfly-like” shapes, and the elastic–plastic boundary is smooth. The size of the plastic zone for orthotropic composites is less at the crack tip for various loading conditions, compared with the case of isotropic materials. Crack inclination angle and loading conditions affect greatly the size and shape of crack tip plastic zone. The mode I crack has a crucial effect on the plastic zone for mixed mode case in plane stress state. The plastic zone for pure mode I crack and pure mode II crack have a symmetrical distribution to the initial crack plane.     

Analytical representation of the non-square-root singular stress field at a finite angle sharp notch

The stress field near the tip of a finite angle sharp notch is singular. However, unlike a crack, the order of the singularity at the notch tip is less than one-half. Under tensile loading, such a singularity is characterized by a generalized stress intensity factor which is analogous to the mode I stress intensity factor used in fracture mechanics, but which has order less than one-half. By using a cohesive zone model for a notional crack emanating from the notch tip, we relate the critical value of the generalized stress intensity factor to the fracture toughness. The results show that this relation depends not only on the notch angle, but also on the maximum stress of the cohesive zone model. As expected the dependence on that maximum stress vanishes as the notch angle approaches zero. The results of this analysis compare very well with a numerical (finite element) analysis in the literature. For mixed-mode loading the limits of applicability of using a mode I failure criterion are explored.     

Plane-strain crack-tip stress field for anisotropic perfectly-plastic materials

Plane-strain crack-tip stress solutions for anisotropic perfectly-plastic materials are presented. These solutions are obtained using the plane-strain slip-line theory developed by Rice (1973). The plastic anisosotropy is described by the Hill quadratic yield condition. The crack-tip stress solutions under symmetric (Mode I) and anti-symmetric (Mode II) conditions agree well with the low-hardening solutions for the corresponding power-law hardening materials. The crack-tip stress solutions under mixed Mode I and II conditions are also presented. All the solutions indicate that the general features of the slip-line field near a crack tip in orthotropic plastic materials with the elliptical yield contours in the Mohr plane are the same as those associated with isotropic plastic materials. However, the angular variations of the crack-tip stress fields for the materials with large plastic orthotropy differ substantially from those for isotropic plastic materials. Modifications due to polygonal yield contours are outlined and implications of solutions to the fracture analysis of ductile composite materials containing macroscopic flaws are discussed.     

Analysis of dynamic mixed-mode isochromatics

The mixed-mode, elastodynamic state of stress in the neighborhood of a constant-velocity crack tip is used to generate numerically unsymmetric isochromatics. Unsymmetry associated with the third-order terms of a mixed-mode stress field, with and without the Mode II singular stress term, is also investigated. In extractingK _I from an unsymmetric isochromatic pattern, errors in the Mode I fracture parameters due to the assumed presence ofK _II in aK _I stress field were found to be significant when data are taken more than 4 mm from the crack tip. Paper was presented at V International Congress on Experimental Mechanics held in Montreal, Quebec, Canada on June 10–15, 1984.     

Near crack line elastic–plastic analysis for a infinite plate loaded by two pairs of point tensile forces

The near crack line analysis method is used to investigate a center crack loaded by two pairs of point tensile forces in an infinite plate in an elastic–perfectly plastic solid, and the analytical solutions are obtained in this paper. These solutions include: the elastic–plastic stress field near the crack line, the law that the length of the plastic zone along the crack line is varied with an external loads and the bearing capacity of an infinite plate with a center crack. The results of this paper are sufficiently precise near the crack line because the assumptions of the small scale yielding theory have not been used and no other assumptions have been taken.     

Mixed-mode I/II fields around a crack with a cohesive zone ahead of the crack tip

Considering a cohesive zone ahead of the crack tip, the linear elastic crack problem under Mode I, Mode II or mixed-mode conditions is formulated in an elliptic coordinate system, so that the cohesive surfaces are conveniently represented by straight line segments. It is shown that the displacement and stress fields around the crack tip and the cohesive zone, expressed in terms of elliptic coordinates, have a simple mathematical form, which does not contain a stress singularity at the crack tip due to the existence of the cohesive zone.