求助大家给我分享一下方自虎老师的子程序

我用abaqus软件做钢筋混凝土剪力墙抗震模拟，滞回曲线效果很差，没有捏缩现象，前期刚度过大，滞回环呈锯齿状分布，多次调试混凝土CDP模型参数也没有效果，想用方自虎老师的子程序来改善一下滞回曲线的捏缩，但是我找不到方自虎老师的子程序，求助大家给我分享一下方自虎老师的子程序。