An explicit finite element-finite difference method for analyzing the effects of two-dimensional visco-elastic localtopography on earthquake ground motion is prOPosed in this paper. In the method, at first, the finite...An explicit finite element-finite difference method for analyzing the effects of two-dimensional visco-elastic localtopography on earthquake ground motion is prOPosed in this paper. In the method, at first, the finite elementdiscrete model is formed by using the artificial boundary and finite element method, and the dynamic equationsof local nodes in the discrete model are obtained according to the theory of the special finite element method similar to the finite difference method, and then the explicit step-by-step integration formulas are presented by usingthe explicit difference method for solving the visco-elastic dynamic equation and Generalized Multi-transmittingBoundary. The method has the advantages of saving computing time and computer memory space, and it is suitable for any case of topography and has high computing accuracy and good computing stability.展开更多
In this paper, we consider two methods, the Second order Central Difference Method (SCDM) and the Finite Element Method (FEM) with P1 triangular elements, for solving two dimensional general linear Elliptic Partial Di...In this paper, we consider two methods, the Second order Central Difference Method (SCDM) and the Finite Element Method (FEM) with P1 triangular elements, for solving two dimensional general linear Elliptic Partial Differential Equations (PDE) with mixed derivatives along with Dirichlet and Neumann boundary conditions. These two methods have almost the same accuracy from theoretical aspect with regular boundaries, but generally Finite Element Method produces better approximations when the boundaries are irregular. In order to investigate which method produces better results from numerical aspect, we apply these methods into specific examples with regular boundaries with constant step-size for both of them. The results which obtained confirm, in most of the cases, the theoretical results.展开更多
Chaotic vibrations of flexible non-linear Euler-Bernoulli beams subjected to harmonic load and with various boundary conditions（symmetric and non-symmetric）are studied in this work.Reliability of the obtained result...Chaotic vibrations of flexible non-linear Euler-Bernoulli beams subjected to harmonic load and with various boundary conditions（symmetric and non-symmetric）are studied in this work.Reliability of the obtained results is verified by the finite difference method（FDM）and the finite element method（FEM）with the Bubnov-Galerkin approximation for various boundary conditions and various dynamic regimes（regular and non-regular）.The influence of boundary conditions on the Euler-Bernoulli beams dynamics is studied mainly,dynamic behavior vs.control parameters { ωp,q0 } is reported,and scenarios of the system transition into chaos are illustrated.展开更多
The asymptotic convergence of the solution of the parabolic equation is proved. By the eigenvalues estimation, we obtain that the approximate solutions by the finite difference method and the finite element method are...The asymptotic convergence of the solution of the parabolic equation is proved. By the eigenvalues estimation, we obtain that the approximate solutions by the finite difference method and the finite element method are asymptotically convergent. Both methods are considered in continnous time.展开更多
The non_stationary natural convection problem is studied. A lowest order finite difference scheme based on mixed finite element method for non_stationary natural convection problem, by the spatial variations discreted...The non_stationary natural convection problem is studied. A lowest order finite difference scheme based on mixed finite element method for non_stationary natural convection problem, by the spatial variations discreted with finite element method and time with finite difference scheme was derived, where the numerical solution of velocity, pressure, and temperature can be found together, and a numerical example to simulate the close square cavity is given, which is of practical importance.展开更多
Finite difference type preconditioners for spectral element discretizations based on Legendre-Gauss-Lobatto points are analyzed. The latter is employed for the approximation of uniformly elliptic partial differential ...Finite difference type preconditioners for spectral element discretizations based on Legendre-Gauss-Lobatto points are analyzed. The latter is employed for the approximation of uniformly elliptic partial differential problems. In this work, it is shown that the condition number of the resulting preconditioned system is bounded independently of both of the polynomial degrees used in the spectral element method and the element sizes. Several numerical tests verify the h-p independence of the proposed preconditioning.展开更多
A nonconforming finite element method of finite difference streamline diffusion type is proposed to solve the time-dependent linearized Navier-Stokes equations. The backward Euler scheme is used for time discretizatio...A nonconforming finite element method of finite difference streamline diffusion type is proposed to solve the time-dependent linearized Navier-Stokes equations. The backward Euler scheme is used for time discretization. Crouzeix-Raviart nonconforming finite element approximation, namely, nonconforming （P1）2 - P0 element, is used for the velocity and pressure fields with the streamline diffusion technique to cope with usual instabilities caused by the convection and time terms. Stability and error estimates are derived with suitable norms.展开更多
We propose a novel symplectic finite element method to solve the structural dynamic responses of linear elastic systems.For the dynamic responses of continuous medium structures,the traditional numerical algorithm is ...We propose a novel symplectic finite element method to solve the structural dynamic responses of linear elastic systems.For the dynamic responses of continuous medium structures,the traditional numerical algorithm is the dissipative algorithm and cannot maintain long-term energy conservation.Thus,a symplectic finite element method with energy conservation is constructed in this paper.A linear elastic system can be discretized into multiple elements,and a Hamiltonian system of each element can be constructed.The single element is discretized by the Galerkin method,and then the Hamiltonian system is constructed into the Birkhoffian system.Finally,all the elements are combined to obtain the vibration equation of the continuous system and solved by the symplectic difference scheme.Through the numerical experiments of the vibration response of the Bernoulli-Euler beam and composite plate,it is found that the vibration response solution and energy obtained with the algorithm are superior to those of the Runge-Kutta algorithm.The results show that the symplectic finite element method can keep energy conservation for a long time and has higher stability in solving the dynamic responses of linear elastic systems.展开更多
The nonlinear finite element（FE） analysis has been widely used in the design and analysis of structural or geotechnical systems.The response sensitivities（or gradients） to the model parameters are of significant i...The nonlinear finite element（FE） analysis has been widely used in the design and analysis of structural or geotechnical systems.The response sensitivities（or gradients） to the model parameters are of significant importance in these realistic engineering problems.However the sensitivity calculation has lagged behind,leaving a gap between advanced FE response analysis and other research hotspots using the response gradient.The response sensitivity analysis is crucial for any gradient-based algorithms,such as reliability analysis,system identification and structural optimization.Among various sensitivity analysis methods,the direct differential method（DDM） has advantages of computing efficiency and accuracy,providing an ideal tool for the response gradient calculation.This paper extended the DDM framework to realistic complicated soil-foundation-structure interaction（SFSI） models by developing the response gradients for various constraints,element and materials involved.The enhanced framework is applied to three-dimensional SFSI system prototypes for a pilesupported bridge pier and a pile-supported reinforced concrete building frame structure,subjected to earthquake loading conditions.The DDM results are verified by forward finite difference method（FFD）.The relative importance（RI） of the various material parameters on the responses of SFSI system are investigated based on the DDM response sensitivity results.The FFD converges asymptotically toward the DDM results,demonstrating the advantages of DDM（e.g.,accurate,efficient,insensitive to numerical noise）.Furthermore,the RI and effects of the model parameters of structure,foundation and soil materials on the responses of SFSI systems are investigated by taking advantage of the sensitivity analysis results.The extension of DDM to SFSI systems greatly broaden the application areas of the d gradient-based algorithms,e.g.FE model updating and nonlinear system identification of complicated SFSI systems.展开更多
To solve the equation for gravity-gyroscopic waves in a rectangular domain, the distinguished algorithm for the solution of the Cauchy problem for a second-order transient equation is proposed. This algorithm is devel...To solve the equation for gravity-gyroscopic waves in a rectangular domain, the distinguished algorithm for the solution of the Cauchy problem for a second-order transient equation is proposed. This algorithm is developed by using the time-varying finite element method. The space derivatives in the gravity-gyroscopic wave equation are approximated with finite differences. The stability and accuracy of the method are estimated. The procedure for the implementation of the method is developed. The calculations were performed for determining the steady-state modes of fluctuations of the solutions of the gravity-gyroscopic wave equation depending on the problem parameters.展开更多
Based on the experimental results of the work-hardening processes of single crystals,the ac- commodation processes of polycrvstal deformation and the assumption of idealized polycrystal,the stress-strain relation of e...Based on the experimental results of the work-hardening processes of single crystals,the ac- commodation processes of polycrvstal deformation and the assumption of idealized polycrystal,the stress-strain relation of elasto-plastic deformation crystal was derived.The effect of orientation difference on the mechanical properties of the bicrvstalline specimens of aluminum was simulated by means of the finite element method(FEM)of elasto-plastic crystal.The results are in good agreement with the experimental results made by Clark and Chalmers in 1954.展开更多
Partial differential equations arise in formulations of problems involving functions of several variables such as the propagation of sound or heat, electrostatics, electrodynamics, fluid flow, and elasticity, etc. The...Partial differential equations arise in formulations of problems involving functions of several variables such as the propagation of sound or heat, electrostatics, electrodynamics, fluid flow, and elasticity, etc. The present paper deals with a general introduction and classification of partial differential equations and the numerical methods available in the literature for the solution of partial differential equations.展开更多
This paper presents the dimension split element-free Galerkin （DSEFG） method for three-dimensional potential problems, and the corresponding formulae are obtained. The main idea of the DSEFG method is that a three-d...This paper presents the dimension split element-free Galerkin （DSEFG） method for three-dimensional potential problems, and the corresponding formulae are obtained. The main idea of the DSEFG method is that a three-dimensional potential problem can be transformed into a series of two-dimensional problems. For these two-dimensional problems, the improved moving least-squares （IMLS） approximation is applied to construct the shape function, which uses an orthogonal function system with a weight function as the basis functions. The Galerkin weak form is applied to obtain a discretized system equation, and the penalty method is employed to impose the essential boundary condition. The finite difference method is selected in the splitting direction. For the purposes of demonstration, some selected numerical examples are solved using the DSEFG method. The convergence study and error analysis of the DSEFG method are presented. The numerical examples show that the DSEFG method has greater computational precision and computational efficiency than the IEFG method.展开更多
This article analyzes the relationship between the water level and the water-tube tilting in Shuangyang lake,based on the differential deformation features reflected by the NS and EW components of the water-tube tiltm...This article analyzes the relationship between the water level and the water-tube tilting in Shuangyang lake,based on the differential deformation features reflected by the NS and EW components of the water-tube tiltmeter.The results show a good spatiotemporal consistency between the variation of water level and the NS tilt component,which is considered to be affected by the magnitude and duration of the water level variation in Shuangyang Lake.The article uses Landsat remote sensing image data to extract the water boundary of Shuan-gyang Lake,and takes advantage of the finite element numerical simulation method to build three-dimensional models for different geological structural conditions of the Shuangyang seismostation.The simulation results show that when the underground medium is granite,the effect of water level variation on the vertical displacement of the surface is non-directional.With a 50-m soil layer in Model 2,the simulated NS tilt variation is equivalent to the actual observed water-tube tiltmeter NS component when the water level variation is 0.44 m and 0.8m.When the variation of water level reaches 2.0m,the simulation result of the NS component is 79.6 ms,which is slightly larger than the observed result of 60.32 ms.However,the simulation results show that the variation of the EW component is significantly smaller than that of the NS one.Due to the fact that the Shuan-gyang lake is long in the NS direction and short in the EW direction,the existence of the soil layer tends to generate ground deformation along the NS direction in the vicinity of the lake after the increase of water level,thereby resulting in the difference of the ground deformation in the two directions.展开更多
We present the finite difference/element method for a two-dimensional modified fractional diffusion equation.The analysis is carried out first for the time semi-discrete scheme,and then for the full discrete scheme.Th...We present the finite difference/element method for a two-dimensional modified fractional diffusion equation.The analysis is carried out first for the time semi-discrete scheme,and then for the full discrete scheme.The time discretization is based on the L1-approximation for the fractional derivative terms and the second-order backward differentiation formula for the classical first order derivative term.We use finite element method for the spatial approximation in full discrete scheme.We show that both the semi-discrete and full discrete schemes are unconditionally stable and convergent.Moreover,the optimal convergence rate is obtained.Finally,some numerical examples are tested in the case of one and two space dimensions and the numerical results confirm our theoretical analysis.展开更多
文摘An explicit finite element-finite difference method for analyzing the effects of two-dimensional visco-elastic localtopography on earthquake ground motion is prOPosed in this paper. In the method, at first, the finite elementdiscrete model is formed by using the artificial boundary and finite element method, and the dynamic equationsof local nodes in the discrete model are obtained according to the theory of the special finite element method similar to the finite difference method, and then the explicit step-by-step integration formulas are presented by usingthe explicit difference method for solving the visco-elastic dynamic equation and Generalized Multi-transmittingBoundary. The method has the advantages of saving computing time and computer memory space, and it is suitable for any case of topography and has high computing accuracy and good computing stability.
文摘In this paper, we consider two methods, the Second order Central Difference Method (SCDM) and the Finite Element Method (FEM) with P1 triangular elements, for solving two dimensional general linear Elliptic Partial Differential Equations (PDE) with mixed derivatives along with Dirichlet and Neumann boundary conditions. These two methods have almost the same accuracy from theoretical aspect with regular boundaries, but generally Finite Element Method produces better approximations when the boundaries are irregular. In order to investigate which method produces better results from numerical aspect, we apply these methods into specific examples with regular boundaries with constant step-size for both of them. The results which obtained confirm, in most of the cases, the theoretical results.
文摘Chaotic vibrations of flexible non-linear Euler-Bernoulli beams subjected to harmonic load and with various boundary conditions（symmetric and non-symmetric）are studied in this work.Reliability of the obtained results is verified by the finite difference method（FDM）and the finite element method（FEM）with the Bubnov-Galerkin approximation for various boundary conditions and various dynamic regimes（regular and non-regular）.The influence of boundary conditions on the Euler-Bernoulli beams dynamics is studied mainly,dynamic behavior vs.control parameters { ωp,q0 } is reported,and scenarios of the system transition into chaos are illustrated.
文摘The asymptotic convergence of the solution of the parabolic equation is proved. By the eigenvalues estimation, we obtain that the approximate solutions by the finite difference method and the finite element method are asymptotically convergent. Both methods are considered in continnous time.
文摘The non_stationary natural convection problem is studied. A lowest order finite difference scheme based on mixed finite element method for non_stationary natural convection problem, by the spatial variations discreted with finite element method and time with finite difference scheme was derived, where the numerical solution of velocity, pressure, and temperature can be found together, and a numerical example to simulate the close square cavity is given, which is of practical importance.
文摘Finite difference type preconditioners for spectral element discretizations based on Legendre-Gauss-Lobatto points are analyzed. The latter is employed for the approximation of uniformly elliptic partial differential problems. In this work, it is shown that the condition number of the resulting preconditioned system is bounded independently of both of the polynomial degrees used in the spectral element method and the element sizes. Several numerical tests verify the h-p independence of the proposed preconditioning.
基金supported by the National Natural Science Foundation of China(No.10771150)the National Basic Research Program of China(No.2005CB321701)+1 种基金the Program for New Century Excellent Talents in University(No.NCET-07-0584)the Natural Science Foundation of Sichuan Province(No.07ZB087)
文摘A nonconforming finite element method of finite difference streamline diffusion type is proposed to solve the time-dependent linearized Navier-Stokes equations. The backward Euler scheme is used for time discretization. Crouzeix-Raviart nonconforming finite element approximation, namely, nonconforming （P1）2 - P0 element, is used for the velocity and pressure fields with the streamline diffusion technique to cope with usual instabilities caused by the convection and time terms. Stability and error estimates are derived with suitable norms.
基金supported by the National Natural Science Foundation of China(Nos.12132001 and 52192632)。
文摘We propose a novel symplectic finite element method to solve the structural dynamic responses of linear elastic systems.For the dynamic responses of continuous medium structures,the traditional numerical algorithm is the dissipative algorithm and cannot maintain long-term energy conservation.Thus,a symplectic finite element method with energy conservation is constructed in this paper.A linear elastic system can be discretized into multiple elements,and a Hamiltonian system of each element can be constructed.The single element is discretized by the Galerkin method,and then the Hamiltonian system is constructed into the Birkhoffian system.Finally,all the elements are combined to obtain the vibration equation of the continuous system and solved by the symplectic difference scheme.Through the numerical experiments of the vibration response of the Bernoulli-Euler beam and composite plate,it is found that the vibration response solution and energy obtained with the algorithm are superior to those of the Runge-Kutta algorithm.The results show that the symplectic finite element method can keep energy conservation for a long time and has higher stability in solving the dynamic responses of linear elastic systems.
基金National Key Research and Development Program of China under Grant No.2016YFC0701106Natural Sciences and Engineering Research Council of Canada via Discovery under Grant No.NSERC RGPIN-2017-05556 Li
文摘The nonlinear finite element（FE） analysis has been widely used in the design and analysis of structural or geotechnical systems.The response sensitivities（or gradients） to the model parameters are of significant importance in these realistic engineering problems.However the sensitivity calculation has lagged behind,leaving a gap between advanced FE response analysis and other research hotspots using the response gradient.The response sensitivity analysis is crucial for any gradient-based algorithms,such as reliability analysis,system identification and structural optimization.Among various sensitivity analysis methods,the direct differential method（DDM） has advantages of computing efficiency and accuracy,providing an ideal tool for the response gradient calculation.This paper extended the DDM framework to realistic complicated soil-foundation-structure interaction（SFSI） models by developing the response gradients for various constraints,element and materials involved.The enhanced framework is applied to three-dimensional SFSI system prototypes for a pilesupported bridge pier and a pile-supported reinforced concrete building frame structure,subjected to earthquake loading conditions.The DDM results are verified by forward finite difference method（FFD）.The relative importance（RI） of the various material parameters on the responses of SFSI system are investigated based on the DDM response sensitivity results.The FFD converges asymptotically toward the DDM results,demonstrating the advantages of DDM（e.g.,accurate,efficient,insensitive to numerical noise）.Furthermore,the RI and effects of the model parameters of structure,foundation and soil materials on the responses of SFSI systems are investigated by taking advantage of the sensitivity analysis results.The extension of DDM to SFSI systems greatly broaden the application areas of the d gradient-based algorithms,e.g.FE model updating and nonlinear system identification of complicated SFSI systems.
文摘To solve the equation for gravity-gyroscopic waves in a rectangular domain, the distinguished algorithm for the solution of the Cauchy problem for a second-order transient equation is proposed. This algorithm is developed by using the time-varying finite element method. The space derivatives in the gravity-gyroscopic wave equation are approximated with finite differences. The stability and accuracy of the method are estimated. The procedure for the implementation of the method is developed. The calculations were performed for determining the steady-state modes of fluctuations of the solutions of the gravity-gyroscopic wave equation depending on the problem parameters.
文摘Based on the experimental results of the work-hardening processes of single crystals,the ac- commodation processes of polycrvstal deformation and the assumption of idealized polycrystal,the stress-strain relation of elasto-plastic deformation crystal was derived.The effect of orientation difference on the mechanical properties of the bicrvstalline specimens of aluminum was simulated by means of the finite element method(FEM)of elasto-plastic crystal.The results are in good agreement with the experimental results made by Clark and Chalmers in 1954.
文摘Partial differential equations arise in formulations of problems involving functions of several variables such as the propagation of sound or heat, electrostatics, electrodynamics, fluid flow, and elasticity, etc. The present paper deals with a general introduction and classification of partial differential equations and the numerical methods available in the literature for the solution of partial differential equations.
基金supported by the National Natural Science Foundation of China (Grants 11571223, 51404160)Shanxi Province Science Foundation for Youths (Grant 2014021025-1)
文摘This paper presents the dimension split element-free Galerkin （DSEFG） method for three-dimensional potential problems, and the corresponding formulae are obtained. The main idea of the DSEFG method is that a three-dimensional potential problem can be transformed into a series of two-dimensional problems. For these two-dimensional problems, the improved moving least-squares （IMLS） approximation is applied to construct the shape function, which uses an orthogonal function system with a weight function as the basis functions. The Galerkin weak form is applied to obtain a discretized system equation, and the penalty method is employed to impose the essential boundary condition. The finite difference method is selected in the splitting direction. For the purposes of demonstration, some selected numerical examples are solved using the DSEFG method. The convergence study and error analysis of the DSEFG method are presented. The numerical examples show that the DSEFG method has greater computational precision and computational efficiency than the IEFG method.
基金sponsored by National Key Research and Development Program of China(2018YFC0807000)The Spark Program of Earthquake Technology of CEA(XH20070Y)The Earthquake Tracking Task of CEA(2021010221).
文摘This article analyzes the relationship between the water level and the water-tube tilting in Shuangyang lake,based on the differential deformation features reflected by the NS and EW components of the water-tube tiltmeter.The results show a good spatiotemporal consistency between the variation of water level and the NS tilt component,which is considered to be affected by the magnitude and duration of the water level variation in Shuangyang Lake.The article uses Landsat remote sensing image data to extract the water boundary of Shuan-gyang Lake,and takes advantage of the finite element numerical simulation method to build three-dimensional models for different geological structural conditions of the Shuangyang seismostation.The simulation results show that when the underground medium is granite,the effect of water level variation on the vertical displacement of the surface is non-directional.With a 50-m soil layer in Model 2,the simulated NS tilt variation is equivalent to the actual observed water-tube tiltmeter NS component when the water level variation is 0.44 m and 0.8m.When the variation of water level reaches 2.0m,the simulation result of the NS component is 79.6 ms,which is slightly larger than the observed result of 60.32 ms.However,the simulation results show that the variation of the EW component is significantly smaller than that of the NS one.Due to the fact that the Shuan-gyang lake is long in the NS direction and short in the EW direction,the existence of the soil layer tends to generate ground deformation along the NS direction in the vicinity of the lake after the increase of water level,thereby resulting in the difference of the ground deformation in the two directions.
基金This research was partly supported by the National Basic Research Program of China973 Program under Grant No.2011CB706903+3 种基金the Program for New Century Excellent Talents in University under Grant No.NCET-09-0438the National Natural Science Foundation of China under Grant No.10801067the Fundamental Research Funds for the Central Universities under Grant No.lzujbky-2010-63No.lzujbky-2012-k26。
文摘We present the finite difference/element method for a two-dimensional modified fractional diffusion equation.The analysis is carried out first for the time semi-discrete scheme,and then for the full discrete scheme.The time discretization is based on the L1-approximation for the fractional derivative terms and the second-order backward differentiation formula for the classical first order derivative term.We use finite element method for the spatial approximation in full discrete scheme.We show that both the semi-discrete and full discrete schemes are unconditionally stable and convergent.Moreover,the optimal convergence rate is obtained.Finally,some numerical examples are tested in the case of one and two space dimensions and the numerical results confirm our theoretical analysis.
