PFC双轴压缩模拟：松散与密实砂样的力学行为对比

1. 项目概述：松散与密实砂样的双轴压缩模拟

在岩土工程和地质力学领域，理解颗粒材料的力学行为对边坡稳定、地基承载力和地震液化等实际问题至关重要。这个PFC案例通过离散元方法，系统研究了四种典型砂土配置在双轴压缩条件下的力学响应差异。作为从业十余年的岩土仿真工程师，我认为这个案例的价值在于它同时考虑了材料密实度和块石夹杂这两个关键因素——这正好对应了实际工程中常遇到的"原生结构"和"施工扰动"两种工况。

松散砂样代表自然沉积或人工回填的初始状态，而密实砂样则对应经过碾压或振动加密处理后的材料。块石的引入则模拟了地质构造中常见的夹杂物或填方工程中的粗骨料。二维双轴压缩试验是评估这类材料剪切特性的黄金标准，可以获取应力-应变关系、体积变化和破坏模式等核心参数。

2. 模型构建与参数设定

2.1 颗粒体系生成逻辑

创建四种试样时，需遵循严格的制备流程以保证结果可比性。对于松散砂样，采用低摩擦系数(μ≈0.1)的颗粒随机生成，随后仅施加重力完成自然堆积；密实砂样则需在重力沉积后施加各向同性压缩，直至孔隙比达到目标值(通常en<0.6)。块石通过生成特定粒径的"簇颗粒"实现，其尺寸一般为砂粒平均直径的5-10倍。

关键参数设置示例：

;砂粒基本参数 set disk_diameter 0.1mm set density 2650kg/m³ set kn 1e8N/m set ks 0.3kn ;块石特殊设置 cluster create id=1 radius_ratio=8.0 set cluster_kn 2e8N/m ;块石刚度提高

2.2 边界条件与加载方案

采用刚性墙模拟加载板，底部固定而顶部施加竖向位移荷载，侧向采用伺服控制保持恒定围压。建议采用阶梯式加载策略：

1. 各向同性固结阶段：同步增加四向压力至目标围压(如100kPa)
2. 剪切阶段：固定侧压，以0.1%/s的轴向应变速率压缩
3. 终止条件：轴向应变达到15%或试样完全破坏

重要提示：伺服控制算法需设置合理的阻尼系数(建议η=0.7)，避免系统振荡影响应力计算精度。

3. 关键结果分析与解读

3.1 应力-应变特性对比

通过监测边界墙受力可计算偏应力q=σ1-σ3，四种试样的典型响应表现为：

含块石试样表现出更显著的应变硬化特性，这是因为块石骨架形成了"桥接效应"，延迟了剪切带的形成。密实砂的强度峰值比松散砂高出约94%，印证了密度对抗剪强度的决定性影响。

3.2 体积变化规律

通过跟踪颗粒体系总体积变化，可观察到：

• 松散砂：持续剪缩，最大体积应变达-5%
• 密实砂：初始轻微剪缩后转为剪胀，最终体积应变+3%
• 含块石试样：体积变化幅度减小约30%，反映块石抑制了颗粒重排列

这种差异可通过配位数变化解释：密实砂初始配位数约4.5，剪切过程中升至5.2；而松散砂配位数始终低于3.8。

4. 微观力学机制探究

4.1 力链网络演化

使用PFC的contact force可视化功能，可以清晰看到：

1. 松散砂：力链稀疏且短暂，主要沿压缩方向排列
2. 密实砂：形成密集的网状力链，存在明显的力链分叉
3. 含块石试样：块石周围产生高强度力链集中区，形成局部"拱结构"

通过统计力链方向分布发现，密实砂在峰值强度时力链各向异性度达0.65，而松散砂仅0.4。

4.2 能量耗散分析

记录各类能量分量随时间变化：

• 弹性应变能：密实砂比松散砂高2-3倍
• 摩擦耗能：占总耗能70%以上
• 动能：仅在剪切带形成时出现短暂峰值

特别发现：含块石试样中，块石-砂粒接触面的摩擦耗能占比达40%，说明界面特性对整体行为影响显著。

5. 工程启示与参数优化建议

5.1 本构模型校准要点

根据模拟结果校准Mohr-Coulomb参数时需注意：

• 内摩擦角φ：密实砂取35-40°，松散砂28-32°
• 剪胀角ψ：密实砂取5-10°，松散砂设为0
• 粘结力c：含块石试样建议增加2-5kPa表观粘结力

5.2 数值建模实用技巧

1. 粒径分布控制：采用Fuller曲线分布可更好反映实际级配

   def fuller_distribution(d_min, d_max, n=0.5): return [d_min + (d_max-d_min)*(i/100)**n for i in range(101)]

2. 阻尼系数调整：剪切阶段建议降低至0.3以提高收敛速度
3. 数据采样频率：至少每0.1%应变保存一次完整快照

6. 常见问题解决方案

6.1 数值不稳定现象处理

6.2 计算结果验证方法

建议通过三重验证确保结果可靠性：

1. 能量平衡检验：总输入功≈弹性应变能+耗散能（误差<5%）
2. 与实验室数据对比：重点校核峰值强度和体积变化曲线
3. 网格敏感性分析：逐步增加颗粒数量直至结果收敛

在实际项目中，我们曾发现当块石含量超过30%时，二维模拟会高估强度约15%，这时需要考虑采用三维模型。另一个容易忽视的细节是颗粒形状效应——通过clump模拟非球形颗粒可使摩擦角提高2-3°。