夜雨聆风PFC FLAC3D 代做，土木工程、水利、岩土、地下工程等各类|||颗粒流程序调试，经验丰富，专注研究对象：土、岩石、土石混合体、混凝土等颗粒材料。
提供服务：
- 单轴抗压/抗拉、巴西劈裂、双轴、三轴等各类静力/动力实验
- 刚性/柔性挡土墙、盾构隧道等工程问题
- 各类前/后处理技术，成样、参数标定等
PFC颗粒流 离散元 耦合 静力触探 ，pfc颗粒流 常规三轴 离散元 模拟，损伤 dfn 损伤效果 裂纹数量 应力应变 pfc3d，pfc颗粒流（PFC5.0 ）参数标定，组合体，单轴，双轴，冲击，pfc 颗粒流三轴单轴建模 可套袋可询 单轴、三轴、直剪、边坡、耦合。
离散元方法仿真，多相流，爆炸冲击lsdyn，pfc颗粒流离散元分析
干湿交替循环和循环荷载作用下边坡稳定性研究
在颗粒流程序中，接触模型起着至关重要的作用，其主要作用是描述和模拟颗粒之间的相互作用力。主要分为接触刚度模型，接触滑动模型和粘结模型。
粘结模型（Bonding Model）设定了颗粒法向力和剪切力的最大合力值。该模型有两种形式：接触粘结模型（Contact-Bond Model）和平行粘结模型（Parallel-Bond Model）。平行粘结模型广泛应用于岩石及其他颗粒材料的力学分析中，该模型假设颗粒之间通过一个平行于接触面的粘结层进行连接，从而模拟了颗粒间的粘结作用，颗粒通过接触力和粘结力相互作用，当颗粒之间的受力超过粘结力时，粘结层会断裂，导致颗粒之间失去粘结关系，从而模拟了材料的破裂和失效过程。在平行粘结模型中，粘结强度、接触半径、弹性模量和摩擦系数等参数共同决定了颗粒间的粘结行为，它够模拟岩石或土体在不同应力状态下的破坏过程和变形行为，能够揭示颗粒材料的力学行为及其破坏机理。平行粘结接触模型及其力学行为示意图如图 5-2 所示。
计算模型建立及参数标定
通过 PFC2D 构建模型边界，确保颗粒在运行过程中不会超出边界范围，在其中生成 6516 个随机分布的颗粒，颗粒密度为 2700kg/m3，消除模型建立过程中产生的悬浮颗粒后，最终得到岩石的计算模型，如图 5-3 所示。
根据文献[114]的研究，在 PFC 的模拟过程中，平行粘结接触模型在模拟岩石材料损伤过程中的力学行为表现出较好的效果，故本文采用平行粘结接触模型作为颗粒之间的细观本构模型，模拟颗粒之间的粘结作用，通过对建立的计算模型施加应力腐蚀模型，模拟岩石在不同围压条件下的应力状态。表 5-1 为平行粘结模型的参数。图 5-3 为PFC2D 计算模型。
对于细观参数的标定，本文采用“试错法”，先初步选择一组细观参数，通过不断调整细观参数改变 PFC 数值试验的曲线，并与实际试验曲线进行对比，直至数值试验曲线与力学试验曲线的差异较小[115]，即为完成细观参数标定，再以同样的方法，对比不同围压条件下的力学试验曲线，获得相应的细观参数。
图5-5 为常规压缩中数值试验曲线和室内试验曲线的对比。由于 PFC 计算过程中，颗粒之间的接触方式于实际环境有所不同，压密阶段在初始静力平衡阶段就已经完成，因此，参照王震[63]的对比方法，在不考虑压密段的情况下，将 PFC 中得到的数值试验曲线平移后与实际力学试验曲线对比，由图 5-4 中平移后 PFC 试验曲线可以看出，细观参数模拟得到的岩样常规压缩应力-应变曲线与室内试验差距较小，最终确定了数值模型的细观参数。
利用标定的细观参数，通过模型计算得到不同围压条件下白云岩的循环加卸载应力-应变曲线，图 5-5 为白云岩在循环加卸载数值试验的应力-应变曲线，表 5-2 为循环加卸载数值试验与室内试验峰值抗压强度的对比。
通过表 5-2 可以发现，不同围压条件下，数值试验与室内试验峰值抗压强度的差值分别为 3.85%、1.37%、0.26%、1.79%和 2.28%，均不超过 4%，表明数值试验中所标定的细观参数具有合理性。在 PFC 模型中，岩石力学特性由颗粒间的粘结决定，水对白云岩颗粒间的水-岩作用是导致其力学性能下降的重要原因，通过“试错法”得到的细观参数，能够反应白云岩在不同干湿交替循环次数条件下的力学特性，具有合理性。
工程背景
以云南磷化集团某露天矿为案例，采用 PFC2D 对经历不同干湿交替循环次数的边坡进行数值模拟，将标定的细观参数应用到边坡模型中，采用改进的重力放大法，探究边坡变形发展的破坏模式、位移特征以及稳定性系数变化。
数值模型的建立
假设工程体上方岩层均为白云岩，细观参数采用上文中标定得到的平行粘结模型细观参数。为了减少计算量，将计算模型简化为二维的边坡模型，图 5-7 为边坡边界及PFC2D 边坡模型，模型具体建立过程如下。
（1）生成颗粒
定义边坡模型的底部长为 50 m，高为 35 m，将边坡的外轮廓确定为墙体，生成颗粒后，模型共有 21647 个颗粒，采用两种大小不同的颗粒；削坡后大约有 14900 个颗粒，坡高 21 m，坡顶为 19.86 m，边坡倾角为 70°。
（2）细观参数的确定
通过 PFC2D 对降雨作用下边坡的渐进破坏进行了研究，接触模型选择了平行粘结模型，同时取得了良好的效果。因此，本节同样采用该模型，并以室内试验标定的细观参数为基础，应用到边坡模型的动力学计算当中。PFC2D 进行常规压缩数值模拟时，因不受重力影响，颗粒的密度对试验结果没有影响。但是，模拟边坡宏观破坏时，由于重力作用，经历过干湿交替循环作用的白云岩受到的重力大于干燥状态下的白云岩，必须考虑颗粒密度的影响，根据第二章基础物理参数测得的干燥与饱和白云岩的质量差，按比例增加干湿交替循环作用后边坡模型的颗粒密度。具体的细观参数值见表5-3，表 5-4 为不同干湿交替循环次数下边坡模型的颗粒密度。
（3）设定边界条件与地应力平衡
边坡模型的边界条件为边坡两侧以及下部的位移为零，边坡颗粒仅在它们自身重量作用下达到平衡（在模型中施加重力引起的地应力），将初始重力加速度设置为 9.8 m/s2。
数值模拟结果分析
通过 PFC2D 模拟得到了在重力放大法下，不同干湿交替循环次数及不同围压对边坡的影响。选取了具有代表性的渐进破坏模式进行分析，图 5-8 展示了边坡渐进破坏模式。总体来看，边坡的渐进破坏模式可以分为以下几个阶段：
（1）裂隙发育阶段
在重力加载的前期，由于边坡受到重力作用，裂隙尖端出现应力集中，引起裂隙的萌发与发育，随着裂隙的扩展，裂隙不仅会向坡面延伸，还会向坡顶扩展。在这一阶段，整个坡体除了坡脚基本没有位移，但由于坡脚区域的应力集中，出现了少量垮塌。
（2）坡脚滑动阶段
随着重力的不断增加，坡体的应力逐渐集中于坡脚区域，当应力超过临界应力时，坡脚处的裂隙和原生裂隙贯通使得产生的块体发生滑动并向下滑出。
（3）坡体滑落阶段
在坡脚处发生滑移后，该区域裂隙发育严重，进而开始产生自下而上的渐进式破坏。上部岩体受到重力作用，裂隙迅速发育并向外扩展，块体呈破碎状向下滑动，最终导致坡体的滑落，造成边坡失稳。
（4）坡顶坍塌阶段
边坡后缘拉裂破坏，坡顶受重力作用导致下滑力高于抗滑力，坍塌现象进一步加剧。
边坡位移特征分析
围压 1 MPa 及不同干湿交替循环次数条件下的边坡位移云图如图 5-9 所示，可以看出不同干湿交替循环次数对边坡的位移有较大的影响，对于干燥状态及干湿交替循环次数较少（0 次和 1 次）的边坡，在重力加载的作用下边坡的整体位移相对较小，裂隙发育较少，仅在原生裂隙尖端出现少量裂隙。随着重力作用的增大，裂隙的发育开始加速，边坡的稳定性逐渐降低。在边坡的潜在滑动面与裂隙的共同作用下，坡脚处裂纹尖端出现应力集中，最终该区域的应力超过临界值，从而导致边坡失稳现象。
对于经历过较多次数干湿交替循环作用（5 次、10 次和 30 次）的边坡，干湿交替循环通过水-岩的物理、化学作用加速了裂隙的发育与扩展，坡体内部着随干湿交替循环次数的增加，强度急剧降低，裂隙发育情况严重，几乎遍布整个坡面，干湿交替循环作用加速了裂隙的发育进程，也加大了边坡失稳的风险。
图 5-10 展示了 30 次干湿交替循环及不同围压条件下的边坡位移云图。可以看到，尽管围压条件有所不同，但无论围压的高低，边坡破坏面积都较为广泛，颗粒滑动距离较大，表明干湿交替循环作用对颗粒之间的粘结能力有明显的削弱。由于重力的增加导致坡脚区域承受较大的应力和剪切力，颗粒之间粘结能力减弱，无法有效抵抗重力作用，裂隙不断扩展，坡脚处首先发生破坏，开始产生位移，并发展成自下而上的渐进式破坏。
围压的增加虽然可以提高颗粒之间的粘结，但其作用无法抵消干湿交替循环导致的劣化作用。因此，不同围压条件下边坡的破坏模式具有相似的特征。
稳定性系数分析
边坡稳定性系数可以直接的反映出边坡的稳定性，采用改进的重力放大法对边坡稳定性系数进行计算，并将计算结果绘制成图 5-11，不同干湿交替循环及不同围压的边坡稳定性系数如图 5-11 所示。
观察图 5-11，可以看出，围压 0 MPa 及不同干湿交替循环次数条件下，稳定性系数由 1.18 降至 1.04，围压 1 MPa 及不同干湿交替循环次数条件下，稳定性系数由 1.16 降至 1.05，围压 2 MPa 及不同干湿交替循环次数条件下，稳定性系数由 1.17 降至 1.07，围压 3 MPa 及不同干湿交替循环次数条件下，稳定性系数由 1.17 降至 1.08，围压 5 MPa及不同干湿交替循环次数条件下，稳定性系数由 1.18 降至 1.08，干燥条件下的边坡稳定性系数高于经历过干湿交替交替循环的边坡稳定性系数，表明干燥条件下的边坡是相对稳定的。随着干湿交替循环次数的增加，边坡稳定性系数均会降低，且降低速率减小，表明颗粒的粘结能力不断下降，干湿交替循环对边坡的累积劣化作用逐渐减缓。
对比 0 次、1 次、5 次、10 次和 30 次干湿交替循环及不同围压条件下边坡的稳定性系数，可以看出，稳定性系数分别由 1.15、1.11、1.08、1.05 和 1.04 增加至 1.18、1.16、1.13、1.11 和 1.08，这表明当白云岩的干湿交替循环次数条件一定时，围压的增加提高了颗粒的粘结能力，边坡体内部的裂隙数量减少，边坡的结构更加稳定，因此，边坡稳定性系数随围压的增加而增大。
（1）通过 PFC2D 得到的数值模拟结果与室内力学试验的对比误差较小，所标定的细观参数具有合理性，以及采用的应力腐蚀模型具有较高的准确度。
（2）岩质边坡的失稳是自下而上依次贯通的渐进破坏过程，按照边坡破坏模式可以分为四个阶段：裂隙发育阶段、坡脚滑动阶段、坡体滑落阶段和坡顶坍塌阶段。
（3）通过改进的重力放大法，计算了不同条件下的边坡稳定性系数，数据显示围压1 MPa 及不同干湿交替循环次数条件下，稳定性系数由 1.16 降至 1.05，30 次干湿交替循环及不同围压条件下，稳定性系数由 1.04 增至 1.08，结果表明围压越大，边坡稳定性系数越高，结构更加稳定，边坡稳定性系数随干湿交替循环次数的增加而降低，呈现负相关性。
