flac动力分析中如何改变材料参数？

dell

希望找人讨论一下，flac2d中进行地震响应计算时，希望将单元剪切模量随剪应变增大而损伤降低的关系曲线通过fish语言实现，函数大概如下：
def  jiangcan
whilestepping

loop i (1,izones)
   loop j (1,jzones)
           shear_mod(i,j) = 8e9 *(1-ln(ssi(i,j)) )
endloop
endloop
end

问题是这个命令在静力计算中可以实现，即剪应变增大，剪切模量降低
但在动力计算中，可以计算，系统没有报错，但是剪切模量始终没有变化，不知何解，有人能帮助我么？