G.S. Li¹, D. Yao², Y.Z. Wang³, H.Y. Jiang²

CMES-Computer Modeling in Engineering & Sciences, Vol.51, No.1, pp. 53-72, 2009, DOI:10.3970/cmes.2009.051.053

Abstract In this paper, an undisturbed soil-column infiltrating experiment is investigated, and a mathematical model describing multi-components solutes transport behaviors in the column is put forward by combing hydro-chemical analysis with advection dispersion mechanisms, which is a group of advection-dispersion-reaction partial differential equations. Since the model involving six reaction coefficients which can not be obtained directly, an optimal perturbation regularization algorithm of determining these parameters is performed, and numerical simulations under different conditions are carried out. Furthermore, the inversion algorithm is applied to solve the real inverse problem by utilizing the measured breakthrough data. The reconstruction data basically coincide with the… More >

Ehsan Jabbari¹, Behrouz Gatmiri^2,3

CMES-Computer Modeling in Engineering & Sciences, Vol.18, No.1, pp. 31-44, 2007, DOI:10.3970/cmes.2007.018.031

Abstract In this paper after a discussion about the evolution of the unsaturated soils' governing differential equations and a brief history of the Green's functions for porous media, the governing equations, i.e., the mathematical model in the presence of heat effects are presented and simplified so as the derivation of the associated Green's functions be in the realm of possibility. The thermal two- and three-dimensional, full- and half-space Green's functions for unsaturated porous media, although in a relatively simplified form, are being introduced for the first time, following the previous works of the authors. The derived Green's functions have been demonstrated… More >

Daoming Shen¹, Hua Si^1,*, Jinhong Xia¹, Shunqun Li²

FDMP-Fluid Dynamics & Materials Processing, Vol.15, No.1, pp. 63-76, 2019, DOI:10.32604/fdmp.2019.04799

Abstract The correct determination of thermal parameters, such as thermal conductivity and specific heat of soil during freezing, is the most important and basic problem for the construction of an appropriate freezing method. In this study, a calculation model of three stages of soil temperature was established. At the unfrozen and frozen stages, the specific temperatures of dry soil, water, and ice are known. According to the principle of superposition, a calculation model of unfrozen and frozen soils can be established. Informed by a laboratory experiment, the latent heat of the adjacent zone was calculated for the freezing stage based on… More >

Chunguang Li^{1, *}, Cuihua Li^{1, 2}, Cong Sun³, Hong Zheng¹

CMC-Computers, Materials & Continua, Vol.54, No.1, pp. 21-41, 2018, DOI:10.3970/cmc.2018.054.021

Abstract Previous approaches can only tackle anisotropic problems with cohesion varying with direction. A novel linearization of the Mohr-Coulomb yield criterion associated with plane strain problem has been achieved by simulating the Mohr’s circle with orientation lines in σ-τ space, which allows for lower bound solution of soils with cohesion and friction coefficient varying with direction. The finite element lower limit analysis formulation using the modified anisotropic yield criterion is then developed. Several examples are given to illustrate the capability and effectiveness of the proposed numerical procedure for computing rigorous lower bounds for anisotropic soils. More >

Yang Pu¹, Xiang Zhihai¹, Hu Xiasong², Li Guorong², Zhu Haili², Mao XiaoqinCen², Zhangzhi^1,3

CMC-Computers, Materials & Continua, Vol.10, No.2, pp. 117-138, 2009, DOI:10.3970/cmc.2009.010.117

Abstract This paper studies the influence of root reinforcement on shallow soil protection by using Finite Element (FE) method. Taking the root-soil composite as a periodic material, the homogenization method is used to construct a Representative Volume Element (RVE) that consists of roots and soil. This RVE is discretized by a two-dimensional (2-D) FE mesh, while special formulation is established so that this model is capable of describing three-dimensional (3-D) deformations when the strain is invariant along the fiber axis. The important effect of debonding on the interface between the fiber and the matrix is also considered by using a special… More >