Research Article  Open Access
Fengde Wang, Wensheng Xiao, Yanan Yao, Qi Liu, Changjiang Li, "An Analytical Procedure to Predict Transverse Vibration Response of JackUp Riser under the Random Wave Load", Shock and Vibration, vol. 2020, Article ID 5072989, 9 pages, 2020. https://doi.org/10.1155/2020/5072989
An Analytical Procedure to Predict Transverse Vibration Response of JackUp Riser under the Random Wave Load
Abstract
Marine riser is a key equipment in offshore drilling operation, and failure of the riser can lead to drilling moratorium; in severe cases, it may cause oil and gas leaks. In this paper, the timedependent boundary conditions of the riser and the randomness of wave load are considered to improve the calculation efficiency and accuracy of the dynamic response of the jackup riser. Based on the Euler–Bernoulli beam theory, an analytical method to determine the response of the jackup riser subjected to the random wave load was established by the Mindlin–Goodman method in the frequency domain, and an experiment was carried out to verify it. The research shows that transverse dynamic response is the main component of the transverse response of the riser, and the method proposed is feasible to calculate the transverse response of the riser.
1. Introduction
In this study, the jackup riser is taken as the research object. As the main equipment of offshore oil exploration and development, the jackup platform plays a leading role in the continental shelf waters. Jackup platform which can decrease the lateral rigidity of the platform is designed to adapt to the growth of oil and gas demand in deeper water. Under the random environment loads, random dynamic response is produced by the platform. Although the marine riser is a key equipment in offshore drilling, it is a weak component. As the working depth of the platform increases, the dynamic response of the riser becomes more and more complicated, which may lead to operational accidents, cause huge economic losses, and even threaten the lives of operators and the marine environment. Therefore, it is necessary to study the dynamic response characteristics of the jackup riser. In the sea area with the depth of 100200 m, subsea blowout preventers (BOP) are used to drill well (Figure 1). The jackup riser can be modelled as an Euler beam because of its large slenderness ratio, the top tensioning system creates axial tension force on the riser, and the lateral vibration of the platform can be treated as the timedependent boundary conditions of the riser. Therefore, the jackup riser can be treated as an axialloaded Euler beam with timedependent boundary conditions.
For the transverse vibration of the marine riser, some scholars have studied the natural vibration characteristics and forced vibration response of the riser. Ignoring the variable tension force, Clauss et al. [1] and Nayfeh et al. [2] investigated the natural vibration characteristics of the marine pipes and beams. Based on the work by Nayfeh et al., Sinir [3] investigated the pseudononlinear dynamic characteristics of buckled pipes by considering the effect of internal velocity. MontoyaHernández et al. [4] developed a numerical algorithm to evaluate the natural frequency of risers under an internal multiphase flow. Klaycham et al. [5] derived the motion equation of a largedisplacement deepwater riser based on the Hamilton variational principle and studied the nonlinearfree vibration characteristics of the marine riser. By the Galerkin method, Zhou et al. [6] conducted research on the effects of the elastic constraints simulating damaged and undamaged boundaries on the natural frequencies and shapes of marine risers with a variable axial tension. Alfosail et al. [7] studied natural frequencies and mode shapes of the linear vibrations of inclined risers and proposed a numerical statespace approach to examine the natural frequencies and critical buckling limits of the marine riser [8]. Alfosail and Younis also carried out research on the twotoone internal resonance [9] and threetoone internal resonance of an inclined marine riser [10]; in their research, the excitation is formulated by the harmonic functions and the boundary conditions were treated as constants. Franzini and Mazzilli [11] studied the lateral motion of a slender and immersed rod subjected to harmonic and axial top motion with the reducedorder model (ROM). The dynamic response of marine risers under the harmonic excitation was also studied by Yin et al. [12]. According to the finitedifference approximation method, Wang et al. [13] put forward an idea of computing the dynamic response of risers under returning working conditions. Considering the interaction of the internal and external fluids on the marine riser, Wang et al. [14] studied lateral vibration characteristics of the riser. Reducedorder models were proposed by Vernizzi et al. [15] to investigate the parametric excitation of a vertical rod. Mao et al. [16] established a dynamic model to investigate the mechanical characteristics of a drilling riser and studied the offset of the platform as a series of determined values. For the nonlinear model, Lei et al. [17] studied the response of a deepwater riser to the platform motion, in which the motion of the platform was also treated as deterministic. Based on the Euler–Bernoulli theory, Liu et al. [18] and Zhao et al. [19] conducted research on the stabilization scheme for the flexible marine riser system, and the axial force is regarded as a constant. According to the linear wave theory, Fan et al. [20] and Wu et al. [21] investigated the vibration characteristics of the drilling riser. Do [22] proposed a constructive design of boundary controllers to stabilize lateral motion of flexible marine risers under random loads. The well posedness of the boundary conditions at the top end of the riser was studied by Guo et al. [23] with the assumption of constant axial force; and a threedimensional simulation of the curved riser was fulfilled by Zhu et al. [24] without considering the randomness of loads. Tang et al. [25] discussed the effects of tension ratio and platform deflection on the lateral displacement, bending moment, and the stress of the drilling riser. Recently, the penalty method was applied by Klaycham et al. [26] to investigate the nonlinear response of the marine riser.
The above literature review presents that most scholars used the numerical methods to investigate the lateral vibration response of marine risers. In their research, the average axial force and harmonic functions are used to establish the mathematical model of risers, and the quasistatic method is mainly used to solve the lateral vibration response of risers. However, neglecting the randomness of load and boundary conditions will reduce the calculation accuracy of the vibration response. Based on the current research achievements, the average axial force is also applied in this paper to establish the mathematical model of a riser; but we consider the timedependent boundary conditions of the riser and the randomness of the wave load in this paper. We propose an analytical method to determine the random lateral vibration response by the Mindlin–Goodman method. The research results can further improve the random vibration theory of beams; in terms of engineering application, this method can provide technical support to evaluate whether a special sea area and its sea conditions are available for the jackup, can adjust the installation sequence of the riser, and improve the service life of the marine riser system. The rest of this paper is organized as follows: an analytical procedure and an experiment are elaborated in Section 2. Then, a case study is presented in section 3. Finally, several conclusions are summarized in Section 4.
2. Analysis
2.1. Random Vibration Response of JackUp
The motion of a jackup is treated as the timedependent boundary condition of the riser in this research, thus it should be studied first. Figure 2 shows the sketch map of the jackup drilling system. According to the requirements of the China classification society, the jackup platform is modelled as a single degree of the freedom system, and the motion of the platform can be formulated by the following equation:where u(t) is the lateral displacement of the platform, m_{e} is the equivalent mass of the platform, c_{e} is the equivalent damping of the jackup leg, k_{e} is the bending stiffness of the jackup leg, and p(t) is the random wave load.
The frequency response function of the platform iswhere ω_{p} is the natural frequency of the platform. In this study, the Pierson–Moskowitz spectrum S() is selected to describe the wave:where H_{s} is the significant wave height of the wave, is the natural frequency of the wave, and the unit of the numerical is m^{2}/s^{4}. Based on the Morison equation and Borgman’s linearization method, the total wave force spectrum of one leg of a jackup can be obtained by the integral along the direction of the water depth, as shown in the following equation:where C_{D} is the drag force coefficient, C_{M} is the inertia force coefficient, k is the wavenumber of a significant wave, ρ is the density of seawater, is the gravitational acceleration, D_{o} is the outer diameter of the leg, z is the operating depth of the platform, and the parameter σ^{2} is the variance of horizontal speed of the water particle, which is the function of the height h. Considering the linearization of the drag force is an approximate treatment, and the variance can be obtained by the following equation:
Substituting equation (5) into equation (4), the approximate expression of the total wave force spectrum of one leg can be derived in the following equation:
According to the definition of power spectral density function and the relationship between autocorrelation function and power spectral density function, the total wave force spectrum of the platform S_{p}() can be obtained as follows:
By using equations (2) and (7), the power spectral density function of the random vibration S_{u}() of the platform can be obtained by the following equation:
2.2. Random Vibration Response of JackUp Riser with Subsea BOP
When a jackup platform uses subsea BOP to drill an exploration well, the top end of the riser is connected to the platform by a ball joint and a tension system, and the bottom end of the riser is connected to the subsea drilling system through a ball joint. In our research, several assumptions are applied as follows: (1) joints of the top end and the bottom end are modelled as a hinge; (2) the geometry characteristics and material properties of the riser are considered constant; (3) linearly varying axial force is replaced by the average axial force. Based on the above assumptions, the Euler–Bernoulli beam theory is adopted to establish the transverse motion equation of the riser as follows (without considering the internal damping of material):where y is the transverse deflection of the axis of the riser, t is the time, x is the coordinate measured along the beam axis, E is the modulus of elasticity, I is the area moment of inertia, T_{m} is the average axial tensile force, m is the mass per unit length, and c is the damping coefficient.
The timedependent boundary conditions of the riser at x = 0 and x = l (l is the length of the riser) are formulated as
By using the Mindlin–Goodman method, the lateral deformation of the riser is decomposed into quasistatic displacement y_{s}(x, t) and transverse dynamic displacement y_{d}(x, t) as follows [27]:
The coefficient r in equation (12) is decided by the boundary conditions of the riser, and in this research, r = 1. φ_{n} (x) is the vibration mode function, (x) is the static influence function, and q_{n} (t) is modal coordinate of ith mode:
The corresponding boundary conditions for the riser are expressed in the following equations:
Substituting equation (11) into equation (9) yields the following equation:
Based on equation (14), the first and second terms on the right side of equation (16) are both zero. Therefore, equation (16) is rewritten as
The substitution of equations (12) and (13) into equation (17) generates the following equation:
The transverse dynamic displacement of the riser satisfies the homogeneous boundary condition, which makes the mode shape function of the riser conform to the following equation:
Substituting equation (19) into equation (20) yields the following equation:
By applying the orthogonality conditions, equation (20) is decoupled to the following equation:where is the weight coefficient of generalized load. According to linear theory, the frequency response function of a modal coordinate versus u(t) can be obtained as follows:where ω is the natural frequency of the riser. By equation (23), the frequency response function between the transverse dynamic displacement and u(t) is generated in the following equation:
Then, the power spectral density functions of the transverse dynamic displacement can be formulated as
The cross spectral density function between u(t) and the transverse dynamic displacement is
By using the definition of autocorrelation function and the Wiener–Khinchin principle, the power spectral density functions of the quasistatic displacement is derived as follows:
Based on the above analysis, the power spectral density function of the beam’s transverse displacement can be obtained by the following equation:
According to the relationship between displacement and stress, the frequency response function of bending stress response σ(x, t) versus u(t) is written aswhere W is the bending modulus of the riser, and it can be calculated by equation (30). In equation (30), D is the outer diameter of the riser and d is the inner diameter of the riser:
On the basis of equations (25) and (29), the power spectral density functions of stress resulted from y_{d}(x, t) is
According to equations (27) and (29), the power spectral density functions of stress resulted from y_{s}(x, t) is
Combining equations (26) and (29), the cross power spectral density functions of the bending stress of the riser is obtained as
According to the summation formula of spectrum, the power spectral density function of the bending stress of the riser can be obtained by combining equations (31), (32), and (33):
2.3. Experimental Study
Based on the Mindlin–Goodman method, we propose an analytical method to obtain the random dynamic response of the Euler beam with axial load and timedependent boundary conditions. Subsequently, an experiment is carried out to determine the effectiveness of our method. Overall scheme of the experiment is (a) solve the dynamic stress response of the Euler beam under the sinusoidal constant frequency excitation based on our proposed method. (b) The stress response of the beam resulted from the above sinusoidal excitation is tested by the experimental system. (c) Make a contrastive analysis of the stress response data obtained by the proposed method and the experimental test. Figure 3 shows the experiment system.
2.3.1. Harmonic Excitation
In order to better observe the experimental phenomena and collect more accurate experimental data, we constantly adjust the frequency and amplitude of the sine constant frequency excitation and then choose the following expression E(t) as the excitation:
The beam in the experiment is fixed at x = 0 and hinged at x = l (l = 0.965 m), and the sinusoidal constant frequency excitation E(t) is given to the hinge support end. The coordinate of the test point is x = 0.03 m, and the transverse motion equation at x = 0.03 m of the beam is
The static influence function of the beam can be obtained by the initial parameter method, and the function of the beam’s pseudostatic displacement can be derived according to equation (12):
Then, the modal coordinate equation of the beam can be derived aswhere is the weight coefficients of generalized load. The modal coordinates of the dynamic displacement of the beam under sinusoidal constant frequency excitation is
The material of the beam is 45 carbon steel, and the diameter of the beam is 6 mm. Since the series in equations have a fast convergence rate, the first five orders of the series are enough for the analysis. Therefore, the first five order natural frequencies of the beam and weight coefficients are calculated first, which are listed in Tables 1 and 2. Subsequently, the first five modal coordinates of the transverse dynamic displacement of the beam are calculated by equation (39).


Based on the above analysis, the total transverse displacement of the beam in the experiment is
From the relationship between displacement and stress, the stress response of the beam can be obtained as
2.3.2. Experimental Results and Discussion
As the beam is cylindrical, its curved surface causes initial deformation of the strain gauge; therefore, this part of stress caused by the deformation should be subtracted from the original data. To ensure a fixation at the bottom, we installed a pin at the connection between the beam and the base, which generates bending stress; this part of stress caused by installation should be removed from the initial stress too. In Figure 4, the stress response curve obtained from theoretical method and experimental data are drawn, respectively. Due to the limitation of the experimental conditions, the foundation stiffness of the exciter and the beam is not large enough, so the extreme response amplitude of the beam is lower in the experiment.
The stress response curve obtained by experiment basically coincides with the theoretical calculation data in the main vibration response time period of the beam, which demonstrates that the proposed method is effective in solving the random response of a jackup riser. Engineers can use our method to conduct a preliminary study of the marine riser, and if they want to accurately obtain the dynamic response of the riser at a certain position, they need to establish the motion equation at that point.
3. Case Study
We take a typical riser with the length of 22.86 m and an external diameter of 1.372 m as the example, and the dynamic response of the riser with subsea BOP under the random wave load is investigated in our research. The parameters of the platform are m_{e} = 6.48 × 10^{6} kg, k_{e} = 4.71 × 10^{6} N/m, c_{e} = 8.77 × 10^{5} N m/s, D_{o} = 3.62 m, EI = 8.24 × 10^{8} N m^{2}, m = 461 kg/m, l = 110 m, T_{m} = 3 × 10^{3} N, H_{s} = 10 m, C_{D} = 2.0, C_{M} = 2.0, z = 100 m, ρ = 1025 kg/m^{3}, and = 9.8 m/s^{2}.
The first five order generalized load coefficients of the riser are calculated and listed in Table 3. The distribution of the stress response spectrums of the riser in the frequency domain and the spatial domain is plotted in Figures 5–7. The results show that the dynamical stress spectrum has two peak frequencies, which are the dominant frequency of the wave force spectrum and the fundamental frequency of the platform, and the total stress response spectrum is a narrow band spectrum; its dominating vibration frequency is the fundamental frequency of the platform. And, the dynamical displacement response of the jackup riser is the main response according to the comparative analysis between Figures 5 and 7.

Based on the above research results, the mean square deviation of the random displacement response and the stress response of the riser are calculated and depicted in Figures 8 and 9. Figure 8 shows that the minimum displacement response occurs in the middle of the riser. The stress response amplitude of the riser is the largest near its top end and bottom end, and the stress response is evenly distributed in the middle area of the riser system. The dynamical response results of our study are available to evaluate whether a certain sea area meets the drilling conditions, and it can also provide technical support for adjusting the installation sequence of risers and improve the service life of the jackup riser.
4. Conclusions
In this paper, a frequency domain analysis method is proposed to solve the random dynamic response of the jackup riser with subsea BOP and timedependent boundary conditions based on the Mindlin–Goodman method. This method can also be used to calculate the lateral response of the Euler beam with other boundary conditions. In addition, an experimental system was established and a response test experiment of the Euler beam was finished. The experimental results show that the method is effective for solving the transverse vibration response of the Euler beam. Through the case study, it is found that the lateral dynamic response of the jackup riser is dominated by the transverse dynamic responses, the minimum displacement response amplitude appears in the middle of the riser system, and the maximum stress response amplitude occurs at both ends of the riser system.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
The authors acknowledge the Ministry of Industry and Information Technology of the PR China for supporting this study through the project “Research on key technologies of integrated dismantling equipment for ultralarge offshore oilfield facilities” with the grant number 0300/05M1903005A.
References
 G. F. Clauss, H. Weede, and A. Saroukh, “Nonlinear static and dynamic analysis of marine pipelines during laying,” Ship Technology, vol. 38, no. 2, pp. 76–108, 1991. View at: Google Scholar
 A. H. Nayfeh, W. Kreider, and T. J. Anderson, “Investigation of natural frequencies and mode shapes of buckled beams,” AIAA Journal, vol. 33, no. 6, pp. 1121–1126, 1995. View at: Publisher Site  Google Scholar
 B. G. Sinir, “Pseudononlinear dynamic analysis of buckled pipes,” Journal of Fluids and Structures, vol. 37, pp. 151–170, 2013. View at: Publisher Site  Google Scholar
 D. J. MontoyaHernández, A. O. VázquezHernández, R. Cuamatzi, and M. A. Hernandez, “Natural frequency analysis of a marine riser considering multiphase internal flow behavior,” Ocean Engineering, vol. 92, pp. 103–113, 2014. View at: Publisher Site  Google Scholar
 K. Klaycham, C. Athisakul, and S. Chucheepsakul, “Nonlinear vibration of marine riser with large displacement,” Journal of Marine Science and Technology, vol. 22, no. 2, pp. 361–375, 2017. View at: Publisher Site  Google Scholar
 X. Zhou, M. Duan, and J. J. Granados, “Precise integration method for natural frequencies and mode shapes of ocean risers with elastic boundary conditions,” Applied Mathematical Modelling, vol. 61, pp. 709–725, 2018. View at: Publisher Site  Google Scholar
 F. K. Alfosail, A. H. Nayfeh, and M. I. Younis, “Natural frequencies and mode shapes of statically deformed inclined risers,” International Journal of Nonlinear Mechanics, vol. 94, pp. 12–19, 2017. View at: Publisher Site  Google Scholar
 F. K. Alfosail, A. H. Nayfeh, and M. I. Younis, “A state space approach for the eigenvalue problem of marine risers,” Meccanica, vol. 53, no. 45, pp. 747–757, 2018. View at: Publisher Site  Google Scholar
 F. K. Alfosail and M. I. Younis, “Twotoone internal resonance of an inclined marine riser under harmonic excitations,” Nonlinear Dynamics, vol. 95, no. 2, pp. 1301–1321, 2019. View at: Publisher Site  Google Scholar
 F. K. Alfosail and M. I. Younis, “Threetoone internal resonance of inclined marine riser,” International Journal of Nonlinear Mechanics, vol. 109, pp. 107–117, 2019. View at: Publisher Site  Google Scholar
 G. R. Franzini and C. E. N. Mazzilli, “Nonlinear reducedorder model for parametric excitation analysis of an immersed vertical slender rod,” International Journal of Nonlinear Mechanics, vol. 80, pp. 29–39, 2016. View at: Publisher Site  Google Scholar
 D. Yin, E. Passano, H. Lie et al., “Experimental and numerical study of a top tensioned riser subjected to vessel motion,” Ocean Engineering, vol. 171, pp. 565–574, 2019. View at: Publisher Site  Google Scholar
 S.w. Wang, X.s. Xu, B.h. Yao, and L. Lian, “A finite difference approximation for dynamic calculation of vertical free hanging slender risers in reentry application,” China Ocean Engineering, vol. 26, no. 4, pp. 637–652, 2012. View at: Publisher Site  Google Scholar
 Y. Wang, D. Gao, and J. Fang, “Study on lateral vibration analysis of marine riser in installationvia variational approach,” Journal of Natural Gas Science and Engineering, vol. 22, pp. 523–529, 2015. View at: Publisher Site  Google Scholar
 G. J. Vernizzi, G. R. Franzini, and S. Lenci, “Reducedorder models for the analysis of a vertical rod under parametric excitation,” International Journal of Mechanical Sciences, vol. 163, p. 105122, 2019. View at: Publisher Site  Google Scholar
 L. Mao, Q. Liu, S. Zhou, G. Wang, and Q. Fu, “Deep water drilling riser mechanical behavior analysis considering actual riser string configuration,” Journal of Natural Gas Science and Engineering, vol. 33, pp. 240–254, 2016. View at: Publisher Site  Google Scholar
 S. Lei, X. Y. Zheng, and D. Kennedy, “Dynamic response of a deepwater riser subjected to combined axial and transverse excitation by the nonlinear coupled model,” International Journal of Nonlinear Mechanics, vol. 97, pp. 68–77, 2017. View at: Publisher Site  Google Scholar
 Y. Liu and F. Guo, “Output feedback boundary control of a flexible marine riser system,” Journal of Vibration and Control, vol. 24, no. 16, pp. 3617–3630, 2017. View at: Publisher Site  Google Scholar
 Z. Zhao, Y. Liu, and F. Guo, “Robust output feedback stabilization for a flexible marine riser system,” ISA Transactions, vol. 78, pp. 130–140, 2018. View at: Publisher Site  Google Scholar
 H. Fan, C. Li, Z. Wang, L. Xu, Y. Wang, and X. Feng, “Dynamic analysis of a hangoff drilling riser considering internal solitary wave and vessel motion,” Journal of Natural Gas Science and Engineering, vol. 37, pp. 512–522, 2017. View at: Publisher Site  Google Scholar
 Z. Wu, C. Xie, G. Mei, and H. Dong, “Dynamic analysis of parametrically excited marine riser under simultaneous stochastic waves and vortex,” Advances in Structural Engineering, vol. 22, no. 1, pp. 268–283, 2019. View at: Publisher Site  Google Scholar
 K. D. Do, “Boundary control of transverse motion of flexible marine risers under stochastic loads,” Ocean Engineering, vol. 155, pp. 156–172, 2018. View at: Publisher Site  Google Scholar
 F. Guo, Y. Liu, Y. Wu, and F. Luo, “Observerbased backstepping boundary control for a flexible riser system,” Mechanical Systems and Signal Processing, vol. 111, pp. 314–330, 2018. View at: Publisher Site  Google Scholar
 H.j. Zhu and P.z. Lin, “Numerical simulation of the vortexinduced vibration of a curved flexible riser in shear flow,” China Ocean Engineering, vol. 32, no. 3, pp. 301–311, 2018. View at: Publisher Site  Google Scholar
 L. Tang, Z. Huang, X. Zhu, Y. Zhou, and B. Li, “Investigation of the mechanical response of a deepwater drilling riser to ocean currents and waves,” Advances in Mechanical Engineering, vol. 11, no. 1, pp. 1–11, 2019. View at: Publisher Site  Google Scholar
 K. Klaycham, C. Athisakul, and S. Chucheepsakul, “Nonlinear response of marine riser with large displacement excited by topend vessel motion using penalty method,” International Journal of Structural Stability and Dynamics, vol. 25, Article ID 2050052, 2020. View at: Publisher Site  Google Scholar
 R. W. Clough and J. Penzien, Dynamics of Structures, McGrawHill, New York, NY, USA, 3rd edition, 1995.
Copyright
Copyright © 2020 Fengde Wang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.