Finite Element Modeling of Crack Tip Blunting for Estimation of Energy Release Rate Component of Mode I Crack near a Strength Mismatched Interface

The paper presents finite element modeling of crack tip blunting for numerical estimation of fracture parameter of a Mode I crack, in weak alloy steel, which is near and normal to the interface of elastically identical but stronger maraging steel. The bimetallic body is subjected to monotonic load in K dominated regime. Crack tip yield zone across the interface, treated as Dugdale’s cohesive zone, is isolated from the bimetallic domain and is modeled alone under the action of respective cohesive stress over yield zones for obtaining the contribution of mismatch between yield strength of the steels in crack energy release rate component, . Effect of far field load on is found separately from a theoretical model. Numerical and theoretical results of are in good agreement. interface J interface J


Introduction
The effect of strength mismatch between two un-identical bodies is felt by the crack tip in parent body as it approaches the interface body [1] due to crack tip plasticity or yield zone spreading over into the interface body.If plasticity is modeled by Dugdale's hypothesis, then the part of yield zone in interface body is subjected to closing cohesive stress different from that acting over the part of yield zone in parent body tip due to strength mismatch between the bodies.As such, the plasticity induced load transfer towards the interface or parent body, depending upon the direction of strength gradient across the interface, changes the stress field around the crack which induces shielding or amplification effects at its tip.This phenomenon has also been confirmed with the help of finite element analysis [2].J integral over the path far away from the crack tip around the interface provides the applied value of J integral, applied , whereas the integral over the path near the crack tip in parent body without crossing the interface results in J integral at the crack tip, J tip J .tip J is less or more than depending upon the crack approaching a stronger or weaker material respectively.Difference in strength between the bodies results in non-homogeneity effect of the mismatched interface that is represented by, , or J integral at the interface, [3].
a J pplied inh c interface Previously, Bhat and Ukadgaonker [4], while adopting a different approach to simplify the analysis, isolated the Dugdale's cohesive zone from the bimetallic domain of elastically identical but strength and plastically mismatched steels viz., Mode I cracked ASTM 4340 alloy steel joined with ultra strong MDN 250 maraging steel and modeled it alone under the effect of cohesive stress over yield zones by finite element method.J integral, 1 J J , over the path in yield zones around the interface was obtained from the numerical solution whereas J integral, 2 J , representing the effect of far field load at the interface was determined separately from a theoretical model.Evaluation of 2 J in such a manner was justified by the fact that finite element analysis of the crack under far field load is a common exercise which previously was undertaken by many researchers with their results matching closely with the theoretical estimations. 1 J and 2 J were superimposed to obtain .The ef-fect of crack tip blunting was however not presented in that work although it was included in the analysis.Since blunting is of significance in the problem, this paper reports the methodology of its inclusion in evaluation of interface for the crack that is near and normal to the interface between the stated steels using Dugdale's plasticity model.Numerical value of interface is found to be in good agreement with that obtained from the theoretical model.

Theoretical Model
Refer Figure 1.As the crack in ductile homogenous body is subjected to monotonic far field load, its tip blunts by opening in load line direction by distance, , from the axis and by extending longitudinally by distance,  . represents the size of, highly stressed, process or fracture zone at the crack tip in which the material degradation or damage occurs by nucleation of voids [5].The parent crack finally merges with the process zone.The value of  is approximately of the order of 1.6 to 2 times the value of crack tip opening displacement (CTOD) or 3.2 to 4 times the value of   . The yield zone surrounds the process zone.The concept holds good for the crack in strength mismatched bi-material as well although the value of  may not be exactly same as that in homogenous body.
Refer   over cohesive zones in parent and interface bodies respectively in small scale yielding (SSY) or K dominated regime are given as follows [6]: applie K d in the case of Dugdale's cohesive zone across the interface is written as [7]: applied by and

Case Study
The stated theoretical model is applied to material combination of weak ASTM 4340 alloy steel, A, and strong MDN 250 maraging steel, B. Refer Table 1.The results presented as Case I and Case II, represent fracture data of alloy steel at different positions of crack near the interface of maraging steel when subjected to monotonic load in SSY regime under plane stress condition.Material and crack data are suitably selected.Refer Figure 3. Stress field at the interface of maraging steel in both the cases, defined conventionally by,     bimetallic higher order terms 2π , reveals nil yielding of maraging steel at and beyond the interface.But on viewing the bimetallic domain in comparison with the homogenous body of alloy steel, load is still transferred elastically to maraging steel due to its higher yield limit than that of alloy steel.Elastically strained zone in maraging steel under stress less than its yield strength is replaced by much smaller cohesive zone under the action of larger cohesive stress to make the application of the theoretical model possible.However, maraging steel shall also yield as the crack grows nearer towards the interface.Before undertaking the finite element analysis, the results of the selected cases are verified in the following manner:   K in the theoretical model due to fracture conditions.The value of l is iteratively assumed in Equation ( 3) till the output applied satisfies Equation (2).Final value of applied equals . Since the crack faces a stronger steel (weak-strong interface),

Finite Element Analysis and Results
The cohesive zones of both the cases are modeled by finite element method.Half of the cohesive zone is only considered in each case due to symmetry.
, at fracture, tapers from at crack tip to at the interface, finally reducing to zero at the tip of the cohesive zone.As cohesive zone has minimal lateral dimensions when compared to its length, its height is assumed constant as over distance, a, in alloy steel and over distance, l, in maraging steel to facilitate modeling. , are applied as pressure in -ve y direction over the boundary nodes of cohesive zones in steels A and B respectively.Plane stress with thickness option is adopted.Since the blunting zone is unsupported and the value of  is not exactly known, nodes ahead of the crack tip, up to the distance of , are initially left unconstrained.Remaining nodes on crack axis are constrained in y direction.The mesh model is displayed in Figure 5.

 
Using the post-processor displacement and stress solutions, the required values are mapped over chosen path, P, around the interface to obtain J integral, J 1 , from the expression, 1 d d , where are the traction components with n x and n y representing unit vectors in x and y directions and u and v as displacements in the stated directions.Since cohesive stress has a closing or compressive effect over the crack that opposes the effect of far field tensile load, the value of J 1 is taken with a negative sense.Sample stress and displacement plots near the crack tip in cohesive zone of alloy steel in Case I and Case II are presented in Figure 6 and Figure 7 respectively.The plots at the interface are available elsewhere [4].Since the material properties employed in the analysis are linear in nature, very high stresses are found to develop at certain top locations in the cohesive zone due to bending caused by the unsupported part.This however can be eliminated by using the actual elastic-plastic properties of alloy steel.However, to avoid these hypothetical values, J integral paths are made to pass only through those areas that are elastically stressed.Different cyclic paths are tried.Appropriate value of 1 J is obtained on path P that reaches right up to the crack tip where displacement values are In Case I, the condition, J and interface are equal to +115.18N/mm and +29.74 N/mm respectively.Value of interface is again in good agreement with the theoretical value of +28.3 N/mm.J J J 1 is not found to be path independent.As mentioned earlier, appropriate values of this integral are obtained over paths reaching up to the crack tip nodes where displacements are higher.On paths terminating at constrained nodes, the displacements are less that result in reduced value of the integral.As the result, is also not path independent.interface J

Conclusion
Blunted crack tip near the interface of elastically identical but strength and plastically mismatched bodies is modeled by finite element method to numerically obtain J integral, J 1 , that represents the contribution of mismatch in yield strength between the bodies on crack energy release rate component, interface .Numerical results are well validated thereby supporting the feasibility of iso-lating the cohesive zone from the bimetallic domain and modeling it alone under the dual action of un-identical cohesive stress over yield zones in parent and interface bodies.The approach is simple and reasonably accurate.J

Acknowledgements
Support received from the School of Mechanical Building Sciences, VIT University, Vellore, India during the course of this work is gratefully acknowledged.

Figure 2 .
The crack is in the parent body, A, with its tip at distance, a, from the interface body, B. Its cohesive zone has penetrated into the interface body by distance, l, with total length of cohesive zone, b, being equal to   a l  .The height of cohesive zone above the crack axis,   v x , and crack tip stress intensity parameter, tip K , under the action of applied stress intensity parame- ter, , due to far field tensile load and cohesive stress appl K ied and A B

Refer Figure 4 .
2D, 8 noded quadrilateral plane 82 elements are chosen for meshing.Cohesive stress, A  and B

Table 1 . Results of case study.
more in Case II than in Case I because the crack is nearer to the interface in Case II.The conservation of energy release rate criterion, higher.The size of  or the number of unconstrained nodes is then increased and J 1 is determined in each state till 1 , is satisfied at and beyond  of 0.304 mm.At  of 0.256 mm, 0.272 mm, 0.288 mm, J 1 is -65.0N/mm, -87.1 N/mm and -104.2N/mm respectively.Average of J 1 values is -85.44 N/mm. 2