A time-domain method is applied to simulate nonlinear wave diffraction around a surface piercing 3-D arbitrary body. The method involves the application of Taylor series expansions and the use of perturbation procedur...A time-domain method is applied to simulate nonlinear wave diffraction around a surface piercing 3-D arbitrary body. The method involves the application of Taylor series expansions and the use of perturbation procedure to establish the corresponding boundary value problems with respect to a time-independent fluid domain. A boundary element method based on B-spline expansion is used to calculate the wave field at each time step, and the free surface boundary condition is satisfied to the second order of wave steepness by a numerical integration in time. An artificial damping layer is adopted on the free surface for the removal of wave reflection from the outer boundary. As an illustration, the method is used to compute the second-order wave forces and run-up on a surface-piercing circular cylinder. The present method is found to be accurate, computationally efficient, and numerically stable.展开更多
The time domain responses of the tunnel element under wave actions during its immersion are investigated based on the linear wave diffraction theory. The integral equation is derived by using the time-domain Green fun...The time domain responses of the tunnel element under wave actions during its immersion are investigated based on the linear wave diffraction theory. The integral equation is derived by using the time-domain Green function that satisfies the free water surface condition in the finite water depth, and is solved by the boundary element method. The motion equations of the tunnel element are solved by the fourth order Runge-Kutta method. A comparison between the computed and measured results reveals that the numerical model can effectively simulate the motion responses of the tunnel element and the cable tensions when the motions of the tunnel element are within some limit. Taking the tunnel element of 100 m in length, 15 m in width and 10 m in height as an example, the computational results of the motion responses of the tunnel element and the cable tensions in different immersing depths are obtained under different incident wave conditions.展开更多
A 3-D iterative Rankine Boundary Element Method （BEM） for seakeeping problem in time domain is developed in the framework of linear potential theory. Waves generated by both submerged and surface-piercing bodies mov...A 3-D iterative Rankine Boundary Element Method （BEM） for seakeeping problem in time domain is developed in the framework of linear potential theory. Waves generated by both submerged and surface-piercing bodies moving at a constant forward speed in otherwise calm water, and the resultant steady wave pattern, wave profile and resistance are computed to validate this newly-developed code. A rectangular computational domain moving with the same forward speed as the body is introduced, in which an artificial damping beach is installed at an outer portion of the free surface except the downstream side for satisfying the radiation condition. The velocity potential on the ship hull and the normal velocity on the free surface are obtained directly by solving the boundary integral equation, with the Rankine source used as the kernel function. An iterative time-marching scheme is employed for updating both kinematic and dynamic free surface boundary conditions to stabilize the calculation. Extensive results including the wave patterns, wave profiles and wave resistances for a submerged spheroid and a Wigley hull with forward speed are presented to validate the efficiency of the proposed 3-D time-domain higher-order approach. Finally, the sensitivity of ship-generated waves to the water depth is investigated. Computed results show satisfactory agreement with the corresponding experimental data and other numerical solutions.展开更多
Environmental effects have an important influence on Offshore Wind Turbine (OWT) power generation efficiency and the structural stability of such turbines. In this study, we use an in-house Boundary Element (BEM)-panM...Environmental effects have an important influence on Offshore Wind Turbine (OWT) power generation efficiency and the structural stability of such turbines. In this study, we use an in-house Boundary Element (BEM)-panMARE code-to simulate the unsteady flow behavior of a full OWT with various combinations of aerodynamic and hydrodynamic loads in the time domain. This code is implemented to simulate potential flows for different applications and is based on a three-dimensional first-order panel method. Three different OWT configurations consisting of a generic 5 MW NREL rotor with three different types of foundations (Monopile, Tripod, and Jacket) are investigated. These three configurations are analyzed using the RANSE solver which is carried out using ANSYS CFX for validating the corresponding results. The simulations are performed under the same environmental atmospheric wind shear and rotor angular velocity, and the wave properties are wave height of 4 m and wave period of 7.16 s. In the present work, wave environmental effects were investigated firstly for the two solvers, and good agreement is achieved. Moreover, pressure distribution in each OWT case is presented, including detailed information about local flow fields. The time history of the forces at inflow direction and its moments around the mudline at each OWT part are presented in a dimensionless form with respect to the mean value of the last three loads and the moment amplitudes obtained from the BEM code, where the contribution of rotor force is lower in the tripod case and higher in the jacket case and the calculated hydrodynamic load that effect on jacket foundation type is lower than other two cases.展开更多
A seismic free field input formulation of the coupling procedure of the finite elelnent(FE)and the scaled boundary finite-element(SBFE)is proposed to perform the unbounded soil-structure interaction analysis in time d...A seismic free field input formulation of the coupling procedure of the finite elelnent(FE)and the scaled boundary finite-element(SBFE)is proposed to perform the unbounded soil-structure interaction analysis in time domain. Based on the substructure technique,seismic excitation of the soil-structure system is represented by the free-field motion of an elastic half-space.To reduce the computational effort,the acceleration unit-impulse response function of the unbounded soil is decomposed into two functions:linear and residual.The latter converges to zero and can be truncated as required. With the prescribed tolerance parameter,the balance between accuracy and effMency of the procedure can be controlled. The validity of the model is verified by the scattering analysis of a hemi-spherical canyon subjected to plane harmonic P,SV and SH wave incidence.Numerical results show that the new procedure is very efficient for seismic problems within a nor- real range of frequency.The coupling procedure presented herein can be applied to linear and nonlinear earthquake re- sponse analysis of practical structures which are built on unbounded soil.展开更多
Developing serpent-type wave generators to generate solitary waves in a 3D-basin was investigated in this study. Based on the Lagrangian description with time-marching procedures and finite differences of the time der...Developing serpent-type wave generators to generate solitary waves in a 3D-basin was investigated in this study. Based on the Lagrangian description with time-marching procedures and finite differences of the time derivative, a 3D multiple directional wave basin with multidirectional piston wave generators was developed to simulate ocean waves by using BEM with quadrilateral elements, and to simulate wave-caused problems with fully nonlinear water surface conditions. The simulations of perpendicular solitary waves were conducted in the first instance to verify this scheme. Furthermore, the comparison of the waveform variations confirms that the estimation of 3D solitary waves is a feasible scheme.展开更多
Based on the modified scale boundary finite element method and continued fraction solution,a high-order doubly asymptotic transmitting boundary(DATB)is derived and extended to the simulation of vector wave propagation...Based on the modified scale boundary finite element method and continued fraction solution,a high-order doubly asymptotic transmitting boundary(DATB)is derived and extended to the simulation of vector wave propagation in complex layered soils.The high-order DATB converges rapidly to the exact solution throughout the entire frequency range and its formulation is local in the time domain,possessing high accuracy and good efficiency.Combining with finite element method,a coupled model is constructed for time-domain analysis of underground station-layered soil interaction.The coupled model is divided into the near and far field by the truncated boundary,of which the near field is modelled by FEM while the far field is modelled by the high-order DATB.The coupled model is implemented in an open source finite element software,OpenSees,in which the DATB is employed as a super element.Numerical examples demonstrate that results of the coupled model are stable,accurate and efficient compared with those of the extended mesh model and the viscous-spring boundary model.Besides,it has also shown the fitness for long-time seismic response analysis of underground station-layered soil interaction.Therefore,it is believed that the coupled model could provide a new approach for seismic analysis of underground station-layered soil interaction and could be further developed for engineering.展开更多
The unsteady flow around a ship advancing at a constant speed was studied in time domain by using a higher-order boundary element method. The Rankine source was employed. The ship hull and free surface were discretize...The unsteady flow around a ship advancing at a constant speed was studied in time domain by using a higher-order boundary element method. The Rankine source was employed. The ship hull and free surface were discretized into 8-node isoparametric elements. A modified Green formula was used to overcome the numerical difficulty arising from singular or quasi-singular integral equation. In the case of τ>0.25, a rigid wall condition was applied to the upstream and a numerical beach was applied to the downstream to absorb the outgoing waves. The numerical results are satisfactory compared with the experimental results.展开更多
To study wave-current actions on 3-D bodies a time-domain numerical model was established using a higher-order boundary element method(HOBEM).By assuming small flow velocities,the velocity potential could be expressed...To study wave-current actions on 3-D bodies a time-domain numerical model was established using a higher-order boundary element method(HOBEM).By assuming small flow velocities,the velocity potential could be expressed for linear and higher order components by perturbation expansion.A 4th-order Runge-Kutta method was applied for time marching.An artificial damping layer was adopted at the outer zone of the free surface mesh to dissipate scattering waves.Validation of the numerical method was carried out on run-up,wave exciting forces,and mean drift forces for wave-currents acting on a bottom-mounted vertical cylinder.The results were in close agreement with the results of a frequency-domain method and a published time-domain method.The model was then applied to compute wave-current forces and run-up on a Seastar mini tension-leg platform.展开更多
The equivalent stress fundamental solution for the elastoplastic dynamic plane strain problem is proposed to transform the virtual work in the third direction to the plane.Subsequently,based on Betti reciprocal theore...The equivalent stress fundamental solution for the elastoplastic dynamic plane strain problem is proposed to transform the virtual work in the third direction to the plane.Subsequently,based on Betti reciprocal theorem,by adopting the time dependent fundamental solutions in terms of displacement,traction and equivalent stress,the boundary integral equations for dynamic elastoplastic analysis for the plane strain problem are established.The establishment procedures for the displacement and the stress boundary integral equations,together with the stress equation at boundary points,are presented in details,while the standard discretization both in time and space under the frame of time domain boundary element method(TD-BEM)and the solution of the algebraic equations are also briefly stated.Two verification examples are presented from different viewpoints,for elastic and elastoplastic analysis,for 1-D and 2-D geometries,and for finite and infinite domains.The TD-BEM formulation for dynamic elastoplastic analysis is presented for the plane strain problem as an example,where the formulation is also applicable for the plane stress problem by properly transforming the elastic constants and adopting the corresponding fundamental solutions.展开更多
基金The project was financially supported by the National Natural Science Foundation of China under the Grant No. 19732004 the National Science Fund for Distinguished Young Scholars under the Grant No. 50029002
文摘A time-domain method is applied to simulate nonlinear wave diffraction around a surface piercing 3-D arbitrary body. The method involves the application of Taylor series expansions and the use of perturbation procedure to establish the corresponding boundary value problems with respect to a time-independent fluid domain. A boundary element method based on B-spline expansion is used to calculate the wave field at each time step, and the free surface boundary condition is satisfied to the second order of wave steepness by a numerical integration in time. An artificial damping layer is adopted on the free surface for the removal of wave reflection from the outer boundary. As an illustration, the method is used to compute the second-order wave forces and run-up on a surface-piercing circular cylinder. The present method is found to be accurate, computationally efficient, and numerically stable.
基金supported by the National Natural Science Foundation of China (Grant No.50439010)the Key Project of the Ministry of Education of China (Grant No.305003)
文摘The time domain responses of the tunnel element under wave actions during its immersion are investigated based on the linear wave diffraction theory. The integral equation is derived by using the time-domain Green function that satisfies the free water surface condition in the finite water depth, and is solved by the boundary element method. The motion equations of the tunnel element are solved by the fourth order Runge-Kutta method. A comparison between the computed and measured results reveals that the numerical model can effectively simulate the motion responses of the tunnel element and the cable tensions when the motions of the tunnel element are within some limit. Taking the tunnel element of 100 m in length, 15 m in width and 10 m in height as an example, the computational results of the motion responses of the tunnel element and the cable tensions in different immersing depths are obtained under different incident wave conditions.
基金sponsored by the Fundamental Research Developing Association for Shipbuilding and Offshore(REDAS)the Special Coordination Funds for Promoting Science and Technology,Ministry of Education,Culture,Sports,Science and Technology(MEXT),Japan
文摘A 3-D iterative Rankine Boundary Element Method （BEM） for seakeeping problem in time domain is developed in the framework of linear potential theory. Waves generated by both submerged and surface-piercing bodies moving at a constant forward speed in otherwise calm water, and the resultant steady wave pattern, wave profile and resistance are computed to validate this newly-developed code. A rectangular computational domain moving with the same forward speed as the body is introduced, in which an artificial damping beach is installed at an outer portion of the free surface except the downstream side for satisfying the radiation condition. The velocity potential on the ship hull and the normal velocity on the free surface are obtained directly by solving the boundary integral equation, with the Rankine source used as the kernel function. An iterative time-marching scheme is employed for updating both kinematic and dynamic free surface boundary conditions to stabilize the calculation. Extensive results including the wave patterns, wave profiles and wave resistances for a submerged spheroid and a Wigley hull with forward speed are presented to validate the efficiency of the proposed 3-D time-domain higher-order approach. Finally, the sensitivity of ship-generated waves to the water depth is investigated. Computed results show satisfactory agreement with the corresponding experimental data and other numerical solutions.
文摘Environmental effects have an important influence on Offshore Wind Turbine (OWT) power generation efficiency and the structural stability of such turbines. In this study, we use an in-house Boundary Element (BEM)-panMARE code-to simulate the unsteady flow behavior of a full OWT with various combinations of aerodynamic and hydrodynamic loads in the time domain. This code is implemented to simulate potential flows for different applications and is based on a three-dimensional first-order panel method. Three different OWT configurations consisting of a generic 5 MW NREL rotor with three different types of foundations (Monopile, Tripod, and Jacket) are investigated. These three configurations are analyzed using the RANSE solver which is carried out using ANSYS CFX for validating the corresponding results. The simulations are performed under the same environmental atmospheric wind shear and rotor angular velocity, and the wave properties are wave height of 4 m and wave period of 7.16 s. In the present work, wave environmental effects were investigated firstly for the two solvers, and good agreement is achieved. Moreover, pressure distribution in each OWT case is presented, including detailed information about local flow fields. The time history of the forces at inflow direction and its moments around the mudline at each OWT part are presented in a dimensionless form with respect to the mean value of the last three loads and the moment amplitudes obtained from the BEM code, where the contribution of rotor force is lower in the tripod case and higher in the jacket case and the calculated hydrodynamic load that effect on jacket foundation type is lower than other two cases.
基金the National Key Basic Research and Development Program under Grant No.2002CB412709
文摘A seismic free field input formulation of the coupling procedure of the finite elelnent(FE)and the scaled boundary finite-element(SBFE)is proposed to perform the unbounded soil-structure interaction analysis in time domain. Based on the substructure technique,seismic excitation of the soil-structure system is represented by the free-field motion of an elastic half-space.To reduce the computational effort,the acceleration unit-impulse response function of the unbounded soil is decomposed into two functions:linear and residual.The latter converges to zero and can be truncated as required. With the prescribed tolerance parameter,the balance between accuracy and effMency of the procedure can be controlled. The validity of the model is verified by the scattering analysis of a hemi-spherical canyon subjected to plane harmonic P,SV and SH wave incidence.Numerical results show that the new procedure is very efficient for seismic problems within a nor- real range of frequency.The coupling procedure presented herein can be applied to linear and nonlinear earthquake re- sponse analysis of practical structures which are built on unbounded soil.
基金supported by the Science Council under the Project Nos.NSC-95-2221-E-019-075-MY3(CRC)andNSC-97-2221-E-236-011-(RSS)
文摘Developing serpent-type wave generators to generate solitary waves in a 3D-basin was investigated in this study. Based on the Lagrangian description with time-marching procedures and finite differences of the time derivative, a 3D multiple directional wave basin with multidirectional piston wave generators was developed to simulate ocean waves by using BEM with quadrilateral elements, and to simulate wave-caused problems with fully nonlinear water surface conditions. The simulations of perpendicular solitary waves were conducted in the first instance to verify this scheme. Furthermore, the comparison of the waveform variations confirms that the estimation of 3D solitary waves is a feasible scheme.
基金This research investigation was supported by the National Natural Science Foundation of China(Grant No.51678248 and Grant No.51878296)the Fundamental Research Funds for the Central Universities.And sincere thanks also to State Key Lab of Subtropical Building Science,South China University of Technology under Grant No.2017KB15 and the Open Research Fund of State Key Laboratory of Simulation and Regulation of Water Cycle in River Basin under Grant No.IWHRSKL-KF201818.
文摘Based on the modified scale boundary finite element method and continued fraction solution,a high-order doubly asymptotic transmitting boundary(DATB)is derived and extended to the simulation of vector wave propagation in complex layered soils.The high-order DATB converges rapidly to the exact solution throughout the entire frequency range and its formulation is local in the time domain,possessing high accuracy and good efficiency.Combining with finite element method,a coupled model is constructed for time-domain analysis of underground station-layered soil interaction.The coupled model is divided into the near and far field by the truncated boundary,of which the near field is modelled by FEM while the far field is modelled by the high-order DATB.The coupled model is implemented in an open source finite element software,OpenSees,in which the DATB is employed as a super element.Numerical examples demonstrate that results of the coupled model are stable,accurate and efficient compared with those of the extended mesh model and the viscous-spring boundary model.Besides,it has also shown the fitness for long-time seismic response analysis of underground station-layered soil interaction.Therefore,it is believed that the coupled model could provide a new approach for seismic analysis of underground station-layered soil interaction and could be further developed for engineering.
基金the National Natural Science Foundation of China and the Chinese State Edu-cation Commission Grant.
文摘The unsteady flow around a ship advancing at a constant speed was studied in time domain by using a higher-order boundary element method. The Rankine source was employed. The ship hull and free surface were discretized into 8-node isoparametric elements. A modified Green formula was used to overcome the numerical difficulty arising from singular or quasi-singular integral equation. In the case of τ>0.25, a rigid wall condition was applied to the upstream and a numerical beach was applied to the downstream to absorb the outgoing waves. The numerical results are satisfactory compared with the experimental results.
基金Supported by the National Natural Science Foundation of China under (Grant No.107 72040,50709005 and 50921001)the Major National Science and Technology Projects of China under (Grant No.2008ZX05026-02)the Open Fund of State Key Laboratory of Ocean Engineering
文摘To study wave-current actions on 3-D bodies a time-domain numerical model was established using a higher-order boundary element method(HOBEM).By assuming small flow velocities,the velocity potential could be expressed for linear and higher order components by perturbation expansion.A 4th-order Runge-Kutta method was applied for time marching.An artificial damping layer was adopted at the outer zone of the free surface mesh to dissipate scattering waves.Validation of the numerical method was carried out on run-up,wave exciting forces,and mean drift forces for wave-currents acting on a bottom-mounted vertical cylinder.The results were in close agreement with the results of a frequency-domain method and a published time-domain method.The model was then applied to compute wave-current forces and run-up on a Seastar mini tension-leg platform.
基金The authors would like to acknowledge the financial support provided by Hebei Education Department(Grant QN2020135)the National Key R&D Program of China(Grants 2019YFC1511105 and 2019YFC1511104)the National Natural Science Foundation of China(Grant 51778193).
文摘The equivalent stress fundamental solution for the elastoplastic dynamic plane strain problem is proposed to transform the virtual work in the third direction to the plane.Subsequently,based on Betti reciprocal theorem,by adopting the time dependent fundamental solutions in terms of displacement,traction and equivalent stress,the boundary integral equations for dynamic elastoplastic analysis for the plane strain problem are established.The establishment procedures for the displacement and the stress boundary integral equations,together with the stress equation at boundary points,are presented in details,while the standard discretization both in time and space under the frame of time domain boundary element method(TD-BEM)and the solution of the algebraic equations are also briefly stated.Two verification examples are presented from different viewpoints,for elastic and elastoplastic analysis,for 1-D and 2-D geometries,and for finite and infinite domains.The TD-BEM formulation for dynamic elastoplastic analysis is presented for the plane strain problem as an example,where the formulation is also applicable for the plane stress problem by properly transforming the elastic constants and adopting the corresponding fundamental solutions.
