Research Article  Open Access
Piguang Wang, Yifu Chang, Mi Zhao, Junyan Han, "Earthquake and Wave Analysis of Circular Cylinder considering WaterStructureSoil Interaction", Advances in Civil Engineering, vol. 2020, Article ID 4271378, 18 pages, 2020. https://doi.org/10.1155/2020/4271378
Earthquake and Wave Analysis of Circular Cylinder considering WaterStructureSoil Interaction
Abstract
Offshore structures in zones of active seismicity are under a potential threat caused by the combined action of earthquakes and waves. Taking a submerged circular cylinder as the prototype and considering waterstructuresoil interaction, the present study is devoted to the investigation of the combined action of earthquakes and waves. Watercylinder interaction and soilstructure interaction are simulated by added mass and rigid circular massless foundations, respectively. Based on the radiation and diffraction wave theory, the scaled boundary finite element method is utilized to determine the earthquakeinduced and waveinduced pressure on a circular cylinder. Then, a closedform expression for the natural frequencies and mode shapes of the system is derived by using the transfer matrix method, where the transfer matrix is obtained based on EulerāBernoulliās beam differential equation. Furthermore, the dynamic response of the system under the combined action of earthquakes and waves is derived by using the mode superposition method. Finally, the effects of the hydrodynamic force, wave force, and soilstructure interaction on the dynamic response of the submerged cylinder are investigated. The results indicate that the wave forces can substantially increase the dynamic responses of the cylinder and that the influence increases as the stiffness ratio increases and the widthdepth ratio decreases. It is necessary to consider the combined action of earthquakes and waves in the seismic design of offshore structures.
1. Introduction
An increasing number of offshore structures, such as seacrossing bridges, artificial islands, and offshore wind turbines, have been constructed around the world in recent years especially in China. Ocean wave loads are the most important environmental load to consider for the design of offshore structures. However, earthquake loads will dominate the design in zones with high seismic intensities, such as the eastern coast of China and the western coast of the USA. This means that an earthquake can coincide with a wave condition [1]. Therefore, to guarantee the reliability of offshore structures, it is necessary to investigate the dynamic response of offshore structures under the combined action of earthquakes and waves. This paper takes a submerged circular cylinder as the prototype to study the earthquake and wave responses of offshore structures including waterstructuresoil interaction.
The dynamic analysis of offshore structures requires special considerations because of the fluidstructure interaction, which does not occur for onshore structures. The fluidstructure interaction produces an additional hydrodynamic force on the structure when it vibrates in water. Studies show that the additional hydrodynamic force can modify the dynamic properties of the structure and result in structural damage [2ā4]. Earthquakeinduced hydrodynamic forces on a circular cylinder have been studied by many researchers [5] who initially investigated the effects of hydrodynamic forces on the earthquake response of cantilever circular cylinders. The results indicated that the water compressibility is negligible for slender cylinders. The surface waves have no influence on earthquakeinduced hydrodynamic forces except at low loading frequencies and will have little consequence on the seismic response of cylinders surrounded by water. The seismic responses of circular cylinders were further investigated by Williams [6] and Tanaka and Hudspeth [7]. Considering water compressibility, Du et al. [8] proposed a simplified formula for the hydrodynamic force on a circular cylinder in the time domain and Wang et al. [9] proposed an accurate timedomain model to replace the 3D infinite water layer in the watercylinder interaction analysis. The effects of the surrounding water can be entirely equivalent to an added mass when water compressibility is neglected. The simplified methods for evaluating the added hydrodynamic mass for circular, elliptical, roundended, and rectangular cylinders were given by Goyal and Chopra [10]; Li and Yang [11]; Yang and Li [12]; Jiang et al. [13]; and Wang et al. [14ā16].
Wave forces acting on vertical cylinders have also been investigated by many researchers. When the wavelength is substantially longer than the cylinder diameter, the Morison equation [17] can be used to estimate the wave forces on slender cylinders. The Morison equation is made up of two components, including a viscous drag force and an inertia force. However, for large offshore structures, the viscous drag force can be neglected. MacCamy and Fuchs [18] developed the diffraction wave theory to evaluate the wave forces on large circular cylinders immersed in water. Furthermore, the analytical solutions to calculate the wave forces on elliptical cylinders and cylinders with an arbitrary smooth cross section were proposed by Chen and Mei [19]; Williams [20]; and Liu et al. [21]. Li et al. [22] developed a semianalytical solution method to calculate the wave forces on cylinders with arbitrary shapes.
Nevertheless, only a few studies have been carried out to investigate the dynamic response of offshore structures subjected to combined earthquake and wavecurrent action. Penzien et al. [23] studied the dynamic responses of fixedbase offshore towers subjected to random waves and earthquakes. Liu et al. [24] conducted an experiment to study the combined action of earthquakes, waves, and currents on a pile group cablestayed bridge tower foundation. Ding et al. [25] performed underwater shaking table tests to investigate the effect of wavecurrent action on the seismic responses of piers. It should be noted that the soilstructure interaction was not taken into account. Studies have indicated that the soilstructure interaction may have a substantial influence on the dynamic response of offshore structures [26]. Considering the waterstructuresoil interaction, Goyal and Chopra [10,27] and Xu and Spyrakos [28] investigated the waterstructure interaction and soilstructure interaction on the seismic responses of intake towers. The results showed that the seismic response was increased by the waterstructure interaction and decreased by the soilstructure interaction. Yamada et al. [29] studied the dynamic response of offshore towers considering soilstructure interaction subjected to random waves and random earthquake ground motions. It was observed that the effects of wave forces on the seismic response became important with increasing tower period or water depth. This paper presents a simple model to compute the dynamic responses of circular cylinders including waterstructuresoil interaction under the combined action of earthquakes and waves. The effects of waterstructure interaction and soilstructure interaction are systematically analyzed.
2. Mathematical Formulation
The waterstructuresoil interaction problem in threedimensional space is shown in Figure 1. A circular cylinder with radius a and height H is situated in water with a mean depth h. The mass of the superstructure is M_{S}. The soilstructure interaction is considered to be a rigid circular massless foundation with radius R_{0} on deep soil strata, where the translational and rotational stiffness of the foundation are denoted by spring constants k_{h} and k_{r}, respectively. The foundation has an acceleration time history of the earthquake motion along the xdirection. The Cartesian coordinate system (x, y, z) is defined with an origin located on the seabed level with the zaxis directed vertically upwards. The center of the cylinder is taken as the origin of the polar coordinate system (r, Īø), where Īø is measured counterclockwise from the positive xaxis. The fluid is assumed to be incompressible and inviscid. The system is initially at rest.
The EulerBernoulli beam theory is adopted for the dynamic analysis of the cylinder. The differential equations of the cylinder subjected to wave and earthquake actions can be expressed as follows:where u and are the transverse deflection and acceleration along the zdirection, is the bending stiffness of the pile, is the density of the cylinder, is the crosssectional area of the cylinder, is the earthquakeinduced hydrodynamic force on the cylinder, and is the wave force on the cylinder.
In the polar coordinate system, the earthquakeinduced hydrodynamic force can be expressed as follows [16]:where and the hydrodynamic pressure is governed by a twodimensional Helmholtz equation with the boundary conditions on the surface of the cylinder as follows:in which is the water density and is the Laplace operator.
Assuming that the wavelength is substantially longer than the cylinder diameter, the wave force can be expressed as [21]where , is the gravitational acceleration, is the wave height, k is the wavenumber, is the wave frequency, is the incident wave pressure, and is the scattered wave pressure. In the polar coordinate system, is governed by a twodimensional Helmholtz equation with the boundary conditions on the surface of the cylinder as follows:
In addition, the radiation pressure waves and propagate away from the structure and decay with the traveling distance. The range of the applications of the present wave model is 2a/L_{0}ā>ā0.2, where L_{0} represents the wavelength.
3. EarthquakeInduced and WaveInduced Pressures
The scaled boundary finite element method (SBFEM) developed by Wolf [30] was used to determine the earthquakeinduced and waveinduced pressures on a circular cylinder. The formulation of the scaled boundary finite element (SBFE) governing equations of the 2D pressure P_{j} or P_{s} is derived by converting (3) or (6) into the scaled boundary coordinates with the scaling center at the center of the cylinder in the present study, as shown in Figure 2. The circumferential coordinate is the anticlockwise direction along the defining curve S, which is closed in this case. The normalized radial coordinate defined as 1 at curve S is a scaling factor with for the infinite water domain. Therefore, the polar coordinates are transformed to the scaled boundary coordinates and with the following scaling equations:
Following the concept of an isoparametric element, the pressure Pā=āP_{j} or Pā=āP_{s} at a point is calculated bywhere is the shape function, the vector is the hydrodynamic pressure along the radial lines and is analytical with respect to .
By performing a scaled boundary transformation, the operator can be written aswhere is related to the determinant of the Jacobian matrix at the boundary:
The infinitesimal volume for any radial coordinate is calculated as
The finite element method requires the weighted residuals of the governing equation in (3) or (6) to be zero. Applying the Galerkin approach, the weighting function can be chosen as the same shape function as in (9):
3.1. Solution of the EarthquakeInduced Pressure
By applying the weighted residual method and performing integration by parts, the solution of the differential in (3) with the boundary condition in (4) is equivalent to the following weak form:where .
Substituting (9), (10), (13), and (14) into (15), after some manipulations, the scaled boundary finite element equation and the boundary condition equation for the earthquakeinduced 2D pressure can be unified aswhere , , , , and are expressed asin which , , is the identity matrix of rank n, and n denotes the number of nodes in curve S.
Using and , premultiplying both sides of (17) by and simplifying, we getwhere . Equation (19) is the matrix form of the modified Bessel differential equation. Considering the radiation condition, it is logical to select as a base solution for (19). Therefore, the solution for can be expressed in the form of a series aswhere are coefficients, are vectors of rank n, and are the modified Bessel functions of the second kind.
Substituting (20) into (19) and using the following properties of the modified Bessel function of the second kind givewhere the prime and the double prime represent the first and second derivatives with respect to argument , respectively, and we obtain
Usually, , so that
Let the variable be the eigenvalues of ; then ; and are the eigenvectors of .
Substituting (20) into (17), we obtain
Using (20) and (24), the pressure can be obtained asHere, the dynamicstiffness matrix is expressed aswhere ādiagā represents a diagonal matrix with the elements in square brackets on the main diagonal.
3.2. Solution of the WaveInduced Pressure
By applying the weighted residual method and performing integration by parts, the solution of the differential equation in (6) with the boundary condition in (7) is equivalent to the following weak form:where .
Substituting (9), (10), (13), and (14) into (27), after some manipulations, the scaled boundary finite element equation and the boundary condition equation for the waveinduced 2D pressure can be unified aswhere . Using and , premultiplying both sides of (28) by and simplifying, we obtainwhere .
Equation (30) is the matrix form of Besselās differential equation. Considering the radiation condition, it is logical to select as a base solution of (30). Therefore, the solution for can be expressed in the form of a series aswhere are the Hankel functions of the first kind.
Substituting (31) into (30) and using the properties of the Hankel function give
We obtain
Substituting (33) into (29), we obtain
Using (34) and (35), the pressure can be obtained aswhere the dynamicstiffness matrix is expressed as
4. Solution of the Equations of Motion
The simplified model of the fluidcylindersoil system is shown in Figure 3. The mode superposition method is utilized to solve the equations of motion in the present study. The equations of motion of the free vibration are first solved, and the results are further used to solve the dynamic responses of the cylinder under the combined action of waves and earthquakes.
The free vibration equation of the cylinder vibrating in water can be written aswhere denotes the added mass on the cylinder due to the watercylinder interaction.
Assuming a harmonic solution and applying the method of the separation of variables, the transverse deflection u can be expressed as
Substituting (39) into (38), we obtain
The four boundary conditions associated with the equations of motion in (40) can be expressed as follows:(1)Bending moment at zā=ā0:(2)Shear force at zā=ā0:(3)Bending moment at zā=ā:(4)Shear force at zā=ā:
Generally, the differential equations of motion in (40) do not possess an exact solution. In addition, the added mass and stiffness of the foundation change along the cylinder and pile. Therefore, the transfer matrix method (TMM) is used to solve the previous equations of motion. During analysis by the TMM, the cylinder and pile are discretized into L continuous beam elements as shown in Figure 3, and the equation of motion for each element is derived. The added mass and stiffness of the foundation are assumed to be constant along each element. The analytical solution of the motion for each element is derived in the following section.
4.1. Free Vibration
The differential equation of motion for the ith element can be written aswhere , is the equivalent added mass, , and is the angular frequency. The general solution form of (42) can be expressed as follows:where .
According to the continuity conditions of the displacement, slope, shear force, and bending moment between the ith and (iā+ā1)th elements, we obtainin which derived in Appendix is the transfer matrix; the undetermined coefficients and are expressed by
By using the TMM, the relationship between and can be given as follows:
According to the boundary conditions in (41a) and (41b), we obtain
According to boundary conditions in (41c) and (41d), we obtain
Substituting (47) into equations ((50a), (50b)), we get
Furthermore, using (49a) into (51), we obtain
The constant vector in (53) cannot be zero. Therefore, the equation governing the natural frequencies of the cylinder is written as
The natural frequencies of the system denoted by can be obtained by solving (54), and the corresponding mode shape is denoted by . It should be noted that (54) is a transcendental equation that needs to be solved numerically. In this study, (54) is solved using MATLAB software.
4.2. Dynamic Responses in the Time Domain
By using the mode superposition method, the deformation of the cylinder can be written as follows:where are the generalized coordinates and the damping ratio is denoted by . Substituting (55) into (1), considering the damping and using the orthogonality of the mode shapes, the generalized coordinates in the equation of motion of the cylinder subjected to earthquake and wave actions are expressed as [5]
Equation (56) is a dynamic equation of a singledegreeoffreedom system. It can be solved by the standard implicit time integration algorithm such as the Newmark or Wilson method. In addition to the implicit algorithm, it can also be solved by the secondorder accurate explicit algorithm proposed by the authors [31].
4.3. Dynamic Responses in the Frequency Domain
The deformation of the cylinder under earthquake action can be written as follows:where is the frequency generalized coordinates, which can be calculated bywhere is the ground acceleration in the frequency domain.
5. Results and Discussion
The material parameters of the flexible cylinder, including Youngās modulus, density, and damping ratio, are 30000āMPa, 2500ākg/m^{3}, and 0.05, respectively. The water depth is equal to the height of the cylinder. The peak acceleration of the ground motion (a_{max}) is equal to 0.1āg. The wave height and wave period are ā=ā3ām and Tā=ā8ās, respectively. The translational and rotational spring constants k_{h} and k_{r} for a rigid circular massless foundation with radius R_{0} on deep soil strata can be evaluated by the following frequencyindependent expressions [32]:where and ā=ā0.3 are the shear modulus and Poissonās ratio of the soil, respectively.
The effects of the fluidstructure interaction (FSI), soilstructure interaction (SSI), and fluidstructuresoil interaction (FSSI) on the watercylindersoil system can be better understood by studying the system dynamic properties in terms of dimensionless parameters. Therefore, four dimensionless parameters, including the widthdepth ratio (), mass ratio (), heightradius ratio (), and stiffness ratio (), are introduced as
5.1. Validation
The solutions for the earthquakeinduced pressure and waveinduced pressure are first verified. Figure 4 shows the comparison of the hydrodynamic pressure P_{j} on a circular cylinder between the present SBFEM and the finite element method (FEM) developed by Wang et al. [33]. Figure 5 shows the comparison of the total wave forces on a circular cylinder computed by the present SBFEM and FEM developed by Wang et al. [33], where . It can be seen that the proposed SBFEM agrees well with the FEM.
The solution of the free vibration is further validated by assuming the foundation rigid. The first natural frequencies of the circular cylinders vibrating in water are compared with those in a previous study [15], where Eā=ā29.4āGPa, Ļ_{c}ā=ā2450ākg/m^{3}, and Hā=āhā=ā20ām. The natural frequencies and the difference between the two studies are tabulated in Table 1. It is evident that good agreement is observed.

The solution of the dynamic response is finally verified by assuming the foundation rigid. Figure 6 shows the comparison of the displacement amplitude at the top of a circular cylinder excited by the harmonic ground motion between the proposed solution and the experimental values [34]. It can be seen that good agreement is observed. Figure 7 shows the comparison of the displacement time history at the top of a circular cylinder excited by wave forces between the proposed solution and the software ABAQUS, where Eā=ā29.4āGPa, Ļ_{c}ā=ā2450ākg/m^{3}, aā=ā2ām, Hā=āhā=ā20ām, ā=ā2ām, and . It can also be seen that good agreement is observed.
5.2. Free Vibration
The fundamental frequencies of the cylinder in air with a rigid foundation, considering the SSI, FSI, and FSSI, are denoted by , , , and .
The effects of the SSI on the fundamental frequency of the cylinder are first investigated. Figure 8 shows the ratio with respect to for different values of and H with ā=ā0 and ā=ā0.5. It can be seen that the ratio decreases as increases, which means that the fundamental frequency of the cylinder decreases when the stiffness of the foundation decreases. It can also be seen that the ratio is independent of and H. Figure 9 shows the ratio with respect to for different values of with ā=ā0.5. It can be seen that the ratio increases as increases which means that the effects of the SSI on the fundamental frequency of the cylinder decrease as increases. Figure 10 shows the ratio with respect to for different values of with ā=ā0. It can be seen that the ratio decreases marginally with increasing .
The effects of the FSI on the fundamental frequency of the cylinder are investigated next. Figure 11 shows the ratio with respect to for different values of and H. It can be seen that the ratio increases as and increase, which means that the effects of the FSI on the fundamental frequency of the cylinder decrease as and increase. It can also be seen that the ratio is independent of H. Figure 12 shows the ratio with respect to for different values of with ā=ā0 and ā=ā0.2. It can be seen that the ratio decreases marginally with the increase in and decrease in .
The effects of the SSI and FSI on the vibration shapes of the cylinder are thirdly investigated, where Hā=ā8ām, ā=ā0.1, ā=ā0, and ā=ā0.5. Figure 13 shows the vibration mode shapes of the cylinder in a vacuum for different values of G with Hā=ā8ām, ā=ā0.1, ā=ā0, and ā=ā0.5. The SSI has a substantial effect on the vibration mode shapes of the cylinder when the soil stiffness is weak, especially for the highorder mode. Figure 14 shows the vibration mode shapes of the cylinder in a vacuum and in water with Hā=ā8ām, Gā=ā50āMPa, ā=ā0.1, ā=ā0, and ā=ā0.5. It can be seen that the FSI has little effect on the fundamental vibration modes, while it has a slight effect on the higher modes.
(a)
(b)
(c)
(a)
(b)
(c)
5.3. Dynamic Responses in the Time Domain
The dominant frequency of the seismic waves is substantially different for the different project sites. The farfield suite of 22 records introduced in the FEMAP695 [35] selected from the NGAWest2 ground motion database of the Pacific Earthquake Engineering Research Center is also considered in the present study. The details of the three ground motions are given in Table 1. The dynamic responses of the pier are investigated when the peak acceleration (a_{max}) is equal to 0.1āg, ā=ā0, ā=ā3ām, and the wave period ā=ā8ās.
The effects of the SSI on the seismic responses of the cylinder are first investigated. The peak displacements on the top of the cylinder with rigid foundation and considering the SSI are denoted by u_{0max} and u_{1max}, respectively. The peak bending moment on the base of the cylinder with a rigid foundation and considering the SSI are denoted by bm_{0max} and bm_{1max}, respectively. Two dimensionless parameters are introduced as Ru_{1}ā=āu_{1max}/u_{0max} and Rm_{1}ā=ābm_{1max}/bm_{0max}. Figures 15 and 16 show the parameters Ru_{1} and Rm_{1} versus for different values of H with ā=ā0.2 and ā=ā0.5. Figures 17 and 18 show the parameters Ru_{1} and Rm_{1} versus for different values of H with ā=ā5 and ā=ā0.5. The SSI increases the displacement of the cylinder and this influence increases as and increase and H decreases, which means that the SSI has greater effect on the dynamic response of the squatty cylinder. The bending moment of the cylinder may decrease due to the SSI when H is large and is small, while it may increase due to the SSI when H is small and is large.
(a)
(b)
(c)
(a)
(b)
(c)
(a)
(b)
(c)
(a)
(b)
(c)
The effects of the FSI on the seismic responses of the cylinder are secondly investigated next, where ā=ā0.5. The peak displacement on the top of the cylinder and the peak bending moment on the base of the cylinder considering the FSSI are denoted by u_{2max} and bm_{2max}, respectively. Two dimensionless parameters are introduced as Ru_{2}ā=āu_{2max}/u_{1max} and Rm_{2}ā=ābm_{2max}/bm_{1max}. Figures 19 and 20 show the parameters Ru_{2} and Rm_{2} versus for different values of H and . It can be seen that the FSI increases the displacement and bending moment of the cylinder. However, the parameters Ru_{2} and Rm_{2} have no obvious tendency as increases.
(a)
(b)
(c)
(a)
(b)
(c)
The effects of wave forces on the dynamic responses of the cylinder are investigated thirdly, where ā=ā0.5. The peak displacements on the top of the cylinder under the combined action of earthquakes and waves are denoted by u_{3max}. A dimensionless parameter is introduced as Ru_{3}ā=āu_{3max}/u_{2max}. Figure 21 shows the history of the relative displacement on top of the cylinder under the combined action of earthquakes and waves under the HECTOR earthquake (ID number 4 in Table 2), where ā=ā0.2, ā=ā0.5, and Hā=ā40ām. Figure 22 shows the parameters Ru_{3} versus for different values of H with ā=ā0.2. Figure 23 shows the parameters Ru_{3} versus for different values of H with ā=ā5. The wave forces substantially increase the dynamic responses of the cylinder and this influence increases as increases and decreases. It can also be seen from Figure 17 that the dynamic responses of the cylinder under the combined action of earthquakes and waves display obvious periodicity.
(a)
(b)

(a)
(b)
(c)
(a)
(b)
(c)
5.4. Dynamic Responses in the Frequency Domain
The effects of the FSI on the frequency responses of the cylinder are investigated in this section, where Hā=ā40ām and ā=ā0.
Figure 24 shows the displacement amplitude of the cylinder on the top with respect to the earthquake frequency in the case of the SSI and FSSI, where ā=ā0.2 and ā=ā0.5. The frequency responses of the cylinder increase due to the FSI. Figure 25 shows the ratio with respect to for different values of with ā=ā0.5. Figure 26 shows the ratio with respect to for different values of with ā=ā0.2. It can be seen that the ratio decreases marginally as increases. In general, the dynamic responses of the cylinder increase due to the FSI by approximately 30%ā¼40%.
(a)
(b)
6. Conclusions
The effects of fluidstructuresoil interaction on the dynamic response of offshore cylinders are investigated by idealizing the cylinder as a uniform beam, replacing the watercylinder interaction with added mass, and replacing the soilstructure interaction with translational and rotational springs. The earthquakeinduced and waveinduced pressure on a circular cylinder are determined by the scaled boundary finite element method. The closedform expression for the natural frequencies and mode shapes of the system is derived by using the transfer matrix method and the dynamic response of the system under the combined action of earthquakes and waves further obtained and derived by the time integration algorithm. From all the results presented in the present study, the following conclusion can be made:(1)The SSI decreases the fundamental frequency of the cylinder and the influence increases as the stiffness ratio and heightradius ratio increase and mass ratio decreases. The FSI decreases the fundamental frequency of the cylinder and the influence increases as the widthdepth ratio and mass ratio decrease. The SSI has a substantial effect on the vibration mode shapes of the cylinder when the soil stiffness is weak and the FSI has a slight effect on the higher modes.(2)The SSI increases the displacement of the cylinder and the influence increases as the stiffness ratio and widthdepth ratio increase, and the cylinder height decreases. The SSI can decrease the bending moment of the cylinder for slender cylinders with large heights, while it can increase the bending moment of the cylinder for squatty cylinders with small heights.(3)The FSI increases the displacement and bending moment of the cylinder and the influence decreases marginally with increasing widthdepth ratio. It should be noted that the effects of the FSI on the dynamic responses of the cylinder have no obvious tendency as the stiffness ratio increases.(4)The wave forces can substantially increase the dynamic responses of the cylinder, and the influence increases as the stiffness ratio increases and the widthdepth ratio decreases. The wave forces may be the dominant loading in the seismic design of offshore structures for soft soil. Consequently, it is necessary to consider the combined action of earthquakes and waves in the dynamic analysis of offshore structures.
Generally, during an earthquake, the wave condition and the soil behaviour will be significantly affected by the ground movement of the earthquake. However, the present model is a simplified model without considering the interaction between the earthquake, wave, and soil. The linear regular wave is only used to investigate the response trends of the circular cylinder subjected to earthquake and wave loadings in the present study. However, a wave subjected to an earthquake cannot be modeled as a linear regular wave in a real scenario. In practical engineering applications, the wave spectrum, such as the JONSWAP spectrum [36] and PiersonMoskowitz spectrum [37], can be used to define a stationary random sea state caused by windgenerated surface waves. It should be noted that the linear wave model is not sufficient to represent the wave condition during an earthquake, which is usually highly nonlinear. Therefore, the range of applications of the present linear wave model is h/L_{0}ā>ā0.2 and //hāā¤ā0.2 [38]. In future work, we will improve the present model by considering the nonlinear wave force on the cylinder.
Appendix
According to the continuity conditions of the displacement, slope, shear force, and bending moment between the ith and (iā+ā1)th elements, we obtainwhere the matrices and are expressed as
Equation (A.1) can be further rewritten aswhere .
Data Availability
No data were used to support this study.
Conflicts of Interest
The authors declare that there are no conflicts of interest regarding the publication of this paper.
Acknowledgments
This work was supported by the National Natural Science Foundation of China (51678015 and 51708010) and Ministry of Education Innovation Team of China (IRT_17R03).
References
 X. Y. Zheng, H. Li, W. Rong, and W. Li, āJoint earthquake and wave action on the monopile wind turbine foundation: an experimental study,ā Marine Structures, vol. 44, pp. 125ā141, 2015. View at: Publisher Site  Google Scholar
 R. P. S. Han and H. Xu, āA simple and accurate added mass model for hydrodynamic fluidstructure interaction analysis,ā Journal of the Franklin Institute, vol. 333, no. 6, pp. 929ā945, 1996. View at: Publisher Site  Google Scholar
 B. Song, F. Zheng, and Y. Li, āStudy on a simplified calculation method for hydrodynamic pressure to slender structures under earthquakes,ā Journal of Earthquake Engineering, vol. 17, no. 5, pp. 720ā735, 2013. View at: Publisher Site  Google Scholar
 K. Wei, W. Yuan, and N. Bouaanani, āExperimental and numerical assessment of the threedimensional modal dynamic response of bridge pile foundations submerged in water,ā Journal of Bridge Engineering, vol. 18, no. 10, pp. 1032ā1041, 2013. View at: Publisher Site  Google Scholar
 C.Y. Liaw and A. K. Chopra, āDynamics of towers surrounded by water,ā Earthquake Engineering & Structural Dynamics, vol. 3, no. 1, pp. 33ā49, 1974. View at: Publisher Site  Google Scholar
 A. N. Williams, āEarthquake response of submerged circular cylinder,ā Ocean Engineering, vol. 13, no. 6, pp. 569ā585, 1986. View at: Publisher Site  Google Scholar
 Y. Tanaka and R. T. Hudspeth, āRestoring forces on vertical circular cylinders forced by earthquakes,ā Earthquake Engineering & Structural Dynamics, vol. 16, no. 1, pp. 99ā119, 1988. View at: Publisher Site  Google Scholar
 X. Du, P. Wang, and M. Zhao, āSimplified formula of hydrodynamic pressure on circular bridge piers in the time domain,ā Ocean Engineering, vol. 85, pp. 44ā53, 2014. View at: Publisher Site  Google Scholar
 P. Wang, M. Zhao, H. Li, and X. Du, āAn accurate and efficient timedomain model for simulating watercylinder dynamic interaction during earthquakes,ā Engineering Structures, vol. 166, pp. 263ā273, 2018. View at: Publisher Site  Google Scholar
 A. Goyal and A. K. Chopra, āSimplified evaluation of added hydrodynamic mass for intake towers,ā Journal of Engineering Mechanics, vol. 115, no. 7, pp. 1393ā1412, 1989. View at: Publisher Site  Google Scholar
 Q. Li and W. Yang, āAn improved method of hydrodynamic pressure calculation for circular hollow piers in deep water under earthquake,ā Ocean Engineering, vol. 72, pp. 241ā256, 2013. View at: Publisher Site  Google Scholar
 W. Yang and Q. Li, āA new added mass method for fluidstructure interaction analysis of deepwater bridge,ā KSCE Journal of Civil Engineering, vol. 17, no. 6, pp. 1413ā1424, 2013. View at: Publisher Site  Google Scholar
 H. Jiang, B. Wang, X. Bai, C. Zeng, and H. Zhang, āSimplified expression of hydrodynamic pressure on deep water cylindrical bridge piers during earthquakes,ā Journal of Bridge Engineering, vol. 22, no. 6, Article ID 04017014, 2017. View at: Publisher Site  Google Scholar
 P. Wang, M. Zhao, and X. Du, āAnalytical solution and simplified formula for earthquake induced hydrodynamic pressure on elliptical hollow cylinders in water,ā Ocean Engineering, vol. 148, pp. 149ā160, 2018. View at: Publisher Site  Google Scholar
 P. Wang, M. Zhao, and X. Du, āA simple added mass model for simulating elliptical cylinder vibrating in water under earthquake action,ā Ocean Engineering, vol. 179, pp. 351ā360, 2019. View at: Publisher Site  Google Scholar
 P. Wang, M. Zhao, and X. Du, āSimplified formula for earthquakeinduced hydrodynamic pressure on roundended and rectangular cylinders surrounded by water,ā Journal of Engineering MechanicsASCE, vol. 145, no. 2, Article ID 04018137, 2019b. View at: Publisher Site  Google Scholar
 J. R. Morison, J. W. Johnson, S. A. Schaaf, and S. A. Schaaf, āThe force exerted by surface waves on piles,ā Journal of Petroleum Technology, vol. 2, no. 5, pp. 149ā154, 1950. View at: Publisher Site  Google Scholar
 R. C. MacCamy and R. A. Fuchs, Wave forces on piles: a diffraction theory. Beach Erosion Board Technical Memorandum U.S. Army Corps of Engineers, Report 69, United States Army Corps of Engineers, Washington, DC, USA, 1954.
 H. S. Chen and C. C. Mei, āWave forces on a stationary platform of elliptical shape,ā Journal of Ship Research, vol. 17, no. 2, pp. 61ā71, 1973. View at: Google Scholar
 A. N. Williams, āWave forces on an elliptic cylinder,ā Journal of Waterway, Port, Coastal, and Ocean Engineering, vol. 111, no. 2, pp. 433ā449, 1985. View at: Publisher Site  Google Scholar
 J. Liu, A. Guo, and H. Li, āAnalytical solution for the linear wave diffraction by a uniform vertical cylinder with an arbitrary smooth crosssection,ā Ocean Engineering, vol. 126, pp. 163ā175, 2016. View at: Publisher Site  Google Scholar
 B. Li, L. Cheng, A. J. Deeks, and M. Zhao, āA semianalytical solution method for twodimensional Helmholtz equation,ā Applied Ocean Research, vol. 28, no. 3, pp. 193ā207, 2006. View at: Publisher Site  Google Scholar
 J. Penzien, M. K. Kaul, and B. Berge, āStochastic response of offshore towers to random sea waves and strong motion earthquakes,ā Computers and Structures, vol. 2, no. 56, pp. 733ā756, 1972. View at: Publisher Site  Google Scholar
 C. Liu, S. Zhang, and E. Hao, āJoint earthquake, wave and current action on the pile group cablestayed bridge tower foundation: an experimental study,ā Applied Ocean Research, vol. 63, pp. 157ā169, 2017. View at: Publisher Site  Google Scholar
 Y. Ding, R. Ma, Y.D. Shi, and Z.X. Li, āUnderwater shaking table tests on bridge pier under combined earthquake and wavecurrent action,ā Marine Structures, vol. 58, pp. 301ā320, 2018. View at: Publisher Site  Google Scholar
 B.F. Chen, āDynamic responses of coastal structures during earthquakes including sedimentseastructure interaction,ā Soil Dynamics and Earthquake Engineering, vol. 20, no. 5ā8, pp. 445ā467, 2000. View at: Publisher Site  Google Scholar
 A. Goyal and A. K. Chopra, āHydrodynamic and foundation interaction effects in dynamics of intake towers: earthquake responses,ā Journal of Structural Engineering, vol. 115, no. 6, pp. 1386ā1395, 1989. View at: Publisher Site  Google Scholar
 C. Xu and C. C. Spyrakos, āSeismic analysis of towers including foundation uplift,ā Engineering Structures, vol. 18, no. 4, pp. 271ā278, 1996. View at: Publisher Site  Google Scholar
 Y. Yamada, H. Iemura, K. Kawano, and K. Venkataramana, āSeismic response of offshore structures in random seas,ā Earthquake Engineering & Structural Dynamics, vol. 18, no. 7, pp. 965ā981, 1989. View at: Publisher Site  Google Scholar
 J. P. Wolf, The Scaled Boundary Finite Element Method, John Wiley & Sons, Chichester, UK, 2003.
 M. Zhao, H. Li, S. Cao, and X. Du, āAn explicit time integration algorithm for linear and nonlinear finite element analyses of dynamic and wave problems,ā Engineering Computations, vol. 36, no. 1, pp. 161ā177, 2019. View at: Publisher Site  Google Scholar
 G. Gazetas, āAnalysis of machine foundations: state of the art,ā Soil Dynamic and Earthquake Engineering, vol. 2, pp. 1ā42, 1983. View at: Publisher Site  Google Scholar
 P. Wang, M. Zhao, X. Du, and X. Cheng, āA finite element solution of earthquakeinduced hydrodynamic forces and wave forces on multiple circular cylinders,ā Ocean Engineering, vol. 189, Article ID 106336, 2019. View at: Publisher Site  Google Scholar
 H. Goto and K. Toki, āVibrational characteristics and a seismic design of submerged bridge piers,ā Proceeding of 3rd World Conference in Earthquake Engineering, pp. 107ā122, 1965. View at: Google Scholar
 ATC/FEMA, āApplied technology council/federal emergency management agancy (ATC/FEMA),ā in FEMA P695: Quantification of Building Seismic Performance Factors, ATC/FEMA, Washington, DC, USA, 2009. View at: Google Scholar
 K. Hasselmann, T. Barnett, E. Bouws, H. Carlson, D. Cartwright, and K. Enke, Measurements of WindWave Growth and Swell Decay during the Joint North Sea Wave Project (JONSWAP), Deutches Hydrographisches Institut, Hamburg, Germany, 1973.
 K.Y. Oh, J.Y. Kim, and J.S. Lee, āPreliminary evaluation of monopile foundation dimensions for an offshore wind turbine by analyzing hydrodynamic load in the frequency domain,ā Renewable Energy, vol. 54, pp. 211ā218, 2013. View at: Publisher Site  Google Scholar
 Y. Zhu, āAnalyses of the range of validity for several wave theories,ā Hai An Gong Cheng, vol. 2, no. 2, pp. 11ā27, 1983, in Chinese. View at: Google Scholar
Copyright
Copyright © 2020 Piguang 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.