体变视角下合肥地区非饱和膨胀土土水特征的多维度解析

体变视角下合肥地区非饱和膨胀土土水特征的多维度解析一、绪论1.1研究背景与意义膨胀土作为一种特殊的黏性土，在我国分布广泛，合肥地区便是典型的膨胀土分布区域。其显著特点是在环境干湿交替的作用下，体积会明显胀缩，性质极不稳定，素有工程“癌症”之称。在合肥地区的工程建设中，膨胀土带来了诸多棘手问题。例如在建筑地基工程中，膨胀土的胀缩特性导致地基不均匀沉降，使建筑物墙体开裂、基础变形，严重影响建筑物的稳定性和安全性。在道路工程方面，膨胀土路基在水分变化时的体积变化，会造成路面起伏、开裂，降低道路的平整度和使用寿命，增加维护成本。在基坑工程中，合肥膨胀土基坑在开挖卸荷、干湿交替等外界因素作用下，土体强度急剧降低，基坑边坡稳定性急剧下降，由于土体内大量裂隙的存在，边坡破坏由表层浅部土块的滑裂变形开始逐渐向深部发展，如合肥某基坑工程就因膨胀土问题发生了边坡坍塌事故，造成了严重的经济损失和工期延误。非饱和膨胀土的土水特征对于深入理解其工程性质起着关键作用。土水特征直接影响着膨胀土的强度特性，基质吸力的变化会改变土体的有效应力，进而影响土体的抗剪强度。当土体含水率增加，基质吸力减小，抗剪强度降低，容易引发工程事故。土水特征与膨胀土的变形特性紧密相关，含水率的变化导致土体膨胀或收缩，而土水特征曲线能反映这种变化关系，为工程设计提供重要依据。在渗透特性方面，土水特征决定了水分在土体中的迁移规律，影响着土体的渗透性，对于工程中的防水、排水设计至关重要。研究合肥地区非饱和膨胀土在体变作用下的土水特征，能够为该地区的工程建设提供准确的参数和理论支持，有效提高工程的稳定性和可靠性，减少工程事故的发生，降低工程成本，具有重要的现实意义。1.2国内外研究现状1.2.1膨胀土体变特性研究进展国外对膨胀土体变特性的研究开展较早。一些学者通过大量试验，研究了不同地区膨胀土在干湿循环作用下的体变规律。发现随着干湿循环次数的增加，膨胀土的累积塑性变形逐渐增大，且体变特性与土的初始含水率、干密度等因素密切相关。在荷载作用下，膨胀土的体变不仅受到自身胀缩特性影响，还与荷载大小、加载方式等有关，高荷载会抑制膨胀土的膨胀变形，但会加剧其收缩变形。国内学者也对膨胀土体变特性进行了深入研究。针对南水北调工程中的膨胀土，通过室内和现场试验，分析了干湿循环和不同上覆荷载条件下的体变特性，得出上覆荷载能有效约束膨胀土的膨胀变形，而干湿循环会导致土体结构破坏，使其体变特性更加复杂的结论。还有学者研究了不同压实度膨胀土的体变规律，发现压实度越高，膨胀土的膨胀量越小，收缩量也相对减小，土体的稳定性更好。1.2.2土水特征曲线研究综述土水特征曲线是描述土壤含水量与土体等效孔隙比值之间关系的重要曲线，它反映了土壤力学及水文学特征的变化规律，在非饱和土力学中具有举足轻重的地位。其定义为土壤水势（土壤水吸力）与土壤水分含量之间的关系曲线，通常土壤含水量以体积百分数表示，土壤吸力以大气压表示。测量方法主要有张力计法、压力膜仪法和离心机法等。张力计可直接测定土壤的水势，但压力范围有限，无法抵达植物的萎蔫系数；压力膜仪是目前默认的仪器设备，国际上一般都承认其测定数据，它通过高压气泵或高压氮气瓶向密封容器中充气加压，从土壤样品中渗出水分来测定含水量，从而得到土水特征曲线；离心机法则是利用离心力使土壤中的水分排出，根据离心力和水分排出量的关系确定土水特征曲线。国外在土水特征曲线研究方面起步早，取得了众多成果。建立了多种理论模型来描述土水特征曲线，如VanGenuchten模型，该模型基于土壤孔隙大小分布理论，能较好地拟合大多数土壤的土水特征曲线，但对于结构复杂的膨胀土，拟合效果存在一定局限性。国内学者也对土水特征曲线进行了大量研究，针对不同地区的非饱和土，通过试验测定土水特征曲线，并对现有模型进行改进和优化，以提高对本土土壤的适用性。1.2.3非饱和膨胀土土水特征模型研究目前，已存在多种非饱和膨胀土土水特征模型。MK模型将非饱和膨胀土看做一种由按层分布的三种孔隙组成的多孔介质，包括内部气体（僵固气体）、外部空隙（可压缩空气）和孔隙水，支持基于物理机制的耦合分析，可以模拟土壤的压缩和蠕变性能，能自然地适应非饱和状态下土壤的行为特性，但该模型涉及方程较多，参数确定较为复杂，在实际应用中受到一定限制。VanGenuchten模型在非饱和土土水特征曲线模拟中应用广泛，它通过数学函数形式较好地描述了土壤水分特征曲线的形状，但对于膨胀土这种具有特殊结构和性质的土壤，在描述其吸力与含水率关系的某些阶段时不够准确。这些模型各自具有优缺点和适用条件。简单模型计算方便，但对复杂膨胀土特性的描述能力有限；复杂模型能更全面地考虑膨胀土的特性，但参数获取困难，计算量大。在实际应用中，需要根据具体工程需求和膨胀土的特性，合理选择和改进土水特征模型。1.2.4试验过程数值模拟研究现状数值模拟在非饱和膨胀土土水特征试验中发挥着重要作用。常用的模拟软件如GeoStudioVADOSE_W，是一款强大的地学模拟软件，专为处理涉及多孔介质流的工程问题而设计，能够模拟水在土壤中的运动，支持各种土壤水分特征曲线的输入，准确计算出水在土壤中的运移过程，还提供多孔介质流的可视化模拟，帮助用户深入理解土壤水分运动机制。国外利用数值模拟软件对非饱和膨胀土的土水特征进行了大量研究，模拟不同工况下水分在膨胀土中的迁移、分布以及土体的变形等，与试验结果对比验证，取得了较好的效果，为工程设计和分析提供了重要参考。国内也广泛开展了相关数值模拟研究，针对具体工程中的膨胀土问题，结合室内外试验数据，运用数值模拟软件进行分析，优化工程方案，提高工程的安全性和可靠性。但在数值模拟过程中，还存在模型简化、参数确定不准确等问题，需要进一步研究和改进。1.3研究内容与方法本研究旨在深入探究合肥地区非饱和膨胀土在体变作用下的土水特征，具体研究内容如下：合肥膨胀土体变特性研究：通过室内试验，研究合肥膨胀土在不同干湿循环次数、不同荷载条件下的体变规律，分析体变与初始含水率、干密度等因素的关系。合肥膨胀土土水特性试验及曲线拟合：采用压力膜仪等设备，测定合肥非饱和膨胀土在不同吸力下的含水率，获取土水特征曲线，并运用合适的数学模型进行曲线拟合，确定拟合参数。吸力-孔隙比-含水率关系研究：分析吸力、孔隙比和含水率之间的内在联系，建立三者之间的定量关系模型，进一步揭示非饱和膨胀土的土水特征。数值模拟研究：利用GeoStudioVADOSE_W等数值模拟软件，对合肥膨胀土在不同工况下的土水特征进行数值模拟，将模拟结果与试验结果对比分析，验证模型的准确性和可靠性，为工程应用提供理论支持。本研究采用试验研究与数值模拟相结合的方法。在试验研究方面，严格按照相关标准和规范进行室内试验，确保试验数据的准确性和可靠性。在数值模拟方面，根据试验结果合理选择和确定模拟参数，对模拟过程进行细致的设置和分析，充分发挥数值模拟的优势，与试验研究相互补充，全面深入地研究合肥地区非饱和膨胀土在体变作用下的土水特征。二、合肥膨胀土的体变特性试验研究及曲线拟合2.1试验基础准备本试验所用土样取自合肥地区典型膨胀土分布区域，在深度2-3m处进行采样，以获取具有代表性的原状土样。通过一系列标准土工试验，对土样的基本土性指标进行了详细测定。结果表明，该膨胀土的液限为58.5%，塑限为26.8%，塑性指数高达31.7，属于高液限、高塑性黏土。其自由膨胀率为56.0%，判定为中等膨胀土。土样的天然含水率为22.5%，天然密度为1.90g/cm³，孔隙比为0.75。这些基本土性指标为后续的体变特性试验研究提供了重要的基础数据。在试样制备过程中，对于原状土样，小心地将其切削成直径为61.8mm、高度为20mm的标准环刀试样，确保试样的完整性和结构不受破坏。对于扰动土样，先将土样风干、碾碎，过2mm筛，去除较大颗粒。然后按照预定的含水率和干密度，采用静压法制备试样。具体操作是将配置好的土样分多次装入模具中，每次装入后用小型千斤顶施加一定压力，使其达到规定的干密度。在量测方法上，采用百分表测量试样的竖向变形，将百分表固定在稳固的支架上，其测量头与试样的上表面垂直接触，确保测量的准确性。为了测量试样的体积变化，采用排水法，将试样放入盛满水的量筒中，测量排出水的体积，从而得到试样的体积。在测量含水率时，使用精度为0.01g的电子天平，采用烘干法，将试样在105-110℃的烘箱中烘干至恒重，通过前后质量差计算含水率。2.2加湿试验研究2.2.1加湿试验原理与仪器加湿试验基于水分迁移原理，利用环境湿度差使水分逐渐进入膨胀土试样。试验在恒温恒湿试验箱中进行，该试验箱型号为TH-225，能够精确控制试验环境的温度和湿度。温度控制范围为10-60℃，精度可达±0.5℃；湿度控制范围为30%-98%RH，精度为±3%RH。操作步骤如下：首先将制备好的膨胀土试样放入试验箱内的试样架上，设置试验箱的湿度为目标加湿湿度，如80%RH，温度设定为25℃，以模拟合肥地区常见的温湿度环境。试验过程中，定期使用电子天平称量试样的质量，记录加湿过程中试样质量随时间的变化，从而计算出试样的含水率变化。同时，每隔一定时间，采用游标卡尺测量试样的尺寸，包括直径和高度，以监测试样在加湿过程中的体积变化情况。2.2.2加湿曲线拟合模型与结果选择常用的指数函数模型对加湿曲线进行拟合，其表达式为：w=w_0+A(1-e^{-Bt})，式中w为加湿后试样的含水率，w_0为试样的初始含水率，A和B为拟合参数，t为加湿时间。利用最小二乘法对试验数据进行拟合，得到拟合参数A和B的值。例如，对于某一试样，拟合得到A=0.12，B=0.05。通过拟合优度R^2来评估拟合效果，R^2越接近1，说明拟合效果越好。经计算，该试样的拟合优度R^2=0.96，表明指数函数模型能够较好地描述该膨胀土试样的加湿过程。为了进一步验证拟合模型的准确性，对比了其他模型，如线性模型w=C+Dt（其中C和D为拟合参数）。线性模型拟合该试样的拟合优度R^2=0.85，明显低于指数函数模型的拟合优度。这表明线性模型在描述膨胀土的加湿过程时存在一定局限性，无法准确反映含水率随时间的非线性变化趋势，而指数函数模型能更准确地体现膨胀土在加湿过程中的特性。2.2.3测量方法精度对比在测量试样体积变化时，对比了游标卡尺法与照相法。游标卡尺法是直接使用精度为0.02mm的游标卡尺测量试样的直径和高度，然后根据圆柱体体积公式V=\pir^2h计算试样体积。该方法操作简单，但测量过程中可能因人为操作因素，如测量位置不准确、卡尺与试样接触力度不均匀等，导致测量误差。而且对于形状不规则的试样，游标卡尺法测量的准确性会受到较大影响。照相法是利用高清数码相机对试样进行拍照，然后通过图像分析软件，如ImageJ，识别试样边界，测量其尺寸并计算体积。该方法的优点是可以非接触式测量，避免了对试样的扰动，能够更全面地获取试样的形状信息，对于形状不规则的试样也能准确测量。但照相法的精度受到相机分辨率和图像分析算法的限制，相机分辨率较低时，可能无法准确识别试样边界，导致测量误差。而且图像分析过程较为复杂，需要一定的专业知识和技能。通过对多个试样的体积测量对比试验，发现对于形状规则的膨胀土试样，游标卡尺法和照相法的测量结果较为接近，相对误差在3%以内。但对于因膨胀而产生明显变形、形状不规则的试样，游标卡尺法的测量误差可达到8%，而照相法的测量误差能控制在5%以内。总体来说，照相法在测量变形较大、形状不规则的膨胀土试样体积时，具有更高的精度和可靠性，但操作相对复杂；游标卡尺法操作简便，适用于形状规则试样的测量。2.3收缩试验研究2.3.1收缩试验原理与仪器收缩试验原理是基于土在自然风干条件下失水导致体积收缩的现象。试验采用的主要仪器为收缩仪，该收缩仪由多孔板、百分表、支架等组成。多孔板上小孔的面积为整板面积的50%以上，以保证水分能够顺利排出。环刀选用直径61.8mm、高度20mm的标准环刀。试验操作流程如下：首先将制备好的试样用环刀取出，轻轻推出环刀，放置在多孔板上。在试样上表面放置百分表，百分表的测量头与试样紧密接触，用于测量试样的竖向收缩变形。将整套装置放置在通风良好、温度恒定（25℃）的室内环境中进行自然风干。每隔1-4h测记百分表读数，同时使用电子天平称取整套装置和试样的质量，准确至0.1g，以计算试样的含水率。2d后，由于收缩速度减慢，每隔6-24h测记百分表读数并称质量，直至两次百分表读数基本不变，表明试样收缩已基本稳定。试验结束后，将试样在105-110℃的烘箱中烘干，称取干土质量，用于后续计算。2.3.2收缩曲线拟合模型与缩限求解选用双曲线模型拟合收缩曲线，其表达式为：\delta=\frac{t}{a+bt}，式中\delta为线缩率，t为时间，a和b为拟合参数。通过对试验数据进行拟合，确定参数a和b的值。例如，对于某一试样，拟合得到a=2.5，b=0.3。利用拟合得到的模型可以预测不同时间的线缩率。缩限的求解采用作图法，以线缩率为纵坐标，含水率为横坐标，绘制关系曲线。在收缩曲线中，第一阶段土体收缩与含水率减少成正比，第二阶段随着含水率的减少，土体收缩率愈来愈小，第三阶段含水率继续减小，但土体不再收缩或收缩甚微。延长第一阶段和第三阶段的直线段至相交，交点所对应的横坐标即为原状土的缩限。经测定，该试样的缩限含水率为16.5%。收缩特性与缩限密切相关，缩限是土体收缩状态转变的关键含水率，当含水率高于缩限时，土体收缩明显；当含水率低于缩限时，土体收缩基本停止，了解这一关系对于分析膨胀土的工程性质具有重要意义。2.3.3失水路径对收缩特性的影响研究了不同失水路径下膨胀土的收缩特性，设置了快速失水和缓慢失水两种工况。快速失水工况是将试样放置在通风良好且温度较高（35℃）的环境中，加快水分蒸发；缓慢失水工况是将试样放置在温度为25℃、相对湿度为50%的恒温恒湿环境中。试验结果表明，快速失水时，膨胀土的收缩速率明显加快，在较短时间内达到较大的收缩量。这是因为快速失水导致土体内部水分迅速散失，颗粒间的吸力快速增大，从而使土体快速收缩。而缓慢失水时，收缩速率较为平缓，收缩量相对较小。失水速率还会影响土体的结构，快速失水容易导致土体内部产生较大的应力集中，使土体出现更多的微裂隙，破坏土体结构，进而影响其工程性质；缓慢失水时，土体结构破坏相对较小。环境湿度也是影响收缩特性的重要因素。在高湿度环境下，水分散失困难，收缩量较小；在低湿度环境下，水分容易散失，收缩量较大。当环境相对湿度从30%增加到70%时，相同时间内试样的收缩量降低了约30%。这说明在实际工程中，环境湿度的变化对膨胀土的收缩变形有显著影响，需要充分考虑环境湿度因素来评估膨胀土的工程稳定性。2.4缩胀试验研究2.4.1缩胀特性研究方法采用控制变量法研究膨胀土的缩胀特性。确定研究的变量包括初始干密度、试样形状、干湿循环次数等，控制因素为试验环境的温度和湿度，保持温度为25℃，相对湿度为60%。对于初始干密度的研究，制备不同初始干密度的试样，如1.60g/cm³、1.70g/cm³、1.80g/cm³，在相同的干湿循环条件下，观察和测量试样的缩胀变形。对于试样形状的研究，制作圆柱体和长方体两种形状的试样，圆柱体试样直径为61.8mm、高度为20mm，长方体试样尺寸为50mm×50mm×20mm，在相同工况下对比它们的缩胀特性。在研究干湿循环次数的影响时，设置干湿循环次数分别为1次、3次、5次，每次干湿循环过程中，先将试样饱和，然后在自然条件下风干至初始含水率，测量每次循环后的缩胀量。2.4.2初始干密度对缩胀特性的影响分析不同初始干密度的膨胀土在缩胀过程中的特性发现，初始干密度越大，膨胀土的膨胀量越小，收缩量也相对减小。当初始干密度从1.60g/cm³增加到1.80g/cm³时，在相同的加湿条件下，膨胀量降低了约40%。这是因为干密度较大时，土体颗粒排列更加紧密，孔隙较小，水分进入土体的难度增加，从而限制了膨胀变形。在收缩过程中，干密度大的土体颗粒间的连接力较强，抵抗收缩变形的能力也较强，所以收缩量相对较小。初始干密度还与缩胀速率密切相关。干密度较小的试样，在干湿变化时，水分的迁移速度较快，缩胀速率较大；而干密度较大的试样，水分迁移相对困难，缩胀速率较慢。在干湿循环过程中，干密度为1.60g/cm³的试样，其体积变化在较短时间内就能达到较大值，而干密度为1.80g/cm³的试样，体积变化相对较为缓慢，需要更长时间才能达到稳定状态。2.4.3试样形状对缩胀特性的影响探讨不同形状试样在缩胀试验中的表现发现，圆柱体试样和长方体试样在缩胀特性上存在一定差异。在相同的干湿条件下，长方体试样的边角部位更容易发生变形，因为边角处的水分散失和吸收速度更快，导致应力集中，从而产生较大的变形。而圆柱体试样由于形状较为规则，应力分布相对均匀，变形相对较为均匀。从体积变化量来看，在膨胀过程中，长方体试样的体积膨胀量略大于圆柱体试样；在收缩过程中，长方体试样的收缩量也相对较大。这是因为长方体试样的表面积与体积之比较大，水分交换面积大，所以在干湿变化时，体积变化更为明显。例如，在相同的加湿条件下，长方体试样的体积膨胀率比圆柱体试样高约10%。这种形状因素对缩胀特性的影响机制表明，在工程实际中，对于不同形状的膨胀土结构体，需要考虑其形状因素对变形的影响，合理进行设计和施工，以确保工程的稳定性。2.5本章小结通过对合肥膨胀土的加湿试验、收缩试验和缩胀试验研究，全面了解了其体变特性。加湿试验中，指数函数模型能较好地拟合加湿曲线，照相法在测量变形试样体积时精度更高；收缩试验确定了双曲线模型拟合收缩曲线的方法及缩限求解方式，失水路径对收缩特性有显著影响；缩胀试验揭示了初始干密度和试样形状对缩胀特性的影响规律。这些体变特性的研究结果对于深入理解合肥膨胀土的工程性质具有重要意义，为后续研究非饱和膨胀土在体变作用下的土水特征提供了关键的基础数据和理论支持，有助于准确评估膨胀土在工程中的变形行为，为工程设计和施工提供科学依据。三、合肥膨胀土的土水特性试验研究及曲线拟合3.1非饱和土力学相关理论基础非饱和土由固相（土颗粒）、液相（孔隙水）和气相（孔隙气）三相组成，其性质与各相的含量、分布以及相互作用密切相关。在非饱和土中，吸力是一个重要的概念，它反映了土中孔隙水的能量状态，是引起水分迁移和土体变形的重要因素。总吸力可分为基质吸力和溶质吸力两部分。基质吸力是由于土中孔隙水与孔隙气之间的弯液面产生的毛细作用而形成的，当土为非饱和时，孔隙中的水与毛细管中的水一样，存在弯液面，其上的水蒸气压力u_v，要小于土处于饱和状态（相同水质）时（即水面水平），水面上方的蒸气压力u_{v0}，表示毛细压力S的存在使相对湿度降低，毛细压力是产生吸力的重要因素，称为基质吸力。溶质吸力则是由于孔隙水中溶解的盐分等溶质而产生的，当水中含盐且水面水平时，水面上方的蒸气压力u_{v1}，又要小于纯净水水面水平时上方的蒸气压力，故土壤水含盐时，相对湿度会降低，即吸力会升高，与溶质种类及浓度有关，称为溶质吸力。对于一般工程问题，溶质吸力相对较小，可忽略不计，主要考虑基质吸力。土水特征曲线是描述非饱和土中基质吸力与含水率之间关系的曲线，它是研究非饱和土土水特性的重要工具。通常，土壤含水量以体积百分数表示，土壤吸力以大气压表示。土水特征曲线能够反映不同土壤的持水和释水特性，也可从中了解给定土类的一些土壤水分常数和特征指标。其曲线的斜率倒数称为比水容量，是用扩散理论求解水分运动时的重要参数。曲线的拐点可反映相应含水量下的土壤水分状态，如当吸力趋于0时，土壤接近饱和，水分状态以毛管重力水为主；吸力稍有增加，含水量急剧减少时，用负压水头表示的吸力值约相当于支持毛管水的上升高度；吸力增加而含水量减少微弱时，以土壤中的毛管悬着水为主，含水量接近于田间持水量；饱和含水量和田间持水量间的差值，可反映土壤给水度等。测量土水特征曲线的方法有多种，常见的有张力计法、压力膜仪法、离心机法等。张力计可直接测定土壤的水势，但压力范围有限，无法抵达植物的萎蔫系数；压力膜仪是目前默认的仪器设备，国际上一般都承认其测定数据，它通过高压气泵或高压氮气瓶向密封容器中充气加压，从土壤样品中渗出水分来测定含水量，从而得到土水特征曲线；离心机法则是利用离心力使土壤中的水分排出，根据离心力和水分排出量的关系确定土水特征曲线。3.2渗析法试验3.2.1渗析法试验原理与率定曲线渗析法试验基于溶质扩散原理来测量土水特征。其原理是利用半透膜的选择透过性，只允许溶剂分子（通常是水）通过，而溶质分子被截留在膜的一侧，从而实现溶质与溶剂的分离。在本试验中，将非饱和膨胀土试样与一定浓度的盐溶液通过半透膜隔开，由于土中孔隙水与盐溶液之间存在浓度差，水分会在浓度梯度的作用下通过半透膜进行扩散，当达到平衡状态时，根据盐溶液浓度的变化以及相关理论公式，就可以计算出土样的基质吸力和含水率。在进行正式试验之前，需要绘制渗析法的率定曲线。选用不同浓度的标准盐溶液，如氯化钠溶液，其浓度分别设置为0.1mol/L、0.2mol/L、0.3mol/L等。将已知含水率的纯净滤纸分别与这些不同浓度的盐溶液进行渗析平衡试验，在达到平衡后，测量滤纸的含水率。以盐溶液浓度为横坐标，滤纸平衡含水率为纵坐标，绘制出率定曲线。通过对率定曲线的分析，可以确定盐溶液浓度与滤纸平衡含水率之间的定量关系。率定曲线在试验中的作用至关重要，它为后续根据盐溶液浓度变化计算土样含水率提供了关键的依据。在实际试验中，当土样与盐溶液达到渗析平衡后，通过测量盐溶液浓度的变化，利用率定曲线就能够准确地确定土样在该吸力状态下的含水率。3.2.2渗析法试验仪器与步骤渗析法试验主要仪器包括渗析装置和高精度电子天平。渗析装置由有机玻璃制成的密封容器、半透膜（如醋酸纤维膜）、试样放置架等部分组成。半透膜将容器分隔为两个腔室，一侧放置土样，另一侧放置盐溶液。高精度电子天平用于准确称量土样和盐溶液的质量，精度可达0.001g。试验操作步骤如下：首先，将采集的合肥膨胀土土样制备成直径为50mm、高度为10mm的试样，用电子天平准确称取其初始质量m_1。将制备好的土样放入渗析装置的土样腔室中，在盐溶液腔室中加入一定浓度C_1、体积为V的盐溶液，并准确称取盐溶液和容器的总质量m_2。密封渗析装置，确保其密封性良好，防止水分蒸发和外界干扰。将渗析装置放置在恒温（25℃）、恒湿（相对湿度50%）的环境中，让土样与盐溶液进行渗析反应。每隔一定时间（如12h），轻轻摇晃渗析装置，使盐溶液浓度均匀，并观察渗析情况。经过一段时间（一般为3-5天），当土样与盐溶液达到渗析平衡时，打开渗析装置，迅速称取盐溶液和容器的总质量m_3，通过质量差计算出盐溶液质量的变化\Deltam=m_3-m_2。同时，取出土样，立即用电子天平称取其质量m_4，然后将土样放入105-110℃的烘箱中烘干至恒重，称取干土质量m_5，从而计算出土样的含水率w=\frac{m_4-m_5}{m_5}\times100\%。根据盐溶液质量的变化和体积，计算出平衡后盐溶液的浓度C_2，利用相关公式计算出土样的基质吸力S。在试验过程中，有诸多注意事项。半透膜的选择和安装至关重要，要确保半透膜无破损且安装紧密，否则会影响渗析效果和试验结果的准确性。试验环境的温度和湿度对渗析过程有显著影响，必须严格控制在设定范围内，以保证试验条件的一致性。在称取土样和盐溶液质量时，要快速、准确，尽量减少水分蒸发对质量测量的影响。3.3滤纸法试验3.3.1滤纸法试验原理与率定曲线滤纸法依据滤纸与土样水分平衡原理来测量吸力。其原理是当滤纸与非饱和土样紧密接触时，土样中的水分会在基质吸力的作用下向滤纸迁移，直到滤纸与土样的基质吸力相等，达到水分平衡状态。此时，滤纸的含水率与土样的基质吸力存在一定的对应关系。滤纸法的率定曲线绘制过程如下：选用定量滤纸，将其剪成直径为50mm的圆形滤纸。准备不同基质吸力的标准土样，可通过压力膜仪等设备制备。将圆形滤纸分别与这些不同吸力的标准土样紧密接触，放置在恒温（25℃）、恒湿（相对湿度50%）的环境中，让滤纸与土样进行水分平衡。达到平衡后，取出滤纸，立即放入密封袋中，防止水分蒸发。用精度为0.001g的电子天平称取滤纸的湿质量m_{湿}，然后将滤纸放入105-110℃的烘箱中烘干至恒重，称取干滤纸质量m_{干}，计算出滤纸的含水率w_{滤纸}=\frac{m_{湿}-m_{干}}{m_{干}}\times100\%。以标准土样的基质吸力为横坐标，滤纸平衡含水率为纵坐标，绘制率定曲线。对率定曲线进行分析可知，滤纸平衡含水率随着基质吸力的增大而减小，且呈现出一定的非线性关系。在低吸力范围内，滤纸含水率变化较为明显；随着吸力增大，滤纸含水率变化逐渐趋于平缓。这是因为在低吸力时，土样中的水分相对较多，向滤纸迁移的水分量较大，导致滤纸含水率变化显著；而在高吸力时，土样中可迁移的水分减少，滤纸含水率变化也就相应减小。通过率定曲线，可以建立起滤纸含水率与基质吸力之间的定量关系，为后续利用滤纸法测量土样吸力提供依据。3.3.2滤纸法试验仪器与步骤滤纸法试验主要用到的仪器和材料有定量滤纸、密封容器和电子天平。定量滤纸选用孔径均匀、吸水性良好的型号，如Whatman42号滤纸。密封容器采用玻璃材质的干燥器，其密封性良好，能够有效防止水分蒸发。电子天平精度为0.001g，用于准确称量滤纸和土样的质量。试验操作流程如下：将采集的合肥膨胀土土样制备成直径为50mm、高度为10mm的试样。在密封容器底部放置一层厚度约为5mm的海绵，海绵预先浸泡适量的水，以保持容器内的相对湿度稳定。将制备好的土样放置在海绵上方的试样架上，在土样上表面紧密覆盖一张圆形滤纸，确保滤纸与土样充分接触。盖上密封容器的盖子，将其放置在恒温（25℃）的环境中，让滤纸与土样进行水分平衡。每隔12h，打开密封容器，快速观察滤纸与土样的接触情况，确保滤纸未发生位移或脱落。经过48-72h，当认为滤纸与土样达到水分平衡时，迅速取出滤纸，放入密封袋中，用电子天平称取其湿质量m_{湿}。然后将滤纸放入105-110℃的烘箱中烘干至恒重，称取干滤纸质量m_{干}，计算出滤纸的含水率w_{滤纸}。根据事先绘制的率定曲线，由滤纸含水率查得对应的基质吸力S。同时，取出土样，立即用电子天平称取其质量m_{土湿}，然后将土样放入105-110℃的烘箱中烘干至恒重，称取干土质量m_{土干}，计算出土样的含水率w_{土}=\frac{m_{土湿}-m_{土干}}{m_{土干}}\times100\%。在数据处理方面，为了提高测量的准确性，对每个土样进行3次平行试验，取3次测量结果的平均值作为最终结果。计算每次测量结果与平均值的相对误差，若相对误差大于5%，则重新进行试验。通过对多个土样的测量和数据处理，得到不同含水率下土样的基质吸力数据，为后续分析土水特征提供数据支持。3.4土水特征拟合模型选择与应用在众多土水特征拟合模型中，选择Fredlund-Xing模型来拟合合肥膨胀土的土水特征曲线。Fredlund-Xing模型是一种广泛应用且适用性较强的模型，其表达式为：\theta=\frac{\theta_s}{(ln(e+\left(\frac{S}{S_0}\right)^a))^n}，式中\theta为体积含水率，\theta_s为饱和体积含水率，S为基质吸力，S_0为进气值吸力，a和n为拟合参数。应用该模型拟合试验数据的过程如下：将渗析法和滤纸法得到的基质吸力S和对应的含水率w数据进行整理，由于模型中使用的是体积含水率，需要将质量含水率w转换为体积含水率\theta，转换公式为\theta=w\times\rho_d/\rho_w，其中\rho_d为土的干密度，\rho_w为水的密度。利用最小二乘法对Fredlund-Xing模型进行参数优化，通过不断调整参数a、n和S_0的值，使模型计算得到的体积含水率\theta与试验测量值之间的误差平方和最小。例如，对于一组试验数据，经过多次迭代计算，最终确定拟合参数a=2.5，n=1.2，S_0=10。将确定的参数代入Fredlund-Xing模型中，得到拟合的土水特征曲线方程。将该方程绘制在坐标系中，与试验数据点进行对比，直观地展示模型的拟合效果。通过计算拟合优度R^2来评估拟合效果，R^2越接近1，说明模型对试验数据的拟合效果越好。经计算，该组数据的拟合优度R^2=0.92，表明Fredlund-Xing模型能够较好地拟合合肥膨胀土的土水特征曲线，能够准确地描述合肥膨胀土在不同吸力下的含水率变化规律。3.5试验结果分析对比渗析法和滤纸法的试验结果，发现两种方法得到的土水特征曲线在趋势上基本一致，都呈现出随着基质吸力的增大，含水率逐渐减小的规律。但在具体数值上存在一定差异，在低吸力范围内，渗析法测得的含水率略高于滤纸法；而在高吸力范围内，滤纸法测得的含水率略高于渗析法。这可能是由于两种方法的测量原理和试验过程存在差异导致的。渗析法是基于溶质扩散原理，在试验过程中，土样与盐溶液之间的扩散平衡可能受到半透膜性能、溶液浓度均匀性等因素的影响；滤纸法是基于水分平衡原理，滤纸与土样的接触状态、水分迁移过程中的阻力等因素可能对测量结果产生影响。分析土水特征曲线的特点，曲线呈现出明显的非线性关系。在低吸力阶段，含水率随吸力的变化较为缓慢，这是因为此时土中的孔隙主要被自由水占据，吸力的变化对自由水的影响较小；随着吸力的增大，进入中吸力阶段，含水率迅速下降，说明土中的毛细水开始大量排出；当吸力进一步增大到高吸力阶段，含水率下降趋势逐渐变缓，此时土中主要为结合水，结合水与土颗粒的结合力较强，较难被吸出。探讨影响土水特征的因素，土的颗粒组成是重要因素之一。合肥膨胀土中黏土矿物含量较高，黏土颗粒细小，比表面积大，对水分的吸附能力强，导致土水特征曲线较为平缓，在相同吸力下，含水率相对较高。初始干密度也对土水特征有显著影响，干密度越大，土颗粒排列越紧密，孔隙越小，水分迁移的通道变窄，在相同吸力下，含水率较低，土水特征曲线向左上方移动。此外，土的结构性也会影响土水特征，原状土具有一定的结构性，其土水特征与重塑土有所不同，原状土的结构能够在一定程度上阻碍水分的迁移，使得在相同吸力下，原状土的含水率相对较高。3.6本章小结通过渗析法和滤纸法对合肥膨胀土的土水特性进行试验研究，并运用Fredlund-Xing模型进行曲线拟合，得到了合肥膨胀土的土水特征曲线。试验结果表明，两种试验方法得到的曲线趋势一致，但存在数值差异，Fredlund-Xing模型能较好地拟合试验数据。土水特征曲线呈现非线性关系，受土的颗粒组成、初始干密度和结构性等因素影响。这些试验结果对于深入理解合肥膨胀土的土水特征具有重要意义，为后续研究非饱和膨胀土在体变作用下的土水特征提供了关键的试验数据和理论支持，有助于准确评估膨胀土在工程中的水理性质和稳定性，为工程设计和施工提供科学依据。四、膨胀土增湿过程中吸力-孔隙比-含水率关系4.1曲线拟合参数与初始孔隙比间关系在研究膨胀土增湿过程中的土水特征时，分析土水特征曲线（SWCC）和体积膨胀曲线（VSC）拟合参数与初始孔隙比的关系至关重要。对于SWCC拟合参数，以常用的Fredlund-Xing模型为例，其表达式为\theta=\frac{\theta_s}{(ln(e+\left(\frac{S}{S_0}\right)^a))^n}，其中a和n为拟合参数。通过对不同初始孔隙比的合肥膨胀土试样进行土水特征试验，发现随着初始孔隙比e_0的增大，参数a呈现出先增大后减小的趋势。当初始孔隙比从0.65增加到0.80时，a值从2.2增大到3.0，而后当e_0继续增大到0.90时，a值减小至2.5。这表明在初始孔隙比较小时，土体孔隙结构相对紧密，随着孔隙比增大，孔隙结构变得更加复杂，对吸力与含水率关系的影响呈现出特定的变化规律。参数n则随着初始孔隙比的增大而逐渐减小，当初始孔隙比从0.65增大到0.90时，n值从1.5减小到1.2。这说明初始孔隙比越大，含水率随吸力变化的速率越慢，土体的持水能力相对增强。在VSC拟合参数方面，采用的拟合模型如e=e_0+k_1w+k_2w^2（其中k_1和k_2为拟合参数）。研究发现，拟合参数k_1和k_2与初始孔隙比e_0也存在密切关系。随着e_0的增大，k_1和k_2的值均呈现出增大的趋势。当初始孔隙比从0.65增大到0.90时，k_1从0.1增加到0.3，k_2从0.05增加到0.12。这意味着初始孔隙比越大，在相同含水率变化下，孔隙比的变化幅度越大，土体的膨胀性越强。初始孔隙比的变化对土水特征有着多方面的影响。它会改变土体的孔隙结构，进而影响水分在土体中的储存和迁移。初始孔隙比的不同还会导致土水特征曲线的形态发生变化，从而影响土体的持水能力和强度特性等工程性质。4.2含水率、吸力及初始孔隙比的关系4.2.1重量含水率与吸力及初始孔隙比的关系通过对合肥膨胀土的一系列试验数据进行深入分析，建立了重量含水率w与吸力\psi及初始孔隙比e_0的数学关系。经研究发现，重量含水率与吸力之间存在着明显的非线性关系，且这种关系受到初始孔隙比的显著影响。基于试验数据，拟合得到的数学表达式为w=P_{11}e_0+P_{12}+\frac{P_{13}}{1+(\frac{\psi}{P_{14}})^{P_{15}}}+P_{16}，其中P_{11}至P_{16}为拟合参数。从该表达式可以看出，当吸力\psi固定时，重量含水率w随着初始孔隙比e_0的增大而增大。例如，当吸力为50kPa时，初始孔隙比从0.7增加到0.8，重量含水率从18%增加到22%。这是因为初始孔隙比越大，土体中的孔隙空间越大，能够容纳更多的水分。而当初始孔隙比e_0固定时，随着吸力\psi的增大，重量含水率w逐渐减小，呈现出指数衰减的趋势。这是由于吸力的增大使得土体对水分的吸附能力增强，水分逐渐被吸出，从而导致重量含水率降低。4.2.2孔隙比与重量含水率（或吸力）及初始孔隙比的关系在增湿过程中，孔隙比e与重量含水率w、吸力\psi和初始孔隙比e_0存在紧密的内在联系。通过试验数据拟合得到孔隙比与重量含水率及初始孔隙比的关系式为e=P_{21}e_0+P_{22}w+\frac{P_{23}}{1+(\frac{w}{P_{24}})^{P_{25}}}+P_{26}，其中P_{21}至P_{26}为拟合参数。从该式可以看出，当重量含水率w增大时，孔隙比e也随之增大，这是因为水分的增加使土体发生膨胀，孔隙体积增大。当初始孔隙比e_0增大时，在相同的重量含水率下，孔隙比e也会增大，这表明初始孔隙比大的土体在增湿过程中更容易发生膨胀。孔隙比与吸力及初始孔隙比也存在一定关系。随着吸力\psi的增大，孔隙比e逐渐减小，这是因为吸力增大导致水分排出，土体收缩，孔隙比降低。当初始孔隙比e_0固定时，吸力对孔隙比的影响更为显著。在低吸力范围内，孔隙比随吸力的变化较为明显；在高吸力范围内，孔隙比的变化逐渐趋于平缓。4.2.3体积含水率与吸力及初始孔隙比的关系为了建立体积含水率\theta与吸力\psi、初始孔隙比e_0的关系模型，首先将重量含水率w转换为体积含水率\theta，转换公式为\theta=w\times\frac{\rho_d}{\rho_w}，其中\rho_d为土的干密度，\rho_w为水的密度。通过对试验数据的分析和拟合，得到体积含水率与吸力及初始孔隙比的关系式为\theta=P_{31}e_0+P_{32}+\frac{P_{33}}{1+(\frac{\psi}{P_{34}})^{P_{35}}}+P_{36}，其中P_{31}至P_{36}为拟合参数。在工程应用中，该关系模型具有重要意义。在道路工程中，了解体积含水率与吸力及初始孔隙比的关系，可以预测路基在不同气候条件下的水分变化情况，从而合理设计路基的排水系统和防护措施，防止路基因水分变化而产生病害。在基坑工程中，利用该模型可以评估基坑周边土体的稳定性，为基坑支护设计提供重要依据。如果基坑周边土体的体积含水率过高，会导致土体强度降低，增加基坑坍塌的风险，通过该模型可以提前预测这种风险，采取相应的加固措施。4.3空间曲面的获取方法及分析获取吸力-孔隙比-含水率空间曲面主要有曲线拟合法和曲面拟合法两种方法。曲线拟合法的步骤较为复杂。以重量含水率w与吸力\psi、初始孔隙比e_0的关系曲面w=f_1(\psi,e_0)为例，首先分别对不同初始孔隙比时的w-\psi试验数据点采用合适的函数（如w=a+\frac{b}{1+(\frac{\psi}{c})^{d}}+e）进行曲线拟合，求得曲线拟合参数a，b，c，d，e。然后分别计算各初始孔隙比情况下拟合参数d及e的平均值，此平均值即为曲面拟合参数P_{15}及P_{16}。接着将d及e设为定值（各e_0的平均值），再分别对不同初始孔隙比时的w-\psi试验数据点采用上述函数进行曲线拟合，求得曲线拟合参数a，b，c。最后参考相关公式，分别拟合曲线a-e_0及b-e_0，得这两条曲线的拟合参数，它们即曲面拟合参数P_{11}至P_{14}。采用曲线拟合法的优点是便于分析各拟合参数与e_0的关系，能够深入了解每个参数随初始孔隙比的变化规律；缺点是需要进行多次曲线拟合，计算工作量偏大，过程较为繁琐。曲面拟合法相对简便。当空间曲面的具体表达式（如w=P_{11}e_0+P_{12}+\frac{P_{13}}{1+(\frac{\psi}{P_{14}})^{P_{15}}}+P_{16}）已经确定时，可直接对各初始孔隙比的所有试验数据点进行曲面拟合，一次性求得曲面拟合参数P_{11}至P_{16}。采用该方法可方便迅速地进行曲面拟合，其计算工作量远小于相应的曲线拟合法，大大提高了工作效率。将两种方法得到的结果进行对比，以w曲面为例，将两种方法得到的曲面绘于同一图形中，可以发现两种方法得到的曲面几乎重合，二者只在曲面的左上部略有不同。这表明两种方法的结果十分相近，但由于曲面拟合法计算简便，在后续分析中通常只取曲面拟合法的结果进行深入研究。4.4结果深入分析对体变曲面特征进行分析，在w-e_0-e坐标系中的体变曲面e=f_2(w,e_0)，可以看作由无数条不同e_0条件下的体积膨胀曲线（VSC）所组成。该曲面可以大致分为两部分，一部分是饱和部分，另一部分是非饱和部分。在饱和部分，随着重量含水率的增加，孔隙比基本保持不变，因为此时土体孔隙已被水充满，水分的继续增加不会导致孔隙比的明显变化；在非饱和部分，随着重量含水率的增加，孔隙比逐渐增大，土体发生膨胀。曲面由非饱和区进入饱和区的转折点对应的重量含水率随着初始孔隙比的增大而增大，转折点对应的吸力随着初始孔隙比的增大而减小。这是因为初始孔隙比越大，土体能够容纳更多的水分才会达到饱和状态，所以转折点对应的重量含水率增大，而吸力与含水率呈反比关系，所以转折点对应的吸力减小。在\psi-e_0-w或\psi-e_0-\theta坐标系中的土水特征曲面，与特定初始孔隙比对应的土水特征曲线是纵坐标恒定的平面曲线。这表明在固定初始孔隙比的情况下，吸力与重量含水率或体积含水率之间存在明确的对应关系，反映了土体在该初始孔隙比条件下的土水特性。在\psi-e-w或\psi-e-\theta坐标系中的土水特征曲面，与特定初始孔隙比对应的土水特征曲线是纵坐标在变化空间曲线，它能同时表示初始孔隙比的影响及试验过程中孔隙比的变化。随着吸力的变化，孔隙比和含水率也会相应改变，且这种变化受到初始孔隙比的影响。在低吸力阶段，孔隙比和含水率的变化相对较小；随着吸力增大，孔隙比逐渐减小，含水率也逐渐降低，且初始孔隙比越大，这种变化的幅度相对越大。这进一步揭示了土水特征与体变的耦合关系，即吸力的变化不仅影响含水率，还通过体变（孔隙比的变化）对土水特征产生综合影响。4.5本章小结通过对膨胀土增湿过程中吸力-孔隙比-含水率关系的研究，明确了SWCC和VSC拟合参数与初始孔隙比的关系，建立了重量含水率、孔隙比、体积含水率分别与吸力及初始孔隙比的定量关系模型。掌握了曲线拟合法和曲面拟合法获取空间曲面的方法，并对两种方法的结果进行了对比分析。深入分析了体变曲面和土水特征曲面的特征，揭示了土水特征与体变的耦合关系。这些研究成果对于深入理解非饱和膨胀土在体变作用下的土水特征具有重要意义，丰富了非饱和膨胀土的理论体系，为后续的工程应用和进一步研究提供了坚实的理论基础。五、土水特征试验过程的数值模拟5.1数值分析软件及数据来源本研究选用GeoStudio软件对土水特征试验过程进行数值模拟。GeoStudio是一款功能强大的地学模拟软件，专为处理涉及多孔介质流的工程问题而设计。其核心功能包括土壤水动力学模拟、地下水流动和热传输分析、以及多相流系统模拟等。在处理非饱和膨胀土的土水特征问题时，该软件能够准确模拟水在土壤中的运动，支持各种土壤水分特征曲线的输入，通过内置的算法和模型，能够精确地计算出水在膨胀土中的运移过程。其提供的多孔介质流的可视化模拟功能，可直观展示水分在土体中的分布和变化情况，帮助深入理解土壤水分运动的机制。数值模拟的数据主要来源于前文所述的试验结果。在体变特性试验中，获取的不同干湿循环次数、不同荷载条件下的体变数据，以及初始含水率、干密度等相关参数，为模拟膨胀土在不同工况下的变形提供了基础数据。土水特性试验中，渗析法和滤纸法得到的基质吸力与含水率数据，用于校准和验证数值模型，确保模型能够准确反映合肥膨胀土的土水特征。土样的基本土性指标，如液限、塑限、塑性指数、自由膨胀率、天然含水率、天然密度、孔隙比等，也是数值模拟中不可或缺的参数，它们决定了膨胀土的基本物理性质，影响着模型中材料参数的设定。5.2渗析法试验过程的数值模拟5.2.1模型几何尺寸与边界条件设定建立渗析法试验的数值模型时，充分考虑试验的实际情况来确定模型的几何尺寸。将模型设置为二维轴对称模型，以简化计算过程。模型的中心轴对应渗析装置的中心，模型的半径根据渗析装置中土样腔室的半径确定，设置为25mm，高度根据土样的高度设置为10mm，以准确模拟土样在渗析装置中的状态。在边界条件设定方面，水分边界条件至关重要。在模型的外边界，设置为定水头边界条件，根据试验中盐溶液的水位高度，将水头值设定为0.05m，以模拟盐溶液与土样之间的水分交换。在模型的顶部和底部，设置为不透水边界条件，防止水分从这些部位流失或进入，确保水分主要在土样与盐溶液之间进行渗析。应力边界条件的设置也不容忽视。在模型的顶部，施加一个与试验中实际压力相等的均布压力，根据试验加载情况，压力值设定为50kPa，以模拟土样在实际工程中可能承受的荷载。在模型的底部，设置为固定约束边界条件，限制土体在垂直方向的位移，模拟土体底部与渗析装置的接触状态。通过合理设定这些边界条件，能够更真实地模拟渗析法试验过程中土体的受力和水分迁移情况。5.2.2材料参数确定根据前面的试验结果，精确确定膨胀土的材料参数。渗透系数是影响水分迁移的关键参数，通过室内渗透试验，测得合肥膨胀土在不同含水率下的渗透系数。在数值模拟中，根据试验土样的初始含水率，选取对应的渗透系数值，一般在10⁻⁸-10⁻⁹m/s之间。弹性模量反映了土体抵抗弹性变形的能力，根据室内压缩试验结果，结合经验公式计算得到膨胀土的弹性模量。在不同的初始干密度和含水率条件下，弹性模量会有所变化，在数值模拟中，根据实际土样的初始状态，取值范围一般在5-15MPa。泊松比表示土体在横向应变与纵向应变之比，通过相关试验和理论分析，确定合肥膨胀土的泊松比取值为0.35。土水特征曲线参数则根据试验数据，利用Fredlund-Xing模型拟合得到，如进气值吸力、拟合参数a和n等，确保模型能够准确描述膨胀土在不同吸力下的含水率变化。5.2.3模拟结果分析将数值模拟结果与试验结果进行详细对比分析。在含水率变化方面，模拟结果显示，随着渗析时间的增加，土样的含水率逐渐增加，这与试验中观察到的现象一致。在渗析初期，含水率增加较快，后期逐渐趋于平缓，模拟结果与试验数据的变化趋势基本相符。通过计算模拟值与试验值之间的相对误差，发现大部分数据点的相对误差在10%以内，表明模拟结果具有较高的准确性。在基质吸力变化方面，模拟结果表明，随着渗析的进行，基质吸力逐渐减小，这与试验结果趋势一致。在低吸力阶段，模拟值与试验值较为接近；在高吸力阶段，由于模型简化和参数不确定性等因素，模拟值与试验值存在一定偏差，但偏差在可接受范围内。数值模拟结果对试验具有重要的补充和验证作用。通过模拟，可以直观地展示渗析过程中水分在土体中的迁移路径和分布情况，这在试验中难以直接观测到。模拟还可以预测不同工况下的试验结果，为进一步优化试验方案提供参考，有助于更深入地理解渗析法试验过程中膨胀土的土水特征变化规律。5.3滤纸法试验过程的数值模拟5.3.1模型几何尺寸与边界条件设定构建滤纸法试验的数值模型时，考虑到滤纸与土样的实际接触情况，将模型设置为二维平面模型。模型的长度根据土样的直径确定，设置为50mm，宽度根据滤纸和土样的总厚度设置为15mm，以准确模拟滤纸与土样的相互作用。在水分边界条件设置上，在模型的四周，设置为相对湿度边界条件。根据试验环境的相对湿度，将边界相对湿度设定为50%，以模拟试验过程中水分在滤纸与土样之间的平衡过程。在模型内部，设置滤纸与土样的界面为水分交换边界，允许水分在两者之间自由迁移。对于应力边界条件，在模型的底部，设置为固定约束边界条件，限制土体在垂直方向的位移，模拟土体底部与试验容器的接触状态。在模型的顶部，设置为自由边界条件，允许滤纸和土样在垂直方向自由变形，以符合试验实际情况。通过合理设置这些边界条件，能够较好地模拟滤纸法试验过程中土体的水分迁移和受力状态。5.3.2材料参数确定确定用于数值模拟的滤纸和膨胀土的材料参数。滤纸的渗透系数根据滤纸的材质和孔径等特性，通过相关资料和试验测定，取值为10⁻⁵m/s左右，以反映滤纸的水分传导能力。滤纸的孔隙率根据滤纸的结构，设定为0.7，用于描述滤纸内部的孔隙空间。膨胀土的材料参数在滤纸法试验模拟中与渗析法试验模拟有所不同。由于滤纸法试验中土体的受力状态和水分迁移机制与渗析法存在差异，在弹性模量方面，根据滤纸法试验的特点和相关研究，取值范围调整为3-10MPa。泊松比根据实际情况，仍取值为0.35。土水特征曲线参数同样根据滤纸法试验数据，利用Fredlund-Xing模型重新拟合确定，以适应滤纸法试验中膨胀土的土水特征。5.3.3模拟结果分析对滤纸法试验的数值模拟结果进行深入分析，将其与试验数据进行详细对比。在基质吸力方面，模拟结果显示，随着试验时间的增加，基质吸力逐渐减小，这与试验结果趋势一致。在低吸力范围内，模拟值与试验值的偏差较小，平均相对误差在8%以内；在高吸力范围内，由于模型对复杂物理过程的简化，模拟值与试验值的偏差略有增大，但仍在可接受范围内，平均相对误差在15%以内。在含水率变化方面，模拟结果表明，滤纸和土样的含水率在试验过程中逐渐趋于平衡，这与试验观察到的现象相符。模