flac3d耦合pfc3d边坡模型位移连续性良好。最近在折腾边坡稳定性分析的时候发现纯连续介质模型遇到碎裂带就拉胯PFC颗粒流倒是能模拟破碎过程但整体位移场总像拼图少了几块。硬是把FLAC3D和PFC3D这两冤家撮合在一起整出个耦合模型效果居然比相亲节目牵手成功还自然。先甩个耦合界面的配置代码镇楼zone face apply velocity-x (funcget_vel_x()) range group interface zone face apply velocity-y (funcget_vel_y()) range group interface # PFC侧颗粒绑定 fish define bind_particles loop foreach local ball ball.list if ball.pos.z interface_z - 0.1 ball.fix 1 ball.group boundary endif endloop end这段代码暗藏乾坤——FLAC用函数指针动态获取颗粒边界速度PFC这边用fish函数锁死界面层颗粒。特别是那个0.1的容差范围调参时差点让我头秃。太大会把活动颗粒误杀成僵尸节点太小又容易漏绑颗粒搞出鬼影界面。耦合最骚的操作在数据交换。来看个应力传递的硬核操作// 自定义耦合模块 void CouplingModule::transferStress() { for (auto zone : flacZones) { Tensor stress zone-getStress(); vectorParticle* neighbors pfcGrid-queryNeighbors(zone-centroid()); for (auto p : neighbors) { double weight calcWeight(zone-centroid(), p-pos); p-addStress(stress * weight); // 应力加权分配 } } }这里用空间网格加速查找周边颗粒比暴力遍历快20倍不止。但权值函数要是用不好应力分配就像撒胡椒粉——要么颗粒扎堆吃撑要么边缘颗粒饿死。实测高斯衰减比线性插值靠谱半径取1.5倍平均粒径时数据跳变最小。flac3d耦合pfc3d边坡模型位移连续性良好。位移连续性验证得整点实在的。跑完模型截取剖面数据% 位移场比对脚本 flac_disp load(flac_disp.dat); pfc_disp load(pfc_disp.dat); smooth_disp movmedian(pfc_disp, 50); % 中值滤波去颗粒抖动 figure; subplot(211) contourf(flac_disp, LineColor,none) title(连续区位移场) subplot(212) scatter(pfc_disp(:,1), pfc_disp(:,2), 10, smooth_disp, filled) colorbar title(离散区位移场)关键在颗粒位移的中值滤波窗口大小——窗口太小滤不干净颗粒碰撞噪声太大又会抹平真实位移梯度。当滑动面形成时连续区的位移突变线会和颗粒流滑带完美衔接就像高速公路突然变成石子路但路面高度愣是没断差。耦合模型最怕出现量子纠缠现象明明参数设的一样跑出来的结果每次都不重样。后来发现是FLAC的显式求解和PFC的时步没同步就像两个人跳舞各踩各的拍子。加个时步耦合控制器才搞定-- 时步同步脚本 while sim.time total_time do local flac_dt flac.solve.step() local pfc_dt pfc.cycle(flac_dt) actual_dt math.min(flac_dt, pfc_dt) flac.set.timestep(actual_dt) pfc.set.timestep(actual_dt) end这招让两者步调一致避免出现FLAC已经算到2030年PFC还在2023年磨叽的魔幻场景。不过收敛判据得放宽点否则容易陷入无限小步长死循环。折腾下来最大的感悟耦合模型像个矫情的对象得顺着它的脾气来。连续区和离散区的参数不能完全独立设置比如PFC颗粒的刚度至少得是FLAC单元刚度的10倍否则会出现颗粒陷入连续体的灵异现象。但真要较真起来这参数匹配问题够写三篇SCI了。
FLAC3D与PFC3D耦合边坡模型，位移连续性优异
flac3d耦合pfc3d边坡模型位移连续性良好。最近在折腾边坡稳定性分析的时候发现纯连续介质模型遇到碎裂带就拉胯PFC颗粒流倒是能模拟破碎过程但整体位移场总像拼图少了几块。硬是把FLAC3D和PFC3D这两冤家撮合在一起整出个耦合模型效果居然比相亲节目牵手成功还自然。先甩个耦合界面的配置代码镇楼zone face apply velocity-x (funcget_vel_x()) range group interface zone face apply velocity-y (funcget_vel_y()) range group interface # PFC侧颗粒绑定 fish define bind_particles loop foreach local ball ball.list if ball.pos.z interface_z - 0.1 ball.fix 1 ball.group boundary endif endloop end这段代码暗藏乾坤——FLAC用函数指针动态获取颗粒边界速度PFC这边用fish函数锁死界面层颗粒。特别是那个0.1的容差范围调参时差点让我头秃。太大会把活动颗粒误杀成僵尸节点太小又容易漏绑颗粒搞出鬼影界面。耦合最骚的操作在数据交换。来看个应力传递的硬核操作// 自定义耦合模块 void CouplingModule::transferStress() { for (auto zone : flacZones) { Tensor stress zone-getStress(); vectorParticle* neighbors pfcGrid-queryNeighbors(zone-centroid()); for (auto p : neighbors) { double weight calcWeight(zone-centroid(), p-pos); p-addStress(stress * weight); // 应力加权分配 } } }这里用空间网格加速查找周边颗粒比暴力遍历快20倍不止。但权值函数要是用不好应力分配就像撒胡椒粉——要么颗粒扎堆吃撑要么边缘颗粒饿死。实测高斯衰减比线性插值靠谱半径取1.5倍平均粒径时数据跳变最小。flac3d耦合pfc3d边坡模型位移连续性良好。位移连续性验证得整点实在的。跑完模型截取剖面数据% 位移场比对脚本 flac_disp load(flac_disp.dat); pfc_disp load(pfc_disp.dat); smooth_disp movmedian(pfc_disp, 50); % 中值滤波去颗粒抖动 figure; subplot(211) contourf(flac_disp, LineColor,none) title(连续区位移场) subplot(212) scatter(pfc_disp(:,1), pfc_disp(:,2), 10, smooth_disp, filled) colorbar title(离散区位移场)关键在颗粒位移的中值滤波窗口大小——窗口太小滤不干净颗粒碰撞噪声太大又会抹平真实位移梯度。当滑动面形成时连续区的位移突变线会和颗粒流滑带完美衔接就像高速公路突然变成石子路但路面高度愣是没断差。耦合模型最怕出现量子纠缠现象明明参数设的一样跑出来的结果每次都不重样。后来发现是FLAC的显式求解和PFC的时步没同步就像两个人跳舞各踩各的拍子。加个时步耦合控制器才搞定-- 时步同步脚本 while sim.time total_time do local flac_dt flac.solve.step() local pfc_dt pfc.cycle(flac_dt) actual_dt math.min(flac_dt, pfc_dt) flac.set.timestep(actual_dt) pfc.set.timestep(actual_dt) end这招让两者步调一致避免出现FLAC已经算到2030年PFC还在2023年磨叽的魔幻场景。不过收敛判据得放宽点否则容易陷入无限小步长死循环。折腾下来最大的感悟耦合模型像个矫情的对象得顺着它的脾气来。连续区和离散区的参数不能完全独立设置比如PFC颗粒的刚度至少得是FLAC单元刚度的10倍否则会出现颗粒陷入连续体的灵异现象。但真要较真起来这参数匹配问题够写三篇SCI了。
