冻融循环下黄土力学性质损伤的试验解析与机制探究

冻融循环下黄土力学性质损伤的试验解析与机制探究一、引言1.1研究背景与意义1.1.1研究背景黄土作为一种在全球广泛分布的特殊土类，主要集中于中纬度干旱或半干旱的大陆性气候地区，涵盖了中欧、北美以及亚洲中部等地。在我国，黄土主要分布于西北和华北地区，总面积约63.5万平方千米，其中原生黄土为38.1万平方千米，次生黄土为25.4万平方千米。黄河中下游的陕西北部、甘肃中部和东部、宁夏南部和山西西部是黄土分布最为集中的区域，这些地区地势较高，共同构成了世界闻名的黄土高原，其黄土厚度最厚可达200m。黄土的颗粒成分以粉土粒级为主，含量超过50%，富含碳酸钙质、钙质结核，有时还含有硫酸岩或氯化物盐类，呈现出肉眼可见孔隙的第四纪陆相沉积物特征。在季节性冻土地区，黄土会受到冻融循环现象的显著影响。随着外界气候的季节性交替变化，土体反复经历冻胀和融沉过程。在冬季，温度下降致使黄土中的水分冻结，水变成冰后体积膨胀，进而导致土体膨胀；而在春季，气温回升，土壤中的冰逐渐融化成水，土体体积收缩。这种周期性的冻融循环过程对黄土的物理力学性质产生了至关重要的影响，如导致黄土内部孔隙结构发生变化、毛细管压力改变，进而引起黄土的强度、变形模量等力学特性改变。在实际工程中，冻融循环作用引发的黄土力学性质变化给工程建设带来了诸多问题。在黄土地区修建的建筑物地基，由于冻融循环导致黄土承载力下降，可能引发建筑物的不均匀沉降、开裂甚至倒塌；黄土边坡在冻融循环作用下，稳定性降低，容易发生滑坡、坍塌等地质灾害，严重威胁到工程安全和人民生命财产安全。在道路工程中，黄土路基受冻融循环影响，会出现路面开裂、鼓包、翻浆等病害，不仅增加了道路维护成本，还影响了道路的使用寿命和行车安全。因此，深入研究冻融循环作用对黄土力学性质的影响具有重要的现实意义。1.1.2研究意义本研究具有重要的理论和实践意义。在理论方面，有助于深入理解黄土在冻融循环作用下的物理力学行为，丰富和完善土力学理论体系。通过探究冻融循环对黄土基本物理性质（如含水率、密度、孔隙比等）、力学性质（如抗剪强度、压缩性、渗透性等）以及微观结构的影响机制，能够揭示黄土在复杂环境条件下的力学响应规律，为黄土力学的进一步发展提供理论支持。在实践方面，研究成果可为黄土地区的工程建设提供科学依据和技术指导。在地基处理工程中，依据研究结果可以合理选择地基处理方法和参数，提高地基的稳定性和承载能力，减少建筑物的不均匀沉降；在边坡防护工程中，能够为边坡的设计和加固提供理论依据，增强边坡在冻融循环条件下的稳定性，预防滑坡等地质灾害的发生；在道路工程中，有助于优化路基设计和施工工艺，提高道路的耐久性和行车舒适性，降低道路维护成本。此外，本研究对于黄土地区的水利工程、隧道工程等其他基础设施建设也具有重要的参考价值，能够有效保障工程的安全运行和长期稳定性，促进黄土地区的经济发展和社会稳定。1.2国内外研究现状国内外学者针对冻融循环作用对黄土力学性质的影响开展了大量研究。在黄土的强度特性方面，众多研究表明冻融循环对黄土的抗剪强度影响显著。有学者通过直剪试验发现，随着冻融循环次数的增加，黄土的抗剪强度呈现先增加后减小的趋势。在冻融循环初期，由于水分结冰产生的冰晶对黄土颗粒间的连接有一定的强化作用，使得颗粒间的连接更加紧密，抗剪强度有所增加；但随着冻融循环次数的进一步增多，冰晶的反复膨胀和融化导致黄土结构逐渐破坏，颗粒间的摩擦力和黏聚力逐渐减弱，抗剪强度逐渐减小。还有学者通过三轴剪切试验研究了冻融黄土黏聚力和内摩擦角与冷端温度、融化温度的关系，发现黏聚力的大小与冷端温度、融化温度有关，而内摩擦角未呈现出规律性的波动。在变形特性研究中，学者们借助一维压缩试验、动三轴试验等手段，分析了冻融循环对黄土压缩性和动力变形特性的影响。研究发现，经过冻融循环之后，黄土的压缩系数增大，表明其压缩性增强。在动力变形方面，随着冻融循环次数的增加，黄土的动弹性模量和动剪切模量逐渐减小，而动阻尼比逐渐增大，这意味着黄土在冻融循环作用下，对地震等动力荷载的响应能力降低。在微观结构研究领域，扫描电镜（SEM）、压汞仪（MIP）和计算机断层扫描（CT）等技术被广泛应用。通过SEM观察发现，在冻融初期，黄土颗粒间的连接变得更加紧密，形成了更为稳固的结构形态；然而，随着冻融循环次数的增加，黄土的微观结构逐渐变得松散，颗粒间的连接逐渐减弱，甚至出现了一些裂缝和孔隙。MIP测试结果表明，冻融循环会导致黄土孔隙结构发生变化，孔隙尺寸分布和孔隙连通性改变。CT扫描技术则能够直观地呈现黄土内部结构在冻融循环过程中的三维变化情况，为深入研究微观结构演变提供了有力支持。尽管国内外学者在该领域已取得了一定的研究成果，但仍存在一些不足之处。现有研究大多集中在特定地区的黄土，对于不同区域黄土的对比研究较少；研究因素相对单一，综合考虑多种因素（如温度、含水率、应力状态等）对黄土力学性质影响的研究还不够深入；在微观结构与宏观力学性质的定量关系研究方面，还需要进一步加强，以更准确地揭示冻融循环作用下黄土力学性质变化的内在机制。1.3研究内容与方法1.3.1研究内容本研究将从多个方面深入探讨冻融循环作用对黄土力学性质的影响。首先，对黄土的基本物理性质展开研究，分析冻融循环过程中黄土的含水率、密度、孔隙比等物理指标的变化规律，明确冻融循环对黄土物理状态的影响。其次，重点研究冻融循环对黄土力学性质的影响，通过室内试验，包括直剪试验、三轴剪切试验、压缩试验等，获取不同冻融循环次数下黄土的抗剪强度、黏聚力、内摩擦角、压缩系数、变形模量等力学参数，分析其变化趋势和影响因素。再者，运用微观测试技术，如扫描电镜（SEM）、压汞仪（MIP）等，研究冻融循环作用下黄土微观结构的变化特征，包括颗粒排列方式、孔隙结构、颗粒间连接等，揭示微观结构变化与宏观力学性质之间的内在联系。此外，基于试验结果，深入探讨冻融循环作用下黄土力学性质损伤的机制，分析水分迁移、冰晶生长、结构破坏等因素对力学性质损伤的影响过程和作用机理。最后，采用数值模拟方法，建立黄土在冻融循环作用下的力学模型，对黄土的力学行为进行数值模拟分析，验证试验结果的准确性，预测不同条件下黄土力学性质的变化，为工程实践提供更具前瞻性的指导。1.3.2研究方法本研究采用多种研究方法相结合，以确保研究的全面性和准确性。在野外取样方面，选择具有代表性的黄土地区进行现场取样，保证所取黄土样品能够真实反映该地区黄土的特性。对采集的样品进行妥善保存和运输，避免样品受到扰动和污染，为后续的室内试验提供可靠的材料。室内试验是本研究的重要方法之一。通过冻融循环试验，模拟自然环境中的冻融过程，控制冷冻温度、融化温度和冻融循环次数等参数，对黄土样品进行周期性的冻融处理。利用直剪试验、三轴剪切试验等力学试验，测定不同冻融循环次数下黄土的抗剪强度等力学参数；通过压缩试验，获取黄土的压缩性指标。在试验过程中，严格按照相关试验标准和操作规程进行操作，确保试验数据的准确性和可靠性。微观测试方法用于研究黄土的微观结构。运用扫描电镜（SEM）对不同冻融循环次数下的黄土样品进行微观形貌观察，分析颗粒的大小、形状、排列方式以及颗粒间的连接情况；采用压汞仪（MIP）测试黄土的孔隙结构参数，包括孔隙大小分布、孔隙体积、孔隙率等，深入了解冻融循环对黄土微观孔隙结构的影响。数值模拟方法则借助专业的岩土工程数值分析软件，建立黄土在冻融循环作用下的力学模型。根据试验获得的物理力学参数，对模型进行参数赋值，模拟不同工况下黄土的力学响应，如应力分布、应变发展、变形情况等。通过数值模拟，不仅能够对试验结果进行验证和补充，还可以预测在复杂条件下黄土的力学行为，为工程设计和分析提供更全面的依据。二、黄土基本物理性质与冻融循环作用概述2.1黄土基本物理性质2.1.1黄土的物质组成黄土的物质组成是决定其物理力学性质的关键因素之一，主要包括矿物成分和颗粒组成。在矿物成分方面，黄土的化学成分主要由SiO₂、Al₂O₃、CaO等组成。其中，SiO₂含量约占50%左右，这主要源于黄土矿物成分中大量的石英。石英硬度高、化学性质稳定，使得黄土具有一定的抗风化能力。Al₂O₃含量在8-15%，主要与黄土中的长石等矿物相关。CaO含量约为10%左右，黄土中含有较多的碳酸盐类矿物，如方解石，是CaO的主要来源。此外，Fe₂O₃含量为4-5%，MgO含量在2-3%，K₂O含量约2%。Fe₂O₃和MgO的含量与黄土中含有的带铁锰成分的不稳定矿物有关，如角闪石、辉石、黑云母等。这些矿物成分相互作用，共同影响着黄土的性质。例如，碳酸盐类矿物在黄土中常起到胶结作用，使黄土颗粒在天然状态下呈团聚体存在，对黄土的结构强度有一定贡献。但当黄土遇水时，可溶盐胶结的团聚体被破坏，可能导致黄土湿陷性增强。从颗粒组成来看，黄土以粉土粒级为主，含量超过50%。黄土的颗粒组成具有较高的均一性，通常可将其颗粒分为5级（以毫米为单位）：0.25-0.1为粗细沙，0.1-0.05为细沙，0.05-0.01为粗粉沙，0.01-0.005为粉沙土，＜0.005为粘土。其中，粗粉粒（0.05-0.01mm）含量在50%以上。这种颗粒组成特征使得黄土具有独特的物理力学性质。粉土颗粒之间的相互作用较弱，导致黄土的孔隙率较大，结构相对疏松，具有较好的透水性，但抗剪强度相对较低。而一定比例的细沙和黏土颗粒的存在，又在一定程度上影响着黄土的黏聚性和压缩性。细沙颗粒可以增加黄土的骨架支撑作用，提高其承载能力；黏土颗粒则具有较强的吸附性和可塑性，能够影响黄土的黏性和变形特性。2.1.2黄土的天然结构特性黄土的天然结构特性对其力学性质有着重要影响，主要体现在孔隙结构和颗粒排列方式等方面。黄土具有独特的孔隙结构，其孔隙度一般在40%-50%，孔隙大小分布较为复杂，包含了大孔隙、中孔隙和小孔隙。大孔隙主要是由黄土颗粒间的架空形成，对黄土的渗透性和透气性影响较大；中孔隙和小孔隙则主要存在于颗粒团聚体内部和颗粒之间的接触部位，对黄土的持水性和力学性质有重要作用。黄土中的孔隙呈现出不规则的形状，且孔隙之间的连通性较好，这使得水分和气体在黄土中能够较为顺畅地运移。在颗粒排列方式上，黄土颗粒之间交错排列，形成了一种多孔结构。黄土中的黏土矿物和细粒矿物相互粘结，将较大的颗粒连接在一起，进一步增强了黄土结构的稳定性。然而，这种结构也具有一定的脆弱性，当受到外界因素如冻融循环、荷载作用等影响时，颗粒间的连接可能会被破坏，导致黄土结构的改变。在冻融循环过程中，水分冻结膨胀产生的冰晶压力会使颗粒间的连接受到破坏，孔隙结构发生变化，进而影响黄土的力学性质。黄土中的氧化铁颗粒和碳酸盐颗粒存在于颗粒表面，并粘结在一起，这也为黄土的结构特性提供了一定的支持。但在某些情况下，如遇水或化学作用，这些粘结物质可能会被溶解或破坏，从而降低黄土的结构强度。2.1.3黄土的初始物理参数测定为了深入研究冻融循环作用对黄土力学性质的影响，需要对黄土的初始物理参数进行准确测定。本研究通过一系列室内试验，测定了黄土的初始干密度、含水率、液塑限等物理参数。采用环刀法测定黄土的初始干密度。具体操作如下：在环刀内壁涂一层薄薄的凡士林油，将其刃口向下放在试样上。用修土刀沿环刀外缘将土样削成略大于环刀直径的土柱，然后慢慢将环刀垂直下压，边压边削，直到土样伸出环刀为止。用刮土刀仔细刮平两端余土，注意刮平时不得使土样扰动或压密。擦净环刀外壁，称量环刀加土的质量，准确至0.1g。通过公式计算得出黄土的初始干密度，本次试验测定的黄土初始干密度为[X]g/cm³。对于含水率的测定，采用烘干法这一室内试验的标准方法。取代表性土样15-30g，放入称量盒内，盖好盒盖，放在天平上称量，准确至0.01g。揭开盒盖，套在盒底，放入烘箱，在温度105-110°C下烘至质量恒定。将烘干后的土样，盖好盒盖收入干燥器内冷却至室温，称干土质量，准确至0.01g。通过公式计算含水率，经测定，黄土的初始含水率为[X]%。液塑限的测定采用液塑限联合测定法。取0.5mm筛下的代表性土样200g，分开放入三个盛土杯皿中，加入不同数量的蒸馏水，使土样的含水率分别控制在液限（a点）、略大于塑限（c点）和两者的中间状态（b点），用调土刀调匀，放置18h以上。将备好的土样搅拌均匀，分层装入盛土杯，压密，使空气逸出。将装好土样的试杯放在联合测定仪的升降座上，待锥尖与土样刚好接触，锥体下落，5s时读取锥入深度h1。改变锥尖与土的接触位置，测得锥入深度h2，h1与h2允许误差0.5mm，否则应重做，取h1、h2平均值作为该点的锥入深度h。去掉锥尖处的凡士林，取10g以上的土样两个，测定其含水率，计算含水率平均值ω。重复上述步骤，对其他两个含水率土样进行试验，测其锥入深度与含水率。在双对数坐标纸上，以含水率为横坐标，锥入深度为纵坐标，点绘a、b、c三点，连此3点应成一条直线，如3点不在同一条直线上，要过a点与b、c两点连两条直线，根据液限在hp-ω1图上查得hp，以此hp再在h-ω图的ab、ac两直线上求出相应的两个含水率，当两个含水率的差值小于2%时，以该两点含水率的平均值与a点连成一条直线，当两个含水率差值大于2%时，应重做试验。在h-ω图上，对于100g锥，查得纵坐标入土深度h=20mm所对应的含水率为液限；若采用76g锥，则在h-ω图上，查得纵坐标入土深度h=17mm所对应的横坐标的含水率为土样的液限。经测定，该黄土的液限为[X]%，塑限为[X]%。2.2冻融循环作用原理及自然环境中的冻融情况2.2.1冻融循环的物理过程冻融循环是一个复杂的物理过程，主要涉及水分的冻结膨胀和融化收缩，这一过程对黄土内部结构产生着显著影响。当温度下降时，黄土中的水分开始冻结。水在冻结过程中，其体积会膨胀约9%。这是因为水分子在低温下形成规则的晶体结构，分子间的距离增大，导致体积膨胀。在黄土孔隙中，冻结的水分形成冰晶，冰晶的生长会对周围的黄土颗粒产生向外的压力。随着冰晶的不断生长，这种压力逐渐增大，当超过黄土颗粒间的连接强度时，就会使颗粒间的连接被破坏，孔隙结构发生改变。冰晶的膨胀可能会使黄土颗粒发生位移，原本紧密排列的颗粒变得松散，大孔隙数量增加，孔隙尺寸增大。当温度回升时，黄土中的冰晶开始融化。融化后的水分会填充到孔隙中，使孔隙中的水分含量增加。由于冰晶融化后体积缩小，黄土内部会出现一定的空隙，导致土体结构变得松弛。水分在融化过程中还可能会携带部分细小颗粒发生迁移，进一步改变黄土的孔隙结构和颗粒分布。这种反复的冻融循环过程，使得黄土的物理结构不断受到破坏和重塑，进而对其力学性质产生重要影响。随着冻融循环次数的增加，黄土的孔隙率逐渐增大，颗粒间的摩擦力和黏聚力逐渐减小，抗剪强度降低，压缩性增大。2.2.2自然环境中黄土冻融循环的特点在自然环境中，黄土的冻融循环具有明显的地域和季节差异，其温度变化和循环频率等特点受多种因素影响。不同地区的气候条件不同，导致黄土冻融循环的温度变化范围存在差异。在高纬度地区或海拔较高的山区，冬季气温较低，黄土冻结深度较大，冻结温度可达到零下十几摄氏度甚至更低。而在中低纬度地区，冬季气温相对较高，黄土冻结深度较浅，冻结温度一般在零下几摄氏度。在我国东北地区，冬季平均气温可达-20°C左右，黄土冻结深度可达1-2m；而在华北地区，冬季平均气温一般在-5°C至-10°C之间，黄土冻结深度通常在0.5-1m。季节变化也使得黄土冻融循环呈现出不同的特点。在冬季，随着气温逐渐降低，黄土开始进入冻结阶段，冻结过程较为缓慢。而在春季，气温回升较快，黄土中的冰迅速融化，融化过程相对较短。这种季节性的冻融循环频率也有所不同，在寒冷地区，一年中可能经历多次冻融循环；而在温暖地区，冻融循环次数相对较少。在青藏高原地区，由于气候寒冷，一年中可能会经历5-10次冻融循环；而在黄土高原部分地区，一年中可能仅经历1-3次冻融循环。此外，自然环境中的降水、日照等因素也会对黄土冻融循环产生影响。降水会增加黄土中的含水率，使得在冻结过程中水分冻结膨胀的作用更加明显，对黄土结构的破坏更大。日照时间和强度会影响黄土的升温速度和降温速度，进而影响冻融循环的进程。在日照时间长、强度大的地区，黄土在白天升温较快，夜间降温也较快，冻融循环的频率可能会相对增加。三、冻融循环作用对黄土力学性质影响的试验研究3.1试验材料与方案设计3.1.1试验材料的选取与制备本试验所用黄土样品采自[具体地点]，该地区属于典型的黄土分布区域，其黄土具有代表性。采样点位于[详细的采样地点描述，如某山坡、某场地等]，此处黄土堆积深厚，且未受到明显的人为扰动和其他地质因素的强烈影响。在采集过程中，采用了专门的原状土采样器，以确保采集的黄土样品保持其天然结构和状态。采样器为[采样器型号]，其内径为[X]mm，高度为[X]mm。首先，在采样点处清理表面杂物，然后将采样器垂直压入黄土中，缓慢旋转并下压，直至采样器完全进入黄土预定深度。取出采样器后，用密封胶带将两端密封，以防止水分散失和结构扰动，并及时运回实验室进行后续处理。在实验室中，对采集的黄土样品进行了初步处理。将样品从采样器中小心取出，去除表面的浮土和杂质。对于用于物理性质测试的样品，采用四分法将其分成若干份，每份样品质量约为[X]g，用于测定初始干密度、含水率、液塑限等物理参数。对于用于冻融循环试验和力学性质测试的样品，根据试验要求制备成相应的尺寸和形状。制备无侧限抗压强度试验样品时，使用环刀（内径[X]mm，高度[X]mm）在原状土样上切取，确保样品的完整性和均匀性，每个环刀样品质量约为[X]g；制备三轴剪切试验样品时，使用专门的制样模具将黄土制成直径为[X]mm，高度为[X]mm的圆柱体，每个样品质量约为[X]g；制备直剪试验样品时，将黄土制成边长为[X]mm的正方体，每个样品质量约为[X]g。制备好的样品均用保鲜膜包裹，放置在保湿缸中备用，以保持其含水率不变。3.1.2冻融循环试验方案设计冻融循环试验在高精度的冻融循环试验箱中进行，该试验箱型号为[试验箱型号]，能够精确控制温度和时间。试验温度控制范围设定为-20°C至20°C，这一温度范围能够较好地模拟自然环境中黄土可能经历的温度变化。在寒冷的冬季，黄土的冻结温度可能会达到-20°C甚至更低，而在春季气温回升时，黄土的融化温度最高可接近20°C。循环次数设置为0次、5次、10次、15次、20次，以研究不同冻融循环次数对黄土力学性质的影响。选择这几个循环次数是基于对实际工程和自然环境的考虑。在实际工程中，黄土可能会在几年到几十年的时间内经历不同次数的冻融循环，通过设置这几个具有代表性的循环次数，可以初步了解冻融循环对黄土力学性质的长期影响趋势。在自然环境中，不同地区的黄土冻融循环次数也有所不同，通过这些循环次数的设置，可以涵盖一定范围内的自然情况。具体的试验步骤如下：将制备好的黄土样品放入冻融循环试验箱中，首先将温度降至-20°C，保持[X]小时，模拟黄土的冻结过程。在冻结过程中，黄土中的水分逐渐结冰，体积膨胀，对黄土结构产生压力。然后将温度升至20°C，保持[X]小时，模拟黄土的融化过程。在融化过程中，冰逐渐融化成水，黄土结构发生变化。如此完成一次冻融循环，按照设定的循环次数重复上述过程。在每次冻融循环过程中，使用高精度温度传感器实时监测样品内部的温度变化，确保温度控制的准确性。3.1.3力学性质测试试验方案在完成冻融循环试验后，对黄土样品进行无侧限抗压强度、三轴剪切、直剪等力学性质测试。无侧限抗压强度测试采用万能材料试验机，型号为[试验机型号]，加载速率设定为[X]mm/min。将经过冻融循环的黄土样品放置在试验机的承压板上，确保样品与承压板紧密接触且中心对齐。启动试验机，以恒定的加载速率对样品施加垂直压力，直至样品破坏，记录破坏时的荷载值。根据样品的尺寸和破坏荷载，计算无侧限抗压强度，计算公式为：q_{u}=\frac{P}{A}，其中q_{u}为无侧限抗压强度，P为破坏荷载，A为样品的横截面积。三轴剪切试验采用三轴剪切仪，型号为[仪器型号]。试验前，将黄土样品用乳胶膜包裹，放入三轴压力室中。首先对样品施加围压\sigma_{3}，围压值分别设定为50kPa、100kPa、150kPa，以模拟不同的实际受力情况。在实际工程中，黄土可能会受到不同大小的侧向压力，通过设置这几个围压值，可以研究围压对黄土力学性质的影响。然后以[X]mm/min的速率施加轴向压力，使样品逐渐发生剪切破坏，记录破坏时的轴向应力\sigma_{1}和对应的轴向应变\varepsilon_{1}。根据试验数据，计算得到黄土的抗剪强度参数，黏聚力c和内摩擦角\varphi，计算公式为：c=\frac{\sigma_{1f}-\sigma_{3f}}{2}\sin\varphi，\varphi=\arcsin\frac{\sigma_{1f}-\sigma_{3f}}{\sigma_{1f}+\sigma_{3f}}，其中\sigma_{1f}和\sigma_{3f}分别为破坏时的大主应力和小主应力。直剪试验使用直剪仪，型号为[直剪仪型号]。将经过冻融循环的黄土样品放入直剪盒中，对样品施加竖向压力\sigma，竖向压力值分别为100kPa、200kPa、300kPa、400kPa。同样，这些竖向压力值是根据实际工程中黄土可能承受的压力范围来确定的。以[X]mm/min的剪切速率对样品进行剪切，记录剪切过程中的剪应力\tau和剪切位移\Deltal，直至样品发生剪切破坏。根据试验数据，绘制剪应力-剪切位移曲线，确定样品的抗剪强度。通过不同竖向压力下的抗剪强度，利用库仑定律计算得到黏聚力c和内摩擦角\varphi，库仑定律表达式为：\tau=c+\sigma\tan\varphi。3.2试验结果与数据分析3.2.1冻融循环次数对黄土无侧限抗压强度的影响随着冻融循环次数的增加，黄土的无侧限抗压强度呈现出先增加后减小的变化趋势。在冻融循环次数为0次时，即原状黄土，其无侧限抗压强度为[X]kPa。当冻融循环次数增加到5次时，无侧限抗压强度增加至[X]kPa，较原状黄土提高了[X]%。这是因为在冻融循环初期，水分冻结形成的冰晶对黄土颗粒间的连接有一定的强化作用，使得颗粒间的接触更加紧密，从而提高了黄土的结构强度，进而增加了无侧限抗压强度。然而，当冻融循环次数继续增加到10次时，无侧限抗压强度开始下降，降至[X]kPa，较5次冻融循环时降低了[X]%。随着冻融循环次数进一步增加到15次和20次，无侧限抗压强度分别降至[X]kPa和[X]kPa。这是由于随着冻融循环次数的增多，冰晶的反复膨胀和融化对黄土结构造成了严重破坏。冰晶膨胀时产生的压力使黄土颗粒间的连接逐渐被破坏，孔隙结构不断恶化，颗粒间的摩擦力和黏聚力逐渐减小，导致黄土的结构强度降低，最终使得无侧限抗压强度逐渐减小。3.2.2冻融循环作用下黄土抗剪强度参数的变化经过冻融循环后，黄土的黏聚力和内摩擦角均发生了明显变化，这些变化对土体稳定性有着重要影响。在三轴剪切试验中，原状黄土的黏聚力c为[X]kPa，内摩擦角\varphi为[X]°。当冻融循环次数达到5次时，黏聚力c增加至[X]kPa，内摩擦角\varphi变化不大，略有增加至[X]°。在冻融循环初期，冰晶对黄土颗粒间的胶结作用使得颗粒间的连接增强，从而增加了黏聚力；而内摩擦角变化不大是因为此时黄土颗粒的形状和排列方式尚未发生显著改变，颗粒间的摩擦特性基本保持稳定。随着冻融循环次数增加到10次，黏聚力c开始下降至[X]kPa，内摩擦角\varphi也有所下降至[X]°。当冻融循环次数达到15次和20次时，黏聚力c分别降至[X]kPa和[X]kPa，内摩擦角\varphi分别降至[X]°和[X]°。随着冻融循环次数的不断增加，黄土结构逐渐被破坏，颗粒间的胶结物质被破坏，颗粒间的连接减弱，导致黏聚力大幅下降；同时，颗粒的排列变得更加松散，颗粒间的咬合作用减弱，使得内摩擦角也逐渐减小。黏聚力和内摩擦角是影响土体抗剪强度的重要参数，它们的减小会导致土体抗剪强度降低，从而降低土体的稳定性。在工程实际中，如黄土边坡，土体抗剪强度的降低可能会导致边坡更容易发生滑动破坏，威胁工程安全。因此，在黄土地区的工程设计和施工中，必须充分考虑冻融循环对黄土抗剪强度参数的影响，采取相应的措施来提高土体的稳定性。3.2.3不同含水率黄土在冻融循环后的力学性质差异为了研究含水率对黄土在冻融循环后力学性质的影响，选取了初始含水率分别为10%、15%、20%、25%的黄土样品进行试验。经过相同次数（10次）的冻融循环后，不同含水率黄土的力学性质表现出明显差异。在无侧限抗压强度方面，初始含水率为10%的黄土样品，无侧限抗压强度为[X]kPa；初始含水率为15%的黄土样品，无侧限抗压强度为[X]kPa；初始含水率为20%的黄土样品，无侧限抗压强度为[X]kPa；初始含水率为25%的黄土样品，无侧限抗压强度为[X]kPa。可以看出，随着初始含水率的增加，黄土的无侧限抗压强度逐渐降低。这是因为含水率较高时，在冻融循环过程中，水分冻结膨胀产生的冰晶压力更大，对黄土结构的破坏更为严重，导致无侧限抗压强度降低。在三轴剪切试验中，初始含水率为10%的黄土样品，黏聚力c为[X]kPa，内摩擦角\varphi为[X]°；初始含水率为15%的黄土样品，黏聚力c为[X]kPa，内摩擦角\varphi为[X]°；初始含水率为20%的黄土样品，黏聚力c为[X]kPa，内摩擦角\varphi为[X]°；初始含水率为25%的黄土样品，黏聚力c为[X]kPa，内摩擦角\varphi为[X]°。随着初始含水率的增加，黏聚力和内摩擦角均呈现下降趋势。含水率的增加使得黄土颗粒间的润滑作用增强，颗粒间的连接减弱，从而导致黏聚力和内摩擦角减小。这些结果表明，含水率是影响黄土在冻融循环后力学性质的重要因素。在工程实践中，对于黄土地区的工程，如地基处理、边坡防护等，应充分考虑黄土的含水率情况，采取相应的措施来控制含水率，以提高工程的稳定性和安全性。3.2.4冻结温度对黄土力学性质的影响分析为了研究不同冻结温度对黄土力学性质的影响，设置了冻结温度分别为-10°C、-15°C、-20°C的试验工况，对黄土样品进行冻融循环试验（循环次数均为10次）后，测试其力学性质。在无侧限抗压强度方面，当冻结温度为-10°C时，黄土的无侧限抗压强度为[X]kPa；当冻结温度为-15°C时，无侧限抗压强度为[X]kPa；当冻结温度为-20°C时，无侧限抗压强度为[X]kPa。随着冻结温度的降低，无侧限抗压强度逐渐减小。这是因为冻结温度越低，水分冻结时产生的冰晶体积膨胀越大，对黄土结构的破坏越严重，从而导致无侧限抗压强度降低。在三轴剪切试验中，当冻结温度为-10°C时，黄土的黏聚力c为[X]kPa，内摩擦角\varphi为[X]°；当冻结温度为-15°C时，黏聚力c为[X]kPa，内摩擦角\varphi为[X]°；当冻结温度为-20°C时，黏聚力c为[X]kPa，内摩擦角\varphi为[X]°。随着冻结温度的降低，黏聚力和内摩擦角均呈现下降趋势。较低的冻结温度使得冰晶对黄土颗粒间的连接破坏更严重，颗粒间的摩擦力和黏聚力减小，导致黏聚力和内摩擦角降低。在变形模量方面，当冻结温度为-10°C时，黄土的变形模量为[X]MPa；当冻结温度为-15°C时，变形模量为[X]MPa；当冻结温度为-20°C时，变形模量为[X]MPa。随着冻结温度的降低，变形模量逐渐减小，说明黄土在较低冻结温度下，抵抗变形的能力减弱。综上所述，冻结温度对黄土的力学性质有着显著影响。在工程实际中，尤其是在寒冷地区的黄土工程中，应充分考虑冻结温度的因素，合理设计工程结构，采取有效的保温措施，以减小冻结温度对黄土力学性质的不利影响，确保工程的安全稳定。四、冻融循环作用下黄土微观结构变化与力学性质损伤机制4.1冻融循环对黄土微观结构的影响4.1.1微观结构观测方法与结果为了深入探究冻融循环对黄土微观结构的影响，本研究采用了扫描电镜（SEM）和压汞仪（MIP）等先进的微观测试手段。扫描电镜能够直观地呈现黄土微观结构的形态特征，而压汞仪则可以精确地测定黄土孔隙结构的参数。在扫描电镜观测中，将经过不同冻融循环次数的黄土样品进行精心制备。首先，对样品进行干燥处理，以去除水分对微观结构观测的干扰。然后，将干燥后的样品进行喷金处理，增强其导电性，以便在扫描电镜下获得清晰的图像。通过扫描电镜观察发现，原状黄土的颗粒排列相对紧密，颗粒间的接触较为紧密，存在一定的胶结物质将颗粒连接在一起，孔隙大小相对均匀，多为小孔径孔隙，孔隙形状不规则，连通性较差。当冻融循环次数达到5次时，黄土颗粒间的连接开始出现松动，部分胶结物质被破坏，孔隙数量略有增加，孔径有所增大，大孔隙的数量开始增多，颗粒的排列也变得相对松散。随着冻融循环次数增加到10次，黄土微观结构的变化更为明显，颗粒间的连接进一步削弱，胶结物质大量减少，孔隙结构变得更加复杂，大孔隙和小孔隙相互连通，形成了较为复杂的孔隙网络。当冻融循环次数达到15次和20次时，黄土的微观结构遭到严重破坏，颗粒间的连接几乎完全丧失，大量大孔隙出现，颗粒呈现出离散状态，孔隙连通性良好，形成了贯通的孔隙通道。相关的微观结构图像如图1所示。[此处插入不同冻融循环次数下黄土微观结构的SEM图像，如原状黄土、冻融5次、10次、15次、20次的图像]利用压汞仪对黄土孔隙结构进行测试，能够得到孔隙大小分布、孔隙体积、孔隙率等重要参数。测试结果表明，原状黄土的孔隙率为[X]%，孔隙体积主要集中在小孔径范围内，孔径分布相对集中。经过5次冻融循环后，黄土的孔隙率增加至[X]%，大孔径孔隙的体积分数有所增加，小孔径孔隙的体积分数略有减小，孔径分布范围变宽。当冻融循环次数达到10次时，孔隙率进一步增加到[X]%，大孔径孔隙的体积分数显著增加，小孔径孔隙的体积分数明显减小，孔隙大小分布更加分散。随着冻融循环次数继续增加到15次和20次，孔隙率分别增加到[X]%和[X]%，大孔径孔隙的体积分数持续增加，小孔径孔隙的体积分数进一步减小，孔隙分布呈现出明显的双峰特征，即大孔径孔隙和小孔径孔隙并存，且大孔径孔隙的占比逐渐增大。相关的孔隙结构参数变化曲线如图2所示。[此处插入不同冻融循环次数下黄土孔隙率、孔隙体积、孔径分布等参数的变化曲线]4.1.2微观结构变化与宏观力学性质的关联黄土微观结构的变化与宏观力学性质之间存在着紧密的内在联系。随着冻融循环次数的增加，黄土微观结构的破坏逐渐加剧，这直接导致了宏观力学性质的劣化。在抗剪强度方面，原状黄土由于颗粒排列紧密，颗粒间的胶结物质和摩擦力能够有效地抵抗剪切力，因此具有较高的抗剪强度。然而，随着冻融循环次数的增加，颗粒间的连接逐渐被破坏，胶结物质减少，颗粒排列变得松散，使得抗剪强度逐渐降低。在三轴剪切试验中，原状黄土的黏聚力为[X]kPa，内摩擦角为[X]°；经过10次冻融循环后，黏聚力降至[X]kPa，内摩擦角降至[X]°。这是因为微观结构的变化使得颗粒间的咬合作用减弱，摩擦力减小，同时胶结物质的减少也降低了颗粒间的黏聚力，从而导致抗剪强度下降。在变形特性方面，微观结构的变化也对黄土的压缩性和变形模量产生了显著影响。原状黄土的孔隙率较低，颗粒间的接触紧密，在受到压力时，颗粒能够有效地传递应力，变形较小，变形模量较大。随着冻融循环次数的增加，孔隙率增大，孔隙结构变得松散，颗粒间的接触点减少，在受到压力时，颗粒容易发生相对位移，导致变形增大，变形模量减小。在压缩试验中，原状黄土的压缩系数为[X]MPa⁻¹，经过10次冻融循环后，压缩系数增大至[X]MPa⁻¹，表明黄土的压缩性增强，抵抗变形的能力减弱。综上所述，黄土微观结构的变化是导致其宏观力学性质损伤的重要原因。通过对微观结构变化与宏观力学性质之间关联的研究，可以更深入地理解冻融循环作用下黄土力学性质变化的内在机制，为黄土地区的工程建设提供更科学的理论依据。4.2黄土力学性质损伤机制分析4.2.1冻融过程中水分迁移与应力变化在冻融循环过程中，水分迁移是导致黄土力学性质损伤的关键因素之一，它与土体内部应力状态的改变密切相关。当温度下降时，黄土中的水分开始冻结。水分冻结过程中，冰晶的生长会产生体积膨胀，这种膨胀力对周围的土体颗粒产生挤压作用。由于黄土中水分分布不均匀，冰晶生长所产生的膨胀力也不均匀，这就导致土体内部产生了不均匀的应力分布。在孔隙较小的区域，冰晶生长受到限制，产生的膨胀力较大；而在孔隙较大的区域，冰晶有更多的生长空间，膨胀力相对较小。这种不均匀的应力分布会使土体颗粒发生位移和变形，从而破坏颗粒间的连接，导致土体结构的损伤。随着温度的回升，黄土中的冰晶开始融化。融化后的水分在重力和毛细管力的作用下发生迁移。水分迁移会改变土体内部的水分分布，使得原本干燥的区域变得湿润，而湿润区域的含水率进一步增加。这种水分分布的变化会导致土体的物理性质发生改变，如密度、孔隙比等。水分迁移还会对土体内部的应力状态产生影响。当水分从高处向低处迁移时，会对土体颗粒产生拖曳力，使得颗粒间的摩擦力和黏聚力发生变化。在水分迁移过程中，土体内部还可能会产生渗流力，进一步改变土体的应力状态。如果渗流力过大，可能会导致土体颗粒的流失，从而破坏土体的结构，降低土体的力学性能。此外，冻融循环过程中，土体的温度变化也会导致热胀冷缩现象的发生。土体在温度变化时，不同部位的膨胀和收缩程度不同，这也会在土体内部产生应力。这种热应力与冰晶膨胀力和水分迁移所产生的应力相互叠加，使得土体内部的应力状态更加复杂。长期的冻融循环作用下，土体内部的应力不断积累，超过土体的承受能力时，就会导致土体结构的破坏，力学性质的劣化。4.2.2颗粒间连接破坏与结构重组冻融循环作用会对黄土颗粒间的连接产生严重破坏，进而引发土体结构的重组，这是导致黄土力学性能劣化的重要机制。在天然状态下，黄土颗粒间存在着多种连接方式，包括化学胶结、静电引力、摩擦力等。这些连接方式使得黄土颗粒能够相互作用，形成相对稳定的结构。化学胶结是通过黄土中的胶结物质，如碳酸钙、铁铝氧化物等，将颗粒粘结在一起，提供了较强的连接强度。静电引力则是由于颗粒表面电荷的存在，使得颗粒之间产生相互吸引的作用力。摩擦力则是颗粒之间相互接触和摩擦所产生的阻力。在冻融循环过程中，水分的冻结和融化对颗粒间的连接产生了强烈的破坏作用。当水分冻结时，冰晶的膨胀会对颗粒间的连接产生挤压和拉伸作用。冰晶的膨胀力可能会超过颗粒间连接的强度，导致化学胶结物质的破裂、静电引力的减弱和摩擦力的降低。随着冻融循环次数的增加，这种破坏作用不断累积，使得颗粒间的连接逐渐丧失。在融化过程中，水分的迁移会带走部分胶结物质，进一步削弱颗粒间的连接。随着颗粒间连接的破坏，黄土的土体结构开始发生重组。原本紧密排列的颗粒变得松散，颗粒之间的相对位置发生改变。在冰晶膨胀和水分迁移的作用下，颗粒可能会发生位移、旋转和重新排列。一些大颗粒可能会被推挤到孔隙中，导致孔隙结构的改变。原本连通性较差的孔隙可能会因为颗粒的移动而相互连通，形成更大的孔隙通道。这种结构重组使得黄土的微观结构变得更加复杂和不稳定，宏观上表现为土体的强度降低、压缩性增大。土体结构的重组还会影响黄土的其他力学性能。例如，孔隙结构的改变会影响土体的渗透性，使得水分在土体中的运移更加容易，进一步加剧了土体的劣化。颗粒排列的变化还会影响土体的各向异性，使得土体在不同方向上的力学性能出现差异。因此，颗粒间连接破坏与结构重组是冻融循环作用下黄土力学性质损伤的重要原因，深入研究这一机制对于理解黄土的力学行为具有重要意义。4.2.3基于损伤力学的黄土力学性质损伤模型构建为了更准确地描述冻融循环作用下黄土力学性质的损伤过程，本研究尝试引入损伤力学理论，建立考虑冻融循环影响的黄土力学性质损伤模型。损伤力学是一门研究材料在外部荷载和环境因素作用下内部损伤演化规律的学科，通过引入损伤变量来描述材料内部结构的劣化程度。在构建黄土损伤模型时，首先确定损伤变量的定义和计算方法。考虑到冻融循环对黄土微观结构和力学性质的影响，选择以孔隙率的变化作为损伤变量的度量。随着冻融循环次数的增加，黄土的孔隙率逐渐增大，孔隙结构不断恶化，这与损伤的概念相符合。损伤变量D的计算公式为：D=\frac{n-n_0}{n_{max}-n_0}，其中n为经过冻融循环后的孔隙率，n_0为原状黄土的孔隙率，n_{max}为孔隙率的最大值，可通过试验或理论分析确定。根据损伤力学理论，建立黄土的应力-应变关系。在弹性阶段，考虑损伤的影响，应力-应变关系可表示为：\sigma=(1-D)E\varepsilon，其中\sigma为应力，\varepsilon为应变，E为原状黄土的弹性模量。在塑性阶段，引入损伤对屈服函数和硬化规律的影响。假设屈服函数为f(\sigma,\varepsilon^p,D)，其中\varepsilon^p为塑性应变。随着损伤的发展，屈服函数会发生变化，屈服面逐渐扩大，反映了土体强度的降低。硬化规律也会受到损伤的影响，塑性应变的发展会加快，导致土体的变形能力增强。为了确定模型中的参数，进行了一系列的室内试验。通过冻融循环试验和力学性质测试，获得不同冻融循环次数下黄土的孔隙率、弹性模量、屈服应力等参数。利用这些试验数据，采用最小二乘法等参数拟合方法，对损伤模型中的参数进行优化，使得模型能够较好地拟合试验结果。经过参数拟合和验证，建立的损伤模型能够较好地描述冻融循环作用下黄土力学性质的损伤过程，为黄土地区的工程分析和设计提供了有力的工具。五、考虑冻融循环的黄土力学性质数值模拟5.1数值模拟方法与模型建立5.1.1数值模拟软件介绍与选择本研究选用有限元软件ABAQUS进行数值模拟。ABAQUS是一款功能强大的通用有限元分析软件，在岩土工程模拟领域具有显著优势。其具备丰富的材料本构模型库，能够提供多种适用于岩土材料的本构模型，如摩尔库仑模型、修正剑桥模型、Drucker-Prager模型等。这些模型能够真实地反映土体在复杂应力状态下的性状，包括土体的剪胀性、屈服性、硬化特性等，为准确模拟黄土在冻融循环作用下的力学行为提供了有力工具。ABAQUS拥有高效稳定的求解器，能够处理复杂的非线性问题和大规模计算。在岩土工程中，土体的力学行为往往呈现出高度的非线性，如材料非线性、几何非线性和接触非线性等。ABAQUS的求解器可以有效地处理这些非线性因素，准确模拟土体在冻融循环过程中的大变形、复杂的接触相互作用以及材料性能的变化，为工程师提供精确的仿真结果。该软件还提供了方便的建模工具和强大的网格划分功能。工程师可以轻松地建立复杂的岩土工程模型，包括隧道、边坡、地基等结构，以及考虑地质条件、施工过程和外部荷载等因素的复杂模型。在模拟黄土冻融循环时，可以方便地对黄土模型进行建模，并根据需要对不同区域进行网格划分，以提高计算精度和效率。ABAQUS的后处理功能强大且直观，工程师可以方便地查看和分析计算结果，包括位移、应力、应变等关键参数的图形化显示，以及敏感性分析、优化设计等高级功能，有助于深入理解黄土在冻融循环作用下的力学响应。5.1.2黄土本构模型的选择与参数确定根据黄土的特性以及冻融循环作用的影响，本研究选择修正剑桥模型作为黄土的本构模型。修正剑桥模型是在剑桥模型的基础上发展而来，它考虑了土体的剪胀性和硬化特性，能够较好地描述正常固结土和弱超固结土在复杂应力路径下的力学行为，对于黄土这种具有一定结构性和特殊力学性质的土体具有较好的适用性。修正剑桥模型的参数主要包括：弹性模量E、泊松比\nu、临界状态线斜率M、压缩指数\lambda、回弹指数\kappa、初始孔隙比e_0等。这些参数的确定对于准确模拟黄土的力学行为至关重要。本研究通过室内试验和相关文献资料来确定模型参数。通过三轴剪切试验获得黄土的应力-应变关系曲线，根据曲线特征和相关理论计算得到弹性模量E和泊松比\nu。通过固结试验确定压缩指数\lambda和回弹指数\kappa。初始孔隙比e_0则根据黄土的初始物理参数测定结果确定。临界状态线斜率M可通过三轴排水剪切试验数据，利用相关公式计算得出。对于冻融循环作用下黄土的参数变化，参考试验研究结果，考虑冻融循环对黄土孔隙结构、颗粒间连接等因素的影响，对参数进行适当调整。在冻融循环次数增加时，黄土的孔隙率增大，颗粒间连接减弱，弹性模量E和剪切模量G会相应减小，压缩指数\lambda会增大，通过试验数据拟合得到参数随冻融循环次数变化的函数关系，从而在数值模拟中能够更准确地反映冻融循环对黄土力学性质的影响。5.1.3冻融循环作用的数值模拟实现方法在数值模型中，通过设置热-力耦合分析步来实现冻融循环过程的模拟。首先，定义黄土的热力学参数，包括热膨胀系数\alpha_T、导热系数k、比热容c_p等。热膨胀系数反映了黄土在温度变化时的体积膨胀或收缩特性，导热系数决定了热量在黄土中的传导速率，比热容则表示单位质量黄土温度升高1°C所需吸收的热量。这些参数可通过相关试验测定或参考已有文献资料确定。在分析步中，设置温度边界条件来模拟冻融循环过程中的温度变化。施加周期性的温度载荷，如将温度从-20°C变化到20°C，模拟一个完整的冻融循环周期。通过控制温度变化的时间历程和循环次数，实现不同冻融循环条件的模拟。在每个时间步中，ABAQUS会根据热传导方程计算黄土内部的温度分布，考虑到水分迁移对热量传输的影响，在模型中引入水分迁移方程，与热传导方程进行耦合求解。水分迁移主要受到温度梯度、孔隙水压力梯度等因素的影响，通过建立相应的数学模型来描述水分在黄土孔隙中的迁移过程。在模拟过程中，还需要考虑冻融过程中黄土的相变问题，即水分的冻结和融化。通过引入相变潜热的概念，在温度变化过程中，当温度达到水的冻结点或融化点时，考虑相变潜热对热量传输的影响，以更准确地模拟冻融循环过程中黄土的热力学行为。通过以上方法，能够在数值模型中较为真实地实现冻融循环作用的模拟，为研究黄土在冻融循环下的力学性质提供基础。5.2数值模拟结果与试验结果对比验证5.2.1数值模拟结果展示通过ABAQUS数值模拟，得到了黄土在冻融循环作用下的应力、应变分布云图以及力学性质随冻融循环次数的变化曲线。在应力分布云图中，可以清晰地看到在冻融循环过程中，黄土内部应力的分布情况。在冻结阶段，由于水分结冰膨胀，黄土内部产生较大的应力，尤其是在孔隙周围和颗粒接触部位，应力集中现象较为明显。随着冻融循环次数的增加，应力集中区域逐渐扩大，且应力值也有所变化。在融化阶段，应力分布相对较为均匀，但由于土体结构的变化，整体应力水平也会发生改变。应变分布云图则直观地展示了黄土在冻融循环过程中的变形情况。在冻结过程中，黄土会发生膨胀变形，应变较大的区域主要集中在孔隙周围和土体的薄弱部位。随着冻融循环次数的增加，土体的变形逐渐累积，应变分布范围扩大，且变形程度也逐渐增大。在融化阶段，土体发生收缩变形，应变方向与冻结阶段相反。力学性质随冻融循环次数的变化曲线表明，随着冻融循环次数的增加，黄土的无侧限抗压强度逐渐降低。在冻融循环初期，无侧限抗压强度下降较为缓慢，但随着循环次数的进一步增加，下降速度加快。抗剪强度参数黏聚力和内摩擦角也呈现出类似的变化趋势，黏聚力和内摩擦角均随冻融循环次数的增加而逐渐减小。这些变化趋势与试验研究结果具有一定的相似性，初步验证了数值模拟的合理性。相关的应力、应变分布云图和力学性质变化曲线如图3-图5所示。[此处插入黄土在冻融循环作用下的应力分布云图、应变分布云图以及力学性质随冻融循环次数的变化曲线]5.2.2模拟结果与试验结果对比分析将数值模拟结果与试验结果进行对比分析，以评估数值模型的准确性和可靠性。在无侧限抗压强度方面，数值模拟得到的无侧限抗压强度随冻融循环次数的变化趋势与试验结果基本一致，均呈现出先略微增加后逐渐减小的趋势。在冻融循环初期，由于冰晶的胶结作用，数值模拟和试验结果中的无侧限抗压强度都有一定程度的增加。但随着冻融循环次数的增加，由于土体结构的破坏，两者的无侧限抗压强度都逐渐降低。然而，数值模拟结果与试验结果在具体数值上存在一定差异。在冻融循环次数为10次时，试验测得的无侧限抗压强度为[X]kPa，而数值模拟结果为[X]kPa，相对误差为[X]%。这种差异可能是由于数值模型中对黄土的本构关系简化、参数取值的不确定性以及试验过程中的误差等因素导致的。在抗剪强度参数方面，数值模拟得到的黏聚力和内摩擦角随冻融循环次数的变化趋势也与试验结果相符，均逐渐减小。但同样在具体数值上存在一定偏差。在冻融