Microsoft Word - numero_41_art_39.docx A.S. Chernyatin et alii, Frattura ed Integrità Strutturale, 41 (2017) 293-298; DOI: 10.3221/IGF-ESIS.41.39 293 Focused on Crack Tip Fields The effect of residual stress on the nonsingular T-stresses A.S. Chernyatin, Yu.G. Matvienko, I.A. Razumovsky Mechanical Engineering Research Institute of the Russian Academy of Sciences, Moscow, Russia cas@inbox.ru, http://orcid.org/0000-0002-2708-8132 matvienko7@yahoo.com, http://orcid.org/0000-0003-2367-0966/ murza45@gmail.com ABSTRACT. The effect of residual stresses on surface fatigue crack propagation and fracture mechanics parameters of intersecting orthogonal cracks is discussed. Mutual influence of intersecting cracks and biaxiality affects the nonsingular parameters Txx and Tzz. At the same time, the effect on the stress intensity factor is negligible. Two-parameter fracture mechanics is employed for an analysis of fatigue crack propagation taking into account residual stresses. The final configuration of the crack front is close to the semi-elliptical configuration. KEYWORDS. Residual Stress; Fatigue Surface Crack; T-stresses; Intersecting Cracks; Finite Element Simulation. Citation: Chernyatin, A.S., Matvienko, Yu.G., Razumovsky, I.A., The effect of residual stress on the nonsingular T-stresses, Frattura ed Integrità Strutturale, 41 (2017) 293-298. Received: 28.02.2017 Accepted: 03.05.2017 Published: 01.07.2017 Copyright: © 2017 This is an open access article under the terms of the CC-BY 4.0, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. INTRODUCTION n analysis of propagation of surface cracks is considered in a number of publications [1-5], where the results of the experimental and numerical analysis of a flat semi-elliptical crack growth under cyclic loading are presented. The Paris formula or some subsequent modifications is used for calculations of fatigue crack increment. For assessment of the crack front configuration, the principle of self-similarity is used in most investigations. In this paper, the numerical method and the appropriate program realization are presented for the analysis of the surface crack growth in weld joint of a pipeline, including the calculation of changes in the crack front configuration due to action of the active cyclic loading, influence of no uniformly distributed residual stresses and the T-stress. It should be noted that the results of the mutual influence of intersecting perpendicular cracks on the parameters of fracture mechanics are not enough presented in the literature [6]. It is clear that the effect of the T-stress can be significant and the two-parameter fracture mechanics should be employed. The analysis of mutual influence intersecting cracks on the values and distribution of singular (the stress intensity factor) and nonsingular (the T-stress) components of the stress field along the crack fronts of intersecting cracks are also presented. A A.S. Chernyatin et alii, Frattura ed Integrità Strutturale, 41 (2017) 293-298; DOI: 10.3221/IGF-ESIS.41.39 294 FATIGUE CRACK PROPAGATION UNDER CYCLIC LOADING AND RESIDUAL STRESS Crack propagation low he method is based on the modified Forman expression which includes the nonsingular T-stress        max 1 mT c C qdl K KdN K (1) where   max minK K K is the stress intensity factor (SIF) range and  min maxR K K is the SIF ratio. The correction function qT can be represented taking into account the T-stress as follows [3]          1 T CT Y Tq e (2) The constants C, m, CT, as well as the critical stress intensity factor Kc and yield strength σY are the properties of the material. It should be noted that in this case, the SIF ratio varies along the crack front (here, it is not a characteristic of external loads) due to the effect of the residual stress (RS). Since, the fatigue crack increment occurs when load is increased to maximum values at the cycle, the T-stress in Eq. (2) should be determined as Tmax. Therefore, Eq. (1) becomes       max max, mdl C K T K dN (3) Numerical simulation of three-dimensional planar cracks can be constructed as a gradual process of crack increments (after finite number of the cycles ΔN) in a limited set of the crack front points, the position of which will be characterized by a dimensionless local coordinate s that runs along the crack front from one external points to other. An infinitely small value dN can be replaced by a finite cycle increment ΔN and dl replaced by Δl in Eq. (1). So, the crack increment at the current crack front configuration (after N cycles) can be calculated according to the follow equation               max max, , , , , m l N s C K N s T N s K N s N (4) For the successful solution of the crack propagation problem under the action of residual stresses and taking into account the constraint along the crack front by means implementation of Eq. (4), the principle of remeshing procedure of finite elements in the vicinity of the crack front is proposed. Adaptive parametric finite elements model with varying crack configuration The parametric finite element model (FE-model) of crack area in the form of a prismatic region is developed in ANSYS software. The model includes the regions with different structures of the element mesh and its sizes. Embedding the crack volume into the FE-model is carried out by means of the macro "Crack" which contains the following features [7]: - configuration of the crack front as a line can be arbitrary, but must have smoothness; - an mapped mesh of the singular elements is created along the crack front. - crack front may pass through several geometric volumes. It is especially convenient for modeling of cracks intersection. Numerical procedure of the SIF and T-stress determination The elastic crack-tip displacement and stress fields of mode I crack can be represented as follows [8] T A.S. Chernyatin et alii, Frattura ed Integrità Strutturale, 41 (2017) 293-298; DOI: 10.3221/IGF-ESIS.41.39 295 (5) where r, θ and x, y are local polar and Cartesian coordinates related to the crack tip, the plane x0y is perpendicular to the crack front line; κ is the parameter of the stress state. It should be noted that  1 2IK a ,  24Т а . The program in the MATLAB environment, so called "Williams", in conjunction with macro "Crack" provides the calculation of the coefficients of stress field expansion terms (Eq. (5) in the vicinity of the crack tip. The "Williams" interacts with the macro in ANSYS software, which enables collection of the displacements and stresses at measuring points surrounding the crack tip in plane x0y. The determination of the coefficients of expansion (5) is carried out by means of least squares method that provides comparison of these data and data calculated via expression (5). These calculations are performed in many planes which are perpendicular to the crack front for estimation of the SIF and T- stress distribution along front. Algorithm of numerical crack propagation simulation Control program (with GUI) that provides collaboration between program elements of ANSYS and MATLAB software is developed in MATLAB. A schematic diagram of its implementation is given in Fig. 1. The data files with special structures are formed at the beginning for initialization of different data, namely, sizes and geometry of the model, the crack, elements, material properties, options of loading (minimum and maximum value of cyclic loading and residual stresses distribution along thickness) and other data. The solution begins with the start of the control program and "Williams". Next step is the launching of ANSYS in the background mode for creation of the FE-model and solution of minimum and maximum values for the loading cycle. The collected data in the analyzed points are transiting to "Williams". Figure 1: The algorithm of numerical solution for the crack propagation problem. After computing the coefficients of William’s expansion for the minimum and maximum loads per cycle (as well as Kmin, Tmax, Kmax), the control program calculates values of the fatigue crack increment at each points of the crack front in accordance with Eq. (4). The new configuration of the crack front is formed and saved in the data file for the "Crack" macro. This specified sequence of procedures in ANSYS and MATLAB are executing repeatedly. A.S. Chernyatin et alii, Frattura ed Integrità Strutturale, 41 (2017) 293-298; DOI: 10.3221/IGF-ESIS.41.39 296 Internal surface crack propagation in weld joint of a pipeline Let consider a welded zone of the pipe with outer diameter 1420 mm and wall thickness 18.7 mm. The pipeline internal pressure is 10.6 MPa and it leads to the active axial and circumference stress 200 MPa and 400 MPa, respectively. Note that axial direction is perpendicular to the plane of the weld joint. In operation, there are additional axial stresses that occur because of bending and thermal loading. The axial stresses in conjunction with the pressure pulsation lead to cyclic loading. The following parameters of the loading cycle are accepted in the work: average stress is 200 MPa, the minimum and maximum values are 133 MPa and 267 MPa, respectively. In addition, there are residual stresses in the zone of welding. Distribution of axial residual stresses (RS) along the thickness was accepted in accordance with the results of well-known experimental study of a full-scale pipeline (Fig. 2). a) b) Figure 2: Loading condition of the pipe (a) and the axial residual stress distribution along pipe thickness (b) Taking into account the small curvature of the pipe, the region with the crack is considered as thick plate. The initial front geometry of the crack located in the weld joint plane on the inner pipe surface is assumed in form of a semicircle with radius of 4 mm. Ten terms in the Williams expansion are kept to estimate the SIF and the T-stress at points of the fatigue crack front after discrete cycles ΔN = 105. The total number of cycles, at which a calculation stop was occurred, (if previously SIF did not reach the critical value) is N=107. The results of calculations in the form of finite configurations of the crack front at various conditions are shown in Fig. 3. Figure 3: Final configuration of the crack front (half) under different loading conditions. It should be noted that independently of the initial crack geometry, final configuration of the crack front is close to the semi-elliptical. As expected, the stress T>0 leads to decreasing of the crack growth rate. This fact allows concluding about the conservatism of classic Paris and Forman formulas (not taking into account the effect of constraint and the residual stress). Otherwise, the residual stress (RS) affects on the rate of crack growth. Stress intensity factor distribution along the crack front strives for a constant value during crack size growing in absent of the RS. Presence of the RS leads to quite complex distribution of the SIF along the crack front. In contrast to the SIF, distribution and values of the T-stress along the crack front are weakly dependent on the RS. A.S. Chernyatin et alii, Frattura ed Integrità Strutturale, 41 (2017) 293-298; DOI: 10.3221/IGF-ESIS.41.39 297 INTERSECTING CRACKS IN PIPELINE WITH RESIDUAL STRESS racture mechanics parameters for intersecting surface cracks in the weld joint of pipeline are analyzed bellow. The pipe has the following dimensions, namely, outer diameter is 325 mm; wall thickness is 16 mm. The material is corrosion-resistant austenitic steel. The circular crack (crack A) with depth of 6 mm is assumed to be formed on the inner side of the pipe in the plane of welding zone. The semi-elliptical crack is located in the meridional section (crack B). These cracks are intersected perpendicular each other (Fig. 4). Dimensions of the elliptical crack are 8 mm (crack depth) and 12 mm (length). Figure 4: Intersecting cracks on inner surface of the pipe. The active circumferential stress is 67.6 MPa and axial stress is 108.8 MPa (taking into account not only pressure, but other loading factors). In addition, there are residual stresses occurred during the welding process. It is assumed that the axial residual stress is linearly varying through the pipe thickness. The residual stress reaches the yield strength on the inside surface near the weld root and a similar value with the opposite sign on the outer surface. The circumferential residual stress is also equal to the yield strength and this stress is constant through the thickness (they are varying along the pipe). a) b) Figure 5: Distribution of KI, Txx, Tzz along the crack front A (a) and crack front B (b) under different loading conditions. The computed results are summarized in Fig. 5. Note that sA and sB are dimensionless local coordinates (Fig. 4) transiting along a front of the crack A and B, respectively. Distributions of the Tzz-stress [9] as well as the SIF and the Txx-stress are determined by means of the program "Williams" for only active external stress ("act") as well as in combination with the residual stress ("res"). An analysis of fracture mechanics parameters allows concluding the following. The axial crack (B) has a great influence on the Txx-stress and the Tzz-stress. The values of Txx and especially Tzz are decreasing (staying more F A.S. Chernyatin et alii, Frattura ed Integrità Strutturale, 41 (2017) 293-298; DOI: 10.3221/IGF-ESIS.41.39 298 negative). Strong influence of the RS on the SIF of the crack A is not observing. However, the residual stress acts on SIF values in more manner then on T-stress values of crack A. The distribution of fracture parameters and the effect of the RS for the crack B are more complex. The SIF has an almost constant value along the crack front, whereas T-stresses have values great varied. It should be noted that there is the influence of the crack A on the state of the crack. CONCLUSIONS he effect of residual stresses on surface fatigue crack propagation and fracture mechanics parameters of intersecting orthogonal cracks in the pipeline is discussed. The method on the basis of modified Foreman equation and program software for realization of gradual remeshing of the finite element model during incremental crack growth is developed. It is implemented for an numerical analysis of fatigue surface crack propagation in welded area of pipeline taking into account nonlinear distributed residual stress and constraint effect in the crack tip employing two-parameter fracture mechanics. In spite of the fact that the residual stress and the nonsingular T-stress have a significant influence on the crack growth rate. Final configuration of the crack front is close to the semi-elliptical configuration. Mutual influence of intersecting surface cracks and biaxiality effects on the nonsingular parameters Txx and Tzz is observed. At the same time, their effect on the stress intensity factor is negligible. ACKNOWLEDGMENTS he authors acknowledge the support of the Russian Science Foundation (Project N 14-19-00383). REFERENCES [1] Zhao, C.Y., Huang, P.Y., Zhou, H., Zheng, X.H., Numerical analysis of KI of semi-elliptical surface crack in steel structure strengthened with FRP under tensile load, J. App. Mech. Mat., 137 (2012) 42–49. [2] Roychowdhury, S., Dodds, R.H., Three-dimensional effects on fatigue crack closure in the small scale yielding regime - a finite element study, J. Fat. Fract. Eng. Mat. Struct., 8 (2003) 663–673. [3] Hamam, R., Pommier, S., Bumbieler, F. Mode I fatigue crack growth under biaxial loading, Int. J. Fatigue, 27 (2005) 1342–1346. [4] Hoshide, T, Tanaka, K., Stress-ratio effect of fatigue crack propagation in a biaxial stress field, J. Fat. Fract. Eng. Mat. Struct., 4 (1981) 355–366. [5] Forman, R. G., Study of fatigue crack initiation from flaws using fracture mechanics theory, J. Eng. Fract. Mech., 2 (1972) 333–345. [6] Chen, Y.Z., Evaluation of T-stresses in multiple crack problems of finite plate, J. Fat. Fract. Eng. Mat. Struct., 35 (2012) 173–184. [7] Razumovskii, I.A., Chernyatin, A.S., Experimental numerical method of loading estimation for structures with surface cracks, J. Machinery Manufacture and Reliability, 38-3 (2009) 247–253. [8] Williams, M.L., On the stress distribution at the base of a stationary crack, J. Appl. Mech., 24 (1957) 109–114. [9] Nakamura, T., Parks, D.M., Determination of elastic T-stress along three-dimensional crack front an interaction integral, Int. J. Solid and Struct., 29 (1992) 1597–1611. T T << /ASCII85EncodePages false /AllowTransparency false /AutoPositionEPSFiles true /AutoRotatePages /None /Binding /Left /CalGrayProfile (Dot Gain 20%) /CalRGBProfile (sRGB IEC61966-2.1) /CalCMYKProfile (U.S. Web Coated \050SWOP\051 v2) /sRGBProfile (sRGB IEC61966-2.1) /CannotEmbedFontPolicy /Error /CompatibilityLevel 1.4 /CompressObjects /Tags /CompressPages true /ConvertImagesToIndexed true /PassThroughJPEGImages true /CreateJobTicket false /DefaultRenderingIntent /Default /DetectBlends true /DetectCurves 0.0000 /ColorConversionStrategy /CMYK /DoThumbnails false /EmbedAllFonts true /EmbedOpenType false /ParseICCProfilesInComments true /EmbedJobOptions true /DSCReportingLevel 0 /EmitDSCWarnings false /EndPage -1 /ImageMemory 1048576 /LockDistillerParams false /MaxSubsetPct 100 /Optimize true /OPM 1 /ParseDSCComments true /ParseDSCCommentsForDocInfo true /PreserveCopyPage true /PreserveDICMYKValues true /PreserveEPSInfo true /PreserveFlatness true /PreserveHalftoneInfo false /PreserveOPIComments true /PreserveOverprintSettings true /StartPage 1 /SubsetFonts true /TransferFunctionInfo /Apply /UCRandBGInfo /Preserve /UsePrologue false /ColorSettingsFile () /AlwaysEmbed [ true ] /NeverEmbed [ true ] /AntiAliasColorImages false /CropColorImages true /ColorImageMinResolution 300 /ColorImageMinResolutionPolicy /OK /DownsampleColorImages true /ColorImageDownsampleType /Bicubic /ColorImageResolution 300 /ColorImageDepth -1 /ColorImageMinDownsampleDepth 1 /ColorImageDownsampleThreshold 1.50000 /EncodeColorImages true /ColorImageFilter /DCTEncode /AutoFilterColorImages true /ColorImageAutoFilterStrategy /JPEG /ColorACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /ColorImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000ColorACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000ColorImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasGrayImages false /CropGrayImages true /GrayImageMinResolution 300 /GrayImageMinResolutionPolicy /OK /DownsampleGrayImages true /GrayImageDownsampleType /Bicubic /GrayImageResolution 300 /GrayImageDepth -1 /GrayImageMinDownsampleDepth 2 /GrayImageDownsampleThreshold 1.50000 /EncodeGrayImages true /GrayImageFilter /DCTEncode /AutoFilterGrayImages true /GrayImageAutoFilterStrategy /JPEG /GrayACSImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /GrayImageDict << /QFactor 0.15 /HSamples [1 1 1 1] /VSamples [1 1 1 1] >> /JPEG2000GrayACSImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /JPEG2000GrayImageDict << /TileWidth 256 /TileHeight 256 /Quality 30 >> /AntiAliasMonoImages false /CropMonoImages true /MonoImageMinResolution 1200 /MonoImageMinResolutionPolicy /OK /DownsampleMonoImages true /MonoImageDownsampleType /Bicubic /MonoImageResolution 1200 /MonoImageDepth -1 /MonoImageDownsampleThreshold 1.50000 /EncodeMonoImages true /MonoImageFilter /CCITTFaxEncode /MonoImageDict << /K -1 >> /AllowPSXObjects false /CheckCompliance [ /None ] /PDFX1aCheck false /PDFX3Check false /PDFXCompliantPDFOnly false /PDFXNoTrimBoxError true /PDFXTrimBoxToMediaBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXSetBleedBoxToMediaBox true /PDFXBleedBoxToTrimBoxOffset [ 0.00000 0.00000 0.00000 0.00000 ] /PDFXOutputIntentProfile () /PDFXOutputConditionIdentifier () /PDFXOutputCondition () /PDFXRegistryName () /PDFXTrapped /False /CreateJDFFile false /Description << /ARA /BGR /CHS /CHT /CZE /DAN /DEU /ESP /ETI /FRA /GRE /HEB /HRV (Za stvaranje Adobe PDF dokumenata najpogodnijih za visokokvalitetni ispis prije tiskanja koristite ove postavke. Stvoreni PDF dokumenti mogu se otvoriti Acrobat i Adobe Reader 5.0 i kasnijim verzijama.) /HUN /ITA /JPN /KOR /LTH /LVI /NLD (Gebruik deze instellingen om Adobe PDF-documenten te maken die zijn geoptimaliseerd voor prepress-afdrukken van hoge kwaliteit. De gemaakte PDF-documenten kunnen worden geopend met Acrobat en Adobe Reader 5.0 en hoger.) /NOR /POL /PTB /RUM /RUS /SKY /SLV /SUO /SVE /TUR /UKR /ENU (Use these settings to create Adobe PDF documents best suited for high-quality prepress printing. Created PDF documents can be opened with Acrobat and Adobe Reader 5.0 and later.) >> /Namespace [ (Adobe) (Common) (1.0) ] /OtherNamespaces [ << /AsReaderSpreads false /CropImagesToFrames true /ErrorControl /WarnAndContinue /FlattenerIgnoreSpreadOverrides false /IncludeGuidesGrids false /IncludeNonPrinting false /IncludeSlug false /Namespace [ (Adobe) (InDesign) (4.0) ] /OmitPlacedBitmaps false /OmitPlacedEPS false /OmitPlacedPDF false /SimulateOverprint /Legacy >> << /AddBleedMarks false /AddColorBars false /AddCropMarks false /AddPageInfo false /AddRegMarks false /ConvertColors /ConvertToCMYK /DestinationProfileName () /DestinationProfileSelector /DocumentCMYK /Downsample16BitImages true /FlattenerPreset << /PresetSelector /MediumResolution >> /FormElements false /GenerateStructure false /IncludeBookmarks false /IncludeHyperlinks false /IncludeInteractive false /IncludeLayers false /IncludeProfiles false /MultimediaHandling /UseObjectSettings /Namespace [ (Adobe) (CreativeSuite) (2.0) ] /PDFXOutputIntentProfileSelector /DocumentCMYK /PreserveEditing true /UntaggedCMYKHandling /LeaveUntagged /UntaggedRGBHandling /UseDocumentProfile /UseDocumentBleed false >> ] >> setdistillerparams << /HWResolution [2400 2400] /PageSize [612.000 792.000] >> setpagedevice