FFP在黏土中贯入过程的CFD模拟

圆锥静力触探法(cone penetration test,CPT)广泛应用于土的抗剪强度测试中.在此基础上发展的自由落体式贯入仪(free fall penetrometer,FFP),依靠自由下落获得的动能和自身重力势能贯入土中,不需要借助外部加载装置,提高其使用的便捷性.但在动力贯入过程中FFP与土的相互作用更为复杂,涉及到土体的率效应和拖曳阻力等的影响.因此,对FFP各项受力和相关参数的准确分析有助于提高其实用性以及测量的准确性.作者采用基于计算流体动力学(computational fluid dynamics,CFD)的分析软件ANSYS CFX 17.0模拟FFP在均质黏土中的贯入过程,借助动网格的大变形分析方法模拟FFP贯入过程中的运动边界问题.提出了薄层单元法模拟FFP与土体的界面摩擦接触行为.在CFD模拟中,土体材料采用非牛顿流体来模拟,其剪切强度受土体切应变率的影响(即土体的率效应).通过模拟贯入仪以不同的速度在黏土中贯入的过程,研究FFP的端部阻力和侧壁阻力与贯入速率、土体强度和密度、界面摩擦系数以及率效应参数之间的关系,建立了端部承载力系数、端部和侧壁率效应参数及拖曳系数的表达式,并提出了土体不排水抗剪强度的预测方法,为FFP测试数据的解析提供依据.     

基于CEL有限元法的大管径钢管桩贯入海基过程的承载力分析

鉴于耦合欧拉-拉格朗日(Coupled Eulerian-Lagrangian, CEL)有限元方法在计算分析接触大变形问题时的优势,建立了大管径钢管桩贯入海基的CEL有限元模型,即在被贯入的海基层采用欧拉有限元法建模,其它层采用拉格朗日有限元法建模,两者之间通过罚函数的接触算法进行耦合．在欧拉有限元模型中,采用算子分裂技术．为兼顾计算模型的精度和效率,本文先对模型的网格密度和贯入速度进行稳定性分析,以此确定合适的网格模型和适宜的贯入速度,同时验证CEL方法在贯入模拟中的可行性．随后针对不同土体类型和不同土体参数变化,探究土体贯入阻力随贯入深度的变化规律．计算结果表明砂土的贯入阻力最大,粉砂土次之,粘土的贯入阻力最小;且贯入阻力对土体内摩擦角的变化最为敏感,对土体粘聚力的变化比较敏感,而对土体模量的变化不敏感．     

土体的本构模型和超重力物理模拟

数值模拟和物理模拟是分析土体沉降和稳定性的主要手段. 本构模型作为描述土体应力应变关系的数学表达式, 是数值模拟的基础. 土体具有碎散性, 这一基本物理特性导致了其具有压硬性、摩擦性和剪胀性, 这是土的力学特性区别于金属的主要特征, 在土体的本构模型中必须反映这3个基本特性. 传统土力学将土体的变形和强度分离考虑, 分别采用弹性理论和基于刚塑性模型的极限平衡理论分析, 虽然应用广泛, 但由于不能全面地反映土的基本力学特性, 计算结果的精度常常难以满足定量分析的需要. 剑桥模型作为第一个全面反映压硬性、摩擦性和剪胀性的弹塑性本构模型, 实现了变形和强度的统一, 能较好地描述饱和正常固结黏土的应力应变关系, 被视为是现代土力学的开端; 统一硬化模型通过引入一个独特的硬化参数进一步发展了剑桥模型, 将适用范围扩大到超固结黏土. 作者认为, 未来岩土体本构模型研究的挑战是: 如何考虑岩土体在受力过程中土骨架相变与多场耦合, 以解决目前本构模型尚无法定量分析的能源、交通、环境和水利相关的重大岩土工程问题. 超重力物理模拟具有缩尺效应和缩时效应, 克服了常重力物理模拟中模型的应力水平低于原型的缺点, 特别适用于大尺度、长历时问题的模拟. 相较数值模拟, 超重力物理模拟的优势在于能够检验本构模型的合理性, 揭示本构模型无法描述的未知特性. 最后, 介绍了采用数值模拟和物理模拟联合分析大直径钢管桩水平受荷特性的工程案例.      

CONSTITUTIVE MODELS AND HYPERGRAVITY PHYSICAL SIMULATION OF SOILS$^{\bf 1)}$

数值模拟和物理模拟是分析土体沉降和稳定性的主要手段. 本构模型作为描述土体应力应变关系的数学表达式, 是数值模拟的基础. 土体具有碎散性, 这一基本物理特性导致了其具有压硬性、摩擦性和剪胀性, 这是土的力学特性区别于金属的主要特征, 在土体的本构模型中必须反映这3个基本特性. 传统土力学将土体的变形和强度分离考虑, 分别采用弹性理论和基于刚塑性模型的极限平衡理论分析, 虽然应用广泛, 但由于不能全面地反映土的基本力学特性, 计算结果的精度常常难以满足定量分析的需要. 剑桥模型作为第一个全面反映压硬性、摩擦性和剪胀性的弹塑性本构模型, 实现了变形和强度的统一, 能较好地描述饱和正常固结黏土的应力应变关系, 被视为是现代土力学的开端; 统一硬化模型通过引入一个独特的硬化参数进一步发展了剑桥模型, 将适用范围扩大到超固结黏土. 作者认为, 未来岩土体本构模型研究的挑战是: 如何考虑岩土体在受力过程中土骨架相变与多场耦合, 以解决目前本构模型尚无法定量分析的能源、交通、环境和水利相关的重大岩土工程问题. 超重力物理模拟具有缩尺效应和缩时效应, 克服了常重力物理模拟中模型的应力水平低于原型的缺点, 特别适用于大尺度、长历时问题的模拟. 相较数值模拟, 超重力物理模拟的优势在于能够检验本构模型的合理性, 揭示本构模型无法描述的未知特性. 最后, 介绍了采用数值模拟和物理模拟联合分析大直径钢管桩水平受荷特性的工程案例.     

饱和黏土中考虑土塞效应的柱孔扩张分析

为研究静压开口管桩的挤土效应,在考虑土塞效应的前提下,采用柱孔扩张理论对管桩压入饱和软黏土时的扩孔过程进行了分析。基于可以合理描述黏性土强度特性的拓展Lade-Duncan屈服准则,对柱孔扩张过程中的桩周土体进行了弹塑性理论推导,得到了桩周土体应力场、位移场、塑性区半径、极限扩孔压力和塑性区外侧边界径向位移的解析解。在此基础上,通过室内模型实验获得了软土中开口管桩的土塞高度,且运用等效替代法验证了理论解答的合理性;并通过引入土塞长度比对参数进行了分析。结果表明:在沉桩过程中随着土塞长度比的增大,塑性区半径、极限扩孔压力和塑性区外侧边界径向位移均逐渐减小;土体的剪切模量和黏聚力对沉桩挤土效应有显著的影响,挤土效应始终随着土塞效应的增强而减弱。该研究成果对开口管桩工程的实际应用具有一定理论指导意义。     

被动桩中土拱效应特征与影响参数研究

采用平面有限元方法,对被动桩在粘性、无粘性土体条件下的土拱形成机理进行了分析。从土拱效应的4个主要方面:土拱形状、桩周土体塑性(拉裂)区的分布、竖向位移等值线及桩后土体残余荷载分担比,探讨了改变桩、土参数,如粘聚力、内摩擦角、弹性模量、泊松比、剪胀角、桩土接触面特性时,土拱效应的变化规律。研究表明,土性不同时,土拱形成过程中土体的受力和破坏形式不同;泊松比、剪胀角、桩土接触特征对土拱效应的影响最明显;强度高、剪胀角大、桩土接触面粗糙或泊松比小的土体,易形成土拱。根据研究结果,指出Tom io Ito排桩塑性绕流土压力理论计算公式尚存不足。     

基于van Genuchten模型的非饱和土非线性抗剪强度研究

为了更直观和全面地揭示非饱和土的抗剪强度变化规律,针对土体剪切破坏面上抗剪强度与法向应力呈非线性关系以及基质吸力关联摩擦角随基质吸力改变呈非线性变化这两种实验现象,结合vanGenuchten土水特征曲线模型,提出了非饱和土的非线性抗剪强度包络壳模型,进而推导得出了基于土水特征曲线(SWCC)相关参数的非饱和土非线性抗剪强度计算公式。研究表明:模型中各参数较易获取且无须增加额外的试验;与Fredlund包络(平)面模型相比,本文模型参数获取过程更为客观,可有效保证计算的唯一性和准确性;通过与已有文献所给试验结果进行对比分析,本文模型的适用性与可靠性得到了验证,相较于著名的Fredlund非饱和土抗剪强度计算公式,由本文方法给出的非线性抗剪强度与土体真实抗剪强度更为接近。     

饱和超固结黏性土的三剪弹塑性本构模型研究

针对饱和超固结黏性土现有下加载面修正剑桥模型中破坏应力比为定值、土体黏聚力为零,以及不能准确反映不同应力状态下土的强度差异这些问题,基于三剪统一强度准则以及应力坐标平移法得到了扩展破坏应力比,其特点是能更好地反映应力状态变化以及土体黏聚力的影响。在此基础上提出了饱和超固结黏性土的三剪弹塑性本构模型,该模型的特点是能描述土体受力时的中间主应力效应,应力区间效应和拉压差影响,同时也能更好地考虑土体黏聚力的影响。基于该模型对ABAQUS软件进行了二次开发,并利用其模拟了饱和超固结黏性土在排水和不排水条件下的真三轴和常规三轴压缩试验特性。对常规三轴压缩条件下土体力学特性作了模拟和试验结果对比。结果表明所提模型能很好地反映不同超固结比下土体的变形、剪胀、孔隙水压力变化特性。     

含石量对土石混合体压桩承载力影响的离散元分析

土石混合体是介于土体和岩体之间的一种非均质、非连续和非线性的特殊工程地质材料,其在压桩贯入过程中的承载力受含石量的影响非常显著。本文分别采用球形颗粒和非规则镶嵌组合颗粒模拟土体颗粒和块石,对不同含石量下压桩贯入过程进行离散元数值分析。计算结果表明,桩柱阻力及其波动规律在不同含石量下有很大的差别。高含石量下的阻力要明显大于低含石量下的阻力,且其波动性也更加明显。通过对土石混合体内部力链结构的细观分析,揭示了压桩贯入过程中承载力随含石量变化的内在机理。以上研究有助于分析土石混合体材料的宏观力学行为,深入研究其在复杂工程条件下的力学特性。     

柠条、沙柳根与土及土与土界面摩擦特性

通过直剪试验研究了垂直荷载、土体干密度、土体含水率对柠条根-土界面、沙柳根-土界面以及土-土界面摩擦特性的影响.结果表明:柠条根-土界面、沙柳根-土界面与土-土界面的剪切特性相似,均服从莫尔-库仑理论;柠条根-土界面、沙柳根-土界面的摩擦系数随着土体干密度的增大而增大;土壤含水率对柠条根-土界面、沙柳根-土界面摩擦系数和黏聚力的影响与素土相同,三者的摩擦系数随含水率增加逐渐减小,而黏聚力随着含水率增大呈先增大后减少的变化.在自然生境下主要根系分布层,柠条、沙柳的根-土界面摩擦系数大于土-土界面,尽管根-土界面黏聚力小于土-土界面黏聚力,但是根-土界面的摩擦系数的增大仍使根-土界面抗剪强度大于土-土界面,说明根-土结合面抵抗位移的能力大于土-土界面;2种植物相比,柠条根系抵抗根-土分离的能力优于沙柳.     

Evaluation of penetrometers for measuring soil strength

Cone index, as determined by a cone penetrometer, is frequently used as a measure of soil strength. The index is a compound parameter involving components of shear, compressive and tensile strength and soil metal friction. In order to assess the effect of soil type and condition on the relative contributions of these components to penetration resistance, the forces required to push blunt and sharp probes into two soils under a range of moisture contents and bulk densities were investigated. The maximum penetration force in homogeneous soil was not uniquely related to dry bulk density or cohesion, but varied with soil moisture content.At high and low moisture contents, the soil tended to interact with the shaft of the penetrometer thus increasing the resistance to penetration. At low moisture content, bodies of compressed soil formed in front of the probe, effectively changing the probe geometry.It was concluded that interpretation of cone index in typical layered field soils is difficult. Even in homogeneous soils, the proportion of shear, compressive and tensile components that the cone index reflects varies with soil condition.     

Cone penetration resistance equation as a function of the clay ratio, soil moisture content and specific weight

The cone penetrometer is a simple versatile device which is widely used to monitor the strength of a soil in terms of its resistance to the penetration of a standard cone. The soil penetration resistance is a function of soil moisture content, soil specific weight and soil type. The soil type is characterised by means of a clay ratio which is the ratio of the clay content of the soil to the content of silt and sand.Based on the classical bearing capacity theories for strip foundations, a general cone penetration resistance equation is developed to represent the variability of cohesion and friction angle by means of soil type and moisture content. The empirical relationship is shown to give an accurate prediction of the cone penetration resistance for a wide range of soils from a loamy sand to a heavy clay (clay ratios 0.10–1.60) and over a wide spectrum of soil moisture contents from 10 to 65% w/w.     

Non-linear finite element analysis of cone penetration in layered sandy loam soil – Considering precompression stress state

Axisymmetric finite element (FE) method was developed to simulate cone penetration process in layered granular soil. The FE was modeled using ABAQUS/Explicit, a commercially available package. Soil was considered as a non-linear elastic plastic material which was modeled using variable elastic parameters of Young’s Modulus and Poisson’s ratio and Drucker–Prager criterion with yield stress dependent material hardening property. The material hardening parameters of the model were estimated from the USDA-ARS National Soil Dynamics Laboratory – Auburn University (NSDL-AU) soil compaction model. The stress–strain relationship in the NSDLAU compaction model was modified to account for the different soil moisture conditions and the influence of precompression stress states of the soil layers. A surface contact pair (‘slave-master’) algorithm in ABAQUS/Explicit was used to simulate the insertion of a rigid cone (RAX2 ABAQUS element) into deformable and layered soil medium (CAX4R ABAQUS element). The FE formulation was verified using cone penetration data collected on a soil chamber of Norfolk sandy loam soil which was prepared in two compaction treatments that varied in bulk density in the hardpan layer of (1) 1.64 Mg m⁻³ and (2) 1.71 Mg m⁻³. The FE model successfully simulated the trend of cone penetration in layered soils indicating the location of the sub-soil compacted (hardpan) layer and peak cone penetration resistance. Modification of the NSDL-AU model to account for the actual soil moisture content and inclusion of the influence of precompression stress into the strain behavior of the NSDL-AU model improved the performance of FE in predicting the peak cone penetration resistance. Modification of the NSDL-AU model resulted in an improvement of about 42% in the finite element-predicted soil cone penetration forces compared with the FE results that used the NSDL-AU ‘virgin’ model.     

A static cone penetrometer for ocean sediments: T.J. Hirst, A.F. Richards and A.L. InderbitzenASTM Spec. Techn. Publ. No. STP 501, 69–80 (1972)

A static cone penetrometer for deep ocean exploration has been developed and successfully tested in soft sea floor sediments at an ocean depth in excess of 1200 m. The penetration resistance of a standard 60-deg cone is detected by appropriately insulated strain gauges whose output may be monitored remotely from the sensing unit. Initial deployment of the device was aboard the submersible Deep Quest in the San Diego Trough off the coast of southern California. Resistances as low as 7 N (1·6 lb) were measured near the surface of the sediment, generally increasing with depth. At the maximum depth of penetration of 1.1 m (restricted by the limited negative buoyancy of the submersible) tip resistance reached 80 N (18 lb). By means of a conventional bearing capacity analysis, the cone resistance was converted to in situ strength and compared to vane shear measurements taken in the immediate vicinity of the penetration tests. Reasonable agreement between vane shear and penetration strengths were obtained for an assumed bearing capacity factor N_c = 9.     

Effect of wetting and drying on soil physical properties

Agricultural soils are subject to seasonal wetting and drying cycles. Effect of drying stress, as influenced by one cycle of wetting and drying, on physical properties of a clay–loam soil was investigated in the laboratory. The physical properties studied were soil bulk density, cone penetration resistance, shear strength, adhesion and aggregate size and stability. Three drying stress treatments were made by wetting air-dried soil of initial moisture content of 12% (on dry weight basis) to three different higher moisture contents, namely 27, 33 and 40%, and then drying each of them back to their original moisture content of 12%. Thus, the soil was subjected to three different degrees of drying stress. The results showed that the soil strength indicated by cone penetration resistance and cohesion, and soil aggregate size, increased with the degree of drying stress. However, the soil bulk density did not change significantly with the drying stress.     

静力触探锥头阻力的近似理论与实验研究进展

锥头阻力在静力触探试验中扮演着十分重要的角色.从不同角度,对触探中锥头阻力的研究进行简要阐述,对承载力理论、空洞膨胀理论、应变路径法及运动点位错法等几种理论分析方法进行了回顾.另外,对数值分析和实验研究的进展情况进行了叙述.并对各种方法的适用性进行了比较.承载力理论虽然简单,但忽略了土的压缩性和探杆周围初始应力的增加,所以不能精确地模拟锥头的深层贯入.空洞膨胀理论提供了一个分析锥头阻力的简单而较精确的方法,它考虑了土的压缩性(或膨胀性)和锥头贯入过程中锥杆周围应力增加的影响.但这种方法是将锥头贯入与空洞膨胀之间做了一个等效模拟,所以不同的模拟方法,得到的结果差别较大.应变路径法能够有效解决饱和粘土中的不排水贯入,但不适用于砂土.运动点位错法因为考虑了部分排水,所以能较好地预测固结系数,但采用了线弹性分析,故位错法在其他方面的应用还需要大量的试验验证.有限元法在处理锥头贯入这类慢侵彻问题时缺乏一种很好的处理技术,导致它在进行破坏荷载计算时有显著的误差和数值计算困难.标定槽试验将在验证和建立锥头阻力与土的性能关系方面继续起到一个重要作用,但其结果需经过校正后才可应用到现场.最后对该领域的研究趋势进行了讨论.      

Cone resistance in sand by incremental method

A new method has been developed for the determination of cone resistance under drained conditions. Numerical methods are used for the solution of the differential equations of plasticity theory for soils and for the determination of the stress states in the soil produced by the penetration of the cone. It is assumed that the stresses produced by the penetration of the cone remain ‘locked in’ the soil and constitute boundary conditions for further penetration. The computation starts with the cone base at the surface and is continued by successively incrementing the depth by a small amount. Charts are given for the computation of cone resistance in sands for various friction angles. The importance of the effect of the shear stresses generated at the surface of the cone and characterized by the interface friction angle, δ, is discussed in detail.     

扁铲侧胀试验计算地基承载力

简单介绍扁铲侧胀试验,首次提出用扁铲侧胀试验计算地基承载力的方法。利用上海地区部分重大工程中的扁铲 侧胀试验数据,在与室内土工试验、静力触探试验及十字板剪切试验的试验参数相关分析的基础上,经数理统计后给出了上 海地区用扁铲侧胀试验计算地基承载力的经验公式,并对其计算精度进行了评价。     

排水粉喷桩(2D工法)加固软土地基对桩周土体强度的影响

排水粉喷桩(简称2D工法)作为一种新的软土地基处理方法,在淮盐公路软基处理中得到试用,效果良好;但是目前对其桩间土强度的变化规律还不是很清楚,本文首次通过大量现场试验,采用静力触探、十字板剪切等试验手段,测得了排水粉喷桩施工后桩间土强度变化的规律。结果显示在排水粉喷桩施工后,桩间土强度先降低,之后随着龄期的增长而提高。     

Micromechanical modeling for behavior of silty sand with influence of fine content

Silty sand is a soil mixture with coarse grains and fine grains. Experimental observations have shown that small amount of fines may reduce the undrained shear strength significantly. The purpose of this paper is to propose a micromechanical model for the stress–strain behavior of silty sand influenced by fines under drained and undrained conditions. The micromechanical stress–strain model accounts for the influence of fines on the density state of the soil mixture, thus consequently affect the critical state friction angle and the amount of sliding between particles. The present model is examined by simulating typical drained and undrained tests in conventional triaxial conditions. The simulated stress–strain curves are compared with the measured results on samples made of Ottawa sand and Foundry sand with various amounts of fines. The predictive ability of the present model for simulating the behavior of silty sand is discussed.