基于PFC6.0泥石流淤积与运移模拟仿真应用（限时赠项目包）

详细信息

作者优秀: 优秀教师/意见领袖/博士学历/特邀专家/独家讲师
平台推荐: 内容稀缺

3月前浏览3153

导读：1月14日20时(明晚)，仿真秀2024新年设计仿真新年报告第四期邀请仿真秀年度优秀讲师Lobby做《离散元PFC仿真入门学习思路》公开课，分享PFC顶流创作者的PFC入门学习的见解（包括PFC现有的知识体系、PFC学习之路、和推荐2024年新版基于6.0创建的课程)。

一、Lobby的PFC6.0入门学习讲义

继PFC顶流：110万+单个视频播放量！愿作离散元池中的“鲇鱼”创作以来，引起了离散元仿真学习者的高度关注。目前在仿真秀单课播放量突破了259万次（全站排名第一）。为满足用户对PFC6.0学习需要，Lobby再次在仿真秀原创首发了《Lobby的PFC6.0入门学习讲义【更新中】》，当前超过200+名用户订阅和学习。

该课程全程配有字幕，方便大家能听懂Lobby的发音，另外查报错的过程也会被简化。下面为课程的目录，欢迎各位支持推荐。

本套课程注重实操，案例丰富，配有字幕。希望同学们可以借助这门课程完成PFC入门。课程更新中，从1月1号开始，剩余章节预计一周更新两章，一个月以内完成所有章节。目前上线课程到第17章节，前20章节的课程已录制结束，代码已在附件中，等待字幕配置。

课程可随时回放，可开具发票

讲师提供vip群知识圈答疑和模型下载

扫码立即试看

二、基于PFC6.0泥石流淤积与运移模拟

接下来，Lobby基于PFC6.0泥石流淤积与运移模拟案例应用。众所周知，泥石流在离散元方面的文章是比较多的，有二维做一些定性研究的，也有三维对真实的工况做模拟。本文以某深切沟的地形为案例，做一个简单的泥石流的运移模拟，在模拟过程中也可以观察到在地形上的淤积。

1、地形生成

本文和之前的一篇文章一样，用Microstation软件生成地形网格【【PFC6.0】三维真实边坡落石轨迹追踪】。Bentley的软件专门为公路设计，对于地形这种长范围的生成技术是比较友好的。一开始直接用等高线生成的网格文件过大，并且质量不高，导致wall上面有很多的空洞。而Microstation软件网格工具中有一个收缩包络网格的方法，可以生成比较均匀的地形网格，并且质量比较高。

等高线直接生成的网格

优化后的网格

生成地形的代码如下，geotools文件在之前的项目里面找就可以了。

2、泥石流生成

在三维空间中指定泥石流的生成范围，这里将后缘泥石流的推进，认为是在地形上方直接生成颗粒。为了模拟后面泥石流的推进过程，在计算过程中每隔15s会补充泥石流可以，所以在计算过程中可以发现颗粒数目是在不断增多的。文件是每5s保存一个sav，所以每保存3个sav就可以进行一次颗粒补充了。

count/3-count/3.0=0 这个是我们课程中提到的利用3和3.0的区别来进行是否可以整除的判断。

整体流程不复杂，代码量也不大。主要是参数的取值要慎重，对于这种动态问题，阻尼系数的影响是很大的。

model restore "dixing"[center=vector(45,63,20)][cubic_size=8]ball generate radius 0.2 0.4 number 4000 box [comp.x(center)-cubic_size] [comp.x(center)+cubic_size] ...                                [comp.y(center)-cubic_size] [comp.y(center)+cubic_size] ...                               [comp.z(center)] [comp.z(center)+cubic_size*2]                               contact cmat default model linear method deform emod 10e7 kratio 1.5 property fric 0.1ball attribute density 3e3 damp 0.2model gravity 9.8[baocunpinlv=5][time_record=mech.time.total-100][count=0]def savefile        if mech.time.total-time_record >= baocunpinlv then        filename=string.build("jieguo%1",count)        command            model save @filename        endcommand        if count/3-count/3.0=0 then             command            ball generate radius 0.2 0.4 number 500 box [comp.x(center)-cubic_size] [comp.x(center)+cubic_size] ...                                [comp.y(center)-cubic_size] [comp.y(center)+cubic_size] ...                               [comp.z(center)] [comp.z(center)+cubic_size*2]            ball attribute density 3e3 damp 0.2        endcommand        endif        time_record=mech.time.total        count +=1    endif    endfish callback add @savefile -1.0model solve time 500