Microsoft Word - numero 25 art 23 C. J. Christopher et alii, Frattura ed Integrità Strutturale, 25 (2013) 161-166; DOI: 10.3221/IGF-ESIS.25.23 161 Special Issue: Characterization of Crack Tip Stress Field Extension of the CJP model to mixed mode I and mode II C. J. Christopher, G. Laboviciute, M. N. James University of Plymouth, Drake Circus, Plymouth PL4 8AA, UK mjames@plymouth.ac.uk E. A. Patterson University of Liverpool, Brownlow Hill, Liverpool, L69 3GH, UK ABSTRACT. The present authors have previously proposed a novel ‘plastic inclusion’ approach for dealing with the local plasticity which occurs at the tip of a growing fatigue crack. This meso-scale model provides a modified set of crack tip stress intensity factors that include the magnitude of plastic wake-induced crack tip shielding and which have the potential to help resolve some long-standing controversies associated with plasticity-induced closure. The present work extends the CJP model to deal with the case of mixed Mode I and Mode II loading and thus opens up enhanced possibilities for testing it on inclined cracks in metallic specimens. This extension requires the addition of only one new force parameter to the model, i.e. an anti-symmetric shear force on either side of the crack. KEYWORDS. Mixed mode fatigue; CJP crack tip stress model; plastic inclusion; crack tip shielding. INTRODUCTION he present authors have previously proposed a novel ‘plastic inclusion’ approach for dealing with the local plasticity which occurs at the tip of a growing fatigue crack [1]. Localised plasticity arises from crack growth mechanisms and essentially blunts the crack, creates a reversed cyclic plastic zone, and induces shear along the crack flanks, along with the possible generation of wake contact stresses which act on the applied elastic stress field at the boundary of the elastic-plastic enclave surrounding the crack. The outcome of this meso-scale model is a modified set of crack tip stress intensity factors that include the magnitude of plastic wake-induced crack tip shielding and which have the potential to help resolve some long-standing controversies associated with plasticity-induced closure. A full-field approach has been developed for stress using photoelasticity and also for displacement using digital image correlation. This model has been termed the CJP model by the authors, and is independent of the mechanisms of plastic deformation and is therefore potentially applicable to a variety of materials. The definition of the forces on the crack allows roughness-induced closure to also be accounted for in the calculated stress intensity factors. The model can also be used to mathematically explore the effect on shape of the crack tip stress field of changes in magnitude of the various parameters; for, example, the effect of variation in magnitude of increasing positive or negative T-stress on crack tip fringe patterns. Under Mode I loading, the new model uses four parameters to characterize the stress fields generated by the forces in Fig. 1; an opening mode stress intensity factor KF, the shear stress intensity factor KS, the retardation stress intensity factor KR, and the T-stress. In applications involving the DIC technique, stress intensity factors in the new four-parameter model can be solved directly from measured displacement fields using Muskhelishvili’s potential functions. T http://dx.medra.org/10.3221/IGF-ESIS.25.23&auth=true http://www.gruppofrattura.it C. J. Christopher, Frattura ed Integrità Strutturale, 25 (2013) 161-166; DOI: 10.3221/IGF-ESIS.25.23 162 Figure 1: Schematic idealisations of forces acting at the interface of the plastic enclave and the surrounding elastic material, where FA is the Mode I applied force generating the crack tip stress field usually characterized by KI, FB is the additional applied force induced by Mode II loading, FT represents the force due to the T-stress shown in this example as positive, FS is the interfacial shear force between the elastic and plastic zones, FC and FP together create the shielding effect. FP is the force generated by the constraint of compatibility on the plastically deformed material and FC is the contact force between the flanks of the crack generated by the interference of the plastic zones along the flanks [1]. The present work extends the CJP model to deal with the case of mixed Mode I and Mode II loading and thus opens up enhanced possibilities for testing it on inclined cracks in metallic specimens. This extension requires the addition of only one new force parameter to the model, shown as FB in Fig. 1, i.e. an anti-symmetric shear force on either side of the crack. THE CJP MODEL UNDER BIAXIAL (MODE I AND MODE II) LOADING he CJP Model [1] is defined as: 1 3 1 3 02 2 2 2 y x xy2i Az z z CB ln(z)z Dz z z ln(z)E            (1) The extension to the model was obtained by making the coefficients A and B complex and making the assumptions BiArA i3 , B BBr i i  , D E 0  . 1 3 1 3 02 2 2 2 y x xy r i r i2i ( A i3B )z iB )z z Cz Dz z z ln(z( )B ln(z) E               (2) The elastic stress fields near the crack tip are given by: Forces on plastic field Forces on elastic field Maximum applied load Minimum applied load FS FS FS FS FS FS FS FS FAy FAy FAy FAy FAx FAx FPy FPy FPx FPx FPy FPy FC FC FC FC FT FT FT FT FT FT Notational boundary of  plastically deformed material Closed crack with flank contact FT Open crack  FT Forces on plastic field Forces on elastic field Maximum applied load Minimum applied load FS FS FS FS FS FS FS FS FAy FAy FAy FAy FAx FAx FPy FPy FPx FPx FPy FPy FC FC FC FC FTFT FTFT FTFT FTFT FTFT FTFT Notational boundary of  plastically deformed material Closed crack with flank contact FTFT Open crack  FTFT T http://dx.medra.org/10.3221/IGF-ESIS.25.23&auth=true http://www.gruppofrattura.it C. J. Christopher et alii, Frattura ed Integrità Strutturale, 25 (2013) 161-166; DOI: 10.3221/IGF-ESIS.25.23 163 1 1 1 2 2 2 y r r r i 1 1 2 2 1 1 1 2 2 2 r r ix r 1 1 5 1 5 ( A cos B B sin 2 2 2 2 2 2 2 1 5 5 E 4B 8E )r r cos r sin r ln(r ) cos sin ) 4B 8E)r r c 5cos 5sin O(r 2 2 2 2 2 1 1 5 1 5 ( A cos B Bos r si 7sin 2 2 n 2 2 2 2 2                                                          1 1 2 2 1 1 2 2 xy r r i 1 1 2 2 1 5 5 E 3cos 3sin C O(r 2 2 2 2 2 1 5 1 5 r B B r cos 3cos 2 2 2 2 2 2 3 3 E ln(r )cos sin O r ln(r ) cos sin ) A sin s (r 2 2 in r sin )                                                            (3) FK is defined from the asymptotic limit of y as x 0  , along y 0 , i.e. towards the crack tip from the front along the crack line: r 0 1 2 F y r rK lim 2Er ln(r ))2 r ( ( A 3B 8E ) 2               (4) RK was obtained by evaluating x in the limit as x 0  , along y 0 , i.e. towards the crack tip from behind along the crack flank: R i r x 0 2 rK lim (2B E 2 4 )            (5) S K is derived from the asymptotic limit of xy as x 0  , along y 0 , i.e. towards the crack tip from behind along the crack flank: S 0 xy r r r 2 rK l A )im ( B 2           (6) The +ve sign indicates y 0 , and a –ve sign that y 0 . The quantity II K characterizes mode II loading, and is derived from the asymptotic limit of xy as x 0  , along y 0 , i.e. towards the crack tip from the front along the crack line: xyI iI r 0 2 r 2 BK lim 2        (7) T-stress is the transverse stress which is added to x as a constant term and is given by T C  (8) The new five-parameter model can be solved for displacement fields: http://dx.medra.org/10.3221/IGF-ESIS.25.23&auth=true http://www.gruppofrattura.it C. J. Christopher, Frattura ed Integrità Strutturale, 25 (2013) 161-166; DOI: 10.3221/IGF-ESIS.25.23 164 1 1 2 2 x y r i 1 1 2 2 r i 1 1 1 2 2 2 r i C (u iu ) 2(B iB )z z 4 C z (B iB 2E)z z ln(z) 4 C ( A i3B )z z ln(z 2 2Ez ln( ) z) E D 2D z z 2                                  (9) Where xu and yu are the horizontal and vertical displacements respectively, E 2(1 )     ; 3 4   (plane strain) or 3 1       (plane stress). xu and yu are shown explicitly below with the assumption D E 0  . 1 1 1 2 2 2 x r r i i r 1 1 2 2 i 2 ( A (2B (B 2E)cos B sin E ln(r ) 3 u r 2B 2E )cos r 3B )sin r 2 2 2 3 3 3 r r (1 2 )cos sin (1 2 )sin 2 2 2 2 2 C r(1 )c 4 co os s                                            (10) 1 1 1 2 2 2 i i r r r 1 1 2 2 i y2 ( 2B ( A (B 2E )sin B cos E ln(r ) 3 u r 3B )cos r 2B 2E )sin r 2 2 2 3 3 3 r r sin (1 2 )sin (1 2 )cos 2 2 2 2 2 C r(3 )sin 4 cos                                             (11) APPLICATION OF THE MODEL TO MIXED MODE CRACK GROWTH hus far, the model has been tested against limited data obtained from an inclined crack in a 2mm thick 2024-T6 aluminium CT specimen with non-standard dimensions, which is shown in Fig. 2. A jeweller’s saw with blade thickness of 0.15 mm was used to extend the notch tip into a 45° slit some 5 mm long, which is schematically shown in Fig. 2. A fatigue crack some 2 mm long was then grown in Mode I collinear with this slit. This was achieved by starting with a larger dimension CT specimen with additional loading holes in a similar fashion to the disk-shaped compact specimen described by Ding et al [2]. The specimen was then machined to final dimensions and approximately 100 cycles of loading was applied in the usual orientation, so as to put the fatigue crack under a combination of mixed Mode I and Mode II loading. The applied load ratio was R = 0.1 and the peak load was 1.2 kN. A Dantec digital image correlation (DIC) system operating in 2D mode was used to measure the crack tip displacement field and to compare the predictions of the CJP model with the measured strain field data. A facet size of 17 pixels with a centre-line pitch of 17 pixels was used with a scale of 107 pixels per millimetre. It is acknowledged that this method of cracking is not optimum, as it produces a Mode I plastic wake, rather than one that reflects mixed mode loading. However, the technique gives data for a mixed mode crack tip strain field and a Mode I crack wake and this is not dissimilar to the data acquired by many other workers, e.g. [2]. Attention will be given in the immediate future to testing cracks of various inclinations, i.e. various mode-mixity ratios, which have been developed under true mixed mode loading, as well as under Mode I loading. Fig. 3 presents a comparison between the measured x and y components of the crack tip displacement and those predicted by the 5-parameter CJP mixed mode model, whilst Fig. 4 show typical stress intensity data measured through a loading half cycle. In Fig. 4 the stress intensity factors are all defined with respect to the local crack plane, i.e. KF and KII values are defined as respectively being perpendicular to, and parallel with, the plane of the inclined crack. The stress intensity factors can be fairly easily re-defined using simple geometry, with Cartesian coordinates defined in terms of the load line T http://dx.medra.org/10.3221/IGF-ESIS.25.23&auth=true http://www.gruppofrattura.it C. J. Christopher et alii, Frattura ed Integrità Strutturale, 25 (2013) 161-166; DOI: 10.3221/IGF-ESIS.25.23 165 and the plane of the original notch in the CT specimen, e.g. KF = (KF+KII)/√2 and KII = (Kf-KII)/√2 (or the negative of this, depending on how the force directions were chosen). However, the belief of the present authors is that it is more useful to define stress intensity factors in terms of the local crack plane, because crack growth rate is governed by these local conditions. : Geometry and dimensions (in millimetres) of the polycarbonate compact tension (CT) specimens used in this work. Note Figure 2 that the notch length (equal to the initial crack length) ai = 20 mm and the width b = W = 72 mm are defined from the centreline of the holes. DIC: y-displacement a) CJP model: y- b) DIC: x-displacement c) CJP model: x- d) : Comparison between DIC and CJP model predictions of the x and y-displacement fields at the tip of a crack under mixed Figure 3 mode fatigue loading. http://dx.medra.org/10.3221/IGF-ESIS.25.23&auth=true http://www.gruppofrattura.it C. J. Christopher, Frattura ed Integrità Strutturale, 25 (2013) 161-166; DOI: 10.3221/IGF-ESIS.25.23 166 : Stress intensity values given by the mixed-mode CJP equation, with KF defined perpendicular to the 45° slant crack and KII Figure 4 defined parallel with the slant crack in the CT specimen. CONCLUSIONS he CJP model has been extended to deal with fatigue cracking under mixed mode I and II conditions, by incorporating an additional anti-symmetric horizontal shear force into the analysis. The model appears able to accurately predict the x and y-displacement fields ahead of the tip of a crack under mixed mode loading. Further work is on-going to explore the accuracy of the predictions, the values of the mixed mode stress intensity factors and their use in predicting fatigue crack growth rates. REFERENCES [1] James, M.N., Christopher, C.J., Lu, Y. and Patterson, E.A., Local crack plasticity and its influences on the global elastic stress field. Int. J. Fatigue, 46 (2013) 4-15. [2] Ding, F, Zhao, T. and Jiang, Y., A study of fatigue crack growth with changing load direction. Engng Fract. Mech., 74 (2007) 2014-2029. T http://dx.medra.org/10.3221/IGF-ESIS.25.23&auth=true http://www.gruppofrattura.it