And by means of Oula difference method, the unified constitutive equations numerical solutions were also obtained.

• 采用欧拉差分法，得到了统一本构方程的数值解。