知识点二：初始地应力场-更改强度参数的弹塑性求解法

model new
zone create brick size 6 8 8
model title 'Tutorial Examples'
zone cmodel assign mohr-coulomb
zone property bulk 1e8 shear 0.3e8 fric 35
zone property cohesion 1e10 tension 1e10
model gravity 0 0 -9.81
zone initialize density 1000
zone gridpoint fix velocity-x range position-x 0.0
zone gridpoint fix velocity-x range position-x 6.0
zone gridpoint fix velocity-y range position-y 0.0
zone gridpoint fix velocity-y range position-y 8.0
zone gridpoint fix velocity-z range position-z 0.0

history interval 5
model history mechanical ratio
zone history displacement-z position 4 4 8

model solve
model save 'trench'
zone property cohesion 1e3 tension 1e3
;初始地应力场的生成
zone cmodel assign null range position-x 2 4 position-y 2 6 position-z 5 10
model largestrain on
zone gridpoint initialize displacement-x 0
zone gridpoint initialize displacement-y 0
zone gridpoint initialize displacement-z 0
model cycle 2000
model save 'trench2'

更改强度参数的弹塑性求解法生成是指求解过程中始终采用塑性模型，但为防止在计算过程中出现屈服区域，将粘聚力和抗拉强度设为大值，计算至平衡后，再将粘聚力和抗拉强度改为分析所采用的值计算至最终平衡状态。