FLAC3D隧道蠕变加支护，需要的人看过来。如若有错，望指...
本帖最后由 706811291 于 2017-4-5 10:37 编辑model elastic
prop bulk 10e9 shear 10e9 dens 2000 range z -80.7836 79.2164
set gravity 0 0 -10
ini sxx -4.8e6 grad 0,0,0.036e6
ini syy -4.8e6 grad 0,0,0.036e6
apply szz -3.6e6 grad 0 0 1.8e4 range z 79.1164 79.3164
fix x range x -79.8631 -79.8831
fix x range x 80.1169 80.1369
fix x y z range z -80.6836 -80.8836
fix y range y -0.01 0.01
fix y range y 69.99 70.01
step 1000
solve
save c1
ini xdis 0 ydis 0 zdis 0
ini xvel 0 yvel 0 zvel 0
ini density 2200 range gr 4
ini density 2200 range gr 1 any gr 2 any gr 3 any
config creep
def modulus
s_mod=Emod/(2.0*(1.0+p_ratio))
b_mod=Emod/(3.0*(1.0-2.0*p_ratio))
end
set Emod=1e9
set p_ratio=0.4
modulus
model cvisc ran gr 4
prop bulk b_mod coh 2.0e5 fric 25 di 2.68 tension 1.0e5 ran gr 4
prop kshear 2.2e8 kviscosity 9.6e16 mshear s_mod mviscosity 2.0e17 ran gr 4
model cvisc range gr 1 any gr 2 any gr 3 any
prop bulk b_mod coh 1.9e5 fric 24 di 2.65 tension 9.9e4 range gr 1 any gr 2 any gr 3 any
prop kshear 2.1e8 kviscosity 9.5e16 mshear s_mod mviscosity 1.9e17 range gr 1 any gr 2 any gr 3 any
set creep mindt=1e-2 maxdt=1
set creep on
his unbal
his gp zdisp 0 0 6.4621
his gp xdisp -7.1910 0 1.6377
his gp xdisp 7.1922 0 1.6395
his crtime
def excav_ms
loop i(1,14)
step_rock1=5*i-5
step_rock2=5*i
command
mo null ran gr 1 y step_rock1 step_rock2
end_command
end_loop
end
excav_ms
step 500
save c2
def ins_cab_gd
loop iidx(1,70)
y=iidx*1.0
command
sel cable id iidx begin -11.0289 y 3.3533 end -7.1910 y 1.6377 nseg 4
sel cable id iidx begin -10.3380 y 4.7328 end -6.7490 y 2.5435 nseg 4
sel cable id iidx begin -9.4698 y 6.0180 end -6.1890 y 3.3868 nseg 4
sel cable id iidx begin -8.4427 y 7.1810 end -5.5242 y 4.1520 nseg 4
sel cable id iidx begin -7.2708 y 8.2041 end -4.7643 y 4.8324 nseg 4
sel cable id iidx begin -5.9739 y 9.0672 end -3.9161 y 5.3978 nseg 4
sel cable id iidx begin -4.5737 y 9.7559 end -2.9968 y 5.8502 nseg 4
sel cable id iidx begin -3.0852 y 10.2577 end -2.0219 y 6.1832 nseg 4
sel cable id iidx begin -1.5614 y 10.5619 end -1.0240 y 6.3873 nseg 4
sel cable id iidx begin -0.0011 y 10.7139 end 0.0000 y 6.462 nseg 4
sel cable id iidx begin 1.5603 y 10.5623 end 1.0240 y 6.3876 nseg 4
sel cable id iidx begin 3.0842 y 10.2585 end 2.0219 y 6.1837 nseg 4
sel cable id iidx begin 4.5728 y 9.7571 end 2.9969 y 5.8510 nseg 4
sel cable id iidx begin 5.9732 y 9.0688 end 3.9164 y 5.3988 nseg 4
sel cable id iidx begin 7.2704 y 8.2060 end 4.7647 y 4.8337 nseg 4
sel cable id iidx begin 8.4425 y 7.1832 end 5.5248 y 4.1534 nseg 4
sel cable id iidx begin 9.4699 y 6.0204 end 6.1898 y 3.3884 nseg 4
sel cable id iidx begin 10.3385 y 4.7355 end 6.7500 y 2.5453 nseg 4
sel cable id iidx begin 11.0297 y 3.3562 end 7.1922 y 1.6395 nseg 4
sel cable prop dens 7800 emod 210e9 xcarea 3.80e-4 gr_per 0.314 gr_fric 25 ytens 40e4 gr_k 17.5e6 gr_coh 20e4 range id iidx
end_command
end_loop
end
ins_cab_gd
step 500
solve
save c3
def excav_mx
loop j(1,14)
step_rock3=5*j-5
step_rock4=5*j
command
mo null ran gr 2 y step_rock3 step_rock4
end_command
end_loop
end
excav_mx
sel shell id=1 ran gr 4 x -7.3084 7.3097 z 0 6.4621 y 0.1 69.9
sel shell prop iso=(30e9,0.2) thick=0.3 dens 2500 range id=1
step 500
save c4 第一次见人用step做蠕变~ 我也是第一次，计算怎么样还得看结果吧 学习一下
页:
[1]