flac3d5.0软件隧道支护与开挖命令流.txt
《flac3d5.0软件隧道支护与开挖命令流.txt》由会员分享，可在线阅读，更多相关《flac3d5.0软件隧道支护与开挖命令流.txt（6页珍藏版）》请在冰豆网上搜索。
;-----------------------------------------------------
;----ExcavationandSupportforaShallowTunnel---
;-----------------------------------------------------
new;�½���Ŀ
setfishautocreateoff
title'ExcavationandSupportforaShallowTunnel';������Ŀ
;generateprimitivecomponentsofgrid
;concreteliner-uppertunnel
genzoncshellp0000p1700p20510p3005.5&
dim5555size25110
groupzone'concreteliner'
;
;uppertunnel
genzoncylinderp0000p1500p20510p3005&
size55110
groupzonetunnelrangegroup'concreteliner'not
;
;lowertunnel&liner
genzonebrickp000-4.5p1add700p2add0510p3add004.5&
size7513
;
;surroundingrock（8primitives）
genzonradcylp0000p12700p20510p30025&
dim75.575.5size551108rat1111.3
;
genzonebrickp070-4.5p1270-15p2add0510p3700&
p42751-15p57510p62700p727510&
size8513ratio1.311
;
genzonebrickp000-15p1add2700p2add0510p300-4.5&
p42751-15p5051-4.5p670-4.5p7751-4.5&
size7518rat110.7692307692307692
;
genzonbrickp00025p1add2700p2add0510p3add0010&
size5512
;
genzonbricp027025p1add1700p2add0510p3add0010&
size2512rat211
;
genzonbricp0270-15p1add1700p2add0510p3add0040&
size2518rat211
;
genzonbricp0270-40p1add1700p2add0510p3add0025&
size2512rat210.5
;
genzonbricp000-40p1add2700p2add0510p3add0025&
size7512rat110.5
;
;assignnamestogroupsofzones
groupzonerockrangegroup'concreteliner'notgrouptunnelnot
;
;assignMohr-Coulombmaterialmodel
modelmechmohr
probulk50e6she18e6fric20coh25e3ten0dil0rangez2535
probulk4e8she1.5e8fric20coh50e3ten5e3dil3rangez-5025
;assignboundaryconditions;ʩ�ӱ߽��������������ֱ����Ϊ���������
fixxrangex-.1.1
fixxrangex43.944.1
fixzrangez-40.1-39.9
fixyrangey-.1.1
fixyrangey50.951.1
;assigninitialstressstate;��ʼӦ��״̬
setgrav00-10
inidensity2200
iniszz-770e3grad0022000
inisxx-770e3grad0022000
inisyy-385e3grad0011000;ʩ�ӳ�ʼӦ��
;monitorvariablesinmodel;ģ���б����ļ��
histaddunbal;��ز�ƽ����
histaddgpzdisp005.5
histaddgpxdisp700
histaddgpzdisp000
histaddgpzdisp0035
histaddgpzdisp0305.5
histaddgpxdisp7300
histaddgpzdisp0300
histaddgpzdisp03035
histaddgpzdisp01235
histaddgpzdisp01835
histaddgpzdisp02435
histaddgpzdisp03635
histaddgpzdisp53035
histaddgpzdisp103035;���������Щ���x��z����λ�Ʊ仯
;
savgeom1
;
defconc_parm;����֧��������parm��parameter����������˼
globalbmc=20.7e9;�������ģ��Ϊȫ�ֱ�����b����bulk��m����modulus��c����concrete
globalsmc=12.6e9;�������ģ��Ϊȫ�ֱ�����s����shear��m����modulus��c����concrete
end
@conc_parm
;
;definethelocationsofcablepatterns1,2and3
;
defcab_parm;����ê�˲���
globalx_b=get_array（4,3）;�������飨4,3������ê�˵�λ��
globalz_b=get_array（4,3）
globaly0=-3;��ê�˵�λ�����������ʽ��ʾ������x,y,z��
x_b（1,1）=0.8
x_b（2,1）=2.1
x_b（3,1）=3.5
x_b（4,1）=5.5
z_b（1,1）=5.5
z_b（2,1）=2.4
z_b（3,1）=4.7
z_b（4,1）=1.5
x_b（1,2）=0.8
x_b（2,2）=0.8
x_b（3,2）=3.5
x_b（4,2）=5.5
z_b（1,2）=0.6
z_b（2,2）=4.0
z_b（3,2）=2.4
z_b（4,2）=0.6
x_b（1,3）=0.8
x_b（2,3）=2.6
x_b（3,3）=5.0
x_b（4,3）=3.5
z_b（1,3）=2.4
z_b（2,3）=4.0
z_b（3,3）=3.0
z_b（4,3）=0.6;��ÿ��ê�˵�λ�ý��и�ֵ
end
definip（iidx）;�����ʼê��λ��
globalx1=x_b（1,iidx）
globalx2=x_b（2,iidx）
globalx3=x_b（3,iidx）
globalx4=x_b（4,iidx）
globalz1=z_b（1,iidx）
globalz2=z_b（2,iidx）
globalz3=z_b（3,iidx）
globalz4=z_b（4,iidx）
end
@cab_parm
;
;installinitialcables;��װ��ʼê��
;
defins_cab;����Ҫ��ʼ��װ��ê�ˣ�����ins_cab��ʾinstallinitialcables
globaliidx;����ȫ�ֱ���iidx
globalcab_seg;����ê�˻��ֵ�Ԫ��
globalcab_seg_m;ê�˳��ȣ��д���һ��ȷ�ϣ�
loopiidx（1,3）
inip（iidx）
cab_seg=cab_seg_m-3*（3-iidx）
globaly1=0.
globaly2=float（cab_seg）
command
selcableid@iidxbegin@x1@y1@z1end@x1@y2@z1nseg@cab_seg
selcableid@iidxbegin@x2@y1@z2end@x2@y2@z2nseg@cab_seg
selcableid@iidxbegin@x3@y1@z3end@x3@y2@z3nseg@cab_seg
selcableid@iidxbegin@x4@y1@z4end@x4@y2@z4nseg@cab_seg
selcableproemod45e9xcarea1.57e-3gr_per1.0&
yten25e4gr_k17.5e6gr_c20e4rangeid@iidx;ʩ����ʼê��
end_command
end_loop
end
set@cab_seg_m15
@ins_cab
;installpre-supportconcrete;Ԥ֧��
;
selshellid10grouprockrangecylend100-1.5end201-1.5rad7.4&
cylend100-1.5end201-1.5rad6.7not&
z-0.16
selshellpropisotropic10.5e9,0.25thickness0.3density2500
defmonit
globalipt_surf=gp_near（0,30,35）;�ر�
globalipt_crown=gp_near（0,30,5.5）;����
globalipt_spring=gp_near（7,30,0）
end
@monit
savem_ini
table1name'groundsurfaceattunnelcenterline'
table2name'tunnelcrown';�������
table3name'tunnelsidewall';����������֣������ǽ
;
;FISHfunctiontocontrolexcavationandsupportsequence
defexcav
y0=y0+3
localcut_i=y0/3+1
globalcut
loopcut（cut_i,16）
localcut_cur=cut
localii=out（'EXCAVATIONSTEP'+string（cut））
y0=3*（cut-1）
y1=y0+3
globalyp0=y0+1
globalyp1=y1+1
globalys0=yp0-3
globalys1=yp1-3
globalyc0=y0-3
globalyc1=y1-3
globalid_=10
;id_=10*（cut+1）;useifshellsunconnected
command
;installpresupportconcrete
selshellid@id_grouprock&
rangecylend10@yp0-1.5end20@yp1-1.5rad7.4&
cylend10@yp0-1.5end20@yp1-1.5rad6.7not&
z-0.16
selshellpropisotropic10.5e9,0.25thickness0.3density2500&
rany@yp0@yp1
;excavatenextcut
modelmechnullrangegrouptunnely@y0@y1
modelmechnullrangegroup'concreteliner'y@y0@y1
;delete-cablesintheexcavatedarea
seldeletecablerangeid1y@y0@y1
seldeletecablerangeid2y@y0@y1
seldeletecablerangeid3y@y0@y1
end_command
localcut_1=cut-1
iidx=int（cut_1-3*（cut_1/3））+1
y2=min（y1+15,51）
inip（iidx）
ii=out（'CABLEBOLTPATTERN'+string（iidx））
command
seldeletecablerangeid@iidx
;installnewcables
selcableid@iidxbegin@x1@y1@z1end@x1@y2@z1nseg@cab_seg_m
selcableid@iidxbegin@x2@y1@z2end@x2@y2@z2nseg@cab_seg_m
selcableid@iidxbegin@x3@y1@z3end@x3@y2@z3nseg@cab_seg_m
selcableid@iidxbegin@x4@y1@z4end@x4@y2@z4nseg@cab_seg_m
selcableproemod45e9xcarea1.57e-3gr_per1.0&
yten25e4gr_k17.5e6gr_c20e4ranid@iidx
;shotcrete
selshellpropisotropic10.5e9,0.25thickness0.5density2500&
rany@ys0@ys1
end_command
ifcut>1then
command
;concreteliner
modelmechelrangegroup'concreteliner'y@yc0@yc1
propbulk@bmcsh@smcrangegroup'concreteliner'y@yc0@yc1
end_command
end_if
command
step3000
end_command
;storedisplacementsintables;��λ�ƴ����ڱ����
xtable（1,cut）=3.0*cut
ytable（1,cut）=gp_zdisp（ipt_surf）
xtable（2,cut）=3.0*cut
ytable（2,cut）=gp_zdisp（ipt_crown）
xtable（3,cut）=3.0*cut
ytable（3,cut）=gp_zdisp（ipt_spring）
command
savem1
end_command
ifcut=5then
command
savem1_15
end_command
end_if
ifcut=9then
command
savem1_27
end_command
end_if
ifcut=10then
command
savem1_30
end_command
end_if
end_loop
end
@excav
return
升级会员 