The pile-soil system interaction computational model in liquefaction-induced lateral spreading ground was established by the finite difference numerical method.Considering an elastic-plastic subgrade reaction method,n...The pile-soil system interaction computational model in liquefaction-induced lateral spreading ground was established by the finite difference numerical method.Considering an elastic-plastic subgrade reaction method,numerical methods involving finite difference approach of pile in liquefaction-induced lateral spreading ground were derived and implemented into a finite difference program.Based on the monotonic loading tests on saturated sand after liquefaction,the liquefaction lateral deformation of the site where group piles are located was predicted.The effects of lateral ground deformation after liquefaction on a group of pile foundations were studied using the fmite difference program mentioned above,and the failure mechanism of group piles in liquefaction-induced lateral spreading ground was obtained.The applicability of the program was preliminarily verified.The results show that the bending moments at the interfaces between liquefied and non-liquefied soil layers are larger than those at the pile＇s top when the pile＇s top is embedded.The value of the additional static bending moment is larger than the peak dynamic bending moment during the earthquake,so in the pile foundation design,more than the superstructure＇s dynamics should be considered and the effect of lateral ground deformation on pile foundations cannot be neglected.展开更多
To predict the heat diffusion in a given region over time, it is often necessary to find the numerical solution for heat equation. However, the computational domain of classical numerical methods are limited to fiat s...To predict the heat diffusion in a given region over time, it is often necessary to find the numerical solution for heat equation. However, the computational domain of classical numerical methods are limited to fiat spacetime. With the techniques of discrete differential calculus, we propose two unconditional stable numerical schemes for simulation heat equation on space manifold and time. The analysis of their stability and error is accomplished by the use of maximum principle.展开更多
To overcome the deficiencies of conventional geosynthetic-reinforced and pile-supported (GRPS) embankment, a new improvement technique, fixed geosynthetic technique of GRPS embankment (FGT embankment), was developed a...To overcome the deficiencies of conventional geosynthetic-reinforced and pile-supported (GRPS) embankment, a new improvement technique, fixed geosynthetic technique of GRPS embankment (FGT embankment), was developed and introduced. Based on the discussion about the load transfer mechanism of FGT embankment, a simplified check method of the requirement of geosynthetic tensile strength and a mechanical model of the FGT embankment were proposed. Two conditions, the pile cap and pile beam conditions are considered in the mechanical model. The finite difference method is used to solve the mechanical model owing to the complexity of the differential equations and the soil strata. Then, the numerical procedure is programmed. Finally, a field test is conducted to verify the mechanical model and the calculated results are in good agreement with field measured data.展开更多
In this paper we consider averaging and finite difference methods for solving the 3-D boundary-value problem in a multilayered domain. We consider the metal concentration in the 3 layered peat blocks. Using experiment...In this paper we consider averaging and finite difference methods for solving the 3-D boundary-value problem in a multilayered domain. We consider the metal concentration in the 3 layered peat blocks. Using experimental data the mathematical model for calculating the concentration of metal at different points in peat layers is developed. A specific feature of these problems is that it is necessary to solve the 3-D boundary-value problems for the partial differential equations （PDEs） of the elliptic type of second order with piece-wise diffusion coefficients in the three layer domain. We develop here a finite-difference method for solving a problem of the above type with the periodical boundary condition in x direction. This procedure allows reducing the 3-D problem to a system of 2-D problems by using a circulant matrix.展开更多
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.展开更多
Optical waveguide is the main element in integrated optics. Therefore many numerical methods are used on these elements of integrated optics. Simulation response of an optical slab waveguide used in integrated optics ...Optical waveguide is the main element in integrated optics. Therefore many numerical methods are used on these elements of integrated optics. Simulation response of an optical slab waveguide used in integrated optics needs such numerical methods. These methods must be precise and useful in terms of memory capacity and time duration. In this paper, we study basic analytical and finite difference methods to determine the effective refractive index of AIGaAs-GaAs slab waveguide. Also, appropriate effective refractive index value is obtained with respect to number of grid points and number of matrix sizes. Finally, the validity of the obtained values by both methods is compared to using waveguide type.展开更多
To guarantee the accuracy of error analysis and evaluate the manufacturing tolerance s influence,anumerical error analysis method for parallel kinematic machines (PKMs) is presented in this paper.Quasi-Newton method a...To guarantee the accuracy of error analysis and evaluate the manufacturing tolerance s influence,anumerical error analysis method for parallel kinematic machines (PKMs) is presented in this paper.Quasi-Newton method and genetic algorithm are introduced for the forward kinematic solution.Based onthe inverse and forward kinematic solutions,the end-effector s error calculation procedure is developed.To solve the accuracy problem caused by the length and angular parameters' different units,a normalizationmethod is proposed based on the manufacturing tolerance.Comparison between the error analysis resultscalculated by the traditional method and the numerical method for a 4RRR PKM shows that,this numericalerror analysis method is more accurate,simpler,and can evaluate the machine s real error basedon the manufacturing tolerance.展开更多
Homogenization is concerned with obtaining the average properties of a material. The problem on its own has no easy solution, except in cases like the periodic case, when it can be obtained in closed form. In this pap...Homogenization is concerned with obtaining the average properties of a material. The problem on its own has no easy solution, except in cases like the periodic case, when it can be obtained in closed form. In this paper we consider a numerical solution of the elliptic homogenization problem for the case of rapidly varying tensor or boundary conditions. The method makes use of an adaptive finite element method to correctly capture the rapid change in the tensor or boundary condition. In the numerical experiments we vary the mesh size and do a posteriori error analysis on test problems.展开更多
The damage identification is made by the numerical simulation analysis of a five-storey-and-two-span RC frame structure, using improved and unimproved direct analytical method respectively; and the fundamental equatio...The damage identification is made by the numerical simulation analysis of a five-storey-and-two-span RC frame structure, using improved and unimproved direct analytical method respectively; and the fundamental equations were solved by the minimal least square method （viz. general inverse method）. It demonstrates that the feasibility and the accuracy of the present approach were impoved significantly, compared with the result of unimproved damage identification.展开更多
Semiparametric regression models and estimating covariance functions are very useful for longitudinal study. To heed the positive-definiteness constraint, we adopt the modified Cholesky decomposition approach to decom...Semiparametric regression models and estimating covariance functions are very useful for longitudinal study. To heed the positive-definiteness constraint, we adopt the modified Cholesky decomposition approach to decompose the covariance structure. Then the covariance structure is fitted by a semiparametric model by imposing parametric within-subject correlation while allowing the nonparametric variation function. We estimate regression functions by using the local linear technique and propose generalized estimating equations for the mean and correlation parameter. Kernel estimators are developed for the estimation of the nonparametric variation function. Asymptotic normality of the the resulting estimators is established. Finally, the simulation study and the real data analysis are used to illustrate the proposed approach.展开更多
The optimal nearly-analytic discrete(ONAD) method is a new numerical method developed in recent years for solving the wave equation.Compared with other methods,such as popularly-used finite-difference methods,the ONAD...The optimal nearly-analytic discrete(ONAD) method is a new numerical method developed in recent years for solving the wave equation.Compared with other methods,such as popularly-used finite-difference methods,the ONAD method can effectively suppress the numerical dispersion when coarse grids are used.In this paper,the ONAD method is extended to solve the 2-dimensional SH-wave equation in the spherical coordinates.To investigate the accuracy and the efficiency of the ONAD method,we compare the numerical results calculated by the ONAD method and other methods for both the homogeneous model and the inhomogeneous IASP91 model.The comparisons indicate that the ONAD method for solving the SH-wave equation in the spherical coordinates has the advantages of less numerical dispersion,small memory requirement for computer codes,and fast calculation.As an application,we use the ONAD method to simulate the SH-wave propagating between the Earth's surface and the core-mantle boundary(CMB).Meanwhile,we investigate the SH-wave propagating in the mantle through analyzing the wave-field snapshots in different times and synthetic seismograms.展开更多
基金Supported by grants from Yunnan Province Science and Technology Department(Grant No.2011BC010)Kunming Science and Technology Bureau(Grant No.10H100102)。
文摘经典天体测量仪器以铅垂线为基准测量本地的天文经纬度,因而能探测到本地铅垂线的偏转.本地铅垂线的偏转代表着测站周围重力场的变化,而这一变化与地下物质的再分布相关,由此有望帮助了解地下物质变化的情况.将本地的铅垂线偏转(Plumb Line Variation,PLV)与地下物质变化联系起来,建立了双质量体模型.在97°E-107°E和21°N-29°N范围内,考虑了3480个包含正、负质量变化(相对地质背景)的质量体组合算例,并利用差分进化(Differential Evolutionary,DE)算法进行了解算.解算得到的质量体相对于模拟值的位置误差小于米量级,质量误差小于10^(11) kg,结果精度较高.
基金Project(51109208)supported by the National Natural Science Foundation of ChinaProject(2013M531688)supported by the Postdoctoral Science Foundation of China+1 种基金Project(Z012009)supported by the Open Research Fund of State Key Laboratory of Geomechanics and Geotechnical Engineering(Institute of Rock and Soil Mechanics,Chinese Academy of Sciences)Project(CKSF2012054)supported by the Foundation of Changjiang River Scientific Research Institute,China
文摘The pile-soil system interaction computational model in liquefaction-induced lateral spreading ground was established by the finite difference numerical method.Considering an elastic-plastic subgrade reaction method,numerical methods involving finite difference approach of pile in liquefaction-induced lateral spreading ground were derived and implemented into a finite difference program.Based on the monotonic loading tests on saturated sand after liquefaction,the liquefaction lateral deformation of the site where group piles are located was predicted.The effects of lateral ground deformation after liquefaction on a group of pile foundations were studied using the fmite difference program mentioned above,and the failure mechanism of group piles in liquefaction-induced lateral spreading ground was obtained.The applicability of the program was preliminarily verified.The results show that the bending moments at the interfaces between liquefied and non-liquefied soil layers are larger than those at the pile＇s top when the pile＇s top is embedded.The value of the additional static bending moment is larger than the peak dynamic bending moment during the earthquake,so in the pile foundation design,more than the superstructure＇s dynamics should be considered and the effect of lateral ground deformation on pile foundations cannot be neglected.
基金Supported by China Postdoctoral Science Foundation under Grant No.20090460102
文摘To predict the heat diffusion in a given region over time, it is often necessary to find the numerical solution for heat equation. However, the computational domain of classical numerical methods are limited to fiat spacetime. With the techniques of discrete differential calculus, we propose two unconditional stable numerical schemes for simulation heat equation on space manifold and time. The analysis of their stability and error is accomplished by the use of maximum principle.
基金Project(51278216) supported by the National Natural Science Foundation of ChinaProject(20091341) supported by the Scientific Research Foundation for Returned Overseas Chinese Scholars,Ministry of Education,ChinaProject(HF-08-01-2011-240) supported by the Graduates’ Innovation Fund of Huazhong University of Science and Technology,China
文摘To overcome the deficiencies of conventional geosynthetic-reinforced and pile-supported (GRPS) embankment, a new improvement technique, fixed geosynthetic technique of GRPS embankment (FGT embankment), was developed and introduced. Based on the discussion about the load transfer mechanism of FGT embankment, a simplified check method of the requirement of geosynthetic tensile strength and a mechanical model of the FGT embankment were proposed. Two conditions, the pile cap and pile beam conditions are considered in the mechanical model. The finite difference method is used to solve the mechanical model owing to the complexity of the differential equations and the soil strata. Then, the numerical procedure is programmed. Finally, a field test is conducted to verify the mechanical model and the calculated results are in good agreement with field measured data.
文摘In this paper we consider averaging and finite difference methods for solving the 3-D boundary-value problem in a multilayered domain. We consider the metal concentration in the 3 layered peat blocks. Using experimental data the mathematical model for calculating the concentration of metal at different points in peat layers is developed. A specific feature of these problems is that it is necessary to solve the 3-D boundary-value problems for the partial differential equations （PDEs） of the elliptic type of second order with piece-wise diffusion coefficients in the three layer domain. We develop here a finite-difference method for solving a problem of the above type with the periodical boundary condition in x direction. This procedure allows reducing the 3-D problem to a system of 2-D problems by using a circulant matrix.
基金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.
文摘Optical waveguide is the main element in integrated optics. Therefore many numerical methods are used on these elements of integrated optics. Simulation response of an optical slab waveguide used in integrated optics needs such numerical methods. These methods must be precise and useful in terms of memory capacity and time duration. In this paper, we study basic analytical and finite difference methods to determine the effective refractive index of AIGaAs-GaAs slab waveguide. Also, appropriate effective refractive index value is obtained with respect to number of grid points and number of matrix sizes. Finally, the validity of the obtained values by both methods is compared to using waveguide type.
基金Supported by the National High Technology Research and Development Programme of China ( No. 2007AA041901 )the National Natural Science Foundation of China ( No. 50775117 )+1 种基金the National S&T Major Project ( No. 2009XZ04001-025 )the Technology Innovation Fund of AVIC ( No.2009E 13224 )
文摘To guarantee the accuracy of error analysis and evaluate the manufacturing tolerance s influence,anumerical error analysis method for parallel kinematic machines (PKMs) is presented in this paper.Quasi-Newton method and genetic algorithm are introduced for the forward kinematic solution.Based onthe inverse and forward kinematic solutions,the end-effector s error calculation procedure is developed.To solve the accuracy problem caused by the length and angular parameters' different units,a normalizationmethod is proposed based on the manufacturing tolerance.Comparison between the error analysis resultscalculated by the traditional method and the numerical method for a 4RRR PKM shows that,this numericalerror analysis method is more accurate,simpler,and can evaluate the machine s real error basedon the manufacturing tolerance.
文摘Homogenization is concerned with obtaining the average properties of a material. The problem on its own has no easy solution, except in cases like the periodic case, when it can be obtained in closed form. In this paper we consider a numerical solution of the elliptic homogenization problem for the case of rapidly varying tensor or boundary conditions. The method makes use of an adaptive finite element method to correctly capture the rapid change in the tensor or boundary condition. In the numerical experiments we vary the mesh size and do a posteriori error analysis on test problems.
文摘The damage identification is made by the numerical simulation analysis of a five-storey-and-two-span RC frame structure, using improved and unimproved direct analytical method respectively; and the fundamental equations were solved by the minimal least square method （viz. general inverse method）. It demonstrates that the feasibility and the accuracy of the present approach were impoved significantly, compared with the result of unimproved damage identification.
基金supported by National Natural Science Foundation of China (GrantNos.10931002,10911120386)
文摘Semiparametric regression models and estimating covariance functions are very useful for longitudinal study. To heed the positive-definiteness constraint, we adopt the modified Cholesky decomposition approach to decompose the covariance structure. Then the covariance structure is fitted by a semiparametric model by imposing parametric within-subject correlation while allowing the nonparametric variation function. We estimate regression functions by using the local linear technique and propose generalized estimating equations for the mean and correlation parameter. Kernel estimators are developed for the estimation of the nonparametric variation function. Asymptotic normality of the the resulting estimators is established. Finally, the simulation study and the real data analysis are used to illustrate the proposed approach.
基金supported by National Science Fund of Distinguished Young Scholars of China(Grant No. 40725012)40821002)National Natural Science Foundation of China (Grant No. 41074073)
文摘The optimal nearly-analytic discrete(ONAD) method is a new numerical method developed in recent years for solving the wave equation.Compared with other methods,such as popularly-used finite-difference methods,the ONAD method can effectively suppress the numerical dispersion when coarse grids are used.In this paper,the ONAD method is extended to solve the 2-dimensional SH-wave equation in the spherical coordinates.To investigate the accuracy and the efficiency of the ONAD method,we compare the numerical results calculated by the ONAD method and other methods for both the homogeneous model and the inhomogeneous IASP91 model.The comparisons indicate that the ONAD method for solving the SH-wave equation in the spherical coordinates has the advantages of less numerical dispersion,small memory requirement for computer codes,and fast calculation.As an application,we use the ONAD method to simulate the SH-wave propagating between the Earth's surface and the core-mantle boundary(CMB).Meanwhile,we investigate the SH-wave propagating in the mantle through analyzing the wave-field snapshots in different times and synthetic seismograms.
