Mathematical Models of Pipelines Alternative Stress States Viktor A. Rukavishnikova and Oleg P. Tkachenkoa a Computing Center of the Far Eastern Branch of the Russian Academy of Sciences, Kim Yu Chen Str., 65, Khabarovsk, 680000, Russia Abstract A pipeline mathematical model as moment shell with singularities caused by the domain geometry is built. The parameters that differentiate the various stress states of the pipeline are determined. Two cases of pipeline geometry are considered, which differ in the numerical parameter we have selected. The analysis of the limiting states of the pipeline is performed. The existence of singularities in solutions is established. A numerical analysis of the stress state is carried out for a weakly bent pipeline and a pipeline with a singularity. Keywords 1 Pipeline, stress-strain state, singularities, computing experiment 1. Introduction Currently, many scientific problems associated with pipeline transportation of oil products remain unresolved. Here we analyse two such problems: 1. The problem of pipeline deviation from the design position [1], associated with the problem of studying the pipeline dynamics in a deformable medium. The curvature of the pipe centerline is small. 2. The problem of stress concentration at pipe joints. The curvature of the pipe centerline tends to infinity. The problem of imposing conjugation conditions on the shells junction was studied in [2]. The aim of the work is numerical study of the features of the stress-strain state of the pipeline under these limiting geometries. Research objectives: (1) defining a numerical criterion for detecting limiting cases; (2) building mathematical models; (3) setting up computational experiments. 2. Problem formulation We investigate a pipeline of length L with a circular cross-section of radius R0 and a wall thickness h (see [3, p. 908]). The pipe centerline is curved along the curve  0 = { x0 , y 0 : x0 = x0 ( s ), y 0 = y 0 ( s )} , where s is the arc length. The pipe is filled with a stationary fluid flow with a velocity  s 0 . We define the curvature parameter  = R0 max  0 ( s ) , (1) where  0 ( s ) is the initial curvature of the line  . Let the basic geometric relation of the theory of elastic cylindrical shells be satisfied: h R0 = 1 . The parameter (1) allows distinguishing two limiting cases of pipe geometry:  = 1 – weakly bent pipe; (2) VI International Conference Information Technologies and High-Performance Computing (ITHPC-2021), September 14–16, 2021, Khabarovsk, Russia EMAIL: vark0102@mail.ru (V.A. Rukavishnikov); olegt1964@gmail.com (O.P. Tkachenko) ORCID: 0000-0002-3702-1126 (V.A. Rukavishnikov); 0000-0003-1806-0274 (O.P. Tkachenko) ©️ 2021 Copyright for this paper by its authors. Use permitted under Creative Commons License Attribution 4.0 International (CC BY 4.0). CEUR Workshop Proceedings (CEUR-WS.org) 173  →  – pipe with singularity. (3) In the case of (2), the pipe is considered as a technical Vlasov shell (see [4]). In the case of (3), the pipe is considered as a moment shell (see [2]). The pipe geometry at implementation (2) was examined in [3, 5]. Let's introduce curvilinear coordinates: s – defined above;  and R are the angle and radius of polar coordinates in the section s . We find for the middle surface of the pipe wall following [6]: A = 1 + R0 ( s , t ) sin  , B = R0 , k1 =  ( s , t ) sin  , k2 = 1 ; (4) (1 +  ( s , t ) R0 sin  ) R0 where k1 and k 2 are the main curvatures of a median surface,  is the axis curvature. The geometry of a pipeline with a kink in the profile was studied in [7, p. 7553]. Let's consider case (3) on the example of two cylindrical pipes connected at right angles. We denote: R i – inner radius of the pipe; R e – outer radius of the pipe; L1 , L 2 – lengths of the first and second pipe sections along the centerline, respectively; R0 = 0.5( Re + Ri ) . The connected pipes are hereinafter referred to as sections (1) and (2). We construct cylindrical coordinate systems on sections (1), (2), and denote these coordinates by si (axial),  i (angular) and  i (radial), where i is the section number. We fix the numbering of curvilinear coordinates: x1 = s , x2 =  , x3 =  . (5) 2.1. Equations of the pipeline statics and dynamics The dynamics of a pipeline is governed by the equations of an elastic body [6]:  t a k =  i ki , (6) where  t is the density of the pipe material, a are the acceleration components,  are the stress- k ki tensor components, and  i is the covariant derivative. In the stationary case, in the absence of external distributed loads, the equilibrium equations are as follows:  i ki =0. (7) We use the equations (6) to describe the case (2), that is, the dynamics of the slow motion of a curved pipeline. We use the equations (7) to describe a pipeline with singular profile, that is, for the case (3). 2.1.1. Dynamics of a bent pipeline In the case (2), the movement of the inner flow is considered as quasi stationary. The Darcy's law of a friction [8] was chosen as the law of hydraulic resistance. Ibid also given the equations of stationary motion of an incompressible fluid taking into account the Darcy friction force Φ ( s 0 ) . We denoted:  f – fluid density,  f – fluid viscosity, p – fluid pressure. For a bent pipe, the equations system is obtained [3]: 1 I 1 −   0  k w  1 − (0) 2 − + (1 −  )  k1 k 2 u − 2 =− X, A s B   A s  Eh 1 I 1 −   0  k w  1 − (0) 2 + + (1 −  )  k1 k 2 v − 1 =− Y, B  A s  B   Eh 1 −     h 2 − ( k1 + k 2 )  I + 2 ABk1 k 2 w + ( Bk 2 u ) + ( Ak1v )  −   w− (0) 2 2 AB  s   12 h2 1 − 2 −  2  ( k12 + k 22 ) w  = − Z; (8) 12   Eh 174 1    B    A     =    +   2  . AB   s  A  s     B     1  u  1 A  1   w   v   2 2 A v I (0) =  + + v cos   +  +  sin   w − 2   +   , A   s R0    A  R0  2 A   s   s   1    0 = ( Bv ) − ( Au )  ; 2 AB   s   1  2u 1 X = − t +  t ( s 0 ) , h t 2 h 1  2v 2  u * cos  Y = − t − , h t 2    eu *  hR0  0.5 − ln R0   4     1 2w 1 Z = − t + ( p − p e ). h t 2 h We denoted: u , v , w – displacements of the median pipe surface along the coordinates s ,  , R ; X , Y , Z – the components of the forces density acting on the shell along the coordinates s ,  , R ; pe – external pressure. A , B are defined by formulas (4). The system of equations (8) is supplemented by the boundary conditions of rigid fixing and by the homogeneous initial conditions. 2.1.2. Static of a pipeline with a singular profile In the case of (3), we use the equations (7). The problem statement within the framework of the moment shells theory is investigated in [7]. The mathematical model equations are derived in the coordinates (5):  u i 1 −   u i 1 +   vi  wi h  wi (1 −  ) h  wi 1 − 2 2 2 2 3 2 3 2 + + + − + + X i = 0,  si 2 R0   i 2 R0 si   i R0 si 12 R0 si 24 R0 si   i 2 2 2 3 3 2 Eh 1 +   ui 1  vi 1 −   vi 1 wi 3 −  wi 1 − 2 2 2 2 3 2 h + + + − + Yi = 0, (9) 2 R0  s i   i R0   i  si R0   i 2 12 R0   i si 2 2 2 2 2 2 2 Eh ui 1  vi h   ui 2  wi 3 −   v i 2 3 2 3 wi − − − +  3 − 3 + −  si R0   i R0 12   si R0   i2 2 R0   i  si2 1 −  3ui  h 2   4 wi 2  wi 4 1  wi  4 1 − 2 −  − R  + +  + R Z i = 0. 2 R02  si   i2  12   si R02  si2   i2 R04   i4  0 4 0 Eh Here X i , Yi , Z i are the external force components acting on the shell; u i , vi , wi are components of the middle surface displacement vector of the i -th section. On the pipes outer ends, we impose the boundary conditions of rigid fastening. We imposed conjugation conditions on the shells connection line. We obtained formulas for the conjugation conditions in [7]. Geometric conjugation conditions: v2 − u 2 cos  = v1 + u1 cos  , u 2 + v2 cos  + w2 sin  = u1 − v1 cos  − w1 sin  , − u 2 sin  − v 2 sin  cos  + w2 (1 + cos 2  ) = (10) = u1 sin  − v1 sin  cos  + w1 (1 + cos 2  ) ,  w1  w2 = .   Conjugation conditions for force factors: 175 M (1) = M (2) ; S (1) = − S (2) ; −Q 1 + 3 cos  + N sin  = N sin  + Q 1 + 3 cos  ; (2) 2 (2) (1) (1) 2 (11) −Q (2) sin  − N (2) 1 + 3 cos  = N 2 (1) 1 + 3 cos  − Q2 (1) sin  . We denoted: M – bending moments, S – shear forces, Q – cutting efforts, N ( i ) – normal (i ) (i ) (i ) efforts, ( i ) stands for section number. The mathematical model consists of: (1) system of equations (9) with fixed constrains as boundary conditions; (2) geometric conditions on the conjugation line (10); (3) interface conditions for force factors (11). 3. Methods of numerical experiments 3.1. Pipe bending problem Let's introduce dimensionless variables into the equations (8):  = s , r = R R0 ,  =  ,  =  t . Here ,  are the characteristic length and frequency of processes in the pipeline. Displacements of the pipe middle surface: u  = u R0 , v  = v R0 , w  = w R0 ; fluid pressure: p  = p p a . We obtain the dimensionless form of the equation (8). Then we represent their solutions in the form: u ( ,  , ) = u 0 ( ) +  u1 ( , ) sin  + O (  2 ) ; v ( ,  , ) =  v1 ( , ) cos  + O (  ) ; 2 (12) w ( ,  , ) = w0 ( ) +  w1 ( , ) sin  + O (  ) ; 2 p  (  ,  , r ) = p (0) (  , r ) +  p (1) (  , r ) sin  +  p (2) (  , r ) cos  + O (  2 ) . Within the scope of this paper, the dynamics of the fluid is considered known and its study is available in [5, 3]. The zero-order solutions (12) are also supposed to be known. The equations of the first approximation (see [3]):  2 u1 1 − 1 +  v1  w1 2 − u1 −  +  +    2 2 2 1 −   2u0 w  +f u 0 − 2 2 +  (1 −  ) 0  −  2  2     w  2 w0  w0  2 w1   w0  2 w0  t R02 2  2 u1 − 3  1 +  + 3 3 f  = ; (13)      2    2  2 2 E* 1 −  2 v1 1 2u1*  1 +  u1 2 − v1 − +  + 2  2 E h * *    u  * 2  R0  0.5 − ln | e R0 |  1  4    3 −  u 0  2  w0  w1  t R0   2 v1 2 2 + w1 + f  w0 −   −  = ; (14)  2     E*  2 h *2  4  w1 2  w1  u 4 2 w1 +  −  2  +  1 − v1 + 12   4    2  u 0  2  w0  w1  2   w0  + f  2 w0 + (1 −  )  −  + f  =      2      1 2 u1  *   R 2 2  2 w = * *   f s 0 f − − 2 t 0 1 . (15) E h    e  u1*  E *  2 R0  0.5 − ln | R0 |    4   176 Thus, the three-dimensional problem (8) is reduced to a one-dimensional formulation. A difference scheme for numerical solution of equations (13)–(15) was constructed in [3]. 3.2. Numerical method for a pipeline with a singular profile In the problem with a kinked pipe, there are points of singularity of the stress field, as indicated in [2]. The problem of calculating such a stress field in a mathematical sense is close to the problem of calculating stresses in the L-shaped domain, see [9]. Therefore, to solve the problem under the condition (3), it is necessary to develop a new computational algorithm. We plan to create this algorithm based on the approach developed in [9, 10, 11, 12]. Numerical experiments on calculating the stress-strain state of a pipeline with a break in the profile were performed in the FreeCAD software package (see [13]) to illustrate the existence of a singularity and estimate the limiting stress values. We created a pipeline modeling algorithm in the FreeCAD software package. It provides for the creation of a solid 3D model, mesh generation by the finite element method, data entry, solver setup, and visualization of calculation results. The CalculiX finite element method package and the NetGen meshing package are used. 4. Numerical results 4.1. The pipeline bending problem Physical and geometric parameters of the test problem: h = 0.005 m ,  e = 1700 kg m 3 ,  = 10000 N  s m 2 ,  t = 7200 kg m 3 , E = 2.07  1011 N m 2 ,  = 0.24 , R0 = 0.3 m ,  f = 0.667 N  s m ,  f = 850 kg m , L = 12000 m ,  s 0 = 1 m s . Calculation time Tend = 691200 s . 2 3 The centerline of the pipeline is described by a fractional rational function: y = 40 (1 − 0.001x ) (1 + 10 x ) , − 6000  x  6000. −6 2 In numerical experiments, the following are found: displacements of the centerline, angular deformations of the walls, coordinates of the centerline x (t , s ) , y ( t , s ) , longitudinal displacements of the first approximation u1 . From (12) it follows that the physical meaning of  u1 R0 is warping of pipe cross sections. Warping of the cross-sections of a cylindrical tube were observed in the experiments of V.S. Vlasov [4]. The coordinates of the centerline at the beginning and end of the calculation are shown in Figure 1 (a). This illustrates the consistency of the numerical solution with the mechanics laws: the profile displacement is directed towards distributed load from the fluid flow. Figure 1 (b) shows the graph of  u1 , repeating the cross-sections warping. Warping has a maximum of about 0.003R0 . In another numerical experiment for a cubic parabola profile, the warping reached 0.02R0 . 4.2. The pipeline with singularity Physical and geometric parameters of the test problem: h = 5 mm, R0 = 47.5 mm, L1 = 250 mm, L2 = 250 mm. Material: S335JO steel,  t = 7800 kg m 3 , E = 210 GPa,  = 0.3 ,  t = 343 MPa. Rigid fixing conditions are imposed on the outer ends. Pressure p = 10 MPa is applied to the inner pipes surface. 177 Figure 1: (a) – profile coordinates at the beginning (dashed line) and at the end of calculation; (b) – longitudinal displacement in the first approximation. Mesh parameters are set to high precision. As a result, the stress distribution was found, see Figure 2. Figure 2 (a) shows the von Mises stresses. Figure 2 (b) shows a histogram of stress distribution by the number of mesh nodes. Limiting stress values: maximum stress 395.5 MPa (reached on the inner side in the reentrant corner of the domain), average stress 89.5 MPa, minimum stress 2.9 MPa. Figure 2: Stress field in the pipe: (a) – graphic isolines on the grid, MPa; (b) – histogram. 5. Conclusion A mathematical model of the pipeline as a moment shell of irregular geometry, with singularities caused by the domain geometry, has been built. The parameters that differentiate the various stress states of the pipeline are determined. A constructive algorithm for the asymptotic solution of the original three-dimensional boundary value problem for a pipeline is created, and on the basis of the obtained asymptotics, reduced mathematical models are constructed. The parameters of the problem are found for which various models are applicable. A numerical criterion has been determined that makes it possible to distinguish between two limiting states of the pipeline geometry. 178 Numerical experiments were performed. In the case (2), the results consistency of numerical experiments on the proposed mathematical model with the mechanics laws is shown. The existence of the pipe cross-sections warping is proved. In the case (3), the pipeline modeling algorithm in the FreeCAD package is created. The stress fields are found by the finite element method. The presence of a singularity in the stress field is proved. 6. Acknowledgements The work has been supported by the Russian Science Foundation grant № 21-11-00039, https://rscf.ru/en/project/21-11-00039/. Computational resources were provided by the Shared Services Center "Data Center of FEB RAS". 7. References [1] I. Towhata, Geotechnical Earthquake Engineering, Springer-Verlag, Berlin, Heidelberg, 2008. [2] A. L. Gol’denveizer, Theory of Elastic Thin Shells, Published for the American Society of Mechanical Engineers by Pergamon Press, 1961. [3] V. A. Rukavishnikov, O. P. Tkachenko, Dynamics of a fluid-filled curvilinear pipeline, Applied Mathematics and Mechanics 39 (2018) 905–922. doi:10.1007/s10483-018-2338-9. [4] V. S. Vlasov, Basic Differential Equations in General Theory of Elastic Shells, Technical Report NACA-TM-1241, National Advisory Committee for Aeronautics, Washington, 1951. URL: https://ntrs.nasa.gov/search.jsp?R=20050028489. [5] V. A. Rukavishnikov, O. P. Tkachenko, Numerical and asymptotic solution of the equations of propagation of hydroelastic vibrations in a curved pipe, Journal of Applied Mechanics and Technical Physics 41 (2000) 1102–1110. doi:10.1023/A:1026619009228. [6] L. I. Sedov, Mechanics of Continuous Media, World Scientific Publishing Company, Singapore, 1997. [7] V. A. Rukavishnikov, O. P. Tkachenko, Mathematical model of the pipeline with angular joint of elements, Mathematical Methods in the Applied Sciences 43 (2020) 7550–7568. doi:10.1002/mma.5751. [8] L. G. Loitsyanskii, Mechanics of Liquids and Gases, Pergamon Press, Oxford-New York, 1966. [9] V. A. Rukavishnikov, S. G. Nikolaev, On the R -generalized solution of the Lamé system with corner singularity, Doklady Mathematics 92 (2015) 421–423. doi:10.1134/S1064562415040080. [10] V. A. Rukavishnikov, On the uniqueness of the R -generalized solution of boundary value problems with noncoordinated degeneration of the initial data, Doklady Mathematics 63 (2001) 68–70. [11] V. A. Rukavishnikov, E. I. Rukavishnikova, The finite difference method for boundary value problem with singularity, Lecture Notes in Computer Science (including subseries Lecture Notes in Artificial Intelligence and Lecture Notes in Bioinformatics) 9045 (2015) 72–83. doi: 10.1007/978-3-319-20239-6_7. [12] V. A. Rukavishnikov, A. O. Mosolapov, E. I. Rukavishnikova, Weighted finite element method for elasticity problem with a crack, Computers and Structures 243 (2021) 106400. doi: 10.1016/j.compstruc.2020.106400. [13] B. Collette, FreeCAD for Inventors: Practical Examples and Clear Descriptions, Purple Squirrel Productions LLC, Gurley, USA, 2018. 179