夜雨聆风层状岩体真三轴加卸载：PFC3D把"层理面"那点事儿掰开了给你看
地下工程开挖，岩体不是被简单压着——而是被"卸"着。真实工况下，中间主应力和最小主应力压根不相等，岩体还带着层理面。这种复杂受力状态，常规三轴试验根本模拟不了。
华东交通大学熊良宵团队用PFC3D干了一件事：把层状岩体在真三轴加卸载条件下的力学行为，从微观颗粒层面完整复现了出来。
先搭模型——两种加载方式，七种层理角度
试样尺寸50×50×100mm（长宽高），在PFC3D里用离散颗粒+平行黏结模型构建。层理面用光滑节理（Smooth Joint）模拟，倾角β从0°到90°，每隔15°取一个工况，共7种。
两种加载方式很关键：A类——σ₁垂直于层理面方向加载；B类——σ₁平行于层理面方向加载。这决定了主应力与层理面的几何关系，直接控制破坏模式。
图：A类（左）和B类（右）加载方式下，主应力方向与层理面倾角β的几何关系示意图。A类加载时σ₁垂直于层理面，B类加载时σ₁平行于层理面。
参数标定——拿单轴压缩曲线"对表"
颗粒微观参数这样取：E_c=6.9GPa，k_n/k_s=1.2，摩擦系数μ=0.5，密度2600kg/m³。平行黏结抗拉、抗压强度均设为42±2MPa。光滑节理参数：k_n=700GPa/m，k_s=200GPa/m，节理摩擦系数0.8，黏聚力C=3MPa，内摩擦角45°。
拿Tien等[12]的室内单轴压缩试验做对标——完整试样模拟出来抗压强度49.2MPa（试验48.1MPa），弹性模量8.6GPa（试验8.9GPa），泊松比0.168（试验0.170），几乎吻合。含节理面的试样也一一对标，45°倾角偏差最大，但整体趋势完全一致。参数可信，往下走。
常规三轴加载——层理倾角一变，力学参数全跟着跳
固定围压，只加σ₁到破坏。结果很直观：
抗压强度随β增大先降后升，45°附近最低——此时节理面最利于滑移，岩体"沿层面"破坏，强度直接腰斩。弹性模量也是同一规律，45°最小。泊松比反过来，45°最大——侧向膨胀最厉害。峰值应变60°最小，意味着此时岩体最"脆"。
围压越大，各向异性越弱。3MPa围压下，β=60°时抗压强度只有约15MPa；围压升到10MPa，各倾角下的强度差距明显缩小。
图：常规三轴加载条件下，抗压强度、弹性模量、泊松比和峰值应变随节理面倾角β的变化规律。围压从0到10MPa，清晰展示各向异性特征。
加载vs卸载——围压相同，强度差不多；但应变差得远
初始围压3MPa时加载到破坏，和"升轴压+卸围压"到破坏——两种路径的抗压强度几乎一样。但卸载路径的峰值应变明显更小，泊松比更大。
换句话说：卸围压让岩体"提前"侧向鼓出，但抗压能力没怎么掉。真到破坏那一刻，加载和卸载的强度殊途同归。
图：初始围压3MPa和8MPa时，常规三轴加载与卸载条件下抗压强度、弹性模量、泊松比和峰值应变的对比。卸载路径峰值应变更小、泊松比更大。
真三轴加载——中间主应力一上，规律变了
真三轴的核心差异：σ₂≠σ₃。A类加载下，固定σ₂=10MPa、σ₃=3MPa，β从0°到90°扫一遍——抗压强度依然先减后增，60°最低。但中间主应力对强度的影响不是万能的：β=45°、60°、75°时，σ₂从5MPa加到9MPa，抗压强度几乎不动。只有在0°、15°、30°、90°这些"非滑移有利"角度，σ₂才明显抬强度。
B类加载的规律类似，但整体强度比A类高——σ₁平行于层理面时，岩体更难沿节理滑移，自然更"扛压"。
图：A类（上）和B类（下）真三轴加载条件下，抗压强度随节理面倾角β的变化。不同σ₂、σ₃组合下，60°附近均为强度低谷。
弹性模量和泊松比的规律跟常规三轴一致：45°弹性模量最低，泊松比最高。峰值应变60°最小——这个"最脆角度"很稳定，不受σ₂、σ₃组合影响。
真三轴加卸载对比——大多数情况下，强度还是差不多
破坏时最小主应力相同、σ₂相同，β取0°、15°、30°、45°、90°时，加载和卸载的抗压强度基本重合。只有β=60°和75°时，A类加载和卸载才有明显差异。
峰值应变依然敏感——加载路径的峰值应变普遍大于卸载路径。
图：A类（左）和B类（右）真三轴加载与卸载条件下，抗压强度随破坏时最小主应力的变化。大多数倾角下两条曲线几乎重叠。
一句话总结
层状岩体的力学行为，核心变量就两个：层理面倾角β，以及主应力与层理面的几何关系。45°~60°是最危险区间——强度最低、变形最大、各向异性最显著。开挖卸荷时，虽然最终抗压强度和加载路径差不多，但峰值应变和泊松比的差异提醒我们：地下工程的变形控制，不能只看强度这一个指标。
用PFC3D的好处就在这——它让你看清每个颗粒怎么转、每条节理怎么滑，而不只是一根应力-应变曲线。
建模到报告的全流程服务，包售后优化
