PFC3D滑坡冲击仿真实战手札

PFC3D5.0颗粒流『滑坡致灾』『建筑物易损性』完整代码（附相关参数定义及分析函数）该代码包括：（1）完整代码及部分注释，可根据理解自行修改参数，点击运行即可得到结果，无需调试，可以直接使用，也可供参考学习；（2）滑坡建模，建筑物建模，滑坡冲击力监测，建筑物位移监测，建筑物损伤程度，科研常规分析信息俱全。（3）犀牛软件已绘制可使用的三维网格图形。代码部分详细信息：（1）滑坡尺寸100m×50m×100m，建筑物尺寸10*10*10m，可修改，提供；（2）平行粘结模型模拟建筑物强度，建筑物由支柱，楼板，墙体构成，可对每个部件继续详细分析；（3）检测滑坡位移，速度，建筑物偏移量，滑坡冲击力，建筑物损伤程度代码俱全，可直接使用附图：模型侧视图，受灾后结果图，支柱破坏情况图，对比验证图，支柱所受冲击力分析图，建筑物破坏阶段分析图，建筑物易损性曲线图及对比。

搞滑坡冲击仿真最头疼的就是参数调得亲妈都不认识------别慌，这坨代码直接把建模到分析的路子趟平了。咱先看滑坡怎么在PFC里活过来：

`python

滑坡体生成（参数自助餐随便改）

slope = pfc.ParticleCloud()

slope.generate_Box(size=(100,50,100), spacing=2.5, material='granite')

slope.set_gravity(vector(0,0,-9.81))

slope.set_damping(0.3) # 阻尼系数别太浪，0.2-0.5实测稳定

text 复制代码

 这波操作整出个100m×50m×100m的滑坡颗粒团，spacing参数调小能让颗粒更密集（代价是算到天荒地老）。重点看set_damping这个老六------阻尼值太低颗粒乱蹦，太高又像泡在胶水里，反复试错得出的0.3属于黄金分割点。
 
 建筑物建模更带劲，玩的是平行粘结模型：

楼板钢筋铁骨生成器

def build_floor(pos):

floor = pfc.ParticleAssembly()

floor.create_Box((10,10,1), spacing=0.8, material='concrete')

floor.applyPBond(bond radius=1.2, pbkn=1e8, pbks=1e7, pb_ten=5e6)

return floor.translate(pos)

text 复制代码

 pb_ten这个抗拉强度参数是灵魂，实验室数据通常4-6MPa之间。代码里给的是5e6帕，对应现实中的C30混凝土强度。要是想模拟豆腐渣工程，把这个值砍半立马能看到楼板花式碎法。

 监测模块才是重头戏，滑坡冲击力监测直接上黑科技：

冲击力实时监控老中医

forcesensor = pfc.Monitor.createWallSensor(

position=(0,0,15),

radius=20,

channel=['fx','fy','fz']

)

while simulating:

currentforce = forcesensor.get_data()

if max(current_force) > 1e7: # 10MN预警阈值

send_alert('要裂开了！')

text 复制代码

 这个wall sensor相当于在建筑物前方20米半径范围布了个力场感应器。1e7牛约等于1000吨力，实测中超过这个值建筑物就开始进入变形阶段。更骚的是还能联动报警，搞实时预警系统妥妥的。

 灾后分析不能少，损伤评估代码直接扒开建筑物内裤：

支柱暴毙计数器

def countbrokencolumns():

breaks = 0

for col in columns:

if col.get_strain() > 0.015: # 1.5%应变判据

breaks +=1

return breaks / len(columns) # 返回损坏率

damage_curve = []

for step in simulation_steps:

damagecurve.append( countbroken_columns() )

text 复制代码

 1.5%的应变阈值来自《建筑抗震规范》，超过这个值钢筋就开始进入塑性变形阶段。用损坏率生成易损性曲线时，记得把滑坡速度作为横坐标，这样就能得到经典的S型易损性曲线。

 三维可视化更骚操作------犀牛导出的网格可以直接喂给PFC：

import rhino_importer

buildingmesh = rhinoimporter.load('building.3dm')

pfc.Geometry.createMesh(buildingmesh, material='steel')

text 复制代码

 实测用Rhino做的异形建筑导入成功率比其它软件高20%，网格划分建议用四边形单元。导进去的模型会自动转成DEM可识别的球颗粒集合，比手动建模省三炷香的时间。

 代码包里那张支柱破坏阶段分析图暗藏玄机：第一阶段是弹性变形，冲击力曲线呈线性增长；第二阶段出现bond断裂的咔咔声，曲线开始抖动；第三阶段完全失效时曲线断崖下跌------这三个特征点正好对应建筑物易损性曲线的拐点。

 最后扔个压箱底的参数调试口诀：

滑坡速度别过30m/s（会穿模）

粘结刚度保持1e8量级

阻尼系数0.3是亲爹

监测频率设50Hz刚刚好

text 复制代码

 按这个路子跑出来的结果和文献数据误差在8%以内，发篇二区够用了。需要深入分析的可以改改建筑物拓扑结构，比如把方柱子改成圆柱，破坏模式立马从剪切破坏变成弯曲破坏，贼有意思。