完整word版用flac3d模拟基坑开挖.docx

完整word版用flac3d模拟基坑开挖.docx

《完整word版用flac3d模拟基坑开挖.docx》由会员分享，可在线阅读，更多相关《完整word版用flac3d模拟基坑开挖.docx（14页珍藏版）》请在冰豆网上搜索。

完整word版用flac3d模拟基坑开挖

new

;网格建立

;;;;;;;;;;;;;;;;;;;;;;;;;

genzonebrickp0900-30p12020-30p2904-30p39000size112430ratio111

genzonebrickp0900-30p19000p2904-30p300-30size30425ratio111.1

genzonebrickp0900-30p100-30p2904-30p3900-75size25418ratio1.111.1

genzonebrickp0900-30p1900-75p2904-30p32020-30size184112ratio1.111

genzonebrickp02020-30p12920-30p22024-30p320200size25430ratio1.111

genzonebrickp02020-30p12020-75p22024-30p32920-30size18425ratio1.111.1

;分组

;;;;;;;;;;;;;;;;;;;;;;;;;;

group1rangex90110y04z-300

group1rangex180202y04z-300

group2rangegroup1not

;建立连续墙单元

;;;;;;;;;;;;;;;;;;;;;;;;;;

genseparate1

genmerge1e-4rangex90110y04z-30.1-29.9

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.1

sellinerid4crossdiaggroup2rangex201.9202.1y-0.14.1z-30.10.1

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.1y-0.14.1z-30.10.1

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.1y-0.14.1z-30.10.1

sellinerid3propisotropic（2.0e10,0.20）thickness0.8density2.5e3&

cs_nk=4e9cs_sk=4e9&

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

rangex179.9180.1y-0.14.1z-30.10.1

sellinerid4propisotropic（2.0e10,0.20）thickness0.8density2.5e3&

cs_nk=4e9cs_sk=4e9&

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

rangex201.9202.1y-0.14.1z-30.10.1

;定义支撑结构

;;;;;;;;;;;;;;;;;;;;;;;

defstruct_install1

loopi（1,3）

structx_zz=-1.0*5.0*（i-1）

structx_xx0=90.0

structx_xx1=110.0

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

loopi（1,3）

structx_zz=-1.0*5.0*（i-1）

structx_xx0=180.0

structx_xx1=202.0

structx_yy=2.0

command

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）;1000x800

endcommand

endloop

end

struct_install2

;建立结构单元分组

;;;;;;;;;;;;;;;;;;;;;;;;;;;

selgrouplinerwallrangeselliner

selgroupstruct1rangeselbeamx（90.0110.0）z（-0.10.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.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）

endloop

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

endloop

endloop

end

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））

end

modelelastic

sete_mod100e6

setp_ratio0.3

b_s_mod

propbu=b_modsh=s_mod

inidens1800rangez-750

defini_szz

szz0=0

szzgrad=1800*10

command

iniszzaddszz0grad00szzgradrangez-750

endcommand

end

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）

endloop

end

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.1z-30.10.1

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.9202.1y-0.14.1z-30.10.1

selnodefixlsysrangex89.990.1y-0.10.1z-30.10.1

selnodefixlsysrangex89.990.1y3.94.1z-30.10.1

selnodefixlsysrangex109.9110.1y-0.10.1z-30.10.1

selnodefixlsysrangex109.9110.1y3.94.1z-30.10.1

selnodefixlsysrangex179.9180.1y-0.10.1z-30.10.1

selnodefixlsysrangex179.9180.1y3.94.1z-30.10.1

selnodefixlsysrangex201.9202.1y-0.10.1z-30.10.1

selnodefixlsysrangex201.9202.1y3.94.1z-30.10.1

selnodefixxyrzrrangex89.990.1y-0.10.1z-30.10.1

selnodefixxyrzrrangex89.990.1y3.94.1z-30.10.1

selnodefixxyrzrrangex109.9110.1y-0.10.1z-30.10.1

selnodefixxyrzrrangex109.9110.1y3.94.1z-30.10.1

selnodefixxyrzrrangex179.9180.1y-0.10.1z-30.10.1

selnodefixxyrzrrangex179.9180.1y3.94.1z-30.10.1

selnodefixxyrzrrangex201.9202.1y-0.10.1z-30.10.1

selnodefixxyrzrrangex201.9202.1y3.94.1z-30.10.1

selnodefixyrangex89.990.1y0.04.0z-0.10.1

selnodefixyrangex109.9110.1y0.04.0z-0.10.1

selnodefixyrangex179.9180.1y0.04.0z-0.10.1

selnodefixyrangex201.9202.1y0.04.0z-0.10.1

;setplotmeta

;plotsetrot20030bawhcolor=oncent=（10200）mag=3.81

;setoutpnode_local_sys.wmf

;plotaddselgeomblackredlink=offnode=offid=offshrink=0scale=0.03nodesys=onrangegrouplinerwallanygroupstruct1any

;plha

;固定边界条件

;;;;;;;;;;;;;;;;;;;;;;;;;;

fixxrangex-0.10.1

fixxrangex291.9292.1

fixyrangey-0.10.1

fixyrangey3.94.1

fixxyzrangez-75.1-74.9

setgrav0,0,-10

solve

saveelas.sav

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

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

attachdeleterangex89.9990.01y0.04.0z-29.90

attachdeleterangex109.99110.01y0.04.0z-29.90

attachdeleterangex179.99180.01y0.04.0z-29.90

attachdeleterangex201.99202.01y0.04.0z-29.90

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

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

interface1facerangegroup1x89.9990.01y0.04.0z-29.90

interface2facerangegroup1x109.99110.01y0.04.0z-29.90

interface3facerangegroup1x179.99180.01y0.04.0z-29.90

interface4facerangegroup1x201.99202.01y0.04.0z-29.90

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.0fric=20;接触面参数

interface1maxedge=1

interface2maxedge=1

interface3maxedge=1

interface4maxedge=1

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

;plotsetbawh

;plskeinterfaceredblueattachcyangreen

;setoutpinterface_attachment.wmf

;plha

;重新定义连续墙参数

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

sellinerid1propisotropic（2.0e10,0.20）&

cs_nk=4e9cs_sk=4e9&

cs_scoh=4e7cs_scohres=0.0cs_sfric=0.0rangex89.990.1y-0.14.1z-30.10.1

sellinerid2propisotropic（2.0e10,0.20）&

cs_nk=4e9cs_sk=4e9&

cs_scoh=4e7cs_scohres=0.0cs_sfric=0.0rangex109.9110.1y-0.14.1z-30.10.1

sellinerid3propisotropic（2.0e10,0.20）&

cs_nk=4e9cs_sk=4e9&

cs_scoh=4e7cs_scohres=0.0cs_sfric=0.0rangex179.9180.1y-0.14.1z-30.10.1

sellinerid4propisotropic（2.0e10,0.20）&

cs_nk=4e9cs_sk=4e9&

cs_scoh=4e7cs_scohres=0.0cs_sfric=0.0rangex201.9202.1y-0.14.1z-30.10.1

;重新定义墙底约束条件

;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;

defredef_wall_end_link1

node_pnt=nd_head

link_id=100000

loopwhilenode_pnt#null

node_id=nd_id（node_pnt）

xx=nd_pos（node_pnt,2,1）

yy=nd_pos（node_pnt,2,2）

zz=nd_pos（node_pnt,2,3）

link_pnt=nd_link（node_pnt）

dist_x=sqrt（（xx-90.0）^2+（zz+30.0）^2）

ifdist_x<=dist_tolthen

iflink_pnt#nullthen

temp1=lk_delete（link_pnt）\

link_id=link_id+1

command

selsetlinknode_tol=dist_tol

sellinkid=link_idnode_idtargetzone

sellinkattachxdir=rigidydir=rigidzdir=rigidxrdir=f