整理用flac3d模拟基坑开挖Word格式文档下载.docx

整理用flac3d模拟基坑开挖Word格式文档下载.docx

《整理用flac3d模拟基坑开挖Word格式文档下载.docx》由会员分享，可在线阅读，更多相关《整理用flac3d模拟基坑开挖Word格式文档下载.docx（15页珍藏版）》请在冰豆网上搜索。

genmerge1e-4rangex180202y04z-30。

1—29。

9

attachfacerangex89.9990。

01y0.04。

0z-29.90

attachfacerangex109.99110。

01y0。

04。

0z-29。

90

attachfacerangex179。

99180.01y0。

04.0z-29.90

attachfacerangex201.99202。

01y0.04.0z-29.90

sellinerid1crossdiaggroup2rangex89.990。

1y—0.14.1z—30。

10.1

sellinerid2crossdiaggroup2rangex109。

9110.1y—0.14。

1z—30.10.1

sellinerid3crossdiaggroup2rangex179.9180。

1y—0。

14。

1z-30.10。

sellinerid4crossdiaggroup2rangex201。

9202。

1y-0。

1z—30。

10。

sellinerid1propisotropic（2。

0e10,0.20）thickness0.8density2.5e3&

cs_nk=4e9cs_sk=4e9&

cs_ncut=4e7cs_scoh=4e7cs_scohres=0cs_sfric=20.0＆

rangex89。

990。

sellinerid2propisotropic（2。

0e10,0.20）thickness0。

8density2.5e3&

cs_nk=4e9cs_sk=4e9＆

cs_ncut=4e7cs_scoh=4e7cs_scohres=0cs_sfric=20。

0&

rangex109.9110。

14.1z-30.10.1

sellinerid3propisotropic（2.0e10，0。

20）thickness0.8density2.5e3&

rangex179。

9180.1y—0.14。

1z—30.10。

sellinerid4propisotropic（2。

0e10，0。

20）thickness0.8density2.5e3＆

cs_ncut=4e7cs_scoh=4e7cs_scohres=0cs_sfric=20.0&

rangex201。

1z-30。

定义支撑结构

defstruct_install1

loopi（1，3）

structx_zz=—1。

0＊5.0＊（i—1）

structx_xx0=90。

0

structx_xx1=110。

structx_yy=2.0

command

selbeamid=2begin（structx_xx0,structx_yy，structx_zz）end（structx_xx1,structx_yy,structx_zz）nseg=10

selbeamid=2propdens=0。

000emod=1.0e—6nu=0.2&

xcarea=0.80xcj=10.94e-2xciy=6。

67e-2xciz=4.27e—2ydirection=（00—1）;

1000x800

endcommand

endloop

end

struct_install1

defstruct_install2

structx_zz=-1。

0＊5。

0＊（i—1）

structx_xx0=180.0

structx_xx1=202。

selbeamid=3begin（structx_xx0，structx_yy,structx_zz）end（structx_xx1，structx_yy，structx_zz）nseg=11

selbeamid=3propdens=0。

000emod=1.0e-6nu=0.2&

xcarea=0.80xcj=10.94e-2xciy=6.67e-2xciz=4.27e-2ydirection=（00-1）;

struct_install2

建立结构单元分组

selgrouplinerwallrangeselliner

selgroupstruct1rangeselbeamx（90.0110.0）z（-0。

1）

selgroupstruct2rangeselbeamx（90。

0110.0）z（—5.1—4。

9）

selgroupstruct3rangeselbeamx（90.0110.0）z（—10.1—9.9）

selgroupstruct4rangeselbeamx（180。

0202.0）z（-0.10.1）

selgroupstruct5rangeselbeamx（180。

0202.0）z（—5.1-4。

9）

selgroupstruct6rangeselbeamx（180.0202。

0）z（-10。

1-9。

删除beam单元的link

seldelelinkrangeselbeamz（—300）

建立liner间的节点间的刚性link

defmerge_link0

node_num=0

node_pnt0=nd_head

loopwhilenode_pnt0＃null；

寻找总节点数，注：

不能自己任生成node,程序缺省的方式为连续生成无不连续

node_num=node_num+1

node_pnt0=nd_next（node_pnt0）

node_num_minus1=node_num-1

link_id=30000

loopii（1，node_num_minus1）

node_pnt1=nd_find（ii）

xxa=nd_pos（node_pnt1，2,1）

yya=nd_pos（node_pnt1,2，2）

zza=nd_pos（node_pnt1，2,3）

ii_plus1=ii+1

loopjj（ii_plus1，node_num）

node_pnt2=nd_find（jj）

xxb=nd_pos（node_pnt2,2,1）

yyb=nd_pos（node_pnt2，2,2）

zzb=nd_pos（node_pnt2，2,3）

node_dist=sqrt（（xxa-xxb）^2+（yya—yyb）^2+（zza-zzb）^2）

dist_tol=1e-1

ifnode_dist<

=dist_tolthen

link_pnt1=nd_link（node_pnt1）

link_pnt2=nd_link（node_pnt2）

；

iflink_pnt1#nullthen

;

temp1=lk_delete（link_pnt1）

endif

iflink_pnt2＃nullthen

temp2=lk_delete（link_pnt2）

endif

link_id=link_id+1

command；

生成新link（6自由度全固结），大的node的id作为targetnode,小的node的id作为sourcenode，需注意不同情况下的灵活调整

selsetlinknode_tol=dist_tol

sellinkid=link_idjjtarget=nodetgt_num=ii；

指定link的ID

sellinkiitarget=nodetgt_num=jj；

不指定link的id，自动生成

sellinkattachxdir=rigidydir=rigidzdir=rigidxrdir=rigidyrdir=rigidzrdir=rigidrangeid=link_id

endcommand

endif

merge_link0

设置土层材料参数

defb_s_mod

b_mod=e_mod/（3.0＊（1.0—2。

0＊p_ratio））

s_mod=e_mod/（2.0*（1.0+p_ratio））

modelelastic

sete_mod100e6

setp_ratio0.3

b_s_mod

propbu=b_modsh=s_mod

inidens1800rangez-750

defini_szz

szz0=0

szzgrad=1800*10

iniszzaddszz0grad00szzgradrangez—750

ini_szz

defini_sxx_syy

pnt=zone_head

loopwhilepnt＃null

val=k0*z_szz（pnt）

z_sxx（pnt）=val

z_syy（pnt）=val

pnt=z_next（pnt）

setk0=0。

50

ini_sxx_syy

定义边界处的结构边界条件

cyc0

selnodelocalxdir=（0,1，0）ydir=（0,0,1）rangex89。

990.1y—0。

14.1z—30.10.1

selnodelocalxdir=（0，1，0）ydir=（0,0,—1）rangex109。

9110。

1y—0.14。

selnodelocalxdir=（0,1，0）ydir=（0,0，1）rangex179.9180.1y—0.14。

1z-30.10.1

selnodelocalxdir=（0，1，0）ydir=（0，0，—1）rangex201。

selnodefixlsysrangex89.990。

10.1z-30。

selnodefixlsysrangex89.990.1y3。

94.1z-30.10.1

selnodefixlsysrangex109.9110。

10.1z—30.10。

1y3.94。

selnodefixlsysrangex179.9180.1y-0.10。

selnodefixlsysrangex179。

9180。

1y3。

94.1z-30.10。

selnodefixlsysrangex201。

9202.1y—0.10.1z-30.10.1

94.1z-30。

selnodefixxyrzrrangex89.990.1y-0.10。

selnodefixxyrzrrangex89.990.1y3.94。

selnodefixxyrzrrangex109。

9110.1y-0。

selnodefixxyrzrrangex109.9110.1y3.94.1z—30。

selnodefixxyrzrrangex179.9180。

94。

selnodefixxyrzrrangex201.9202。

1y—0.10。

selnodefixyrangex89。

1y0.04.0z-0.10.1

selnodefixyrangex109。

9110.1y0.04.0z-0.10.1

selnodefixyrangex179。

1y0。

04.0z—0。

selnodefixyrangex201.9202.1y0.04.0z-0.10。

setplotmeta

plotsetrot20030bawhcolor=oncent=（10200）mag=3.81

setoutpnode_local_sys。

wmf

plotaddselgeomblackredlink=offnode=offid=offshrink=0scale=0。

03nodesys=onrangegrouplinerwallanygroupstruct1any

plha

固定边界条件

fixxrangex-0。

fixxrangex291。

9292.1

fixyrangey—0。

fixyrangey3。

fixxyzrangez—75.1-74.9

setgrav0，0，-10

solve

saveelas.sav

删除侧面内外土体间的连接约束

attachdeleterangex89。

9990.01y0。

04.0z-29。

attachdeleterangex109.99110.01y0。

0z—29。

attachdeleterangex179。

99180.01y0.04。

attachdeleterangex201.99202.01y0。

0z—29.90

在墙内土体的外侧建立接触面

interface1facerangegroup1x89。

9990。

01y0.04.0z-29。

interface2facerangegroup1x109.99110。

interface3facerangegroup1x179。

99180。

interface4facerangegroup1x201。

99202。

interface1propkn=4e8ks=4e8tens=5e3coh=0。

0fric=20；

接触面参数

interface2propkn=4e8ks=4e8tens=5e3coh=0.0fric=20；

interface3propkn=4e8ks=4e8tens=5e3coh=0。

0fric=20;

interface4propkn=4e8ks=4e8tens=5e3coh=0。

interface1maxedge=1

interface2maxedge=1

interface3maxedge=1

interface4maxedge=1

interface1propkn=4e8ks=4e8tens=1e10sbratio=100

plotsetbawh

plskeinterfaceredblueattachcyangreen

setoutpinterface_attachment。

重新定义连续墙参数

0e10,0。

20）＆

cs_scoh=4e7cs_scohres=0.0cs_sfric=0.0rangex89。

sellinerid2propisotropic（2.0e10,0.20）＆

cs_scoh=4e7cs_scohres=0.0cs_sfric=0.0rangex109.9110.1y-0。

14.1z-30。

sellinerid3propisotropic（2。

0e10，0.20）＆

cs_scoh=4e7cs_scohres=0.0cs_sfric=0.0rangex179。

9180.1y-0.14.1z—30.10.1

sellinerid4propisotropic（2.0e10,0.20）&

cs_scoh=4e7cs_scohres=0.0cs_sfric=0。

0rangex201。

9202.1y-0。

重新定义墙底约束条件