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.展开更多
This study presents an effective hybrid simulation approach for simulating broadband ground motion in complex near-fault locations.The approach utilizes a deterministic approach based on the spectral element method(SE...This study presents an effective hybrid simulation approach for simulating broadband ground motion in complex near-fault locations.The approach utilizes a deterministic approach based on the spectral element method(SEM),which is used to simulate low-frequency ground motion(f<1 Hz)by incorporating an innovative efficient discontinuous Galerkin(DG)method for grid division to accurately model basin sedimentary layers at reduced costs.It also introduces a comprehensive hybrid source model for high-frequency random scattering and a nonlinear analysis module for basin sedimentary layers.Deterministic outcomes are combined with modified three-dimensional stochastic finite fault method(3D-EXSIM)simulations of high-frequency ground motion(f>1 Hz).A fourth-order Butterworth filter with zero phase shift is employed for time-domain filtering of low-and high-frequency time series at a crossover frequency of 1 Hz,merging the low and high-frequency ground motions into a broadband time series.Taking an Ms 6.8 Luding earthquake,as an example,this hybrid method was used for a rapid and efficient simulation analysis of broadband ground motion in the region.The accuracy and efficiency of this hybrid method were verified through comparisons with actually observed station data and empirical attenuation curves.Deterministic method simulation results revealed the effects of mountainous topography,basin effects,nonlinear effects within the basin’s sedimentary layers,and a coupling interaction between the basin and the mountains.The findings are consistent with similar studies,showing that near-fault sedimentary basins significantly focus and amplify strong ground motion,and the soil’s nonlinear behavior in the basin influences ground motion to varying extents at different distances from the fault.The mountainous topography impacts the basin’s response to ground motion,leading to barrier effects.This research provides a scientific foundation for seismic zoning,urban planning,and seismic design in nearfault mountain basin regions.展开更多
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.展开更多
With the acceleration of the global aging process and the increase of cardiovascular ancerebrovascular diseases,more and more patients are paralyzed due to accidents,so theexoskeleton robot began to appear in people&#...With the acceleration of the global aging process and the increase of cardiovascular ancerebrovascular diseases,more and more patients are paralyzed due to accidents,so theexoskeleton robot began to appear in people's sight,and the lower limb exoskeleton robot withrehabilitation training is also favored by more and more people.In this paper,the structural designand analysis of the lower limb exoskeleton robot are carried out in view of the patients'expectation ofnormal walking.First,gait analysis and structural design of lower limb exoskeleton robot.Based onthe analysis of the walking gait of normal people,the freedom of the three key joints of the lower limbexoskeleton robot hip joint,knee joint and ankle joint is determined.at the same time,according tothe structuralcharacteristics of each joint,the three key joints are modeled respectively,and theoverall model assembly of the lower limb exoskeleton robot is completed.Secondly,the kinematicsanalysis of the lower limb exoskeleton robot was carried out to obtain the relationship between thelinear displacement,linear speed and acceleration of each joint,so as to ensure the coordination ofthe model with the human lower limb movement.Thirdly,the static analysis of typical gait of hipjoint,knee joint and ankle joint is carried out to verify the safety of the design model under thepremise of ensuring the structural strength requirements.Finally,the parts of the model were 3Dprinted,and the rationality of the design was further verified in the process of assembling the model.展开更多
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.展开更多
In this paper, using the theory of stochastic analysis of the response to earthquake load, a stochastic analysis method of the response of piled platforms to earthquake load has been established. In the method, the st...In this paper, using the theory of stochastic analysis of the response to earthquake load, a stochastic analysis method of the response of piled platforms to earthquake load has been established. In the method, the strong ground motion is considered as three dimensional stationary white noise process and the pile-soil interaction and water-structure interaction are considered. The stochastic response of a typical platform to earthquake load has been computed with this method and the results compared with those obtained with the response spectrum analysis method. The comparison shows that the stochastic analysis method of the response of piled platforms to earthquake load is suitable for this kind of analysis.展开更多
文摘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.
基金National Natural Science Foundation of China under Grant Nos.U2139208 and 52278516Key Laboratory of Earthquake Engineering and Engineering Vibration,China Earthquake Administration under Grant No.2024D15Key Laboratory of Soft Soil Characteristic and Engineering Environment,Tianjin Chengjian University under Grant No.2022SCEEKL003。
文摘This study presents an effective hybrid simulation approach for simulating broadband ground motion in complex near-fault locations.The approach utilizes a deterministic approach based on the spectral element method(SEM),which is used to simulate low-frequency ground motion(f<1 Hz)by incorporating an innovative efficient discontinuous Galerkin(DG)method for grid division to accurately model basin sedimentary layers at reduced costs.It also introduces a comprehensive hybrid source model for high-frequency random scattering and a nonlinear analysis module for basin sedimentary layers.Deterministic outcomes are combined with modified three-dimensional stochastic finite fault method(3D-EXSIM)simulations of high-frequency ground motion(f>1 Hz).A fourth-order Butterworth filter with zero phase shift is employed for time-domain filtering of low-and high-frequency time series at a crossover frequency of 1 Hz,merging the low and high-frequency ground motions into a broadband time series.Taking an Ms 6.8 Luding earthquake,as an example,this hybrid method was used for a rapid and efficient simulation analysis of broadband ground motion in the region.The accuracy and efficiency of this hybrid method were verified through comparisons with actually observed station data and empirical attenuation curves.Deterministic method simulation results revealed the effects of mountainous topography,basin effects,nonlinear effects within the basin’s sedimentary layers,and a coupling interaction between the basin and the mountains.The findings are consistent with similar studies,showing that near-fault sedimentary basins significantly focus and amplify strong ground motion,and the soil’s nonlinear behavior in the basin influences ground motion to varying extents at different distances from the fault.The mountainous topography impacts the basin’s response to ground motion,leading to barrier effects.This research provides a scientific foundation for seismic zoning,urban planning,and seismic design in nearfault mountain basin regions.
文摘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.
基金College Student Innovation andEntrepreneurship Project(Grant No.:S202414435026)ingkou Institute of Technology campus level research project——Development of food additive supercritical extraction equipment and fluid transmission systemresearch(Grant No.HX202427).
文摘With the acceleration of the global aging process and the increase of cardiovascular ancerebrovascular diseases,more and more patients are paralyzed due to accidents,so theexoskeleton robot began to appear in people's sight,and the lower limb exoskeleton robot withrehabilitation training is also favored by more and more people.In this paper,the structural designand analysis of the lower limb exoskeleton robot are carried out in view of the patients'expectation ofnormal walking.First,gait analysis and structural design of lower limb exoskeleton robot.Based onthe analysis of the walking gait of normal people,the freedom of the three key joints of the lower limbexoskeleton robot hip joint,knee joint and ankle joint is determined.at the same time,according tothe structuralcharacteristics of each joint,the three key joints are modeled respectively,and theoverall model assembly of the lower limb exoskeleton robot is completed.Secondly,the kinematicsanalysis of the lower limb exoskeleton robot was carried out to obtain the relationship between thelinear displacement,linear speed and acceleration of each joint,so as to ensure the coordination ofthe model with the human lower limb movement.Thirdly,the static analysis of typical gait of hipjoint,knee joint and ankle joint is carried out to verify the safety of the design model under thepremise of ensuring the structural strength requirements.Finally,the parts of the model were 3Dprinted,and the rationality of the design was further verified in the process of assembling the model.
文摘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.
文摘In this paper, using the theory of stochastic analysis of the response to earthquake load, a stochastic analysis method of the response of piled platforms to earthquake load has been established. In the method, the strong ground motion is considered as three dimensional stationary white noise process and the pile-soil interaction and water-structure interaction are considered. The stochastic response of a typical platform to earthquake load has been computed with this method and the results compared with those obtained with the response spectrum analysis method. The comparison shows that the stochastic analysis method of the response of piled platforms to earthquake load is suitable for this kind of analysis.
