基于动态空腔膨胀模型的土壤Mohr-Coulomb参数反演方法

基于动态空腔膨胀模型的土壤Mohr-Coulomb参数反演方法 Inversion method of soil Mohr-Coulomb parameters based on dynamic cavity expansion model 1 2 first-author 邹 欧（1998-），男，湖南娄底人，硕士研究生，2020年于中国矿业大学（北京）获得学士学位，主要从事车辆地面力学相关研究。E-mail： zouou20@mails.ucas.ac.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41464330&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465404&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465399&type= 22.01333237 32.00399780 邹 欧 （1998-），男，湖南娄底人，硕士研究生，2020年于中国矿业大学（北京）获得学士学位，主要从事车辆地面力学相关研究。E-mail： zouou20@mails.ucas.ac.cn zouou20@mails.ucas.ac.cn 3 2 2 3 2 corresp E-mail： liyuqiong@imech.ac.cn E-mail： liyuqiong@imech.ac.cn 李玉琼（1982-），男，湖南娄底人，研究员，2010年于北京理工大学获得博士学位，主要从事车辆地面力学，星壤力学探测与模型构建研究。E-mail：liyuqiong@imech.ac.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465428&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465436&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465432&type= 22.01333618 32.00399780 李玉琼 （1982-），男，湖南娄底人，研究员，2010年于北京理工大学获得博士学位，主要从事车辆地面力学，星壤力学探测与模型构建研究。E-mail： liyuqiong@imech.ac.cn liyuqiong@imech.ac.cn 中国科学院大学 工程科学学院，北京 100049 中国科学院大学 工程科学学院，北京 100049 College of Engineering Science， University of Chinese Academy of Sciences， Beijing 100049， China College of Engineering Science， University of Chinese Academy of Sciences， Beijing 100049， China 中国科学院 力学研究所 流固耦合系统力学重点实验室，北京 100190 中国科学院 力学研究所 流固耦合系统力学重点实验室，北京 100190 Key Laboratory for Mechanics in Fluid Solid Coupling Systems， Institute of Mechanics， Chinese Academy of Sciences， Beijing 100190， China Key Laboratory for Mechanics in Fluid Solid Coupling Systems， Institute of Mechanics， Chinese Academy of Sciences， Beijing 100190， China 中国科学院 力学研究所 非线性力学国家重点实验室，北京 100190 中国科学院 力学研究所 非线性力学国家重点实验室，北京 100190 State Key Laboratory of Nonlinear Mechanics， Institute of Mechanics， Chinese Academy of Sciences， Beijing 100190， China State Key Laboratory of Nonlinear Mechanics， Institute of Mechanics， Chinese Academy of Sciences， Beijing 100190， China 为实现土壤原位力学参数的精准、快速解译，定量评估军事装备机动性，本文利用冲击贯入过程中测得的阻力-速度曲线、加速度时程曲线、速度时程曲线以及贯入深度-初速度关系，基于动态空腔膨胀模型建立了土壤Mohr-Coulomb参数的反演方法。本文分析表明准确拟合贯入阻力系数是反演Mohr-Coulomb参数的关键，而利用冲击贯入过程中的速度时程曲线进行拟合能够对Mohr-Coulomb参数进行准确的反演。针对Forrestal等人报道的实验结果，利用本文方法对土壤内聚力进行反演的相对误差为2.14%，对土壤内摩擦角反演的相对误差为9.77%。本文进一步给出了动态空腔膨胀模型下Mohr-Coulomb参数反演的可解域，讨论了参数敏感性。本文提出的反演方法突破了传统半经验解译方法中物理图像不清晰、经验参数依赖性强的问题，可为复杂地质环境下土壤Mohr-Coulomb参数快速确定及土壤承载力评估提供一个新的方法和途径。 To achieve accurate and rapid interpretation of in-situ soil mechanical parameters and quantitatively evaluate military equipment mobility， an inversion method for soil Mohr-Coulomb parameters is established based on the dynamic cavity expansion model using the resistance-velocity relationship， acceleration and velocity curves， and the relationship between depth and initial velocity measured during the impact penetration process. It indicates that accurate fitting of the resistance coefficients is critical for interpreting Mohr-Coulomb parameters. Precise results could be derived using the measured velocity curve. Based on the experimental results reported by Forrestal et al . using the proposed method， the relative errors of soil cohesion and internal friction angle interpretation are 2.14% and 9.77%， respectively. Furthermore， the solution domain of the Mohr-Coulomb parameters under the dynamic cavity expansion model is revealed， and parameter sensitivity is investigated. The proposed inversion method resolves the problems of unclear physical images and strong parameter dependency in traditional semi-empirical interpretation methods. It could provide a new approach for rapidly determining soil Mohr-Coulomb parameters and assessing soil bearing capacity in complex geological environments. 动态空腔膨胀模型 土壤强度 Mohr-Coulomb准则 参数反演 dynamic cavity expansion model soil strength Mohr-Coulomb criterion parametric inversion 1 引 言 1 土壤强度参数解译在定量评估军事装备机动性和指导装备轻量化结构设计等方面具有重要作用，在地质勘测、装备全域机动、地外行星科考等重大需求领域中有广泛的应用背景 ［ 1 1 - 5 5 1-5 ］ 。传统“取样-室内实验”是典型的非原位测试方法，取样造成的土体扰动——特别是原位应力的解除——使得所采集试样的力学性质不能真实反映土壤原本的力学特征，从而影响对土壤承载性能及装备车辆通行性的可靠评估 ［ 6 6 ］ 。因此，基于现场原位贯入触探测试并结合相关理论实现土壤原位力学参数的精准、快速反演的需求日益紧迫。目前，土壤原位贯入测试方法主要有静力触探（Cone Penetration Test，CPT）和自由落体贯入触探。静力触探能获得土壤强度参数随深度的变化情况，一般需要外载荷持续作用以保持匀速贯入（2.5 cm/s），设备安装耗时长，操作费用高 ［ 7 7 - 8 8 7-8 ］ ，难以满足复杂地质环境和大范围区域内土壤力学参数的测量，限制了其在装备机动性评估方面的应用。自由落体贯入触探技术是在静力触探基础上发展起来的一种新型原位测试技术，它依靠自由落体贯入仪（Free Fall Penetrometer，FFP）以一定初始速度贯入土壤，通过传感器实时记录贯入过程中贯入仪的加速度和受到的阻力，来反演土壤强度参数。FFP结构轻便、操作简单且较为经济，能够实现复杂环境和较大区域内土壤强度参数的快速测量 ［ 7 7 ， 9 9 ］ 。 目前，在海洋工程中，已有相关研究利用FFP获取土壤强度参数。Aubeny等人 ［ 10 10 ］ 提供了基于数值模拟建立的不排水抗剪强度计算方法，给出了承载力系数的经验公式；Stoll等人 ［ 11 11 ］ 通过对贯入实验所得数据进行最小二乘拟合，得到了贯入阻力和最大减速度之间的关系。研究表明贯入过程中的阻力可表示为一个与贯入速度有关的阻力函数 ［ 12 12 - 14 14 12-14 ］ 。通过确定土壤强度参数与阻力函数之间的关系，结合贯入仪上传感器记录的加速度数据可反演土壤强度参数。True ［ 15 15 ］ 给出了一个描述FFP贯入海床土过程的半经验阻力函数，其中考虑了端部阻力、侧壁摩擦阻力、土壤拖曳阻力以及土体浮力的作用。Chow等人 ［ 9 9 ］ 利用离心机实验中传感器测量的加速度数据，反演了饱和软黏土的不排水抗剪强度，并考虑了率效应对软黏土强度的影响。刘君等人 ［ 16 16 ］ 采用数值仿真的方法模拟了贯入仪在软黏土中的运动，研究其所受阻力与贯入速度、土壤强度和界面摩擦系数之间的关系，提出了一种反演软黏土不排水抗剪强度的方法。以上反演方法均基于True提出的半经验阻力函数，其中包括承载力系数，界面摩擦系数，拖曳阻力系数等许多经验参数，当缺少场地土壤条件前期资料时，经验参数取值会对反演方法的准确性产生重要影响。此外，以上方法仅能针对海洋工程中常见饱和软黏土的不排水抗剪强度进行反演。这表明，由于实际地质环境的复杂性，基于自由落体贯入技术的冲击贯入机理还有待更深入的研究，参数反演方法还需完善，其适用性还需拓展。 在贯入阻力的理论研究方面，现有方法主要包括承载力理论 ［ 17 17 ］ ，应变路径理论 ［ 18 18 ］ 和空腔膨胀理论 ［ 19 19 - 26 26 19-26 ］ 。承载力理论仅考虑土体破坏区域的整体平衡，而破坏区域内各点并不处于平衡状态。Janbu等人 ［ 27 27 ］ 通过黏性土的静力触探实验表明，利用承载力理论求得的贯入阻力相比于实测值低约40%。应变路径法只能近似求解饱和黏性土的情况 ［ 18 18 ］ ，对于砂土，由于贯入时的流场难以构造，该方法并不适用。 空腔膨胀理论是研究侵彻问题的重要模型。由于FFP的刚度远大于土壤，在贯入过程中变形较小，因此FFP贯入土壤的过程可视作刚性弹侵彻靶体的过程。若将贯入仪看作由一系列与其内切的球形空腔组成，则贯入仪与土体之间的相互作用就可看作这些球形空腔在无限大介质中均匀膨胀时与土体介质间的相互挤压 ［ 28 28 ］ 。通过求解介质中空腔膨胀的径向应力分布，就可得到弹体贯入过程中受到的贯入阻力。Bishop ［ 19 19 ］ 首先给出了准静态柱/球形空腔膨胀理论的控制方程，并利用它求解了尖锥体对金属靶的压入阻力问题。Chadwick ［ 20 20 ］ 将这个理论应用到理想土壤材料，结合Mohr-Coulomb屈服准则，得到了土壤径向位移的表达式。Cao等人 ［ 21 21 ］ 假定土的本构为修正剑桥模型，通过求解准静态空腔膨胀理论控制方程，得到了土体内的超孔隙水压力和总应力。Hill ［ 22 22 ］ 则首先给出动态空腔膨胀理论控制方程，得到了不可压缩理想弹塑性靶体对弹体的阻力。基于动态空腔膨胀理论，Forrestal等 ［ 23 23 - 25 25 23-25 ］ 开展了针对于不同靶体材料（如岩石、混凝土及土壤等）的侵彻研究，大大拓宽了这一理论的应用范围。Shi等人 ［ 26 26 ］ 针对干砂的冲击贯入，建立了考虑压缩性的动态空腔膨胀理论。动态空腔膨胀模型准确地刻画了稳定贯入过程中贯入阻力和贯入速度平方间的线性关系，展现了清晰的物理图像，而土壤强度参数则反映在阻力系数当中，这为土壤强度参数反演提供了理论基础。 本文基于冲击贯入触探技术和动态空腔膨胀模型，开展了土壤Mohr-Coulomb参数（内聚力和内摩擦角）反演方法的研究。利用土壤冲击贯入实验中测得的阻力-速度曲线、加速度时程曲线、速度时程曲线以及贯入深度-初速度关系，对贯入阻力系数进行了拟合，进而对Mohr-Coulomb参数进行了反演。最后，对参数可解域和参数敏感性进行了讨论。 2 动态球形空腔膨胀模型 1 一般来讲，在对土壤的冲击贯入过程中，贯入仪可视作刚体。此时，动态球形空腔膨胀模型为计算弹-靶相互作用提供了一个简单的物理图像。这个模型将贯入仪看作是由一系列与其内切的球形空腔组成，空腔球心在贯入仪的对称轴上 ［ 28 28 ］ 。贯入仪与土体的相互作用可看作是这些球形空腔在无限大介质中均匀膨胀时与土体介质间的相互挤压。这样一来，贯入仪侧表面上每一点的法向应力可近似为球形空腔膨胀时内表面的径向应力，它与空腔膨胀的速度有关，而空腔膨胀速度为贯入仪侧表面上该点的法向速度分量 ［ 25 25 ， 28 28 ］ 。由于球对称性的引入，经过以上简化后，贯入仪侧表面上的压力分布就可转化为一个一维问题进行求解。贯入仪图示及球形空腔膨胀模型如 图1 图1 。 球形空腔的膨胀将导致土体介质出现相应的动态响应，形成弹塑性分区 ［ 25 25 ］ （如 图1 图1 （b））。记空腔膨胀速度为 V 、弹塑性边界扩张速度为 c 、初始时刻空腔的半径为零，则在 t 时刻塑性区的范围为（ Vt ， ct ），弹性区的范围为（ ct ，+ ∞ ）。在弹性区内，土壤的本构关系采用Hooke定律来描述，并且假定弹性区内土壤不可压缩。在塑性区内，采用体积锁定的假设 ［ 25 25 ］ 来描述土壤的应力应变关系，即发生塑性变形后土壤的体积应变为常数 η * ，不随平均应力的增大而变化。在这个假定下，静水压力将不影响土壤的后继屈服。记 ρ 0 和 ρ * 分别是土壤的初始密度和锁定密度，则： <math display="block" id="M1"><msup><mrow><mi>η</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup><mo>=</mo><mn mathvariant="normal">1</mn><mo>-</mo><msub><mrow><mi>ρ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>/</mo><msup><mrow><mi>ρ</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465462&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465457&type= 22.26733208 5.07999992 . 土壤的屈服准则采用Mohr-Coulomb准则，在球坐标系下可写为 ［ 25 25 ］ ： <math display="block" id="M2"><msub><mrow><mi>σ</mi></mrow><mrow><mi>r</mi></mrow></msub><mo>-</mo><msub><mrow><mi>σ</mi></mrow><mrow><mi>θ</mi></mrow></msub><mo>=</mo><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>+</mo><mi>λ</mi><mi>p</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465473&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465469&type= 26.24666595 4.57200003 ， 其中： σ r ， σ θ 分别为径向应力和环向应力（以压为正）； p =（ σ r +2 σ θ ）/3为平均应力； τ 0 ， λ 是与土壤内聚力 C 和内摩擦角 φ 有关的常数，且有 ［ 29 29 ］ ： <math display="block" id="M3"><mtable columnalign="left"><mtr><mtd><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>=</mo><mfrac><mrow><mn mathvariant="normal">6</mn><mi>C</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mi>φ</mi></mrow><mrow><mn mathvariant="normal">3</mn><mo>-</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mi>φ</mi></mrow></mfrac></mtd></mtr><mtr><mtd><mi>λ</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">6</mn><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mi>φ</mi></mrow><mrow><mn mathvariant="normal">3</mn><mo>-</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mi>φ</mi></mrow></mfrac></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465484&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465478&type= 22.26733208 19.47333336 . 冲击贯入过程中，质量、动量及能量守恒定律均需得到满足。对于贯入阻力和贯入深度的解答，仅利用质量守恒和动量守恒关系即可得到。在球坐标系下，拉格朗日坐标表示的质量守恒方程和动量守恒方程可写为： <math display="block" id="M4"><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">3</mn></mrow></mfrac><mfrac><mrow><mo>∂</mo></mrow><mrow><mo>∂</mo><mi>r</mi></mrow></mfrac><msup><mrow><mo stretchy="false">(</mo><mi>r</mi><mo>+</mo><mi>u</mi><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup><mo>=</mo><mfrac><mrow><msub><mrow><mi>ρ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><msup><mrow><mi>ρ</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup></mrow></mfrac><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465496&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465492&type= 33.10466385 8.97466660 ， <math display="block" id="M5"><mtable><mtr><mtd><msup><mrow><mo stretchy="false">(</mo><mi>r</mi><mo>+</mo><mi>u</mi><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mfrac><mrow><mo>∂</mo><msub><mrow><mi>σ</mi></mrow><mrow><mi>r</mi></mrow></msub></mrow><mrow><mo>∂</mo><mi>r</mi></mrow></mfrac><mo>+</mo><mn mathvariant="normal">2</mn><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn><mo>+</mo><mfrac><mrow><mo>∂</mo><mi>u</mi></mrow><mrow><mo>∂</mo><mi>r</mi></mrow></mfrac></mrow></mfenced><mo stretchy="false">(</mo><mi>r</mi><mo>+</mo><mi>u</mi><mo stretchy="false">)</mo><mo stretchy="false">(</mo><msub><mrow><mi>σ</mi></mrow><mrow><mi>r</mi></mrow></msub><mo>-</mo><msub><mrow><mi>σ</mi></mrow><mrow><mi>θ</mi></mrow></msub><mo stretchy="false">)</mo><mo>+</mo></mtd></mtr><mtr><mtd><msub><mrow><mi>ρ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mfrac><mrow><msup><mrow><mo>∂</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mi>u</mi></mrow><mrow><mo>∂</mo><msup><mrow><mi>t</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac><mo>=</mo><mn mathvariant="normal">0</mn><mo>,</mo></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465505&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465500&type= 69.17266846 19.38866615 其中： r 是拉格朗日坐标， u 为径向位移（以向外为正）。空腔表面处的边界条件为： <math display="block" id="M6"><mi>u</mi><mo stretchy="false">(</mo><mi>r</mi><mo>=</mo><mn mathvariant="normal">0</mn><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mo>=</mo><mi>V</mi><mi>t</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465517&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465512&type= 25.56933212 3.55599999 . 弹塑性界面处的Hugoniot条件可表示为 ［ 24 24 ］ ： <math display="block" id="M7"><msup><mrow><mi>ρ</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup><mfenced separators="|"><mrow><msub><mrow><mi>v</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>-</mo><mi>c</mi></mrow></mfenced><mo>=</mo><msub><mrow><mi>ρ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mfenced separators="|"><mrow><msub><mrow><mi>v</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>-</mo><mi>c</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465527&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465522&type= 36.49133301 3.97933316 ， <math display="block" id="M8"><msub><mrow><mi>σ</mi></mrow><mrow><mi>r</mi><mn mathvariant="normal">2</mn></mrow></msub><mo>+</mo><msup><mrow><mi>ρ</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup><msub><mrow><mi>v</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mfenced separators="|"><mrow><msub><mrow><mi>v</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>-</mo><mi>c</mi></mrow></mfenced><mo>=</mo><msub><mrow><mi>σ</mi></mrow><mrow><mi>r</mi><mn mathvariant="normal">1</mn></mrow></msub><mo>+</mo><msub><mrow><mi>ρ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><msub><mrow><mi>v</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mfenced separators="|"><mrow><msub><mrow><mi>v</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>-</mo><mi>c</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465537&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465532&type= 59.85933304 3.97933316 ， 其中， v 是介质的径向速度，下标1、2分别代表弹塑性界面处的弹性区和塑性区。 利用以上本构关系和边界条件求解控制方程（4）和（5）可解析地得到贯入仪端部表面的法向应力分布为 ［ 25 25 ］ ： <math display="block" id="M9"><msub><mrow><mi>σ</mi></mrow><mrow><mi mathvariant="normal">n</mi></mrow></msub><mo>=</mo><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mi>A</mi><mo>+</mo><msub><mrow><mi>ρ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mi>B</mi><msup><mrow><mfenced separators="|"><mrow><msub><mrow><mi>V</mi></mrow><mrow><mi>z</mi></mrow></msub><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mi>ϕ</mi></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465548&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465543&type= 42.07933044 5.50333309 ， 其中： V z 是贯入仪在竖直方向的速度； A ， B 是与 τ 0 / E ， λ 和 η * 有关的参数。根据 λ 的不同取值， A ， B 的表达式存在如下三种不同的情况 ［ 25 25 ］ ： （1） λ ≠0且 λ ≠3/4： <math display="block" id="M10"><mi>A</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mi>α</mi><msup><mrow><mi>γ</mi></mrow><mrow><mn mathvariant="normal">2</mn><mi>α</mi></mrow></msup></mrow></mfrac><msup><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn><mo>+</mo><mfrac><mrow><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><mn mathvariant="normal">2</mn><mi>E</mi></mrow></mfrac></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn><mi>α</mi></mrow></msup><mo>-</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mi>λ</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465559&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465554&type= 40.21666718 10.07533360 ， <math display="block" id="M11"><mtable><mtr><mtd><mi>B</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">3</mn></mrow><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn><mo>-</mo><msup><mrow><mi>η</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup></mrow></mfenced><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn><mo>-</mo><mn mathvariant="normal">2</mn><mi>α</mi></mrow></mfenced><mfenced separators="|"><mrow><mn mathvariant="normal">2</mn><mo>-</mo><mi>α</mi></mrow></mfenced></mrow></mfrac><mo>+</mo></mtd></mtr><mtr><mtd><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><msup><mrow><mi>γ</mi></mrow><mrow><mn mathvariant="normal">2</mn><mi>α</mi><mo>+</mo><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac><msup><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn><mo>+</mo><mfrac><mrow><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><mn mathvariant="normal">2</mn><mi>E</mi></mrow></mfrac></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn><mi>α</mi></mrow></msup><mfenced open="{" close="" separators="|"><mrow><mfrac><mrow><mn mathvariant="normal">3</mn><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><mi>E</mi></mrow></mfrac><mo>+</mo><msup><mrow><mi>η</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup><msup><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn><mo>-</mo><mfrac><mrow><mn mathvariant="normal">3</mn><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><mn mathvariant="normal">2</mn><mi>E</mi></mrow></mfrac></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfenced><mo>-</mo></mtd></mtr><mtr><mtd><mfenced open="" close="}" separators="|"><mrow><mo> </mo><mo> </mo><mo> </mo><mo> </mo><mfrac><mrow><msup><mrow><mi>γ</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup><mfenced open="[" close="]" separators="|"><mrow><mn mathvariant="normal">2</mn><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn><mo>-</mo><msup><mrow><mi>η</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup></mrow></mfenced><mfenced separators="|"><mrow><mn mathvariant="normal">2</mn><mo>-</mo><mi>α</mi></mrow></mfenced><mo>+</mo><mn mathvariant="normal">3</mn><msup><mrow><mi>γ</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup></mrow></mfenced></mrow><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn><mo>-</mo><msup><mrow><mi>η</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup></mrow></mfenced><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn><mo>-</mo><mn mathvariant="normal">2</mn><mi>α</mi></mrow></mfenced><mfenced separators="|"><mrow><mn mathvariant="normal">2</mn><mo>-</mo><mi>α</mi></mrow></mfenced><msup><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn><mo>+</mo><mrow><mrow><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mo>/</mo><mrow><mn mathvariant="normal">2</mn><mi>E</mi></mrow></mrow></mrow></mfenced></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msup></mrow></mfrac></mrow></mfenced><mo>,</mo></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465568&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465565&type= 65.95532990 38.52333450 其中： <math display="block" id="M12"><mtable><mtr><mtd><mi>α</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">3</mn><mi>λ</mi></mrow><mrow><mn mathvariant="normal">3</mn><mo>+</mo><mn mathvariant="normal">2</mn><mi>λ</mi></mrow></mfrac></mtd></mtr><mtr><mtd><mo> </mo><mo> </mo><mi>γ</mi><mo>=</mo><mfrac><mrow><mi>V</mi></mrow><mrow><mi>c</mi></mrow></mfrac><mo>=</mo><msup><mrow><mfenced open="[" close="]" separators="|"><mrow><msup><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn><mo>+</mo><mfrac><mrow><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><mn mathvariant="normal">2</mn><mi>E</mi></mrow></mfrac></mrow></mfenced></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup><mo>-</mo><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn><mo>-</mo><msup><mrow><mi>η</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup></mrow></mfenced></mrow></mfenced></mrow><mrow><mrow><mrow><mn mathvariant="normal">1</mn></mrow><mo>/</mo><mrow><mn mathvariant="normal">3</mn></mrow></mrow></mrow></msup></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465575&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465572&type= 52.66266251 21.92866707 . （2） λ= 3/4： <math display="block" id="M13"><mi>A</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">2</mn></mrow><mrow><mi>γ</mi></mrow></mfrac><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn><mo>+</mo><mfrac><mrow><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><mn mathvariant="normal">2</mn><mi>E</mi></mrow></mfrac></mrow></mfenced><mo>-</mo><mfrac><mrow><mn mathvariant="normal">4</mn></mrow><mrow><mn mathvariant="normal">3</mn></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465600&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465595&type= 33.61266708 9.22866726 ， <math display="block" id="M14"><mtable><mtr><mtd><mi>B</mi><mo>=</mo><mfrac><mrow><mo>-</mo><mn mathvariant="normal">2</mn><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mi>γ</mi></mrow><mrow><mn mathvariant="normal">1</mn><mo>-</mo><msup><mrow><mi>η</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup></mrow></mfrac><mo>+</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><msup><mrow><mi>γ</mi></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup></mrow></mfrac><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn><mo>+</mo><mfrac><mrow><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><mn mathvariant="normal">2</mn><mi>E</mi></mrow></mfrac></mrow></mfenced><mfenced open="[" close="]" separators="|"><mrow><mfrac><mrow><mn mathvariant="normal">3</mn><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><mi>E</mi></mrow></mfrac><mo>+</mo><msup><mrow><mi>η</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup><msup><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn><mo>-</mo><mfrac><mrow><mn mathvariant="normal">3</mn><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><mn mathvariant="normal">2</mn><mi>E</mi></mrow></mfrac></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfenced><mo>-</mo></mtd></mtr><mtr><mtd><mtext>          </mtext><mfrac><mrow><mn mathvariant="normal">2</mn></mrow><mrow><mn mathvariant="normal">3</mn></mrow></mfrac><mfenced open="[" close="]" separators="|"><mrow><msup><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn><mo>+</mo><mfrac><mrow><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><mn mathvariant="normal">2</mn><mi>E</mi></mrow></mfrac></mrow></mfenced></mrow><mrow><mo>-</mo><mn mathvariant="normal">3</mn></mrow></msup><mo>-</mo><mfrac><mrow><mn mathvariant="normal">3</mn><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn><mo>+</mo><mrow><mrow><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mo>/</mo><mrow><mn mathvariant="normal">2</mn><mi>E</mi></mrow></mrow></mrow></mfenced></mrow><mrow><mn mathvariant="normal">1</mn><mo>-</mo><msup><mrow><mi>η</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup></mrow></mfrac></mrow></mfenced><mo>.</mo><mtext>          </mtext><mo stretchy="false">(</mo><mn mathvariant="normal">14</mn><mo stretchy="false">)</mo></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465609&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465605&type= 76.36933136 24.72266769 （3） λ= 0： <math display="block" id="M15"><mi>A</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">2</mn></mrow><mrow><mn mathvariant="normal">3</mn></mrow></mfrac><mo>-</mo><mn mathvariant="normal">2</mn><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mfenced open="[" close="]" separators="|"><mrow><mi>γ</mi><msup><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn><mo>+</mo><mrow><mrow><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mo>/</mo><mrow><mn mathvariant="normal">2</mn><mi>E</mi></mrow></mrow></mrow></mfenced></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465619&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465614&type= 47.83666611 7.87400007 ， <math display="block" id="M16"><mtable><mtr><mtd><mi>B</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">3</mn></mrow><mrow><mn mathvariant="normal">2</mn><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn><mo>-</mo><msup><mrow><mi>η</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup></mrow></mfenced></mrow></mfrac><mo>+</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><msup><mrow><mi>γ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac><mfenced open="[" close="]" separators="|"><mrow><mfrac><mrow><mn mathvariant="normal">3</mn><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><mi>E</mi></mrow></mfrac><mo>+</mo><msup><mrow><mi>η</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup><msup><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn><mo>-</mo><mfrac><mrow><mn mathvariant="normal">3</mn><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><mn mathvariant="normal">2</mn><mi>E</mi></mrow></mfrac></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfenced><mo>-</mo></mtd></mtr><mtr><mtd><mfrac><mrow><mi>γ</mi></mrow><mrow><mn mathvariant="normal">2</mn><msup><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn><mo>+</mo><mrow><mrow><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mo>/</mo><mrow><mn mathvariant="normal">2</mn><mi>E</mi></mrow></mrow></mrow></mfenced></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msup></mrow></mfrac><mfenced open="[" close="]" separators="|"><mrow><mn mathvariant="normal">1</mn><mo>+</mo><mfrac><mrow><mn mathvariant="normal">3</mn><msup><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn><mo>+</mo><mrow><mrow><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mo>/</mo><mrow><mn mathvariant="normal">2</mn><mi>E</mi></mrow></mrow></mrow></mfenced></mrow><mrow><mn mathvariant="normal">3</mn></mrow></msup></mrow><mrow><mn mathvariant="normal">1</mn><mo>-</mo><msup><mrow><mi>η</mi></mrow><mrow><mi mathvariant="normal">*</mi></mrow></msup></mrow></mfrac></mrow></mfenced><mo>.</mo></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465630&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465625&type= 69.17266846 26.16200066 式（12） 式（12） 中的第2式表明，在贯入过程中空腔膨胀速度 V 与弹塑性边界扩张速度 c 的比值始终保持为常数。随着贯入深度的增加，贯入速度减小，空腔膨胀速度 V 也逐渐减小，对应的弹塑性区域的扩张逐渐变慢。 对于贯入仪相对土壤垂直贯入的情况，将 式（9） 式（9） 在贯入仪端部积分，则得到贯入仪所受到的竖直方向阻力 F z 为 ［ 25 25 ］ <math display="block" id="M17"><msub><mrow><mi>F</mi></mrow><mrow><mi>z</mi></mrow></msub><mo>=</mo><msub><mrow><mi>α</mi></mrow><mrow><mi>s</mi></mrow></msub><mo>+</mo><msub><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow></msub><msubsup><mrow><mi>V</mi></mrow><mrow><mi>z</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465637&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465633&type= 23.36800003 4.82600021 ， 其中： <math display="block" id="M18"><msub><mrow><mi>α</mi></mrow><mrow><mi>s</mi></mrow></msub><mo>=</mo><mi>π</mi><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mi>A</mi><mfenced open="[" close="]" separators="|"><mrow><mn mathvariant="normal">1</mn><mo>+</mo><mn mathvariant="normal">4</mn><mi>μ</mi><msup><mrow><mi>ψ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mfenced separators="|"><mrow><mrow><mrow><mi>π</mi></mrow><mo>/</mo><mrow><mn mathvariant="normal">2</mn></mrow></mrow><mo>-</mo><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfenced><mo>-</mo><mi>μ</mi><mfenced separators="|"><mrow><mn mathvariant="normal">2</mn><mi>ψ</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></mfenced><msup><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">4</mn><mi>ψ</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></mfenced></mrow><mrow><mrow><mrow><mn mathvariant="normal">1</mn></mrow><mo>/</mo><mrow><mn mathvariant="normal">2</mn></mrow></mrow></mrow></msup></mrow></mfenced><mo>,</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465649&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465645&type= 92.37133789 6.68866634 <math display="block" id="M19"><mtable columnalign="left"><mtr><mtd><msub><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow></msub><mo>=</mo><mi>π</mi><msup><mrow><mi>r</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><msub><mrow><mi>ρ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mi>B</mi><mfenced open="[" close="" separators="|"><mrow><mfrac><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">8</mn><mi>ψ</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></mfenced></mrow><mrow><mn mathvariant="normal">24</mn><msup><mrow><mi>ψ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac><mo>+</mo><mi>μ</mi><msup><mrow><mi>ψ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mfenced separators="|"><mrow><mrow><mrow><mi>π</mi></mrow><mo>/</mo><mrow><mn mathvariant="normal">2</mn></mrow></mrow><mo>-</mo><msub><mrow><mi>θ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfenced></mrow></mfenced><mo>-</mo></mtd></mtr><mtr><mtd><mfenced open="" close="]" separators="|"><mrow><mi>μ</mi><mfrac><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">2</mn><mi>ψ</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></mfenced><mfenced separators="|"><mrow><mn mathvariant="normal">6</mn><msup><mrow><mi>ψ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><mn mathvariant="normal">4</mn><mi>ψ</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></mfenced><msup><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">4</mn><mi>ψ</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></mfenced></mrow><mrow><mrow><mrow><mn mathvariant="normal">1</mn></mrow><mo>/</mo><mrow><mn mathvariant="normal">2</mn></mrow></mrow></mrow></msup></mrow><mrow><mn mathvariant="normal">24</mn><msup><mrow><mi>ψ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac></mrow></mfenced><mo>.</mo></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465659&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465655&type= 65.44733429 25.06133461 <math id="M20"><msub><mrow><mi>α</mi></mrow><mrow><mi>s</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465669&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465665&type= 2.79399991 3.72533321 和 <math id="M21"><msub><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465678&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465675&type= 3.04800010 4.57200003 是与土壤力学特性有关的系数。其中 ψ = s /（2 r ）是贯入仪头部的曲径比，并且有 θ 0 =arcsin［（2 ψ -1）/（2 ψ ）］（见 图1 图1 （a））； μ 为贯入仪与土壤间的滑动摩擦系数，二者之间的切应力可表示为 σ τ = μσ n 。如果忽略贯入仪与土体间的摩擦（ μ =0），并令： <math display="block" id="M22"><mi>Φ</mi><mo>=</mo><mfrac><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">8</mn><mi>ψ</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></mfenced></mrow><mrow><mn mathvariant="normal">24</mn><msup><mrow><mi>ψ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465686&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465683&type= 21.92866707 9.73666668 ， 式（18） 式（18） 和 式（19） 式（19） 可进一步简化为： <math display="block" id="M23"><msub><mrow><mi>α</mi></mrow><mrow><mi>s</mi></mrow></msub><mo>=</mo><mi mathvariant="normal">π</mi><msup><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mi>A</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465696&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465691&type= 18.79599953 3.97933316 ， <math display="block" id="M24"><msub><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow></msub><mo>=</mo><mi mathvariant="normal">π</mi><msup><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mi>Φ</mi><msub><mrow><mi>ρ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mi>B</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465706&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465701&type= 21.84399986 4.82600021 . 由 式（9） 式（9） 和 式（17） 式（17） 可以看到，在当前模型假设下，贯入仪端部的应力分布及贯入阻力由两部分贡献，一项是静阻力项，另一项是流动阻力项（即速度平方项），而黏性效应和附加质量效应 ［ 14 14 ］ 未被考虑。 贯入仪在竖直方向的贯入速度 V z 可由牛顿第二定律计算得到，注意到 式（17） 式（17） ，有： <math display="block" id="M25"><mi>m</mi><mfrac><mrow><mi mathvariant="normal">d</mi><msub><mrow><mi>V</mi></mrow><mrow><mi>z</mi></mrow></msub></mrow><mrow><mi mathvariant="normal">d</mi><mi>t</mi></mrow></mfrac><mo>=</mo><mo>-</mo><mo stretchy="false">(</mo><msub><mrow><mi>α</mi></mrow><mrow><mi>s</mi></mrow></msub><mo>+</mo><msub><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow></msub><msubsup><mrow><mi>V</mi></mrow><mrow><mi>z</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465716&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465712&type= 36.15266418 7.95866632 ， 其中， m 为贯入仪的质量。积分上式即得到： <math display="block" id="M26"><msub><mrow><mi>V</mi></mrow><mrow><mi>z</mi></mrow></msub><mo>=</mo><msup><mrow><mfenced separators="|"><mrow><mfrac><mrow><msub><mrow><mi>α</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow><mrow><msub><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow></mfrac></mrow></mfenced></mrow><mrow><mrow><mrow><mn mathvariant="normal">1</mn></mrow><mo>/</mo><mrow><mn mathvariant="normal">2</mn></mrow></mrow></mrow></msup><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mfenced open="{" close="}" separators="|"><mrow><mi mathvariant="normal">a</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mfenced open="[" close="]" separators="|"><mrow><msup><mrow><mfenced separators="|"><mrow><mfrac><mrow><msub><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow><mrow><msub><mrow><mi>α</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow></mfrac></mrow></mfenced></mrow><mrow><mrow><mrow><mn mathvariant="normal">1</mn></mrow><mo>/</mo><mrow><mn mathvariant="normal">2</mn></mrow></mrow></mrow></msup><msub><mrow><mi>V</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfenced><mo>-</mo><mfrac><mrow><msup><mrow><mfenced separators="|"><mrow><msub><mrow><mi>α</mi></mrow><mrow><mi>s</mi></mrow></msub><msub><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow></mfenced></mrow><mrow><mrow><mrow><mn mathvariant="normal">1</mn></mrow><mo>/</mo><mrow><mn mathvariant="normal">2</mn></mrow></mrow></mrow></msup><mi>t</mi></mrow><mrow><mi>m</mi></mrow></mfrac></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465727&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465721&type= 74.93000031 12.86933327 ， 其中， V 0 = V z （ t =0）是贯入仪的初始贯入速度。贯入过程中贯入仪加速度大小为： <math display="block" id="M27"><msub><mrow><mi>a</mi></mrow><mrow><mi>z</mi></mrow></msub><mo>=</mo><mfrac><mrow><mrow><mrow><msub><mrow><mi>α</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow><mo>/</mo><mrow><mi>m</mi></mrow></mrow></mrow><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><msup><mrow><mi mathvariant="normal">s</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mfenced open="{" close="}" separators="|"><mrow><mi mathvariant="normal">a</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mfenced open="[" close="]" separators="|"><mrow><msup><mrow><mfenced separators="|"><mrow><mfrac><mrow><msub><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow><mrow><msub><mrow><mi>α</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow></mfrac></mrow></mfenced></mrow><mrow><mrow><mrow><mn mathvariant="normal">1</mn></mrow><mo>/</mo><mrow><mn mathvariant="normal">2</mn></mrow></mrow></mrow></msup><msub><mrow><mi>V</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfenced><mo>-</mo><mfrac><mrow><msup><mrow><mfenced separators="|"><mrow><msub><mrow><mi>α</mi></mrow><mrow><mi>s</mi></mrow></msub><msub><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow></mfenced></mrow><mrow><mrow><mrow><mn mathvariant="normal">1</mn></mrow><mo>/</mo><mrow><mn mathvariant="normal">2</mn></mrow></mrow></mrow></msup><mi>t</mi></mrow><mrow><mi>m</mi></mrow></mfrac></mrow></mfenced></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465735&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465730&type= 64.76999664 18.03400040 . 最终贯入深度 P 为： <math display="block" id="M28"><mi>P</mi><mo>=</mo><mfrac><mrow><mi>m</mi></mrow><mrow><mn mathvariant="normal">2</mn><msub><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow></mfrac><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mfenced separators="|"><mrow><mn mathvariant="normal">1</mn><mo>+</mo><mfrac><mrow><msub><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow></msub><msubsup><mrow><mi>V</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow><mrow><msub><mrow><mi>α</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow></mfrac></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465744&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465740&type= 34.54399872 9.65200043 . 3 土壤Mohr-Coulomb参数的反演 1 3.1　Mohr-Coulomb参数 C ， φ 的反演方法 2 土壤Mohr-Coulomb参数即土壤的内聚力 C 和内摩擦角 φ 。由于动态空腔膨胀模型刻画了土壤力学特性参数对弹靶相互作用的影响，因此对 C 和 φ 的反演可以建立在动态空腔膨胀模型的基础上。由 式（3） 式（3） 可知： <math display="block" id="M29"><mtable><mtr><mtd><mi>C</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">3</mn><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><mn mathvariant="normal">2</mn><mroot><mrow><mfenced separators="|"><mrow><mn mathvariant="normal">3</mn><mo>+</mo><mn mathvariant="normal">2</mn><mi>λ</mi></mrow></mfenced><mfenced separators="|"><mrow><mn mathvariant="normal">3</mn><mo>-</mo><mi>λ</mi></mrow></mfenced></mrow><mrow/></mroot></mrow></mfrac></mtd></mtr><mtr><mtd><mi>φ</mi><mo>=</mo><mi mathvariant="normal">a</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mfrac><mrow><mn mathvariant="normal">3</mn><mi>λ</mi></mrow><mrow><mn mathvariant="normal">6</mn><mo>+</mo><mi>λ</mi></mrow></mfrac></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465753&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465750&type= 38.94666672 20.40466690 . 上式表明，要反演内聚力 C 和内摩擦角 φ ，实际上就是要求得参数 τ 0 和 λ 的值。由 式（18） 式（18） 、 式（19） 式（19） （或对于不考虑摩擦时的方程（21）和（22））以及 A 、 B 的表达式可见，参数 τ 0 、 λ 包含在阻力系数 α s 、 β s 中。一旦求得 α s 和 β s ，联立方程（18）和（19）（或对于不考虑摩擦的情况，联立方程（21）和（22）），即可求得参数 τ 0 和 λ 的值。 在联立方程（18）和（19）（或方程（21）和（22））求解参数 τ 0 和 λ 时，需针对 A 、 B 取值的三种情况进行讨论。由于参数 A 、 B 的表达式过于复杂，一般不能得到 τ 0 和 λ 的解析表达式。此时，可采用如下几何化方法进行求解。为简便起见，忽略掉贯入仪与土体间的摩擦（即令 μ =0），不影响对方法的讨论。 （1）如果 λ ≠0且 λ ≠3/4： 此时，由 式（10） 式（10） 和 式（11） 式（11） 可见，参数 A 、 B 均为（ τ 0 ， λ ）的函数。根据 式（21） 式（21） 和 式（22） 式（22） ， α s 、 β s 也均是（ τ 0 ， λ ）的函数。因此，由 式（21） 式（21） 和 式（22） 式（22） 可定义如下曲面： <math display="block" id="M30"><msub><mrow><mi>Π</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mfenced separators="|"><mrow><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>,</mo><mi>λ</mi></mrow></mfenced><mo>=</mo><mfrac><mrow><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><mi>E</mi></mrow></mfrac><mi>A</mi><mo>-</mo><mfrac><mrow><msub><mrow><mi>α</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow><mrow><mi mathvariant="normal">π</mi><msup><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mi>E</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465763&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465756&type= 39.03133011 8.21266651 ， <math display="block" id="M31"><msub><mrow><mi>Π</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mfenced separators="|"><mrow><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>,</mo><mi>λ</mi></mrow></mfenced><mo>=</mo><mi>B</mi><mo>-</mo><mfrac><mrow><msub><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow><mrow><mi mathvariant="normal">π</mi><msup><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mi>Φ</mi><msub><mrow><mi>ρ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465789&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465768&type= 37.42266846 8.97466660 . 此种情况下若存在 τ 0 和 λ 的解，那么参数 τ 0 和 λ 必是曲面 Π 1 和 Π 2 的公共零点。 （2）如果 λ =3/4或 λ =0： 此时，由式（13）~ 式（16） 式（16） 可见，参数 A 、 B 仅是 τ 0 的函数。因此 α s 、 β s 也仅是 τ 0 的函数。于是，根据 式（21） 式（21） 和 式（22） 式（22） 可定义如下曲线： <math display="block" id="M32"><msub><mrow><mi>L</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mfenced separators="|"><mrow><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfenced><mo>=</mo><mfrac><mrow><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow><mrow><mi>E</mi></mrow></mfrac><mi>A</mi><mo>-</mo><mfrac><mrow><msub><mrow><mi>α</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow><mrow><mi mathvariant="normal">π</mi><msup><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mi>E</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465799&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465794&type= 35.22133255 8.21266651 ， <math display="block" id="M33"><msub><mrow><mi>L</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mfenced separators="|"><mrow><msub><mrow><mi>τ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfenced><mo>=</mo><mi>B</mi><mo>-</mo><mfrac><mrow><msub><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow><mrow><mi mathvariant="normal">π</mi><msup><mrow><mi>a</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mi>Φ</mi><msub><mrow><mi>ρ</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465811&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465804&type= 33.69733429 8.97466660 . 此种情况下若存在 τ 0 的解，则 τ 0 必是曲线 L 1 和 L 2 的公共零点。 至此，对 C 、 φ 的反演就归结为对阻力系数 α s 和 β s 的计算，它们可通过贯入阻力、贯入速度、贯入加速度以及最终贯入深度的测量结果计算得到。下面就针对这些情况进行具体说明。 3.2　基于阻力-速度关系进行的反演 2 在冲击贯入过程中，最重要的物理量是贯入仪受到的贯入阻力 F z ，它综合反映了靶体介质的物理力学性质以及贯入仪与靶体材料之间的相互作用。动态空腔膨胀模型的核心就是揭示了贯入过程中的稳定阶段贯入阻力 F z （ t ）与速度平方 <math id="M34"><msubsup><mrow><mi>V</mi></mrow><mrow><mi>z</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mfenced separators="|"><mrow><mi>t</mi></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465820&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465816&type= 9.22866726 3.97933316 间的线性关系（ 式（17） 式（17） ）。因此，我们首先利用这个关系来反演土壤的Mohr-Coulomb参数。 实际当中，贯入阻力 F z （ t ）可由加速度传感器测量得到；而贯入速度 V z （ t ）则可由加速度时程曲线积分得出。因此，利用测量数据集｛ F z ， i ， V z ， i ｝并基于 式（17） 式（17） 进行线性拟合，即可得到参数 α s 和 β s 的数值，并进而求得内聚力 C 和内摩擦角 φ 。为此，构造误差函数： <math display="block" id="M35"><msub><mrow><mi>Q</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>=</mo><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><msup><mrow><mfenced separators="|"><mrow><msub><mrow><mi>F</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><mo>-</mo><msub><mrow><mi>α</mi></mrow><mrow><mi>s</mi></mrow></msub><mo>-</mo><msub><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow></msub><msubsup><mrow><mi>V</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465832&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465827&type= 44.19599915 8.55133343 . 令 <math display="block" id="M36"><mfrac><mrow><mo>∂</mo><msub><mrow><mi>Q</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow><mrow><mo>∂</mo><msub><mrow><mi>α</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow></mfrac><mo>=</mo><mn mathvariant="normal">0</mn><mo>,</mo><mo> </mo><mfrac><mrow><mo>∂</mo><msub><mrow><mi>Q</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></mrow><mrow><mo>∂</mo><msub><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow></mfrac><mo>=</mo><mn mathvariant="normal">0</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465843&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465837&type= 28.87133408 8.80533314 . 可解得： <math display="block" id="M37"><msub><mrow><mi>α</mi></mrow><mrow><mi>s</mi></mrow></msub><mo>=</mo><mfrac><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><msub><mrow><mi>F</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub></mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><msubsup><mrow><mi>V</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msubsup></mrow><mo>-</mo><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><msub><mrow><mi>F</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><msubsup><mrow><mi>V</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><msubsup><mrow><mi>V</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow></mrow><mrow><mi>N</mi><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><msubsup><mrow><mi>V</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msubsup></mrow><mo>-</mo><msup><mrow><mfenced separators="|"><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><msubsup><mrow><mi>V</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465853&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465848&type= 60.19800186 20.31999969 ， <math display="block" id="M38"><msub><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow></msub><mo>=</mo><mfrac><mrow><mi>N</mi><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><msub><mrow><mi>F</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub><msubsup><mrow><mi>V</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow><mo>-</mo><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><msub><mrow><mi>F</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow></msub></mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><msubsup><mrow><mi>V</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow></mrow><mrow><mi>N</mi><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><msubsup><mrow><mi>V</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow><mrow><mn mathvariant="normal">4</mn></mrow></msubsup></mrow><mo>-</mo><msup><mrow><mfenced separators="|"><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><msubsup><mrow><mi>V</mi></mrow><mrow><mi>z</mi><mo>,</mo><mi>i</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465861&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465857&type= 53.34000015 20.31999969 ， 其中， N 是数据采样点数目。得到 α s ， β s 的数值后，即可按照上节所述步骤进行 C ， φ 值的反演。下面采用Forrestal等人 ［ 25 25 ］ 报道的实验结果来对这个反演方法进行详细说明。 文献［ 25 25 ］中，通过对不同采样地点取样的土壤进行三轴压缩及单轴压缩标定实验，测得的土壤参数如 表1 表1 所示。需要说明的是采样地点和现场贯入实验地点的土层按照深度可分为三层 ［ 25 25 ］ ，第一层（Ⅰ：0~1.8 m）的弹性模量为120 MPa，锁定体积应变为0.17；第二层（Ⅱ：1.8~3.0 m）的弹性模量为150 MPa，锁定体积应变为0.10；第三层（Ⅲ：3.0~6.1 m）的弹性模量为210 MPa，锁定体积应变为0.13。为了简化模型计算，文献［ 25 25 ］忽略了土壤的分层，将其视作均匀介质，并取三层土壤的平均弹性模量（160 MPa）和平均锁定体积应变（0.13）作为土壤基本参数。本文按照同样的假定进行计算。 密度 ρ 0 /（g·cm -3 ） 平均弹性模量 E /MPa τ 0 /MPa λ 平均锁定体积应变 η * 1.86 160 10 0 0.13 以文献［ 25 25 ］中的贯入实验3为例，贯入仪质量为23.1 kg，弹体半径 r =47.6 mm，弹头曲径比 ψ =3。忽略贯入过程中贯入仪与土体间的摩擦。由实验测得的加速度数据计算得到贯入阻力和贯入速度后，给出的阻力-速度平方关系如 图2 图2 所示。由 图2 图2 可见，除初始阶段贯入仪与土体碰撞造成的震荡以及贯入结束阶段的土体扰动外，在中间的稳定贯入阶段，贯入阻力与贯入速度的平方之间呈线性关系，这表明动态空腔膨胀模型揭示了贯入过程的关键特征。根据 式（17） 式（17） ，对稳定贯入阶段的阻力-速度关系进行线性拟合，结果如 图2 图2 中红线所示，拟合得到的阻力系数为 α s =85.39 kN， β s =2.52 kg/m。图中 A 、 B 两点分别是拟合范围的起点和终点，黑色圆点是参与拟合的数据。关于拟合范围的选取，将在第4.2节中进行详细讨论。 由拟合得到的阻力系数 α s 、 β s ，可按照 式（28） 式（28） 和 式（29） 式（29） 画出曲面 Π 1 和 Π 2 的图像，如 图3 图3 。由 图3 图3 （a）、 图3 图3 （b）可见，曲面 Π 1 与坐标平面 Π 1 =0有交线 l 1 （曲面 Π 1 零点的集合）；曲面 Π 2 与坐标平面 Π 2 =0有交线 l 2 （曲面 Π 2 零点的集合）。将交线 l 1 和 l 2 画在同一坐标系下，如 图4 图4 所示，可见二者有公共交点（ τ 0 =5.52 MPa、 λ =0.58），显然，这个交点就是曲面 Π 1 和 Π 2 的公共零点。按照第3.1节的分析，此时的 τ 0 和 λ 就是所要求的解。此外，对于 λ =3/4和 λ =0的情况也应当进行校验。按照 式（30） 式（30） 和 式（31） 式（31） 分别画出曲线 L 1 和 L 2 的图像（如 图5 图5 （a）和 图5 图5 （b）所示），可以看到，对于这两种情况， L 1 和 L 2 都没有公共零点，因此不存在 τ 0 和 λ 的解。综合以上分析，将 τ 0 =5.52 MPa和 λ =0.58代入 式（27） 式（27） ，就得到 C 、 φ 的反演结果分别为： 1;2; <math display="block" id="M39"><mi>C</mi><mo>=</mo><mn mathvariant="normal">2.61</mn><mtext> </mtext><mi mathvariant="normal">M</mi><mi mathvariant="normal">P</mi><mi mathvariant="normal">a</mi><mo>,</mo><mi>φ</mi><mo>=</mo><mn mathvariant="normal">15.33</mn><mo>°</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465942&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465940&type= 42.58732986 3.89466691 . 对于实验中的土壤试样，Forrestal等人 ［ 25 25 ］ 给出了三轴实验下剪切强度 τ 与平均应力 p 之间的关系，如 图6 图6 ，此处，剪切强度 τ 定义为试样屈服时第一主应力 σ 1 与第三主应力 σ 3 的差 ［ 25 25 ］ （ τ = σ 1 -σ 3 ）。注意到空腔膨胀模型的球对称性，第一主应力 σ 1 就是空腔膨胀时的径向应力 σ r ，第三主应力 σ 3 就是空腔膨胀时的环向应力 σ θ （此时 σ 2 = σ 3 = σ θ ），于是Mohr-Coulomb准则（ 式（2） 式（2） ）可以写为： <math display="block" id="M40"><mi>τ</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">6</mn><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mi>φ</mi></mrow><mrow><mn mathvariant="normal">3</mn><mo>-</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mi>φ</mi></mrow></mfrac><mi>p</mi><mo>+</mo><mfrac><mrow><mn mathvariant="normal">6</mn><mi>C</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mi>φ</mi></mrow><mrow><mn mathvariant="normal">3</mn><mo>-</mo><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mi>φ</mi></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465968&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465965&type= 43.51866913 8.80533314 . 根据上式，对 图6 图6 中的数据进行线性拟合（如 图6 图6 中红线所示），可解得相应的 C 、 φ 值分别为 C =3.73 MPa、 φ =3.48°。 对比 式（36） 式（36） 所示的反演结果可见，对于内聚力 C 的反演结果与三轴实验的结果较为接近，而对内摩擦角 φ 的反演结果却与实验值差别较大。这一定程度上与采用的拟合方法有关。接下来，我们分别基于贯入加速度数据和贯入速度数据对阻力系数进行拟合，并进一步讨论拟合结果对方法的依赖。 3.3　基于加速度时程曲线和速度时程曲线进行的反演 2 贯入加速度与贯入阻力间仅相差 m 倍，原则上讲，由贯入加速度进行阻力系数 α s 、 β s 的拟合与采用贯入阻力进行拟合没有本质差别。但是，在利用阻力-速度关系进行拟合时，贯入速度是由加速度的数值积分获得，由于进行了数值积分操作，数据处理过程与直接利用加速度时程曲线进行拟合是不同的，因此拟合结果也将产生差别。此外，由 式（17） 式（17） 可见，贯入阻力 F z 与贯入速度的平方 V z 2 间呈线性关系，阻力系数 α s 和 β s 通过线性最小二乘拟合即可得到；而由 式（25） 式（25） 可见，贯入加速度 a z 与时间 t 呈现非线性关系， α s 和 β s 只能通过非线性最小二乘拟合得到。 基于加速度时程曲线 ［ 25 25 ］ 和速度时程曲线对阻力系数的拟合如 图7 图7 。根据贯入加速度表达 式（25） 式（25） ，采用Levenberg-Marquardt方法 ［ 30 30 ］ 对文献［ 25 25 ］中贯入实验3给出的加速度数据进行非线性最小二乘拟合，如 图7 图7 （a）中红线所示。拟合的初值设置为根据阻力-速度关系拟合得到的数值，即 α s =85.39 kN， β s =2.52 kg/m。由加速度曲线拟合得到的阻力系数分别为 α s =89.22 kN、 β s =2.69 kg/m。由此得到 τ 0 =5.61 MPa、 λ =0.72，Mohr-Coulomb参数为： <math display="block" id="M41"><mi>C</mi><mo>=</mo><mn mathvariant="normal">2.64</mn><mtext> </mtext><mi mathvariant="normal">M</mi><mi mathvariant="normal">P</mi><mi mathvariant="normal">a</mi><mo>,</mo><mi>φ</mi><mo>=</mo><mn mathvariant="normal">18.75</mn><mo>°</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465985&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41465981&type= 42.58732986 3.89466691 . 1,2; 与三轴实验测量值相比， C 值的相对误差为29.2%，与上一节中基于阻力-速度关系反演得到的结果接近， φ 值的反演误差仍然较大。此外， 图7 图7 （a）给出了按照文献［ 25 25 ］中取 τ 0 =10 MPa、 λ =0时计算得到的加速度时程曲线（ 图7 图7 （a）中蓝线），将其与红线相比较可以看到， τ 0 =10 MPa和 λ =0并不是测量数据的最小二乘解。 按照相同的非线性最小二乘拟合方法，利用贯入速度表达 式（24） 式（24） 亦可对阻力系数 α s 和 β s 进行拟合，如 图7 图7 （b）所示。此时，同样将拟合的初值设置为根据阻力-速度关系拟合得到的数值，即 α s =85.39 kN， β s =2.52 kg/m。由速度曲线拟合得到的阻力系数为 α s =101.48 kN、 β s =1.87 kg/m。相应的 τ 0 =7.76 MPa、 λ =0.11，Mohr-Coulomb参数为： <math display="block" id="M42"><mi>C</mi><mo>=</mo><mn mathvariant="normal">3.81</mn><mtext> </mtext><mi mathvariant="normal">M</mi><mi mathvariant="normal">P</mi><mi mathvariant="normal">a</mi><mo>,</mo><mi>φ</mi><mo>=</mo><mn mathvariant="normal">3.14</mn><mo>°</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41466013&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41466011&type= 40.72467041 3.89466691 . 将上式与三轴实验测量结果对比， C 值的相对误差为2.14%，较前两种方法的反演结果误差大大减小， φ 值的相对误差也减小为9.77%。 总结起来，对比 式（36） 式（36） 、 式（38） 式（38） 和 式（39） 式（39） 的反演结果可见，利用阻力-速度关系进行拟合的结果与利用加速度直接进行拟合的结果相近，而与利用贯入速度进行拟合的结果差别较大。更重要的是，相较于前两者的拟合结果，利用贯入速度拟合得到的内聚力 C 和内摩擦角 φ 与三轴实验给出的结果更加接近。以上结论并不是偶然的，它反映了贯入阻力系数对加速度、速度等物理量有不同的影响，这将在第4.2节中进行更进一步的讨论。此外，以上结论为土壤Mohr-Coulomb参数反演方法的实际应用提供了重要参照。 3.4　基于深度-初速度关系进行的反演 2 阻力系数 α s 、 β s 的数值还可以通过对同一状态的土壤进行多次贯入实验后测得的贯入深度 P 和初始贯入速度 V 0 来拟合。为此，根据 式（26） 式（26） 构造误差函数： <math display="block" id="M43"><msub><mrow><mi>Q</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>=</mo><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><msup><mrow><mfenced open="[" close="]" separators="|"><mrow><msub><mrow><mi>P</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>-</mo><mfrac><mrow><mi>m</mi></mrow><mrow><mn mathvariant="normal">2</mn><msub><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow></mfrac><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo>+</mo><msubsup><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow><mrow><mi>'</mi></mrow></msubsup><msubsup><mrow><mi>V</mi></mrow><mrow><mn mathvariant="normal">0</mn><mo>,</mo><mi>i</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo stretchy="false">)</mo></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41466016&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41466015&type= 53.93266296 10.58333302 ， 其中， N 是贯入实验次数。 <math display="block" id="M44"><msubsup><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow><mrow><mi>'</mi></mrow></msubsup><mo>=</mo><mrow><mrow><msub><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow><mo>/</mo><mrow><msub><mrow><mi>α</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41466022&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41466019&type= 14.90133286 4.14866638 . 令 <math display="block" id="M45"><mfrac><mrow><mo>∂</mo><msub><mrow><mi>Q</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow><mrow><mo>∂</mo><msub><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow></mfrac><mo>=</mo><mn mathvariant="normal">0</mn><mo>,</mo><mo> </mo><mfrac><mrow><mo>∂</mo><msub><mrow><mi>Q</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></mrow><mrow><mo>∂</mo><msubsup><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow><mrow><mi>'</mi></mrow></msubsup></mrow></mfrac><mo>=</mo><mn mathvariant="normal">0</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41466030&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41466026&type= 28.87133408 8.80533314 . 可解得： <math display="block" id="M46"><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><mfenced separators="|"><mrow><msub><mrow><mi>P</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>-</mo><mfrac><mrow><mi>m</mi></mrow><mrow><mn mathvariant="normal">2</mn><msub><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow></msub></mrow></mfrac><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo>+</mo><msubsup><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow><mrow><mi>'</mi></mrow></msubsup><msubsup><mrow><mi>V</mi></mrow><mrow><mn mathvariant="normal">0</mn><mo>,</mo><mi>i</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo stretchy="false">)</mo></mrow></mfenced><mfrac><mrow><msubsup><mrow><mi>V</mi></mrow><mrow><mn mathvariant="normal">0</mn><mo>,</mo><mi>i</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow><mrow><mn mathvariant="normal">1</mn><mo>+</mo><msubsup><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow><mrow><mi>'</mi></mrow></msubsup><msubsup><mrow><mi>V</mi></mrow><mrow><mn mathvariant="normal">0</mn><mo>,</mo><mi>i</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow></mfrac></mrow><mo>=</mo><mn mathvariant="normal">0</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41466034&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41466032&type= 66.97133636 9.90600014 ， <math display="block" id="M47"><msub><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow></msub><mo>=</mo><mfrac><mrow><mi>m</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mfrac><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><mi mathvariant="normal">l</mi><msup><mrow><mi mathvariant="normal">n</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo>+</mo><msubsup><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow><mrow><mi>'</mi></mrow></msubsup><msubsup><mrow><mi>V</mi></mrow><mrow><mn mathvariant="normal">0</mn><mo>,</mo><mi>i</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo stretchy="false">)</mo></mrow></mrow><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><msub><mrow><mi>P</mi></mrow><mrow><mi>i</mi></mrow></msub><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo>+</mo><msubsup><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow><mrow><mi>'</mi></mrow></msubsup><msubsup><mrow><mi>V</mi></mrow><mrow><mn mathvariant="normal">0</mn><mo>,</mo><mi>i</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo stretchy="false">)</mo></mrow></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41466038&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41466036&type= 43.34933472 18.54199982 . 将 式（44） 式（44） 代入 式（43） 式（43） 并令： <math display="block" id="M48"><mi>f</mi><mfenced separators="|"><mrow><msubsup><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow><mrow><mi>'</mi></mrow></msubsup></mrow></mfenced><mo>=</mo><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><mfenced separators="|"><mrow><msub><mrow><mi>P</mi></mrow><mrow><mi>i</mi></mrow></msub><mo>-</mo><mfrac><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><msub><mrow><mi>P</mi></mrow><mrow><mi>i</mi></mrow></msub><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo>+</mo><msubsup><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow><mrow><mi>'</mi></mrow></msubsup><msubsup><mrow><mi>V</mi></mrow><mrow><mn mathvariant="normal">0</mn><mo>,</mo><mi>i</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo stretchy="false">)</mo></mrow></mrow><mrow><mstyle displaystyle="true"><munderover><mo largeop="true">∑</mo><mrow><mi>i</mi><mo>=</mo><mn mathvariant="normal">1</mn></mrow><mrow><mi>N</mi></mrow></munderover></mstyle><mrow><mi mathvariant="normal">l</mi><msup><mrow><mi mathvariant="normal">n</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo>+</mo><msubsup><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow><mrow><mi>'</mi></mrow></msubsup><msubsup><mrow><mi>V</mi></mrow><mrow><mn mathvariant="normal">0</mn><mo>,</mo><mi>i</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo stretchy="false">)</mo></mrow></mrow></mfrac><mi mathvariant="normal">l</mi><mi mathvariant="normal">n</mi><mo stretchy="false">(</mo><mn mathvariant="normal">1</mn><mo>+</mo><msubsup><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow><mrow><mi>'</mi></mrow></msubsup><msubsup><mrow><mi>V</mi></mrow><mrow><mn mathvariant="normal">0</mn><mo>,</mo><mi>i</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup><mo stretchy="false">)</mo></mrow></mfenced><mfrac><mrow><msubsup><mrow><mi>V</mi></mrow><mrow><mn mathvariant="normal">0</mn><mo>,</mo><mi>i</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow><mrow><mn mathvariant="normal">1</mn><mo>+</mo><msubsup><mrow><mi>β</mi></mrow><mrow><mi>s</mi></mrow><mrow><mi>'</mi></mrow></msubsup><msubsup><mrow><mi>V</mi></mrow><mrow><mn mathvariant="normal">0</mn><mo>,</mo><mi>i</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow></mfrac></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41466042&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41466040&type= 101.09200287 19.72733307 . 函数 f （ β' s ）的零点即是方程（43）的解。求得 β' s 后，由 式（41） 式（41） 和 式（44） 式（44） 即可得到 α s 和 β s ，并进而算得 C 、 φ 的反演结果。下面以何翔等人 ［ 31 31 ］ 的弹靶实验为例，来具体说明这个计算过程。 实验中，弹体质量为286 g，弹头长度为140 mm，曲径比 ψ =8.54 ［ 31 31 - 32 32 31-32 ］ 。土壤的基本参数如 表2 表2 所示。 密度 ρ 0 /（g·cm -3 ） 弹性模量 E /MPa 内聚力 C /kPa 内摩擦角/（°） 1.86 36.0 61 30 表3 表3 中列出了同一刚性弹以不同初速度多次贯入同一状态土壤后测得的贯入深度结果。 图8 图8 为 f （ β' s ）的图像及贯入深度-贯入初速度关系。将 表3 表3 中的9个实验数据代入 式（45） 式（45） 可画出函数 f （ β' s ）的图像（ 图8 图8 （a）），并得到 f （ β' s ）的零点为 β' s =1.09×10 -3 s 2 /m 2 ，联立 式（41） 式（41） 和 式（44） 式（44） 可解得 α s =276.92 N， β s =0.30 kg/m。取锁定体积应变 η * =0.13，可得到 τ 0 =52.51 kPa， λ =1.16，Mohr-Coulomb参数为： 序号 初始贯入速度 V 0 /（m·s -1 ） 贯入深度 P /m 1 152.0 1.58 2 150.0 1.61 3 149.0 1.43 4 203.0 1.83 5 201.0 1.79 6 203.0 1.85 7 253.0 2.11 8 251.0 1.97 9 248.0 1.95 1;2; <math display="block" id="M49"><mi>C</mi><mo>=</mo><mn mathvariant="normal">25.17</mn><mtext> </mtext><mi mathvariant="normal">k</mi><mi mathvariant="normal">P</mi><mi mathvariant="normal">a</mi><mo>,</mo><mi>φ</mi><mo>=</mo><mn mathvariant="normal">29.08</mn><mo>°</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41466087&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=41466084&type= 42.92599869 3.89466691 . 反演得到的内聚力 C 的相对误差为58.74%，内摩擦角 φ 的相对误差为3.07%。可见，利用9个实验数据整体拟合反演得到的内聚力 C 与三轴实验测量结果相差较大，而得到的内摩擦角 φ 与三轴实验测量结果一致。 此外，对 表3 表3 中所列数据进行分组拟合。根据文献［ 31 31 ］， 表3 表3 中所列实验结果对应三种不同的贯入初速度 V 0 ，分别为150 m/s（序号1~3），200 m/s（序号4~6）和250 m/s（序号7~9），在每个速度等级下取1个数据点，总计3个数据点作为一组数据（总共27组），来拟合 α s 、 β s 的值，结果见 表4 表4 。分组拟合及反演得到的内聚力 C 的平均值及标准差为29.45±20.69 kPa，内摩擦角 φ 的平均值及标准差为25.00±17.88°。此时，分组拟合并反演得到的 C 的平均值较整体拟合结果有所提高，但相比61 kPa仍差别较大；而对于 φ 的反演，整体拟合结果与实验测得的30°更加接近。 分组 编号 对应表2的实验序号 反演得到的内聚力 C /kPa 反演得到的内摩擦角 φ /（°） 1 1，4，7 33.45 22.27 2 1，4，8 14.02 43.77 3 1，4，9 13.10 46.44 4 1，5，7 34.32 22.21 5 1，5，8 14.11 44.48 6 1，5，9 13.12 47.52 7 1，6，7 33.24 22.16 8 1，6，8 13.98 43.31 9 1，6，9 13.10 45.80 10 2，4，7 23.62 28.08 11 2，4，8 10.31 57.86 12 2，4，9 10.41 63.33 13 2，5，7 24.18 28.05 14 2，5，8 10.57 59.23 15 2，5，9 11.01 65.70 16 2，6，7 23.55 27.87 17 2，6，8 10.19 56.90 18 2，6，9 10.15 61.88 19 3，4，7 75.24 10.19 20 3，4，8 42.77 20.89 21 3，4，9 41.67 21.45 22 3，5，7 76.28 10.20 23 3，5，8 42.67 21.28 24 3，5，9 41.35 21.95 25 3，6，7 74.60 10.16 26 3，6，8 42.60 20.73 27 3，6，9 41.62 21.24 以上分析表明，利用贯入深度-初速度关系进行反演时，内摩擦角 φ 的反演结果更加可靠。 4 讨 论 1 4.1　参数的可解域 2 土壤Mohr-Coulomb参数 C 、 φ 并不是任意的，具有一定范围，相应的阻力系数 α s 和 β s 也对应了一定范围，称为参数的可解域。只有拟合得到的 α s 和 β s 落在可解域范围内，才能在动态空腔膨胀模型框架下进行 C 、 φ 值的反演。 图9 图9 （a）和 图9 图9 （b）分别给出了Forrestal ［ 25 25 ］ 和何翔 ［ 31 31 ］ 等人报道的实验结果对应的参数可解域。 图9 图9 （a）中给出了利用阻力-速度关系、加速度时程曲线和速度时程曲线拟合得到的 α s 、 β s 值，其均落在参数的可解域内。 图9 图9 （b）中的分组拟合值为第3.4节中利用贯入深度-初速度关系对27组数据拟合得到的 α s 、 β s 值；整体拟合值为利用9个实验数据拟合得到的 α s 、 β s 值。这些拟合结果也都在参数的可解域内，因此可以利用动态空腔膨胀理论来进行 C 、 φ 值的反演。 1;2; 实际的地质环境往往十分复杂，数据噪声也难以避免，根据测量数据拟合得到的 α s 和 β s 并不总能反演得到Mohr-Coulomb参数 C 、 φ 的数值。此时，就需要对偏离模型的过程进行更细致的分析 参数的可解域是土壤Mohr-Coulomb参数对阻力系数作出的限制，它反映了土壤物理力学特性对贯入阻力的影响。对于给定土壤，Mohr-Coulomb参数总是确定的，因此在理想情况下，阻力系数的实验值必定落在可解域内。但在实际当中，由于噪声的影响，或简化的模型未能涵盖贯入过程的复杂细节，阻力系数的拟合值可能超出参数的可解域。此时，就需要对冲击贯入的物理过程进行更加详细的分析，并对相应的数据处理手段进行优化。 4.2　拟合区间的选取及对反演结果的进一步讨论 2 由第3节的分析可以知道，要准确地反演 C 、 φ 的数值，关键是要准确地拟合阻力系数 α s 和 β s 。第3.1节至3.3节利用贯入时程曲线进行拟合时，拟合数据是稳定贯入阶段贯入阻力与贯入速度平方呈线性关系的部分（ 图2 图2 中 A 、 B 之间的数据）。数据拟合范围对阻力系数的拟合及Mohr-Coulomb参数的反演有重要影响，在这一节中就来对数据拟合范围的选取进行讨论。 为此选取拟合窗口起点为 图2 图2 中的点 A （对应时刻为 t =4.84 ms），改变拟合窗口终点 t end 进行阻力系数 α s 和 β s 的拟合以及Mohr-Coulomb参数的反演，结果如 图10 图10 所示。由图可见，不论拟合窗口的大小如何改变，基于阻力-速度关系反演的结果与基于加速度曲线反演的结果接近，而与基于速度曲线的反演有较大的差别；并且，对于内摩擦角 φ ，前两者的反演结果始终大于后者。 1;2; 图10 图10 （a）显示了两个关键拟合范围。第一个范围对应三条拟合曲线的交点（ t end = t 2 ）。对于内聚力 C 的反演，当 t end < t 2 时，基于阻力-速度关系和加速度曲线反演得到的结果小于利用速度曲线反演的结果；当 t end > t 2 时，基于阻力-速度关系和加速度曲线反演得到的结果大于利用速度曲线反演的结果。在 t end = t 2 的拟合范围下，三种拟合方法得到的 C 的反演结果一致。此时反演得到的数值（ C =4 MPa）比三轴实验结果（ C =3.73 MPa）大7.24%。另一个关键拟合范围对应于 图10 图10 （a）中速度曲线拟合结果的最小值点 t end = t 1 （这一点也是 图10 图10 （b）中速度曲线拟合结果的最大值点），此时，由速度曲线拟合并反演得到的内聚力 C =3.81 MPa，内摩擦角为 φ =3.14°，与三轴试验结果一致。这一范围就是在第3.1节至3.3节中拟合时所采用的范围， 图10 图10 中时刻 t end = t 1 就对应 图2 图2 和 图7 图7 中的点 B 。 以上分析为反演方法的实际应用提供了一个可行方案：对于冲击贯入实验测得的数据，用速度曲线进行参数 C （ φ ）的拟合，改变区间的大小进行搜索，使得 C （ φ ）的拟合值最小（大）时对应的区间即为最终的拟合区间，相应的 C 、 φ 值即为最终的反演结果。 由图 2 2 、 7 7 、 10 10 可以看到，拟合区间终点 t end 一般并不能覆盖整个贯入过程的数据段， 图10 图10 中能够反演得到Mohr-Coulomb参数的最大 t end 为27.80 ms，这是由于参数可解域的限制。当 t end 大于27.80 ms后，拟合得到的 α s 和 β s 已经超出了参数的可解域，这表明超过27.80 ms后的贯入过程可能存在模型忽略掉的细节。在利用动态空腔膨胀模型进行参数反演时，假定了土壤介质的初始状态是均匀的，这一假设保证了模型的球对称性，而在现场贯入实验中，土壤实际上是分层的。 图11 图11 给出了文献［ 25 25 ］中贯入实验3的加速度时程曲线并根据贯入仪进入不同土层（Ⅰ、Ⅱ、Ⅲ）的时刻进行了划分。由第3.2节的说明可知，这三层土壤的弹性模量和锁定体积应变是不同的，特别是第三层土壤的弹性模量明显高于前两层；此外，第三层土层的厚度也几乎是前两层的两倍，贯入仪大部分的运动时间在第三层土壤中。这表明，第三层土壤对贯入过程的贡献应当比前两层土壤的大，然而平均化的假定抹平了这些差别。当拟合区间为 图11 图11 中的 AB 范围时，贯入仪在三个土层中的运动时间相差不大，此时，平均化的假定基本适用；而当拟合区间范围逐渐扩大，比如对于 AC 段或 AD 段，此时第三个土层的贡献逐渐增大（但仍与前两层可比拟），平均化的假定将产生误差，此时可能需要考虑土壤分层的影响。总之，准确地反演土壤的力学参数，需要模型能够准确地把握贯入过程的关键细节。 4.3　参数的敏感性 2 在求解冲击贯入问题时，对于贯入阻力和贯入深度的求解属于正问题。此时，弹性模量和锁定体积应变等参数对于计算结果的影响不大 ［ 25 25 ］ 。然而，对于通过测量贯入阻力求解靶体介质力学特性参数这样的反问题，反演结果往往具有参数敏感性。这种敏感性不仅与材料参数本身有关还与阻力系数有关。 图12 图12 和 图13 图13 分别给出了文献［ 25 25 ］和文献［ 31 31 ］中，土壤弹性模量 E 和锁定体积应变 η * 对 C 、 φ 值反演的影响。由图可见，随着弹性模量 E 的增大，反演得到的 C 值减小、 φ 值增大；随着锁定体积应变 η * 的增大，反演得到的 C 、 φ 值均增大。 1;2; 1;2; 由 图12 图12 （a）和 图13 图13 （a）可见，弹性模量 E 对于 C 、 φ 值反演结果的影响较小。在 图12 图12 （a）中，当弹性模量 E 的变化范围为135 MPa时， C 值的变化大约仅有1 MPa，而 φ 值的变化约3°。在 图13 图13 （a）中，当弹性模量 E 的变化范围为45 MPa时， C 值的变化不到2 kPa，而 φ 值的变化小于1°。这就是说， C 、 φ 值的反演结果对弹性模量 E 不敏感。相反地，由 图12 图12 （b）和 图13 图13 （b）可见，锁定体积应变 η * 对于 C 、 φ 值的反演结果却有很大影响，即 C 、 φ 值的反演结果对参数 η * 是敏感的。 此外，从 图12 图12 （b）和 图13 图13 （b）还可以看到，不同阻力系数导致根据相同范围内的 η * 反演得到的 C 、 φ 值也不相同。这说明反演结果的参数敏感性也与阻力系数 α s 、 β s 有关。 5 结 论 1 本文基于动态空腔膨胀模型，利用土壤冲击贯入实验中测得的阻力-速度关系、加速度时程曲线、速度时程曲线，以及贯入深度-初速度关系，建立了对土壤Mohr-Coulomb参数（内聚力 C 和内摩擦角 φ ）进行反演的方法。其中，基于前三者建立的反演方法是利用单次贯入过程中的时程曲线进行的反演；而基于贯入深度-初速度关系建立的反演方法则是在对同一状态的土壤进行多次贯入时采用。 在贯入过程中的稳定贯入阶段，贯入阻力与贯入速度平方之间呈线性关系，准确地拟合线性系数（阻力系数） α s 和 β s 是反演Mohr-Coulomb参数的关键。本文揭示了利用动态空腔膨胀模型进行Mohr-Coulomb参数反演时存在一个关于参数的可解域。对于测得的实验数据，只有首先保证拟合得到的阻力系数落在可解域范围内， C 、 φ 值的反演才能在动态空腔膨胀模型的框架下进行。 本文分析表明利用阻力-速度关系进行反演的结果与利用加速度时程曲线进行反演的结果相近，而这两者与利用速度时程曲线进行反演的结果差别较大。这本质上是由于阻力系数 α s 、 β s 分别对加速度和速度的贡献不同。因此，这三种反演方法对 C 、 φ 值反演的准确性也是不同的。其中，利用速度时程曲线进行参数反演的结果最为准确。利用Forrestal等人 ［ 25 25 ］ 报道的土壤原位贯入实验结果进行的方法验证表明，利用速度时程曲线对土壤内聚力反演的相对误差为2.14%，对土壤内摩擦角反演的相对误差为9.77%。 此外，Mohr-Coulomb参数的反演具有参数敏感性，这不仅与材料参数本身有关还与阻力系数 α s 和 β s 有关。 C 、 φ 值的反演结果对锁定体积应变 η * 敏感，而对弹性模量 E 不敏感。 本文建立的反演方法可为复杂地质环境下土壤Mohr-Coulomb参数的快速确定提供基础，并为大范围内的地质力学信息勘测提供一个新的途径。 参考文献 WASFY T ， JAYAKUMAR P . Next-generation NATO reference mobility model complex terramechanics - Part 1： definition and literature review ［J］. Journal of Terramechanics ， 2021 ， 96 ： 45 - 57 . doi: 10.1016/j.jterra.2021.02.002 http://dx.doi.org/10.1016/j.jterra.2021.02.002 Next-generation NATO reference mobility model complex terramechanics - Part 1： definition and literature review WASFY T M ， JAYAKUMAR P . Next-generation NATO reference mobility model complex terramechanics - Part 2： requirements and prototype ［J］. Journal of Terramechanics ， 2021 ， 96 ： 59 - 79 . doi: 10.1016/j.jterra.2021.02.007 http://dx.doi.org/10.1016/j.jterra.2021.02.007 Next-generation NATO reference mobility model complex terramechanics - Part 2： requirements and prototype 陈薪硕 ， 李守定 ， 张晓静 ， 等 . 行星地质工程原位测试方法 ［J］. 工程地质学报 ， 2021 ， 29 （ 5 ）： 1525 - 1544 . 行星地质工程原位测试方法 CHEN X S ， LI S D ， ZHANG X J ， et al . In-situ testing method for planetary geological engineering ［J］. Journal of Engineering Geology ， 2021 ， 29 （ 5 ）： 1525 - 1544 . （in Chinese） In-situ testing method for planetary geological engineering 华琛 ， 牛润新 ， 余彪 . 地面车辆机动性评估方法与应用 ［J］. 吉林大学学报（工学版） ， 2022 ， 52 （ 6 ）： 1229 - 1244 . 地面车辆机动性评估方法与应用 HUA C ， NIU R X ， YU B . Methods and applications of ground vehicle mobility evaluation ［J］. Journal of Jilin University （Engineering and Technology Edition） ， 2022 ， 52 （ 6 ）： 1229 - 1244 . （in Chinese） Methods and applications of ground vehicle mobility evaluation 张克健 . 车辆地面力学 ［M］. 北京 ： 国防工业出版社 ， 2002 . 车辆地面力学 ZHANG K J . Vehicle-Terramechanics ［M］. Beijing ： National Defense Industry Press ， 2002 . （in Chinese） Vehicle-Terramechanics 何健 ， 马吉胜 ， 吴大林 . 车辆地面力学研究综述 ［J］. 火炮发射与控制学报 ， 2017 ， 38 （ 2 ）： 89 - 93 . 车辆地面力学研究综述 HE J ， MA J S ， WU D L . A research review of vehicle terramechanics ［J］. Journal of Gun Launch ＆ Control ， 2017 ， 38 （ 2 ）： 89 - 93 . （in Chinese） A research review of vehicle terramechanics ROSKODEN R . Processing advancements of free fall penetrometer data and their ground proving in regional studies in New Zealand ［D］. German ： University of Bremen ， 2020 . Processing advancements of free fall penetrometer data and their ground proving in regional studies in New Zealand 孟醒 ， 王栋 . 自由落体贯入试验获得土体强度的简化分析 ［J］. 土木与环境工程学报（中英文） ， 2019 ， 41 （ 1 ）： 22 - 28 . 自由落体贯入试验获得土体强度的简化分析 MENG X ， WANG D . Simplified analysis of determining undrained shear strength from FFP tests ［J］. Journal of Civil and Environmental Engineering ， 2019 ， 41 （ 1 ）： 22 - 28 . （in Chinese） Simplified analysis of determining undrained shear strength from FFP tests CHOW S H ， O'LOUGHLIN C D ， WHITE D J ， et al . An extended interpretation of the free-fall piezocone test in clay ［J］. Géotechnique ， 2017 ： 1 - 14 . doi: 10.1680/jgeot.16.p.220 http://dx.doi.org/10.1680/jgeot.16.p.220 An extended interpretation of the free-fall piezocone test in clay AUBENY C P ， SHI H . Interpretation of impact penetration measurements in soft clays ［J］. Journal of Geotechnical and Geoenvironmental Engineering ， 2006 ， 132 （ 6 ）： 770 - 777 . doi: 10.1061/(asce)1090-0241(2006)132:6(770) http://dx.doi.org/10.1061/(asce)1090-0241(2006)132:6(770) Interpretation of impact penetration measurements in soft clays AKAL T ， STOLL R D . A dynamic penetrometer for rapid assessment of seafloor parameters ［J］. The Journal of the Acoustical Society of America ， 1994 ， 96 （ 5 ）： 3219 . doi: 10.1121/1.411178 http://dx.doi.org/10.1121/1.411178 A dynamic penetrometer for rapid assessment of seafloor parameters BENDOR G ， DUBINSKY A ， ELPERIN T . Ballistic impact： recent advances in analytical modeling of plate penetration dynamics–A review ［J］. Applied Mechanics Reviews ， 2005 ， 58 （ 6 ）： 355 - 371 . doi: 10.1115/1.2048626 http://dx.doi.org/10.1115/1.2048626 Ballistic impact： recent advances in analytical modeling of plate penetration dynamics–A review BOGUSLAVSKII Y ， DRABKIN S ， SALMAN A . Analysis of vertical projectile penetration in granular soils ［J］. Journal of Physics D： Applied Physics ， 1996 ， 29 （ 3 ）： 905 - 916 . doi: 10.1088/0022-3727/29/3/059 http://dx.doi.org/10.1088/0022-3727/29/3/059 Analysis of vertical projectile penetration in granular soils 陈小伟 . 穿甲/侵彻问题的若干工程研究进展 ［J］. 力学进展 ， 2009 ， 39 （ 3 ）： 316 - 351 . doi: 10.6052/1000-0992-2009-3-J2007-090 http://dx.doi.org/10.6052/1000-0992-2009-3-J2007-090 穿甲/侵彻问题的若干工程研究进展 CHEN X W . Advances in the penetration/perforation of rigid projectiles ［J］. Advances in Mechanics ， 2009 ， 39 （ 3 ）： 316 - 351 . （in Chinese） . doi: 10.6052/1000-0992-2009-3-J2007-090 http://dx.doi.org/10.6052/1000-0992-2009-3-J2007-090 Advances in the penetration/perforation of rigid projectiles TRUE D . Undrained Vertical Penetration Into ocean Bottom Soils ［D］. Berkeley ， USA ， University of California ， 1976 Undrained Vertical Penetration Into ocean Bottom Soils 刘君 ， 张雨勤 . FFP在黏土中贯入过程的CFD模拟 ［J］. 力学学报 ， 2018 ， 50 （ 1 ）： 167 - 176 . doi: 10.6052/0459-1879-17-284 http://dx.doi.org/10.6052/0459-1879-17-284 FFP在黏土中贯入过程的CFD模拟 LIU J ， ZHANG Y Q . Cfd simulation on the penetration of ffp into uniform clay ［J］. Chinese Journal of Theoretical and Applied Mechanics ， 2018 ， 50 （ 1 ）： 167 - 176 . （in Chinese） . doi: 10.6052/0459-1879-17-284 http://dx.doi.org/10.6052/0459-1879-17-284 Cfd simulation on the penetration of ffp into uniform clay KARL T . Theoretical Soil Mechanics ［M］. Hoboken， NJ， USA ： John Wiley & Sons， Inc. ， 1943 . Theoretical Soil Mechanics BALIGH M M . Strain path method ［J］. Journal of Geotechnical Engineering ， 1985 ， 111 （ 9 ）： 1108 - 1136 . doi: 10.1061/(asce)0733-9410(1985)111:9(1108) http://dx.doi.org/10.1061/(asce)0733-9410(1985)111:9(1108) Strain path method BISHOP R F ， HILL R ， MOTT N F . The theory of indentation and hardness tests ［J］. Proceedings of the Physical Society ， 1945 ， 57 （ 3 ）： 147 - 159 . doi: 10.1088/0959-5309/57/3/301 http://dx.doi.org/10.1088/0959-5309/57/3/301 The theory of indentation and hardness tests CHADWICK P . The quasi-static expansion of a spherical cavity in metals and ideal soils ［J］. The Quarterly Journal of Mechanics and Applied Mathematics ， 1959 ， 12 （ 1 ）： 52 - 71 . doi: 10.1093/qjmam/12.1.52 http://dx.doi.org/10.1093/qjmam/12.1.52 The quasi-static expansion of a spherical cavity in metals and ideal soils CAO L F ， TEH C I ， CHANG M F . Analysis of undrained cavity expansion in elasto-plastic soils with non-linear elasticity ［J］. International Journal for Numerical and Analytical Methods in Geomechanics ， 2002 ， 26 （ 1 ）： 25 - 52 . doi: 10.1002/nag.189 http://dx.doi.org/10.1002/nag.189 Analysis of undrained cavity expansion in elasto-plastic soils with non-linear elasticity HILL R . The Mathematical Theory of Plasticity ［M］. Oxford University Press ， 1998 . The Mathematical Theory of Plasticity FORRESTAL M J . Penetration into dry porous rock ［J］. International Journal of Solids and Structures ， 1986 ， 22 （ 12 ）： 1485 - 1500 . doi: 10.1016/0020-7683(86)90057-0 http://dx.doi.org/10.1016/0020-7683(86)90057-0 Penetration into dry porous rock LUK V K ， FORRESTAL M J . Penetration into semi-infinite reinforced-concrete targets with spherical and ogival nose projectiles ［J］. International Journal of Impact Engineering ， 1987 ， 6 （ 4 ）： 291 - 301 . doi: 10.1016/0734-743x(87)90096-0 http://dx.doi.org/10.1016/0734-743x(87)90096-0 Penetration into semi-infinite reinforced-concrete targets with spherical and ogival nose projectiles FORRESTAL M J ， LUK V K . Penetration into soil targets ［J］. International Journal of Impact Engineering ， 1992 ， 12 （ 3 ）： 427 - 444 . doi: 10.1016/0734-743x(92)90167-r http://dx.doi.org/10.1016/0734-743x(92)90167-r Penetration into soil targets SHI C C ， WANG M Y ， LI J ， et al . A model of depth calculation for projectile penetration into dry sand and comparison with experiments ［J］. International Journal of Impact Engineering ， 2014 ， 73 ： 112 - 122 . doi: 10.1016/j.ijimpeng.2014.06.010 http://dx.doi.org/10.1016/j.ijimpeng.2014.06.010 A model of depth calculation for projectile penetration into dry sand and comparison with experiments JANBU N ， SENNESET K . Effective stress interpretation of in situ static penetration tests ［C］. Proceedings of the 1st European symposium on penetration testing ， ESOPT， Stockholm ， 1974 ， 2 ： 181 - 93 . Effective stress interpretation of in situ static penetration tests 钱七虎 ， 王明洋 . 高等防护结构计算理论 ［M］. 南京 ： 江苏科学技术出版社 ， 2009 . 高等防护结构计算理论 QIAN Q H ， WANG M Y . Calculation theory for advanced protective structure ［M］. Nanjing ： Phoenix Science Press ， 2009 . （in Chinese） Calculation theory for advanced protective structure 吴华杰 ， 王明洋 ， 吴春梅 ， 等 . 钻地武器侵彻岩石滑动摩擦系数影响及计算研究 ［J］. 岩土工程学报 ， 2009 ， 31 （ 12 ）： 1845 - 1850 . doi: 10.3321/j.issn:1000-4548.2009.12.007 http://dx.doi.org/10.3321/j.issn:1000-4548.2009.12.007 钻地武器侵彻岩石滑动摩擦系数影响及计算研究 WU H J ， WANG M Y ， WU C M ， et al . Influence and calculation of friction coefficient of sliding of projectile penetration into rock ［J］. Chinese Journal of Geotechnical Engineering ， 2009 ， 31 （ 12 ）： 1845 - 1850 . （in Chinese） . doi: 10.3321/j.issn:1000-4548.2009.12.007 http://dx.doi.org/10.3321/j.issn:1000-4548.2009.12.007 Influence and calculation of friction coefficient of sliding of projectile penetration into rock MARQUARDT D W . An algorithm for least-squares estimation of nonlinear parameters ［J］. Journal of the Society for Industrial and Applied Mathematics ， 1963 ， 11 （ 2 ）： 431 - 441 . doi: 10.1137/0111030 http://dx.doi.org/10.1137/0111030 An algorithm for least-squares estimation of nonlinear parameters 何翔 ， 吴祥云 ， 曲建波 . 细长弹丸侵彻土壤介质的实验研究 ［C］. 第二届全国爆炸力学实验技术交流会论文集，安徽，中国 ： 安徽省力学学会 ， 2002 ： 92 - 100 . 细长弹丸侵彻土壤介质的实验研究 HE X ， WU X Y ， QU J B . Experimental study on the penetration of slender projectile into the soil medium ［C］. Proceedings of the Second National Explosion Mechanics Experimental Technology Exchange Conference， Anhui， China ： ASTAM ， 2002 ： 92 - 100 . （in Chinese） Experimental study on the penetration of slender projectile into the soil medium 武海军 ， 曾费隐 ， 李金柱 ， 等 . 考虑含水率干密度的土壤动态球形空腔膨胀理论 ［J］. 北京理工大学学报 ， 2021 ， 41 （ 5 ）： 459 - 464 . 考虑含水率干密度的土壤动态球形空腔膨胀理论 WU H J ， ZENG F Y ， LI J Z ， et al . Dynamic spherical cavity expansion theory of soil considering moisture content and dry density ［J］. Transactions of Beijing Institute of Technology ， 2021 ， 41 （ 5 ）： 459 - 464 . （in Chinese） Dynamic spherical cavity expansion theory of soil considering moisture content and dry density 10.37188/OPE.20233105.0565 TP391.4 A 1004-924X（2023）05-0565-16 邹欧,薛健,李娜等.基于动态空腔膨胀模型的土壤Mohr-Coulomb参数反演方法[J].光学精密工程,2023,31(05):565-580. 邹欧,薛健,李娜等.基于动态空腔膨胀模型的土壤Mohr-Coulomb参数反演方法[J].光学精密工程,2023,31(05):565-580. ZOU Ou,XUE Jian,LI Na,et al.Inversion method of soil Mohr-Coulomb parameters based on dynamic cavity expansion model[J].Optics and Precision Engineering,2023,31(05):565-580. ZOU Ou,XUE Jian,LI Na,et al.Inversion method of soil Mohr-Coulomb parameters based on dynamic cavity expansion model[J].Optics and Precision Engineering,2023,31(05):565-580. 中国科学院青年创新促进会资助项目（No.2018024）；国家自然科学基金资助项目（No.61575209）；钱学森实验室太空探索实验培育项目（No.TKTSPY-2020-05-01） 2022-08-26 2022-10-14 2023-03-10 2023-03-22 光学精密工程 5 31