颗粒材料双轴压缩试验的光弹测试

郑 虎， 牛文清， 毛无卫， 黄 雨

（1. 同济大学 土木工程学院，上海 200092；2. 同济大学 岩土及地下工程教育部重点实验室，上海 200092）

颗粒材料是由大量离散的固体颗粒所组成的宏观体系，也是自然界中存在最为广泛的物质之一（如砂土、砾石等），约有75%的原材料属于颗粒物质，其宏观力学行为主要是由颗粒的形状、材质、大小、排布等因素决定［1-2］。颗粒材料的离散性决定了其力学行为的复杂性，许多自然灾害现象和地质灾害的本质即为岩土颗粒材料变形破坏［3-5］。岩土颗粒材料的变形破坏是通过细观力链和应变局部化演化耦合作用后的宏观表现，因此，深入研究颗粒体系应变局部化演化特征［6］，揭示颗粒物质失稳的物理机制，对预防此类地质灾害的发生有着重要的意义。

颗粒物质的非仿射位移是表征颗粒材料发生剪切变形时颗粒尺度上各颗粒间相对位移的有效参量，非仿射位移场可以通过从连续位移场中减去单颗粒位移而得到［7-8］。非仿射位移与局部区域内变形相关［9］，被认为是导致局部不可逆塑性或剪切带产生的原因［10-11］。当外部荷载过大会导致非仿射位移累积成宏观尺度上的剪切带，剪切带演化的本质是细观尺度上的力链断裂与重建和颗粒尺度上颗粒空间位置发生非仿射位移后重新排布［12-13］。通过分析非仿射位移场来研究剪切带的演化特征，将有助于揭示岩土颗粒材料变形失稳的内在机制。

光弹试验对揭示颗粒系统内部细观尺度力链结构分布特征有着极大的优势［14-17］，可以获得剪切带产生过程中颗粒尺度上空间位置的分布、转动，细观尺度上的力链网络等数字图像信息［18-21］，对从细观尺度上揭示岩土体失稳过程中剪切带演化特征的研究有着重要的应用前景。Wang 等［22］通过光弹试验研究了颗粒系统剪切膨胀过程中形成阻塞状态的结构特征。Zheng 等［23］基于光弹测试技术研究了颗粒材料中产生剪切带时的细观尺度力链演化特征，发现颗粒系统局部密度变化和颗粒旋转与剪切带有很强的相关性。因此，基于光弹试验从细观尺度上识别力链演化特征对于理解颗粒材料的破坏力学行为起到了关键作用。

本文将基于光弹测试原理对二维光弹颗粒系统进行双轴压缩试验，结合粒子跟踪测速 （PTV）技术对试样中各个颗粒在时间和空间上进行连续追踪，得到颗粒位移的位移场。通过分析颗粒体系细观力链演化规律和宏观位移场之间的关系，揭示颗粒材料压缩变形过程中细观尺度上非均质力链网络演化机制。

1 试验技术与试验装置

基于光弹测试技术进行颗粒材料相关力学性质的研究已在物理学软物质领域里应用比较成熟。光弹测试技术是基于偏振光穿透暂时（人工）光学各向异性体时发生双折射现象后所产生的相位差来计算光弹材料所受到力的大小，所谓的相位差表现出来的现象即为在颗粒材料表面出现由于频率相同的光波发生干涉后的明暗条纹［24-25］，结合光学应力定理［19］便可以定量地算出颗粒所有的应力张量。因此，可以通过提供相对均匀并稳定的偏振光场，再根据具体的试验需要设计出相应的试验方案。除此之外，还可以根据具体的研究情况研制出不同形状、不同强度的试验材料，同时颗粒与颗粒间的摩擦力也可以定量地标定。三轴试验是研究岩土体力学性质的经典且常规的试验方法。但目前采用光弹测试方法定量地分析三维的光弹试验数据还不是很成熟，已有的试验方法基本都是通过将试验简化到二维平面内进行研究分析［26］。

图1是通过结合三轴压缩试验方法设计的一个二维双轴压缩光弹试验装置示意图。光源提供稳定均匀的自然白光，由放置在光源上方的圆形偏振片将自然光转变成均匀且垂直于偏振片的圆偏振光，当圆偏振光穿过受力状态的光弹颗粒后，颗粒材料内部所有处于受力状态的颗粒均会产生双折射现象从而连接成网络结构，即为力链。进一步通过控制数码相机前端的另一个对应的圆偏振片实现同时捕捉光弹颗粒在自然光照和偏振光场中的图片信息。试样的两侧以两根细绳为边界，通过控制细绳的张力来控制试样的边界压力，并以分布式马达来实现对颗粒试样的压缩剪切。采用拟静态的加载方案，每压缩一步，ΔS=1 mm，对应的轴向应变约为0.25%。通过对偏振光场的控制，并在每两个压缩步之间通过数码相机对试样进行数据采集。

图1 二维双轴压缩光弹试验装置示意图Fig. 1 Biaxial setup based on photoelastic technique

图2中矩形区域为本次试验所用的光弹颗粒材料。光弹颗粒材料质地为环氧树脂，单个颗粒材料的弹性模量为4.8 MPa，泊松比为0.4～0.5，应变光学系数C≈0.002 MPa-1，组成试样的大小为200 mm×395 mm。为了保证试样更具有代表性，试验过程中采用两种不同大小的颗粒制样，以避免试样中出现晶状堆积的现象。两种颗粒的直径分别为8.6 、7.4 mm，厚度为6.4 mm，且以1： 3 的比例制样。

图2 颗粒样品的俯视图Fig. 2 Top view of granular samples

基于Hough 变换等数字图像处理手段，自主编制图像处理程序对所采集的数字图像进行分析，实现颗粒边界和位置的智能识别，识别的准确率可以达到100%，并通过数字图像粒子追踪技术在时间域内对单个颗粒进行追踪和编号定位。结合各颗粒的位置和其光弹效应，基于平均灰度值梯度法［19］计算各颗粒的受力信息，并重构颗粒系统细观尺度上的力链网络。

2 试验结果及分析

2.1 基于Hough变换的粒子位置追踪

颗粒位置追踪是基于相邻两个压缩步颗粒的位置变化特征来确定的，例如，在第n压缩步时某颗粒A的坐标为（xA(n)，yA(n)），那么在第n+1 压缩步中寻找距离A最近的颗粒，即

式中：Xi(n+1)和Yi(n+1)分别为所有颗粒在第n+1 压缩步时的横坐标与纵坐标。依此类推处理每一压缩步。这样便可对每个颗粒进行编号，在时间上进行颗粒追踪（particle tracking velocimetry，PTV），确保每个压缩步时同一个颗粒拥有唯一一个编号，从而实现对试样中各个颗粒在时间上和空间上的连续追踪。

图3 为试样在应变γ=25%时各颗粒的位移场。其中图3a为各颗粒实际位移的位移场，可以看出在试样底部颗粒几乎没有发生位移，位移的大小沿着压缩的方向逐渐减小。而颗粒沿着压缩方向的位移大小主要由颗粒离压缩端的位置决定，其位移场如图3b 所示。虽然试验过程中只是对试样施加了沿轴向的压缩，但颗粒在被挤压的过程中会发生与压缩方向垂直的剪胀，如图3c所示，在靠近压缩端和底部固定端位置的颗粒并没有发生（或者只是发生了相对很小）沿着垂直于压缩方向的位移。颗粒自试样中部开始向两侧发生剪胀，越靠近边界位置的颗粒位移越大。由图3c的构型可以看出，该构型与传统的土体以液压为围压的三轴试验压缩后的截面极为相似，从而也验证了该二维柔性边界试验的可靠性。

图3 试样在γ=25%时的位移场Fig. 3 Displacement field of granular sample at γ=25%

鉴于颗粒材料在固定端附近沿着压缩方向位移的差异性，从一个方向对试样进行压缩，分析细/微观尺度上颗粒的位移特性，往往只能直接得到以压缩端为主的位移场。因此，需要对颗粒位移进行重新归一化处理，进而分析试样在压缩后颗粒变形的本质。

2.2 颗粒位移场归一化处理

通过对颗粒材料试样的一端以应变控制进行压缩剪切，对试样沿着压缩方向的位移场进行归一化处理，其中第i个颗粒在压缩第n步时的位移可由下式进行计算：

式中：yi0和yin分别为试样中第i个颗粒在压缩端初始位置以及在压缩了第n步后的位置；δ和n分别为每压缩剪切步的压缩量以及压缩步数；ye和y0分别为试验开始时试样的固定端以及压缩端的位置。

试样归一化处理后的位移场如图4a 所示，可以看出，试样中心位置的颗粒并没有发生相对位移，在试样的固定端和压缩端分别形成了似锥形的压缩块体向试样内部延伸。而在试样的中部颗粒主要是以沿着与压缩方向垂直的位移为主。通过图4b 可以更清楚地看到类似锥形的压缩块体，而在试样中部沿压缩方向的颗粒位移几乎为零。

图4 试样在γ=25%时归一化后的位移场Fig. 4 Normalized displacement field of sample at γ=25%

颗粒材料位移的概率密度分布（probability density function， PDF）如图5所示，可以看出颗粒在沿压缩方向的位移分布基本均匀，最大位移量分别为在固定端和压缩端附近。而颗粒在垂直于压缩方向的位移主要集中在零附近，最大位移量与压缩方向相比较小，最大位移的颗粒数目也较少。

图5 试样中颗粒位移沿着压缩方向和垂直于压缩方向的概率分布Fig. 5 PDF of granular displacement along compression direction and perpendicular to compression direction

试样中颗粒总位移量的PDF 如图6 所示，颗粒的绝对位移大部分超过10 mm，只有很少一部分颗粒发生了小于3 mm 的位移，结合图4 可以看出，这部分颗粒主要集中在试样的中心位置。

图6 试样中颗粒总位移量的概率分布Fig. 6 PDF of granular total displacement

2.3 试样内部力链网络分布

基于光弹测试技术进行颗粒材料力学特性相关试验研究的最大优势是可以直观地观测到试样在压缩变形后内部应力的网络结构。力链的明暗程度代表着颗粒间接触力的强弱，当压缩量较小时颗粒材料内部会出现明显的应力诱导各向异性的力链网络结构［19］。

图7 为试样在γ=25%时的力链网络结构照片，颗粒材料中力链容易发生弯曲［27］，当存在侧限时会形成与压缩方向垂直的力链，由于试样此时的压缩量相对较大，可以看出试样的力链方向除了沿压缩方向外，在垂直于压缩方向也可以看到力链的分布。颗粒材料受力后的特点与固体材料或者液体不同，并不是每个颗粒都参与外部荷载的分担。从图7右局部放大图可以看出，若干个颗粒形成了环形的强力链网络，在环形力链中的颗粒几乎没有承担任何力。因此，即使移除这些颗粒对试样的应力平衡状态或颗粒试样内部力链网络结构影响也较小。试样的光弹效应照片除了能够定性地描述颗粒力链传递特征之外，还可以通过颗粒内部灰度值梯度平方的平均值＜G2＞，定性计算各个颗粒的压力大小，照片中每个像素点位置的灰度梯度平方∇I2i，j可以通过下式计算：

图7 试样在γ=25%时的力链网络结构Fig. 7 Force chain network at γ=25%

式中：Ii+1，j、Ii-1，j、Ii，j+1、Ii，j-1、Ii+1，j+1、Ii-1，j-1、Ii-1，j+1和Ii+1，j-1分别为待计算像素点相邻像素点的灰度值。而每个颗粒的平均灰度值平方可以由式（4）计算。

图8 为计算所得试样中各颗粒的＜G2＞分布图，图中圆盘代表颗粒的位置，其中白色的圆盘表示颗粒所受压力最大，黑色表示颗粒所受压力相对较小。黑色圆盘和白色圆盘的分布特征直接反映了颗粒材料该状态下的内部受力分布，力链网络的分布沿着纵向相对于水平方向较强。值得注意的是，材料的中心位置颗粒所受的压力并不是很大，但是可以清楚地看到该部分的颗粒分布相对均匀，颗粒的局部密度相对较小，此处的颗粒相对位移比较大。

图8 试样在 γ=25%时＜G2 ＞的分布Fig. 8 ＜G2 ＞ distribution of granular system at γ=25%

3 结论

通过自主设计的柔性边界试验装置对由两种不同粒径组成的光弹颗粒混合系统进行了应变控制式的双轴压缩试验。结合试验数据，利用Hough 变换等数字图像处理手段对每个颗粒进行定位，并通过PTV技术对各颗粒进行编号，实现了对每个颗粒在时间和空间上的连续定位。通过分析颗粒归一化位移场发现，在应变为25%时颗粒材料中心位置颗粒的相对位移非常小（小于3 mm），试样主要通过固定端和压缩端附近以及试样中部靠近边界处的颗粒移动实现变形，而在固定端和压缩端主要以沿压缩方向的位移为主，试样中部靠近变形处的颗粒主要以垂直与压缩方向的位移为主。颗粒中心处的力链也呈现出相对各向同性的性质，密度也相对较小，正是剪切带所穿越的地方。颗粒材料的光弹效应清晰地揭示了颗粒在压缩过程中局部颗粒几乎没有受力的现象。

作者贡献声明：

郑 虎：论文的选题、指导、修改。

牛文清：具体工作的开展和论文撰写。

毛无卫：论文的指导、修改。

黄 雨：论文的指导、修改。