After a brief review of studies on artificial boundaries in dynamic soil-structure interaction, a three-dimensional viscous-spring artificial boundary （VSAB） in the time domain is developed in this paper. First, the...After a brief review of studies on artificial boundaries in dynamic soil-structure interaction, a three-dimensional viscous-spring artificial boundary （VSAB） in the time domain is developed in this paper. First, the 3D VSAB equations in the normal and tangential directions are derived based on the elastic wave motion theory. Secondly, a numerical simulation technique of wave motion equations along with the VSAB condition in the time domain is studied. Finally, numerical examples of some classical elastic wave motion problems are presented and the results are compared with the associated theoretical solutions, demonstrating that high precision and adequate stability can be achieved by using the proposed 3D VSAB. The proposed 3D VSAB can be conveniently incorporated in the general finite element program, which is commonly used to study dynamic soil-structure interaction problems.展开更多
The finite element artificial transmitting boundary method is employed here to treat the near field scattering of a cylindrical wave from an irregular cylinder. A comparison is made between this method and the analy...The finite element artificial transmitting boundary method is employed here to treat the near field scattering of a cylindrical wave from an irregular cylinder. A comparison is made between this method and the analytical one. And then examples are given to demonstrate the solution of several problems of the irregular object scattering. The method can not only produce clear physical pictures, but can efficiently handle many complicated scattering problems.展开更多
An improved version of the regular boundary element method, the artificial boundary node approach, is derived. A simple contact algorithm is designed and implemented into the direct boundary element, regular boundary ...An improved version of the regular boundary element method, the artificial boundary node approach, is derived. A simple contact algorithm is designed and implemented into the direct boundary element, regular boundary element and artificial boundary node approaches. The exisiting and derived approaches are tested using some case studies. The results of the artificial boundary node approach are compared with those of the existing boundary element program, the regular element approach, ANSYS and analytical solution whenever possible. The results show the effectiveness of the artificial boundary node approach for a wider range of boundary offsets.展开更多
The theory of perfectly matched layer （PML） artificial boundary condition （ABC）, which is characterized by absorption any wave motions with arbitrary frequency and arbitrarily incident angle, is introduced. The co...The theory of perfectly matched layer （PML） artificial boundary condition （ABC）, which is characterized by absorption any wave motions with arbitrary frequency and arbitrarily incident angle, is introduced. The construction process of PML boundary based on elastodynamic partial differential equation （PDE） system is developed. Combining with velocity-stress hybrid finite element formulation, the applicability of PML boundary is investigated and the numerical reflection of PML boundary is estimated. The reflectivity of PML and multi-transmitting formula （MTF） boundary is then compared based on body wave and surface wave simulations. The results show that although PML boundary yields some reflection, its absorption performance is superior to MTF boundary in the numerical simulations of near-fault wave propagation, especially in comer and large angle grazing incidence situations. The PML boundary does not arise any unstable phenomenon and the stability of PML boundary is better than MTF boundary in hybrid finite element method. For a specified problem and analysis tolerance, the computational efficiency of PML boundary is only a little lower than MTF boundary.展开更多
Anew artificial boundary model based on multi-directional transmitting and viscous-spring artificial boundary theories is proposed to absorb stress waves in a saturated soil foundation in dynamic analysis. Since shear...Anew artificial boundary model based on multi-directional transmitting and viscous-spring artificial boundary theories is proposed to absorb stress waves in a saturated soil foundation in dynamic analysis. Since shear waves （S-waves） are the same in a saturated soil foundation and a single-phase medium foundation, a tangential visco-elastic boundary condition for a single-phase medium foundation can also be used for saturated soil foundations. Thus, the purpose of the artificial boundary proposed in this paper is primarily to absorb two types of P-waves in a saturated soil foundation. The main idea is that the stress of the P-waves in the saturated soil foundation is decomposed into two types. The first type of stress, δra＇ is absorbed by the first artificial boundary. The second type of stress, δrb, is balanced by the stress generated by the second artificial boundary. Ultimately, both types of P-waves （fast-P-waves and slow-P-waves） are absorbed by the artificial boundary model proposed in this paper. In particular, note that the fast-P-waves and slow-P-waves are absorbed at the position of the first boundary. Thus, the artificial boundary model proposed herein can simultaneously absorb P-fast waves, P-slow waves and shear waves. Finally, a numerical example is given to examine the proposed artificial boundary model, and the results show that it is very accurate.展开更多
An explicit finite element-finite difference method for analyzing the effects of two-dimensional visco-elastic localtopography on earthquake ground motion is prOPosed in this paper. In the method, at first, the finite...An explicit finite element-finite difference method for analyzing the effects of two-dimensional visco-elastic localtopography on earthquake ground motion is prOPosed in this paper. In the method, at first, the finite elementdiscrete model is formed by using the artificial boundary and finite element method, and the dynamic equationsof local nodes in the discrete model are obtained according to the theory of the special finite element method similar to the finite difference method, and then the explicit step-by-step integration formulas are presented by usingthe explicit difference method for solving the visco-elastic dynamic equation and Generalized Multi-transmittingBoundary. The method has the advantages of saving computing time and computer memory space, and it is suitable for any case of topography and has high computing accuracy and good computing stability.展开更多
Based on presumed active fault and corresponding model, this paper predicted the near-fault ground motion filed of a scenario earthquake （Mw=6 3/4 ） in an active fault by the explicit finite element method in combin...Based on presumed active fault and corresponding model, this paper predicted the near-fault ground motion filed of a scenario earthquake （Mw=6 3/4 ） in an active fault by the explicit finite element method in combination with the source time function with improved transmitting artificial boundary and with high-frequency vibration contained. The results indicate that the improved artificial boundary is stable in numerical computation and the predicted strong ground motion has a consistent characteristic with the observed motion.展开更多
The mixed finite element method is used to solve the exterior Poisson equations with higher-order local artificial boundary conditions in 3-D space. New unknowns are introduced to reduce the order of the derivatives ...The mixed finite element method is used to solve the exterior Poisson equations with higher-order local artificial boundary conditions in 3-D space. New unknowns are introduced to reduce the order of the derivatives of the unknown to two. The result is an equivalent mixed variational problem which was solved using bilinear finite elements. The primary advantage is that special finite elements are not needed on the adjacent layer of the artificial boundary for the higher-order derivatives. Error estimates are obtained for some local artificial boundary conditions with prescibed orders. A numerical example demonstrates the effectiveness of this method.展开更多
An a posteriori error estimator is obtained for a nonconforming finite element approximation of a linear elliptic problem, which is derived from a corresponding unbounded domain problem by applying a nonlocal approxim...An a posteriori error estimator is obtained for a nonconforming finite element approximation of a linear elliptic problem, which is derived from a corresponding unbounded domain problem by applying a nonlocal approximate artificial boundary condition. Our method can be easily extended to obtain a class of a posteriori error estimators for various conforming and nonconforming finite element approximations of problems with different artificial boundary conditions. The reliability and efficiency of our a posteriori error estimator are rigorously proved and are verified by numerical examples.展开更多
The present paper focuses on the wave radiation by an oscillating body with six degrees of freedom by using the DtN artifi-cial boundary condition.The artificial boundary is usually selected as a circle or spherical s...The present paper focuses on the wave radiation by an oscillating body with six degrees of freedom by using the DtN artifi-cial boundary condition.The artificial boundary is usually selected as a circle or spherical surface to solve various types of fields,such as sound waves or electromagnetic waves,provided that the considered domain is infinite or unbounded in all directions.However,the substantial wave motion is considered in water of finite depth,that is,the fluid domain is bounded vertically but unbounded horizon-tally.Thus,the DtN boundary condition is given on an artificial cylindrical surface,which divides the water domain into an interior and exterior region.The boundary integral equation is adopted to implement the present model.In the case of a floating cylinder,the results of hydrodynamic coefficients of a chamfer box are discussed.展开更多
In earlier approach, the 2-D acoustical field profiles on the substrate region are often calculated with BPM. In this paper, we present a new approach based on the finite element -artificial transmitting boundary meth...In earlier approach, the 2-D acoustical field profiles on the substrate region are often calculated with BPM. In this paper, we present a new approach based on the finite element -artificial transmitting boundary method and calculate acoustical field on the substrate region.展开更多
We consider the linearized incompressible Navier-Stokes (Oseen) equations in a flat channel. A sequence of approximations to the exact boundary condition at an artificial boundary is derived. Then the original problem...We consider the linearized incompressible Navier-Stokes (Oseen) equations in a flat channel. A sequence of approximations to the exact boundary condition at an artificial boundary is derived. Then the original problem is reduced to a boundary value problem in a bounded domain, which is well-posed. A finite element approximation on the bounded domain is given, furthermore the error estimate of the finite element approximation is obtained. Numerical example shows that our artificial boundary conditions are very effective.展开更多
We consider the numerical approximations of the three-dimensional steady potential flow around a body moving in a liquid of finite constant depth at constant speed and distance below a free surface in a channel. One v...We consider the numerical approximations of the three-dimensional steady potential flow around a body moving in a liquid of finite constant depth at constant speed and distance below a free surface in a channel. One vertical side is introduced as the upstream artificial boundary and two vertical sides are introduced as the downstream artificial boundaries. On the artificial boundaries, a sequence of high-order global artificial boundary conditions are given. Then the original problem is reduced to a problem defined on a finite computational domain, which is equivalent to a variational problem. After solving the variational problem by the finite element method, we obtain the numerical approximation of the original problem. The numerical examples show that the artificial boundary conditions given in this paper are very effective.展开更多
The calculation of wave resistance for a ship moving at constant speed near a free surface is considered. This wave resistance is calculated with a linearized steady potential model. To deal with the unboundedness of...The calculation of wave resistance for a ship moving at constant speed near a free surface is considered. This wave resistance is calculated with a linearized steady potential model. To deal with the unboundedness of the physical domain in the potential flow problem, we introduce one vertical side as an artificial upstream boundary and two vertical sides as the artificial downstream boundaries. On the artificial boundaries, a sequence of high-order global artificial boundary conditions are given. Then the potential flow problem is reduced to a problem defined on a finite computational domain, which is equivalent to a variational problem. The solution of the variational problem by the finite element method gives the numerical approximation of the potential flow around the ship, which was used to calculate the wave resistance. The numerical examples show the accuracy and efficiency of the proposed numerical scheme.展开更多
The artificial boundary method is applied to solve three-dimensional exterior problems. Two kind of rotating ellipsoids are chosen as the artificial boundaries and the exact artificial boundary conditions are derived ...The artificial boundary method is applied to solve three-dimensional exterior problems. Two kind of rotating ellipsoids are chosen as the artificial boundaries and the exact artificial boundary conditions are derived explicitly in terms of an infinite series. Then the well-posedness of the coupled variational problem is obtained. It is found that error estimates derived depend on the mesh size, truncation term and the location of the artificial boundary. Three numerical examples are presented to demonstrate the effectiveness and accuracy of the proposed method.展开更多
基金National Natural Science Foundation of ChinaUnder Grant No.50478014Special Funds for Major State Basic Research Project Under Grant No.2002CB412706Research Funds from National Civil Defense Oficce of Chinafor the Tenth Five-year Plan。
文摘After a brief review of studies on artificial boundaries in dynamic soil-structure interaction, a three-dimensional viscous-spring artificial boundary （VSAB） in the time domain is developed in this paper. First, the 3D VSAB equations in the normal and tangential directions are derived based on the elastic wave motion theory. Secondly, a numerical simulation technique of wave motion equations along with the VSAB condition in the time domain is studied. Finally, numerical examples of some classical elastic wave motion problems are presented and the results are compared with the associated theoretical solutions, demonstrating that high precision and adequate stability can be achieved by using the proposed 3D VSAB. The proposed 3D VSAB can be conveniently incorporated in the general finite element program, which is commonly used to study dynamic soil-structure interaction problems.
文摘The finite element artificial transmitting boundary method is employed here to treat the near field scattering of a cylindrical wave from an irregular cylinder. A comparison is made between this method and the analytical one. And then examples are given to demonstrate the solution of several problems of the irregular object scattering. The method can not only produce clear physical pictures, but can efficiently handle many complicated scattering problems.
文摘An improved version of the regular boundary element method, the artificial boundary node approach, is derived. A simple contact algorithm is designed and implemented into the direct boundary element, regular boundary element and artificial boundary node approaches. The exisiting and derived approaches are tested using some case studies. The results of the artificial boundary node approach are compared with those of the existing boundary element program, the regular element approach, ANSYS and analytical solution whenever possible. The results show the effectiveness of the artificial boundary node approach for a wider range of boundary offsets.
基金National Natural Science Foundation of China (50608024 and 50538050).
文摘The theory of perfectly matched layer （PML） artificial boundary condition （ABC）, which is characterized by absorption any wave motions with arbitrary frequency and arbitrarily incident angle, is introduced. The construction process of PML boundary based on elastodynamic partial differential equation （PDE） system is developed. Combining with velocity-stress hybrid finite element formulation, the applicability of PML boundary is investigated and the numerical reflection of PML boundary is estimated. The reflectivity of PML and multi-transmitting formula （MTF） boundary is then compared based on body wave and surface wave simulations. The results show that although PML boundary yields some reflection, its absorption performance is superior to MTF boundary in the numerical simulations of near-fault wave propagation, especially in comer and large angle grazing incidence situations. The PML boundary does not arise any unstable phenomenon and the stability of PML boundary is better than MTF boundary in hybrid finite element method. For a specified problem and analysis tolerance, the computational efficiency of PML boundary is only a little lower than MTF boundary.
基金National Natural Science Foundation of China Under Grant Nos.51109029,51178081,51138001,51009020China Postdoctoral Science Foundation Under Grant No. 20110491535
文摘Anew artificial boundary model based on multi-directional transmitting and viscous-spring artificial boundary theories is proposed to absorb stress waves in a saturated soil foundation in dynamic analysis. Since shear waves （S-waves） are the same in a saturated soil foundation and a single-phase medium foundation, a tangential visco-elastic boundary condition for a single-phase medium foundation can also be used for saturated soil foundations. Thus, the purpose of the artificial boundary proposed in this paper is primarily to absorb two types of P-waves in a saturated soil foundation. The main idea is that the stress of the P-waves in the saturated soil foundation is decomposed into two types. The first type of stress, δra＇ is absorbed by the first artificial boundary. The second type of stress, δrb, is balanced by the stress generated by the second artificial boundary. Ultimately, both types of P-waves （fast-P-waves and slow-P-waves） are absorbed by the artificial boundary model proposed in this paper. In particular, note that the fast-P-waves and slow-P-waves are absorbed at the position of the first boundary. Thus, the artificial boundary model proposed herein can simultaneously absorb P-fast waves, P-slow waves and shear waves. Finally, a numerical example is given to examine the proposed artificial boundary model, and the results show that it is very accurate.
文摘An explicit finite element-finite difference method for analyzing the effects of two-dimensional visco-elastic localtopography on earthquake ground motion is prOPosed in this paper. In the method, at first, the finite elementdiscrete model is formed by using the artificial boundary and finite element method, and the dynamic equationsof local nodes in the discrete model are obtained according to the theory of the special finite element method similar to the finite difference method, and then the explicit step-by-step integration formulas are presented by usingthe explicit difference method for solving the visco-elastic dynamic equation and Generalized Multi-transmittingBoundary. The method has the advantages of saving computing time and computer memory space, and it is suitable for any case of topography and has high computing accuracy and good computing stability.
基金Heilongjiang Province Postdoctoral Science Foundation and China Earthquake Administration’s Tenth Five-year Plan Project
文摘Based on presumed active fault and corresponding model, this paper predicted the near-fault ground motion filed of a scenario earthquake （Mw=6 3/4 ） in an active fault by the explicit finite element method in combination with the source time function with improved transmitting artificial boundary and with high-frequency vibration contained. The results indicate that the improved artificial boundary is stable in numerical computation and the predicted strong ground motion has a consistent characteristic with the observed motion.
基金Supported by the National Natural Science Foundationof China(No.19772 0 2 2 )
文摘The mixed finite element method is used to solve the exterior Poisson equations with higher-order local artificial boundary conditions in 3-D space. New unknowns are introduced to reduce the order of the derivatives of the unknown to two. The result is an equivalent mixed variational problem which was solved using bilinear finite elements. The primary advantage is that special finite elements are not needed on the adjacent layer of the artificial boundary for the higher-order derivatives. Error estimates are obtained for some local artificial boundary conditions with prescibed orders. A numerical example demonstrates the effectiveness of this method.
基金supported by the Special Funds for Major State Basic Research Projects(2005CB321701)NSFC(10431050, 10571006 and 10528102)RFDP of China
文摘An a posteriori error estimator is obtained for a nonconforming finite element approximation of a linear elliptic problem, which is derived from a corresponding unbounded domain problem by applying a nonlocal approximate artificial boundary condition. Our method can be easily extended to obtain a class of a posteriori error estimators for various conforming and nonconforming finite element approximations of problems with different artificial boundary conditions. The reliability and efficiency of our a posteriori error estimator are rigorously proved and are verified by numerical examples.
文摘The present paper focuses on the wave radiation by an oscillating body with six degrees of freedom by using the DtN artifi-cial boundary condition.The artificial boundary is usually selected as a circle or spherical surface to solve various types of fields,such as sound waves or electromagnetic waves,provided that the considered domain is infinite or unbounded in all directions.However,the substantial wave motion is considered in water of finite depth,that is,the fluid domain is bounded vertically but unbounded horizon-tally.Thus,the DtN boundary condition is given on an artificial cylindrical surface,which divides the water domain into an interior and exterior region.The boundary integral equation is adopted to implement the present model.In the case of a floating cylinder,the results of hydrodynamic coefficients of a chamfer box are discussed.
基金This work was supported by the National Natural Science Fund of China (10084001) the Action Plan for Booming Education of Tianjin University.
文摘In earlier approach, the 2-D acoustical field profiles on the substrate region are often calculated with BPM. In this paper, we present a new approach based on the finite element -artificial transmitting boundary method and calculate acoustical field on the substrate region.
基金This work was supported by the Climbing Program of National Key Project of Foundation andDoctoral Program foundation of Instit
文摘We consider the linearized incompressible Navier-Stokes (Oseen) equations in a flat channel. A sequence of approximations to the exact boundary condition at an artificial boundary is derived. Then the original problem is reduced to a boundary value problem in a bounded domain, which is well-posed. A finite element approximation on the bounded domain is given, furthermore the error estimate of the finite element approximation is obtained. Numerical example shows that our artificial boundary conditions are very effective.
基金This work was supported partly by the Special Funds for Major State Basic Research Projects of China and the National Natural Science Foundation of China.Computation was supported by the State Key Lab.of the Scientific nd Engineering Computing in China
文摘We consider the numerical approximations of the three-dimensional steady potential flow around a body moving in a liquid of finite constant depth at constant speed and distance below a free surface in a channel. One vertical side is introduced as the upstream artificial boundary and two vertical sides are introduced as the downstream artificial boundaries. On the artificial boundaries, a sequence of high-order global artificial boundary conditions are given. Then the original problem is reduced to a problem defined on a finite computational domain, which is equivalent to a variational problem. After solving the variational problem by the finite element method, we obtain the numerical approximation of the original problem. The numerical examples show that the artificial boundary conditions given in this paper are very effective.
基金Supported partly by the Special Funds for Major StateBasic Research Projects of China(No.G19990 32 80 2 )and the National Natural Science Foundation of China(No. 19772 0 2 2 )
文摘The calculation of wave resistance for a ship moving at constant speed near a free surface is considered. This wave resistance is calculated with a linearized steady potential model. To deal with the unboundedness of the physical domain in the potential flow problem, we introduce one vertical side as an artificial upstream boundary and two vertical sides as the artificial downstream boundaries. On the artificial boundaries, a sequence of high-order global artificial boundary conditions are given. Then the potential flow problem is reduced to a problem defined on a finite computational domain, which is equivalent to a variational problem. The solution of the variational problem by the finite element method gives the numerical approximation of the potential flow around the ship, which was used to calculate the wave resistance. The numerical examples show the accuracy and efficiency of the proposed numerical scheme.
基金subsidized by the National Basic Research Program of China under the grant 2005CB321701the National Natural Science Foundation of China under the grant 10531080the Beijing Natural Science Foundation under the grant 1072009 and the Research Project of Zhejiang Ocean University (X08M013,X08Z04)
文摘The artificial boundary method is applied to solve three-dimensional exterior problems. Two kind of rotating ellipsoids are chosen as the artificial boundaries and the exact artificial boundary conditions are derived explicitly in terms of an infinite series. Then the well-posedness of the coupled variational problem is obtained. It is found that error estimates derived depend on the mesh size, truncation term and the location of the artificial boundary. Three numerical examples are presented to demonstrate the effectiveness and accuracy of the proposed method.
