Microsoft Word - numero_51_art_34_2679 Y. Dubyk et alii, Frattura ed Integrità Strutturale, 51 (2020) 459-466; DOI: 10.3221/IGF-ESIS.51.34 459 Stress-strain assessment of plain dents in gas pipelines Yaroslav Dubyk IPP-Centre LLC, Ukraine dubykir@gmail.com, http://orcid.org/0000-0002-4327-412X Iryna Seliverstova G. S. Pysarenko Institute for Problem of Strength, National Academy of Sciences of Ukraine, Ukraine seliverstovairyna24@gmail.com, http://orcid.org/0000-0003-1166-792X ABSTRACT. This paper presents the analytical solution of the stress-strain state for a dented pipeline, based on the method of equivalent loads. Firstly, a solution for a harmonic imperfection was found, then using Fourier series expansion a semi-analytical procedure was proposed to assess a single dent. A comparison between analytical and numerical results for the axial force and pressure load were given. The influence of the dent dimensions, shell radius to thickness ratio and initial loading to stress concentration factor were discussed. KEYWORDS. Cylindrical shell; Plane dent; Approximate solution; Citation: Dubyk Y., Seliverstova I., Stress- strain assessment of plain dents in gas pipelines, Frattura ed Integrità Strutturale, 51 (2020) 459-466. Received: 11.11.2019 Accepted: 06.12.2019 Published: 01.01.2020 Copyright: © 2020 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 runk pipelines are one of the most common means used in the transportation of various kinds of energy resources. Ensuring trouble-free operation of pipelines is a pledge of both the environmental and energy security of the country through which they pass. The most common cause of accidents on trunk pipelines is the presence of mechanical damage (dents and gauges) and corrosion [1]. Therefore, the key point in the study of the technical condition and the extension of the life of trunk pipelines is the determination of their stress state, taking into account these local damages. Dents are one of the most common types of defects arising due to external interaction – damage by heavy machinery (excavators, tractors, pipe layers), falling stones on top of a pipeline or resting on supports [2]. Local imperfections like dents on pipelines are quite complicated for engineering analysis. To assess this problem, contact mechanics and nonlinearities in both material and geometry have to be considered in all phases of the dent life cycle. Recently these problems are addressed using FEA only [3, 4]. Finite Element method is the most largely used method for solving problems of engineering but since it is a numerical study it does not give the total accurate understanding of the problem. To investigate dent problems, it would be best to develop a closed-form analytical model would be the wisest choice to investigate dent problems but it is almost impossible due to the complicated geometry and it doesn't give considerable advantages comparing to FEM analysis. Instead a semi-analytical numerical procedure is proposed in the current study based on the equivalent load method. T http://www.gruppofrattura.it/VA/51/2679.mp4 Y. Dubyk et alii, Frattura ed Integrità Strutturale, 51 (2020) 459-466; DOI: 10.3221/IGF-ESIS.51.34 460 Equivalent loads method was proposed by [5], to study the stress behavior of thin–shell structures containing geometric imperfections. The stress behavior of a geometrically imperfect shell under load is equivalent to the sum of two stress fields [6]. The first stress field is produced by the original load acting on the geometrically perfect version of the shell. The second stress field results from applying to the perfect version of the shell a load pattern that results in perturbation stress behavior equivalent to that which would be induced by the geometric imperfection. It should be noted that the equivalent load does not produce the correct dent shape but instead produce a stress field that is similar to that which would have been produced by the dent. Equivalent load is given usually in the form [6]: 2x xx x xp N N N        (1) , ,x xN N N   membrane forces, that arise from the action of the initial load on the curved shell; , ,x x     change curves in two main directions and the change curvature in twist. Eq. (1) is standard in the solution of non-ideal shell problems and is used by [7], who proposed to improve its accuracy by considering the additional terms for imperfections amplitudes larger than the shell thickness [8]. Modern regulatory documents [9] in the damage risk analysis consider only the permissible depth of such defects, neglecting other geometrical parameters and the load level. The main idea of a current study is to develop effective semi-analytical procedure for dent analysis based on equivalent loads method with second order terms to perform the convergency analysis and to investigate the dent shape effect and loading level on the stress concentrator factor. ANALYTICAL BACKGROUND General equations he theory of the thin-shells is quite well-known and effective for pipelines application. General equation of cylindrical shell equilibrium under external forces:    [ ] L u F (2)    Tu v wu is the displacement vector; ,u v and w are the orthogonal components of displacements in the ,x  and radial directions;   2 T x rR H p p p   F is the external loading vector; [ ]L is a matrix differential operator [10]: 1 [ ] [ ] [ ] [ ]D M MOD INIk H   L L L L (3) Differential matrix operator [ ]D M L corresponds to Donnel-Mushtari shell theory, the simplest one; [ ]MOD L modifying operator, for a different shell theory, further we will use Flugge shell theory; [ ]INI L differential operator which implies the initial stress influence. These operators are written in explicit form below: 2 2 2 2 2 2 2 2 2 2 4 1 1 2 2 1 1 [ ] 2 2 1 D M s ss s s k s                                                  L (4) T Y. Dubyk et alii, Frattura ed Integrità Strutturale, 51 (2020) 459-466; DOI: 10.3221/IGF-ESIS.51.34 461 2 3 3 2 3 2 2 3 2 2 3 3 3 2 3 2 2 2 1 1 0 2 2 1 1 [ ] 0 3 3 2 2 1 1 3 1 2 2 2 MOD s s s s s s s                                                          L (5) 2 2 2 2 2 2 2 2 3 2 2 2 2 2 2 2 2 2 0 1 [ ] 0 2 3 2 2 2 x x INI x x x x x N N N N x ss N N N xs s N N N N N N s x xs                                                                   L (6) The following notations are used: 2 2 2 2 ; ; ; 12 1 h Eh k H D H R R         (7) R  shell radius, h  shell thickness, E  Young’s modulus,   Poisson ratio. The complete explicit solution of Eq. (2) can be found only for simplest geometries, and loadings, for a single dent it can’t be obtained, thus a numerical procedure is developed below based on the accurate solution for the harmonic imperfection and Fourier series expansion. Harmonic imperfection A harmonic imperfection was considered as a base for further solution and the displacements representation can be found using:  cos sinumnu C n x R         ,  sin cosvmnv C n x R         ,  cos coswmnw C n x R         (8) /m R l  , ,n m  wave number in circumferential and axial directions, l  length of the shell, , ,u v wmn mn mnC C C  modes coefficients for corresponded directions. Substituting representations (8) in Eq. (2), we can get a simple algebraic set of equations:   2 2 2 2 2 2 22 2 2 2 1 1 2 2 2 0 2 2 22 u v w mn mn mnx k n kn n Nn C C C k n k knN N n Hn H                                   (9)   2 2 2 2 2 22 2 2 3 3 2 1 1 3 2 2 0 2 2 22 u v w mn mn mnx k k n Nn n C C C k kN N n H H                                   (10) Y. Dubyk et alii, Frattura ed Integrità Strutturale, 51 (2020) 459-466; DOI: 10.3221/IGF-ESIS.51.34 462 2 2 2 2 2 2 2 2 2 4 2 2 4 2 2 2 2 3 2 2 2 2 2 2 1 2 u v mn mn xw mn xx N Nn C k n k kn C k k H H N N n n Rl n R m C k k n kn kn k N N H m l                                                     (11) This set of equations can be easily solved using any mathematical software, like Maple or Mathematica, but its solution is too long to be given here. The solution of a harmonic imperfection is itself useful, as it can be used to approximate the long dent, but the solution for a single dent can be further obtained using Fourier series expansion. Solution for a single dent Dents can have different profiles, which depend on the mechanism of dent formation and working conditions. To work with a specific dent, the profile should be measured directly and afterwards approximated by a Fourier expansion but, to test the analytical procedure, the following smooth shape function was used to describe the dent profile: 2 2 0 0 1 1 ( , ) exp exp 2 2 x R x R x                               (12) 0x - length of the dent in the axial direction, 0 is the length in circumferential direction,  - dent depth. A very common way to deal with single imperfections is to use Fourier expansion. Thus, to find a solution for a dent, load coefficients (1) for every mode, using double integrals, have to be found: 2 2 2 , 0 2 0 0 00 1 1 1 1 1 exp exp 2 2 l m n x l x x p N dxd l x xx                                                     (13)   2 2 2 , 2 0 0 00 1 1 1 1 1 exp exp cos cos 2 2 l m n xx x l x x p N n x dxd l x x Rx                                                           (14)   2 2 2 , 2 0 0 00 1 1 1 1 1 exp exp cos cos 2 2 l m n l x p N n x dxd l x R                                                               (15) These integrals can be computed numerically, so it's possible to deal with any dent form, thus, in contrast to other analytical solutions, the dent form influence can be further investigated. Also, it is very useful for analysing the dents found during the inspection. For a complete solution, internal forces and moments with a simple summing for each mode must be found:  2 cos cosu v wx mn mn mn n m H x N C C n C k n R R            (16)  21 cos cosu v wmn mn mn n m H x N C C n C kn k n R R           (17)  2 22 cos cosu v wx mn mn mn n m D x M C C n C n n RR             (18) Y. Dubyk et alii, Frattura ed Integrità Strutturale, 51 (2020) 459-466; DOI: 10.3221/IGF-ESIS.51.34 463  2 22 1 cos coswmn n m D x M C n n RR         (19) Stresses were obtained using the following formulas of the shell theory, and were divided into membrane and bending part, which is very common for nuclear, oil and gas industries:   2 6 mx bx x x x N M h h      (20)   2 6 m b N M h h           (21) In the following section, this semi-analytical procedure will be tested in detail, and some properties will be discussed RESULTS AND DISCUSSION o analyze the dent behavior, the stress concentration factor will be introduced, which is the ratio of a local stress caused by the dent to the nominal stress level in a perfect pipe: mx bxforce nom SCF      (22) m b pressure nom SCF       (23) Here nom is a nominal stress level in a shell without a dent. In this study, two main loadings for oil and gas pipelines – internal pressure (P) and axial force (F) were considered: 2 x nom force N F h Rh     (24) nom pressure N PR h h    (25) First of all, the convergence criteria will be considered to the analysis, i.e. how many modes should be included in Fourier expansion. On Fig. 1, the convergence solution can be seen with respect to the number of modes involved (number of terms in Fourier expansion) and it is possible to conclude that 50th-60th modes were quite enough for dent assessment. Also, for internal pressure loading, it's possible to note that the proposed solution converges faster. On Fig. 2 and Fig. 3, a comparison with numerical results obtained by the commercial FEA ANSYS can be seen. As an applied load axial force (see Fig. 2) and pressure (Fig. 3) were used, the nominal stress level was equal to 1 MPa, note that the stress categorization according to Eqns. (20) and (21) was used. However, an equivalent load method underestimates the peak stress concentrated value, it predicts quite well the stress profile and from Fig. 1b it can be seen that 20-25 modes will be enough for profile prediction. The influence of the nominal stress level was analyzed in Fig. 4. It can be seen the continuous reduction of the SCF with the increasing of the stress level, this effect is introduced in the model with the initial stress matrix (6), without it horizontal lines would be obtained. The initial stress influence can be explained by the fact that high level forces, especially pressure, tries to ‘fix’ the dent and bring the shell to a perfect condition. In addition, from Fig. 4, the effect of non-proportional dent dimensions can be observed and conclude that a dent with a longer dimension perpendicular to a nominal stress is more dangerous. Thus, in dent assessment not only dent depth should be considered, but also length and width have to be taken into account. From Figs. 2-4, it's possible to conclude that internal pressure leads to bigger stress concentration factors with the same nominal stress level. T Y. Dubyk et alii, Frattura ed Integrità Strutturale, 51 (2020) 459-466; DOI: 10.3221/IGF-ESIS.51.34 464 Figure 1: Convergence of results with respect to Fourier terms of expansion: a - stress concentrating factor for pressure and force loading; b – profiles of circumferential bending stress for pressure loading. Figure 2: Comparison with FEM solution for axial force loading: R=400mm, h =10mm, ξ=0.5h, 0R = 0x =60mm, analytical solution; FEM solution. Figure 3: Comparison with FEM solution for internal pressure loading; R=400mm, h =10mm, ξ=0.5h, 0R = 0x =60mm, analytical solution; FEM solution. Y. Dubyk et alii, Frattura ed Integrità Strutturale, 51 (2020) 459-466; DOI: 10.3221/IGF-ESIS.51.34 465 Figure 4: Effect of initial stress: a) Axial force loading, b) Internal pressure loading; R=400mm, h =10mm: 0R = 0x =60mm; 0R =90mm, 0x =60mm; 0R =60mm, 0x =90mm. In Fig. 5, the effect of shell radius to thickness ratio was analyzed and again it is possible to see the different behavior for internal pressure and axial force loading. It can be clearly seen from Fig. 5 that with increasing R to h ratio a higher SCF was obtained and the stress concentration can be up to 40 times. In oil and gas industries typical trunk pipelines have R/h=15…40, so the results will be in the lower part of Fig. 5 and with the increasing of the nominal stress level SCF factor will decrease (see Fig. 4) and this effect will be more significant for flexible shells. Nevertheless, if using high strength steels for piping, the thickness can be reduced and must be aware of the increasing in SCF for dents in very flexible shells. Figure 5: Effect of radius to thickness ration: a) Axial force loading, b) Internal pressure loading; σnom=1MPa, h =10mm: 0R = 0x =60mm; 0R =90mm, 0x =60mm; 0R =60mm, 0x =90mm; CONCLUSIONS simple engineering approach was developed for analysis of stress-strain state of dented pipelines using equivalent load method. Using the harmonic imperfection model, an explicit analytical solution for a single dent was received and for a reliable accurate solution 50 terms in Fourie expansion must be used. Comparison of analytical stresses A Y. Dubyk et alii, Frattura ed Integrità Strutturale, 51 (2020) 459-466; DOI: 10.3221/IGF-ESIS.51.34 466 with the results of numerical simulation shows the effectiveness of the proposed procedures. However, equivalent load method underestimates the peak stress concentration value, but it predicts quite good the stress profile. The peak stress value can be adjusted by a simple semi-analytical correction. Obtained semi-analytical solution allowed to conduct parametric studies for different dent shapes, radius to thickness ratio at different nominal stress level:  for small dents with increasing depth the stresses concentration increases and for small dents, this growth was more significant  the influence of pressure, for unproportional dent was more significant than in case of axial force loading;  concentration due to the pressure load was more significant, and must be taken into account, especially for short dents, since it can exceed the nominal stress level by an excessive value;  increasing in nominal stress level gives a reduction in stress concentration factor, and for pressure loading this decreasing was more significant;  increasing in R to h ratio leads to increasing in stress concentration factor. For pressure and axial force loading this effect is different and should be accounted for very flexible pipelines. The proposed semi-analytical solution is a powerful method for dents assessment, it is especially useful for parametric analysis, and it can be used to understand the factors that influence the pipeline safe operation. REFERENCES [1] EGIG. (2015). Gas Pipeline Incidents. 9th Report of the European Gas Pipeline Incident Data Group. Available at: https://www.egig.eu/reports/$97/$155. [2] Kec, J. and Cerny, I. (2017). Stress-strain assessment of dents in wall of high pressure gas pipeline, Procedia Structural Integrity, 5, pp. 340–346. DOI: 10.1016/j.prostr.2017.07.180. [3] Li, C., and Dang, C. (2017). Plastic damage analysis of oil and gas pipelines with unconstrained and constrained dents. Engineering Failure Analysis, 77, pp. 39–49. DOI: 10.1016/j.engfailanal.2017.02.009. [4] Ying Wu, Y., Jiewen Xiao, J. and Peng Zhang, P. (2016). The analysis of damage degree of oil and gas pipeline with type II plain dent, Engineering Failure Analysis, 66, pp. 212–222. DOI: 10.1016/j.engfailanal.2016.04.004. [5] Calladine, C. R. (1972). Structural Consequences of Small Imperfections in Elastic Thin Shells of Revolution, International Journal of Solids and Structures, 8(5), pp. 679–697. DOI: 10.1016/0020-7683(72)90036-4. [6] Rinehart, A. J. and Keating, P. B. (2007). Stress Concentration Solution for a 2D Dent in an Internally Pressurized Cylinder, Journal of Engineering Mechanics, 133(7), pp. 792–800 DOI: 10.1061/(ASCE)0733-9399(2007)133:7(792). [7] Godoy, L. A. (1996). Thin-Walled Structures with Structural Imperfections, New York, Pergamon. DOI: 10.1016/B978-0-08-042266-4.X5000-3. [8] Godoy, L. A. (1993). On loads equivalent to Geometrical imperfections in shells, J. Eng. Mech, 119(1), pp. 186-190. DOI: 10.1061/(ASCE)0733-9399(1993)119:1(186). [9] API 579-1/ASME FFS-1. (2016). Fitness-For-Service. [10] Leissa, A. W. (1973). Vibration of Shells, NASA-SP-288. << /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