The 2.5D finite/infinite element approach is adopted to study wave propagation problems caused by underground moving trains. The irregularities of the near field, including the tunnel structure and parts of the soil, ...The 2.5D finite/infinite element approach is adopted to study wave propagation problems caused by underground moving trains. The irregularities of the near field, including the tunnel structure and parts of the soil, are modeled by the finite elements, and the wave propagation properties of the far field extending to infinity are modeled by the infinite elements. One particular feature of the 2.5D approach is that it enables the computation of the three-dimensional response of the half-space, taking into account the load-moving effect, using only a two-dimensional profile. Although the 2.5D finite/infinite element approach shows a great advantage in studying the wave propagation caused by moving trains, attention should be given to the calculation aspects, such as the rules for mesh establishment, in order to avoid producing inaccurate or erroneous results. In this paper, some essential points for consideration in analysis are highlighted, along with techniques to enhance the speed of the calculations. All these observations should prove useful in making the 2.5D finite/infinite element approach an effective one.展开更多
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 application of the finite layer & triangular prism element method to the 3D ground subsidence and stress analysis caused by mining is presented. The layer elements and the triangular prism elements have been a...The application of the finite layer & triangular prism element method to the 3D ground subsidence and stress analysis caused by mining is presented. The layer elements and the triangular prism elements have been alternatively used in the numerical simulation system, the displacement pattern, strain matrix, elastic matrix, stiffness matrix, load matrix and the stress matrix of the layer element and triangular prism element have been presented. By means of the Fortran90 programming language, a numerical simulation system based on finite layer & triangular prism element have been built up, and this system is suitable for subsidence prediction and stress analysis of all mining condition and mining methods. Comparing with the infinite element method, this approach dramatically reduces the size of the set of equations that need to be solved, and greatly reduces the amount of data preparation required. It not only saves the internal storage, and the computation time, but also decreases the cost.展开更多
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.展开更多
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.展开更多
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.展开更多
Ground water samples are collected from south West Bank/Palestine and analyzed for different rare elements (Rb, Zr, U, P, Ti, V), rare earth elements (La, Ce, and Nd), and other common trace metals (Li, Na, Mg, Ca, Sr...Ground water samples are collected from south West Bank/Palestine and analyzed for different rare elements (Rb, Zr, U, P, Ti, V), rare earth elements (La, Ce, and Nd), and other common trace metals (Li, Na, Mg, Ca, Sr, Ba, K, Bi) that most of them usually have no maximum acceptable limits as either they are considered not to be toxic to human health or there is no sufficient data about their toxicity to human health. This study was conducted to determine the water quality of ground water which is used for drinking in the study area. Water samples from ten groundwater wells were obtained in three different dates of the year (November 2012, March 2013, and April 2013). Three water samples were obtained from each well for each sampling date;so a total of 90 water samples were collected from the ten wells. The results obtained from this study suggest a possible risk to the population of the study area given the high concentration of some metals that have no maximum allowed concentration, and the fact that for many people in the study area, ground water is a main source of their water supply.展开更多
An element-free Galerkin method(EFGM) is used to solve the two-dimensional(2D) ground penetrating radar(GPR)modelling problems, due to its simple pre-processing, the absence of elements and high accuracy. Different fr...An element-free Galerkin method(EFGM) is used to solve the two-dimensional(2D) ground penetrating radar(GPR)modelling problems, due to its simple pre-processing, the absence of elements and high accuracy. Different from element-based numerical methods, this approach makes nodes free from the elemental restraint and avoids the explicit mesh discretization. First, we derived the boundary value problem for the 2D GPR simulation problems. Second, a penalty function approach and a boundary condition truncated method were used to enforce the essential and the absorbing boundary conditions, respectively. A three-layered GPR model was used to verify our element-free approach. The numerical solutions show that our solutions have an excellent agreement with solutions of a finite element method(FEM). Then, we used the EFGM to simulate one more complex model to show its capability and limitations. Simulation results show that one obvious advantage of EFGM is the absence of element mesh, which makes the method very flexible. Due to the use of MLS fitting, a key feature of EFM, is that both the dependent variable and its gradient are continuous and have high precision.展开更多
As a highly efficient absorbing boundary condition, Perfectly Matched Layer (PML) has been widely used in Finite Difference Time Domain (FDTD) simulation of Ground Penetrating Radar (GPR) based on the first order elec...As a highly efficient absorbing boundary condition, Perfectly Matched Layer (PML) has been widely used in Finite Difference Time Domain (FDTD) simulation of Ground Penetrating Radar (GPR) based on the first order electromagnetic wave equation. However, the PML boundary condition is difficult to apply in GPR Finite Element Time Domain (FETD) simulation based on the second order electromagnetic wave equation. This paper developed a non-split perfectly matched layer (NPML) boundary condition for GPR FETD simulation based on the second order electromagnetic wave equation. Taking two-dimensional TM wave equation as an example, the second order frequency domain equation of GPR was derived according to the definition of complex extending coordinate transformation. Then it transformed into time domain by means of auxiliary differential equation method, and its FETD equation is derived based on Galerkin method. On this basis, a GPR FETD forward program based on NPML boundary condition is developed. The merits of NPML boundary condition are certified by compared with wave field snapshots, signal and reflection errors of homogeneous medium model with split and non-split PML boundary conditions. The comparison demonstrated that the NPML algorithm can reduce memory occupation and improve calculation efficiency. Furthermore, numerical simulation of a complex model verifies the good absorption effects of the NPML boundary condition in complex structures.展开更多
基金Science Council Under Grant No.NSC 89-2211-E-002-020
文摘The 2.5D finite/infinite element approach is adopted to study wave propagation problems caused by underground moving trains. The irregularities of the near field, including the tunnel structure and parts of the soil, are modeled by the finite elements, and the wave propagation properties of the far field extending to infinity are modeled by the infinite elements. One particular feature of the 2.5D approach is that it enables the computation of the three-dimensional response of the half-space, taking into account the load-moving effect, using only a two-dimensional profile. Although the 2.5D finite/infinite element approach shows a great advantage in studying the wave propagation caused by moving trains, attention should be given to the calculation aspects, such as the rules for mesh establishment, in order to avoid producing inaccurate or erroneous results. In this paper, some essential points for consideration in analysis are highlighted, along with techniques to enhance the speed of the calculations. All these observations should prove useful in making the 2.5D finite/infinite element approach an effective one.
基金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 application of the finite layer & triangular prism element method to the 3D ground subsidence and stress analysis caused by mining is presented. The layer elements and the triangular prism elements have been alternatively used in the numerical simulation system, the displacement pattern, strain matrix, elastic matrix, stiffness matrix, load matrix and the stress matrix of the layer element and triangular prism element have been presented. By means of the Fortran90 programming language, a numerical simulation system based on finite layer & triangular prism element have been built up, and this system is suitable for subsidence prediction and stress analysis of all mining condition and mining methods. Comparing with the infinite element method, this approach dramatically reduces the size of the set of equations that need to be solved, and greatly reduces the amount of data preparation required. It not only saves the internal storage, and the computation time, but also decreases the cost.
文摘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.
文摘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.
文摘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.
文摘Ground water samples are collected from south West Bank/Palestine and analyzed for different rare elements (Rb, Zr, U, P, Ti, V), rare earth elements (La, Ce, and Nd), and other common trace metals (Li, Na, Mg, Ca, Sr, Ba, K, Bi) that most of them usually have no maximum acceptable limits as either they are considered not to be toxic to human health or there is no sufficient data about their toxicity to human health. This study was conducted to determine the water quality of ground water which is used for drinking in the study area. Water samples from ten groundwater wells were obtained in three different dates of the year (November 2012, March 2013, and April 2013). Three water samples were obtained from each well for each sampling date;so a total of 90 water samples were collected from the ten wells. The results obtained from this study suggest a possible risk to the population of the study area given the high concentration of some metals that have no maximum allowed concentration, and the fact that for many people in the study area, ground water is a main source of their water supply.
基金Project(41074085)supported by the National Natural Science Foundation of ChinaProject(NCET-12-0551)supported by the Funds for New Century Excellent Talents in University,ChinaProject supported by Shenghua Yuying Program of Central South University,China
文摘An element-free Galerkin method(EFGM) is used to solve the two-dimensional(2D) ground penetrating radar(GPR)modelling problems, due to its simple pre-processing, the absence of elements and high accuracy. Different from element-based numerical methods, this approach makes nodes free from the elemental restraint and avoids the explicit mesh discretization. First, we derived the boundary value problem for the 2D GPR simulation problems. Second, a penalty function approach and a boundary condition truncated method were used to enforce the essential and the absorbing boundary conditions, respectively. A three-layered GPR model was used to verify our element-free approach. The numerical solutions show that our solutions have an excellent agreement with solutions of a finite element method(FEM). Then, we used the EFGM to simulate one more complex model to show its capability and limitations. Simulation results show that one obvious advantage of EFGM is the absence of element mesh, which makes the method very flexible. Due to the use of MLS fitting, a key feature of EFM, is that both the dependent variable and its gradient are continuous and have high precision.
文摘As a highly efficient absorbing boundary condition, Perfectly Matched Layer (PML) has been widely used in Finite Difference Time Domain (FDTD) simulation of Ground Penetrating Radar (GPR) based on the first order electromagnetic wave equation. However, the PML boundary condition is difficult to apply in GPR Finite Element Time Domain (FETD) simulation based on the second order electromagnetic wave equation. This paper developed a non-split perfectly matched layer (NPML) boundary condition for GPR FETD simulation based on the second order electromagnetic wave equation. Taking two-dimensional TM wave equation as an example, the second order frequency domain equation of GPR was derived according to the definition of complex extending coordinate transformation. Then it transformed into time domain by means of auxiliary differential equation method, and its FETD equation is derived based on Galerkin method. On this basis, a GPR FETD forward program based on NPML boundary condition is developed. The merits of NPML boundary condition are certified by compared with wave field snapshots, signal and reflection errors of homogeneous medium model with split and non-split PML boundary conditions. The comparison demonstrated that the NPML algorithm can reduce memory occupation and improve calculation efficiency. Furthermore, numerical simulation of a complex model verifies the good absorption effects of the NPML boundary condition in complex structures.