高速3D-DIC这名字听起来有点绕但对搞岩石力学试验的人来说它基本算得上近几年最值得上手的观测工具之一。我最初接触这个组合是因为手里的凝灰岩试件在干湿循环下实在太“不老实”——表面开裂、掉渣、局部鼓起传统应变片根本贴不住引伸计也只能测一个方向的平均变形等看到破坏的时候整个过程早就过去了。后来把高速3D-DIC接进试验系统才算把凝灰岩从“开始损伤”到“全面破坏”这一段过程完整地看清楚了。这篇就把我这套试验从方案设计、试件制备、系统搭到数据处理的全过程整理出来重点讲为什么选高速3D-DIC、参数怎么定、有哪些坑必须躲给也准备做类似试验的同仁一个可以直接照着走的参考。1. 研究思路与试验方案设计1.1 为什么盯着凝灰岩和干湿循环不放凝灰岩是火山碎屑岩的一种典型的软岩代表。它里面含有大量火山灰和岩屑孔隙率普遍不低部分层位还夹杂蒙脱石、沸石这类遇水敏感矿物。这类岩石在干燥状态下看着挺硬朗一泡水就原形毕露吸水膨胀、强度软化、表面剥落。但凡工程上遇到凝灰岩最头疼的就是库岸边坡、渠道边坡这类长期处于水位变动区的部位。水位上升时岩石泡水水位下降时岩石被晒干这种“湿了又干、干了又湿”的反复过程我们叫它干湿循环。干湿循环是模拟自然环境风化最常用的加速试验手段。相比于长期野外监测室内干湿循环可以在几周到几个月内把几十年的湿度交替效应压缩出来。对凝灰岩来说每一次吸水和脱水都会在矿物颗粒尺度上产生膨胀与收缩这些微小的体积变化日积月累就会在岩石内部累积微裂纹最终表现为宏观力学性能的劣化。我把注意力放在损伤演化而不是单纯的强度衰减上是因为干湿循环对岩石的影响不是一步到位的。循环早期可能只是孔隙增多到了中期微裂纹开始贯通后期才显现出肉眼可见的宏观裂缝。如果只测循环后的一次性抗压强度等于只看到了最终结果中间这些过程全丢了。而损伤演化的研究价值恰恰在于搞清楚它从什么时候开始加速、在哪个区域最先萌生裂纹、破坏前有没有可观测的前兆。1.2 高速3D-DIC在这个场景里解决了什么核心痛点传统手段测岩石损伤大概有这么几种应变片、引伸计、声发射、工业相机加二维DIC。每种都有各自的问题。应变片是最常见的选择但在凝灰岩上非常难受。凝灰岩表面粗糙孔隙多应变片贴上去胶水先渗进孔里实际传递的变形量打了折扣。更麻烦的是当试件表面开始剥落或出现宏观裂缝时应变片不是跟着坏就是直接脱落偏偏这个时候正是最需要数据的阶段。引伸计类似只能测标距范围内的平均应变无法给出空间分布——它根本定位不了裂纹出现在哪里。声发射能捕捉岩石内部的微破裂信号对于“什么时候开始损伤”这个问题回答得比较好但是它不擅长回答“损伤在哪个空间位置、以什么形态发展”。声发射定位精度通常在毫米级别对试件内部小裂纹的分布描述比较模糊。二维DIC解决了全场测量的问题不接触、无损伤、能给出整个表面的位移场和应变场。但它有个根本性的短板它只认面内变形默认试件表面不发生面外位移。岩石试验恰恰是个强烈三维变形的过程。试件受压膨胀时表面会鼓出来裂纹张开时两侧会错动端部摩擦会造成局部隆起这些离面位移在加载后期占比相当大。二维DIC算出的应变场在变形大了以后误差会迅速放大结果不可靠。高速3D-DIC就是冲这几个痛点去的。它用两台高速相机同步拍摄模拟人的双目视觉通过三维重构同时获取试件表面的三维坐标场和应变场。面内、面外变形都能测而且是非接触全场测量。加上“高速”这个前缀意味着每秒可以采集几百到几千帧图像配合材料失稳破坏的极短时间窗口能捕捉到常规相机根本拍不到的瞬时破裂过程。1.3 整套试验思路的搭建逻辑我设计的试验主线是这样先制备一批标准化凝灰岩试件随机分组后分别经历不同次数的干湿循环处理比如0次、5次、10次、15次、20次五个梯度然后在每级循环后取出对应试件在万能试验机上进行单轴压缩加载。加载全程用高速3D-DIC同步拍摄同时记录载荷和位移最后把光学测量数据和力学数据结合起来做损伤定量分析。这个思路有两点特别设计过。第一干湿循环处理和力学加载分开进行循环在烘箱和浸泡箱里完成加载在试验机上进行避免了把整台试验机长期占着做循环的尴尬。第二在每级循环完成后增加一个称重和超声波波速测试环节用非破坏手段先初步评估损伤程度再上试验机做破坏性测试。这样把同批次试件在相近循环次数下的离散性尽量控制住数据前后也能互相对照。试验的总体目标不是只看强度掉了多少而是要回答三个问题干湿循环在细观尺度上怎么改造了凝灰岩的内部结构这种结构改造如何体现在宏观应力-应变行为和表面变形场上损伤变量能不能用3D-DIC的数据建立起来并且跟循环次数呈现单调的对应关系。2. 高速3D-DIC核心原理与参数选型2.1 三维DIC的匹配逻辑DIC的基本思路说穿了很简单试件表面喷上随机散斑相机连续拍照对比变形前后图像中每个小区域的灰度图案找到它在下一帧图像中移动到了哪里这个位移量就是该点的变形。散斑就像给岩石表面做了无数个微小的“指纹标记”算法追踪这些指纹在图像序列中的位置变化。二维DIC只用一个相机假设所有变形都发生在相机成像平面内。三维DIC用两个相机从不同角度同时拍摄同一区域利用双相机之间的空间关系通过标定获得对每个散斑点在左右两张图像中做立体匹配。一旦同一个物理点在左右图中的像素位置都确定了就可以通过三角测量原理算出它的三维空间坐标。对每一帧都做这个操作就得到整个表面的三维位移场再通过数值微分计算应变场。这里有个概念很关键DIC给的是试件表面或者说是可见表面的位移场不是岩石内部的。内部损伤只能靠声发射、CT这类手段去推测但表面变形是内部损伤的“投影”裂纹一旦扩展到表面DIC能立刻看到应变集中的位置和形态。干湿循环导致的损伤恰恰最先从表面开始发展所以表面全场测量在这个场景中非常契合。2.2 “高速”二字到底要多快高速的界定取决于你观测的物理过程的时间尺度。岩石单轴压缩的破坏过程从开始出现可见裂纹到完全失去承载力通常只有几十到几百毫秒。而岩石破坏前往往存在一个短暂的“准静态阶段”之后突然加速进入动态失稳。如果相机只有30帧每秒破裂过程可能就落在两帧之间等于什么都没拍到。我设计加载速率时尽量让破坏过程持续0.5到2秒左右这样高速相机在500帧每秒左右的采样率下能拍到约500到1000帧图像足以分辨裂纹从萌生、扩展到贯穿的完整过程。如果加载速率很快比如冲击加载那就要上几千帧每秒的设备如果只是观察蠕变或非常缓慢的损伤累积一两百帧每秒也够。关键是让“相机时间分辨率”匹配“破坏事件的时间尺度”。分辨率是另一个要权衡的参数。高速相机分辨率越高单个散斑点占的像素越多测量精度越好但帧率和内存会受限。我的经验是在保证帧率不低于200帧每秒的前提下尽量选择高分辨率。如果目标只是捕捉最终破裂瞬间可以降低采样帧率换分辨率如果目标还包括破坏前的应变局部化过程帧率优先级更高。2.3 凝灰岩岩性给测试带来的特殊约束凝灰岩和常见的花岗岩、大理岩在DIC测试上差别很大。花岗岩表面致密喷漆后散斑质量高图像对比度好。凝灰岩表面粗糙又多孔喷漆时漆容易渗进孔隙导致散斑边界模糊等干湿循环处理完试件表面可能已经泡得发酥一碰就掉砂散斑层也跟着遭殃。这些岩性特点直接影响两类决策。第一散斑制作工艺。对凝灰岩来说底漆和面漆都不能喷太厚否则散斑层会作为一个“壳”与岩石表面脱开测的是漆的变形而不完全是岩石的变形。但太薄又覆盖不住深色岩面图像对比度不足。要在这之间找到平衡点通常需要先做几个废试件试喷观察干燥后的散斑边界清晰度。第二干湿循环与散斑制作的时间顺序。这个问题我在后面的实操章节会详细讲但提前说结论建议先完成全部干湿循环处理再做散斑。3. 试件制备、散斑制作与干湿循环处理3.1 取样与端面加工试件的质量决定了实验的可靠程度这一步不值得省时间。首先取样要控制方向。凝灰岩如果存在层理或原生裂隙加载方向和层理的角度关系会显著影响力学行为。我在取样时就统一要求钻孔方向垂直于主要层理面保证所有试件的结构各向异性取向一致。加工精度同样关键。试件做成标准圆柱体直径50毫米高径比2比1这是在岩石力学试验规程里最常规的尺寸。端面不平整会造成加载时应力分布不均局部应力集中会直接影响表面应变场形态。我要求端面平整度控制在0.02毫米以内上下端面的平行度误差不超过0.05毫米。新加工的试件我会先检查外观有肉眼可见裂纹的直接淘汰。加工完的试件还有个细节需要处理侧面尽量保持天然干燥状态不要用油冷或者水冷过度浸泡。某些金刚石取芯钻头加工时冷却液会渗入试件改变初始含水状态。实际操作时我会在加工后把试件放进烘箱在60摄氏度下烘48小时让初始状态统一。3.2 散斑喷涂从废件试喷到成品散斑喷涂是DIC测试里看着简单实际最容易翻车的环节。散斑质量差后续计算相关系数低、位移场噪声大、甚至直接算不出来。散斑的基本要求是随机、高对比、边界清晰、尺寸合适。我给凝灰岩推荐的做法是先在试件表面均匀喷一层哑光白底漆覆盖掉岩石原本的颜色差异形成均匀的底色。待底漆完全干透后再用哑光黑漆在表面喷出细小斑点。黑斑点的大小直接决定了DIC子区匹配的精度一般要求每个散斑颗粒在图像上占据约3到5个像素。以目前的相机分辨率为例视场宽度覆盖试件直径约50毫米时单个散斑颗粒的实际物理尺寸大约在0.5到1毫米之间。喷涂手法有讲究。喷黑漆时喷枪离试件远一点短促点喷让漆雾颗粒随机落在表面而不是喷成连续的漆膜。我习惯分两到三次薄喷每次喷完等几分钟避免一次喷太厚导致漆面流淌。整个过程在通风环境下操作喷完的试件至少静置24小时让漆层充分固化。这里是我踩过最惨的一次坑第一次做干湿循环试验时我先喷了散斑再去做循环结果泡水后散斑层大面积起泡脱落剩下的部分也严重变形DIC计算出来的位移场全线飘移。从那以后我彻底改了流程一律先循环后散斑。循环处理完的试件烘干冷却后表面状态基本稳定再制作散斑这样得到的散斑层在整个加载过程中都能保持完好。3.3 干湿循环方案与控制参数干湿循环方案要贴近工程实际又要保证室内可操作性和可重复性。我采用的方案是在常温清水中浸泡24小时然后取出擦干表面水分放入烘箱在60摄氏度下烘24小时这样一次完整的干湿循环是48小时。循环次数梯度分别为0次、5次、10次、15次、20次。几个参数需要认真斟酌。烘箱温度不要超过60摄氏度。温度过高会改变岩石内部的矿物结构比如某些黏土矿物在高温下会脱去结晶水这种化学变化不是我们要研究的“干湿损伤”属于额外变量。60摄氏度配合24小时足够把试件内部水分基本蒸干再延长烘烤时间边际效应很小。浸泡阶段同样用常温清水不需要加热加速浸泡。加热浸泡虽然能加快水分渗入但也会带来温度梯度和热应力干扰损伤机制。每级循环完成后我会在试件烘干冷却后立即测量质量和纵波波速。质量变化反映干湿循环导致的矿物溶解和颗粒剥落波速变化则直接反映内部微裂纹的累积程度。这两个参数虽然不能直接给出损伤的空间分布但可以作为DIC数据的交叉验证。如果实验数据表明波速下降了DIC表面应变场中也应该能看到相应程度的应变集中。有个容易被忽略的步骤每次浸泡结束后试件拿出来要在室温下沥干表面水分再用湿布轻轻擦掉表面黏附的松散颗粒避免这些颗粒在烘干时黏结在表面形成假壳。这个假壳会影响后续散斑制作时的表面形貌。4. 加载系统搭建与高速采集参数配置4.1 相机布置、标定与光照整套系统布局需要提前想清楚因为一旦开机加载测试过程不可重复现场调整的时间非常有限。双相机的布置逻辑很有讲究。两台高速相机分别布置在试件前方的左右两侧光轴指向试件中心。两相机光轴之间的夹角也就是立体角我设定在20到25度之间。这个角度太小三维重构的深度分辨率不足离面位移测不准角度太大试件表面在左右图像中的形貌差异过大立体匹配难度显著增加。20到25度是一个平衡点既能提供足够的深度分辨率又不会让匹配算法过于吃力。相机到试件的距离取决于焦距和视场需求。我用的是定焦镜头先将视场范围调成刚好覆盖试件全高再留出约15%的余量防止试件压缩变形后超出视场。光圈适当缩小以增加景深因为试件两侧边缘距离两个相机的物距其实略有差异充足的景深能保证整个表面都在合焦范围内。标定是三维DIC精度的基石。高精度标定板按厂家推荐的姿态拍摄15到20张覆盖视场不同位置和倾斜角度。标定结果我会重点看重投影误差一般要求小于0.05像素才算合格。标定不合格就重新拍不要带着误差继续做实验。光源方面我用两个LED常亮灯从试件两侧45度方向照射目的是消除阴影并保证散斑对比度均匀。光照在加载过程中需要保持稳定这个稳定不仅指亮度不变还包括色温和频闪特征。有些LED灯存在肉眼不可见的高频闪烁高速相机一拍就会出现明暗条纹好在绝大多数DIC专用的实验室光源都经过防频闪设计。4.2 加载速率与采集帧率的匹配加载系统选的是电液伺服万能试验机可以进行位移控制。单轴压缩试验我采用位移控制模式加载速率设定为0.5毫米每分钟。这个速率在岩石力学试验中属于常规“准静态”范围既能保证加载过程平稳也不会让破坏来得过于突然。更重要的是加载速率和高速相机帧率的匹配。0.5毫米每分钟换算下来是每秒0.0083毫米的位移速率。对于50毫米高的试件应变速率约为每秒0.00017。在这种速率下岩石试件的破坏过程大约集中在最后0.5到1秒内。我设置高速相机在加载开始后连续采集帧率定为200到500帧每秒确保在破坏加速阶段仍然不会漏掉关键帧。如果你的试验机加载速率提升到每秒0.1毫米以上破坏过程可能只有几十毫秒此时帧率至少要提升到2000帧每秒以上。一般来说宁可把帧率设置偏高一点也不要在破坏阶段因为样本帧不足而后悔。4.3 同步触发让力学数据和光学数据对齐载荷、位移数据记录频率通常在几十到几百赫兹高速相机在几百帧每秒拍摄两套系统不在同一个时钟上数据对齐就成了关键环节。我的方法是采用统一的外部触发源。试验机的数据采集系统每记录一个数据点就输出一个脉冲信号这个信号同时触发高速相机抓拍一帧。这样每一帧图像都与一个力学数据点严格对应。实际做下来这个同步方案在时间轴上可以做到亚毫秒级别的对齐精度。如果没有硬件同步条件也可以用软件同步在相机画面里放一个LED指示灯试验开始时通过信号发生器点亮LED同时触发加载和数据采集后期通过比对LED亮起的时间点来对齐数据。这个方法精度稍低但可作为备用方案。5. 数据处理流程与损伤演化指标提取5.1 DIC计算参数调优原始图像采集完成后进入后处理阶段。DIC计算参数设置是否合理直接影响位移场和应变场的质量。最重要的参数是子区尺寸。子区是相关性匹配时使用的一个局部图像区域需要包含足够多的散斑特征才能被正确追踪。如果子区太小包含的散斑特征太少匹配结果容易受噪声干扰如果子区太大空间分辨率下降应变集中区会被平滑掉。我用的规则是子区边长覆盖3到5个散斑颗粒。以50毫米视场宽度和散斑颗粒1毫米估算子区边长大约设在15到25像素之间。步长决定位移场计算的网格密度。步长等于子区尺寸时就是“无重叠计算”步长小于子区尺寸则有重叠。我通常把步长设为子区尺寸的1/4到1/2在保证计算效率的同时提高位移场空间分辨率。应变计算需要设立应变窗口。应变窗口是位移场中用来计算应变的空间范围窗口太大则应变集中被过度平滑窗口太小则应变场容易出现“椒盐式”噪声。我一般从应变窗口为子区尺寸的两到三倍开始试观察噪声水平如果噪声明显则适当增大窗口如果应变局部化带变得模糊则适当减小。还有几个细节值得提。相关计算一般选零均值归一化最小平方距离函数它对光照变化有一定的容忍度。迭代格式选择双三次样条插值对位移场做亚像素精度计算。图像滤波方面我不建议在计算前做大幅平滑处理DIC计算本身就有低通滤波的效果。5.2 损伤定量参数怎么选从DIC数据里提取损伤演化指标常用的有这么几个最大主应变场是最直观的指标。岩石内部微裂纹扩展到表面时表面会产生强烈的应变集中对应最大主应变值显著升高的局部区域。通过设定阈值可以把应变集中区从背景应变场中分离出来识别出所谓的高应变带。等效应变场可以综合反映面内拉伸、压缩和剪切变形。在岩石压缩破坏中损伤区往往同时包含多种变形模式单一的应变分量不能全面描述。等效应变作为一个标量综合指标在识别损伤区范围时比单向应变更可靠。离面位移场是三维DIC独有的优势数据。岩石在压缩过程中表面会向外鼓起这个鼓起量和位置直接反映了内部损伤的扩展深度。我处理数据时会单独画离面位移云图观察它和表面应变集中区是否一致。在很多试件中离面位移集中的位置与宏观裂纹出现的位置高度吻合甚至比应变场更早出现异常。在具体定义损伤变量时我用的是高应变区域面积占比设定一个应变阈值通常取试件破坏时应变幅值的某个比例计算超过阈值的像素面积占总面积的百分比。这个指标随加载过程变化明显在早期基本接近零到了临近破坏时急剧上升很适合做损伤演化曲线。也可以用某个代表性区域的平均应变值作为区域损伤指标但在非均匀变形场中平均值容易掩盖局部集中不如面积占比灵敏。5.3 从应变场里读破坏前兆我整理多组数据后发现不同循环次数下的凝灰岩破坏模式有明显的规律性变化。未经干湿循环的试件应变场相对均匀破坏前应变集中区主要出现在试件中上部然后快速形成一条贯穿型裂缝整个破坏过程非常快从应变集中到完全失稳基本不到一秒。经历干湿循环后试件表面在加载前就已经存在一定的初始损伤表现为应变场局部有不均匀的分布。加载过程中这些初始损伤区域会率先发展为应变集中区。循环次数越多应变集中区出现的时间越早分布范围越广破坏时往往不是单条裂缝而是多条裂缝同时扩展。利用逐帧图像做破坏前兆分析时我会重点关注一个特征应变局部化带是否开始贯通。这个阶段通常对应应力-应变曲线上的应力降开始之前。我发现干湿循环试件的应变局部化带贯通时间点往往早于应力峰值对应的时刻说明可以从DIC时序数据中提前识别出破坏前兆。6. 试验中的高频问题与避坑清单6.1 散斑脱落与试件表面崩落我做这类试验不下几十次踩过的坑比较有代表性值得集中列一下。散斑脱落是干湿循环试件做DIC最经典的问题。前面提到的“先循环后散斑”是第一个解决办法。但加载过程中当试件临近破坏时表面局部会崩落小块连带着散斑层一起飞掉。这个现象基本无法完全避免能做的就是尽量把加载速率放慢让破坏前兆充分发展避免试件积蓄过多弹性能后突然猛烈崩落。如果加装围压或者使用柔性加载垫片也能减少端部碎裂。6.2 光照变化与频闪干扰实验室的灯光并非完全恒定。高压钠灯和部分荧光灯存在工频闪烁肉眼看不出来但高速相机一拍就是几百帧这些闪烁会导致图像亮度周期变化DIC相关系数跟着周期性波动最终应变场出现规律的条纹噪声。解决方法是使用直流驱动的LED光源并在布置时确认高速相机画面中没有可见的亮度波动。更隐蔽的问题是环境光变化。试验期间如果有人在实验室走动或开关门自然光变化会导致图像亮度缓慢漂移。虽然DIC算法对线性光照变化有一定鲁棒性但严重时仍然影响匹配精度。我把整个试验区的窗帘拉上试验期间禁止人员进出光源固定后不再移动。6.3 烘干收缩、端部摩擦与边界效应干湿循环处理后的试件在加载时存在几个干扰因素。一是烘干后试件体积收缩可能导致其直径略小于压头压板尺寸造成端部接触不均匀。二是端部摩擦产生的环箍效应会让试件两端的面外位移受到约束导致计算得到的应变场在端部区域偏低形成人为边界层效应。处理办法加载前检查试件端面与压板贴合情况必要时在端面垫一层薄橡胶垫片或使用带球座的压板自适应调节。分析数据时端部区域上下高度方向各10%的范围的数据只作参考不纳入损伤定量统计。6.4 数据存储与后处理算力高速相机数据量相当可观。500帧每秒、双相机同时采集一分钟就能产生几十GB原始图像。如果不提前规划存储空间和后处理流程试验做一半就可能跑满硬盘。我的做法是给高速相机配万兆网卡连接到独立存储阵列试验前确认剩余空间至少能容纳两倍预估数据量。后处理阶段DIC计算在高性能工作站上仍然耗时一块几十GB的数据要算几个小时。建议根据试验阶段分批处理先处理低速阶段的少量帧确认整体趋势再处理破坏阶段的密集帧获取细节。这里踩过的坑是刚开始试验时全序列处理等了两天发现参数设置错了白白浪费了大量时间。6.5 试件分组与数据离散性控制天然凝灰岩的离散性比人造材料大得多。即使是同一块岩心钻取的相邻试件力学指标也可能相差一成以上。为了控制离散性我在分组时采用了“层位配对”的方法将取芯位置相近的试件平均分配到不同循环次数组尽量保证每组试件的初始质量相近。数据处理时每个循环梯度至少用三个有效试件最终结果取平均并计算标准偏差。如果某个试件的破坏模式明显偏离组内其他试件比如出现了不必要的端部劈裂我会标记并考虑是否剔除。这里没有绝对的对错标准但试验记录中必须写明避免后期数据解读时被异常点误导。7. 试验完成的后续扩展思路如果条件允许这套试验还可以往两个方向扩展。一是把声发射系统和高速3D-DIC同步使用声发射给出内部微破裂的时间和能量信息DIC给出表面位移场两者结合可以构建“内部损伤-表面响应”的对应关系。二是在加载前对不同循环状态的试件做CT扫描获得内部初始损伤的三维分布加载后再做一次CT对比两次扫描结果就可以直接看到内部裂纹的三维扩展路径弥补DIC只能测表面的局限。如果做数值模拟DIC数据可以作为很好的验证基准。把干湿循环后的试件初始损伤场导入数值模型对照DIC测得的表面应变场演化过程进行参数标定模拟的精度会明显高于只用宏观应力-应变曲线校准的模型。最后再分享一个一直沿用的小技巧每次试验前用同一组标定参数对标准刚体比如一个金属方块做一次DIC测量验证系统精度。如果刚体表面算出的应变量级超过100个微应变说明系统标定或参数设置有潜在问题需要排查后再开始正式试验。我在实际操作中深刻体会到这类检查虽然多花半小时但能拦住很多后面才暴露的大问题。测凝灰岩的损伤演化数据质量比数据量重要得多宁可慢一点、稳一点也不要在后期的处理堆里重新返工。
