(19)国家知识产权局 (12)发明 专利申请 (10)申请公布号 (43)申请公布日 (21)申请 号 20221081096 0.9 (22)申请日 2022.07.11 (71)申请人 贵州正业工程技术投资有限公司 地址 550000 贵州省贵阳市南明区双龙航 空港经济区机场路9号太升国际A2栋 14层1402- 34号 申请人 贵州正建 兴业工程质量检测有限公 司 (72)发明人 沈志平　刘欢　付君宜　孙洪　 许飞洲　张晋　杨欣　 (74)专利代理 机构 北京安瑞克专利代理事务所 (特殊普通 合伙) 11902 专利代理师 焦丽 (51)Int.Cl. G06F 30/17(2020.01)G06F 30/20(2020.01) G06F 119/14(2020.01) (54)发明名称 强夯碎石土PFC颗粒滚动阻力摩 擦系数标定 方法 (57)摘要 本发明提供了一种强夯碎石土PFC颗粒滚动 阻力摩擦系数标定方法， 所述方法包括： 获取强 夯后碎石土地基土体重度、 孔隙率、 粘聚力及内 摩擦角， 计算直立边坡的临界高度hcr； 利用筛分 法获取地基土样的颗粒级配， 基于获取的颗粒级 配、 重度及孔隙率， 利用PFC2D颗粒流软件创建强 夯碎石土填方地基模型； 在保持其他细观参数不 变的情况下， 通过调整所述滚动阻力摩擦系数的 大小， 得到不同的地基模型； 以所述直立边坡临 界高度hcr为开挖高度， 依据所得的不同地基模 型进行直立开挖， 开挖坑壁保持稳定而不坍塌时 所使用的地基模型对应的滚动阻力摩擦系数即 为所确定的与实际碎石土地基相适应的细观参 数。 本发明为PFC模拟强夯碎石土地基提供了快 速便捷的参数标定方法。 权利要求书1页 说明书5页 附图1页 CN 115169041 A 2022.10.11 CN 115169041 A 1.一种强夯碎石土PFC 颗粒滚动阻力摩擦系数 标定方法， 其特 征在于， 所述方法包括： 步骤S1， 获取强夯后碎石土地基土体重度γ、 孔隙率n、 粘聚力c及内摩擦角 并利用所 获取的所述土体重度γ、 孔隙率n、 粘聚力c及内摩擦角 根据卡尔曼法计算直立边坡的临 界高度hcr： 步骤S2， 利用筛分法获取地基土样的颗粒级配， 并基于获取的所述地基土样的颗粒级 配、 重度及孔隙率， 利用PFC2D颗粒流软件创建地基模 型， 其中， 通过调整滚动阻力摩 擦系数 的大小， 获得不同的地基模型； 步骤S3， 以所述直立边坡临界高度hcr为开挖高度， 依据所得的不同地基模型进行直立 开挖， 开挖坑壁保持稳定而不坍塌时所使用的地基模型对应的滚动阻力摩擦系数即为所确 定的与实际碎石土地基相适应的细观参数。 2.根据权利要求1所述的一种强夯碎石土PFC颗粒滚动阻力摩擦系数标定方法， 其特征 在于， 创建所述地基模型过程中所使用的颗粒间接触模型选取抗滚动线性接触模型， 颗粒 与墙体间接触模型选取线性接触模型， 所述抗滚动线性接触模型 的细观参数包括： 接触法 向刚度、 接触切向刚度、 摩 擦系数、 滚动阻力 摩擦系数； 所述线性接触模 型的细观参数包括： 接触法向刚度、 接触切向刚度、 摩擦系数， 其值与所述抗滚动线性接触模型的细观参数一 致。 3.根据权利要求2所述的一种强夯碎石土PFC颗粒滚动阻力摩擦系数标定方法， 其特征 在于， 所述摩擦系数 取值范围为1.0～1.2。 4.根据权利要求2所述的一种强夯碎石土PFC颗粒滚动阻力摩擦系数标定方法， 其特征 在于， 所述 步骤S1还 包括： 获取强夯后碎石土的弹性模量E， 所述法向刚度kn和切向刚度ks通过如下公式计算获得： 其中， E*为有效模量， 其根据所述弹性模量E获得， K*为刚度比， 其取值范围为1.0～1.5， A＝2r， r为颗粒半径， 对于 两个接触颗粒 取其中较小值， L 为两个接触颗粒圆心之间的距离 。 5.根据权利要求1所述的一种强夯碎石土PFC颗粒滚动阻力摩擦系数标定方法， 其特征 在于， 所述 通过调整滚动阻力摩擦系数的大小， 得到不同的地基模型包括： 令所述滚动阻力摩擦系数 μr初始值为0.1， 按照每次增加0.1依次计算， 建立不同的地基 模型。权　利　要　求　书 1/1 页 2 CN 115169041 A 2强夯碎石土PFC颗粒滚动阻力摩擦系数标定方 法 技术领域 [0001]本发明属于PFC数值模拟分析领域， 具体涉及一种强夯碎石土PFC颗粒滚动阻力摩 擦系数标定方法。 背景技术 [0002]PFC(颗粒流程序Particle  Follow Code)作为离散元数值模拟分析软件， 能很好 地从细观上反映出碎石土的不连续性， 对动力作用下的颗粒流动和大变形问题可以进 行很 好的模拟。 [0003]由于使用PFC软件建模时， 不能直接指定材料的宏观物理性质， 必须在物理试验结 果的基础上对参数进行标定， 以获取相应的细观参数， 使所建模型得到与实际材料宏观性 质一致的力学响应。 对于碎石土， 常规标定方法一般采用现场试验得到其粘聚力和内摩擦 角， 再通过PFC建立双轴或三轴数值试验模型， 进行不同围压条件下的压缩试验， 得出偏应 力—轴向应变曲线。 根据莫尔—库伦强度理论， 将计算出 的粘聚力和内摩擦角与实际宏观 参数相比较， 从而得 出粘聚力和内摩擦角与现场试验一 致时的细观参数。 [0004]上述方法具有一定 的盲目性， 且步骤繁琐， 可操作性不强。 因此， 需要提供一种更 加便捷可靠的PFC 颗粒滚动阻力摩擦系数 标定方法。 发明内容 [0005]鉴于现有技术中的上述缺陷或不足， 本发明旨在提供一种强夯碎石土PFC颗粒滚 动阻力摩擦系数 标定方法， 使碎石土的PFC细观参数 标定更加便捷可靠 。 [0006]为了实现上述目的， 本发明实施例采用如下技 术方案： [0007]本发明实施例提供一种强夯碎石土PFC颗粒滚动阻力摩擦系数标定方法， 所述方 法包括： [0008]步骤S1， 获取强夯后碎石 土地基土体重度γ、 孔隙率n、 粘聚力c及内摩擦角 并利 用所获取的所述土体重度γ、 孔隙率n、 粘聚力c及内摩擦角 根据卡尔曼法计算直立边坡 的临界高度hcr： [0009] [0010]步骤S2， 利用筛分法获取地基土样的颗粒级配， 并基于获取的所述地基土样的颗 粒级配、 重度及孔隙率， 利用PFC2D颗粒流软件创建地基模 型， 其中， 通过调整滚动阻力 摩擦 系数的大小， 获得不同的地基模型； [0011]步骤S3， 以所述直立边坡临界高度hcr为开挖高度， 依据所得的不同地基模型进行 直立开挖， 开挖坑壁保持稳定而不坍塌时所使用的地基模型对应的滚动阻力摩擦系数即为 所确定的与实际碎石土地基相适应的细观参数。 [0012]可选的， 创建所述地基模型过程中所使用的颗粒间接触模型选取抗滚动线性接触 模型， 颗粒与墙体间接触模型选取线性接触模型， 所述抗滚动线性接触模型 的细观参数包说　明　书 1/5 页 3 CN 115169041 A 3