试样为 φ39.1 mm × 80 mm 饱和冻融黏土圆柱体，我需要基于分数阶西原模型开发 Abaqus/Fortran UMAT，模拟循环交通载荷下土体完整三阶段流变响应。 重点研究变量为冻融循环次数、围压、动应力幅值和频率。 程序最终要能输出与动三轴试验一致的结果，包括轴向应力—应变曲线、滞回曲线和累积塑性变形曲线，并支持不同变量条件下的参数对比分析以及对应输出位移、应力、应变变化云图