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 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.展开更多
Dynamic contact theory is applied to simulate the sliding of surface fault. Finite element method is used to analyze the effect of surface fault to site ground motions. Calculated results indicate that amplification e...Dynamic contact theory is applied to simulate the sliding of surface fault. Finite element method is used to analyze the effect of surface fault to site ground motions. Calculated results indicate that amplification effect is obvious in the area near surface fault, especially on the site that is in the downside fault. The results show that the effect of surface fault should be considered when important structure is constructed in the site with surface fault.展开更多
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.展开更多
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.展开更多
The aim of this paper is to present a finite element modeling of the dynamic motion of a turbine rotor and its controller design with the mass unbalance under a crack on a rotating shaft. This process is an advanced m...The aim of this paper is to present a finite element modeling of the dynamic motion of a turbine rotor and its controller design with the mass unbalance under a crack on a rotating shaft. This process is an advanced method to the mathematical description of a system including an influence of a mass unbalance and a crack on the rotor shaft. As the first step, the shaft is physically modeled with a finite element method and the dynamic mathematical model is derived by using the Hamilton principle;thus, the system is represented by various subsystems. The equation of motion of a shaft with a mass unbalance and a crack is established by adapting the local mass unbalance and stiffness change through breathing and gaping from the existence of a crack. This is a reference system for the given system. Based on a fictitious model for transient behavior induced from vibration phenomena measured at the bearings, an elementary estimator is designed for the safety control and detection of a mass unbalance on the shaft. Using the state estimator, a bank of an estimator is established to get the diagnosis and the system data for a controller.展开更多
In this paper, we will use the explicit finite element to compute ground motion due to Tangshan earthquake. The explicit finite-element method uses one integration point and an hourglass control scheme. We implement t...In this paper, we will use the explicit finite element to compute ground motion due to Tangshan earthquake. The explicit finite-element method uses one integration point and an hourglass control scheme. We implement the coarse-grain method in a structured finite-element mesh straightforwardly. At the same time, we also apply the coarse-grain method to a widely used, slightly unstructured finite-element mesh, where unstructured finite elements are only used in the vertical velocity transition zones. By the finite-element methods, we can compute the ground velocity with some distance to the seismogenic fault in Tangshan earthquake. Through the computation, we can find the main character of ground motion for the strike slip earthquake events and the high frequency vibration motion of ground motion.展开更多
In the present study, finite element dynamic analysis or time history analysis of two-span beams subjected to asynchronous multi-support motions is carried out by using the moving support finite element. The elemental...In the present study, finite element dynamic analysis or time history analysis of two-span beams subjected to asynchronous multi-support motions is carried out by using the moving support finite element. The elemental equation of the element is based on total displacements and is derived under the concept of the quasi-static displacement decomposition. The use of moving support element shows that the element is very simple and convenient to represent continuous beam moving, deforming and vibrating simultaneously due to support motions. The comparison between the numerical results and analytical solutions indicates that the FE result agrees with the analytical solution.展开更多
General purpose computing on GPU for scientific computing has been rapidly growing in recent years. We investigate the applicability of GPU to discrete element method (DEM) often used in particle motion simulation. NV...General purpose computing on GPU for scientific computing has been rapidly growing in recent years. We investigate the applicability of GPU to discrete element method (DEM) often used in particle motion simulation. NVIDIA provides a sample code for this type of simulation, which obtained superior performance than CPU in computational time. A computational model of the contact force in NVIDIA’s sample code is, however, too simple to use in practice. This paper modifies the NVIDIA’s simple model by replacing it with the practical model. The computing speed of the practical model on GPU is compared with the simple one on GPU and with the practical one on CPU in numerical experiments. The result shows that the practical model on GPU obtains the computing speed 6 times faster than the practical one on CPU while 7 times slower than that of the simple one on GPU. The effects of the GPU architectures on the computing speed are analyzed.展开更多
As the continuation study on amplification of in-plane seismic ground motion by underground group cavities in layered half-space, this study extends to the case of poroelastic half-space with dry poroelastic and satur...As the continuation study on amplification of in-plane seismic ground motion by underground group cavities in layered half-space, this study extends to the case of poroelastic half-space with dry poroelastic and saturated poroelastic soil layers. The influence of poroelastic layers on the amplification of seismic ground motion is studied both in frequency domain and time domain using indirect boundary element method (IBEM). It is shown that for the example of a saturated poroelastic site in Tianjin under the excitation of Taft wave and E1 Centro wave, the amplification of seismic ground motion in poroelastic case is slightly smaller than that in the elastic case, and the amplification of PGA (peak ground acceleration) and its PRS (peak response spectrum).. can be increased up to 38.8% and 64.6%; the predominant period of response spectra in poroelastic case becomes shorter to some extent compared with that in the elastic case. It is suggested that the effect of underground group cavities in poroelastic half-space on design seismic ground motion should be considered.展开更多
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 context,a new boundary element algorithm based on the time-convoluted traction kernels is employed to evaluate the spatially varying earthquake ground motions of the Pacoima dam in the USA subjected to SH,SV a...In this context,a new boundary element algorithm based on the time-convoluted traction kernels is employed to evaluate the spatially varying earthquake ground motions of the Pacoima dam in the USA subjected to SH,SV and P incident waves.An accurate three-dimensional (3D) model of the dam canyon is implemented into the computer code BEMSA to investigate the seismic response of the dam.The analyses are performed in time domain with a linearly elastic constitutive model for the medium.This modeling procedure has been validated by the results reported in the existing literature.According to the results of this study,the response of the dam to earthquake waves is generally influenced by predominant frequency of the incident motion,surface topography,relative distance of observation points,and type of the incident seismic wave.For the cases considered,the incident SV wave has led to the maximum amplification of incident motions,especially at the left side of the dam.The results indicate that the proposed procedure can be employed for accurate prediction of a dam response during an earthquake.展开更多
Amplification of in-plane seismic ground motion by underground group cavities in layered half-space is studied both in frequency domain and time domain by using indirect boundary element method (IBEM), and the effec...Amplification of in-plane seismic ground motion by underground group cavities in layered half-space is studied both in frequency domain and time domain by using indirect boundary element method (IBEM), and the effect of cavity interval and spectrum of incident waves on the amplification are studied by numerical examples. It is shown that there may be large interaction between cavities, and group cavities with certain intervals may have significant amplification to seismic ground motion. The amplification of PGA (peak ground acceleration) and its PRS (peak response spectrum) can be increased up to 45.2% and 84.4%, for an example site in Tianjin, under the excitation of Taft wave and E1 Centro wave; and group cavities may also affect the spectra of the seismic ground motion. It is suggested that the effect of underground group cavities on design seismic ground motion should be considered.展开更多
An integrated dynamic model of China's deep ocean mining system is developed and the fast simulation analysis of its longitudinal reciprocating motion operation processes is achieved. The seafloor tracked miner is bu...An integrated dynamic model of China's deep ocean mining system is developed and the fast simulation analysis of its longitudinal reciprocating motion operation processes is achieved. The seafloor tracked miner is built as a three-dimensional single-body model with six-degree-of-freedom. The track-terrain interaction is modeled by partitioning the track-terrain interface into a certain number of mesh elements with three mutually perpendicular forces, including the normal force, the longitudinal shear force and the lateral shear force, acting on the center point of each mesh element. The hydrodynamic force of the miner is considered and applied. By considering the operational safety and collection efficiency, two new mining paths for the miner on the seafloor are proposed, which can be simulated with the established single-body dynamic model of the miner. The pipeline subsystem is built as a three-dimensional multi-body discrete element model, which is divided into rigid elements linked by flexible connectors. The flexible connector without mass is represented by six spring-damper elements. The external hydrodynamic forces of the ocean current from the longitudinal and lateral directions are both considered and modeled based on the Morison formula and applied to the mass center of each corresponding discrete rigid element. The mining ship is simplified and represented by a general kinematic point, whose heave motion induced by the ocean waves and the longitudinal and lateral towing motions are considered and applied. By integrating the single-body dynamic model of the miner and the multi-body discrete element dynamic model of the pipeline, and defining the kinematic equations of the mining ship, the integrated dynamic model of the total deep ocean mining system is formed. The longitudinal reciprocating motion operation modes of the total mining system, which combine the active straight-line and turning motions of the miner and the ship, and the passive towed motions of the pipeline, are proposed and simulated with the developed 3D dynamic model. Some critical simulation results are obtained and analyzed, such as the motion trajectories of key subsystems, the velocities of the buoyancy modules and the interaction forces between subsystems, which in a way can provide important theoretical basis and useful technical reference for the practical deep ocean mining system analysis, operation and control.展开更多
Based on analysis of regularity of stacking coal,discrete element simultaneous simulation is adopted to predict the process of unloading coal,which is proved to be effcient in the prediction of ash content.The results...Based on analysis of regularity of stacking coal,discrete element simultaneous simulation is adopted to predict the process of unloading coal,which is proved to be effcient in the prediction of ash content.The results show that the altitude of new irregular coal is equal to the income coal volume divided by area of cabin.The distribution of infnitesimal flow velocity helps to induce the motion equation of infnitesimal element,which provides the mathematical model for computer simulation.Swarm,a computer programming language,is utilized in this study.Adaptive infnitesimal stacking algorithm helps settle the diffculties in attainment of infnitesimal elements.The result of simulation is similar to the actual situation,which can accurately predict the ash contents of current time and cumulative time.Coal movement in the cabin is a new project,the result of which can also be applied to other solid particles and the widespread of the result will be highly valued.展开更多
基金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.
基金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.
文摘Dynamic contact theory is applied to simulate the sliding of surface fault. Finite element method is used to analyze the effect of surface fault to site ground motions. Calculated results indicate that amplification effect is obvious in the area near surface fault, especially on the site that is in the downside fault. The results show that the effect of surface fault should be considered when important structure is constructed in the site with surface fault.
文摘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.
文摘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.
文摘The aim of this paper is to present a finite element modeling of the dynamic motion of a turbine rotor and its controller design with the mass unbalance under a crack on a rotating shaft. This process is an advanced method to the mathematical description of a system including an influence of a mass unbalance and a crack on the rotor shaft. As the first step, the shaft is physically modeled with a finite element method and the dynamic mathematical model is derived by using the Hamilton principle;thus, the system is represented by various subsystems. The equation of motion of a shaft with a mass unbalance and a crack is established by adapting the local mass unbalance and stiffness change through breathing and gaping from the existence of a crack. This is a reference system for the given system. Based on a fictitious model for transient behavior induced from vibration phenomena measured at the bearings, an elementary estimator is designed for the safety control and detection of a mass unbalance on the shaft. Using the state estimator, a bank of an estimator is established to get the diagnosis and the system data for a controller.
文摘In this paper, we will use the explicit finite element to compute ground motion due to Tangshan earthquake. The explicit finite-element method uses one integration point and an hourglass control scheme. We implement the coarse-grain method in a structured finite-element mesh straightforwardly. At the same time, we also apply the coarse-grain method to a widely used, slightly unstructured finite-element mesh, where unstructured finite elements are only used in the vertical velocity transition zones. By the finite-element methods, we can compute the ground velocity with some distance to the seismogenic fault in Tangshan earthquake. Through the computation, we can find the main character of ground motion for the strike slip earthquake events and the high frequency vibration motion of ground motion.
文摘In the present study, finite element dynamic analysis or time history analysis of two-span beams subjected to asynchronous multi-support motions is carried out by using the moving support finite element. The elemental equation of the element is based on total displacements and is derived under the concept of the quasi-static displacement decomposition. The use of moving support element shows that the element is very simple and convenient to represent continuous beam moving, deforming and vibrating simultaneously due to support motions. The comparison between the numerical results and analytical solutions indicates that the FE result agrees with the analytical solution.
文摘General purpose computing on GPU for scientific computing has been rapidly growing in recent years. We investigate the applicability of GPU to discrete element method (DEM) often used in particle motion simulation. NVIDIA provides a sample code for this type of simulation, which obtained superior performance than CPU in computational time. A computational model of the contact force in NVIDIA’s sample code is, however, too simple to use in practice. This paper modifies the NVIDIA’s simple model by replacing it with the practical model. The computing speed of the practical model on GPU is compared with the simple one on GPU and with the practical one on CPU in numerical experiments. The result shows that the practical model on GPU obtains the computing speed 6 times faster than the practical one on CPU while 7 times slower than that of the simple one on GPU. The effects of the GPU architectures on the computing speed are analyzed.
基金supported by National Natural Science Foundation of China under grant No. 50978183Key Project for Applied Basic Research of Tianjin Municipality under Grant No. 12JCZDJC29000
文摘As the continuation study on amplification of in-plane seismic ground motion by underground group cavities in layered half-space, this study extends to the case of poroelastic half-space with dry poroelastic and saturated poroelastic soil layers. The influence of poroelastic layers on the amplification of seismic ground motion is studied both in frequency domain and time domain using indirect boundary element method (IBEM). It is shown that for the example of a saturated poroelastic site in Tianjin under the excitation of Taft wave and E1 Centro wave, the amplification of seismic ground motion in poroelastic case is slightly smaller than that in the elastic case, and the amplification of PGA (peak ground acceleration) and its PRS (peak response spectrum).. can be increased up to 38.8% and 64.6%; the predominant period of response spectra in poroelastic case becomes shorter to some extent compared with that in the elastic case. It is suggested that the effect of underground group cavities in poroelastic half-space on design seismic ground motion should be considered.
基金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 context,a new boundary element algorithm based on the time-convoluted traction kernels is employed to evaluate the spatially varying earthquake ground motions of the Pacoima dam in the USA subjected to SH,SV and P incident waves.An accurate three-dimensional (3D) model of the dam canyon is implemented into the computer code BEMSA to investigate the seismic response of the dam.The analyses are performed in time domain with a linearly elastic constitutive model for the medium.This modeling procedure has been validated by the results reported in the existing literature.According to the results of this study,the response of the dam to earthquake waves is generally influenced by predominant frequency of the incident motion,surface topography,relative distance of observation points,and type of the incident seismic wave.For the cases considered,the incident SV wave has led to the maximum amplification of incident motions,especially at the left side of the dam.The results indicate that the proposed procedure can be employed for accurate prediction of a dam response during an earthquake.
基金supported by National Natural Science Foundation of China under grant No. 50978183Tianjin Key Project for Applied Basic Research under grant No. 12JCZDJC29000
文摘Amplification of in-plane seismic ground motion by underground group cavities in layered half-space is studied both in frequency domain and time domain by using indirect boundary element method (IBEM), and the effect of cavity interval and spectrum of incident waves on the amplification are studied by numerical examples. It is shown that there may be large interaction between cavities, and group cavities with certain intervals may have significant amplification to seismic ground motion. The amplification of PGA (peak ground acceleration) and its PRS (peak response spectrum) can be increased up to 45.2% and 84.4%, for an example site in Tianjin, under the excitation of Taft wave and E1 Centro wave; and group cavities may also affect the spectra of the seismic ground motion. It is suggested that the effect of underground group cavities on design seismic ground motion should be considered.
基金supported by the National Natural Science Foundation of China(Grant No.51105386)the National Deep-Sea Technology Project of Development and Research(Grant No.DYXM-115-04-02-01)the Fundamental Research Funds for the Central Universities(Grant No.2011QNZT058)
文摘An integrated dynamic model of China's deep ocean mining system is developed and the fast simulation analysis of its longitudinal reciprocating motion operation processes is achieved. The seafloor tracked miner is built as a three-dimensional single-body model with six-degree-of-freedom. The track-terrain interaction is modeled by partitioning the track-terrain interface into a certain number of mesh elements with three mutually perpendicular forces, including the normal force, the longitudinal shear force and the lateral shear force, acting on the center point of each mesh element. The hydrodynamic force of the miner is considered and applied. By considering the operational safety and collection efficiency, two new mining paths for the miner on the seafloor are proposed, which can be simulated with the established single-body dynamic model of the miner. The pipeline subsystem is built as a three-dimensional multi-body discrete element model, which is divided into rigid elements linked by flexible connectors. The flexible connector without mass is represented by six spring-damper elements. The external hydrodynamic forces of the ocean current from the longitudinal and lateral directions are both considered and modeled based on the Morison formula and applied to the mass center of each corresponding discrete rigid element. The mining ship is simplified and represented by a general kinematic point, whose heave motion induced by the ocean waves and the longitudinal and lateral towing motions are considered and applied. By integrating the single-body dynamic model of the miner and the multi-body discrete element dynamic model of the pipeline, and defining the kinematic equations of the mining ship, the integrated dynamic model of the total deep ocean mining system is formed. The longitudinal reciprocating motion operation modes of the total mining system, which combine the active straight-line and turning motions of the miner and the ship, and the passive towed motions of the pipeline, are proposed and simulated with the developed 3D dynamic model. Some critical simulation results are obtained and analyzed, such as the motion trajectories of key subsystems, the velocities of the buoyancy modules and the interaction forces between subsystems, which in a way can provide important theoretical basis and useful technical reference for the practical deep ocean mining system analysis, operation and control.
基金the financial support provided by the National Natural Science Foundation of China(No.51174202)Jiangsu Natural Science Foundation of China(No.20100095110013)
文摘Based on analysis of regularity of stacking coal,discrete element simultaneous simulation is adopted to predict the process of unloading coal,which is proved to be effcient in the prediction of ash content.The results show that the altitude of new irregular coal is equal to the income coal volume divided by area of cabin.The distribution of infnitesimal flow velocity helps to induce the motion equation of infnitesimal element,which provides the mathematical model for computer simulation.Swarm,a computer programming language,is utilized in this study.Adaptive infnitesimal stacking algorithm helps settle the diffculties in attainment of infnitesimal elements.The result of simulation is similar to the actual situation,which can accurately predict the ash contents of current time and cumulative time.Coal movement in the cabin is a new project,the result of which can also be applied to other solid particles and the widespread of the result will be highly valued.