In 3D frequency domain seismic forward and inversion calculation,the huge amount of calculation and storage is one of the main factors that restrict the processing speed and calculation efficiency.The frequency domain...In 3D frequency domain seismic forward and inversion calculation,the huge amount of calculation and storage is one of the main factors that restrict the processing speed and calculation efficiency.The frequency domain finite-difference forward simulation algorithm based on the acoustic wave equation establishes a large bandwidth complex matrix according to the discretized acoustic wave equation,and then the frequency domain wave field value is obtained by solving the matrix equation.In this study,the predecessor's optimized five-point method is extended to a 3D seven-point finite-difference scheme,and then a perfectly matched layer absorbing boundary condition(PML)is added to establish the corresponding matrix equation.In order to solve the complex matrix,we transform it to the equivalent real number domain to expand the solvable range of the matrix,and establish two objective functions to transform the matrix solving problem into an optimization problem that can be solved using gradient methods,and then use conjugate gradient algorithm to solve the problem.Previous studies have shown that in the conjugate gradient algorithm,the product of the matrix and the vector is the main factor that affects the calculation efficiency.Therefore,this study proposes a method that transform bandwidth matrix and vector product problem into some equivalent vector and vector product algorithm,thereby reducing the amount of calculation and storage.展开更多
A numerical framework was proposed for the seismic analysis of underground structures in layered ground under inclined P-SV waves.The free-field responses are first obtained using the stiffness matrix method based on ...A numerical framework was proposed for the seismic analysis of underground structures in layered ground under inclined P-SV waves.The free-field responses are first obtained using the stiffness matrix method based on plane-wave assumptions.Then,the domain reduction method was employed to reproduce the wavefield in the numerical model of the soil–structure system.The proposed numerical framework was verified by providing comparisons with analytical solutions for cases involving free-field responses of homogeneous ground,layered ground,and pressure-dependent heterogeneous ground,as well as for an example of a soil–structure interaction simulation.Compared with the viscous and viscous-spring boundary methods adopted in previous studies,the proposed framework exhibits the advantage of incorporating oblique incident waves in a nonlinear heterogeneous ground.Numerical results show that SV-waves are more destructive to underground structures than P-waves,and the responses of underground structures are significantly affected by the incident angles.展开更多
In this paper, the vibration characteristics of the structure in the finite fluid domain are analyzed using a coupled finite element method. The added mass matrix is calculated with finite element method （FEM） by 8-...In this paper, the vibration characteristics of the structure in the finite fluid domain are analyzed using a coupled finite element method. The added mass matrix is calculated with finite element method （FEM） by 8-node acoustic fluid elements. The vibration characteristics of the structure in the finite fluid domain are calculated combining structure FEM mass matrix. By writing relevant programs, the numerical analysis on vibration characteristics of a submerged cantilever rectangular plate in finite fluid domain and loaded ship model is performed. A modal identification experiment for the loaded ship model in air and in water is conducted and the experiment results verify the reliability of the numerical analysis. The numerical method can be used for further research on vibration characteristics and acoustic radiation problems of the structure in the finite fluid domain.展开更多
基金supported by the National Natural Science Foundation of China(Project U1901602&41790465)Key Special Project for Introduced Talents Team of Southern Marine Science and Engineering Guangdong Laboratory(Guangzhou)(GML2019ZD0203)+2 种基金Shenzhen Key Laboratory of Deep Offshore Oil and Gas Exploration Technology(Grant No.ZDSYS20190902093007855)Shenzhen Science and Technology Program(Grant No.KQTD20170810111725321)the leading talents of Guangdong province program(Grant No.2016LJ06N652).
文摘In 3D frequency domain seismic forward and inversion calculation,the huge amount of calculation and storage is one of the main factors that restrict the processing speed and calculation efficiency.The frequency domain finite-difference forward simulation algorithm based on the acoustic wave equation establishes a large bandwidth complex matrix according to the discretized acoustic wave equation,and then the frequency domain wave field value is obtained by solving the matrix equation.In this study,the predecessor's optimized five-point method is extended to a 3D seven-point finite-difference scheme,and then a perfectly matched layer absorbing boundary condition(PML)is added to establish the corresponding matrix equation.In order to solve the complex matrix,we transform it to the equivalent real number domain to expand the solvable range of the matrix,and establish two objective functions to transform the matrix solving problem into an optimization problem that can be solved using gradient methods,and then use conjugate gradient algorithm to solve the problem.Previous studies have shown that in the conjugate gradient algorithm,the product of the matrix and the vector is the main factor that affects the calculation efficiency.Therefore,this study proposes a method that transform bandwidth matrix and vector product problem into some equivalent vector and vector product algorithm,thereby reducing the amount of calculation and storage.
基金supported by the National Natural Science Foundation of China(Grant Nos.41922059,42177134,and 51778487)Fundamental Research Funds for the Central Universities,CHD(300102262506)Top Discipline Plan of Shanghai Universities-Class I.
文摘A numerical framework was proposed for the seismic analysis of underground structures in layered ground under inclined P-SV waves.The free-field responses are first obtained using the stiffness matrix method based on plane-wave assumptions.Then,the domain reduction method was employed to reproduce the wavefield in the numerical model of the soil–structure system.The proposed numerical framework was verified by providing comparisons with analytical solutions for cases involving free-field responses of homogeneous ground,layered ground,and pressure-dependent heterogeneous ground,as well as for an example of a soil–structure interaction simulation.Compared with the viscous and viscous-spring boundary methods adopted in previous studies,the proposed framework exhibits the advantage of incorporating oblique incident waves in a nonlinear heterogeneous ground.Numerical results show that SV-waves are more destructive to underground structures than P-waves,and the responses of underground structures are significantly affected by the incident angles.
基金Supported by the National Natural Science Foundation of China （No. 51079027）.
文摘In this paper, the vibration characteristics of the structure in the finite fluid domain are analyzed using a coupled finite element method. The added mass matrix is calculated with finite element method （FEM） by 8-node acoustic fluid elements. The vibration characteristics of the structure in the finite fluid domain are calculated combining structure FEM mass matrix. By writing relevant programs, the numerical analysis on vibration characteristics of a submerged cantilever rectangular plate in finite fluid domain and loaded ship model is performed. A modal identification experiment for the loaded ship model in air and in water is conducted and the experiment results verify the reliability of the numerical analysis. The numerical method can be used for further research on vibration characteristics and acoustic radiation problems of the structure in the finite fluid domain.
