颗粒流环形剪切实验：用代码扒开土体的秘密

PFC3D5.0颗粒流『颗粒材料/土体材料环形剪切实验』完整代码 该代码包括 1完整代码及适量注释可以参考学习也可直接使用无需调试 2环形剪切实验的建模全过程剪切应力法向应力样品体积变化能量等完整fish函数。 附图模型侧视图俯视图加载过程中图像表观摩擦系数曲线剪切应力除以法向应力法向应力曲线恒压加载样品体积变化曲线。玩过堆沙子没颗粒材料在剪切力作用下的反应就跟沙堆塌了似的。环形剪切实验就是专门研究这个的今天咱们用PFC3D5.0整一套能直接拿来用的代码手把手教你搓出颗粒材料的剪切特性曲线。先看看怎么搭模型建模核心就三板斧造墙、撒颗粒、加约束。这里用圆柱形剪切盒下盘固定上盘转。初始化时记得把颗粒的摩擦系数设小点方便后续压实。;造墙 wall generate id 1 cylinder ... ;底下那个固定墙 wall generate id 2 cylinder ... ;上面要转的墙 wall attr damp 0.7 ;墙的阻尼别太小 ;撒芝麻 ball distribute ... ball attr density 2650 ;石英砂密度颗粒生成后要压实到目标孔隙比这步用伺服控制最稳。注意伺服过程中要实时监测接触力别压爆了。压实完成后把摩擦系数调回正常值毕竟实验时要模拟真实摩擦状态。这些参数你肯定得调剪切速率千万别设太高颗粒流最怕数值不稳定。建议用准静态条件剪切速率控制在1e-5到1e-4这个量级。法向应力加载用应力控制模式比位移控制更贴近实际工况。;加载控制 def apply_confining_stress loop local i (1,wall.num) if wall.id(i) 2 wall.force.app -target_stress * area ;压力换算成力 endif endloop end灵魂Fish函数大揭秘剪切过程中要实时记录五个关键指标剪切应力、法向位移、体积变化、动能变化、摩擦系数。这里有个取巧的方法——用墙的旋转角度换算剪切应变。;数据记录 def calc_shear_stress shear_force 0.0 loop foreach cnt wall.contact.list if wall.isin(cnt.end1) ;只统计和活动墙的接触 shear_force cnt.force.shear endif endloop current_shear shear_force / area ;换算成应力 end体积变化计算要注意排除墙位移的影响建议用颗粒的 centroid 坐标统计有效区域内的体积。动能监控能帮我们发现何时系统达到稳定状态这对确定数据采集时机至关重要。结果怎么看才专业当表观摩擦系数曲线剪切应力/法向应力出现稳定平台时说明材料进入临界状态。这时体积变化曲线应该趋于平稳要特别注意剪胀现象——密砂会膨胀松砂可能收缩。PFC3D5.0颗粒流『颗粒材料/土体材料环形剪切实验』完整代码 该代码包括 1完整代码及适量注释可以参考学习也可直接使用无需调试 2环形剪切实验的建模全过程剪切应力法向应力样品体积变化能量等完整fish函数。 附图模型侧视图俯视图加载过程中图像表观摩擦系数曲线剪切应力除以法向应力法向应力曲线恒压加载样品体积变化曲线。!法向应力保持恒定是检验伺服控制是否达标的关键。如果出现应力波动超过5%建议检查墙的刚度和阻尼参数设置。能量曲线里动能占比应该始终小于1%否则要怀疑是不是加载太快了。拿去就能用的骚操作代码包里已经预置了三组典型工况密砂/松砂/过渡状态。改这三个参数就能切换初始孔隙比0.5-0.8颗粒级配系数0.1-0.3围压大小50kPa-1MPa跑完自动生成四个数据文件应力应变曲线.dat体积变化.dat能量记录.dat实时参数.log遇到颗粒飞溅别慌八成是阻尼设小了。建议剪切阶段把局部阻尼调到0.3左右既能保证收敛速度又不影响力学响应。记住颗粒流仿真就像熬汤——火候时步到了自然出结果。