Research Article Parametric instability analysis of the top-tensioned riser in consideration of complex pre-stress distribution Advances in Mechanical Engineering 218, Vol. 1(1) 1 11 Ó The Author(s) 218 DOI: 1.1177/168781417753894 journals.sagepub.com/home/ade eixin i 1,2 and uyun Chen 2 Abstract The parametric vibration instability of a riser is studied in consideration of a complex pre-stress distribution. Differential equations of the riser are derived according to Euler Bernoulli beam theory, and a method to solve the differential equations is proposed. With the parametric vibration of a top-tensioned riser as an example, the effects of the amplitude and direction of complex pre-stress on frequency, mode shapes, and instability characteristics are investigated. Results show that welding residual stress influences the dynamic response of the riser structure. A new approach to eliminate the complex loading of the riser is obtained. Keywords Riser structure, complex pre-stress, top-tensioned riser, parametric vibration, Mathieu instability Date received: 4 September 217; accepted: 12 December 217 Handling Editor: Jianqiao Ye Introduction With the growing demand for crude oil and gas in recent years, deep-water exploitation has become the focus of offshore equipment. The riser structure, as the key equipment linking the platform and wellhead at the sea base, has become a popular issue in engineering design. Aside from the influence of gravity, riser structures are subjected to wave- and current-induced loading and high-pressure oil and gas. Thus, several dynamic responses, such as parametric vibration, vortex-induced vibration (VIV), and collision, occur in riser structures. The longitudinal vibration (i.e. parametric vibration) of risers is caused by the heaving of a floating platform, which easily leads to the destruction of risers. Many studies that involved theoretical analysis, numerical calculation, and experimental investigation have been conducted to investigate the fundamental mechanism of parametric vibration. The non-linear resonance arising from parametric excitation problems has been discussed, and closed-form solutions for a riser have been obtained on the basis of first and second modes through extensive mathematical manipulations. 1 The parametric vibration caused by the wave-induced motions of a floating platform is practically important because it can destroy risers. 2 A stability analysis of the riser structure was performed in a previous study using yapunov stability theory, in which the top boundary controller is considered. 3,4 The finite element method has also been employed to analyze the influence of water depths, environmental conditions, and vessel motions under combined parametric and forcing excitations. 5 Structural natural frequency is an important structural dynamic property of riser structure. However, the natural frequencies of riser structure possess a low modal, and these 1 China Ship Development and Design Center, Wuhan, China 2 State Key aboratory of Ocean Engineering, Shanghai Jiao Tong University, Shanghai, China Corresponding author: uyun Chen, State Key aboratory of Ocean Engineering, Shanghai Jiao Tong University, Shanghai 224, China. Email: cluyun@163.com Creative Commons CC BY: This article is distributed under the terms of the Creative Commons Attribution 4. icense (http://www.creativecommons.org/licenses/by/4./) which permits any use, reproduction and distribution of the work without further permission provided the original work is attributed as specified on the SAGE and Open Access pages (https://us.sagepub.com/en-us/nam/ open-access-at-sage).
2 Advances in Mechanical Engineering frequencies are close to one another because of the structural stiffness problem in risers. A closed-form solution for the natural frequencies and associated mode shapes of axial loading has been deduced with Timoshenko beam theory. 6 The tension force in a riser structure consists of constant and varying parts. 7 According to Spark s theory, the influence on bending stiffness is defined as tension force; natural frequency and mode are deduced with the segmentation method. 8 Simple trigonometric functions are employed as approximations for vibration modes in Galerkin s method, in which the riser assumes a Bernoulli Euler beam model and small rotations. 9 If the frequency of external excitation and the natural frequency of a riser satisfy a certain condition, then the parametric instability problem is defined as Mathieu instability. Hsu 1 was one of the first researchers who analyzed parametric resonance in offshore cable applications and the instability region problem. The asymptotic solution of the Mathieu instability problem under stochastic parametric excitation and non-linear damping has been studied with the combination of linear and non-linear power-law damping. 11 Any possibly undesirable phenomena that occur in the transverse vibration of the riser structure caused by this fluctuation need to be investigated to develop a frequency domain method for linear systems with general time-varying parameters. 12 A suitable mathematical model to explore the stability of a submerged floating pipeline between two floating structures under vortex and parametric excitations has been presented and discussed. 13 Pre-stress (initial stress) often exists in complex structures and may be caused by welding residual stress, structural manufacturing defects, material thermal effects, static external loading, and so on. Prestress can resist or aid in structural deformation and alter the static and dynamic characteristics of structures. Pre-stress exerts a significant influence on local and global stiffness matrices, and the natural frequencies of a structure increase or decrease with pre-stress distribution. Welding residual stress must be considered when analyzing the influence of pre-stress on the dynamic characteristics of risers. The effects of uniform pre-stress distribution on natural frequencies and dynamic responses have been investigated. 14 With regard to uniform Euler Bernoulli beams under linearly varying fully tensile stress, a structure s natural frequencies may increase or decrease; parameters change the forbidden frequencies of the mechanical system after pre-stress is considered. 15 The vibration equation was modified in previous studies to analyze the influence of pre-stress, and the effects of variation in flow velocities and hydrostatic pressure on the dynamic behavior of fluid-conveying shells were defined as prestress. 16,17 Pre-stress includes linear and non-linear parts, and the superposition principle applies to the linear parts. 18 The effects of welding residual stress on the added virtual mass and the quality factor of the diaphragm have been presented. 19 In addition, the differential equation of the vibration of a cylindrical shell with welding residual stress has been derived, and a theoretical solution has been presented. 2 The analytic expression of the influence of complex pre-stress force on a riser structure was derived in a previous study, and the VIV response was compared. 21 The influence of complex pre-stress on the mechanical response of risers, especially on their structural dynamic characteristics, requires further study. However, only a few studies have considered the effects of complex pre-stress on the dynamic response of risers. The main objectives of this study are to investigate the influence of complex pre-stress (welding residual stress) on the natural frequency, modal shape, and size of instability domains for the parametric vibration problem of a top-tensioned riser (TTR). The developed analytical method can analyze the dynamic behavior of riser structures with/without pre-stress distribution, local area or overall pre-stress distribution, and even non-uniform pre-stress distribution. The outline of this article is as follows. A model of complex pre-stress force theory is presented in section Pre-stress force model of a beam structure. A brief description of governing equations for a riser structure with a complex pre-stress distribution is given in section Governing equation of a riser with complex pre-stress distribution. A vibration-free analytical solution of a riser structure with complex pre-stress is introduced in section Free vibration of the riser with complex pre-stress distribution. The parametric vibration instability problem of riser structures is presented in section Parametric vibration instability of the riser structure. A numerical analysis is implemented in section Numerical results and discussion, and the conclusions are presented in section Conclusion. Pre-stress force model of a beam structure Pre-stress force (initial stress) usually exists in continuum structures. The existence of pre-stress force significantly affects the local and global stiffness matrices of riser structures. Existing pre-stress force can resist or aid in the deformation of beam structures. The influence of complex pre-stress force on the static and dynamic characteristics of riser structures, especially on their structural dynamic characteristics, is worthy of studying. In this study, only welding residual stress, structural stress caused by deep-water pressure, and axial tension are discussed. These three types of prestress force can be defined as a complex pre-stress force that is unaffected by external dynamic excitation force.
i and Chen 3 If complex pre-stress force satisfies the linear superposition principle, then it can be expressed as follows s r = s T + s H + s R ð1þ where s r is the complex pre-stress force in the structure. In this study, only the axial direction of complex pre-stress force is considered. Then, complex pre-stress force can be written as s r = s r, z. s T is the pre-stress caused by axial riser tension force, s H is the structural stress caused by deep-water pressure, and s R is the welding residual stress. In riser structures, complex prestress force s r is a non-uniform distribution stress. Governing equation of a riser with complex pre-stress distribution In this section, the governing equation of riser structure with complex pre-stress distribution is established. Modeling of TTR To reveal the parametric vibration problem, a TTR in a complex ocean environment is discussed in consideration of the heaving of the platform and the motion of the tension ring, and the riser structure is subject to current. The riser can be regarded as a long, continuous, tubular member that is straight and vertical, and the boundary conditions at the two ends are known. The heaving of a floating platform, which induces axial tension fluctuation in the riser, is considered. The top of the riser is connected to the main body of the platform through a compensator, and the heave compensator can be simplified as an equivalent spring with stiffness K. Several hypotheses are defined as follows: (1) the material and mechanical properties are uniform along the overall riser structure; (2) the tension variation along the riser length varies linearly with depth; (3) cross-flow vibration is considered, whereas in-line vibration is excluded; and (4) the effect of shear strain is small and can be disregarded. In addition, the pipe wall behaves elastically, that is, no internal damping is considered. A rectangular Cartesian coordinate system is introduced to establish deformation in the riser structure. The riser moves on the plane, as shown in Figure 1. The sea surface is set as the origin of the coordinate system, the x-axis is parallel to the flow velocity, and the z-axis is measured from the top of the riser. The motion equations on the two principal vertical planes are identical and can be derived independently for each plane because of the symmetry of the riser cross section. The riser structure can move only on the xoy plane because the lateral deflection (w(z, t)) is considered small, and every cross section remains perpendicular to the axis. The riser structure can be modeled Figure 1. Mechanical model and reference frame of the riser. as a Euler Bernoulli beam. This equation disregards the effects of rotational inertia. Differential equations of the riser structure With the assumption on the riser structure, the governing motion equation of the lateral deflection w(z, t) of the riser structure by vibration can be written as EI 4 w(z, t) 4 w(z, t) T(z, t) +(m r + m f + m a ) 2 w(z, t) w(z, t) t 2 + c S = f (z, t) ð2þ t where the first item is the bending stiffness of the riser structure and the second term is the axial riser tension force. The third term is the inertia of the riser structure and includes riser structure, internal fluid, and fluid addition masses. f (z, t) is the external force, including hydrodynamic force, and so on. EI represents the bending stiffness of the riser structure, and T(z, t) stands for the effective axial tension of the riser. w(z, t) corresponds to the displacement vertical to the riser structure axis, z is the axial position, and t is the time parameter. In the partial differential equation, m r = 1 4 r sp(d 2 D2 1 ) is the mass per unit length of the riser, m f = 1 4 r f pd 2 1 is the mass per unit length of the internal fluid, and m a = 1 4 C ar w pd 2 is the mass per unit length of the added mass. The influence of the motion of internal fluid in the riser structure is omitted. D is the outer diameter of the riser, D 1 is the inner diameter
4 Advances in Mechanical Engineering of the riser, r s is the riser structural density, r f is water density, r w is internal fluid density, and C a is the added mass coefficient. c S is the damping parameter of the riser structure, and it is defined as in this study. Differential equations of the riser structure with complex pre-stress distribution According to Wu and Zhou, 19 axial riser tension force can be defined as a type of pre-stress force. Then, the partial differential equation of a deep-water riser structure with a complex pre-stress distribution can be written as follows EI 4 w(z, t) 4 s w(z, t) r, za +(m r + m f + m a ) 2 w(z, t) w(z, t) t 2 + c S = f (z, t) ð3þ t where s r, z is the complex pre-stress distribution in the riser stricture and can be written as s r, z = (½T U m z + Km sin (tt)š=a)+s R. In equation (3), A is the cross section of the riser structure; U m is the submerged weight of the riser per unit length, which is equal to m r + m f + m a ; K = U m =l is the stiffness of the heave compensator; parameters m and t correspond to the amplitude and frequency of the heaving of the platform, respectively; l depicts a parameter of the heave compensator; s R is caused by welding residual stress; T = Ag(r s r w )f top is the static tension on the top of the riser; stands for the length of the riser structure; f top is the top tension coefficient (defined as 1.3 in this study); and g corresponds to gravitational acceleration. Equation (3) is a modified function of the riser structure that considers complex pre-stress force distribution. By comparing this equation with the traditional differential equation of deep-water riser structures, the parameter of pre-stress =½s r, z S( w(z, t)=)š is obtained. This parameter can comprehensively describe the stress conditions of the riser structure. If axial tension and welding residual stress are defined, then the structural modes and structural modal of the riser and the solution of the dynamic response can be obtained. displacement at each end is zero, and the boundary conditions can be written as follows w(z, t)j z = =, 2 w(z, t) 2 =, w(z, t)j z = =, z = 2 w(z, t) 2 = ð4þ z = Modal decomposition is based on the assumption that the mode of the riser may be expressed as a sum of eigenmodes or eigenfunctions at any point in time. The solution to this equation can be obtained with a form of power series expansion. To obtain the solution of equation (3), Galerkin s procedure is employed. The solution of linearized equation (3) under the simple support boundary can be expressed as follows w(z, t)=f(z)e ivt = XN n = 1 B n sin (nbz)e ivt ð5þ where z is the axial coordinate, is the length of the riser, t is the time parameter, w n (t) is the modal weight function, and f n (z) is the mode shape function, where n = 1, 2, 3,... The mode shapes can be defined as sinusoid functions as follows: f n = sin (nbz), where b = p= and v is the angular frequency. The free vibration differential formulation of the riser with complex pre-stress can be obtained by introducing orthogonal series sin (jbz) and substituting equation (5) into equation (3). Using the Ð orthogonal property of the trigonometric series, sin (mpz=) sin (npz=)= exists because m 6¼ n integrates the function from z = to z =. Then, equation (3) can be written as follows EI(nb) 4 B n 2 ð s r, z A w sin (jbz)dz = mv 2 B n ð6þ where R = (2=) Ð ( =)(s r, za( w=)) sin (jbz)dz, which is defined as the integration item of the complex pre-stress and dynamic vibration displacement. Free vibration of the riser with complex pre-stress distribution The structural modes and modal problem of the riser with a complex pre-stress distribution are also discussed. Definition of boundary conditions The physical boundary condition at both ends of the riser can be modeled as simple support. The Solution of the beam structure The free vibration of the riser with a complex pre-stress distribution is analyzed, and the analytical solution for equation (6) is discussed as follows: 1. If the complex pre-stress satisfies the equation s r, z =, then equation (6) can be expressed as a vibration function of a simple beam. The solution equation can be written as EI(nb) 4 = rav 2.
i and Chen 5 2. If the complex pre-stress is in a uniform distribution form, then the pre-stress can be written as s r, z A = T, where A is the cross-sectional area. By substituting the complex pre-stress into equation (6), the solution equation can be written as EI(nb) 4 + T (nb) 2 = rav 2. If parameter s r, z is the tension stress, then the natural frequency of the beam structure increases. If parameter s r, z is a compressive stress, then the natural frequency of the beam structure decreases. 3. If the complex pre-stress is a function of variable z, then the complex pre-stress is a one-dimensional (1D) complex pre-stress distribution problem. Its distribution can be fitted by a trigonometric function and can be expressed as follows s r, z = s r, z cos (gbz) ð7þ where s r, z is the amplitude of the complex pre-stress and g is an integral number that is not less than 1. By separating the variables and substituting the complex pre-stress function into equation (6), function R can be obtained as follows R = 2s ð r, za = 2s r, za X N n = 1 w cos (gbz) B n (nb) ð sin (jbz)dzr ½cos (nbz) cos (gbz)š sin (jbz)dz ð8þ According to the trigonometric function, function R can be written as follows R = s r, za ð X N n = 1 B n (nb) fcos½(n + g)bzš + cos½(n g)bzšg sin (jbz)dz 8 = s < r, za X N ð B n (nb 2 )(n + g) sin½(n + g)bzšsin (jbz)dz : n = 1 9 ð + XN = B n (nb 2 )(n g) sin½(n g)bzšsin (jbz)dz ; n = 1 B n (nb 2 )(n g) ð sin½(n g)bzšsin (jbz)dzg ð9þ According to the orthogonal characteristic of the trigonometric function, the solution of equation (9) can be written as follows ð ð sin½(n + g)bzšsin (jbz)dz = =2 n + g = j n + g 6¼ j ð1aþ 8 < =2 n g = j sin½(n g)bzšsin (jbz)dz = =2 n g = j : jn gj 6¼ j ð1bþ By substituting equation (1) into equation (9), function R can be written as follows R = 8 >< >: s r, za 2 s r, za 2 P N g = 1 P N g = 1 B n g (nb 2 )(n g)+b n + g (nb 2 )(n + g) B jn gj (nb 2 )(n g)+b n + g (nb 2 )(n + g) n.g n\g ð11þ According to equation (11), N number functions exist and can be written in a matrix form as follows ( + R g )X = ð12þ where X = ½w 1 (z, t), w 2 (z, t),..., w N (z, t)š T and =diag fei(nb) 4 rav 2 g is a diagonal matrix. The element parameter in the diagonal matrix can be written as EI(nb) 4 rav 2. R g is a sparse matrix, and the value of function R is as follows If p = n and q = jn gj, then the following equation is deduced as 8 >< R pq = >: s r, z A (nb 2 )(n g), 2 s r, za (nb 2 )(n g), 2 n.g n\g ð13aþ If p = n and q = n + g, then the following equation is deduced as R pq = s r, za (nb 2 )(n + g) ð13bþ 2 4. If the distribution of the complex pre-stress has high complexity, then the complex pre-stress is a function of variable z and a 1D complex prestress distribution problem. The distribution of the complex pre-stress can be fitted by a trigonometric function and expressed as follows s r, z = XG g = 1 s r, zg cos (gbz) ð14þ where s r, zg (g = 1, 2,..., J 1, J) is the amplitude of the complex pre-stress force.
6 Advances in Mechanical Engineering By substituting the complex pre-stress function into equation (6), N number functions exist and can be written in a matrix form as follows ( + R)X = ð15þ where is a diagonal matrix. An increase in the series leads to a change in complex pre-stress influence matrix R because the complex pre-stress expression is highly complex. Matrix R is not a sparse matrix. In addition, the series of complex pre-stress follows the linear superposition principle and can be expressed as follows Modal analyses R = XJ R g g = 1 ð16þ According to the complex pre-stress function, the dynamic response function of the beam structure with complex pre-stress can be obtained because the structural complex pre-stress can be expressed as a trigonometric function. For a beam structure with/without a complex pre-stress distribution, the equations are linear, and the determinant factor to the characteristic equation is zero, which can be expressed as follows j + Rj = ð17þ In equation (17), if no complex pre-stress force distribution (R = ) exists, then equation (17) can be written as a classical elastic beam free vibration equation. If complex pre-stress distribution (R 6¼ ) exists, then the matrix + R is not diagonal. Therefore, the characteristic equation and structural mode of the free vibration beam with complex pre-stress need to be modified. Parametric vibration instability of the riser structure The diameter ratio of the riser is relatively large and causes the modal frequencies to be very close. If the vibration frequency of the structure is equal to half of the external excitation frequency, then parametric resonance could easily occur. Parametric vibration equation The parametric excitation problem is caused by the motion of the platform, which leads to the variation of the axial tension force of the riser. Parametric excitation depends on the axial tension force of the riser structure. The axial tension force of the riser includes static and dynamic axial tension forces. The amplitude and frequency of dynamic axial tension force depend on the motion of the platform heaving, which can be simplified as a harmonic force. External excitation, structure damping, and hydrodynamic damping are omitted. Then, the differential equations of the riser with complex pre-stress distribution are deduced as follows EI 4 w(z, t) 4 (T U m z + Km sin (tt)) w(z, t) + (m r + m f + m a ) 2 w(z, t) t 2 = ð18þ where T is the static axial tension force, U m z is the influence of gravity, Km is the amplitude of dynamic axial tension force, and t is the frequency of platform heaving. In equation (17), the solution form can be written as w(z, t)=f(z) q n (t)= P N n = 1 B n sin (nbz) q n (t), in which q n (t) is the unknown time function. Solution of the unstable regions According to equation (3), the Mathieu equation of the riser can be obtained as follows q(t)+v 2 Km cos (tt) nq(t) (m r + m f + m a ) ~ f mn q(t)= n =, 1, 2, 3,... ð19þ where ~ f mn = Ð u n u dz= Ð u2 n dz, v n is the Nth order modal frequency of the riser and t is the parametric excitation frequency. Parameter ~e = (Km=v 2 n (m r + m f + m a ))( Ð u n u dz= Ð u2 ndz) is defined as a parameter excitation value. If ~e\\1, then the solution of equation (19) for the analysis of the unstable region for the Mathieu equation is obtained. Thus, the value of ~e is important in dealing with the parametric instability problem. In equation (19), the riser will resonate if the frequency of platform heaving satisfies t =, 2v, 2v, which easily leads to the destruction of the riser. Then, the resonance t = 2v, 2v of the riser is defined as first-order resonance, that is, the main resonance. A large number of studies have indicated that the parameters of offshore platforms are unstable in the first order, and first-order resonance of the riser can be ensured. In this study, the first-order instability of the vertical tube structure is reviewed. Numerical results and discussion Model description The design parameters of the riser structure for the numerical analysis of parametric vibration are shown in Table 1. The effects of complex pre-stress on the riser s natural frequency and mode are compared.
i and Chen 7 Table 1. Design parameters of the model system. Parameter Symbol Values Unit ength 16 m Outer diameter D.3 m Thickness t.26 m Young s modulus E 2.1E11 Pa Material density r s 785 kg/m 3 Density of water r w 125 kg/m 3 Density of oil r f 8 kg/m 3 Equivalent coefficient l 1 m Top-tension coefficient f top 1.3 Addition mass coefficient C a 1. Stiffness of compensator K 32, N/m Heave amplitude m 5 m Heave frequency T s Figure 4. Distribution model III of welding residual stress. Figure 2. Distribution model I of welding residual stress. Figure 3. Distribution model II of welding residual stress. Modeling of welding residual stress To analyze the influence of welding residual stress on the dynamic characteristics of the riser, three types of welding residual stress distribution, namely, two types of tensile residual stress and a compressive residual stress, are compared. The three types of welding residual stress distribution are as follows: distribution model I is tensile stress; distribution model II is tensile stress, but its maximum amplitude is smaller than that of model I; and distribution model III is defined as compressive stress. The welding residual stress in the riser follows a nonuniform distribution. According to the common welding technology for risers, the peak value of stress is distributed periodically by a certain distance. In this study, every 2 m is defined as peak stress, and the welding residual stress is zero at the end of the riser structure. The welding residual stress in the riser can be fitted by a trigonometric series. With the characteristic of the periodic distribution of welding residual stress considered, five peaks are shown in figures to describe the welding residual stress distribution. Figure 2 shows distribution model I of welding residual stress; a positive value denotes tension stress. Figure 3 shows distribution model II of welding residual stress, and a positive value also denotes tension stress. In comparison with that in Figure 2, the peak value of welding residual stress in Figure 3 exhibits a decrease. Figure 4 shows distribution model III of welding residual stress. A negative value denotes compressive stress. In comparison with that in Figures 2 and 3, the welding residual stress in Figure 4 is compressive stress. The absolute value is similar to that in model I. Effect of welding residual stress on vibration frequency According to the function of welding residual stress, the complex pre-stress in the riser can be written as
8 Advances in Mechanical Engineering Table 2. Comparison of natural frequencies of the riser. Order Axial force, neglecting the sole weight Axial force and sole weight Welding residual stress mode I Welding residual stress mode II Welding residual stress mode III 1.296.14295.15871.15768.12998 2.4193.28597.31747.31538.266 3.629.42913.47635.4732.3935 4.8388.57249.63538.63119.5296 5 1.489.71613.79463.78939.65197 6 1.2591.861.95415.94787.78346 7 1.4696 1.45 1.114 1.167.91553 8 1.684 1.1493 1.2742 1.2658 1.482 9 1.8915 1.2947 1.4348 1.4254 1.1817 1 2.13 1.447 1.5959 1.5855 1.3159 s r, z =(½T U m z + Km sin (tt)š=a)+s R. The vibration characteristics of the riser can then be obtained. Constant axial loads. Constant axial loading is also considered. The complex pre-stress in the riser structure can be written as s r, z = T =A. The mode shapes of the riser structure can be written as EI(df 4 (z)=dz 4 ) T (df 2 (z)=dz 2 ) (m r + m f + m a )v 2 f (z)=. Using the method of separation of variables, the mode function of the vertical tube under fixed axial force can be calculated. Varying axial loads. Vertical pipe weight and internal tension are then considered. The complex pre-stress in the riser structure can be written as s r, z =(T W a z)=a. Using the method of separation of variables, the mode function of the vertical tube under fixed axial force can be calculated, and the modes are not standard sine functions. Varying axial loading and welding residual stress. In this case, welding residual stress, vertical pipe weight, and internal tension are considered. The values of welding residual stress are obtained from the measurement data of the riser structure or through a numerical simulation analysis. The cosine function is used to fit the welding residual stress parameter. The welding residual stress in the riser can be written as s r, z = P 2 g = 1 s r, zg cos (gbz), and the distribution of welding residual stress is shown in Figures 2 4. The complex pre-stress in the riser can be written as s r, z =(T U m z=a)+ P 2 g = 1 s r, zg cos (gbz). By substituting the complex pre-stress into equation (3), the differential formulation can be written as follows " # EI df4 (z) dz 4 T U m z + A X2 g = 1 s r, zg cos (gbz) df(z) + W a (m r + m f + m a )v 2 f(z)= dz df 2 (z) dz 2 Table 2 shows the natural frequency of the riser structure with variation from the first order to the tenth order. Table 2 shows that complex pre-stress exerts a significant effect on the natural frequency of the riser. If the welding residual stress is positive, then the natural frequency of the riser will increase. If the welding residual stress is negative, then the natural frequency of the riser will decrease. In addition, with the increase in the order, the influence of complex pre-stress on natural frequency becomes increasingly significant. These results show that the peak value of welding residual stress distribution is not too evident in natural frequency. A reason is that the range of the riser in the peak area is relatively narrow. According to Table 2, analyzing the influence of welding residual stress is necessary in riser design. Modes of the riser structure Figure 5 shows the first-, second-, third-, fourth-, fifth-, and sixth-order modes of the riser structure. The modes are not standard sine functions. As shown in Figure 5, the maximum amplitude of the mode is moved to the bottom when the sole weight is considered. Welding residual stress exerts a great effect on the mode shapes of the riser structure. The mode shapes become increasingly complex when welding residual stress is considered. The influence of welding residual stress on the riser structure vibration mode is evident. When welding residual stress exists, the mode shape of the riser becomes a rough curve with mutation. The mutation direction depends on the direction of welding residual stress. The reason is that the welding residual stress had changed the local stiffness of the riser. Instability analysis for the riser structure Stability analysis is important for a linear system under parametric excitation. The stability chart is a well-
i and Chen 9 Figure 5. First six mode shapes in different loading cases: (a) first mode, (b) second mode, (c) third mode, (d) fourth mode, (e) fifth mode, and (f) sixth mode. known tool to examine the properties of parametric instability. The chart consists of massive transition curves above which are unstable zones and below which are stable zones. Using the small parameter -P method, the boundary line expression of the stable and unstable regions of the standard Mathieu equation is obtained, as shown in Figure 6. The stability chart is obtained according to the relation between the stable region and unstable dividing line for p or 2p periodic solution. Considering engineering conditions, the motion cycle of platform heaving not less than 4 s, and a mode
1 Advances in Mechanical Engineering Figure 6. Comparison of stability charts for the riser structure (the shaded areas are unstable): (a) axial force with the sole weight neglected, (b) axial force and sole weight, (c) welding residual stress distribution mode I, (d) welding residual stress distribution mode II, and (e) welding residual stress distribution mode III. number of more than six orders can be neglected. In this study, only the first-order instability region of the riser is discussed, as shown in Figure 6. In the figure, d denotes the time parameter, and e denotes displacement: d =4v n =O 2 and e =2~ev 2 n =O 2. The stability charts are computed by considering the first five modes (N = 5). The narrow zones are related to either combined or high-order parametric resonance conditions. Several conclusions are obtained from Figure 6. The riser parametric excited instability area significantly
i and Chen 11 increases when the gravity influence is considered. The unstable region migrates when the welding residual stress is considered. Tensile stress causes the unstable region to migrate to the left direction, and the unstable region is slightly reduced. Compressive stress causes the unstable region to migrate to the right direction, and the unstable region is slightly increased. The reason for this change is that the existence of welding residual stress had changed the local and overall stiffness. Conclusion In this study, to investigate the influence of welding residual stress on the natural frequency, modal shape, and size of instability domains, a new approach to analyze the dynamic characteristics of the riser structure is proposed. A corresponding differential equation is established. The numerical results show that complex prestress force exerts a significant influence on parametric vibration. The distribution of complex pre-stress causes the resonance point of the riser to migrate, and the migration direction corresponds to the complex prestress direction. The mode shape is no longer a smooth curve, and the distortion direction depends on the amplitude of complex pre-stress. In addition, the unstable region moves when complex pre-stress is considered. In this study, pre-stress (welding residual stress) possesses periodic distribution. The aperiodicity problem of prestress distribution will be considered in the future. Declaration of conflicting interests The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article. Funding The author(s) disclosed receipt of the following financial support for the research, authorship, and/or publication of this article: This work was partially supported by the Fund of State Key aboratory of Ocean Engineering (grant no: 157). References 1. Chatjigeorgiou IK and Mavrakos SA. Nonlinear resonances of parametrically excited risers numerical and analytic investigation for O =2v 1. Comput Struct 25; 83: 56 573. 2. Patel MH and Park HI. Combined axial and lateral responses of tensioned buoyant platform tethers. Eng Struct 1995; 17: 687 695. 3. He W, Sun C and Ge SS. Top tension control of a flexible marine riser by using integral-barrier yapunov function. IEEE/ASME T Mech 215; 2: 497 55. 4. He W, He XY and Ge SS. Vibration control of flexible marine riser systems with input saturation. IEEE/ASME TMech216; 21: 254 265. 5. Park H and Jung DH. A finite element method for dynamic analysis of long slender marine structures under combined parametric and forcing excitations. Ocean Eng 22; 29: 1313 1325. 6. Chen YF, Chai YH, i X, et al. An extraction of the natural frequencies and mode shapes of marine risers by the method of differential transformation. Comput Struct 29; 87: 1384 1393. 7. Sparks CP. Transverse modal vibrations of vertical tensioned risers. Oil Gas Sci Technol 22; 57: 71 86. 8. Senjanovic I, justina AM and Parunov J. Natural vibration analysis of tensioned risers by segmentation method. Oil Gas Sci Technol 26; 61: 647 659. 9. Franzini GR and Mazzilli CEN. Non-linear reducedorder model for parametric excitation analysis of an immersed vertical slender rod. Int J Nonlin Mech 216; 8: 29 39. 1. Hsu CS. The response of a parametrically excited hanging string in fluid. J Sound Vib 1975; 39: 35 316. 11. Brouwers JJH. Asymptotic solutions for Mathieu instability under random parametric excitation and nonlinear damping. Physica D 211; 24: 99 1. 12. ei S, Zhang WS, in JH, et al. Frequency domain response of a parametrically excited riser under random wave forces. J Sound Vib 214; 333: 485 498. 13. Wang ZQ and Yang HZ. Parametric instability of a submerged floating pipeline between two floating structures under combined vortex excitations. Appl Ocean Res 216; 59: 265 273. 14. Hamed E and Frostig Y. Natural frequencies of bonded and unbonded prestressed beams prestress force effects. J Sound Vib 26; 295: 28 39. 15. Ashwear N and Eriksson A. Natural frequencies describe the pre-stress in tensegrity structures. Comput Struct 214; 138: 162 171. 16. Zhang Y, Gorman DG and Reese JM. Vibration of prestressed thin cylindrical shells conveying fluid. Thin Wall Struct 23; 41: 113 1127. 17. Zeinoddini M, Harding JE and Parke GAR. Dynamic behaviour of axially pre-loaded tubular steel members of offshore structures subjected to impact damage. Ocean Eng 1999; 26: 963 978. 18. Bakhtiari-Nejad F and Bideleh SMM. Nonlinear free vibration analysis of prestressed circular cylindrical shells on the Winkler/Pasternak foundation. Thin Wall Struct 212; 53: 26 39. 19. Wu HY and Zhou SJ. Dynamic response of sensor diaphragm with residual stress in contact with liquids. J Sound Vib 214; 333: 372 378. 2. iu Y and Chen Y. The effect of weld residual stress on the free vibrational characteristics of cylindrical shell through the analytical method. J Vibroeng 216; 18: 334 348. 21. Chen Y, i X and Yang N. Analysis of vortexinduced vibration of marine riser under the consideration of pre-stress force distribution. J Shanghai Jiaotong Univ 217; 51: 476 484 (in Chinese).