flac3d5.0软件隧道支护与开挖命令流.txt资料文档下载

flac3d5.0软件隧道支护与开挖命令流.txt资料文档下载

《flac3d5.0软件隧道支护与开挖命令流.txt资料文档下载》由会员分享，可在线阅读，更多相关《flac3d5.0软件隧道支护与开挖命令流.txt资料文档下载（6页珍藏版）》请在冰豆网上搜索。

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'

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;

��ÿ

�˵�λ�ý

��и�ֵ

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）

@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

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）

@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

;

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'

y@y0@y1

delete-cablesintheexcavatedarea

seldeletecablerangeid1y@y0@y1

seldeletecablerangeid2y@y0@y1

seldeletecablerangeid3y@y0@y1

localcut_1=cut-1

iidx=int（cut_1-3*（cut_1/3））+1

y2=min（y1+15,51）

ii=out（'

CABLEBOLTPATTERN'

+string（iidx））

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

yten25e4gr_k17.5e6gr_c20e4ranid@iidx

shotcrete

selshellpropisotropic10.5e9,0.25thickness0.5density2500&

rany@ys0@ys1

ifcut>

1then

command

;

concreteliner

modelmechelrangegroup'

y@yc0@yc1

propbulk@bmcsh@smcrangegroup'

end_command

end_if

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）

savem1

ifcut=5then

command

savem1_15

end_command

ifcut=9then

savem1_27

ifcut=10then

savem1_30

@excav

return