知识点十：生成初始应力场

1.弹性求解：适用于浅埋的地表工程

model new
zone create brick size 4 6 8
zone face skin break-angle 45 
zone cmodel assign elastic
zone property bulk 3e7 shear 1e7
zone face apply velocity-normal 0 range group "East" or "West"
zone face apply velocity-normal 0 range group "Worth" or "South"
zone face apply velocity-normal 0 range group "Bottom"
zone initialize density 2000
model gravity 0 0 -10
model solve
model save 'initial'
;;;;;;;;;;达到初始平衡后，重置位移为零，只保留应力;;;;;;;;;
zone gridpoint initialize displacement-x 0
zone gridpoint initialize displacement-y 0
zone gridpoint initialize displacement-z 0
zone gridpoint initialize velocity-x 0
zone gridpoint initialize velocity-y 0
zone gridpoint initialize velocity-z 0

 

1

2.弹塑性求解：第一次把coh和ten最大，第二次把coh和ten改成正常值；

model new
zone create brick size 1 1 2
zone face skin break-angle 45 
zone cmodel assign mohr-coulomb
zone property bulk 3e7 shear 1e7 cohesion 10e10 friction 15 tension 1e10
zone face apply velocity-normal 0 range group "East" or "West"
zone face apply velocity-normal 0 range group "Worth" or "South"
zone face apply velocity-normal 0 range group "Bottom"
zone initialize density 2000
model gravity 0 0 -10
model solve
zone property bulk 3e7 shear 1e7 cohesion 10e3 friction 15 tension 0
model solve
model save 'initial'
;;;;;;;;;;达到初始平衡后，重置位移为零，只保留应力;;;;;;;;;
zone gridpoint initialize displacement-x 0
zone gridpoint initialize displacement-y 0
zone gridpoint initialize displacement-z 0
zone gridpoint initialize velocity-x 0
zone gridpoint initialize velocity-y 0
zone gridpoint initialize velocity-z 0

1

1

3.ini+弹塑性：ini赋值，可以加速平衡，效果一样；需要人为设置K0；

model new
zone create brick size 1 1 2
zone face skin break-angle 45 
zone cmodel assign mohr-coulomb
zone property bulk 3e7 shear 1e7 cohesion 10e3 friction 15 tension 0
zone property bulk 3e7 shear 1e7 fri 15 tension 1e10
zone face apply velocity-normal 0 range group "East" or "West"
zone face apply velocity-normal 0 range group "Worth" or "South"
zone face apply velocity-normal 0 range group "Bottom"
zone initialize density 2000
zone initialize stress-zz -40e3 gradient 0 0 20e3 range position-z 0 2
zone initialize stress-yy -20e3 gradient 0 0 10e3 range position-z 0 2
zone initialize stress-xx -20e3 gradient 0 0 10e3 range position-z 0 2
model gravity 0 0 -10
model solve
model save 'initial'
;;;;;;;;;;达到初始平衡后，重置位移为零，只保留应力;;;;;;;;;
zone gridpoint initialize displacement-x 0
zone gridpoint initialize displacement-y 0
zone gridpoint initialize displacement-z 0
zone gridpoint initialize velocity-x 0
zone gridpoint initialize velocity-y 0
zone gridpoint initialize velocity-z 0

我感觉地应力平衡的效果不太好，按照ABAQUS的地应力平衡标准，位移结果因为得很小~我也不太懂啥原因！！！

---

---

4.initialize-stresses：人为设置侧压力系数，默认侧压力系数为1，后面加关键词ratio设置；

model new
zone create brick size 15 15 15 edge=11
zone face skin
zone cmodel assign elastic
zone property bulk=5e9 shear=1e9 density=2000
zone face apply velocity-normal 0 range group 'East' or 'West'
zone face apply velocity-normal 0 range group 'North' or 'South'
zone face apply velocity-normal 0 range group 'Bottom'
model gravity 10
zone initialize-stresses ratio 0.3
model solve ratio-local 1e-4
model save 'initial'
;;;;;;;;;;达到初始平衡后，重置位移为零，只保留应力;;;;;;;;;
zone gridpoint initialize displacement-x 0
zone gridpoint initialize displacement-y 0
zone gridpoint initialize displacement-z 0
zone gridpoint initialize velocity-x 0
zone gridpoint initialize velocity-y 0
zone gridpoint initialize velocity-z 0

1

1

5.起伏地形（initialize-stresses）：zone initialize-stresses ratio 0.3 0.5

---

---

---