基于DIC方法洞察单轴压缩湿砂样主应变轴偏转规律的实验剖析

基于DIC方法洞察单轴压缩湿砂样主应变轴偏转规律的实验剖析一、绪论1.1研究背景与意义在各类岩土工程实践中，砂土作为一种常见的地质体材料，其力学性质对工程的稳定性与安全性起着关键作用。无论是道路与桥梁建设中的地基处理，还是高层建筑基础的设计与施工，砂土的力学特性都直接关系到工程结构的承载能力与长期稳定性。在水利工程中，堤坝、水闸等设施的基础常建于砂土地基之上，砂土的抗剪强度、压缩性等力学指标决定了这些水利设施能否承受水压力、渗透力等复杂荷载作用，确保水利工程的安全运行。而在地下工程，如地铁隧道、地下停车场等建设中，砂土的力学性质影响着隧道的开挖稳定性、支护结构的设计以及地面沉降控制等关键问题。因此，深入研究砂土的力学性质，对于保障各类工程的安全、稳定与可持续发展具有至关重要的意义。主应变轴偏转规律是砂土力学性质研究中的一个重要方面，它反映了砂土在复杂应力状态下的变形特性。在实际工程中，砂土往往受到多种复杂应力的作用，主应力方向不断变化，导致主应变轴发生偏转。车辆行驶过程中对道路地基的反复荷载作用，海浪对海岸堤坝地基的冲击，地震时地震波引起的土体振动等，都会使砂土的主应力方向发生改变，进而引发主应变轴的偏转。这种主应变轴的偏转对砂土的力学行为产生显著影响，可能导致砂土的强度降低、变形增大，甚至引发工程结构的破坏。在地震作用下，砂土主应变轴的偏转会使其抗剪强度下降，增加地基液化的风险，对建筑物的安全构成严重威胁。从理论研究的角度来看，深入探究主应变轴偏转规律有助于完善砂土的本构模型。现有的砂土本构模型大多基于简单的应力-应变关系，难以准确描述砂土在复杂应力状态下的力学行为。通过研究主应变轴偏转规律，可以更全面地了解砂土在不同应力路径下的变形机制，为建立更精确、更符合实际的砂土本构模型提供理论依据。这不仅能够丰富土力学的理论体系，推动土力学学科的发展，还能为岩土工程的数值模拟分析提供更可靠的工具，提高工程设计的准确性和科学性。在实际工程应用方面，掌握主应变轴偏转规律对岩土工程的设计与施工具有重要的指导意义。在地基处理设计中，考虑主应变轴偏转的影响，可以更合理地选择地基处理方法和参数，提高地基的承载能力和稳定性。在深基坑支护设计中，准确把握砂土在复杂应力下的主应变轴偏转规律，有助于优化支护结构的设计，确保基坑开挖过程中土体的稳定，减少周边环境的影响。对于地震区的工程建设，了解主应变轴偏转对砂土液化特性的影响，能够制定更有效的抗震措施，降低地震灾害对工程结构的破坏风险，保障人民生命财产安全。1.2研究现状1.2.1主应力轴偏转条件下土体变形特性研究概述主应力轴偏转条件下土体变形特性的研究一直是土力学领域的关键内容。早期研究主要集中在简单应力路径下的土体变形分析，随着对工程实际中复杂应力状态认识的加深，主应力轴偏转对土体变形的影响逐渐受到关注。学者们通过室内试验、理论分析和数值模拟等多种手段展开研究，取得了一系列成果。在室内试验方面，空心圆柱扭剪仪等先进设备的应用，使得精确模拟主应力轴偏转条件下的土体受力状态成为可能。通过这类试验，研究人员发现主应力轴偏转时，土体的应力-应变关系呈现出与传统加载方式不同的特征。土体的剪切强度、压缩模量和剪切模量等力学参数会随着主应力轴的偏转发生变化，且这种变化往往具有非线性和各向异性的特点。在主应力轴连续旋转过程中，砂土的应变量和固结指数不断变化，塑性变形显著增加，其剪切强度随主应力轴旋转角度的增大而增大，且增长速率逐渐加快，表现出明显的非线性特征。在理论分析层面，一些学者基于连续介质力学、塑性力学等理论，建立了考虑主应力轴偏转影响的本构模型。这些模型试图从理论上描述土体在复杂应力状态下的变形机制，为工程计算提供理论依据。然而，由于土体材料的复杂性和多样性，目前的本构模型仍存在一定的局限性，难以全面准确地反映主应力轴偏转条件下土体的力学行为。数值模拟方法也被广泛应用于该领域的研究。有限元、离散元等数值方法能够对土体在主应力轴偏转下的变形过程进行数值模拟，直观地展示土体内部的应力、应变分布情况。通过数值模拟，可以深入研究不同因素对土体变形特性的影响，如土体的初始状态、加载速率、主应力轴偏转幅度和频率等。数值模拟结果与试验数据的对比验证，也有助于改进和完善理论模型。然而，当前的研究仍存在一些不足之处。部分研究主要针对特定类型的土体或特定的试验条件，所得结论的普适性有待进一步验证。对于土体在复杂应力路径下的变形机制，尤其是主应力轴快速、大幅度偏转时的变形机理，尚未完全明确。不同研究方法之间的结果存在一定差异，缺乏统一的理论框架来协调和整合这些研究成果。1.2.2土体非共轴特性研究概述土体的非共轴特性是指主应变方向与主应力方向不一致的现象，这一特性在土体力学行为中起着重要作用，近年来受到了广泛关注。研究表明，土体的非共轴特性普遍存在，尤其在主应力轴发生偏转时更为明显。在实验研究方面，众多学者通过各种先进的试验技术，如数字图像相关技术（DIC）、扫描电子显微镜（SEM）等，对土体的非共轴特性进行了深入探究。利用DIC技术可以精确测量土体表面的位移和应变场，从而获取主应变方向的信息；SEM则能够观察土体微观结构的变化，揭示非共轴特性的微观机制。通过这些试验手段，发现土体在受力过程中，其内部颗粒的排列和相互作用会发生改变，导致主应变方向与主应力方向出现偏差。在定向加载的应力状态下，随着应力的加载，非共轴现象会逐渐减弱，但在主应力轴纯旋转的过程中，非共轴现象基本上一直存在。从理论研究角度，为了解释土体的非共轴特性，学者们提出了多种理论和模型。一些基于微观力学的理论，从土体颗粒的接触力学和相互作用出发，建立了描述非共轴特性的微观模型；另一些则在宏观本构模型中引入非共轴参数，以考虑主应变方向与主应力方向的不一致性。这些理论和模型在一定程度上能够解释土体的非共轴现象，但仍存在一些问题和挑战。微观模型虽然能够反映土体的微观结构特征，但计算过程较为复杂，难以应用于实际工程；宏观本构模型中的非共轴参数确定较为困难，且模型的适用性受到一定限制。土体的非共轴特性对主应变轴偏转有着显著影响。由于非共轴特性的存在，主应变轴在土体变形过程中的偏转规律变得更加复杂。这种复杂性不仅增加了对土体力学行为理解的难度，也给岩土工程的设计和分析带来了挑战。在进行地基承载力计算、边坡稳定性分析等工程问题时，如果忽略土体的非共轴特性，可能会导致计算结果与实际情况存在较大偏差，从而影响工程的安全性和可靠性。1.2.3主轴偏转与剪切带关系研究概述主轴偏转与剪切带之间的关系是土力学研究中的一个重要课题，对于理解土体的破坏机制和工程稳定性具有关键意义。剪切带是土体在受力过程中出现的局部化变形区域，其形成与发展往往伴随着主应力轴和主应变轴的偏转。在早期的研究中，学者们主要通过宏观试验观察和理论分析来探讨主轴偏转与剪切带的关系。通过对土体进行单轴压缩、三轴压缩等试验，发现当土体达到一定的应力水平时，会出现明显的剪切带，同时主应力轴和主应变轴也会发生相应的偏转。一些理论分析认为，剪切带的形成是由于土体内部应力分布不均匀，导致局部区域的剪应力集中，从而引发材料的屈服和塑性流动，进而促使主应力轴和主应变轴发生偏转。随着研究的深入，数值模拟方法在该领域得到了广泛应用。有限元、离散元等数值方法能够模拟土体在复杂应力条件下的变形过程，直观地展示剪切带的形成与发展以及主轴的偏转情况。通过数值模拟，可以分析不同因素对主轴偏转和剪切带关系的影响，如土体的初始密度、含水率、加载速率等。研究发现，土体的初始密度对剪切带的形态和主轴偏转角度有显著影响，初始密度较大的土体，其剪切带更窄，主轴偏转角度相对较小；加载速率的变化也会影响剪切带的发展和主轴的偏转，加载速率较快时，剪切带的形成和发展更为迅速，主轴偏转也更为明显。尽管目前在主轴偏转与剪切带关系的研究方面取得了一定的成果，但仍存在一些问题需要进一步解决。对于剪切带内部的细观力学机制，如颗粒的运动和相互作用、孔隙水压力的变化等，尚未完全明确。不同研究方法得到的结果之间存在一定的差异，缺乏统一的认识和解释。在实际工程应用中，如何准确地考虑主轴偏转与剪切带关系对土体力学行为的影响，仍然是一个亟待解决的问题。1.2.4数字图像相关方法的应用概述数字图像相关（DIC）方法作为一种非接触式的光学测量技术，自上世纪80年代后期发展以来，在力学研究领域得到了广泛应用，并逐渐在岩土工程领域展现出独特的优势和重要作用。在岩土工程领域，DIC方法的应用主要体现在土体变形测量方面。传统的土体变形测量方法，如应变片测量、位移计测量等，往往存在接触式测量的局限性，如对土体表面造成扰动、测量点有限、无法获取全场变形信息等。而DIC方法通过对土体表面散斑图像的采集和分析，能够实现对土体表面位移和应变的全场、非接触式测量。在土体的三轴压缩试验、直剪试验、单轴压缩试验等各类力学试验中，DIC方法可以精确测量土体在加载过程中的变形情况，获取土体表面的位移场和应变场信息，包括主应变的大小和方向，从而为研究土体的力学行为提供丰富的数据支持。DIC方法在研究主应变轴偏转规律方面具有重要价值。通过DIC技术测量土体在复杂应力状态下的主应变方向变化，可以直观地观察到主应变轴的偏转过程，分析其偏转规律与应力状态、土体性质等因素之间的关系。在主应力轴偏转条件下的砂土试验中，利用DIC方法能够实时监测砂土表面的变形，准确确定主应变轴的偏转角，为深入研究砂土在复杂应力下的变形特性提供了有效的手段。除了土体变形测量，DIC方法还在岩土工程的其他方面得到应用。在岩石力学试验中，DIC方法可以用于分析岩石的破坏过程，研究岩石裂纹的扩展规律和强度特性；在边坡稳定性研究中，DIC方法能够对边坡表面的变形进行监测，及时发现潜在的滑坡隐患。DIC方法还可以与其他技术相结合，如CT扫描技术、数值模拟技术等，实现对岩土体的多尺度、多维度研究，进一步深化对岩土工程问题的认识。尽管DIC方法在岩土工程领域取得了广泛应用，但在实际应用中仍面临一些挑战。图像采集过程中可能受到环境因素（如光照变化、灰尘干扰等）的影响，导致图像质量下降，从而影响测量精度；对于复杂形状和变形较大的岩土体，图像匹配和分析的难度较大，需要进一步改进算法和技术；DIC方法目前主要适用于表面变形测量，对于岩土体内部的变形信息获取还存在一定的局限性。1.3研究不足与问题尽管目前在砂土主应变轴偏转规律及相关领域已取得了一定的研究成果，但仍存在诸多不足之处，这些问题限制了对砂土力学行为的深入理解和工程应用的准确性。在主应变轴偏转规律研究方面，现有研究对砂土在复杂应力路径下主应变轴偏转的完整过程和内在机制的认识还不够全面和深入。部分研究仅关注了特定加载条件下的主应变轴偏转情况，对于多因素耦合作用，如不同加载速率、复杂应力组合以及土体初始状态差异等对主应变轴偏转规律的影响，缺乏系统的研究。在实际工程中，砂土往往受到多种复杂因素的共同作用，这些因素的相互影响可能导致主应变轴偏转规律与单一因素作用下的情况存在显著差异。在地震作用下，砂土不仅承受着地震波引起的动荷载，还可能受到地下水渗流、土体自身结构性变化等因素的影响，这些因素如何共同作用于主应变轴的偏转，目前尚未有清晰的定论。关于土体非共轴特性对主应变轴偏转的影响，虽然已认识到两者之间存在密切联系，但在具体的量化关系和作用机制方面，研究还相对薄弱。目前用于描述非共轴特性的参数和模型在实际应用中存在一定的局限性，难以准确地反映土体在复杂应力状态下非共轴特性的变化规律，以及这种变化对主应变轴偏转的具体影响程度。不同类型的土体其非共轴特性和主应变轴偏转之间的关系可能存在差异，而现有的研究未能充分考虑土体的多样性，所得结论的普适性有待进一步提高。在主轴偏转与剪切带关系的研究中，虽然对两者的宏观关联有了一定的了解，但对于剪切带内部细观结构变化如何影响主应变轴偏转，以及主应变轴偏转又如何反馈作用于剪切带的发展和演化，相关研究还不够细致和深入。剪切带内部的颗粒破碎、重排列以及孔隙结构变化等细观过程对主应变轴偏转的影响机制尚不完全清楚，这使得在从细观层面解释砂土的力学行为和破坏机制时存在一定的困难。在实际工程中，准确把握主轴偏转与剪切带关系对于预测砂土的破坏模式和工程稳定性至关重要，而目前的研究成果在这方面的应用还存在一定的差距。数字图像相关（DIC）方法在岩土工程中的应用虽然为研究主应变轴偏转提供了新的手段，但在实际应用中仍面临一些挑战。图像采集过程易受环境因素干扰，如光照不均匀、外界振动等，可能导致采集的图像质量下降，从而影响DIC方法测量的准确性和可靠性。对于大变形和复杂形状的砂土试样，图像匹配和应变计算的难度较大，容易出现误差，如何提高DIC方法在这些情况下的测量精度，是需要进一步解决的问题。此外，DIC方法目前主要侧重于表面变形测量，对于砂土内部的应变分布和主应变轴偏转情况，难以直接获取，如何实现DIC方法与其他技术的有效结合，以获取砂土内部的变形信息，也是未来研究的一个重要方向。综上所述，目前砂土主应变轴偏转规律的研究仍存在诸多需要完善和深入的地方。开展基于DIC方法的单轴压缩湿砂样主应变轴偏转规律实验研究，有望在一定程度上弥补现有研究的不足，深入揭示砂土在复杂应力条件下的变形特性，为砂土力学理论的发展和岩土工程的实际应用提供更坚实的基础。1.4研究内容与方法1.4.1研究内容本研究聚焦于基于DIC方法的单轴压缩湿砂样主应变轴偏转规律，主要研究内容涵盖以下几个关键方面：基于DIC方法的单轴压缩湿砂样试验系统搭建：精心挑选适宜的试验设备，如高精度的压力试验机，以确保能够稳定、精确地施加单轴压缩荷载，满足试验过程中对荷载控制的严格要求。同时，合理配置DIC测量系统，包括高分辨率相机、合适的镜头以及专业的图像采集与分析软件，保证能够清晰、准确地采集湿砂样表面的散斑图像，并对其进行高效、精确的处理和分析。此外，还需对试验系统进行全面、细致的标定和调试，严格控制试验环境条件，如温度、湿度等，以确保试验数据的准确性和可靠性，为后续的试验研究奠定坚实的基础。不同工况下单轴压缩湿砂样主应变轴偏转规律分析：系统地开展不同含水率、干密度和加载速率等工况下的单轴压缩湿砂样试验。在试验过程中，运用DIC技术对湿砂样表面的变形进行实时、全面的监测，精确获取主应变轴的偏转角随时间、荷载等因素的变化数据。深入分析不同工况对主应变轴偏转规律的影响，探究含水率的变化如何改变砂土颗粒间的润滑作用和孔隙水压力分布，进而影响主应变轴的偏转；研究干密度的差异怎样导致砂土颗粒排列紧密程度的不同，从而对主应变轴偏转产生作用；分析加载速率的快慢如何影响砂土的变形速率和能量耗散，最终影响主应变轴的偏转特性。考虑非共轴特性的主应变轴偏转机制探讨：深入研究土体非共轴特性对主应变轴偏转的影响机制，通过对试验数据的细致分析和理论推导，揭示主应变方向与主应力方向不一致的内在原因和变化规律。尝试引入合适的非共轴参数，建立考虑非共轴特性的主应变轴偏转理论模型，通过与试验结果的对比验证，不断优化和完善模型，使其能够更准确地描述和预测主应变轴在复杂应力状态下的偏转行为。主应变轴偏转与剪切带形成关系研究：密切关注单轴压缩过程中湿砂样剪切带的形成与发展过程，利用DIC技术精确测量剪切带区域的主应变轴偏转情况，深入分析主应变轴偏转与剪切带形成之间的相互关系。探究主应变轴的偏转如何促使剪切带的萌生和扩展，以及剪切带的形成又如何反过来影响主应变轴的进一步偏转，从细观层面揭示砂土在单轴压缩下的变形破坏机制。主应变轴偏转规律的工程应用分析：将研究得到的主应变轴偏转规律与实际岩土工程案例相结合，如地基沉降分析、边坡稳定性评估等。通过数值模拟和工程实例分析，评估主应变轴偏转对岩土工程稳定性和变形的影响，为岩土工程的设计和施工提供科学、合理的建议和指导，提高工程的安全性和可靠性。1.4.2研究方法本研究综合运用多种研究方法，相互补充、相互验证，以确保研究的全面性、深入性和准确性：室内试验方法：采用单轴压缩试验作为主要的试验手段，对不同工况下的湿砂样进行加载试验。在试验过程中，运用DIC技术实时监测湿砂样表面的变形情况，获取主应变轴的偏转角等关键数据。同时，结合传统的测量方法，如位移计测量、压力传感器测量等，对试验数据进行对比和验证，提高数据的可靠性。理论分析方法：基于土力学、弹性力学、塑性力学等相关理论，对主应变轴偏转的机制进行深入分析。建立考虑非共轴特性的理论模型，推导主应变轴偏转角的计算公式，从理论层面解释主应变轴在不同应力状态下的偏转规律。数值模拟方法：运用有限元软件，如ABAQUS、ANSYS等，建立单轴压缩湿砂样的数值模型。通过模拟不同工况下的试验过程，分析砂土内部的应力、应变分布情况，以及主应变轴的偏转过程。将数值模拟结果与试验数据进行对比，验证数值模型的准确性，进一步深入研究主应变轴偏转的影响因素和变化规律。数据分析方法：对试验数据和数值模拟结果进行全面、系统的分析，运用统计学方法、曲线拟合方法等，揭示主应变轴偏转角与各影响因素之间的定量关系。通过绘制图表、曲线等方式，直观地展示主应变轴偏转规律，为研究结论的得出提供有力的数据支持。二、数字图像相关方法（DIC）2.1DIC方法基本原理数字图像相关（DIC）方法是一种基于计算机视觉和图像处理技术的非接触式光学测量方法，在材料力学、生物力学、岩土工程等众多领域中得到了广泛应用。其核心原理是通过跟踪物体表面变形前后散斑图像中同一像素点的位置变化，来获取该像素点的位移向量，进而计算出物体表面的全场位移和应变分布。DIC方法的实现依赖于物体表面随机分布的散斑图案，这些散斑图案如同物体表面的“指纹”，具有唯一性和稳定性。在实验开始前，首先使用高分辨率相机对物体表面进行拍摄，获取初始状态下的散斑图像，该图像被定义为参考图像。当物体受到外力作用发生变形后，再次使用相机对物体表面进行拍摄，得到变形后的散斑图像，即变形图像。在获取参考图像和变形图像后，DIC方法通过一系列复杂的图像处理和计算步骤来实现位移和应变的测量。将参考图像划分为众多互不重叠的小子区域，每个子区域都被视为一个独立的分析单元。对于每个子区域，在变形图像中通过特定的相关算法进行搜索，寻找与该子区域相关性最高的区域，这个过程类似于在茫茫图像中寻找与目标子区域最为匹配的“孪生区域”。一旦找到匹配区域，就可以确定该子区域在变形前后的位置变化，从而得到该子区域中心点的位移向量。通过对所有子区域的位移向量进行计算和分析，就能够构建出物体表面的全场位移场。在计算应变时，通常采用差分法。对于二维平面问题，假设已知某点在x和y方向的位移分量u和v，则该点在x方向的正应变\varepsilon_{xx}可通过对u关于x求偏导数得到，即\varepsilon_{xx}=\frac{\partialu}{\partialx}；在y方向的正应变\varepsilon_{yy}为\varepsilon_{yy}=\frac{\partialv}{\partialy}；切应变\gamma_{xy}则通过\gamma_{xy}=\frac{\partialu}{\partialy}+\frac{\partialv}{\partialx}计算得出。通过这些公式，将位移场转化为应变场，从而全面了解物体表面的变形情况。以简单的拉伸试验为例，假设在试件表面制作了散斑图案，在拉伸过程中，试件会发生伸长变形。通过DIC系统拍摄拉伸前后的散斑图像，经过处理后，能够清晰地看到散斑图案的移动和变形。通过计算，可以得到试件表面各点在拉伸方向上的位移逐渐增大，垂直于拉伸方向的位移基本不变，进而计算出拉伸方向的正应变和切应变，直观地展示出试件在拉伸作用下的变形特性。在岩土工程的土体单轴压缩试验中，利用DIC方法可以测量土体表面在压缩过程中的位移和应变，观察到土体在加载过程中不同部位的变形差异，为研究土体的力学行为提供了丰富的数据支持。2.2相关运算公式在数字图像相关（DIC）方法中，位移和应变的计算涉及一系列严谨且关键的数学公式，这些公式是准确获取物体表面变形信息的核心工具。在二维平面问题中，假设在参考图像中某点的坐标为(x,y)，变形后该点移动到(x+u,y+v)，其中u和v分别为该点在x方向和y方向的位移分量。对于图像中的一个子区域，通过DIC算法计算得到该子区域中心点在x方向和y方向的位移，即为该点的u和v值。应变的计算基于位移的导数。正应变是指物体在某一方向上的相对伸长或缩短，对于二维问题，x方向的正应变\varepsilon_{xx}可表示为：\varepsilon_{xx}=\frac{\partialu}{\partialx}它反映了在x方向上单位长度的变形量。同理，y方向的正应变\varepsilon_{yy}为：\varepsilon_{yy}=\frac{\partialv}{\partialy}切应变则描述了物体在剪切作用下角度的变化。在二维平面中，切应变\gamma_{xy}的计算公式为：\gamma_{xy}=\frac{\partialu}{\partialy}+\frac{\partialv}{\partialx}它综合考虑了x方向位移在y方向的变化以及y方向位移在x方向的变化，体现了物体在平面内的剪切变形程度。在实际计算中，由于图像是离散的像素点构成，导数的计算通常采用差分法进行近似求解。对于x方向的正应变\varepsilon_{xx}，假设相邻像素点在x方向的间距为\Deltax，则可通过以下差分公式近似计算：\varepsilon_{xx}\approx\frac{u(x+\Deltax,y)-u(x,y)}{\Deltax}同理，对于y方向的正应变\varepsilon_{yy}，假设相邻像素点在y方向的间距为\Deltay，其近似计算公式为：\varepsilon_{yy}\approx\frac{v(x,y+\Deltay)-v(x,y)}{\Deltay}对于切应变\gamma_{xy}，其差分近似计算式为：\gamma_{xy}\approx\frac{u(x,y+\Deltay)-u(x,y)}{\Deltay}+\frac{v(x+\Deltax,y)-v(x,y)}{\Deltax}这些差分公式在DIC方法中起着至关重要的作用，通过对离散像素点位移的计算和处理，能够较为准确地估算出物体表面各点的应变值，从而实现对物体变形的全面分析。在对金属薄板进行拉伸实验时，利用DIC技术结合上述公式，可以精确计算出薄板表面在拉伸过程中不同位置的正应变和切应变，清晰地展示出薄板的变形分布情况。在岩土工程的土体单轴压缩实验中，通过DIC测量系统获取土体表面散斑图像的位移信息，运用这些公式计算应变，能够深入研究土体在压缩过程中的变形特性，为土力学理论研究和工程应用提供有力的数据支持。2.3相关搜索方法在数字图像相关（DIC）方法中，准确地在变形后图像中寻找子区域的对应位置是获取位移信息的关键步骤，这依赖于高效且精确的搜索算法。目前，常用的搜索算法主要包括基于灰度的搜索算法和基于特征的搜索算法，它们各自具有独特的原理和应用场景。基于灰度的搜索算法是DIC方法中最为基础和常用的一类算法，其核心思想是利用子区域内像素的灰度信息来度量图像之间的相似性。这类算法假设在物体变形过程中，子区域内的灰度分布保持相对稳定，通过计算参考图像子区域与变形后图像中不同位置子区域的灰度相关性，寻找相关性最高的区域作为对应位置。归一化互相关（NCC）算法是基于灰度搜索算法的典型代表。该算法通过计算参考子区域与候选子区域的归一化互相关系数来衡量两者的相似程度。假设参考子区域的灰度矩阵为I(x,y)，候选子区域的灰度矩阵为J(x,y)，其中(x,y)表示像素点的坐标，归一化互相关系数C的计算公式为：C=\frac{\sum_{x,y}(I(x,y)-\overline{I})(J(x,y)-\overline{J})}{\sqrt{\sum_{x,y}(I(x,y)-\overline{I})^2\sum_{x,y}(J(x,y)-\overline{J})^2}}其中，\overline{I}和\overline{J}分别为参考子区域和候选子区域的平均灰度值。在搜索过程中，遍历变形后图像的各个位置，计算每个位置的候选子区域与参考子区域的归一化互相关系数，系数最大值所对应的位置即为参考子区域在变形后图像中的对应位置。NCC算法对光照变化具有一定的鲁棒性，在图像灰度变化较为均匀的情况下，能够准确地找到子区域的对应位置，广泛应用于各种DIC测量场景。梯度相关算法也是基于灰度的搜索算法之一，它利用图像的梯度信息来提高搜索的准确性和效率。该算法认为图像的梯度能够更敏感地反映图像的局部特征变化，通过计算参考子区域和候选子区域的梯度相关性来确定对应位置。梯度相关算法在处理具有明显边缘或纹理特征的图像时表现出色，能够更快地收敛到正确的对应位置，减少搜索时间。基于特征的搜索算法则侧重于提取图像中的显著特征，如角点、边缘等，通过匹配这些特征来确定子区域的对应关系。这类算法首先使用特定的特征提取算子，如SIFT（尺度不变特征变换）、SURF（加速稳健特征）等，在参考图像和变形后图像中提取特征点。然后，根据特征点的描述子，如SIFT特征点的128维描述子，通过某种匹配准则，如最近邻匹配、KD树匹配等，找到两幅图像中特征点的对应关系。最后，利用这些对应特征点来确定子区域的对应位置。SIFT算法是一种非常经典的基于特征的搜索算法，它具有尺度不变性、旋转不变性和光照不变性等优点。SIFT算法首先在不同尺度空间上检测特征点，通过高斯差分金字塔来实现尺度空间的构建。然后，对每个特征点计算其主方向，生成特征点的描述子。在匹配阶段，通过计算参考图像和变形后图像中特征点描述子的欧氏距离，寻找最近邻和次近邻匹配点，根据最近邻距离与次近邻距离的比值来判断匹配的可靠性。SIFT算法在处理大变形、旋转和光照变化较大的图像时具有较高的准确性，但计算复杂度较高，处理时间较长。在实际应用中，选择合适的搜索算法需要综合考虑多种因素。对于变形较小、图像质量较好且灰度分布较为均匀的情况，基于灰度的搜索算法，如NCC算法，通常能够满足精度要求，且计算效率较高；而当图像存在较大变形、旋转或光照变化时，基于特征的搜索算法，如SIFT算法，能够提供更可靠的匹配结果，但可能需要付出更多的计算时间和资源。在岩土工程的土体单轴压缩试验中，如果土体表面散斑图像的变形相对较小，光照条件稳定，使用NCC算法可以快速准确地计算出土体表面的位移和应变；但如果在复杂的现场环境下，土体表面图像受到光照不均、灰尘干扰等影响，且土体发生较大变形时，采用SIFT算法可能更有助于准确获取土体的变形信息。2.4基于N-R迭代算法与粒子群优化算法的DIC方法2.4.1基本思想数字图像相关（DIC）方法在测量物体表面位移和应变时，面临着如何提高计算精度和效率的问题。传统的DIC方法在处理复杂变形和噪声干扰时，其计算精度和稳定性可能会受到影响。为了克服这些局限性，将N-R（牛顿-拉夫逊）迭代算法与粒子群优化算法相结合，对DIC方法进行改进，以提升其在复杂情况下的测量性能。N-R迭代算法是一种经典的数值迭代算法，常用于求解非线性方程。在DIC方法中，它通过不断迭代逼近，寻找使目标函数（如灰度相关函数）达到最优值的位移参数。其基本原理是基于泰勒级数展开，将非线性函数在当前迭代点附近进行线性化近似，从而构建迭代公式来逐步更新位移参数。在每次迭代中，根据当前的位移估计值，计算目标函数的梯度和海森矩阵，利用这些信息来确定下一次迭代的位移增量，使目标函数的值不断减小，直至满足收敛条件。N-R迭代算法具有收敛速度快的优点，能够在较少的迭代次数内接近最优解，但它对初始值的选取较为敏感，如果初始值与真实解相差较大，可能会导致迭代过程陷入局部最优解，无法收敛到全局最优。粒子群优化算法是一种基于群体智能的优化算法，模拟了鸟群觅食等群体行为。在该算法中，每个优化问题的解被看作是搜索空间中的一个“粒子”，所有粒子都有一个由目标函数决定的适应值，并且每个粒子都有自己的速度和位置。粒子通过跟踪两个“极值”来更新自己的速度和位置：一个是粒子自身所找到的最优解，即个体极值pBest；另一个是整个种群目前找到的最优解，即全局极值gBest。粒子根据自身经验（pBest）和群体经验（gBest）来调整飞行方向和距离，从而在解空间中不断搜索更优解。粒子群优化算法具有全局搜索能力强、对初始值不敏感等优点，能够在较大的解空间中快速找到较优解，但在后期收敛速度较慢，容易出现搜索停滞的情况。将N-R迭代算法与粒子群优化算法相结合应用于DIC方法，旨在充分发挥两者的优势。利用粒子群优化算法的全局搜索能力，在较大的解空间中快速搜索到一个接近全局最优解的区域，为N-R迭代算法提供一个较好的初始值。然后，借助N-R迭代算法收敛速度快的特点，在粒子群优化算法找到的较优区域内进行精细搜索，快速收敛到全局最优解，从而提高DIC方法计算位移和应变的精度和效率。在对复杂形状的岩土体进行变形测量时，首先通过粒子群优化算法在大范围的位移参数空间中进行搜索，找到大致的位移范围，再利用N-R迭代算法在这个范围内进行精确计算，能够更准确地获取岩土体表面的位移和应变信息。2.4.2算法流程改进后的基于N-R迭代算法与粒子群优化算法的DIC方法的具体计算步骤和流程如下：图像采集与预处理：使用高分辨率相机采集物体变形前后的散斑图像。对采集到的图像进行预处理，包括灰度化、去噪、增强对比度等操作，以提高图像质量，减少噪声对后续计算的干扰。粒子群优化算法初始化：在DIC方法中，将位移参数（如x方向和y方向的位移分量u和v）作为粒子的位置变量。随机生成一定数量的粒子，每个粒子代表一组可能的位移参数解。初始化粒子的速度和位置，速度决定了粒子在解空间中的移动方向和步长，位置则表示当前的位移参数估计值。同时，将每个粒子的个体极值pBest初始化为其当前位置，将全局极值gBest初始化为所有粒子中适应值最优的粒子位置。粒子群优化算法迭代：在每一代的迭代中，首先计算每个粒子的适应值。适应值通过计算粒子所代表的位移参数对应的参考图像子区域与变形后图像中相应子区域的灰度相关性来确定，常用的灰度相关函数如归一化互相关函数。将当前粒子的适应值与个体极值pBest的适应值进行比较，如果当前适应值更优，则更新个体极值pBest为当前粒子的位置和适应值。然后，将所有粒子的个体极值pBest进行比较，找出其中适应值最优的粒子，将其位置更新为全局极值gBest。根据粒子群优化算法的速度和位置更新公式，对每个粒子的速度和位置进行更新。速度更新公式为：v_{ij}(t+1)=w\timesv_{ij}(t)+c_1r_1(pBest_{ij}(t)-p_{ij}(t))+c_2r_2(gBest_{gj}(t)-p_{ij}(t))其中，v_{ij}(t+1)表示第t+1代时第i个粒子在第j维的速度；w为惯性权重系数，用于平衡粒子的全局搜索和局部搜索能力；c_1和c_2为学习因子，通常在0到2之间取值，分别表示粒子对自身经验和群体经验的学习程度；r_1和r_2为[0,1]区间上均匀分布的随机数；pBest_{ij}(t)表示第t代时第i个粒子在第j维所寻找到的最佳位置坐标值；gBest_{gj}(t)表示第t代时整个粒子群在第j维所寻找到的最佳位置坐标值；p_{ij}(t)表示第t代时第i个粒子在第j维的位置坐标值。位置更新公式为：位置更新公式为：p_{ij}(t+1)=p_{ij}(t)+v_{ij}(t+1)其中，p_{ij}(t+1)表示第t+1代时第i个粒子在第j维的新位置坐标值。通过不断迭代，粒子群逐渐向全局最优解靠近。当满足迭代终止条件（如达到最大迭代次数或适应值收敛）时，粒子群优化算法结束，得到一组较为接近全局最优解的位移参数估计值。通过不断迭代，粒子群逐渐向全局最优解靠近。当满足迭代终止条件（如达到最大迭代次数或适应值收敛）时，粒子群优化算法结束，得到一组较为接近全局最优解的位移参数估计值。N-R迭代算法初始化：将粒子群优化算法得到的全局极值gBest作为N-R迭代算法的初始值，即初始的位移参数估计值。同时，计算参考图像子区域与变形后图像中相应子区域在当前位移参数下的灰度相关函数，并对其进行泰勒级数展开，得到线性化的近似函数。N-R迭代算法迭代：在每次N-R迭代中，根据当前的位移参数估计值，计算灰度相关函数关于位移参数的梯度和海森矩阵。梯度反映了函数在当前点的变化率，海森矩阵则描述了函数的曲率信息。利用梯度和海森矩阵构建迭代公式，计算位移参数的增量：\Deltap=-H^{-1}g其中，\Deltap为位移参数的增量向量；H为海森矩阵；g为梯度向量。根据位移参数的增量更新当前的位移参数估计值：根据位移参数的增量更新当前的位移参数估计值：p_{k+1}=p_k+\Deltap其中，p_{k+1}为下一次迭代的位移参数估计值；p_k为当前迭代的位移参数估计值。重复上述步骤，不断迭代更新位移参数，直到满足收敛条件（如位移参数的增量小于设定的阈值或灰度相关函数的变化小于设定值），此时得到的位移参数即为最终的精确解。重复上述步骤，不断迭代更新位移参数，直到满足收敛条件（如位移参数的增量小于设定的阈值或灰度相关函数的变化小于设定值），此时得到的位移参数即为最终的精确解。应变计算：根据计算得到的位移参数，利用应变计算公式计算物体表面各点的应变。对于二维问题，x方向的正应变\varepsilon_{xx}、y方向的正应变\varepsilon_{yy}和切应变\gamma_{xy}的计算公式如前文所述。结果输出与分析：将计算得到的位移和应变结果进行可视化处理，如绘制位移云图、应变云图等，直观地展示物体表面的变形情况。同时，对结果进行分析，与理论值或其他测量方法得到的数据进行对比，评估改进后DIC方法的测量精度和可靠性。通过上述算法流程，将N-R迭代算法与粒子群优化算法相结合，实现了对DIC方法的改进，提高了其在测量物体表面位移和应变时的精度和效率，能够更好地满足复杂工程应用的需求。在对金属材料进行拉伸试验时，利用改进后的DIC方法，可以更准确地测量材料在拉伸过程中的位移和应变变化，为材料力学性能的研究提供更可靠的数据支持。2.5本章小结数字图像相关（DIC）方法作为一种非接触式的光学测量技术，在材料力学、岩土工程等领域具有广泛的应用前景。其基本原理是基于物体表面散斑图像的像素匹配，通过跟踪变形前后散斑图像中同一像素点的位置变化来获取物体表面的位移向量，进而计算出全场位移和应变分布。这一原理使其能够在不接触物体的情况下，精确地测量物体表面的变形情况，避免了传统接触式测量方法对物体表面的扰动和破坏。在DIC方法中，相关运算公式是实现位移和应变计算的核心。通过对位移分量进行偏导数运算，能够准确地计算出正应变和切应变，从而全面了解物体的变形状态。在二维平面问题中，通过对x和y方向位移分量u和v的偏导数计算，得到x方向正应变\varepsilon_{xx}、y方向正应变\varepsilon_{yy}以及切应变\gamma_{xy}，这些应变参数为分析物体的力学行为提供了关键数据。相关搜索方法是DIC技术的关键环节之一，不同的搜索算法适用于不同的场景。基于灰度的搜索算法，如归一化互相关（NCC）算法，利用子区域内像素的灰度信息来度量图像之间的相似性，对光照变化具有一定的鲁棒性，在图像灰度变化较为均匀的情况下能够准确地找到子区域的对应位置；基于特征的搜索算法，如SIFT算法，侧重于提取图像中的显著特征，在处理大变形、旋转和光照变化较大的图像时具有较高的准确性。将N-R迭代算法与粒子群优化算法相结合应用于DIC方法，充分发挥了两者的优势。粒子群优化算法的全局搜索能力能够在较大的解空间中快速搜索到接近全局最优解的区域，为N-R迭代算法提供较好的初始值；N-R迭代算法收敛速度快的特点则能够在粒子群优化算法找到的较优区域内进行精细搜索，快速收敛到全局最优解，从而提高DIC方法计算位移和应变的精度和效率。在实际应用中，这种改进后的DIC方法能够更准确地测量物体表面的位移和应变，为科学研究和工程实践提供更可靠的数据支持。在本研究中，DIC方法对于单轴压缩湿砂样主应变轴偏转规律的研究具有高度的适用性。它能够实时、全面地监测湿砂样表面在单轴压缩过程中的变形情况，精确获取主应变轴的偏转角等关键数据。与传统测量方法相比，DIC方法不仅能够提供全场变形信息，还能避免对湿砂样的干扰，确保实验数据的准确性和可靠性。通过DIC方法，能够直观地观察到湿砂样在不同工况下主应变轴的偏转过程，深入分析其与含水率、干密度、加载速率等因素之间的关系，为揭示砂土的力学行为和变形机制提供了有力的技术手段。三、主应变轴偏转角的DIC计算3.1主应变轴偏转角的理论公式在材料力学中，对于二维平面应力状态，主应变与切应变、正应变之间存在特定的关系，主应变轴偏转角正是基于这些关系推导得出。假设在某一平面内，已知x方向正应变\varepsilon_{xx}、y方向正应变\varepsilon_{yy}以及切应变\gamma_{xy}，主应变轴偏转角\theta的计算公式推导如下。首先，根据材料力学理论，主应变\varepsilon_1和\varepsilon_2（\varepsilon_1\geq\varepsilon_2）可通过以下公式计算：\varepsilon_{1,2}=\frac{\varepsilon_{xx}+\varepsilon_{yy}}{2}\pm\sqrt{(\frac{\varepsilon_{xx}-\varepsilon_{yy}}{2})^2+(\frac{\gamma_{xy}}{2})^2}主应变轴偏转角\theta是指主应变方向与x轴正方向之间的夹角，其计算公式为：\tan2\theta=\frac{\gamma_{xy}}{\varepsilon_{xx}-\varepsilon_{yy}}由此可得：\theta=\frac{1}{2}\arctan(\frac{\gamma_{xy}}{\varepsilon_{xx}-\varepsilon_{yy}})该公式表明，主应变轴偏转角\theta取决于切应变\gamma_{xy}以及x方向和y方向正应变的差值\varepsilon_{xx}-\varepsilon_{yy}。当\gamma_{xy}=0且\varepsilon_{xx}\neq\varepsilon_{yy}时，\tan2\theta=0，此时\theta=0或\theta=90^{\circ}，表示主应变方向与坐标轴方向一致；当\gamma_{xy}\neq0时，主应变轴将发生偏转，偏转角\theta由上述公式确定。在实际应用中，通过数字图像相关（DIC）方法测量得到湿砂样表面各点的\varepsilon_{xx}、\varepsilon_{yy}和\gamma_{xy}，代入上述公式即可计算出主应变轴偏转角\theta，从而分析湿砂样在单轴压缩过程中主应变轴的偏转规律。例如，在某一时刻，通过DIC测量得到湿砂样表面某点的\varepsilon_{xx}=0.005，\varepsilon_{yy}=0.003，\gamma_{xy}=0.002，将这些值代入公式计算可得：\tan2\theta=\frac{0.002}{0.005-0.003}=1则2\theta=45^{\circ}，\theta=22.5^{\circ}，即该点主应变轴相对于x轴正方向偏转了22.5^{\circ}。通过对湿砂样表面多个点的主应变轴偏转角计算，可以全面了解湿砂样在单轴压缩过程中主应变轴的偏转情况，为深入研究砂土的力学行为提供重要依据。3.2湿砂样应变场的DIC计算在单轴压缩湿砂样实验中，利用数字图像相关（DIC）方法计算湿砂样的应变场是研究主应变轴偏转规律的关键环节。其计算过程基于DIC方法的基本原理，通过对湿砂样表面散斑图像的处理和分析来实现。实验开始前，需在湿砂样表面制作随机分布的散斑图案。这一过程通常采用喷涂的方式，先在砂样表面均匀喷涂一层白色底漆，待底漆完全干燥后，再以适当的距离和角度喷涂黑色斑点，形成对比度清晰、分布均匀的散斑场。这样的散斑图案能够为DIC测量提供丰富的变形信息载体。采用高分辨率相机对湿砂样进行拍摄。在单轴压缩加载过程中，按照一定的时间间隔或位移增量，连续拍摄湿砂样表面的散斑图像，以记录湿砂样在不同加载阶段的变形状态。为了确保拍摄的准确性和稳定性，相机需固定在稳定的支架上，并调整好拍摄角度和焦距，使湿砂样表面完全处于相机的视场范围内，且图像清晰、无模糊和畸变。获取散斑图像后，利用DIC分析软件对图像进行处理。软件首先对参考图像（加载前的初始图像）和变形图像（加载过程中不同阶段的图像）进行预处理，包括灰度化、去噪、增强对比度等操作，以提高图像质量，减少噪声对后续计算的干扰。将参考图像划分为众多互不重叠的小子区域，每个子区域通常为正方形或矩形，其大小根据实验精度要求和图像分辨率进行合理设置。对于每个子区域，在变形图像中通过相关搜索算法寻找与之匹配的区域，常用的搜索算法如归一化互相关（NCC）算法。NCC算法通过计算参考子区域与候选子区域的归一化互相关系数来衡量两者的相似程度，遍历变形图像的各个位置，找到互相关系数最大的区域，即为参考子区域在变形图像中的对应位置。通过这种方式，确定每个子区域在变形前后的位置变化，从而得到该子区域中心点在x方向和y方向的位移分量u和v。根据位移分量计算应变。对于二维平面问题，x方向的正应变\varepsilon_{xx}、y方向的正应变\varepsilon_{yy}和切应变\gamma_{xy}的计算如前文所述，通过对位移分量u和v求偏导数得到。在实际计算中，由于图像是离散的像素点构成，通常采用差分法进行近似求解。假设相邻像素点在x方向的间距为\Deltax，在y方向的间距为\Deltay，则x方向正应变\varepsilon_{xx}的近似计算公式为\varepsilon_{xx}\approx\frac{u(x+\Deltax,y)-u(x,y)}{\Deltax}；y方向正应变\varepsilon_{yy}的近似计算公式为\varepsilon_{yy}\approx\frac{v(x,y+\Deltay)-v(x,y)}{\Deltay}；切应变\gamma_{xy}的近似计算公式为\gamma_{xy}\approx\frac{u(x,y+\Deltay)-u(x,y)}{\Deltay}+\frac{v(x+\Deltax,y)-v(x,y)}{\Deltax}。通过这些公式，对每个子区域的位移信息进行计算，得到湿砂样表面各点的应变值，从而构建出湿砂样的应变场。在计算过程中，为了提高计算精度，可采用亚像素插值算法。该算法能够将位移计算精度提高到亚像素级别，从而更精确地反映湿砂样表面的变形情况。对计算得到的应变场数据进行滤波处理，去除可能存在的噪声和异常值，进一步提高数据的可靠性。通过上述步骤，利用DIC方法实现了对单轴压缩湿砂样应变场的准确计算，为后续分析主应变轴偏转角和研究主应变轴偏转规律提供了基础数据。3.3特殊位置测点的布置及双三次样条插值3.3.1特殊位置测点的布置原则在湿砂样表面布置测点时，需遵循特定的原则以确保能够准确获取关键信息，为研究主应变轴偏转规律提供有力支持。测点的布置应充分考虑湿砂样在单轴压缩过程中的变形特点和可能出现的特殊区域。为了准确捕捉剪切带的形成与发展过程，在可能出现剪切带的区域，如湿砂样的中部及边缘部分，应加密布置测点。根据土力学理论和相关研究，在单轴压缩下，砂土试样通常在与加载方向成一定角度的位置形成剪切带，这个角度一般在45°左右。因此，在湿砂样表面以加载方向为基准，在45°方向附近合理布置测点，能够更有效地监测剪切带内的应变变化和主应变轴偏转情况。在湿砂样的中部，由于该区域是剪切带形成的核心区域，应力和应变集中现象较为明显，布置足够数量的测点可以精确记录该区域的变形信息，为研究剪切带内主应变轴的偏转机制提供详细数据。在湿砂样的边缘部分，由于边界条件的影响，其变形情况与内部有所不同，也需要适当布置测点。边缘区域可能会受到试验机加载板与砂样之间的摩擦力、砂样与周围环境的相互作用等因素的影响，导致变形不均匀。通过在边缘区域布置测点，可以监测这些特殊因素对主应变轴偏转的影响，全面了解湿砂样在单轴压缩下的变形特性。除了关注剪切带和边缘区域，还应在湿砂样表面均匀布置一定数量的测点，以获取整体的变形信息。这些均匀分布的测点可以反映湿砂样在不同位置的变形趋势，与特殊区域的测点数据相互补充，有助于分析主应变轴在整个湿砂样表面的偏转规律。在湿砂样表面划分网格，将测点布置在网格的节点上，既能保证测点的均匀分布，又便于后续的数据处理和分析。为了提高测量精度和可靠性，测点的布置应具有一定的冗余性。在关键位置布置多个测点，当某个测点的数据出现异常时，可以参考相邻测点的数据进行修正或验证，确保测量结果的准确性。在可能出现较大变形梯度的区域，如剪切带与非剪切带的交界处，布置多个测点，能够更准确地捕捉应变的变化情况，减少测量误差对研究结果的影响。3.3.2双三次样条插值方法在获取湿砂样表面测点的应变数据后，由于测点数量有限，无法全面、连续地描述湿砂样表面的应变分布，因此采用双三次样条插值方法对测点数据进行处理，以得到更连续、更精确的应变分布信息。双三次样条插值是一种在二维平面上进行数据插值的方法，它使用局部的三次多项式来逼近函数值，通过保证插值函数在节点处的连续性和光滑性，生成新的数据点，从而实现对离散数据的平滑插值。假设在湿砂样表面有一系列离散的测点，其坐标为(x_i,y_j)，对应的应变值为z_{ij}，其中i=1,2,\cdots,n，j=1,2,\cdots,m。双三次样条插值的基本步骤如下：建立网格：以测点的坐标(x_i,y_j)为节点，在湿砂样表面建立一个二维网格，将整个表面划分为多个小矩形区域。计算控制点：对于每个小矩形区域，需要计算出其四个角点处的插值多项式的系数，这些系数构成了控制点。通过已知的测点应变值z_{ij}以及一些边界条件（如插值函数在边界上的一阶导数和二阶导数等），利用三次样条函数的性质和相关数学公式，求解线性方程组，得到每个小矩形区域四个角点处的控制点。应用双三次插值：在每个小矩形区域内，根据计算得到的控制点，使用双三次样条插值公式计算任意位置(x,y)处的应变值z(x,y)。双三次样条插值公式通常基于Lagrange多项式或其他相关的三次多项式形式，通过对控制点进行加权组合，得到插值点的函数值。假设在某个小矩形区域内，四个角点的坐标分别为(x_1,y_1)，(x_1,y_2)，(x_2,y_1)，(x_2,y_2)，对应的控制点分别为c_{11}，c_{12}，c_{21}，c_{22}，则该区域内任意点(x,y)处的插值公式可以表示为：z(x,y)=\sum_{i=1}^{2}\sum_{j=1}^{2}c_{ij}B_{i}(x)B_{j}(y)其中，B_i(x)和B_j(y)是关于x和y的三次样条基函数，它们根据小矩形区域的坐标范围和插值条件确定。通过这种方式，在每个小矩形区域内实现了对离散测点数据的平滑插值，得到了该区域内连续的应变分布。生成插值结果：对湿砂样表面的所有小矩形区域进行双三次样条插值计算，最终得到整个湿砂样表面连续的应变分布场。将插值结果可视化，如绘制应变云图，能够直观地展示湿砂样表面应变的分布情况，便于分析主应变轴的偏转规律。在实际应用中，为了提高双三次样条插值的计算效率和精度，可以采用一些优化算法和技巧。在计算控制点时，使用高效的线性方程组求解算法，如LU分解法、共轭梯度法等，以减少计算时间。对测点数据进行预处理，去除噪声和异常值，提高数据的质量，从而提升插值结果的准确性。通过双三次样条插值方法，能够在有限的测点数据基础上，更准确地描述湿砂样表面的应变分布，为深入研究主应变轴偏转规律提供更丰富、更精确的数据支持。3.4本章小结主应变轴偏转角的理论公式基于二维平面应力状态下主应变与切应变、正应变的关系推导得出，为后续通过实验数据计算主应变轴偏转角提供了理论基础。公式中，主应变轴偏转角取决于切应变以及x方向和y方向正应变的差值，当切应变不为零时，主应变轴将发生偏转。利用数字图像相关（DIC）方法计算湿砂样应变场，通过在湿砂样表面制作散斑图案，使用高分辨率相机拍摄不同加载阶段的散斑图像，经DIC分析软件处理，运用相关搜索算法确定子区域位置变化，进而计算出位移分量，再通过差分法近似求解得到应变值，从而构建湿砂样的应变场。在湿砂样表面布置测点时，遵循在可能出现剪切带区域和边缘部分加密布置，同时均匀布置一定数量测点以获取整体变形信息，且保证测点布置具有冗余性的原则，确保能够准确获取关键信息。采用双三次样条插值方法对测点数据进行处理，通过建立网格、计算控制点、应用双三次插值公式，得到整个湿砂样表面连续的应变分布场，为研究主应变轴偏转规律提供更精确的数据。通过本章内容，成功建立了基于DIC方法计算主应变轴偏转角的体系，为后续开展不同工况下单轴压缩湿砂样主应变轴偏转规律的实验研究奠定了坚实的基础。四、实验过程与方法4.1实验材料与准备本实验所选用的砂土取自[具体地点]，该地区的砂土具有一定的代表性，广泛分布于当地的各类地质工程现场，对其力学性质的研究具有重要的实际意义。为全面了解该砂土的物理性质，进行了一系列详细的测试分析。通过筛分试验测定其颗粒级配，结果显示该砂土中粒径大于0.5mm的颗粒质量占总质量的[X1]%，粒径大于0.25mm的颗粒质量占总质量的[X2]%，粒径大于0.075mm的颗粒质量占总质量的[X3]%，根据砂土的分类标准，该砂土属于[具体砂土类别，如中砂]。采用比重瓶法测定其比重为[具体比重数值]，利用烘干法测得其天然含水率为[具体含水率数值]，通过环刀法测定其天然密度为[具体密度数值]。在制作湿砂样模型时，严格遵循科学的步骤，以确保砂样的质量和实验结果的准确性。首先，将采集到的砂土自然风干，去除其中的杂质和大颗粒，使砂土颗粒均匀。然后，按照预定的含水率要求，计算所需添加的水量，并将水均匀喷洒在砂土上。为使砂土与水充分混合，采用人工搅拌和机械搅拌相结合的方式，搅拌时间持续[具体搅拌时间]，确保砂土的含水率均匀一致。将混合好的湿砂土分多次装入定制的圆柱形模具中，模具内径为[具体内径数值]mm，高度为[具体高度数值]mm，材料选用高强度的有机玻璃，以保证在加载过程中模具不会发生变形，从而确保砂样受力均匀。在装样过程中，使用捣棒分层夯实砂土，每层砂土的厚度控制在[具体厚度数值]mm左右，每层夯实次数为[具体夯实次数]次，以保证砂样的干密度均匀。装样完成后，用刮刀将砂样表面刮平，使其与模具顶部平齐。将装有砂样的模具放置在压力机上，施加[具体预压荷载数值]kN的预压荷载，保持[具体预压时间]min，以消除砂样内部的初始孔隙和应力不均匀性，使砂样在后续实验中能够更加稳定地受力。预压完成后，小心地将砂样从模具中取出，得到成型的湿砂样。为了利用数字图像相关（DIC）方法准确测量湿砂样在单轴压缩过程中的变形，在湿砂样表面制作散斑图案。首先，在砂样表面均匀喷涂一层白色底漆，采用专业的喷枪设备，控制喷涂压力为[具体喷涂压力数值]MPa，喷枪与砂样表面的距离保持在[具体距离数值]cm，以确保底漆均匀覆盖，厚度约为[具体底漆厚度数值]mm。待白色底漆完全干燥后，再以适当的角度和距离喷涂黑色斑点，形成对比度清晰、分布均匀的散斑场。黑色斑点的平均直径控制在[具体斑点直径数值]mm左右，斑点之间的平均间距为[具体间距数值]mm，这样的散斑图案能够为DIC测量提供丰富的变形信息，确保在实验过程中能够准确追踪砂样表面的位移和应变变化。4.2砂土模型制作及DIC实验观测在完成实验材料准备后，着手制作砂土模型。将混合均匀且含水率达标的湿砂土分多层装入特制的模具中。模具选用高强度有机玻璃材质，呈内径[X]mm、高度[X]mm的圆柱形，其尺寸经过精心设计，既能满足实验对砂样尺寸的要求，又便于DIC系统进行图像采集和分析。在装样过程中，严格控制每层砂土的厚度为[X]mm左右，使用捣棒以[X]次/层的频率夯实，确保砂样干密度均匀。装样完成后，用刮刀仔细刮平砂样表面，使其与模具顶部平齐，保证砂样在后续实验中受力均匀。将装有砂样的模具放置在高精度压力机上，施加[X]kN的预压荷载，保持[X]min，以消除砂样内部的初始孔隙和应力不均匀性，使砂样在后续实验中能够更加稳定地受力。预压完成后，小心地将砂样从模具中取出，得到成型的湿砂样。为利用数字图像相关（DIC）方法准确测量湿砂样在单轴压缩过程中的变形，在湿砂样表面制作散斑图案。首先，使用专业喷枪在砂样表面均匀喷涂一层白色底漆，喷涂压力控制在[X]MPa，喷枪与砂样表面的距离保持在[X]cm，确保底漆均匀覆盖，厚度约为[X]mm。待白色底漆完全干燥后，再以适当的角度和距离喷涂黑色斑点，形成对比度清晰、分布均匀的散斑场。黑色斑点的平均直径控制在[X]mm左右，斑点之间的平均间距为[X]mm，这样的散斑图案能够为DIC测量提供丰富的变形信息，确保在实验过程中能够准确追踪砂样表面的位移和应变变化。在进行DIC实验观测时，搭建了一套完善的DIC测量系统。该系统主要由高分辨率相机、稳定的相机支架、专业图像采集软件以及数据分析软件组成。将制作好散斑图案的湿砂样放置在单轴压缩试验机的加载平台上，调整湿砂样的位置，使其中心与加载平台中心重合，确保加载过程中湿砂样受力均匀。将高分辨率相机固定在稳定的支架上，调整相机的位置和角度，使湿砂样表面完全处于相机的视场范围内，并且保证相机光轴与湿砂样表面垂直，以减少图像畸变。根据实验要求，设置相机的参数，如分辨率、帧率、曝光时间等。本实验中，相机分辨率设置为[具体分辨率数值]，帧率为[具体帧率数值]fps，曝光时间根据实验环境光照条件调整为[具体曝光时间数值]s，以确保能够清晰、准确地采集湿砂样表面的散斑图像。在单轴压缩加载过程中，按照一定的时间间隔或位移增量，通过图像采集软件控制相机连续拍摄湿砂样表面的散斑图像。本实验设定每隔[具体时间间隔数值]s或每加载[具体位移增量数值]mm拍摄一张图像，以记录湿砂样在不同加载阶段的变形状态。同时，利用单轴压缩试验机的控制系统，精确控制加载速率，分别设置加载速率为[具体加载速率数值1]mm/min、[具体加载速率数值2]mm/min、[具体加载速率数值3]mm/min等，研究不同加载速率对湿砂样主应变轴偏转规律的影响。在整个实验过程中，密切关注实验设备的运行状态和湿砂样的变形情况，确保实验的顺利进行。4.3实验加载方案本实验采用位移控制加载方式，通过高精度位移控制装置，以恒定的速率对湿砂样施加轴向位移。加载速率分别设置为0.1mm/min、0.5mm/min和1mm/min三个等级。选择这三个加载速率，是基于对砂土力学行为研究的经验以及实际工程应用的参考。在岩土工程中，不同的施工过程和自然作用对砂土的加载速率差异较大，如基础施工过程中对地基的加载速率相对较慢，而地震等自然灾害对砂土的加载则较为快速。0.1mm/min的加载速率模拟了较为缓慢的加载过程，类似于长期的地基沉降过程；0.5mm/min的加载速率处于中等水平，可代表一些常规的工程加载情况；1mm/min的加载速率则模拟了相对快速的加载，类似于地震等短时间内的快速加载情况。通过设置这三个不同等级的加载速率，可以全面研究加载速率对湿砂样主应变轴偏转规律的影响。在加载过程中，当湿砂样的轴向应变达到15%时，停止加载。选择15%的轴向应变作为加载终止条件，是综合考虑了砂土的力学特性和实验目的。在砂土的单轴压缩过程中，当轴向应变达到一定程度后，砂土的变形趋于稳定，其力学行为基本不再发生显著变化。根据相关研究和前期预实验结果，15%的轴向应变能够使湿砂样充分展现其在不同加载条件下的变形特性，同时避免因过度加载导致砂样破坏严重，影响数据的准确性和分析结果。在加载过程中，实时监测湿砂样的变形情况和荷载变化，确保加载过程的稳定性和数据采集的准确性。4.4计算方式在实验数据处理阶段，利用数字图像相关（DIC）方法获取的湿砂样表面位移数据，依据材料力学中的应变计算公式来计算应变。对于二维平面问题，假设在参考图像中某点的坐标为(x,y)，变形后该点移动到(x+u,y+v)，其中u和v分别为该点在x方向和y方向的位移分量。x方向的正应变\varepsilon_{xx}通过对u关于x求偏导数得到，即：\varepsilon_{xx}=\frac{\partialu}{\partialx}在实际计算中，由于图像是离散的像素点构成，采用差分法进行近似求解，假设相邻像素点在x方向的间距为\Deltax，则近似计算公式为：\varepsilon_{xx}\approx\frac{u(x+\Deltax,y)-u(x,y)}{\Deltax}同理，y方向的正应变\varepsilon_{yy}为：\varepsilon_{yy}=\frac{\partialv}{\partialy}其差分近似计算公式为：\varepsilon_{yy}\approx\frac{v(x,y+\Deltay)-v(x,y)}{\Deltay}切应变\gamma_{xy}的计算公式为：\gamma_{xy}=\frac{\partialu}{\partialy}+\frac{\partialv}{\partialx}差分近似计算式为：\gamma_{xy}\approx\frac{u(x,y+\Deltay)-u(x,y)}{\Deltay}+\frac{v(x+\Deltax,y)-v(x,y)}{\Deltax}在得到x方向正应变\varepsilon_{xx}、y方向正应变\varepsilon_{yy}以及切应变\gamma_{xy}后，根据主应变轴偏转角的理论公式计算主应变轴偏转角\theta。主应变轴偏转角\theta是指主应变方向与x轴正方向之间的夹角，计算公式为：\theta=\frac{1}{2}\arctan(\frac{\gamma_{xy}}{\varepsilon_{xx}-\varepsilon_{yy}})在计算过程中，为确保计算结果的准确性，对位移数据进行了滤波处理，采用高斯滤波等方法去除噪声干扰，避免因噪声导致的计算误差。同时，为提高计算精度，对DIC测量系统进行了严格的标定，考虑了相机镜头的畸变、像素尺寸等因素对测量结果的影响，并通过标定参数对计算结果进行修正。4.5本章小结在本次实验中，选用了具有代表性的[具体地点]砂土，通过筛分试验、比重瓶法、烘干法、环刀法等多种方法对其物理性质进行了全面测定，明确了该砂土的颗粒级配、比重、天然含水率和天然密度等关键参数，为后续实验提供了重要的材料基础。在制作湿砂样模型时，严格按照预定的含水率要求，通过均匀喷洒水、充分搅拌等步骤，确保砂土含水率均匀。装样过程中，分层夯实砂土，控制每层厚度和夯实次数，保证砂样干密度均匀。预压处理消除了砂样内部的初始孔隙和应力不均匀性，使砂样在实验中受力更稳定。在砂样表面制作散斑图案，通过精确控制喷涂工艺，形成了对比度清晰、分布均匀的散斑场，为DIC测量提供了良好的条件。搭建了完善的DIC测量系统，将制作好散斑图案的湿砂样放置在单轴压缩试验机加载平台上，调整好相机位置、角度和参数，在加载过程中按照设定的时间间隔或位移增量拍摄散斑图像，同时精确控制加载速率，为研究不同加载速率对湿砂样主应变轴偏转规律的影响提供了数据支持。实验采用位移控制加载方式，设置0.1mm/min、0.5mm/min和1mm/min三个加载速率等级，以模拟不同的工程加载情况。当湿砂样轴向应变达到15%时停止加载，确保能够充分展现湿砂样的变形特性，同时避免过度加载对数据准确性的影响。在数据处理阶段，利用DIC方法获取的位移数据，依据材料力学应变计算公式计算应变，采用差分法近似求解，并对位移数据进行滤波处理和系统标定，以提高计算精度。根据主应变轴偏转角理论公式计算主应变轴偏转角，为后续分析主应变轴偏转规律提供了数据基础。五、实验结果分析5.1观测平面内应变场的分布规律通过数字图像相关（DIC）方法，对不同加载阶段湿砂样观测平面内的应变场进行了精确测量与深入分析。在加载初期，湿砂样整体变形较为均匀，观测平面内的应变分布也相对均匀，各位置处的应变值差异较小。此时，由于荷载较小，砂土颗粒之间的相对位移较小，主要表现为弹性变形阶段。随着加载的持续进行，应变场的分布逐渐发生变化。在湿砂样的中部区域，应变值开始逐渐增大，且增长速率较快，呈现出明显的应变集中现象。这是因为在单轴压缩作用下，湿砂样中部受到的应力较为集中，砂土颗粒之间的相互作用增强，导致该区域的变形加剧。与此同时，在湿砂样的边缘部分，应变值相对较小，增长速率也较为缓慢。这是由于边缘部分受到试验机加载板与砂样之间的摩擦力以及砂样与周围环境相互作用的影响，其受力状态与内部有所不同，使得变形相对较小。在加载接近尾声，即湿砂样的轴向应变达到15%左右时，应变场的分布呈现出更为明显的不均匀性。在湿砂样的中部，应变值达到较大数值，且在某些局部区域出现了应变突变的情况，这些区域往往是剪切带的萌生和发展区域。剪切带的形成使得砂土颗粒之间的相对位移急剧增大，导致应变集中程度进一步加剧。而在湿砂样的边缘部分，虽然应变值仍然相对较小，但在靠近剪切带的区域，应变值也出现了一定程度的增大，这表明剪切带的影响范围逐渐扩大。为了更直观地展示观测平面内应变场的分布规律，绘制了不同加载阶段的应变云图。从应变云图中可以清晰地看到，随着加载的进行，应变高值区域逐渐从湿砂样的中部向周围扩展，形成了以中部为中心的应变集中区域。同时，应变云图的颜色变化也反映了应变值的大小和分布情况，颜色越深表示应变值越大，颜色越浅表示应变值越小。通过对应变云图的分析，可以更准确地把握湿砂样在单轴压缩过程中观测平面内应变场的分布特征和变化趋势。不同加载速率对观测平面内应变场的分布也有一定影响。在加载速率较低时，如0.1mm/min，湿砂样有足够的时间来调整内部结构以适应荷载变化，应变场的分布相对较为均匀，应变集中现象相对不明显。而在加载速率较高时，如1mm/min，湿砂样内部结构来不及充分调整，应变集中现象更为显著，应变场的不均匀性也更加突出。加载速率的变化还会影响剪切带的形成和发展速度，加载速率越快，剪切带的形成和扩展速度也越快，导致应变场的变化更加迅速。5.2观测平面内主应变轴偏转角的分布规律通过对不同加载阶段观测平面内主应变轴偏转角的详细分析，发现其分布呈现出特定的规律。在加载初期，湿砂样整体变形较小，主应变轴偏转角的分布相对较为均匀，大部分区域的偏转角较小，集中在较小的角度范围内，如-5°至5°之间。这是因为在加载初期，砂土颗