Numerical analyses of earthquake effects on the deformation, stability, and load transfer of a slope covered by deposits are traditionally based on the assumption that the slope is a continuum. It would be problem...Numerical analyses of earthquake effects on the deformation, stability, and load transfer of a slope covered by deposits are traditionally based on the assumption that the slope is a continuum. It would be problematic, however, to extend these approaches to the simulation of the slide, collapse and disintegration of the deposits under seismic loading. Contrary to this, a discrete element method (DEM) provides a means to consider large displacement and rotation of the non-continuum. To take the advantages of both methods of continuum and non- continuum analyses, seismic responses of a slope covered by deposits are studied by coupling a twodimensional (a-D) finite difference method and a 2-D DEM, with the bedrock being modelled by the finite difference grids and the deposits being represented by disks. A smooth transition across the boundaries of the continuous/discontinuous domains is obtained by imposing the compatibility condition and equilibrium condition along their interfaces. In the course of computation, the same time-step value is chosen for both continuous and discontinuous domains. The free-field boundaries are adopted for lateral grids of bedrock domain to eliminate the radiation damping effect. When the static equilibrium under gravity load is obtained, dynamic calculation begins under excitation of the seismic wave input from the continuum model bottom. In this way, responses to the earthquake of a slope covered by deposits are analyzed dynamically. Combined with field monitoring data, deformation and stability of the slope are discussed. The effects of the relevant parameters of spectrum characteristic, duration, andpeak acceleration of seismic waves are further investigated and explained from the simulations.展开更多
A procedure of kinematic analysis is presented in this study to assess the reinforcement force of geosynthetics required under seismic loadings, particularly for steep slopes which are hardly able to maintain its stab...A procedure of kinematic analysis is presented in this study to assess the reinforcement force of geosynthetics required under seismic loadings, particularly for steep slopes which are hardly able to maintain its stability. Note that curved sloping surfaces widely exist in natural slopes, but existing literatures were mainly focusing on a planar surface in theoretical derivation, due to complicated calculations. Moreover, the non-uniform soil properties cannot be accounted for in conventional upper bound analysis. Pseudo-dynamic approach is used to represent horizontal and vertical accelerations which vary with time and space. In an effort to resolve the above problems, the discretization technique is developed to generate a discretized failure mechanism, decomposing the whole failure block into various components. An elementary analysis permits calculations of rates of work done by external and internal forces. Finally, the upper bound solution of the required reinforcement force is formulated based on the work rate-based balance equation. A parametric study is carried out to give insights on the implication of influential factors on the performance of geosynthetic-reinforced steep slopes.展开更多
This paper presents an approach to calculate dispersion penalty for VSR-1 optical links.Based on parameters of a specific VSR-1 link,dispersion penalties are computed for various modal dispersion bandwidths respective...This paper presents an approach to calculate dispersion penalty for VSR-1 optical links.Based on parameters of a specific VSR-1 link,dispersion penalties are computed for various modal dispersion bandwidths respectively.The worst-case eye closure is expressed numerically by using the signal waveform at time 0,and the signal waveform is obtained in frequency domain through FFT algorithm.By this approach,the dispersion penalty is determined by the shape of transfer functions of the various components in the links.To simplify the derivation of multimode fiber link transfer function,a Gaussian form of normalized impulse response is used.This calculation approach can be used to estimate the worst-case dispersion penalty of VSR-1 links in the link budget analysis.展开更多
The quantitative evaluation of errors involved in a particular numerical modelling is of prime importance for the effectiveness and reliability of the method. Errors in Distinct Element Modelling are generated mainly ...The quantitative evaluation of errors involved in a particular numerical modelling is of prime importance for the effectiveness and reliability of the method. Errors in Distinct Element Modelling are generated mainly through three resources as simplification of physical model, determination of parameters and boundary conditions. A measure of errors which represent the degree of numerical solution 'close to true value' is proposed through fuzzy probability in this paper. The main objective of this paper is to estimate the reliability of Distinct Element Method in rock engineering practice by varying the parameters and boundary conditions. The accumulation laws of standard errors induced by improper determination of parameters and boundary conditions are discussed in delails. Furthermore, numerical experiments are given to illustrate the estimation of fuzzy reliability. Example shows that fuzzy reliability falls between 75%-98% when the relative standard errors of input data is under 10 %.展开更多
In this paper, we apply homotopy analysis method to solve discrete mKdV equation and successfully obtain the bell-shaped solitary solution to mKdV equation. Comparison between our solution and the exact solution shows...In this paper, we apply homotopy analysis method to solve discrete mKdV equation and successfully obtain the bell-shaped solitary solution to mKdV equation. Comparison between our solution and the exact solution shows that homotopy analysis method is effective and validity in solving hybrid nonlinear problems, including solitary solution of difference-differential equation.展开更多
A method combining the pseudo-dynamic approach and discretization technique is carried out for computing the active earth pressure.Instead of using a presupposed failure mechanism,discretization technique is introduce...A method combining the pseudo-dynamic approach and discretization technique is carried out for computing the active earth pressure.Instead of using a presupposed failure mechanism,discretization technique is introduced to generate the potential failure surface,which is applicable to the case that soil strength parameters have spatial variability.For the purpose of analyzing the effect of earthquake,pseudo-dynamic approach is adopted to introduce the seismic forces,which can take into account the dynamic properties of seismic acceleration.A new type of micro-element is used to calculate the rate of work of external forces and the rate of internal energy dissipation.The analytical expression of seismic active earth pressure coefficient is deduced in the light of upper bound theorem and the corresponding upper bound solutions are obtained through numerical optimization.The method is validated by comparing the results of this paper with those reported in literatures.The parametric analysis is finally presented to further expound the effect of diverse parameters on active earth pressure under non-uniform soil.展开更多
In physics,the Klein-Gordon equation describes the motion of a quantum scalar or pseudoscalar field.Itis important to find actual values of its solutions in general timespace manifold.The paper deals with description ...In physics,the Klein-Gordon equation describes the motion of a quantum scalar or pseudoscalar field.Itis important to find actual values of its solutions in general timespace manifold.The paper deals with description ofdiscrete exterior calculus method for solving this equation numerically on space manifold and the time.The analysis ofstable condition and error for this method is also accomplished.展开更多
In order to obtain an expected numerical solution, a fully discrete discontinuous Galerkin method is applied to a kind of reactive transport problems in two dimension. That is to say, the space variable is discretized...In order to obtain an expected numerical solution, a fully discrete discontinuous Galerkin method is applied to a kind of reactive transport problems in two dimension. That is to say, the space variable is discretized with the symmetric interior penalty Calerkin method (SIPG), and the time variable is done with the backward Euler method. Making use of the duality technique, hp approximation properties and the interpolation theory, a residual-type a posteriori error estimation is achieved, which can be used for adaptivity. Compared with the analyses of semi-discretization, the current presentation is more challenging and more significant.展开更多
基金the National Basic Research Program of China (Grant No. 2008CB425802)
文摘Numerical analyses of earthquake effects on the deformation, stability, and load transfer of a slope covered by deposits are traditionally based on the assumption that the slope is a continuum. It would be problematic, however, to extend these approaches to the simulation of the slide, collapse and disintegration of the deposits under seismic loading. Contrary to this, a discrete element method (DEM) provides a means to consider large displacement and rotation of the non-continuum. To take the advantages of both methods of continuum and non- continuum analyses, seismic responses of a slope covered by deposits are studied by coupling a twodimensional (a-D) finite difference method and a 2-D DEM, with the bedrock being modelled by the finite difference grids and the deposits being represented by disks. A smooth transition across the boundaries of the continuous/discontinuous domains is obtained by imposing the compatibility condition and equilibrium condition along their interfaces. In the course of computation, the same time-step value is chosen for both continuous and discontinuous domains. The free-field boundaries are adopted for lateral grids of bedrock domain to eliminate the radiation damping effect. When the static equilibrium under gravity load is obtained, dynamic calculation begins under excitation of the seismic wave input from the continuum model bottom. In this way, responses to the earthquake of a slope covered by deposits are analyzed dynamically. Combined with field monitoring data, deformation and stability of the slope are discussed. The effects of the relevant parameters of spectrum characteristic, duration, andpeak acceleration of seismic waves are further investigated and explained from the simulations.
文摘A procedure of kinematic analysis is presented in this study to assess the reinforcement force of geosynthetics required under seismic loadings, particularly for steep slopes which are hardly able to maintain its stability. Note that curved sloping surfaces widely exist in natural slopes, but existing literatures were mainly focusing on a planar surface in theoretical derivation, due to complicated calculations. Moreover, the non-uniform soil properties cannot be accounted for in conventional upper bound analysis. Pseudo-dynamic approach is used to represent horizontal and vertical accelerations which vary with time and space. In an effort to resolve the above problems, the discretization technique is developed to generate a discretized failure mechanism, decomposing the whole failure block into various components. An elementary analysis permits calculations of rates of work done by external and internal forces. Finally, the upper bound solution of the required reinforcement force is formulated based on the work rate-based balance equation. A parametric study is carried out to give insights on the implication of influential factors on the performance of geosynthetic-reinforced steep slopes.
基金Supported by"863"Hi-Tech Research and Development Program(2005AA311030) and the National Natural Science Foundationof China(Grant No.60502005)
文摘This paper presents an approach to calculate dispersion penalty for VSR-1 optical links.Based on parameters of a specific VSR-1 link,dispersion penalties are computed for various modal dispersion bandwidths respectively.The worst-case eye closure is expressed numerically by using the signal waveform at time 0,and the signal waveform is obtained in frequency domain through FFT algorithm.By this approach,the dispersion penalty is determined by the shape of transfer functions of the various components in the links.To simplify the derivation of multimode fiber link transfer function,a Gaussian form of normalized impulse response is used.This calculation approach can be used to estimate the worst-case dispersion penalty of VSR-1 links in the link budget analysis.
文摘The quantitative evaluation of errors involved in a particular numerical modelling is of prime importance for the effectiveness and reliability of the method. Errors in Distinct Element Modelling are generated mainly through three resources as simplification of physical model, determination of parameters and boundary conditions. A measure of errors which represent the degree of numerical solution 'close to true value' is proposed through fuzzy probability in this paper. The main objective of this paper is to estimate the reliability of Distinct Element Method in rock engineering practice by varying the parameters and boundary conditions. The accumulation laws of standard errors induced by improper determination of parameters and boundary conditions are discussed in delails. Furthermore, numerical experiments are given to illustrate the estimation of fuzzy reliability. Example shows that fuzzy reliability falls between 75%-98% when the relative standard errors of input data is under 10 %.
基金the State Key Basic Research Program of China under Grant No.2004CB318000
文摘In this paper, we apply homotopy analysis method to solve discrete mKdV equation and successfully obtain the bell-shaped solitary solution to mKdV equation. Comparison between our solution and the exact solution shows that homotopy analysis method is effective and validity in solving hybrid nonlinear problems, including solitary solution of difference-differential equation.
基金Projects(51908557,51378510)supported by the National Natural Science Foundation of China。
文摘A method combining the pseudo-dynamic approach and discretization technique is carried out for computing the active earth pressure.Instead of using a presupposed failure mechanism,discretization technique is introduced to generate the potential failure surface,which is applicable to the case that soil strength parameters have spatial variability.For the purpose of analyzing the effect of earthquake,pseudo-dynamic approach is adopted to introduce the seismic forces,which can take into account the dynamic properties of seismic acceleration.A new type of micro-element is used to calculate the rate of work of external forces and the rate of internal energy dissipation.The analytical expression of seismic active earth pressure coefficient is deduced in the light of upper bound theorem and the corresponding upper bound solutions are obtained through numerical optimization.The method is validated by comparing the results of this paper with those reported in literatures.The parametric analysis is finally presented to further expound the effect of diverse parameters on active earth pressure under non-uniform soil.
基金Supported by China Postdoctoral Science Foundation under Grant No.20090460102 Zhejiang Province Postdoctoral Science Foundation,National Key Basic Research Program of China under Grant No.2004CB318000 National Natural Science Foundation of China under Grant No.10871170
文摘In physics,the Klein-Gordon equation describes the motion of a quantum scalar or pseudoscalar field.Itis important to find actual values of its solutions in general timespace manifold.The paper deals with description ofdiscrete exterior calculus method for solving this equation numerically on space manifold and the time.The analysis ofstable condition and error for this method is also accomplished.
基金supported by Hunan Provincial Natural Science Foundation of China under Grant No. 10JJ3021Scientific Research Fund of Hunan Provincial Education Department under Grant No.11B032the Planned Science and Technology Project of Hunan Province and Aid program for Science and Technology Innovative Research Team in Higher Educational Institutions of Hunan Province
文摘In order to obtain an expected numerical solution, a fully discrete discontinuous Galerkin method is applied to a kind of reactive transport problems in two dimension. That is to say, the space variable is discretized with the symmetric interior penalty Calerkin method (SIPG), and the time variable is done with the backward Euler method. Making use of the duality technique, hp approximation properties and the interpolation theory, a residual-type a posteriori error estimation is achieved, which can be used for adaptivity. Compared with the analyses of semi-discretization, the current presentation is more challenging and more significant.