求助：FLAC2D的模拟对摩尔库伦模型抗拉强度赋值后，没有任何命令更改材料特性，可是在初始平衡以开始求解时抗拉强度自动变化，但是体积和剪切模量都没有变化。命令如下，
请高手指正：
config
grid 496 177
m m
gen 0 0,0 1.768,4.96 1.768, 4.96 0
table 1 0.137 0.257,0 0.257,0 0,0.48 0, 0.48 0.05
gen table 1
gen arc 5.259,-2.772 1.413,1.768 19.15
table 2 4.48 0.05,4.48 0,4.96 0,4.96 0.257,4.823 0.257
gen table 2
gen arc -0.299,-2.772 4.823,0.257 19.15
model null region 38 131
model null region 446 117
table 3 0.902 0.666 1.218 0.666 1.385 0.845 2.48 1.025 3.575 &
0.845 3.742 0.666 4.058 0.666
gen table 3
table 4 1.521 1.33 1.629 1.222 1.749 1.149 1.89 1.138 2.33 1.206 &
2.63 1.206 3.07 1.138 3.211 1.149 3.331 1.222 3.439 1.33
gen table 4
gen arc 5.259,-2.772 1.672,1.462 19.18
gen arc -0.299,-2.772 4.48,0.05 19.18
unmark i 91 j 69
unmark i 407 j 69
unmark i 92 j 69
unmark i 93 j 69
unmark i 406 j 69
unmark i 405 j 69
model null region 248 2
mark i 407 j 69
mark i 406 j 69
mark i 405 j 69
mark i 93 j 69
mark i 92 j 69
mark i 91 j 69
gen line 1.413,1.768 1.672,1.462
gen line 3.288,1.462 3.547 1.768
model null region 248 176
；形成网格
group 'C:251' region 72 90
group 'C:252' region 438 78
group 'C:30' region 250 108
prop density=2400.0 bulk=1.42643E10 shear=1.22108E10 cohesion=550000.0 &
friction=30.0 tension=1800000.0 group 'C:251'
prop density=2400.0 bulk=1.42643E10 shear=1.22108E10 cohesion=550000.0 &
friction=30.0 dilation=0.0 tension=1800000.0 group 'C:252'
prop density=2400.0 bulk=1.5015E10 shear=1.28535E10 cohesion=550000.0 &
friction=30.0 tension=5.4E7 group 'C:30'
；区域分为三组，各自赋予材料特性参数
struct node 1 1.1702,1.3133
struct node 2 1.3329,1.1507
struct node 3 1.4733,1.0127
struct node 4 0.8866,1.0083
struct node 5 1.0599,0.8571
struct node 6 1.2077,0.7282
struct node 7 3.7898,1.3133
struct node 8 3.6271,1.1507
struct node 9 3.4867,1.0127
struct node 10 4.0734,1.0083
struct node 11 3.9,0.8571
struct node 12 3.7523,0.7282
struct rockbolt begin node 1 end node 2 seg 5 prop 4001
struct rockbolt begin node 2 end node 3 seg 5 prop 4002
struct rockbolt begin node 4 end node 5 seg 5 prop 4001
struct rockbolt begin node 5 end node 6 seg 5 prop 4002
struct rockbolt begin node 7 end node 8 seg 5 prop 4001
struct rockbolt begin node 8 end node 9 seg 5 prop 4002
struct rockbolt begin node 10 end node 11 seg 5 prop 4001
struct rockbolt begin node 11 end node 12 seg 5 prop 4002
struct prop 4001 e=210e9 area=3.14e-4 per=0.06 yield=4.7e4 ycomp=1.88e4 &
spacing=0.25 cs_scoh=4.0e5 cs_sstiff=3.5e7
struct prop 4002 e=210e9 area=3.14e-4 per=0.06 yield=2.3e4 ycomp=1.88e4 &
spacing=0.25 cs_scoh=2.0e5 cs_sstiff=1.75e7
；以锚杆单元模拟植筋联结
set gravity 9.81
fix y j=1
fix x i=1
fix x i=497
；设置重力和边界条件
step 20000
