Research Article  Open Access
Kanhua Su, Stephen Butt, Jianming Yang, Hongyuan Qiu, "Coupled Dynamic Analysis for the RiserConductor of Deepwater Surface BOP Drilling System", Shock and Vibration, vol. 2018, Article ID 6568537, 15 pages, 2018. https://doi.org/10.1155/2018/6568537
Coupled Dynamic Analysis for the RiserConductor of Deepwater Surface BOP Drilling System
Abstract
Deepwater surface BOP (surface blowout prevention, SBOP) drilling differs from conventional riser drilling system. To analyze the dynamic response of this system, the riserconductor was considered as a beam with varied crosssections subjected to loads throughout its length; then an equation of motion and free vibration of the riserconductor string for SBOP was developed. The finite difference method was used to solve the equation of motion in time domain and a semianalytical approach based on the concept of section division and continuation was proposed to analyze free vibration. Case simulation results show that the method established for SBOP system natural frequency analysis is reasonable. The mode shapes of the riserconductor are different between coupled and decoupled methods. The soil types surrounding the conductor under mudline have tiny effect on the natural frequency. Given that some papers have discussed the response of the SBOP riser, this work focused on the comparison of the dynamic responses on the wellhead and conductor with variable conditions. The dynamic lateral displacement, the bending moment, and the parameters’ sensitivity of the wellhead and the conductor were analyzed.
1. Introduction
Several operators have developed surface BOP (surface blowout prevention, SBOP) drilling technology for deepwater drilling. SBOP drilling differs from a conventional riser drilling system: the BOP stack is located at the surface below the drill floor of the platform, not at the seabed. Another key difference is that the riser of the SBOP drilling system is designed to contain wellbore pressure, whereas a conventional drilling riser does not contain pressure. The SBOP drilling system can use smaller 2nd or 3rdgeneration semisubmersible rigs for operation. Moreover, it has illustrated a considerable amount of day rate saving over traditional drilling methods using subsea BOP [1]. As presented in the IADC (International Association of Drilling Contractors) guidelines for SBOP from floating MODUs [2], the key components of the SBOP drilling system include the surface BOP stack, upper transition joint, casing riser, lower transition joint, seabed isolation device (SID), subsea wellhead, conductor, and casing strings.
Surface BOP is not a new concept, but it was extended into deepwater only a few years ago. Based on successful drilling campaigns in Asia, Shell extended its SBOP technology to the more demanding offshore operations in Brazil with the implementation of a SID. Unocal and Transocean pioneered the application of the SBOP from floating drilling units, which began in early 1996 in the relatively benign environment of Southeast Asia (Kozicz, 2006) [3]. The deepwater well 1SHEL14RJS (block BMC10 in the Campos Basin, offshore Brazil) was successfully drilled in water depth of 2887 m by Shell in 2003. The 339.7 mm size, P110 grade, and 35 Mpa casing riser was adopted in this well (Brander et al., 2004) [4]. Total used the SBOP technique in 2000 m water depth on its Donggala block in Indonesia. For this SBOP drilling system, a 339.7 mm casing riser complete with specially designed fatigue enhanced connections was used (Simondin et al., 2004) [5]. In 2008, SBOP technology was implemented for both drilling and completion operations for the Parque das Conchas deepwater development in block BC10 offshore of Brazil. A dedicated 406.4 mm, X80 grade, and 41 Mpa high pressure riser replaced the marine riser in the SBOP system (Tarr et al., 2009) [6]. A compact deepwater drilling ship for SBOP operations was built (Claassen et al., 2010) [7].
Typical configurations of the casing riser include diameters of 273.1 mm, 339.7 mm, and 406.4 mm. At the top and bottom of the casing riser, heavy walled transition joints are required to distribute stresses [2]. Due to the small size of the casing riser, a 762.0 mm or smaller conductor is adopted to be jetted into the seabed as a support for the wellhead and other casing strings. Similar to a conventional riser, the high pressure casing riser is under complex forces imposed by vessel motions, waves, and current and soil interaction, which are transferred to the subsea wellhead and conductor.
Specialized computer programs are generally used to predict conventional riser’s behavior under the designed conditions [8]. The motion equation of the riser is usually converted to a system of finite length elements using either a finite difference or a finite element technique. Botke (1975) used a derivation of the riser equations and finite difference method (FDM) of solution [9]. Gardner and Kotch (1976) described finite element method (FEM) applied to the riser [10]. Subsequently, many researchers have studied the dynamics of the deepwater drilling riser.
For the SBOP drilling system, the design loads include bending loads, coming from the riser, stress joint and wellhead above the mudline, and soil reaction below the mudline. IADC identified that the casing riser and conductor analysis should be conducted in a coupled manner [2]. As the loads on the bottom of the riser will be transferred to the wellhead and conductor, the coupled concept of riser and conductor is concerned [11]. King et al. (1993) developed a new approach to analyze the behavior of a drilling riser and conductor as a complete entity, and the FEM was used in his research [12]. Although FEM is able to represent more details of the riser including the connection points, the mesh should be carefully implemented with good modeling practices. Su et al. (2008) studied the stability of the subsea wellhead and the conductor bearing capacity in deepwater drilling and proposed a wellhead stability analysis method with FDM [13]. Yan et al. (2015) analyzed the wellhead stability problem in deepwater drilling using the pile element and nonlinear spring element of ANSYS [14]. Results show that the coupled method to analyze the conventional riserconductor system is reasonable.
However, the casing riser analysis for the SBOP drilling system is conducted in accordance with API RP 16Q for the conventional riser. Some literatures discussed the dynamic behavior of the SBOP drilling riser. Morooka et al. (2008) presented a numerical simulation to estimate the riser behavior for a drilling system with surface BOP, but the research did not consider the coupling effect between the wellhead and conductor [15]. Dib et al. (2009) analyzed the fatigue lifetime of the SBOP riser in various modes of operation and also ignored the effect of the conductor [16]. For the SBOP riserconductor system, as the transition joints replace the flex/ball joint, this allows for the consideration of the riserconductor as a beam with varied crosssections subjected to loads throughout its length. Therefore, it is feasible to model the system considering the coupling between the riser and the conductor for the SBOP system and to solve the equations with FDM and other semianalytical methods.
In this paper, a coupled timedomain dynamic FDM method for the riserconductor of the SBOP system is derived, and a semianalytic method is developed for the free vibration of riserconductor. These methods are more convenient for the analysis coupled with the riser, SID, wellhead, and conductor of the deepwater SBOP drilling system.
2. Equation of Motion of RiserConductor in Transverse Vibration for SBOP Drilling
As the transition joints connect to the surface BOP and SID, there is no rotary joint on the riserconductor for the SBOP system. Thus, the riserconductor can be considered as a beam with varied crosssections subjected to loads. The riserconductor of the SBOP system is modeled as a variable section EulerBernoulli beam undergoing transverse vibration under axial force, as is shown in Figure 1(a). The forces acting on an element of the riserconductor system of length dx are shown in Figures 1(b) and 1(c).
(a)
(b)
(c)
2.1. Equation of Motion of RiserConductor System
By assuming (1) the riserconductor is a EulerBernoulli beam, (2) the riser and casing string are both linear elastic, (3) the drilling string has no effect on its bending rigidity, and (4) the vessel, wave, current, and riser all move in a plane, then the riser can be modeled as a beam subjected to loads throughout its length with boundary conditions at the top and bottom ends.
According to the EulerBernoulli theory [17, 18], the equation of motion for the forced lateral vibration of the riserconductor system can be obtained as inwhere is the lateral displacement of the riserconductor system, m; is the depth coordinate as shown, m; is the riserconductor system mass per unit length, kg; is the riserconductor system damping coefficient, N·s/m; is the riserconductor system stiffness, N/s; is Young’s modulus of the riserconductor system, N/m^{2}; is the area moment of inertia of crosssection, m^{4}; is the effective tension or compressive force on the riserconductor, N; and is the force on the riserconductor system, N.
The external forces on the riser, SID, and wellhead can be computed using the Morison equation (American Petroleum Institute, 2001) [8]; this equation has been widely used in riser dynamic analysis. As the riser transfers its forces to the conductor and casing strings, the lateral force on the conductor and casing strings is zero.
The bending stiffness of the riser can be calculated easily; however, the strings under the mudline are much more complicated containing the conductor, cement, and surface casing. Su et al. (2008) described a method to obtain the equivalent bending rigidity [13]; the equivalent bending rigidity of the casing string can be determined bywhere is the elastic modulus of steel, Pa; is the elastic modulus of the cement sheath between the conductor and the surface casing, Pa; , are the outer and inside diameters, respectively, of the riser, m; are the outer and inside diameters of the transition joint, m; and are the equivalent outer and inside diameters of the SID, m; , are the outer and inside diameters of the conductor, m; , are the outer and inside diameters of the surface casing string, m; , , , , are the lengths of the upper transition joint, riser, lower transition joint, SID, and the conductor above the mudline, m; and , are the lengths of the conductor under the mudline and surface casing string, respectively, m.
The riser mass per unit length should include the mass of the riser itself and the internal drilling mud [8]. The conductor and casing string mass per unit length include the mass of the conductor, the surface casing, and the cement sheath between them. The axial force of the riserconductor can be obtained as in (3). For the riser, the axial force is its effective tension [8, 13], while, in general, the axial force on the SID may be tensile and the force on the conductor and casing string is compressive.where is the vertical force along the riser, N; is the vertical force along the SID, conductor, and surface casing string, N; is the top tension of riser, N; is the riser density, kg/m^{3}; is the drilling mud density, kg/m^{3}; is the seawater density, kg/m^{3}; is the conductor density, kg/m^{3}; is the surface casing string density, kg/m^{3}; is the cement sheath density, kg/m^{3}; is the acceleration of gravity, m/s^{2}; is the weight per length of the casing string, N/m; and is the outside friction of the conductor and casing string, N.
The stiffness of the conductor and casing string can be derived from py curves under period loads according to the actual soil considerations [19]. The damping constant of the conductor and casing string can be expressed with the material damping and the radiation damping (Gazetas and Dobry, 1984) [20].
2.2. Boundary and Initial Conditions
As (1) is a fourthorder equation, four boundary conditions are needed. For the riser, the lateral displacement and bending moment of the upper transition joint are taken as two boundary conditions [21]. For the conductor and casing string, the shear force and the bending moment at their bottom are two lower boundary conditions and can be assumed to be zero. These four boundary conditions can be represented by
The initial condition of the equation of motion iswhere is the bending moment, N·m^{2}; is the shear force, N; is the lateral displacement of the static riserconductor system, m; is the upper rotational stiffness of the transition joint, N·m/rad; and is the upper rotation angle of the transition joint, rad.
is the horizontal deviation of drilling platform motion from its initial location, and it sums the mean offset and the platform drift responding to random waves (Sexton and Agbezuge, 1976) [22]. It can be expressed as follows: where is the mean offset of the platform, m; is the drift amplitude of the platform, m; is the drift period of the platform, s; and , , , , are the wave amplitude, wave number, circular frequency, initial phase, and phase difference of the wave n, respectively. These parameters can be obtained by wave theory and response curves in [21, 22].
2.3. Equation Solution Using FDM
It is difficult to solve the equations analytically; therefore, numerical simulation with the finite difference method was adopted in this paper.
The riserconductor string is divided into equal portions and the length of each section is . By using the threepoint difference format to replace the first and secondorder derivative schemes, the fivepoint difference format takes place of the fourthorder derivative scheme subsequently in (1). Then the finite differential equations of the riserconductor string can be obtained, which can be shown as follows:
According to the difference scheme, the differential equations of the upper boundary condition are expressed as in (8) and the lower boundary condition are shown as in (9):
When , removing the time items in the equations above, the initial conditions can be obtained by the difference equations. The static lateral deformation of the riserconductor system can be easily solved using a matrix or the Glesser method.
Starting from the initial conditions, the responses at a series of discrete time instants can be obtained through direct integration. MATLAB was employed to solve the model by time step. Through iterative calculation, the displacement, offset angle, bending moment, shear force, and soil reaction force at each node and any time were calculated.
3. Free Vibration Equations and Solution for SBOP RiserConductor System
For the free vibration of the riserconductor system, (1) reduces to
Assuming that the system is a uniform beam, (11) reduces to
The solution of (12) can be calculated according to the beam theory; then the natural frequency and the corresponding natural mode shape of the riserconductor system can be expressed as follows: where , , , and are constants that can be found from the initial conditions; is the natural frequency; and is the corresponding natural mode shape of the system.
The riserconductor of the deepwater SBOP system consists of several sections of different diameters shown in Figure 2; therefore, (11) cannot be directly used to solve the problem. Based on the concept of section division and continuation (Cui et al., 2012) [23], a semianalytical approach for analyzing free vibration of the SBOP riserconductor with variable crosssection is proposed. However, each section of the conductor system is with constant crosssection and can be treated as a uniform beam. So the natural frequency and the mode shape of each segment can be solved with (13).
For the segment , the natural mode shape is
Let .
Then (14) becomes
Therefore, the th natural mode shape of the segment is expressed as follows:
Since the deflection, slope, moment, and shear force of the segment and the th segment at node are equal, assuming ,
By substituting (15) and (16) into (17), the following is obtained:whereLetThen (18) becomesFrom (21),where
As (23) is the function of the natural frequency of the riserconductor, the relationship of the undetermined coefficient between the segment and the segment is established.
Assuming , follows (4); from the boundary conditions, we get (24):
Solving (24) using the secant iteration method by MATLAB to obtain the natural frequencies , the iteration formation is shown as in
Then, (22) is used to get the constants , and the corresponding natural mode shapes of the riserconductor.
4. Case Study and Discussion
A case study with the parameters given in Table 1 is carried out.

Seabed soil conditions vary substantially around the world. However, to simplify the calculation process and compare the results, the soil type below the mudline 0–60 m is assumed as all clay or sand layer. The six soil type properties are listed in Table 2.

4.1. Natural Frequencies and Mode Shapes of the SBOP RiserConductor System
The natural frequencies with the data given in Tables 1 and 2 with the proposed method are listed in Table 3. And Figure 3 shows the first six mode shapes of the riserconductor of the SBOP drilling system. As seen especially in Figure 3, the amplitudes of the mode shapes on the SID, wellhead, and conductor are very small because the soil reaction causes stiffness in these sections.

The mode shape comparison results for 4 situations (, coupled; , coupled; , coupled; , decoupled) are shown in Figure 4. From these figures, it can be observed that the mode shapes have some difference between the coupled and decoupled method; however, there is no obvious difference when the mode number is greater than 3. In addition, axial force has great effect on the mode shape and natural frequency, because axial force directly affects the stiffness of bending.
Although the TTR has little effect on the mode shape, its effect on the natural frequencies for the modes is obvious as shown in Figure 5. It can be seen that with the increase of TTR the natural frequency of the riserconductor increases as well. However, it would not have an effect on the natural frequency much for the SBOP riserconductor. It was also observed that the soil types surrounding the conductor under the mudline have very tiny effect on the natural frequency for the riserconductor for SBOP drilling system.
4.2. The Dynamic Response of the SBOP RiserConductor
To analyze the dynamic response of the riserconductor for the SBOP drilling system, lateral displacement, bending moment, and soil reactions at the different positions of the riserconductor string are compared. Given that some papers have discussed the response of the SBOP riser, this work focuses on the comparison of the dynamic responses on the wellhead and the conductor with variable conditions.
4.2.1. Dynamic Response of the Platform
The PM spectrum has been employed to calculate the motion of the platform and the simulation results are shown in Figure 6. From this figure, the motion amplitude of the platform is relatively small without considering offset and drift. This motion considering more conditions will obviously cause the riser response.
4.2.2. Lateral Displacement at Different Positions of the RiserConductor String
For the SBOP drilling system, some key joint points are critical to the drilling operation. This work focuses on 4 positions on the riser, namely, the bottom of the upper transition joint (USJ), the elevation 300 m under the mean water level (MWL), the elevation 800 m under the mean water level (MWL), and the top of the lower transition (stress) joint (LSJ), and 4 positions on the conductor, including the subsea wellhead (WH), the mudline (ML), −5 m under the ML, and −10 m under the ML.
The first four pictures in Figure 7 show the lateral displacements at different locations on the riser and the last four pictures show the displacements on the wellhead or the conductor. The lateral displacement amplitude at the elevation of 800 m of the riser under the MWL is the highest in the first four pictures; the displacement amplitude at the wellhead is more than other places on the conductor under the ML. Their periods are the same in the time domain, following the platform motion.
4.2.3. Dynamic Bending Moment on the Riser, the Wellhead, and the Conductor
In Figure 8, the bending moment’s variations at different positions on the riser, the wellhead, and the conductor are presented. The largest bending moment of the riser focuses on the LSJ, and the moment becomes smaller at the place of the conductor under the ML −5 m.
4.2.4. MomentDisplacement Curves
Comparing the bending moment and the lateral displacement in the same position simultaneously in Figure 9, the change of displacement is found to have a certain delay, compared with the bending moment on the conductor, and there is no delay on the riser. The reason is that the nonlinear soil reaction acted on the conductor under the mudline.
It is also found that the displacement varies more at the bottom of the USJ than at the top of the LSJ. The displacementmoment curves change on the negative axis under the mudline −10 m, because the displacement of the conductor under a certain depth will be negative or equal zero.
4.2.5. Deformation and Stresses of the Conductor in a Drift Period
The deformation and stress of the casing string below the mudline and the soil reaction force around the casing string are analyzed in a slowdrift period of a drilling platform. When the vibration reaches the steady state, 8 time points (t = 200 s, 232.5 s, 260.5 s, 297.4 s, 330.0 s, 362.4 s, 394.9 s, and 427.3 s) are taken in one period. The lateral displacement, bending moment, and soil reaction along the conductor are shown in Figure 10.
It can be seen from Figure 10 that the deformation, stress, and surrounding soil reaction of the casing string change with time in the slowdrift period. In a period, the maximum lateral displacement is along the column; the relative changes in bending moment and maximum soil reaction are 7.7%, 6.6%, and 3.2%, respectively.
4.2.6. Sensitivity Analysis of Parameter Changes to Wellhead and Conductor
(1) Platform Offset and Conductor Size. The different platform offsets of 0%, 1%, 2%, 3%, 4%, 5%, 6%, 7%, and 8% water depth and different conductor sizes (OD = 762.0 mm, 660.4 mm, and 508.0 mm, no conductor) are chosen to compare the results, as shown in Figure 11. The displacement and the bending moment are raised with the offset increase and those have a rapid increase in the condition when there is no conductor.
(2) Platform Surge Amplitude and Riser Size. Figure 12 displays the sensitivity of the lateral displacement and the bending moment to the variable surge amplitudes of platform and the riser’s outer diameters. In Figure 12, the lateral displacement and the bending moment of the wellhead increase with the greater surge amplitude. The 406.4 mm diameter riser will transfer more force to the wellhead and the conductor, so the 339.7 mm diameter riser is the best for the SBOP drilling system.
(3) Soil Properties. As the seabed soil supports the conductor, the soil property is very critical for the transverse stress and deformation of the conductor and the wellhead. Figure 13 compares the displacement and the bending moment of the wellhead and the position of the conductor at the mudline and 5 m below the mudline, using 6 soil types (see Table 2). The parameter of undrained shear strength (USS) is more sensitive to clay than the submerged unit weight. A lower USS of clay will cause weak conductor support. The parameter of the angle of internal friction is also more sensitive to sand than the submerged unit weight. Also, the bending moment does not change more than the lateral displacement of the wellhead and the conductor.
5. Conclusion
For the SBOP drilling system, the riserconductor is considered as combination sections with different crosssections subjected to loads throughout its length, and a FDM solution is derived in time domain. Results show that the displacement amplitude at the wellhead is more than in other places of the conductor under mudline. The largest bending moment of the riser focuses on the LSJ, and the moment becomes smaller at the place of the conductor under the ML −5 m. The deformation, stress, and surrounding soil reaction of the casing string change with time in the slowdrift period.
Based on the concept of section division and continuation, a semianalytical approach for analyzing free vibration of the SBOP riserconductor with variable crosssection is proposed, which can actually be applied to any variable crosssection. And this method established for SBOP system natural frequency analysis is reasonable. Results show that the mode shapes have some difference between the coupled and decoupled method. The natural frequencies at diverse modes have little variation with variable TTR. The soil types surrounding the conductor under mudline have very tiny effect on the natural frequency of the riserconductor system.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this article.
Acknowledgments
This work was done at the Drilling Technology Laboratory (DTL) at the Memorial University of Newfoundland, Canada. This work was supported by the National Natural Science Foundation of China (Grant no. 51004119), the Chongqing Research Program of Basic Research and Frontier Technology (Grant no. CSTC2015jcyjA90021), and the Academician Led Special Project of Chongqing Science and Technology Commission (Grant no. cstc2017zdcyyszxX0009).
References
 P. Azancot, E. Magne, and J. Zhang, Surface BOP—Management System Design Guidelines. Society of Petroleum Engineers, 2002.
 “Guidelines for Surface BOP Drilling from Floating MODUs,” in Proceedings of the International Association of Drilling Contractors (IADC), 2015. View at: Google Scholar
 J. R. Kozicz, “Surface BOP—Recent Experience and Future Opportunities,” in Proceedings of Society of Petroleum Engineers, Mumbai, India, 2006. View at: Publisher Site  Google Scholar
 G. Brander, E. Magne, and T. Newman, “Surface BOP technology steps into deeper water with DP vessels,” Oil Gas Journal, vol. 102, no. 17, p. 65, 2004. View at: Google Scholar
 A. Simondin, D. MacPherson, N. Touboul, and G. Ragnes, “A deepwater well construction alternative: surface BOP drilling concept using environmental safe guard,” in Proceedings of the IADC/SPE Drilling Conference, pp. 153–160, Society of Petroleum Engineers, TX, USA, March 2004. View at: Google Scholar
 B. Tarr, T. Taklo, L. A. Olijnik et al., “Surface BOP system operational experience offshore Brazil in 1,900 m of water,” in Proceedings of the SPE/IADC Drilling Conference and Exhibition, Amsterdam, The Netherlands, 2009. View at: Publisher Site  Google Scholar
 L. C. Claassen, S. M. Hendriks, and G. H. T. Zijderveld, “Newbuild compact deepwater drillship designed for surface BOP system,” in Proceedings of the Offshore Technology Conference 2010, OTC 2010, pp. 1391–1400, USA, May 2010. View at: Google Scholar
 American Petroleum Institute, API RP 16Q Recommended Practice for Design, Selecfion, Operafion and Maintenance of Marine Drilling Riser Systems. Washington, DC, USA: American Petroleum Institute, 2001.
 J. C. Botke, “An Analysis of the Dynamics of Marine Risers,” Tech. Rep., Delco Electronics, 1975. View at: Google Scholar
 T. N. Gardner and M. A. Kotch, “Dynamic Analysis of Risers and Caissons by the Element Method,” in Proceedings of the Offshore Technology Conference, Houston, TX, USA. View at: Publisher Site  Google Scholar
 J. M. Shaughnessy, W. T. Daugherty, R. L. Graff, and T. Durkee, “More Ultradeepwater Drilling Problems,” in Proceedings of the SPE/IADC Drilling Conference, Feburary 2007. View at: Publisher Site  Google Scholar
 G. W. King, K. Burton, and T. Hodgson, “A coupled analysis approach to the assessment of marine drilling systems,” SPE Drilling & Completion, vol. 8, no. 2, pp. 131–137, 1993. View at: Publisher Site  Google Scholar
 K. Su, Z. Guan, and Y. Su, “Mechanical Stability Analysis of Subsea Wellhead for Deepwater Drilling,” Oil Drilling Production Technology, vol. 30, no. 6, pp. 1–4, 2008. View at: Google Scholar
 W. Yan, Z.J. Chen, J.G. Deng, H.Y. Zhu, F.C. Deng, and Z.L. Liu, “Numerical method for subsea wellhead stability analysis in deepwater drilling,” Ocean Engineering, vol. 98, pp. 50–56, 2015. View at: Publisher Site  Google Scholar
 C. K. Morooka, R. I. Tsukada, and D. M. Brandt, “Numerical simulations of ocean drilling system behavior with a surface or a subsea BOP in waves and current,” in Proceedings of the 27th International Conference on Offshore Mechanics and Arctic Engineering, OMAE, pp. 469–476, Germany, June 2008. View at: Publisher Site  Google Scholar
 M. W. Dib, J. Lou, L. Zhu, and M. Bassey, “SBOP Drilling Enables Efficient Drilling in Extreme Water Depths,” in Proceedings of the Offshore Technology Conference, Houston, TX, USA, 2009. View at: Publisher Site  Google Scholar
 S. S. Rao, Vibration of Continuous Systems, John Wiley Sons, Inc, Hoboken, NJ, USA, 2007.
 Tsinghua University, Mechanical Vibration. China Mechine Press,1980.
 M. Novak, “Dynamic stiffness and damping of piles,” Canadian Geotechnical Journal, vol. 11, no. 4, pp. 574–598, 1974. View at: Publisher Site  Google Scholar
 G. Gazetas and R. Dobry, “Horizontal response of piles in layered soils,” Journal of Geotechnical Engineering, vol. 110, no. 1, pp. 20–40, 1984. View at: Publisher Site  Google Scholar
 T. Wang, X. Zhang, and W. Zhu, “Vessel Motion Effects on Nonlinear Dynamics of Deepwater Drilling Riser,” Journal of Ship Mechanics, vol. 14, no. 6, pp. 606–617, 2010. View at: Google Scholar
 R. M. Sexton and L. k. Agbezuge, “Random Wave and Vessel Motion Effects on Drilling Riser Dynamics,” in Proceedings of the Offshore Technology Conference, Houston, Texas, 1976. View at: Publisher Site  Google Scholar
 C. Cui, H. Jiang, and Y.H. Li, “Semianalytical method for calculating vibration characteristics of variable crosssection beam,” Journal of Vibration and Shock, vol. 31, no. 14, pp. 85–88, 2012. View at: Google Scholar
Copyright
Copyright © 2018 Kanhua Su 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.