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这里有几个骚操作：
- 用wall速度直接控制加载方向（vy负值下压）
- 通过应力阈值触发卸载（max_stress变号实现方向切换）
- 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。记得调参时备好枸杞茶，这玩意比写论文还费头发。
