This study proposes a general imperfect thermal contact model to predict the thermal contact resistance at the interface among multi-layered composite structures.Based on the Green-Lindsay(GL)thermoelastic theory,semi...This study proposes a general imperfect thermal contact model to predict the thermal contact resistance at the interface among multi-layered composite structures.Based on the Green-Lindsay(GL)thermoelastic theory,semi analytical solutions of temperature increment and displacement of multi-layered composite structures are obtained by using the Laplace transform method,upon which the effects of thermal resistance coefficient,partition coefficient,thermal conductivity ratio and heat capacity ratio on the responses are studied.The results show that the generalized imperfect thermal contact model can realistically describe the imperfect thermal contact problem.Accordingly,it may degenerate into other thermal contact models by adjusting the thermal resistance coefficient and partition coefficient.展开更多
Considering the fact that in some complex cases,plate anchors are buried in multi-layered geotechnical materials,the ultimate dynamic analysis was performed to investigate the uplift capacity and failure mechanism of ...Considering the fact that in some complex cases,plate anchors are buried in multi-layered geotechnical materials,the ultimate dynamic analysis was performed to investigate the uplift capacity and failure mechanism of shallow strips and circular plate anchors in multi-layered soils.The nonlinear strength criterion and non-associated flow rule of geotechnical materials were introduced to investigate the influence of nonuniformity on the pullout performance and failure mechanism of shallow plate anchors.The expressions of the detaching curves or surfaces were obtained to reflect the failure mechanism,which can be used to figure out the ultimate uplift capacity and failure range.The results are generally in agreement with the numerical simulations and previous research.The effects of various parameters on the ultimate uplift capacity and failure mechanism of plate anchors in multi-layered soils were investigated,and it is found that the ultimate uplift capacity and failure range of shallow anchors increase with the increase of initial cohesion and dilatancy coefficient,but decrease with the unit weight,axial tensile stress and nonlinear coefficient.展开更多
Two modifications for the basic Barcelona model(BBM) are present. One is the replacement of the net stress by the average skeleton stress in unsaturated soil modeling, and the other is the adoption of an expression fo...Two modifications for the basic Barcelona model(BBM) are present. One is the replacement of the net stress by the average skeleton stress in unsaturated soil modeling, and the other is the adoption of an expression for the load-collapse(LC) yield surface that can match flexibly the normal compression lines at different suctions. The predictions of the modified BBM for the controlled-suction triaxial test on the unsaturated compacted clay are presented and compared with the experimental results. A good agreement between the predicted and experimental results demonstrates the reasonability of the modified BBM. On this basis, the coupled processes of groundwater flow and soil deformation in a homogeneous soil slope under a long heavy rainfall are simulated with the proposed elasto-plastic model. The numerical results reveal that the failure of a slope under rainfall infiltration is due to both the reduction of soil suction and the significant rise in groundwater table. The evolution of the displacements is greatly related to the change of suction. The maximum collapse deformation happens near the surface of slope where infiltrated rainwater can quickly reach. The results may provide a helpful reference for hazard assessment and control of rainfall-induced landslides.展开更多
A new approach is proposed to analyze the settlement behavior for single pile embedded in layered soils. Firstly, soil layers surrounding pile shaft are simulated by using distributed Voigt model, and finite soil laye...A new approach is proposed to analyze the settlement behavior for single pile embedded in layered soils. Firstly, soil layers surrounding pile shaft are simulated by using distributed Voigt model, and finite soil layers under the pile end are assumed to be virtual soil-pile whose cross-section area is the same as that of the pile shaft. Then, by means of Laplace transform and impedance function transfer method to solve the static equilibrium equation of pile, the analytical solution of the displacement impedance fimction at the pile head is derived. Furthermore, the analytical solution of the settlement at the head of single pile is theoretically derived by virtue of convolution theorem. Based on these solutions, the influences of parameters of soil-pile system on the settlement behavior for single pile are analyzed. Also, comparison of the load-settlement response for two well-instrumented field tests in multilayered soils is given to demonstrate the effectiveness and accuracy of the proposed approach. It can be noted that the presented solution can be used to calculate the settlement of single pile for the preliminary design of pile foundation.展开更多
Cone penetration test(CPT)is an appropriate technique for quickly determining the geotechnical properties of lunar soil,which is valuable for in situ lunar exploration.Utilizing a typical coupling method recently deve...Cone penetration test(CPT)is an appropriate technique for quickly determining the geotechnical properties of lunar soil,which is valuable for in situ lunar exploration.Utilizing a typical coupling method recently developed by the authors,a finite element method(FEM)-discrete element method(DEM)coupled model of CPTs is obtained.A series of CPTs in lunar soil are simulated to qualitatively reveal the flow of particles and the development of resistance throughout the penetration process.In addition,the effects of major factors,such as penetration velocity,penetration depth,cone tip angle,and the low gravity on the Moon surface are investigated.展开更多
The weak intercalated soils in redbed soft rocks of Badong formation have obvious creep characters. In order to predict the unsaturated creep behaviors of weak intercalated soils, an unsaturated creep model was establ...The weak intercalated soils in redbed soft rocks of Badong formation have obvious creep characters. In order to predict the unsaturated creep behaviors of weak intercalated soils, an unsaturated creep model was established based on the unsaturated creep tests of weak intercalated soils by using GDS triaxial apparatus. The results show that the creep behaviors of intercalated soils are apparent and significantly affected by matric suction. Based on this, an empirical Mesri creep model for intercalated soils under varying matric suctions was built. The fitting results show that the parameters Ed and m of this model are in good power relations with matric suction s and stress level Dr, respectively. An improved Mesri creep model was established involving stress-matric suction-strain-time, which is more precise than the Mesri creep model in predicting the unsaturated creep behaviors of weak intercalated soils.展开更多
In view of the cumbersome and often untimely process of manual collection and observation of frozen soil data parameters,and the damage caused to dams by frost heaving of frozen soil,a remote monitoring and an early w...In view of the cumbersome and often untimely process of manual collection and observation of frozen soil data parameters,and the damage caused to dams by frost heaving of frozen soil,a remote monitoring and an early warning model for frozen soil in dam areas was presented.The Pt100 temperature sensors and JM seam gauges were used as measurement tools in the system.The sensor layout was designed,based on the actual situation in the monitoring area.A 4G network was used for wireless transmission to monitor frozen soil data in real time.BP neural network was used to predict the parameters of frozen soil.After analysis,four factors including the average temperature of frozen soil,the type of frozen soil,the artificial upper limit of frozen soil and the building construction time were selected to establish an early warning model using fuzzy reasoning.The experimental results showed that the early warning model could reflect the influence on dam buildings of frost heaving and sinking of frozen soil,and provided technical support for predicting the hazard level.展开更多
This article investigates the characteristics of shock wave overpressure generated by multi-layer composite charge under different detonation modes.Combining dimensional analysis and the explosion mechanism of the cha...This article investigates the characteristics of shock wave overpressure generated by multi-layer composite charge under different detonation modes.Combining dimensional analysis and the explosion mechanism of the charge,a peak overpressure prediction model for the composite charge under singlepoint detonation and simultaneous detonation was established.The effects of the charge structure and initiation method on the overpressure field characteristics were investigated in AUTODYN simulation.The accuracy of the prediction model and the reliability of the numerical simulation method were subsequently verified in a series of static explosion experiments.The results reveal that the mass of the inner charge was the key factor determining the peak overpressure of the composite charge under single-point detonation.The peak overpressure in the radial direction improved apparently with an increase in the aspect ratio of the charge.The overpressure curves in the axial direction exhibited a multi-peak phenomenon,and the secondary peak overpressure even exceeded the primary peak at distances of 30D and 40D(where D is the charge diameter).The difference in peak overpressure among azimuth angles of 0-90°gradually decreased with an increase in the propagation distance of the shock wave.The coupled effect of the detonation energy of the inner and outer charge under simultaneous detonation improved the overpressure in both radial and axial directions.The difference in peak overpressure obtained from model prediction and experimental measurements was less than 16.4%.展开更多
The creep behaviors of granite residual soil with pre-stress of 100 kPa was investigated by a series of small size creep tests. Three different types of strain curves were obtained at different stress levels. Based on...The creep behaviors of granite residual soil with pre-stress of 100 kPa was investigated by a series of small size creep tests. Three different types of strain curves were obtained at different stress levels. Based on creep characteristics of the granite residual soil under different stress levels, a creep model of the granite residual soil was established by rheological theory, and related parameters of the model were determined according to the experimental data at the same time. Further on, based on the established creep model, a theoretical model of dynamic stress accumulation in the granite residual soil under cyclic loading was deduced. It is found that there is a threshold of dynamic stress accumulation in this theoretical model. The dynamic stress accumulation laws of the granite residual soil are different under different cyclic loading stress. Finally, with the dynamic stress accumulation laws in the small-size samples of granite residual soil under different cycle loading studied and the experimental results comparing with the theoretical results, it verifies the validity of the theoretical model.展开更多
With the high-quality development of urban buildings,higher requirements are come up with for lateral bearing capacity of laterally loaded piles.Consequently,a more accurate analysis to predict the lateral response of...With the high-quality development of urban buildings,higher requirements are come up with for lateral bearing capacity of laterally loaded piles.Consequently,a more accurate analysis to predict the lateral response of the pile within an allowable displacement is an important issue.However,the current p-y curve methods cannot fully take into account the pile-soil interaction,which will lead to a large calculation difference.In this paper,a new analytical p-y curve is established and a finite difference method for determining the lateral response of pile is proposed,which can consider the separation effect of pile-soil interface and the coefficient of circumferential friction resistance.In particular,an analytical expression is developed to determine the compressive soil pressure by dividing the compressive soil pressure into two parts:initial compressive soil pressure and increment of compressive soil pressure.In addition,the relationship between compressive soil pressure and horizontal displacement of the pile is established based on the reasonable assumption.The correctness of the proposed method is verified through four examples.Based on the verified method,a parametric analysis is also conducted to investigate the influences of factors on lateral response of the pile,including internal friction angle,pile length and elastic modulus of pile.展开更多
On the basis of Artificial Neural Network theory, a back propagation neural network with one middle layer is building in this paper, and its algorithms is also given, Using this BP network model, study the case of Mal...On the basis of Artificial Neural Network theory, a back propagation neural network with one middle layer is building in this paper, and its algorithms is also given, Using this BP network model, study the case of Malian-River basin. The results by calculating show that the solution based on BP algorithms are consis- tent with those based multiple - variables linear regression model. They also indicate that BP model in this paper is reasonable and BP algorithms are feasible.展开更多
The characteristics of a disturbed soil-structure interface were studied based on the variation regularities of the disturbed soil within its mining subsidence area using direct shear tests.The effects of the initial ...The characteristics of a disturbed soil-structure interface were studied based on the variation regularities of the disturbed soil within its mining subsidence area using direct shear tests.The effects of the initial moisture content on the shear strength parameters of the soil-structure interfaces were analyzed.The results indicate that the cohesion of the interface initially increased and then decreased as the initial moisture content increased.In addition,the friction angle of the interface decreased as the initial moisture content increased.A constitutive model of the disturbed soil-structure interface,a rigid-plastic model based on the initial void ratio and saturability(VSRP) model,was established based on the results.In order to validate this model,a finite element analysis of DRS-1 direct shear tests was conducted.The finite element model calculations coincided with the results of the DRS-1 direct shear tests.The proposed model also reflected the nonlinear features of the soil-structure interface.展开更多
土壤是具有高度异质性的复合体。早期的数字土壤制图研究主要关注水平方向的土壤空间变异和制图,对垂直方向空间变异和土壤三维制图考虑较少。近年来,三维地理信息技术和对地观测与探测技术的快速发展,极大地促进了土壤三维空间数据获...土壤是具有高度异质性的复合体。早期的数字土壤制图研究主要关注水平方向的土壤空间变异和制图,对垂直方向空间变异和土壤三维制图考虑较少。近年来,三维地理信息技术和对地观测与探测技术的快速发展,极大地促进了土壤三维空间数据获取、三维空间推测、三维数据模型、三维模型构建和可视化方法等方面的研究。本文对三维空间土壤推测与土壤模型构建的已有方法进行梳理和评述,以期为三维数字土壤制图的应用和发展提供建议。以三维土壤制图、三维GIS、三维数据模型、三维地质建模、三维可视化、土壤空间变异、空间推测、克里格插值、土壤-景观分析、深度函数、机器学习、地统计学、随机模拟等为关键词检索Web of Science数据库,基于相关度、引用率和文献来源等因素进一步筛选出重点文献进行分析。归纳整理了土壤空间变异性、三维空间土壤推测、三维空间数据模型和三维模型构建等关键技术的现有研究体系,对各种三维推测和建模方法的优缺点和适用场景作出评价。针对目前研究中存在的垂直方向土壤数据稀少、土壤三维推测精度低、三维模型质量待提高等问题,提出一些可行的研究思路。展开更多
基金Projects(42477162,52108347,52178371,52168046,52178321,52308383)supported by the National Natural Science Foundation of ChinaProjects(2023C03143,2022C01099,2024C01219,2022C03151)supported by the Zhejiang Key Research and Development Plan,China+6 种基金Project(LQ22E080010)supported by the Exploring Youth Project of Zhejiang Natural Science Foundation,ChinaProject(LR21E080005)supported by the Outstanding Youth Project of Natural Science Foundation of Zhejiang Province,ChinaProject(2022M712964)supported by the Postdoctoral Science Foundation of ChinaProject(2023AFB008)supported by the Natural Science Foundation of Hubei Province for Youth,ChinaProject(202203)supported by Engineering Research Centre of Rock-Soil Drilling&Excavation and Protection,Ministry of Education,ChinaProject(202305-2)supported by the Science and Technology Project of Zhejiang Provincial Communication Department,ChinaProject(2021K256)supported by the Construction Research Founds of Department of Housing and Urban-Rural Development of Zhejiang Province,China。
文摘This study proposes a general imperfect thermal contact model to predict the thermal contact resistance at the interface among multi-layered composite structures.Based on the Green-Lindsay(GL)thermoelastic theory,semi analytical solutions of temperature increment and displacement of multi-layered composite structures are obtained by using the Laplace transform method,upon which the effects of thermal resistance coefficient,partition coefficient,thermal conductivity ratio and heat capacity ratio on the responses are studied.The results show that the generalized imperfect thermal contact model can realistically describe the imperfect thermal contact problem.Accordingly,it may degenerate into other thermal contact models by adjusting the thermal resistance coefficient and partition coefficient.
基金Project(51874202) supported by the National Natural Science Foundation of ChinaProject(2017JQ0003) supported by the Sichuan Youth Fund,China。
文摘Considering the fact that in some complex cases,plate anchors are buried in multi-layered geotechnical materials,the ultimate dynamic analysis was performed to investigate the uplift capacity and failure mechanism of shallow strips and circular plate anchors in multi-layered soils.The nonlinear strength criterion and non-associated flow rule of geotechnical materials were introduced to investigate the influence of nonuniformity on the pullout performance and failure mechanism of shallow plate anchors.The expressions of the detaching curves or surfaces were obtained to reflect the failure mechanism,which can be used to figure out the ultimate uplift capacity and failure range.The results are generally in agreement with the numerical simulations and previous research.The effects of various parameters on the ultimate uplift capacity and failure mechanism of plate anchors in multi-layered soils were investigated,and it is found that the ultimate uplift capacity and failure range of shallow anchors increase with the increase of initial cohesion and dilatancy coefficient,but decrease with the unit weight,axial tensile stress and nonlinear coefficient.
基金Project(1301015A)supported by the Post-doctoral Research Fund of Jiangsu Province,ChinaProject Funded by the Priority Academic Program of Jiangsu Higher Education Institution,China+1 种基金Project(2014M561566)supported by China Postdoctoral Science FoundationProject(YK913007)supported by Key Laboratory of Earth-Rock Dam Failure Mechanism and Safety Control Technologies,China
文摘Two modifications for the basic Barcelona model(BBM) are present. One is the replacement of the net stress by the average skeleton stress in unsaturated soil modeling, and the other is the adoption of an expression for the load-collapse(LC) yield surface that can match flexibly the normal compression lines at different suctions. The predictions of the modified BBM for the controlled-suction triaxial test on the unsaturated compacted clay are presented and compared with the experimental results. A good agreement between the predicted and experimental results demonstrates the reasonability of the modified BBM. On this basis, the coupled processes of groundwater flow and soil deformation in a homogeneous soil slope under a long heavy rainfall are simulated with the proposed elasto-plastic model. The numerical results reveal that the failure of a slope under rainfall infiltration is due to both the reduction of soil suction and the significant rise in groundwater table. The evolution of the displacements is greatly related to the change of suction. The maximum collapse deformation happens near the surface of slope where infiltrated rainwater can quickly reach. The results may provide a helpful reference for hazard assessment and control of rainfall-induced landslides.
基金Project(50879077) supported by the National Natural Science Foundation of China
文摘A new approach is proposed to analyze the settlement behavior for single pile embedded in layered soils. Firstly, soil layers surrounding pile shaft are simulated by using distributed Voigt model, and finite soil layers under the pile end are assumed to be virtual soil-pile whose cross-section area is the same as that of the pile shaft. Then, by means of Laplace transform and impedance function transfer method to solve the static equilibrium equation of pile, the analytical solution of the displacement impedance fimction at the pile head is derived. Furthermore, the analytical solution of the settlement at the head of single pile is theoretically derived by virtue of convolution theorem. Based on these solutions, the influences of parameters of soil-pile system on the settlement behavior for single pile are analyzed. Also, comparison of the load-settlement response for two well-instrumented field tests in multilayered soils is given to demonstrate the effectiveness and accuracy of the proposed approach. It can be noted that the presented solution can be used to calculate the settlement of single pile for the preliminary design of pile foundation.
基金Project(51278451) supported by the National Natural Science Foundation of ChinaProject(LZ12E09001) supported by the Zhejiang Natural Science Foundation,China
文摘Cone penetration test(CPT)is an appropriate technique for quickly determining the geotechnical properties of lunar soil,which is valuable for in situ lunar exploration.Utilizing a typical coupling method recently developed by the authors,a finite element method(FEM)-discrete element method(DEM)coupled model of CPTs is obtained.A series of CPTs in lunar soil are simulated to qualitatively reveal the flow of particles and the development of resistance throughout the penetration process.In addition,the effects of major factors,such as penetration velocity,penetration depth,cone tip angle,and the low gravity on the Moon surface are investigated.
基金Project supported by Science&Technology Program of Hubei Traffic and Transport Office,ChinaProject(41272377)supported by the National Natural Science Foundation of China
文摘The weak intercalated soils in redbed soft rocks of Badong formation have obvious creep characters. In order to predict the unsaturated creep behaviors of weak intercalated soils, an unsaturated creep model was established based on the unsaturated creep tests of weak intercalated soils by using GDS triaxial apparatus. The results show that the creep behaviors of intercalated soils are apparent and significantly affected by matric suction. Based on this, an empirical Mesri creep model for intercalated soils under varying matric suctions was built. The fitting results show that the parameters Ed and m of this model are in good power relations with matric suction s and stress level Dr, respectively. An improved Mesri creep model was established involving stress-matric suction-strain-time, which is more precise than the Mesri creep model in predicting the unsaturated creep behaviors of weak intercalated soils.
基金Supported by the Application Technology Research and Development Plan Project of Heilongjiang Province(GY2014ZB0011)the 13th Five-year National Key R&D Program(2016YFD0300610)
文摘In view of the cumbersome and often untimely process of manual collection and observation of frozen soil data parameters,and the damage caused to dams by frost heaving of frozen soil,a remote monitoring and an early warning model for frozen soil in dam areas was presented.The Pt100 temperature sensors and JM seam gauges were used as measurement tools in the system.The sensor layout was designed,based on the actual situation in the monitoring area.A 4G network was used for wireless transmission to monitor frozen soil data in real time.BP neural network was used to predict the parameters of frozen soil.After analysis,four factors including the average temperature of frozen soil,the type of frozen soil,the artificial upper limit of frozen soil and the building construction time were selected to establish an early warning model using fuzzy reasoning.The experimental results showed that the early warning model could reflect the influence on dam buildings of frost heaving and sinking of frozen soil,and provided technical support for predicting the hazard level.
基金funded by the National Natural Science Foundation of China(Grant No.11972018,No.12002336)China Postdoctoral Science Foundation(Grant No.2021M701710)。
文摘This article investigates the characteristics of shock wave overpressure generated by multi-layer composite charge under different detonation modes.Combining dimensional analysis and the explosion mechanism of the charge,a peak overpressure prediction model for the composite charge under singlepoint detonation and simultaneous detonation was established.The effects of the charge structure and initiation method on the overpressure field characteristics were investigated in AUTODYN simulation.The accuracy of the prediction model and the reliability of the numerical simulation method were subsequently verified in a series of static explosion experiments.The results reveal that the mass of the inner charge was the key factor determining the peak overpressure of the composite charge under single-point detonation.The peak overpressure in the radial direction improved apparently with an increase in the aspect ratio of the charge.The overpressure curves in the axial direction exhibited a multi-peak phenomenon,and the secondary peak overpressure even exceeded the primary peak at distances of 30D and 40D(where D is the charge diameter).The difference in peak overpressure among azimuth angles of 0-90°gradually decreased with an increase in the propagation distance of the shock wave.The coupled effect of the detonation energy of the inner and outer charge under simultaneous detonation improved the overpressure in both radial and axial directions.The difference in peak overpressure obtained from model prediction and experimental measurements was less than 16.4%.
基金Projects(41572277,41877229)supported by the National Natural Science Foundation of ChinaProject(2018B030311066)supported by the Natural Science Foundation of Guangdong Province,ChinaProject(201607010023)supported by the Science and Technology Program of Guangzhou,China
文摘The creep behaviors of granite residual soil with pre-stress of 100 kPa was investigated by a series of small size creep tests. Three different types of strain curves were obtained at different stress levels. Based on creep characteristics of the granite residual soil under different stress levels, a creep model of the granite residual soil was established by rheological theory, and related parameters of the model were determined according to the experimental data at the same time. Further on, based on the established creep model, a theoretical model of dynamic stress accumulation in the granite residual soil under cyclic loading was deduced. It is found that there is a threshold of dynamic stress accumulation in this theoretical model. The dynamic stress accumulation laws of the granite residual soil are different under different cyclic loading stress. Finally, with the dynamic stress accumulation laws in the small-size samples of granite residual soil under different cycle loading studied and the experimental results comparing with the theoretical results, it verifies the validity of the theoretical model.
基金Project(52068004)supported by the National Natural Science Foundation of ChinaProject(2018JJA160134)supported by the Natural Science Foundation of Guangxi Province,ChinaProject(AB19245018)supported by Key Research Projects of Guangxi Province,China。
文摘With the high-quality development of urban buildings,higher requirements are come up with for lateral bearing capacity of laterally loaded piles.Consequently,a more accurate analysis to predict the lateral response of the pile within an allowable displacement is an important issue.However,the current p-y curve methods cannot fully take into account the pile-soil interaction,which will lead to a large calculation difference.In this paper,a new analytical p-y curve is established and a finite difference method for determining the lateral response of pile is proposed,which can consider the separation effect of pile-soil interface and the coefficient of circumferential friction resistance.In particular,an analytical expression is developed to determine the compressive soil pressure by dividing the compressive soil pressure into two parts:initial compressive soil pressure and increment of compressive soil pressure.In addition,the relationship between compressive soil pressure and horizontal displacement of the pile is established based on the reasonable assumption.The correctness of the proposed method is verified through four examples.Based on the verified method,a parametric analysis is also conducted to investigate the influences of factors on lateral response of the pile,including internal friction angle,pile length and elastic modulus of pile.
基金Supported by Brilliant Youth Fund in Hebei Province
文摘On the basis of Artificial Neural Network theory, a back propagation neural network with one middle layer is building in this paper, and its algorithms is also given, Using this BP network model, study the case of Malian-River basin. The results by calculating show that the solution based on BP algorithms are consis- tent with those based multiple - variables linear regression model. They also indicate that BP model in this paper is reasonable and BP algorithms are feasible.
基金Project(51274192)supported by National Natural Science Foundation of China
文摘The characteristics of a disturbed soil-structure interface were studied based on the variation regularities of the disturbed soil within its mining subsidence area using direct shear tests.The effects of the initial moisture content on the shear strength parameters of the soil-structure interfaces were analyzed.The results indicate that the cohesion of the interface initially increased and then decreased as the initial moisture content increased.In addition,the friction angle of the interface decreased as the initial moisture content increased.A constitutive model of the disturbed soil-structure interface,a rigid-plastic model based on the initial void ratio and saturability(VSRP) model,was established based on the results.In order to validate this model,a finite element analysis of DRS-1 direct shear tests was conducted.The finite element model calculations coincided with the results of the DRS-1 direct shear tests.The proposed model also reflected the nonlinear features of the soil-structure interface.
文摘土壤是具有高度异质性的复合体。早期的数字土壤制图研究主要关注水平方向的土壤空间变异和制图,对垂直方向空间变异和土壤三维制图考虑较少。近年来,三维地理信息技术和对地观测与探测技术的快速发展,极大地促进了土壤三维空间数据获取、三维空间推测、三维数据模型、三维模型构建和可视化方法等方面的研究。本文对三维空间土壤推测与土壤模型构建的已有方法进行梳理和评述,以期为三维数字土壤制图的应用和发展提供建议。以三维土壤制图、三维GIS、三维数据模型、三维地质建模、三维可视化、土壤空间变异、空间推测、克里格插值、土壤-景观分析、深度函数、机器学习、地统计学、随机模拟等为关键词检索Web of Science数据库,基于相关度、引用率和文献来源等因素进一步筛选出重点文献进行分析。归纳整理了土壤空间变异性、三维空间土壤推测、三维空间数据模型和三维模型构建等关键技术的现有研究体系,对各种三维推测和建模方法的优缺点和适用场景作出评价。针对目前研究中存在的垂直方向土壤数据稀少、土壤三维推测精度低、三维模型质量待提高等问题,提出一些可行的研究思路。