new
gen zone brick size 12 8 5
model mech mohr
prop bulk 1e8 shear 0.3e8 fric 35
prop coh 1e10 tens 1e10
set grav 0 0 -10
ini dens 2000
fix x y range x -0.1 0.1
fix x y range x 11.5 12.5
fix x y range y -0.1 0.1
fix x y range y 7.5 8.5
fix x y z range z -0.1 0.1
solve
;
pause
prop coh 1e3 tens 1e3
model mech null range x 3 9 y 3 5 z 3 5
set large
ini xdis 0 ydis 0 zdis 0
sel cable id 1 beg 3.5 3 5 end 3.5 1 4 nseg 4
sel cable id 1 beg 4.5 3 5 end 4.5 1 4 nseg 4
sel cable id 1 beg 5.5 3 5 end 5.5 1 4 nseg 4
sel cable id 1 beg 6.5 3 5 end 6.5 1 4 nseg 4
sel cable id 1 beg 7.5 3 5 end 7.5 1 4 nseg 4
sel cable id 1 beg 8.5 3 5 end 8.5 1 4 nseg 4
sel cable id 1 beg 3 3.5 5 end 2 3.5 4 nseg 4
sel cable id 1 beg 3 4.5 5 end 2 4.5 4 nseg 4
sel cable id 1 prop xcarea=5.02e-3 emod=200e9 ytens=1e10 &
gr_k=7e6 gr_coh=1e2 gr_per 0.0875 range selid 1 4
sel cable id 1 prop xcarea=5.02e-3 emod=200e9 ytens=1e10 &
gr_k=7e6 gr_coh=1e2 range selid 5 9
sel pile id 2 beg 3 3 5 end 3 3 2.5 nseg 4
sel pile id 2 beg 4 3 5 end 4 3 2.5 nseg 4
sel pile id 2 beg 5 3 5 end 5 3 2.5 nseg 4
sel pile id 2 beg 6 3 5 end 6 3 2.5 nseg 4
sel pile id 2 beg 7 3 5 end 7 3 2.5 nseg 4
sel pile id 2 beg 8 3 5 end 8 3 2.5 nseg 4
sel pile id 2 beg 9 3 5 end 9 3 2.5 nseg 4
sel pile id 2 beg 3 4 5 end 3 4 2.5 nseg 4
sel pile id 2 beg 3 5 5 end 3 5 2.5 nseg 4
sel pile prop Emod=8.0e10 Nu=0.30 XCArea=0.5024 &
XCJ=0.0 XCIy=0.0 XCIz=0.0 &
Per=3.14 &
CS_sK=1.3e11 CS_sCoh=0.0 CS_sFric=30.0 &
CS_nK=1.3e11 CS_nCoh=0.0 CS_nFric= 0.0 CS_nGap=off
sel beam id 3 begin=( 3, 3, 5) end=( 9,3,5) nseg=6
sel beam id 4 begin=( 3, 3, 5) end=( 3,5,5) nseg=2
sel beam prop emod=2.0e11 nu=0.30
sel beam prop XCArea=0.6 XCIz=0.018 XCIy=0.018 XCJ=0.0
SEL JOIN Cable ID 1 beam ID 3 RANG Z 4.9 5.1
SEL JOIN pile ID 2 beam ID 3 RANG Z 4.9 5.1
SEL JOIN beam ID 3 beam ID 4 rang Z 4.9 5.1
SEL JOIN pile ID 2 beam ID 4 RANG Z 4.9 5.1
hist unbal
hist id=5 gp ydis 3 3 5
hist id=6 gp ydis 4 3 5
hist id=7 gp ydis 5 3 5
hist id=8 gp ydis 6 3 5
hist id=9 gp ydis 7 3 5
hist id=10 gp ydis 8 3 5
hist id=11 gp ydis 9 3 5
hist id=12 gp dis 3 4 5
hist id=13 gp dis 3 5 5
solve
刚开始学习flac 3d。以上命令是模拟一个基坑的开挖,但是在求解过程中最大不平衡力不能平衡,请问各位专家怎么解决该问题?怎么输出支护桩和冠梁的弯矩、剪力图(不要云图)?