完整word版用flac3d模拟基坑开挖Word下载.docx
- 文档编号:19732234
- 上传时间:2023-01-09
- 格式:DOCX
- 页数:14
- 大小:18.26KB
完整word版用flac3d模拟基坑开挖Word下载.docx
《完整word版用flac3d模拟基坑开挖Word下载.docx》由会员分享,可在线阅读,更多相关《完整word版用flac3d模拟基坑开挖Word下载.docx(14页珍藏版)》请在冰豆网上搜索。
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&
rangex109.9110.1y-0.14.1z-30.10.1
sellinerid3propisotropic(2.0e10,0.20)thickness0.8density2.5e3&
rangex179.9180.1y-0.14.1z-30.10.1
sellinerid4propisotropic(2.0e10,0.20)thickness0.8density2.5e3&
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
structx_xx0=180.0
structx_xx1=202.0
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&
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)
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.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
重新定义连续墙参数
sellinerid1propisotropic(2.0e10,0.20)&
cs_scoh=4e7cs_scohres=0.0cs_sfric=0.0rangex89.990.1y-0.14.1z-30.10.1
sellinerid2propisotropic(2.0e10,0.20)&
cs_scoh=4e7cs_scohres=0.0cs_sfric=0.0rangex109.9110.1y-0.14.1z-30.10.1
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.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
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 完整 word 版用 flac3d 模拟 基坑 开挖