The transmitting boundary condition is written in a compact form, which can be direct-ly incorporated into finite elements. Basic characteristics of discretization are analyzed throughstudies on wave motion in a one-d...The transmitting boundary condition is written in a compact form, which can be direct-ly incorporated into finite elements. Basic characteristics of discretization are analyzed throughstudies on wave motion in a one-dimensional discrete model and their differences from those in thecorresponding continuum. Tbe analysis leads to identifying a frequency band within which thesimulation is possible, and to a suggestion of using the lumped-mass finite element model forthe simulation. Mechanism of the oscillation instability is then illuminated in the frequencydomain by amplification at the artificial boundary and multi-reflection of wave motion in afinite discrete model. Based on understanding of the mechanism, a modified transmittingboundary condition is devised for eliminating the instability. The special stability criterion forthe modified boundary is finally presented for the one-dimensional model.展开更多
A discrete local transmitting boundary is combined with a lumped-mass finite element technique for simulating steady SH wave motion in a layered and unbounded medium.The combination decouples any node of the finite el...A discrete local transmitting boundary is combined with a lumped-mass finite element technique for simulating steady SH wave motion in a layered and unbounded medium.The combination decouples any node of the finite element model from the others except those in its neighborhood,and an effective algorithm of the Gauss elimination in implementation of the method is then devised to reduce the main memory and computing time.The method and its accuracy are first analysed via simulating steady SH wave motion in layered elastic media,the effective algorithm is then introduced and illustrated by simple examples.展开更多
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.展开更多
Based on Biot theory of two-phase anisotropic media and Hamilton theory about dynamic problem,finite element equations of elastic wave propagation in two-phase anisotropic media are derived in this paper.Numerical sol...Based on Biot theory of two-phase anisotropic media and Hamilton theory about dynamic problem,finite element equations of elastic wave propagation in two-phase anisotropic media are derived in this paper.Numerical solution of finite element equations is given.Finally,Properties of elastic wave propagation are observed and analyzed through FEM modeling.展开更多
This paper briefly reviews the characteristics and major processes of the explicit finite element method in modeling the near-fault ground motion field. The emphasis is on the finite element-related problems in the fi...This paper briefly reviews the characteristics and major processes of the explicit finite element method in modeling the near-fault ground motion field. The emphasis is on the finite element-related problems in the finite fault source modeling. A modified kinematic source model is presented, in which vibration with some high frequency components is introduced into the traditional slip time function to ensure that the source and ground motion include sufficient high frequency components. The model presented is verified through a simple modeling example. It is shown that the predicted near-fault ground motion field exhibits similar characteristics to those observed in strong motion records, such as the hanging wall effect, vertical effect, fling step effect and velocity pulse effect, etc.展开更多
The analysis method of lattice dynamics in classical physics is extended to study the properties of in-plane wave motion in the hybrid-mass finite element model in this paper. The dispersion equations of P and SV wave...The analysis method of lattice dynamics in classical physics is extended to study the properties of in-plane wave motion in the hybrid-mass finite element model in this paper. The dispersion equations of P and SV waves in the discrete model are first obtained by means of separating the characteristic equation of the motion equation, and then used to analyse the properties of P-and SV-homogeneous, inhomogeneous waves and other types of motion in the model. The dispersion characters, cut-off frequencies of P and SV waves, the polarization drift and appendent anisotropic property of wave motion caused by the discretization are finally discussed.展开更多
In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform ...In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform cubic grids, the approach used to establish stable formulas with 2M-order accuracy is discussed in detail, with M being a positive integer, and is illustrated by establishing second order (M=1) recursion formulas. The theoretical results presented in this paper are demonstrated through numerical testing.展开更多
In this paper, a method to develop a hierarchy of explicit recursion formulas for numerical simulation in an irregular grid for scalar wave equations is presented and its accuracy is illustrated via 2-D and 1-D models...In this paper, a method to develop a hierarchy of explicit recursion formulas for numerical simulation in an irregular grid for scalar wave equations is presented and its accuracy is illustrated via 2-D and 1-D models. Approaches to develop the stable formulas which are of 2M-order accuracy in both time and space with Mbeing a positive integer for regular grids are discussed and illustrated by constructing the second order (M= 1) and the fourth order (M = 2) recursion formulas.展开更多
A three-step finite element method (FEM) together with Large Eddy Simulation (LES) is applied to incompressible turbulent flow around seabed pipelines at relatively high Reynolds numbers. Both two dimensional and thre...A three-step finite element method (FEM) together with Large Eddy Simulation (LES) is applied to incompressible turbulent flow around seabed pipelines at relatively high Reynolds numbers. Both two dimensional and three-dimensional numerical simulation is carried out to determine the three-dimensional effect. The results of numerical simulation agree quite well with the wave forces acting on pipeline models measured in physical model test.展开更多
A two-dimensional numerical model based on the Navier-Stokes equations and computational Lagrangian-Eulerian advection remap-volume of fluid (CLEAR-VOF) method was developed to simulate wave and flow problems. The N...A two-dimensional numerical model based on the Navier-Stokes equations and computational Lagrangian-Eulerian advection remap-volume of fluid (CLEAR-VOF) method was developed to simulate wave and flow problems. The Navier-Stokes equations were discretized with a three-step finite element method that has a third-order accuracy. In the CLEAR-VOF method, the VOF function F was calculated in the Lagrangian manner and allowed the complicated free surface to be accurately captured. The propagation of regular waves and solitary waves over a flat bottom, and shoaling and breaking of solitary waves on two different slopes were simulated with this model, and the numerical results agreed with experimental data and theoretical solutions. A benchmark test of dam-collapse flow was also simulated with an unstructured mesh, and the capability of the present model for wave and flow simulations with unstructured meshes, was verified. The results show that the model is effective for numerical simulation of wave and flow problems with both structured and unstructured meshes.展开更多
The cross-flow(CF)vortex-induced vibration(VIV)of a deepwater steep wave riser(SWR)subjected to uniform or shear flow loads is investigated numerically.The model is based on a three-dimensional(3D)nonlinear elastic ro...The cross-flow(CF)vortex-induced vibration(VIV)of a deepwater steep wave riser(SWR)subjected to uniform or shear flow loads is investigated numerically.The model is based on a three-dimensional(3D)nonlinear elastic rod theory coupled with a wake oscillator model.In this numerical simulation,the nonlinear motion equations of the riser with large deformation features are established in a global coordinate system to avoid the transformation between global and local coordinate systems,and are discretized with the time-domain finite element method(FEM).A wakeoscillator model is employed to study the vortex shedding,and the lift force generated by the wake flow is described in a van der Pol equation.A Newmark-βiterative scheme is used to solve their coupling equation for the VIV response of the SWR.The developed model is validated against the existing experimental results for the VIV response of the top-tension riser(TTR).Then,the numerical simulations are executed to determine VIV characteristics of the SWR.The effects of both flow velocity and the spanwise length of the flow field on the drag coefficient in the inline(IL)direction and the lift coefficient in the CF direction are investigated systematically.The results indicate that compared with TTR,the low frequency and multi-modal vibration are the main components of the SWR due to the large deformation and flexible characteristics.For shear flow,the multi-frequency resonance dominates the VIV response of the SWR,especially at the hang-off segment.展开更多
In order to analyze the possibility of detecting defects in bend pipe using low-frequency ultrasonic guided wave, the propagation of T(0,1) mode and L(0,2) mode through straight-curved-straight pipe sections was studi...In order to analyze the possibility of detecting defects in bend pipe using low-frequency ultrasonic guided wave, the propagation of T(0,1) mode and L(0,2) mode through straight-curved-straight pipe sections was studied. FE(finite element) models of bend pipe without defects and those with defects were introduced to analyze energy distribution, mode transition and defect detection of ultrasonic guided wave. FE simulation results were validated by experiments of four different bend pipes with circumferential defects in different positions. It is shown that most energy of T(0,1) mode or L(0,2) mode focuses on extrados of bend but little passes through intrados of bend, and T(0,1) mode or L(0,2) mode is converted to other possible non-axisymmetric modes when propagating through the bend and the defect after bend respectively. Furthermore, L(0,2) mode is more sensitive to circumferential notch than T(0,1) mode. The results of this work are beneficial for practical testing of pipes.展开更多
It is well-known that artificial boundary conditions are crucial for the efficient and accurate computations of wavefields on unbounded domains. In this paper, we investigate stability analysis for the wave equation c...It is well-known that artificial boundary conditions are crucial for the efficient and accurate computations of wavefields on unbounded domains. In this paper, we investigate stability analysis for the wave equation coupled with the first and the second order absorbing boundary conditions. The computational scheme is also developed. The approach allows the absorbing boundary conditions to be naturally imposed, which makes it easier for us to construct high order schemes for the absorbing boundary conditions. A thirdorder Lagrange finite element method with mass lumping is applied to obtain the spatial discretization of the wave equation. The resulting scheme is stable and is very efficient since no matrix inversion is needed at each time step. Moreover, we have shown both abstract and explicit conditional stability results for the fully-discrete schemes. The results are helpful for designing computational parameters in computations. Numerical computations are illustrated to show the efficiency and accuracy of our method. In particular, essentially no boundary reflection is seen at the artificial boundaries.展开更多
If a traditional explicit numerical integration algorithm is used to solve motion equation in the finite element simulation of wave motion, the time-step used by numerical integration is the smallest time-step restric...If a traditional explicit numerical integration algorithm is used to solve motion equation in the finite element simulation of wave motion, the time-step used by numerical integration is the smallest time-step restricted by the stability criterion in computational region. However, the excessively small time-step is usually unnecessary for a large portion of computational region. In this paper, a varying time-step explicit numerical integration algorithm is introduced, and its basic idea is to use different time-step restricted by the stability criterion in different computational region. Finally, the feasibility of the algorithm and its effect on calculating precision are verified by numerical test.展开更多
文摘The transmitting boundary condition is written in a compact form, which can be direct-ly incorporated into finite elements. Basic characteristics of discretization are analyzed throughstudies on wave motion in a one-dimensional discrete model and their differences from those in thecorresponding continuum. Tbe analysis leads to identifying a frequency band within which thesimulation is possible, and to a suggestion of using the lumped-mass finite element model forthe simulation. Mechanism of the oscillation instability is then illuminated in the frequencydomain by amplification at the artificial boundary and multi-reflection of wave motion in afinite discrete model. Based on understanding of the mechanism, a modified transmittingboundary condition is devised for eliminating the instability. The special stability criterion forthe modified boundary is finally presented for the one-dimensional model.
文摘A discrete local transmitting boundary is combined with a lumped-mass finite element technique for simulating steady SH wave motion in a layered and unbounded medium.The combination decouples any node of the finite element model from the others except those in its neighborhood,and an effective algorithm of the Gauss elimination in implementation of the method is then devised to reduce the main memory and computing time.The method and its accuracy are first analysed via simulating steady SH wave motion in layered elastic media,the effective algorithm is then introduced and illustrated by simple examples.
基金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.
文摘Based on Biot theory of two-phase anisotropic media and Hamilton theory about dynamic problem,finite element equations of elastic wave propagation in two-phase anisotropic media are derived in this paper.Numerical solution of finite element equations is given.Finally,Properties of elastic wave propagation are observed and analyzed through FEM modeling.
文摘This paper briefly reviews the characteristics and major processes of the explicit finite element method in modeling the near-fault ground motion field. The emphasis is on the finite element-related problems in the finite fault source modeling. A modified kinematic source model is presented, in which vibration with some high frequency components is introduced into the traditional slip time function to ensure that the source and ground motion include sufficient high frequency components. The model presented is verified through a simple modeling example. It is shown that the predicted near-fault ground motion field exhibits similar characteristics to those observed in strong motion records, such as the hanging wall effect, vertical effect, fling step effect and velocity pulse effect, etc.
基金The project sponsored by the Earthquake Science Foundation under Contract No. 90141
文摘The analysis method of lattice dynamics in classical physics is extended to study the properties of in-plane wave motion in the hybrid-mass finite element model in this paper. The dispersion equations of P and SV waves in the discrete model are first obtained by means of separating the characteristic equation of the motion equation, and then used to analyse the properties of P-and SV-homogeneous, inhomogeneous waves and other types of motion in the model. The dispersion characters, cut-off frequencies of P and SV waves, the polarization drift and appendent anisotropic property of wave motion caused by the discretization are finally discussed.
基金China Postdoctoral Science Foundation Under Grant No.20100480321National Basic Research Program of China Under Grant No. 2007CB714200
文摘In this paper, an explicit method is generalized from 1D and 2D models to a 3D model for numerical simulation of wave motion, and the corresponding recursion formulas are developed for 3D irregular grids. For uniform cubic grids, the approach used to establish stable formulas with 2M-order accuracy is discussed in detail, with M being a positive integer, and is illustrated by establishing second order (M=1) recursion formulas. The theoretical results presented in this paper are demonstrated through numerical testing.
基金National Basic Research Program of China Under Grant No. 2007CB714200National Natural Science Foundation of China Under Grant No. 90715038
文摘In this paper, a method to develop a hierarchy of explicit recursion formulas for numerical simulation in an irregular grid for scalar wave equations is presented and its accuracy is illustrated via 2-D and 1-D models. Approaches to develop the stable formulas which are of 2M-order accuracy in both time and space with Mbeing a positive integer for regular grids are discussed and illustrated by constructing the second order (M= 1) and the fourth order (M = 2) recursion formulas.
基金Part of results of a project financially supported by the State Key Laboratory of Coastal and Offshore Engineering of Dalian University of Technology
文摘A three-step finite element method (FEM) together with Large Eddy Simulation (LES) is applied to incompressible turbulent flow around seabed pipelines at relatively high Reynolds numbers. Both two dimensional and three-dimensional numerical simulation is carried out to determine the three-dimensional effect. The results of numerical simulation agree quite well with the wave forces acting on pipeline models measured in physical model test.
基金supported by the National Natural Science Foundation of China (Grant No. 50679008)
文摘A two-dimensional numerical model based on the Navier-Stokes equations and computational Lagrangian-Eulerian advection remap-volume of fluid (CLEAR-VOF) method was developed to simulate wave and flow problems. The Navier-Stokes equations were discretized with a three-step finite element method that has a third-order accuracy. In the CLEAR-VOF method, the VOF function F was calculated in the Lagrangian manner and allowed the complicated free surface to be accurately captured. The propagation of regular waves and solitary waves over a flat bottom, and shoaling and breaking of solitary waves on two different slopes were simulated with this model, and the numerical results agreed with experimental data and theoretical solutions. A benchmark test of dam-collapse flow was also simulated with an unstructured mesh, and the capability of the present model for wave and flow simulations with unstructured meshes, was verified. The results show that the model is effective for numerical simulation of wave and flow problems with both structured and unstructured meshes.
基金financially supported by the National Natural Science Foundation of China(Grant Nos.52111530137 and 52025112)the Natural Science Found of Jiangsu Province(Grant No.BK20160556)the Jiangsu Provincial Higher Education Natural Science Research Major Project(Grant No.18KJA580003)。
文摘The cross-flow(CF)vortex-induced vibration(VIV)of a deepwater steep wave riser(SWR)subjected to uniform or shear flow loads is investigated numerically.The model is based on a three-dimensional(3D)nonlinear elastic rod theory coupled with a wake oscillator model.In this numerical simulation,the nonlinear motion equations of the riser with large deformation features are established in a global coordinate system to avoid the transformation between global and local coordinate systems,and are discretized with the time-domain finite element method(FEM).A wakeoscillator model is employed to study the vortex shedding,and the lift force generated by the wake flow is described in a van der Pol equation.A Newmark-βiterative scheme is used to solve their coupling equation for the VIV response of the SWR.The developed model is validated against the existing experimental results for the VIV response of the top-tension riser(TTR).Then,the numerical simulations are executed to determine VIV characteristics of the SWR.The effects of both flow velocity and the spanwise length of the flow field on the drag coefficient in the inline(IL)direction and the lift coefficient in the CF direction are investigated systematically.The results indicate that compared with TTR,the low frequency and multi-modal vibration are the main components of the SWR due to the large deformation and flexible characteristics.For shear flow,the multi-frequency resonance dominates the VIV response of the SWR,especially at the hang-off segment.
基金Project(51265044)supported by the National Natural Science Foundation of ChinaProject(2013TT2028)supported by the Science and Technology Project of Hunan Province of ChinaProject(2012QK162)supported by the Science and Technology Project of General Administration of Quality Supervision,Inspection and Quarantine of China
文摘In order to analyze the possibility of detecting defects in bend pipe using low-frequency ultrasonic guided wave, the propagation of T(0,1) mode and L(0,2) mode through straight-curved-straight pipe sections was studied. FE(finite element) models of bend pipe without defects and those with defects were introduced to analyze energy distribution, mode transition and defect detection of ultrasonic guided wave. FE simulation results were validated by experiments of four different bend pipes with circumferential defects in different positions. It is shown that most energy of T(0,1) mode or L(0,2) mode focuses on extrados of bend but little passes through intrados of bend, and T(0,1) mode or L(0,2) mode is converted to other possible non-axisymmetric modes when propagating through the bend and the defect after bend respectively. Furthermore, L(0,2) mode is more sensitive to circumferential notch than T(0,1) mode. The results of this work are beneficial for practical testing of pipes.
文摘It is well-known that artificial boundary conditions are crucial for the efficient and accurate computations of wavefields on unbounded domains. In this paper, we investigate stability analysis for the wave equation coupled with the first and the second order absorbing boundary conditions. The computational scheme is also developed. The approach allows the absorbing boundary conditions to be naturally imposed, which makes it easier for us to construct high order schemes for the absorbing boundary conditions. A thirdorder Lagrange finite element method with mass lumping is applied to obtain the spatial discretization of the wave equation. The resulting scheme is stable and is very efficient since no matrix inversion is needed at each time step. Moreover, we have shown both abstract and explicit conditional stability results for the fully-discrete schemes. The results are helpful for designing computational parameters in computations. Numerical computations are illustrated to show the efficiency and accuracy of our method. In particular, essentially no boundary reflection is seen at the artificial boundaries.
基金National Natural Science Foundation of China (50178065), 973 Program (2002CB412706), National Social Com-monweal Research Foundation (2002DIB30076) and Joint Seismological Science Foundation (101066).
文摘If a traditional explicit numerical integration algorithm is used to solve motion equation in the finite element simulation of wave motion, the time-step used by numerical integration is the smallest time-step restricted by the stability criterion in computational region. However, the excessively small time-step is usually unnecessary for a large portion of computational region. In this paper, a varying time-step explicit numerical integration algorithm is introduced, and its basic idea is to use different time-step restricted by the stability criterion in different computational region. Finally, the feasibility of the algorithm and its effect on calculating precision are verified by numerical test.
文摘目的:构建椎管重建内固定术的有限元模型,并分析椎管重建内固定术对脊柱稳定性的影响,验证椎管重建内固定术在椎管内手术中的有效性和可靠性。方法:筛选1名30岁男性,身高172 cm,体重75 kg的健康志愿者并采集其腰椎CT资料,建立正常腰椎L3-L5的有限元模型,并与体外实体结果和已发表的有限元分析结果进行比较,以验证该模型的有效性。根据处理方式不同分为正常组、椎板切除组和椎管重建组。在相同边界固定和生理载荷条件下,实施前屈、后伸、左弯、右弯、左旋和右旋6种工况活动,分析6种工况活动下L3-L4和L4-L5观察节段的活动度(range of motion,ROM)和L3-L5整体最大ROM的变化情况。结果:构建的L3-L5有限元模型各节段ROM位移与体外实体结果和既往文献数据结果吻合,确认了该模型的有效性。在L3-L4中,椎管重建组仅在后伸时较正常组增加,ROM变化百分比>5%;其他工况下的ROM与正常组相近,变化百分比<5%;而椎板切除组在前屈、后伸、左旋和右旋时的ROM较正常组和椎管重建组增加,变化百分比>5%。在L4-L5中,椎管重建组,在各工况下的ROM与正常组相近,变化百分比<5%,而椎板切除组在6种工况下的ROM均大于正常组和椎管重建组,变化百分比>5%。在L3-L5的整体最大ROM中,椎管重建组仅在后伸时超过正常组,变化百分比>5%;而椎板切除在前屈、后伸、左旋和右旋时大于正常组和椎管重建组,变化百分比>5%。L3-L5各节段ROM及总体ROM的变化趋势为:椎板切除组>椎管重建组>正常组。结论:椎板切除会严重影响脊柱生物力学的稳定性,而应用椎管重建内固定的方式可有效减少脊柱责任节段的ROM位移并维持其生物力学稳定性。