非饱和粉砂临界状态特性及试验研究：理论、方法与工程应用

非饱和粉砂临界状态特性及试验研究：理论、方法与工程应用一、引言1.1研究背景与意义在各类工程建设中，非饱和粉砂是极为常见的一种土体类型，广泛分布于河流、湖泊、海岸等区域。在道路工程中，路基的填筑材料常涉及非饱和粉砂；在水利工程里，大坝、堤防等的基础也可能处于非饱和粉砂地层之上。此外，在城市建设中的基坑开挖、地下工程等，也频繁遇到非饱和粉砂的相关问题。然而，非饱和粉砂的力学性质较为复杂，受到多种因素的综合影响。与饱和土相比，非饱和土中存在气相，使得其内部应力状态更为复杂，存在基质吸力等特殊力学现象。基质吸力是指土中孔隙水压力与孔隙气压力之差，它对非饱和土的强度和变形特性有着显著影响。当土体含水率发生变化时，基质吸力也会随之改变，进而影响土体的力学性能。临界状态是土体力学中的一个重要概念，指土体在特定的应力和变形条件下，达到一种稳定的流动状态。在临界状态下，土体的变形持续发展，而应力比保持恒定，此时土体的性质和行为发生显著变化。对于非饱和粉砂而言，研究其临界状态具有至关重要的意义。在工程设计方面，准确掌握非饱和粉砂的临界状态，能够为地基承载力的计算提供科学依据。通过研究临界状态下非饱和粉砂的强度特性，可以合理确定地基的承载能力，避免因地基承载力不足而导致建筑物的沉降、倾斜甚至破坏等问题。在边坡稳定性分析中，了解非饱和粉砂在临界状态下的力学行为，有助于准确评估边坡的稳定性，为边坡的防护和加固设计提供可靠的理论支持，防止边坡失稳引发的地质灾害。在地震工程中，研究非饱和粉砂在地震荷载作用下达到临界状态的条件和响应，能够有效预测地震对建筑物和基础的影响，为抗震设计提供关键参数，提高建筑物的抗震能力。尽管当前针对非饱和土的研究已取得一定成果，但对于非饱和粉砂临界状态的研究仍有待深入。不同地区的非饱和粉砂，其颗粒级配、矿物成分等存在差异，导致力学性质不尽相同。现有的研究在考虑非饱和粉砂的复杂应力状态、多因素耦合作用等方面还存在不足，难以全面准确地描述其临界状态下的力学行为。因此，开展非饱和粉砂临界状态的试验研究具有重要的理论和实际意义。通过试验研究，能够深入揭示非饱和粉砂在不同应力路径、吸力条件下的力学特性，进一步完善非饱和土力学理论体系。同时，为工程实践中涉及非饱和粉砂的地基处理、边坡防护、地下工程等提供更为可靠的设计依据和技术支持，保障工程的安全稳定运行。1.2研究目的与内容本研究旨在通过系统的试验研究，深入揭示非饱和粉砂在不同应力路径和吸力条件下的临界状态规律，为非饱和土力学理论的完善和工程实践提供坚实的理论支持与可靠的数据依据。研究内容主要涵盖以下几个方面：首先，深入研究非饱和粉砂临界状态的相关理论基础。全面梳理非饱和土的基本力学特性，包括吸力、应力状态变量、强度与体变特性、水力特性等方面的研究进展，明确非饱和土临界状态的概念、特征及现有研究中存在的问题。对非饱和土三轴试验研究类型进行归纳总结，分析不同试验类型的特点和适用范围。同时，深入探讨非饱和土临界状态方程的发展历程和研究现状，为后续试验研究提供理论指导。其次，对试验技术和设备展开深入研究。吸力控制技术是研究非饱和粉砂临界状态的关键技术之一，详细研究轴平移技术、渗透技术和湿度控制技术的原理、优缺点及适用条件，为试验中吸力的精确控制提供技术支持。体变测量技术对于准确获取非饱和粉砂在试验过程中的变形特性至关重要，研究流体法测非饱和土体变、孔隙气相液相体积直接测量法和土样尺寸直接测量法等体变测量技术，对比分析各种方法的精度和可靠性。此外，对非饱和土三轴试验系统（GDS）进行全面介绍，包括三轴压力室、体变测量装置、数字型压力/体积控制器、吸力的测定与控制以及饱和陶土板等组成部分的工作原理和性能特点，确保试验设备能够满足研究需求。同时，对试验中使用的传感器进行校正，保证试验数据的准确性和可靠性。再者，开展试验材料及试验方案的研究。对试验材料的物理参数、颗粒级配和压实特性进行详细测定和分析，明确试验用非饱和粉砂的基本性质。制定全面合理的试验方案，包括非饱和土三轴剪切试验和饱和土三轴剪切试验。在非饱和土三轴剪切试验中，设置不同的吸力和围压条件，研究非饱和粉砂在不同应力路径下的力学特性；饱和土三轴剪切试验作为对比，用于分析饱和状态和非饱和状态下粉砂力学性质的差异。此外，进行非饱和土土-水特征曲线试验，测定非饱和粉砂在不同吸力下的含水率变化，为研究非饱和粉砂的水力特性提供数据支持。最后，对试验结果进行深入分析，并探讨非饱和粉砂临界状态在工程中的应用。对饱和粉砂和非饱和粉砂的试验结果进行详细分析，包括偏应力-轴向应变关系、体应变-轴向应变关系、孔隙水比体积-轴向应变关系、比体积-轴向应变关系和饱和度-轴向应变关系等，总结非饱和粉砂在不同应力路径和吸力条件下的力学特性变化规律。描述非饱和土在q-p'平面、v-lnp'平面和v_w-lnp'平面内的临界状态线，分析临界状态下q-s关系、v-s关系和v_w-s关系，建立非饱和粉砂临界状态的数学模型。基于试验结果，探讨非饱和粉砂临界状态在工程中的应用，如地基承载力计算、边坡稳定性分析和地震工程等方面，为实际工程提供理论依据和技术支持。1.3研究方法与创新点本研究综合运用多种研究方法，力求全面深入地探究非饱和粉砂的临界状态。在文献研究方面，全面梳理国内外相关文献资料，对非饱和土的基本力学特性、三轴试验研究类型、临界状态方程以及以往的临界状态试验研究等内容进行系统总结与分析。通过对已有研究成果的归纳，明确研究现状与存在的问题，为后续研究提供坚实的理论基础和研究方向指引。室内试验是本研究的关键环节。采用先进的非饱和土三轴试验系统（GDS），严格按照相关标准和规范进行试验操作。在试验过程中，精确控制吸力、围压等试验条件，设置多组不同的试验工况，以获取丰富且准确的试验数据。通过对试验数据的详细记录与分析，深入研究非饱和粉砂在不同应力路径和吸力条件下的力学特性，包括偏应力-轴向应变关系、体应变-轴向应变关系等。在理论分析上，基于试验结果，深入剖析非饱和粉砂的力学行为，探讨其临界状态的特征和规律。运用数学方法，建立非饱和粉砂临界状态的数学模型，对临界状态下的力学参数进行量化分析，进一步完善非饱和土力学理论体系。数值模拟方法也被应用于本研究中。借助专业的数值模拟软件，建立非饱和粉砂的数值模型，模拟其在不同工况下的力学响应。将数值模拟结果与试验结果进行对比分析，验证数值模型的准确性和可靠性，同时进一步拓展研究的深度和广度，为工程应用提供更具前瞻性的理论支持。本研究在试验技术和理论分析方面具有一定的创新点。在试验技术上，针对吸力控制和体变测量这两个关键技术难题，进行了深入研究和创新应用。在吸力控制技术方面，对比分析轴平移技术、渗透技术和湿度控制技术的优缺点及适用条件，根据试验需求选择合适的技术，并对其进行优化改进，实现了对吸力的高精度控制，为研究非饱和粉砂在不同吸力条件下的力学特性提供了可靠保障。在体变测量技术方面，研究并应用多种体变测量方法，如流体法测非饱和土体变、孔隙气相液相体积直接测量法和土样尺寸直接测量法等，通过对不同方法的对比分析，选择最适合本试验的测量方法，并对测量装置进行优化设计，提高了体变测量的精度和可靠性，为准确获取非饱和粉砂的变形特性提供了技术支持。在理论分析方面，基于试验结果，对非饱和粉砂在q-p'平面、v-lnp'平面和v_w-lnp'平面内的临界状态线进行了详细描述和分析，建立了非饱和粉砂临界状态的数学模型。该模型充分考虑了非饱和粉砂的复杂应力状态和多因素耦合作用，能够更准确地描述其临界状态下的力学行为，为非饱和土力学理论的发展做出了贡献。同时，将非饱和粉砂临界状态的研究成果应用于实际工程中，如地基承载力计算、边坡稳定性分析和地震工程等，提出了基于临界状态理论的工程设计方法和技术措施，为解决实际工程问题提供了新的思路和方法。二、非饱和土及临界状态理论基础2.1非饱和土基本力学特性2.1.1吸力吸力是理解非饱和土力学性质的关键因素，它反映了土中水的自由能状态。土中吸力通常由基质吸力和溶质吸力两部分构成，总吸力表达式为\psi=(u_a-u_w)+\pi，其中\pi表示渗透吸力，即溶质吸力。基质吸力是由于土中孔隙水与空气交界面处存在弯液面，弯液面的表面张力使得孔隙水产生负压力，从而形成吸力。在实际工程中，对于一般的非饱和土，孔隙水中的盐分含量相对较低，溶质吸力相较于基质吸力往往较小，对土体力学性质的影响也相对较弱，因此在多数情况下可以忽略不计，主要考虑基质吸力的作用。基质吸力的形成机制与毛细现象密切相关。非饱和土的孔隙中同时存在水和空气，在水-气分界面（收缩膜）处，由于水分子受力不平衡，内部水分子承受各向等值的力作用，而收缩膜内的水分子有一指向水体内部的不平衡力作用，为保持平衡，收缩膜内必须产生张力，这就导致了基质吸力的产生。当土中的含水率发生变化时，水-气分界面的形状和位置也会相应改变，进而引起基质吸力的变化。当土体含水率降低时，孔隙中的水体积减小，水-气分界面的曲率增大，根据拉普拉斯公式p_c=\frac{2\sigma\cos\theta}{r}（其中p_c为毛细压力，即基质吸力，\sigma为水的表面张力，\theta为接触角，r为弯液面的曲率半径），曲率半径r减小，基质吸力增大。测量非饱和土吸力的方法众多，每种方法都有其独特的原理、适用范围和优缺点。压力板法是一种较为常用的测量方法，其操作过程为将土样放置在饱和的陶瓷板上，通过在土样上方施加气体压力，使土样中的水分在压力差的作用下逐渐排出。当土样中的水分排出达到平衡状态时，所施加的气体压力与土样中的基质吸力相等，此时通过测量气体压力即可得到基质吸力。该方法操作相对简单，能够适用于不同类型土壤的吸力测量，但测量精度会受到土体密度、含水率等因素的影响。吸力杯法通过将土壤样品置于密封的吸力杯中，测量杯中气体压力变化来确定土壤吸力。该方法测量精度较高且稳定性较好，但需要使用精密的压力测量仪器，操作过程相对复杂，对操作人员的技术要求也较高。张力计法是通过测量土壤样品两侧的压力差来确定土壤吸力，适用于不同类型土壤的吸力测量，测量精度较高。然而，该方法的操作较为复杂，需要确保土壤样品的两端压力稳定，在实际应用中对实验条件的要求较为严格。吸力对非饱和土的力学性质有着显著影响。在强度方面，基质吸力的存在使得土颗粒间的有效应力增加，从而增强了土体的抗剪强度。许多学者通过试验研究发现，随着基质吸力的增大，非饱和土的抗剪强度呈现出上升的趋势。Fredlund和Mogerstern提出的双应力状态变量公式\tau_f=(\sigma-u_a)\tan\varphi'+(u_a-u_w)\tan\varphi^b就明确体现了基质吸力对非饱和土抗剪强度的贡献。在变形特性方面，吸力的变化会导致土体的体积发生改变。当吸力增大时，土颗粒间的吸引力增强，土体有收缩的趋势；反之，当吸力减小时，土体可能会发生膨胀。这种体积变化对非饱和土的工程性质有着重要影响，在地基沉降、边坡稳定性等工程问题中需要予以充分考虑。2.1.2应力状态变量非饱和土由于其三相体系的复杂性，其应力状态变量的描述相较于饱和土更为复杂。目前，学术界普遍认为非饱和土的应力状态不能仅用单一的有效应力来描述，而需要采用多个应力变量的组合。从非饱和土的土骨架平衡方程出发，运用热力学和连续介质力学的应力理论，可以推导出一系列应力状态变量。研究表明，非饱和土的应力状态变量存在多种组合形式，数量在5组以上，并非唯一。在众多应力状态变量组合中，净应力(\sigma-u_a)和吸力(u_a-u_w)的组合在工程实践中应用最为广泛。在大多数实际工程问题中，孔隙气压力往往等于大气压力，即压力表压力为零。采用这种组合方式能够清晰地区分总法向应力变化和孔隙水压力变化对土体力学性状的影响。净应力反映了总应力与孔隙气压力的差值，它对土体的骨架结构产生作用，影响土体的变形和强度。而吸力则体现了孔隙气压力与孔隙水压力的差值，其对土颗粒间的相互作用力和土体的体积变形有着重要影响。Bishop于1959年提出了非饱和土的有效应力公式\sigma'=\sigma+\chis=\sigma-u_w+\chis，其中\sigma'为非饱和土的有效应力，\sigma为总应力，u_a为孔隙气压力，s为吸力（s=u_a-u_w），\chi为与饱和度有关的参数。Bishop试图用这一单一的有效应力来描述非饱和土的强度和变形，但Jennings和Burland在1962年指出，该公式无法解释非饱和土中由浸水引起的湿陷现象。按照Bishop有效应力公式，当浸水湿化导致吸力减小时，有效应力也会减小，根据该理论计算得到的体积应该是膨胀的，但实际上大多数非膨胀性非饱和土在浸水湿化时体积是缩小的。这表明Bishop有效应力公式存在一定的局限性，不能全面准确地描述非饱和土的力学行为。Fredlund等在1977年运用连续介质力学方法证明，非饱和土的力学问题应采用净应力(\sigma-u_a)和吸力s作为其应力状态量，并提出了非饱和土的强度公式\tau_f=(\sigma-u_a)\tan\varphi'+(u_a-u_w)\tan\varphi^b，其中\tan\varphi'为法向净应力的内摩擦系数，\tan\varphi^b为吸力的内摩擦系数。该双应力状态变量理论考虑了非饱和土中孔隙气压力和孔隙水压力的共同作用，能够更合理地解释非饱和土的强度特性。然而，该理论也存在一定的局限性，它忽略了与孔隙率及饱和度的关联性，在描述非饱和土的某些力学现象时可能不够准确。不同应力状态变量在描述非饱和土力学行为时各有优劣。单一有效应力描述虽然形式简单，但无法准确解释非饱和土的复杂力学现象，如湿陷变形等。双应力状态变量描述考虑了孔隙气压力和孔隙水压力的影响，在解释非饱和土的强度特性方面具有一定优势，但在考虑孔隙率和饱和度等因素时存在不足。多应力状态变量描述虽然能够更全面地考虑非饱和土的各种因素，但由于变量众多，在实际应用中会增加计算的复杂性和难度。在选择应力状态变量时，需要综合考虑理论的合理性、应用的便利性以及逻辑关系的正确性等原则。同时，还应通过实验验证等方式，确保所选择的应力状态变量能够准确描述非饱和土的力学行为。2.1.3强度与体变特性非饱和土的强度特性是其力学性质的重要方面，受到多种因素的综合影响。基质吸力作为非饱和土特有的力学参数，对其强度有着显著影响。如前所述，基质吸力的存在使得土颗粒间的有效应力增加，从而增强了土体的抗剪强度。许多学者通过大量的室内试验和现场测试，深入研究了基质吸力与非饱和土抗剪强度之间的关系。研究结果表明，随着基质吸力的增大，非饱和土的抗剪强度呈现出明显的上升趋势。Fredlund和Mogerstern提出的双应力状态变量公式\tau_f=(\sigma-u_a)\tan\varphi'+(u_a-u_w)\tan\varphi^b，清晰地体现了基质吸力对非饱和土抗剪强度的贡献。土颗粒的性质，如颗粒大小、形状、级配等，对非饱和土的强度也有着重要影响。颗粒较大、形状不规则且级配良好的土体，其颗粒间的摩擦力和咬合力较大，能够提供更高的抗剪强度。而颗粒细小、级配不良的土体，其抗剪强度相对较低。土体的含水率同样是影响非饱和土强度的关键因素之一。当含水率发生变化时，基质吸力也会相应改变，进而影响土体的强度。在一定范围内，随着含水率的增加，基质吸力减小，非饱和土的抗剪强度降低；当含水率超过某一临界值后，土体可能会趋近饱和状态，此时抗剪强度的变化规律将与饱和土相似。非饱和土的体变特性同样复杂，受到多种因素的共同作用。应力路径对非饱和土体变有着重要影响。在不同的应力路径下，如等向压缩、剪切等，土体的体变响应各不相同。在等向压缩过程中，随着压力的增加，土体孔隙体积减小，发生压缩变形。而在剪切过程中，土体不仅会发生体积变化，还可能出现剪胀或剪缩现象。当土体受到剪切作用时，如果土颗粒之间能够相互错动并重新排列，使得孔隙体积减小，就会发生剪缩现象；反之，如果土颗粒之间相互挤密，导致孔隙体积增大，则会出现剪胀现象。吸力的变化也会引起非饱和土体变。当吸力增大时，土颗粒间的吸引力增强，土体有收缩的趋势，表现为体积减小；当吸力减小时，土体可能会发生膨胀，体积增大。这种由于吸力变化引起的体变现象在非饱和土的工程应用中需要特别关注，如在地基处理、填方工程等中，吸力的变化可能导致土体的体积不稳定，进而影响工程的稳定性和安全性。目前，描述非饱和土强度和体变的理论和模型众多。Mohr-Coulomb准则是经典的土的强度准则，在非饱和土强度研究中也有广泛应用。该准则认为，土体的抗剪强度由两部分组成，即黏聚力和摩擦力，表达式为\tau_f=c+\sigma_n\tan\varphi，其中\tau_f为抗剪强度，c为黏聚力，\sigma_n为法向应力，\varphi为内摩擦角。然而，Mohr-Coulomb准则在描述非饱和土强度时，需要对其参数进行修正，以考虑基质吸力等因素的影响。一些学者针对非饱和土的特性，提出了更为复杂的本构模型，如巴塞罗那基本模型（BBM）。该模型采用净应力和吸力作为应力状态量，能够描述非膨胀性非饱和土的基本力学特性，如屈服应力随吸力增加而变大、因湿化而引起的湿陷变形等。BBM模型的一个重要部分是荷载湿陷屈服曲线（LC屈服线），它反映了非饱和土的屈服应力随吸力增大而变大的特点，通过该曲线可以预测非饱和土最重要的变形特性——湿陷变形。然而，这些复杂的本构模型在实际应用中也存在一些局限性，如模型参数较多，确定过程较为复杂，且模型的适用范围有限，对于某些特殊的非饱和土或复杂的工程条件，可能无法准确描述其力学行为。2.1.4水力特性非饱和土的水力特性是其区别于饱和土的重要特性之一，对其工程性质有着重要影响。非饱和土的水力特性主要包括水分运移和渗透性两个方面。水分运移在非饱和土中是一个复杂的过程，受到多种因素的影响。基质吸力在水分运移中起着关键作用，它决定了水分在土体中的流动方向和速度。由于基质吸力的存在，水分会从吸力较小的区域向吸力较大的区域移动。当土体中某一区域的含水率较高，吸力较小时，水分会向含水率较低、吸力较大的区域流动，以达到能量平衡状态。土的孔隙结构也对水分运移有着重要影响。孔隙大小、形状和连通性等因素决定了水分在孔隙中的流动路径和阻力。孔隙较大且连通性良好的土体，水分运移相对容易，速度较快；而孔隙细小、连通性差的土体，水分运移则较为困难，速度较慢。温度对非饱和土的水分运移也有一定影响。温度升高会使水分的蒸发和扩散速度加快，从而影响水分在土体中的分布和运移。在炎热的夏季，土体表面的水分蒸发较快，导致表面吸力增大，水分会从土体内部向表面运移。非饱和土的渗透性与饱和土有着显著差异。在饱和土中，孔隙被水完全充满，渗透性主要取决于孔隙的大小和连通性。而在非饱和土中，由于孔隙中存在气体，气体的存在阻碍了水分的流动，使得非饱和土的渗透性明显低于饱和土。非饱和土的渗透性还与含水率密切相关。随着含水率的降低，土体中连通的孔隙水通道减少，渗透性急剧下降。当含水率降低到一定程度时，土体中的水分主要以薄膜水和吸附水的形式存在，几乎不发生渗透现象。研究非饱和土水力特性的方法主要有室内试验和数值模拟两种。室内试验方法包括压汞法、毛细压力计法、干度法等。压汞法是一种高精度的测量非饱和土构成的孔隙系统的方法，通过使用压力和压力变化来测量非饱和土样品内毛细孔的孔径和孔隙度。毛细压力计法是一种常用的测定土壤含水量和气水界面位置的方法，其原理基于毛细作用的基本概念，是测定非饱和土体积含水量和气水相平衡位置的可靠方法。干度法是一种非饱和土渗透性的测定方法，利用非饱和土样品在特定压力条件下排出水分的时间和流量来确定其渗透性系数。数值模拟方法则是利用计算机软件对非饱和土的水力特性进行模拟分析。常用的数值模拟软件如HYDRUS、SEEP/W等，这些软件基于达西定律和土-水特征曲线等理论，能够模拟非饱和土中水分的运移和渗透过程。通过数值模拟，可以快速、准确地获取非饱和土在不同条件下的水力参数，为工程设计和分析提供依据。然而，数值模拟结果的准确性依赖于输入参数的准确性和模型的合理性，因此在应用数值模拟方法时，需要对模型进行验证和校准。非饱和土的水力特性与临界状态之间存在着密切的联系。在临界状态下，土体的结构和孔隙特征会发生变化，进而影响其水力特性。当土体达到临界状态时，孔隙结构可能会发生重塑，孔隙大小和连通性改变，这将直接影响水分的运移和渗透性。在边坡失稳等接近临界状态的情况下，土体的渗透性可能会发生突变，导致水分快速聚集或散失，进一步加剧边坡的失稳。因此，在研究非饱和土的临界状态时，需要充分考虑其水力特性的变化。2.2临界状态理论概述临界状态理论的发展历程可以追溯到20世纪中叶，它是土力学领域中的一个重要理论，为深入理解土体的力学行为提供了关键的框架。早在1958年，以Roscoe为代表的剑桥学派在土力学研究中取得了重要突破，他们通过一系列开创性的研究工作，首次提出了临界状态的概念。这一概念的提出，犹如一颗启明星，为后续土力学的发展指明了方向。Roscoe及其团队通过对土体在不同应力条件下的变形和强度特性进行深入研究，发现土体在受到外力作用时，其变形和应力状态会逐渐发展变化，当达到某一特定状态时，土体的变形将持续发展，而应力比保持恒定，此时土体进入了一种特殊的稳定流动状态，即临界状态。这一发现揭示了土体力学行为的重要规律，为土力学的理论研究和工程应用奠定了坚实的基础。1963年，Roscoe等人进一步深化了对临界状态的研究，他们详细阐述了土体在临界状态下的各种特性。研究表明，在临界状态下，土体的应力-应变关系呈现出独特的特征，土体的体积变形和剪切变形之间存在着紧密的联系。通过对大量试验数据的分析和总结，他们提出了描述临界状态的数学表达式，为定量研究土体的临界状态提供了有力的工具。这些研究成果在土力学领域引起了广泛的关注和深入的讨论，推动了临界状态理论的不断发展和完善。随着时间的推移，临界状态理论在土力学领域得到了广泛的应用和深入的研究。众多学者在Roscoe等人的研究基础上，从不同角度对临界状态理论进行了拓展和深化。一些学者通过改进试验技术和方法，获取了更精确的试验数据，进一步验证和完善了临界状态理论。他们对不同类型的土体进行了大量的试验研究，包括砂土、粉土、黏土等，深入探讨了土体的颗粒组成、孔隙结构、初始应力状态等因素对临界状态的影响。另一些学者则从理论分析的角度出发，运用数学和力学方法，建立了更加完善的临界状态本构模型。这些模型能够更准确地描述土体在复杂应力条件下的力学行为，为工程实践提供了更可靠的理论支持。在工程实践中，临界状态理论被广泛应用于解决各类实际问题。在地基承载力计算方面，临界状态理论为确定地基的承载能力提供了科学的依据。通过研究地基土在临界状态下的强度和变形特性，可以合理评估地基的承载能力，确保建筑物的安全稳定。在边坡稳定性分析中，临界状态理论有助于准确判断边坡土体是否处于临界状态，预测边坡的失稳可能性。通过对边坡土体的应力状态和变形特性进行分析，运用临界状态理论可以制定合理的边坡防护和加固措施，防止边坡失稳引发的地质灾害。在地下工程中，如隧道、基坑等，临界状态理论可以用于分析土体在开挖过程中的力学响应，指导工程的设计和施工。通过预测土体在临界状态下的变形和破坏模式，可以合理选择支护结构和施工方法，确保地下工程的安全顺利进行。临界状态是指土体在外部荷载作用下，经过一系列的变形和应力调整后，达到的一种特殊状态。在这种状态下，土体的变形持续发展，而应力比保持恒定。具体来说，当土体达到临界状态时，其偏应力与平均有效应力的比值（即应力比）不再随变形的增加而变化，保持一个相对稳定的值。同时，土体的体积变形也会达到一个相对稳定的状态，不再发生明显的变化。在三轴试验中，当土体达到临界状态时，偏应力-轴向应变曲线会趋于平缓，体应变-轴向应变曲线也会达到一个相对稳定的阶段。临界状态的特征主要包括以下几个方面：首先，应力比恒定是临界状态的重要特征之一。在临界状态下，土体的应力比保持不变，这意味着土体的抗剪强度达到了一个相对稳定的值。这种应力比的恒定特性使得临界状态成为土体力学行为的一个重要转折点，对于理解土体的强度和变形特性具有重要意义。其次，变形持续发展也是临界状态的显著特征。尽管应力比保持恒定，但土体的变形会持续进行，这表明土体在临界状态下处于一种不稳定的流动状态。这种变形持续发展的特性与土体的内部结构和颗粒间的相互作用密切相关。此外，临界状态下土体的孔隙比和含水率等物理参数也会达到一个相对稳定的值。这些物理参数的稳定反映了土体内部结构的相对稳定性，同时也对土体的力学性质产生重要影响。在土力学中，临界状态理论有着广泛的应用。在地基承载力计算中，临界状态理论为确定地基的极限承载力提供了重要的理论依据。根据临界状态理论，当基础底面的压力达到一定值时，地基土会进入临界状态，此时地基的承载能力达到极限。通过对地基土在临界状态下的力学特性进行分析，可以准确计算地基的极限承载力，为基础设计提供可靠的参数。在边坡稳定性分析方面，临界状态理论可以用于判断边坡土体是否处于临界状态，预测边坡的稳定性。当边坡土体的应力状态接近临界状态时，边坡的稳定性会显著降低，容易发生失稳破坏。通过运用临界状态理论对边坡土体的应力和变形进行分析，可以及时发现潜在的边坡失稳风险，并采取相应的防护和加固措施。在数值模拟中，临界状态理论为建立准确的土体本构模型提供了基础。土体本构模型是描述土体力学行为的数学模型，它对于预测土体在不同荷载条件下的变形和强度具有重要作用。基于临界状态理论建立的本构模型能够更好地反映土体的实际力学行为，提高数值模拟的准确性和可靠性。2.3非饱和土临界状态研究现状国外学者在非饱和土临界状态研究方面开展了大量工作，并取得了一系列重要成果。早在20世纪60年代，Roscoe等人就对饱和黏土的临界状态进行了深入研究，提出了剑桥模型，为后续非饱和土临界状态的研究奠定了基础。随着研究的不断深入，学者们逐渐将研究重点转向非饱和土。Alonso等在1990年提出了巴塞罗那基本模型（BBM），该模型采用净应力和吸力作为应力状态量，能够描述非膨胀性非饱和土的基本力学特性，如屈服应力随吸力增加而变大、因湿化而引起的湿陷变形等。BBM模型的提出，推动了非饱和土临界状态研究的发展，成为非饱和土本构模型研究的经典。Wheeler和Sivakumar在1995年通过吸力控制的非饱和土等向压缩和三轴试验，对BBM模型的合理性进行了验证。他们的研究结果表明，BBM模型能够较好地描述非饱和土在不同吸力条件下的力学行为，为非饱和土临界状态的研究提供了重要的实验依据。随后，许多学者在BBM模型的基础上进行了改进和拓展，提出了各种不同的非饱和土本构模型。这些模型在考虑非饱和土的复杂应力状态、多因素耦合作用等方面取得了一定进展，能够更准确地描述非饱和土的临界状态。国内学者在非饱和土临界状态研究方面也取得了显著成果。陈正汉等采用类似饱和土的Duncan-Chang模型的手法，建立了非饱和土的非线性弹性本构模型。该模型在一定程度上能够描述非饱和土的力学特性，但在考虑非饱和土的剪胀剪缩性、湿陷变形等方面还存在不足。沈珠江提出了一个适用于非饱和土的统一弹塑性本构模型，该模型基于能量原理，能够较好地描述非饱和土的各种力学特性，包括强度、变形和湿陷等。模型在理论上具有一定的创新性，但在实际应用中，由于模型参数较多，确定过程较为复杂，限制了其广泛应用。尽管国内外学者在非饱和土临界状态研究方面取得了一定成果，但目前仍存在一些问题和不足。不同学者提出的临界状态模型和理论往往基于特定的试验条件和假设，其适用范围存在一定局限性。许多模型在描述非饱和土的复杂应力状态和多因素耦合作用时还不够完善，难以全面准确地反映非饱和土的临界状态。在试验研究方面，由于非饱和土的力学性质受到多种因素的影响，试验条件的控制较为困难，导致试验结果的离散性较大，影响了研究结果的准确性和可靠性。此外，目前对于非饱和土临界状态的研究主要集中在室内试验和理论分析方面，现场试验研究相对较少，这使得研究成果在实际工程中的应用受到一定限制。针对上述问题，未来的研究可以从以下几个方向展开：进一步完善非饱和土临界状态的理论和模型，充分考虑非饱和土的复杂应力状态、多因素耦合作用以及土体的结构性等因素，提高模型的适用性和准确性。加强试验技术的研究和改进，提高试验条件的控制精度，减少试验结果的离散性，为理论研究提供更可靠的试验数据。加大现场试验研究的力度，通过现场试验验证和完善理论研究成果，提高研究成果在实际工程中的应用价值。开展多学科交叉研究，将非饱和土临界状态的研究与岩土工程、环境科学、地震工程等学科相结合，拓展研究的广度和深度，为解决实际工程问题提供更全面的理论支持。三、非饱和粉砂临界状态试验技术与设备3.1试验技术3.1.1吸力控制技术吸力控制技术是研究非饱和粉砂临界状态的关键技术之一，它对于准确模拟非饱和粉砂在实际工程中的受力状态和水分条件具有重要意义。目前，常用的吸力控制技术主要包括轴平移技术、渗透技术和湿度控制技术，每种技术都有其独特的原理、优缺点及适用条件。轴平移技术是应用最为广泛的吸力控制技术之一，其原理基于土-水特征曲线。该技术通过高进气值陶土板来实现吸力的控制。高进气值陶土板具有许多均匀的微细孔，当陶土板完全饱和时，在其表面形成收缩膜，收缩膜产生的表面张力能够阻止空气通过陶土板。通过控制孔隙气压力和孔隙水压力，使得二者的差值等于目标吸力。具体操作时，将土样放置在饱和的高进气值陶土板上，土样中的孔隙水通过陶土板与孔隙水压力量测系统相连，孔隙气压力则作用在土样表面。通过调整孔隙气压力和孔隙水压力，即可实现对土样吸力的精确控制。轴平移技术的优点是能够在一定范围内精确控制吸力，操作相对简便，适用于大多数非饱和土的试验研究。然而，该技术也存在一定的局限性，高进气值陶土板的进气值限制了其可控制的吸力范围，一般适用于吸力小于1.5MPa的情况。此外，陶土板的使用可能会对土样的初始状态产生一定的扰动，影响试验结果的准确性。渗透技术是利用渗透吸力来控制土样的吸力。该技术通过在土样周围设置渗析溶液，利用渗析溶液的渗析力对土样施加渗透吸力。不同浓度的渗析溶液会产生不同大小的渗透吸力，当土样中的基质吸力与渗析溶液的渗透吸力达到平衡时，土样的水分状态保持稳定，从而实现对吸力的控制。渗透技术的优点是可以实现较高吸力的控制，适用于研究高吸力条件下非饱和土的力学性质。而且，该技术对土样的扰动较小，能够更好地保持土样的原始结构。但渗透技术的缺点是试验装置较为复杂，需要精确控制渗析溶液的浓度和流量，试验操作难度较大。此外，渗析过程中可能会发生溶质的扩散和交换，影响试验结果的稳定性。湿度控制技术是通过控制土样周围环境的相对湿度来间接控制土样的吸力。根据土-水特征曲线，土样的吸力与周围环境的相对湿度存在一定的关系。通过调节环境的相对湿度，可以使土样达到相应的吸力状态。湿度控制技术的优点是可以实现连续的吸力控制，适用于研究吸力变化对非饱和土力学性质的影响。该技术对土样的扰动较小，能够模拟土样在自然环境中的湿度变化。然而，湿度控制技术的缺点是试验周期较长，因为土样与周围环境达到湿度平衡需要一定的时间。而且，湿度控制的精度受到环境条件的影响较大，在实际应用中需要严格控制环境的温度和湿度。在本研究中，根据试验的具体需求和条件，选择轴平移技术作为主要的吸力控制技术。轴平移技术能够满足本试验对吸力控制范围和精度的要求，且操作相对简便，便于在试验过程中进行控制和调整。同时，为了减少陶土板对土样的扰动，在试验前对陶土板进行了充分的饱和处理，并在安装土样时小心操作，尽量保持土样的初始状态。在试验过程中，密切监测孔隙气压力和孔隙水压力的变化，确保吸力的控制精度。3.1.2体变测量技术体变测量技术对于准确获取非饱和粉砂在试验过程中的变形特性至关重要，它能够为研究非饱和粉砂的力学行为提供关键数据。目前，常用的体变测量技术主要包括流体法测非饱和土体变、孔隙气相液相体积直接测量法和土样尺寸直接测量法，每种方法都有其独特的原理、优缺点。流体法测非饱和土体变是一种较为常用的体变测量方法，其原理基于质量守恒定律。该方法通过测量试验过程中流入或流出土样的流体体积来间接计算土样的体积变化。在三轴试验中，通常使用水或气体作为流体介质。当土样发生体积变化时，会导致流体的流动，通过测量流体的流量和时间，可以计算出流入或流出土样的流体体积，从而得到土样的体变。流体法测非饱和土体变的优点是测量原理简单，操作相对容易，能够满足一般试验的精度要求。该方法在工程实践中应用广泛，积累了丰富的经验。然而，该方法也存在一定的局限性，测量结果容易受到流体的压缩性、试验装置的密封性等因素的影响。如果流体存在压缩性，在测量过程中会导致测量误差；试验装置的密封性不好，会使流体泄漏，影响测量的准确性。孔隙气相液相体积直接测量法是直接测量土样中孔隙气和孔隙水的体积变化，从而计算出土样的总体积变化。该方法通常使用压力-体积控制器来实现对孔隙气和孔隙水体积的测量。在试验过程中，通过控制孔隙气压力和孔隙水压力，测量相应的体积变化。然后，根据孔隙气和孔隙水的体积变化，计算出土样的总体积变化。孔隙气相液相体积直接测量法的优点是测量精度较高，能够准确反映土样的体变情况。该方法对于研究非饱和土的微观结构变化和力学特性具有重要意义。但该方法的缺点是试验装置较为复杂，需要高精度的压力-体积控制器和良好的密封性能。试验过程中，需要对孔隙气和孔隙水的压力和体积进行精确控制和测量，操作难度较大。土样尺寸直接测量法是通过直接测量土样在试验过程中的尺寸变化来计算体变。该方法通常使用位移传感器或应变片等测量装置，直接测量土样的轴向和径向变形。然后，根据土样的初始尺寸和测量得到的变形量，计算出土样的体积变化。土样尺寸直接测量法的优点是测量直观，能够直接反映土样的变形情况。该方法对于研究土样的变形模式和变形机理具有重要作用。然而，该方法的缺点是测量精度受到测量装置的精度和土样表面平整度的影响。如果测量装置的精度不够高，或者土样表面不平整，会导致测量误差增大。在本研究中，综合考虑各种体变测量技术的优缺点和试验的具体需求，选择孔隙气相液相体积直接测量法作为主要的体变测量技术。该方法能够满足本试验对体变测量精度的要求，且能够准确反映非饱和粉砂在不同应力路径和吸力条件下的体变特性。为了确保测量的准确性，在试验前对压力-体积控制器进行了校准，保证其测量精度。在试验过程中，严格控制试验装置的密封性，避免孔隙气和孔隙水的泄漏。同时，对测量数据进行多次采集和处理，减小测量误差。3.2试验设备本研究采用GDS非饱和土三轴试验系统进行试验，该系统是英国GDS公司生产的一款高端非饱和土试验仪器，是对传统三轴试验的升级和扩充，主要用来研究地下水位以上的土的特性，可以模拟现场的应力状态和饱和情况。它在非饱和土力学特性研究领域具有重要地位，能够为准确揭示非饱和粉砂的临界状态提供可靠的试验平台。GDS非饱和土三轴试验系统主要由三轴压力室、体变测量装置、数字型压力/体积控制器、吸力的测定与控制以及饱和陶土板等部分组成。三轴压力室是试验系统的核心部件之一，它主要由压力室底座、有机玻璃筒、压力室顶盖、轴向加载杆和围压管路等组成。压力室底座用于支撑整个压力室和放置土样，有机玻璃筒则为试验提供了一个封闭的空间，便于控制和观察试验过程。压力室顶盖与有机玻璃筒紧密连接，确保压力室的密封性。轴向加载杆穿过压力室顶盖，与土样顶部相连，用于施加轴向荷载。围压管路环绕在压力室周围，通过向压力室内充入液体，实现对土样围压的施加和控制。在试验过程中，围压可以根据试验需求进行精确调整，模拟不同的现场应力状态。例如，在研究非饱和粉砂在不同深度的力学特性时，可以通过调整围压来模拟相应的上覆压力。体变测量装置是获取非饱和粉砂体变数据的关键部分，本试验系统采用孔隙气相液相体积直接测量法来测量体变。该装置主要由压力-体积控制器、管路和传感器等组成。压力-体积控制器通过管路与土样相连，能够精确测量土样中孔隙气和孔隙水的体积变化。在试验过程中，当土样受到荷载作用发生体积变化时，孔隙气和孔隙水的体积也会相应改变。压力-体积控制器通过监测管路中气体和液体的压力变化，根据相关原理计算出孔隙气和孔隙水的体积变化量，从而得到土样的总体积变化。这种测量方法具有较高的精度，能够准确反映非饱和粉砂在不同应力路径和吸力条件下的体变特性。数字型压力/体积控制器是试验系统的重要控制部件，它具有高精度的压力和体积控制能力。该控制器采用先进的数字控制技术，能够精确控制孔隙气压力、孔隙水压力、围压和轴压等参数。在试验前，根据试验方案在控制器上设置好各项参数，如目标吸力、围压大小、轴向加载速率等。在试验过程中，控制器会实时监测压力和体积的变化，并根据预设的参数自动进行调整，确保试验条件的稳定性和准确性。例如，在进行不同吸力条件下的三轴剪切试验时，数字型压力/体积控制器能够根据设定的吸力值，精确控制孔隙气压力和孔隙水压力，使土样达到所需的吸力状态。同时，它还能实时采集压力和体积数据，并将数据传输给计算机进行处理和分析。吸力的测定与控制是研究非饱和粉砂临界状态的关键环节，本试验系统采用轴平移技术来实现吸力的控制。通过高进气值陶土板来分隔孔隙气和孔隙水，高进气值陶土板具有许多均匀的微细孔，当陶土板完全饱和时，在其表面形成收缩膜，收缩膜产生的表面张力能够阻止空气通过陶土板。土样放置在饱和的高进气值陶土板上，土样中的孔隙水通过陶土板与孔隙水压力量测系统相连，孔隙气压力则作用在土样表面。通过控制孔隙气压力和孔隙水压力，使得二者的差值等于目标吸力。在试验过程中，通过数字型压力/体积控制器精确调节孔隙气压力和孔隙水压力，实现对吸力的精确控制。同时，利用高精度的孔隙水压力传感器和孔隙气压力传感器实时监测孔隙水压力和孔隙气压力的变化，确保吸力的控制精度。饱和陶土板在试验中起着重要的作用，它是实现轴平移技术控制吸力的关键部件。GDS试验系统中所用陶土板进气值为1500kPa，可以满足本试验对吸力控制范围的要求。在试验前，需要对陶土板进行饱和处理，以确保其能够正常工作。饱和处理的方法是将连接陶土板底部的细塑料管充满水，减少连接管路中的空气。然后将与反压控制器相连的管路和与孔隙水压力传感器相连的管路分别连接到陶土板底部。将与孔压传感器相连的管路出口浸没入装有无气水的器皿中，利用反压控制器施加30kPa的压力于陶土板底部（注意压力不能超过40kPa，否则会因弯矩过大而折断陶土板）。保持30kPa的压力直至陶土板上表面出现一薄层水膜，同时观察装有无气水的器皿中长时间没有气泡溢出，这个过程约需24小时。饱和后的陶土板能够有效地分隔孔隙气和孔隙水，为准确控制吸力提供保障。3.3传感器校正在试验过程中，传感器的准确性直接影响到试验数据的可靠性和试验结果的准确性。因此，在使用GDS非饱和土三轴试验系统进行试验之前，对系统中的传感器进行校正至关重要。本试验中主要涉及压力室体变测量装置中的压力-体积控制器相关传感器以及其他用于测量孔隙气压力、孔隙水压力和轴向荷载等的传感器。对于压力室体变测量装置中的压力-体积控制器相关传感器，采用高精度标准压力源和体积标准器进行校正。在零点校准阶段，将标准压力源设置为零压力，体积标准器设置为零体积变化，此时记录压力-体积控制器传感器的输出值。若输出值偏离零点，通过调整传感器内部的零点调节机构或软件设置，使输出值回到设定范围。例如，当压力传感器的输出值在零压力时显示为0.05kPa，而设定的零点范围为±0.01kPa，此时就需要通过调节机构将输出值调整到规定范围内。在量程校准阶段，将标准压力源设置为满量程压力，体积标准器设置为满量程体积变化，观察并记录传感器的输出值。若输出值偏离满量程，通过调整传感器内部的量程调节机构或软件设置，使输出值达到设定范围。比如，压力传感器满量程为2MPa，当标准压力源施加2MPa压力时，传感器输出值应为2MPa，但实际输出为1.98MPa，这时就需要进行调整。对于测量孔隙气压力、孔隙水压力和轴向荷载等的传感器，同样使用高精度的标准压力源和荷载标准器进行校正。在零点校准和量程校准的基础上，还需要进行线性度校验。在线性度校验时，在量程范围内选择多个测试点，如对于孔隙气压力传感器，在其量程0-1MPa范围内，选择0.2MPa、0.4MPa、0.6MPa、0.8MPa等测试点。使用标准压力源逐一施加压力，记录每个测试点传感器的输出值。通过计算各测试点的误差，评估传感器的线性度。若线性度不符合要求，可能需要进一步调整或维修。例如，计算得到某测试点的误差超出了允许范围，就需要对传感器进行检查和调整，确保其线性度满足试验要求。对于压力控制器，在使用前需进行严格的校准和检查。根据设备的接口规格，选择合适的连接方式和密封材料，避免气体或液体泄漏。安装完成后，接通电源，按照操作手册进行初步的设置和校准。设定压力值时，需根据实际应用场景和工艺要求来确定。通过控制面板或相关软件，可以方便地设定目标压力值。在设定过程中，要注意单位的统一，避免出现误差。定期对压力控制器进行检查和维护，确保其正常工作。例如，每隔一段时间检查压力控制器的显示是否正常，控制参数是否稳定等。在传感器校正过程中，需严格按照相关标准和操作规程进行操作，确保校正的准确性和可靠性。同时，详细记录校准过程中的各项数据，包括校准时间、校准环境、校准前后的传感器输出值、调整参数等信息。这些记录不仅可以作为传感器校准的依据，还能为后续的试验数据分析和设备维护提供重要参考。通过对传感器进行精确校正，可以有效提高试验数据的质量，为非饱和粉砂临界状态的研究提供可靠的数据支持。四、非饱和粉砂临界状态试验方案设计4.1试验材料本试验所用土样取自[具体地点]，该区域的粉砂具有典型的非饱和特性，广泛分布于当地的各类工程场地中，对其进行研究具有重要的工程实际意义。在土样采集过程中，严格遵循相关标准和规范，确保采集的土样具有代表性。使用专业的取土设备，在不同深度和位置进行多点采样，然后将采集到的土样混合均匀，以消除局部差异对试验结果的影响。对土样的物理参数进行了详细测定，测定结果如表1所示。土样的液限为[X]%，塑限为[X]%，塑性指数为[X]，表明该土样具有一定的可塑性，但可塑性相对较弱。最大干密度为[X]g/cm³，最小干密度为[X]g/cm³，这两个参数对于评估土样的压实性能和工程应用具有重要参考价值。比重为[X]，反映了土颗粒的密度特性，与土样的矿物成分密切相关。通过对这些物理参数的测定，可以初步了解土样的基本性质，为后续试验提供基础数据。[此处插入表1：土样物理参数测定结果][此处插入表1：土样物理参数测定结果]土样的颗粒级配是影响其力学性质的重要因素之一，因此对土样进行了颗粒级配分析。采用筛分法和沉降法相结合的方式，对土样中不同粒径的颗粒含量进行测定。具体操作过程中，首先使用标准筛对土样进行筛分，将土样按照粒径大小分为不同的粒组。对于粒径小于0.075mm的细颗粒部分，采用沉降法进行分析，根据斯托克斯定律，通过测量颗粒在液体中的沉降速度，计算出颗粒的粒径分布。将两种方法得到的数据进行综合分析，绘制出土样的颗粒级配曲线，如图1所示。[此处插入图1：土样颗粒级配曲线][此处插入图1：土样颗粒级配曲线]从颗粒级配曲线可以看出，土样中粒径在0.075-0.25mm之间的颗粒含量较高，约占总质量的[X]%，这表明该土样主要由粉砂颗粒组成。不均匀系数Cu为[X]，曲率系数Cc为[X]，根据相关标准判断，该土样的级配不良。不均匀系数反映了土样中不同粒径颗粒的分布均匀程度，不均匀系数越大，说明土样中颗粒大小差异越大。曲率系数则反映了土样颗粒级配曲线的形状，曲率系数在1-3之间时，土样的级配较为良好。而本试验土样的不均匀系数和曲率系数表明，其颗粒分布相对集中，缺乏中间粒径的颗粒，级配不良。这种级配特性会对土样的压实性能和力学性质产生重要影响，在后续试验和工程应用中需要予以充分考虑。为了研究土样的压实特性，进行了压实试验。采用重型击实试验方法，按照《土工试验方法标准》（GB/T50123-2019）中的相关规定进行操作。试验过程中，将土样按照不同的含水率制备成多个试样，每个试样在规定的击实功下进行击实。通过测量击实后土样的干密度，绘制出干密度与含水率的关系曲线，即击实曲线，如图2所示。[此处插入图2：土样击实曲线][此处插入图2：土样击实曲线]从击实曲线可以看出，随着含水率的增加，土样的干密度先增大后减小，存在一个最优含水率。在最优含水率为[X]%时，土样的最大干密度达到[X]g/cm³。这是因为在含水率较低时，土颗粒之间的摩擦力较大，土样难以压实。随着含水率的增加，水分在土颗粒之间起到润滑作用，减小了颗粒间的摩擦力，使得土样更容易被压实，干密度逐渐增大。当含水率超过最优含水率后，过多的水分占据了土颗粒之间的孔隙，导致土样的干密度反而减小。了解土样的压实特性，对于工程施工中的填方作业具有重要指导意义，在填方时应控制土样的含水率接近最优含水率，以达到最佳的压实效果，提高土体的强度和稳定性。4.2试验内容与方案4.2.1非饱和土三轴剪切试验非饱和土三轴剪切试验的目的在于深入探究非饱和粉砂在不同吸力和围压条件下的力学特性，揭示其应力-应变关系、强度特性以及体变特性等，为建立非饱和粉砂临界状态理论提供关键的试验数据支持。本试验采用GDS非饱和土三轴试验系统，该系统能够精确控制试验过程中的吸力、围压等参数，确保试验结果的准确性和可靠性。在试验前，需对土样进行精心制备。根据压实试验得到的最优含水率和最大干密度，将土样制备成具有特定干密度和含水率的圆柱形试样，以保证试验结果的一致性和可比性。采用静压法将土样压制成高度为80mm、直径为39.1mm的标准试样。在压制过程中，严格控制压力和压实时间，确保土样的压实均匀性。试验方案设计如下：共设计6组不同的试验工况，每组试验设置3个平行试样，以减小试验误差。具体试验工况见表2。在每组试验中，通过轴平移技术精确控制吸力，利用数字型压力/体积控制器设定目标吸力值，并实时监测孔隙气压力和孔隙水压力，确保吸力控制的精度。围压则通过向三轴压力室内充入液体来施加，同样由数字型压力/体积控制器进行精确控制。在试验过程中，按照设定的加载速率，通过轴向加载杆缓慢施加轴向荷载，同时实时采集偏应力、轴向应变、体应变等数据。[此处插入表2：非饱和土三轴剪切试验工况][此处插入表2：非饱和土三轴剪切试验工况]试验步骤如下：首先，将制备好的土样安装在三轴压力室中，在土样底部放置饱和的高进气值陶土板，土样顶部放置透水石，然后用橡皮膜将土样包裹紧密，确保土样与外界环境隔离。向三轴压力室内充入液体，施加一定的围压，使土样在围压作用下达到稳定状态。利用轴平移技术，通过控制孔隙气压力和孔隙水压力，使土样达到设定的吸力值，并保持稳定一段时间，确保土样内部的水分分布均匀。按照设定的加载速率，通过轴向加载杆缓慢施加轴向荷载，同时实时采集偏应力、轴向应变、体应变等数据。当偏应力达到峰值或轴向应变达到15%-20%时，认为土样达到破坏状态，停止加载。试验结束后，拆除土样，对土样的破坏形态进行详细观察和记录。4.2.2饱和土三轴剪切试验饱和土三轴剪切试验的主要目的是获取饱和粉砂的力学特性参数，为与非饱和粉砂的试验结果进行对比分析提供基础数据。通过对比饱和状态和非饱和状态下粉砂的力学性质差异，能够更深入地理解吸力对非饱和粉砂力学行为的影响机制，进一步完善非饱和粉砂临界状态理论。本试验同样采用GDS非饱和土三轴试验系统，利用其高精度的控制和测量功能，确保试验结果的准确性。试验土样的制备方法与非饱和土三轴剪切试验相同，均制备成高度为80mm、直径为39.1mm的标准圆柱形试样。为使土样达到饱和状态，采用反压饱和法。将土样安装在三轴压力室中，通过向土样内部施加反压力，使土样中的孔隙水压力逐渐增大，将土样中的气体排出，从而使土样达到饱和状态。在饱和过程中，密切监测孔隙水压力和体变的变化，当体变不再发生明显变化且孔隙水压力达到稳定值时，认为土样已达到饱和状态。试验方案设置了3种不同的围压条件，分别为50kPa、100kPa和200kPa，每种围压条件下设置3个平行试样。试验过程中，按照设定的加载速率，通过轴向加载杆缓慢施加轴向荷载，同时实时采集偏应力、轴向应变、体应变等数据。当偏应力达到峰值或轴向应变达到15%-20%时，认为土样达到破坏状态，停止加载。试验结束后，对试验数据进行整理和分析，绘制偏应力-轴向应变曲线、体应变-轴向应变曲线等，获取饱和粉砂的抗剪强度参数、变形特性等信息。饱和土三轴剪切试验与非饱和土三轴剪切试验的对比意义在于，通过对比不同饱和度下粉砂的力学性质，可以清晰地观察到吸力对非饱和粉砂强度和变形特性的影响。在饱和土中，不存在基质吸力，其力学性质主要取决于土颗粒间的摩擦力和黏聚力。而在非饱和土中，基质吸力的存在增加了土颗粒间的有效应力，从而改变了土体的力学性质。通过对比试验，可以定量分析吸力对非饱和粉砂抗剪强度、体变特性等的影响程度，为非饱和粉砂临界状态的研究提供更全面的认识。4.2.3非饱和土土-水特征曲线试验非饱和土土-水特征曲线试验的目的是测定非饱和粉砂在不同吸力下的含水率变化，建立土-水特征曲线，从而深入了解非饱和粉砂的水力特性。土-水特征曲线反映了非饱和土的基质吸力与含水率之间的关系，是研究非饱和土水分运移、渗透性等水力特性的重要依据。通过土-水特征曲线，能够预测非饱和土在不同吸力条件下的含水率分布，为分析非饱和粉砂在实际工程中的水分变化情况提供理论支持。本试验采用压力板仪法来测定非饱和粉砂的土-水特征曲线。压力板仪法是一种常用的测定非饱和土土-水特征曲线的方法，其原理是利用高进气值陶土板的毛细作用，将土样中的水分在不同吸力下逐渐排出，通过测量排出的水分体积来计算土样的含水率。在试验前，将土样制备成高度为20mm、直径为50mm的圆形试样。将土样放置在饱和的高进气值陶土板上，然后将土样和陶土板一起放入压力板仪中。通过向压力板仪中施加不同的气体压力，使土样中的水分在压力差的作用下逐渐排出。在排水过程中，每隔一定时间测量一次排出的水分体积，根据排出的水分体积和土样的初始质量，计算出土样在不同吸力下的含水率。土-水特征曲线对临界状态研究的作用主要体现在以下几个方面：首先，土-水特征曲线能够反映非饱和粉砂的孔隙结构和水分分布情况，而孔隙结构和水分分布是影响非饱和粉砂力学性质的重要因素。在临界状态下，非饱和粉砂的孔隙结构和水分分布会发生变化，通过土-水特征曲线可以了解这些变化对临界状态的影响。其次，土-水特征曲线与非饱和粉砂的渗透性密切相关。在临界状态下，非饱和粉砂的渗透性可能会发生突变，通过土-水特征曲线可以预测渗透性的变化，为分析临界状态下非饱和粉砂的渗流特性提供依据。土-水特征曲线还可以用于建立非饱和粉砂的本构模型。在建立本构模型时，需要考虑非饱和粉砂的水力特性，土-水特征曲线为描述非饱和粉砂的水力特性提供了重要的参数，有助于建立更准确的本构模型，从而更好地描述非饱和粉砂在临界状态下的力学行为。4.3土样制备土样制备是试验的关键环节，其质量直接影响试验结果的准确性和可靠性。为确保试验用土样的质量，采取了一系列严格的制备步骤和质量控制措施。首先，对采集的土样进行风干处理。将土样置于通风良好、阳光充足的地方，自然风干至含水率较低且均匀的状态。在风干过程中，定期翻动土样，使其各部分水分均匀散失，避免出现局部干湿不均的情况。风干后的土样便于后续的粉碎和过筛操作，同时也能减少土样中水分对试验结果的影响。接着，对风干后的土样进行粉碎和过筛。使用粉碎设备将土样粉碎成细小颗粒，以便通过筛网。采用标准筛对土样进行过筛，根据试验要求，选择合适孔径的筛网。本试验选用孔径为2mm的筛网，将土样过筛后，去除粒径大于2mm的颗粒，确保土样的颗粒均匀性。过筛后的土样颗粒级配更加均匀，有利于保证试验结果的一致性。然后，根据压实试验得到的最优含水率和最大干密度，进行配样。计算所需土样和水的质量，将过筛后的土样放入容器中，按照计算好的水量加入适量的水，充分搅拌均匀，使土样达到设定的含水率。在搅拌过程中，确保水与土样充分混合，避免出现水分聚集或土样干湿不均的现象。配样完成后，将土样装入密封袋中，放置一段时间，让土样中的水分充分均匀分布，一般放置时间不少于24小时。在进行三轴试验前，采用静压法将配好的土样制备成标准试样。使用专门的制样模具，将土样分层填入模具中，每层土样都进行均匀压实。在压实过程中，严格控制压实压力和压实次数，确保土样的压实度均匀一致。采用压力机对土样进行静压，压力控制在[X]MPa，压实次数为[X]次。将土样压制成高度为80mm、直径为39.1mm的圆柱形试样，以满足三轴试验的要求。土样制备过程中的质量控制措施至关重要。在风干过程中，密切关注土样的含水率变化，确保风干后的土样含水率符合要求。在粉碎和过筛过程中，检查筛网的孔径是否符合标准，避免出现筛网破损或孔径不准确的情况，影响土样的颗粒级配。在配样过程中，严格按照计算好的土样和水的质量进行配制，使用高精度的天平进行称量，确保配样的准确性。在制样过程中，定期检查制样模具的尺寸和形状，保证制样的精度。对制备好的试样进行外观检查，确保试样表面光滑、无裂缝、无气泡等缺陷。通过以上质量控制措施，保证了土样制备的质量，为试验结果的准确性和可靠性提供了有力保障。五、非饱和粉砂临界状态试验结果与分析5.1饱和粉砂的临界状态对饱和粉砂进行三轴剪切试验，得到不同围压下的偏应力-轴向应变关系曲线、体应变-轴向应变关系曲线等，试验结果如表3和图3所示。从偏应力-轴向应变关系曲线可以看出，随着轴向应变的增加，偏应力逐渐增大，当轴向应变达到一定值后，偏应力达到峰值，随后偏应力略有下降并趋于稳定。围压越大，偏应力峰值越高，表明饱和粉砂的抗剪强度随着围压的增大而增大。在体应变-轴向应变关系方面，随着轴向应变的增加，体应变先减小后增大，呈现出剪缩-剪胀的特性。围压较小时，剪胀现象更为明显，这是因为围压较小时，土颗粒之间的约束力较小，在剪切过程中更容易发生错动和膨胀。[此处插入表3：饱和粉砂三轴剪切试验结果][此处插入图3：饱和粉砂偏应力-轴向应变、体应变-轴向应变关系曲线][此处插入表3：饱和粉砂三轴剪切试验结果][此处插入图3：饱和粉砂偏应力-轴向应变、体应变-轴向应变关系曲线][此处插入图3：饱和粉砂偏应力-轴向应变、体应变-轴向应变关系曲线]在临界状态下，饱和粉砂的应力比M（M=q/p'，其中q为偏应力，p'为平均有效应力）和孔隙比e存在一定的关系。通过对试验数据的分析，得到饱和粉砂在临界状态下的应力比M约为[X]，孔隙比e与平均有效应力p'的关系可表示为e=[a]-[b]\lnp'，其中[a]和[b]为试验拟合参数。饱和砂土的临界状态线（CSL）可以在e-\lnp'平面和q-p'平面内进行描述。在e-\lnp'平面内，临界状态线方程为e_{cs}=\lambda_0-\kappa\lnp'，其中\lambda_0为正常固结线在e-\lnp'平面上的截距，\kappa为回弹指数。通过对试验数据的拟合，得到\lambda_0=[X]，\kappa=[X]。在q-p'平面内，临界状态线方程为q=Mp'，其中M为应力比，通过试验得到M=[X]。将本试验得到的临界状态线方程与已有研究成果进行对比，发现本试验结果与[已有研究文献]中关于饱和砂土临界状态线的研究结果基本一致，验证了本试验结果的可靠性。在[已有研究文献]中，通过对不同地区饱和砂土的试验研究，得到的临界状态线方程在形式上与本试验结果相似，且应力比M和相关参数的值也较为接近。这表明不同地区的饱和砂土在临界状态下的力学行为具有一定的共性，本试验结果具有一定的代表性。5.2非饱和粉砂土-水特征曲线进行非饱和土土-水特征曲线试验，采用压力板仪法测定非饱和粉砂在不同吸力下的含水率变化。试验结果如表4和图4所示。从试验结果可以看出，随着吸力的增加，非饱和粉砂的含水率逐渐降低。当吸力较小时，含水率下降较为缓慢；当吸力超过一定值后，含水率下降速度加快。这是因为在吸力较小时，土中孔隙水主要以结合水的形式存在，结合水与土颗粒之间的作用力较强，不易被排出。随着吸力的增大，孔隙水逐渐被排出，土中孔隙水主要以自由水的形式存在，自由水与土颗粒之间的作用力较弱，容易被排出，导致含水率下降速度加快。[此处插入表4：非饱和粉砂土-水特征曲线试验结果][此处插入图4：非饱和粉砂土-水特征曲线][此处插入表4：非饱和粉砂土-水特征曲线试验结果][此处插入图4：非饱和粉砂土-水特征曲线][此处插入图4：非饱和粉砂土-水特征曲线]为了更准确地描述非饱和粉砂的土-水特征曲线，采用VanGenuchten模型对试验结果进行拟合。VanGenuchten模型的表达式为：\theta=\theta_r+\frac{\theta_s-\theta_r}{(1+(\alphas)^n)^{1-\frac{1}{n}}}其中，\theta为体积含水率，\theta_r为残余含水率，\theta_s为饱和含水率，\alpha、n为模型参数，s为吸力。通过对试验数据的拟合，得到VanGenuchten模型的参数\alpha=[X]，n=[X]，\theta_r=[X]，\theta_s=[X]。将拟合曲线与试验数据进行对比，如图5所示。可以看出，VanGenuchten模型能够较好地拟合非饱和粉砂的土-水特征曲线，拟合曲线与试验数据点吻合较好，说明该模型能够准确地描述非饱和粉砂的基质吸力与含水率之间的关系。[此处插入图5：VanGenuchten模型拟合曲线与试验数据对比][此处插入图5：VanGenuchten模型拟合曲线与试验数据对比]非饱和粉砂土-水特征曲线的参数受到多种因素的影响。土样的颗粒级配是影响土-水特征曲线的重要因素之一。颗粒级配良好的土样，其孔隙分布较为均匀，土-水特征曲线相对平缓；而颗粒级配不良的土样，其孔隙分布不均匀，土-水特征曲线变化较为剧烈。本试验所用土样的颗粒级配不良，导致其土-水特征曲线在吸力较大时，含水率下降速度较快。土样的矿物成分也会对土-水特征曲线产生影响。不同矿物成分的土颗粒对水的吸附能力不同，从而影响土-水特征曲线的形状和参数。土样的初始含水率同样会影响土-水特征曲线。初始含水率较高的土样，在相同吸力作用下，含水率下降幅度较大。在实际工程中，需要充分考虑这些因素对土-水特征曲线的影响，以便更准确地预测非饱和粉砂的水力特性。5.3吸力平衡过程分析在吸力平衡试验中，对孔隙水和试样体积变化进行了详细监测。试验结果如表5所示。从孔隙水体积变化来看，随着吸力平衡时间的增加，孔隙水体积逐渐减小。这是因为在吸力作用下，土样中的孔隙水被逐渐排出，导致孔隙水体积减小。在吸力为50kPa时，初始孔隙水体积为[X]cm³，经过12小时的吸力平衡，孔隙水体积减小至[X]cm³。试样体积也随着吸力平衡时间的增加而逐渐减小，这是由于孔隙水的排出导致土颗粒间的孔隙减小，从而使试样体积收缩。在吸力为50kPa时，初始试样体积为[X]cm³，12小时后试样体积减小至[X]cm³。[此处插入表5：孔隙水和试样体积变化试验结果][此处插入表5：孔隙水和试样体积变化试验结果]吸力平衡前后比体积、孔隙水比体积和饱和度的演化规律对于理解非饱和粉砂的力学性质具有重要意义。比体积v（v=1+e，其中e为孔隙比）随着吸力的增加而减小。在吸力平衡前，比体积为[X]，当吸力达到50kPa并平衡后，比体积减小至[X]。这是因为吸力的增加使得土颗粒间的有效应力增大，土颗粒相互靠拢，孔隙比减小，从而导致比体积减小。孔隙水比体积v_w（v_w=\frac{V_w}{V_s}，其中V_w为孔隙水体积，V_s为土颗粒体积）也随着吸力的增加而减小。在吸力平衡前，孔隙水比体积为[X]，当吸力达到50kPa并平衡后，孔隙水比体积减小至[X]。这是由于吸力的作用使孔隙水体积减小，而土颗粒体积不变，导致孔隙水比体积减小。饱和度S_r（S_r=\frac{V_w}{V_v}，其中V_v为孔隙体积）随着吸力的增加而减小。在吸力平衡前，饱和度为[X]%，当吸力达到50kPa并平衡后，饱和度减小至[X]%。这是因为吸力增加使孔隙水体积减小，而孔隙体积不变或略有减小，从而导致饱和度降低。将本试验结果与已有研究成果进行对比，发现本试验中比体积、孔隙水比体积和饱和度随吸力的变化趋势与[已有研究文献]中的结果基本一致。在[已有研究文献]中，通过对不同类型非饱和土的研究，也得出了比体积、孔隙水比体积和饱和度随吸力增加而减小的结论。这进一步验证了本试验结果的可靠性，同时也表明不同类型的非饱和土在吸力平衡过程中，比体积、孔隙水比体积和饱和度的演化规律具有一定的共性。5.4等吸力固结特性在等吸力固结过程中，对非饱和粉砂的比体积、孔隙水比体积和饱和度变化进行了详细监测，试验结果如表6所示。从比体积变化来看，随着等吸力固结压力的增加，比体积逐渐减小。在吸力为50kPa，固结压力从50kPa增加到200kPa时，比体积从[X]减小至[X]。这是因为在等吸力固结过程中，土颗粒在压力作用下逐渐靠拢，孔隙比减小，从而导致比体积减小。[此处插入表6：等吸力固结过程中比体积、孔隙水比体积和饱和度变化试验结果][此处插入表6：等吸力固结过程中比体积、孔隙水比体积和饱和度变化试验结果]孔隙水比体积也随着等吸力固结压力的增加而减小。在吸力为50kPa，固结压力从50kPa增加到200kPa时，孔隙水比体积从[X]减小至[X]。这是由于在固结压力作用下，土样中的孔隙水被挤出，孔隙水体积减小，而土颗粒体积不变，导致孔隙水比体积减小。饱和度同样随着等吸力固结压力的增加而减小。在吸力为50kPa，固结压力从50kPa增加到200kPa时，饱和度从[X]%减小至[X]%。这是因为随着固结压力的增加，孔隙水体积减小，而孔隙体积不变或略有减小，从而导致饱和度降低。等吸力固结前后比体积、孔隙水比体积和饱和度的演化规律表明，等吸力固结过程对非饱和粉砂的物理性质产生了显著影响。在等吸力固结前，非饱和粉砂的比体积、孔隙水比体积和饱和度相对较大。