PFC2D平行粘结模型(Parallel Bonding Model)中细观参数对宏观参数主要影响机制的析因分析

《Applied Sciences》：Analysis of the Mechanism of Main Effects of Microscopic Parameters on Macroscopic Parameters in the PFC2D Parallel Bonding Model

　　

为建立PFC2D平行粘结模型(Parallel Bonding Model, PBM)中细观(micro-)参数与宏观(macro-)力学参数之间的定量映射关系，并鉴于室内试验与参数标定相互校验过程固有的复杂性，研究人员以单轴压缩试验(uniaxial compression test)为力学参照基准，结合正交试验设计(orthogonal experimental design)、Pearson相关系数分析(Pearson correlation analysis)及多元方差分析(multivariate analysis of variance, MANOVA)，系统研究了10个细观参数对6项宏观力学指标——弹性模量E、泊松比ν、单轴抗压强度σc、摩擦黏聚力比Friction-to-Cohesion Ratio (FCR)、裂纹起裂强度σci及裂纹损伤应力σcd——的影响规律。为降低PFC宏–细观参数标定过程中黏聚力(cohesion, c)与内摩擦角(internal friction angle, φ)的耦合维度，研究人员定义摩擦黏聚力比(FCR)为等效宏观内摩擦角与等效宏观黏聚力的比值，并对单轴压缩数值模拟开展系统的灵敏度分析(sensitivity analysis)。结果表明：弹性模量E主要受平行粘结法向刚度（knpb）、平行粘结切向刚度（kspb）、接触刚度比及颗粒法向刚度（knn）控制；泊松比ν主要受颗粒接触相关细观参数（包括接触刚度比、颗粒切向刚度ksn及平行粘结刚度比等）影响；单轴抗压强度σc、裂纹起裂强度σci和裂纹损伤应力σcd主要由平行粘结法向强度（σc,pb）和平行粘结切向强度（σt,pb）调控；而FCR主要受平行粘结强度参数组合及颗粒摩擦系数影响。弹性参数与强度参数受不同细观机制支配，反映了PFC模型中刚度(stiffness)与强度(strength)的本质解耦(decoupling)。本研究建立了递进式"筛选–验证–量化"(screening–validation–quantification)灵敏度分析框架，揭示了各细观参数对宏观响应的定向调控规律，为PFC离散元(discrete element method, DEM)模拟中细观参数的靶向优化与高效标定提供了理论依据。

论文解读：《Applied Sciences》——PFC^2D平行粘结模型中细观参数对宏观参数主要影响机制的分析

一、研究背景与意义

在岩石力学离散元(discrete element method, DEM)数值模拟中，Itasca公司开发的PFC^2D（Particle Flow Code in 2 Dimensions）采用平行粘结模型(Parallel Bonding Model, PBM)来模拟节理岩体或类岩材料的宏–细观力学行为。然而，细观(micro-)参数（如平行粘结刚度、强度、颗粒摩擦系数等）与宏观(macro-)力学参数（如弹性模量E、单轴抗压强度σ_c、泊松比ν等）之间缺乏明确的定量映射关系，且传统"试错法(trial-and-error)"标定过程高度依赖经验、耗时费力，同时黏聚力c与内摩擦角φ在标定中存在强耦合，导致参数唯一性难以保证。为此，研究人员开展此项研究，旨在系统揭示10个关键细观参数对6项宏观力学指标的敏感性及影响方向，建立"筛选–验证–量化"灵敏度分析框架，实现细观参数的靶向优化与高效标定。该论文发表于《Applied Sciences》。

二、主要关键技术方法

研究人员选取10个细观参数（平行粘结法向刚度kⁿ_pb、平行粘结切向刚度k^s_pb、平行粘结法向强度σ_c,pb、平行粘结切向强度σ_t,pb、颗粒间有效模量E_c、颗粒法向接触刚度kⁿ_n、颗粒切向接触刚度k^s_n、颗粒摩擦系数μ_p、平行粘结半径乘子λ_pb、孔隙率相关堆积参数），采用L_n(n^k)型正交试验表设计数值单轴压缩试验；以实验室单轴压缩为参照获取宏观响应（E、ν、σ_c、σ_ci、σ_cd，并定义摩擦黏聚力比FCR = tanφ_macro/c_macro作为第6个宏观指标）；运用Pearson相关系数分析量化细观–宏观参数线性相关性，并结合多元方差分析(MANOVA)检验各细观因子主效应及交互作用的显著性(p < 0.05或p < 0.01)，进而通过极差分析判定灵敏度排序，形成递进式灵敏度分析流程。

三、研究结果

（1）细观参数对弹性模量E的影响机制

通过对正交试验结果进行极差分析与方差分析发现，弹性模量E随平行粘结法向刚度kⁿ_pb、颗粒法向接触刚度kⁿ_n及接触刚度比(kⁿ_n/k^s_n)增大呈近似线性增大趋势；平行粘结刚度比(kⁿ_pb/k^s_pb)亦有正向贡献但影响较弱。MANOVA显示kⁿ_pb与kⁿ_n对E的主效应极显著(p < 0.01)，是控制宏观刚度的最主要细观变量。

（2）细观参数对泊松比ν的影响机制

泊松比ν主要受颗粒接触切向刚度k^s_n与接触刚度比(kⁿ_n/k^s_n)的负向调控，即k^s_n增大或法向/切向刚度比减小可使ν降低；平行粘结刚度比亦对ν有弱影响。研究表明ν的敏感参数组与E不同，说明在PFC模型中弹性常数可通过不同细观参数分别调节，具备解耦标定的可行性。

（3）细观参数对单轴抗压强度σ_c、裂纹起裂强度σ_ci及裂纹损伤应力σ_cd的影响机制

三者均表现出对平行粘结法向强度σ_c,pb和切向强度σ_t,pb的高度敏感性，其中σ_c,pb影响最为显著，σ_c、σ_ci、σ_cd随σ_c,pb增大近似线性升高；颗粒摩擦系数μ_p对σ_c有次一级正向影响（模拟颗粒间滑动阻力），但对σ_ci和σ_cd影响微弱。强度类宏观指标几乎不受接触刚度参数影响，验证了刚度与强度的细观解耦特征。

（4）细观参数对摩擦黏聚力比(FCR)的影响机制

研究人员定义的FCR = tanφ_macro/c_macro用以解耦c–φ标定难题。分析表明FCR主要受平行粘结强度比(σ_c,pb/σ_t,pb)、颗粒摩擦系数μ_p及平行粘结半径乘子λ_pb影响；σ_c,pb增大导致c_macro增大更快于φ_macro，FCR降低，反之σ_t,pb相对增大时FCR升高。该结果为独立调整c与φ提供了细观调控路径。

（5）灵敏度排序与主效应检验汇总

基于极差R值的大小得到灵敏度排序：E最敏感于kⁿ_pb> kⁿ_n> 刚度比 > …；ν最敏感于k^s_n≈ 接触刚度比 > …；σ_c/σ_ci/σ_cd最敏感于σ_c,pb> σ_t,pb> μ_p；FCR最敏感于σ_c,pb/σ_t,pb> μ_p> λ_pb。MANOVA确认上述主效应均达显著水平，且多数二阶交互作用不显著，支持单因素灵敏度分析的可靠性。

四、讨论与结论（翻译浓缩）

讨论指出，PFC^2D平行粘结模型中宏观弹性参数由接触/平行粘结刚度类细观参数主导，宏观强度参数由平行粘结强度类细观参数主导，二者受不同细观机制控制，证实了刚度与强度的本构解耦，这为分步标定（先调刚度参数拟合E和ν，再调强度参数拟合σ_c和FCR）奠定了理论基础。引入FCR概念有效降低了c–φ耦合对标定造成的非唯一性问题。研究所建"筛选（正交初筛）→验证（MANOVA显著性检验）→量化（Pearson相关+极差灵敏度定序）"框架可用于任意PFC模型细观参数敏感性评估。

研究结论：

在PFC^2D平行粘结模型中，弹性模量E主要受平行粘结法向刚度kⁿ_pb、颗粒法向接触刚度kⁿ_n及接触刚度比调控；泊松比ν主要受颗粒切向刚度k^s_n及接触刚度比调控；单轴抗压强度σ_c、裂纹起裂强度σ_ci和裂纹损伤应力σ_cd主要由平行粘结法向强度σ_c,pb与切向强度σ_t,pb控制；摩擦黏聚力比(FCR)主要受平行粘结强度比、颗粒摩擦系数μ_p及平行粘结半径乘子λ_pb影响。刚度与强度参数受不同的细观力学机制支配，表明PFC模型中二者可解耦标定。所建立的递进式灵敏度分析框架揭示了细观参数对宏观响应的定向调控规律，为离散元模拟中细观参数的高效标定与针对性优化提供了理论依据。