PFC循环加卸载实战：单轴双轴不同姿势解锁

pfc5.0/6.0单轴双轴不同应力路径循环加卸载程序

岩土佬们搞循环加载实验，最头疼的就是控制不同应力路径。PFC5.0/6.0这货虽然强大，但面对复杂应力路径时还是得讲究操作手法。今天咱们直接上干货，手把手教你怎么玩转单轴双轴的循环加载。

单轴加载：从入门到放弃

先看个基础的单轴压缩循环代码：

wall delete wall generate box ball distribute ... ...（试样生成代码略） fish define cyclic_load loop_cnt = 0 max_stress = 1e6 while loop_cnt < 10 wall.top.vy = -0.01 ; 加载阶段 command solve ratio 1e-5 end_command stress = wall.top.force.y/wall.top.area if stress > max_stress wall.top.vy = 0.01 ; 卸载反转 max_stress = -max_stress loop_cnt += 1 endif endwhile end @cyclic_load

这里有几个骚操作：

1. 用wall速度直接控制加载方向（vy负值下压）
2. 通过应力阈值触发卸载（max_stress变号实现方向切换）
3. solve ratio控制求解精度，建议不要低于1e-5（别问我是怎么知道的）

特别注意vstrainrate参数设置，新手最容易在这里翻车。建议初始值设在0.001~0.01/s之间，太大会导致颗粒"穿模"。

双轴加载：左右横跳的艺术

双轴才是真考验，得同时伺候x、y两个方向。上硬菜：

fish define biaxial_servo target_stress = 2e6 servogain = 1e-3 loop_cnt = 0 while loop_cnt < 5 ; X方向伺服 sx = (wall.left.force.x - wall.right.force.x)/area_x wall.left.vx = servogain*(target_stress - sx) ; Y方向循环加载 if loop_cnt % 2 == 0 wall.top.vy = -0.005 else wall.top.vy = 0.005 endif command solve ratio 1e-5 cycles 5000 end_command loop_cnt += 1 endwhile end

这段代码的亮点在于：

• 用servo增益系数实现应力控制（类似现实中的三轴仪）
• X方向用伺服控制保持恒定围压
• Y方向交替加载形成循环路径
• solve cycles限定计算步数避免无限循环

遇到过试样扭曲的问题？试试在初始化时加上：

ball attribute displacement multiply 0 ; 清空残余位移 contact method deform ; 启用变形接触模型

路径切换的骚操作

想要实现复杂应力路径（比如从q-p路径切换到εv-εq路径），记得在循环内动态修改target参数：

if loop_cnt > 5 target_q = 1.5*sin(loop_cnt*0.5) ; 正弦波加载 target_p = 2e6 + loop_cnt*1e5 call adjust_stress(target_q, target_p) ; 自定义应力调整函数 endif

这里推荐用傅里叶级数生成复杂波形，比单纯三角波有意思多了。不过要当心数值震荡，建议配合自适应步长算法：

fish define auto_step current_ratio = mech.ratio if current_ratio > 1e-4 mech.set timestep scale 0.8 else mech.set timestep scale 1.2 endif end

最后给个忠告：跑循环前务必保存初始状态！见过太多人一晚上白跑的惨剧。用这个保命命令：

save cycle_0.sav

搞定这些，你的循环加载应该能跑得比隔壁组的师弟还6。记得调参时备好枸杞茶，这玩意比写论文还费头发。