def kaiwa_7_40
loop ii(7,8)
filename_1='step_'+string(ii)+'.sav'
temp1=(ii-1)*1.5
temp2=ii*1.5
temp3=temp2+9
temp4=(ii-7)*1.5
temp5=(ii-6)*1.5
temp6=ii-6
temp7=100+temp6
temp8=200+temp6
command
mod null range y temp1 temp2 group tunnel
apply nstress 55e3 range y temp2 group tunnel
;shield machine
;sel delete shell id=temp6
sel shell id=ii range cyl end1 0 temp1 0 end2 0 temp2 0 rad 4.5
sel shell id=ii prop isotropic=(200e9,0.25) dens=7840 thick=0.3
;injection grout
sel shell id=temp7 range cyl end1 0 temp4 0 end2 0 temp5 0 rad 4.5
sel shell id=temp7 prop isotropic=(1e9,0.25) dens=1200 thick=0.15
;lining rings segment
sel shell id=temp8 range cyl end1 0 temp4 0 end2 0 temp5 0 rad 4.35
sel shell id=temp8 prop isotropic=(3.45e10,0.3) dens=2450 thick=0.35 (这里会出现问题,提示shell 某id号 does not exist)
step 1500
apply remove nstress range y temp2 group tunnel
save filename_1
endcommand
endloop
end
kaiwa_7_40