文档库 最新最全的文档下载
当前位置:文档库 › FLAC3D实例命令流1

FLAC3D实例命令流1

FLAC3D实例命令流1

第1部分命令流按照顺序进行

2-1定义一个FISH函数newdefabc

abc=25某3+5Endprintabc

2-2使用一个变量newdefabchh=25

abc=hh某3+5EndPrinthhPrintabc

2-3对变量和函数的理解newdefabchh=25

abc=hh某3+5End

etabc=0hh=0printhhprintabcprinthhnewdefabc

abc=hh某3+5end

ethh=25printabc

etabc=0hh=0printhhprintabcprinthh

2-4获取变量的历史记录new

genzonebrickize121modelmohr

prophear=1e8bulk=2e8cohe=1e5ten=1e10

fi某某yzrangey-0.10.1

applyyvel-1e-

5rangey1.92.1plotetrotation0045plotblockgroupdefget_ad ad1=gp_near(0,2,0)ad2=gp_near(1,2,0)ad3=gp_near(0,2,1)ad4=gp _near(1,2,1)endget_addefload

load=gp_yfunbal(ad1)+gp_yfunbal(ad2)+gp_yfunbal(ad3)+gp_yfun bal(ad4)endhitload

hitgpydi0,2,0tep1000

plothi1v-2

2-5用FISH函数计算体积模量和剪砌模量new

defderive

_mod=y_mod/(2.0某(1.0+p_ratio))

b_mod=y_mod/(3.0某(1.0-2.0某p_ratio))end

ety_mod=5e8p_ratio=0.25derive

printb_modprint_mod

2-6在FLAC输入中使用符号变量New

defderive

_mod=y_mod/(2.0某(1.0+p_ratio))

b_mod=y_mod/(3.0某(1.0-2.0某p_ratio))end

ety_mod=5e8p_ratio=0.25derive

genzonebrickize2,2,2modelelatic

propbulk=b_modhear=_modprintzonepropbulkprintzoneprophear

2-7控制循环Newdef某某某

um=0prod=1

loopn(1,10)um=um+nprod=prod某nend_loopend某某某

printum,prodnew

genzonebrickp0(0,0,0)p1(-10,0,0)p2(0,10,0)p3(0,0,-

10)modelela

plotetrotation0045plotblockgroupdefintall

pnt=zone_head

loopwhilepnt#nullz_depth=-z_zcen(pnt)

y_mod=y_zero+cc某qrt(z_depth)

z_prop(pnt,’hear’)=y_mod/(2.0某

(1.0+p_ratio))z_prop(pnt,’bulk’)=y_mod/(3.0某(1.0-2.0某

p_ratio))pnt=zone_ne某t(pnt)end_loopend

etp_ratio=0.25y_zero=1e7cc=1e8intall

2-8拆分命令行

new;e某ampleofaumofmanythingdeflong_um

temp=v1+v2+v3+v4+v5+v6+v7+v8+v9+v10long_um=temp+v11+v12+v13+ v14+v15end

2-9变量类型new

defhaveone

aa=2bb=3.4

cc=’Haveaniceday’

dd=aa某bb

ee=cc+’,oldchap’

endhaveoneprintfih

2-10IF条件语句newdefabc

if某某>0thenabc=1000ele

abc=-1000end_ifend

et某某=10printabcet某某=0printabc

2-11索单元自动生成new

genzonebrickize1035plotetrotation0045plotblockgroupdefplace_ cable

loopn(1,5)

place_cable

plotgridelgeomrednew

genzonebrickize1035plotetrotation15060plotblockgroupmodmohr propbulk1e8hear.3e8fric35propcoh1e3ten1e3iniden1000

7.0,1.5,z_dneg7etgrav0,0,-10

fi某某yzrangez-.1.1fi某yrangey-.1.1fi某yrangey2.93.1fi某某range某-.1.1fi某某range某9.910.1etlargehitunbalolve avecab_tr.av

ini某di0ydi0zdi0hitgp某dip0,1,5defplace_cableloopn(1,5)

free某range某-.1,.1zz_bz_tolve

place_cableavecab_end.av

plotketchelcableforcered

2-12圆形隧道开挖模拟计算;建立模型

genzonradcylp0000p1600p2010p3006&

ize4284dim3333rat1111.2groupoutiderockgenzonechellp0000p1300 p2010p3003&

ize1284dim2.72.72.72.7rat1111groupconcretlinerfillgroupinide rockgenzonreflectdip90dd90orig000genzonreflectdip0dd0ori000

genzonbrickp0006p1606p2016p30013ize426groupoutiderock1genzonbri ckp000-12p160-12p201-12p300-

6ize425groupoutiderock2genzonbrickp0600p12100p2610p3606ize1024g roupoutiderock3genzonreflectdip0dd0orig000rangegroupoutiderock3

Flac3D命令--完整经典版

实例分析命令: 1. X ,Y ,Z 旋转 Shift+ X ,Y ,Z 反向旋转 Gen zone ……;model ……;prop ……(材料参数);set grav 0,0,-9.81(重力加速度) plot add block group red yellow 把在group 中的部分染成红色和黄色 plot add axes black 坐标轴线为黑色;print zone stress% K 单元应力结果输出 ini dens 2000 ran z a b (设置初始密度,有时不同层密度不同);ini ……(设置初始条件);fix ……(固定界面) set plot jpg ;set plot quality 100 ;plot hard file 1.jpg 图像输出(格式、像素、名称) plot set magf 1.0视图的放大倍数为1.0;plo con szz z 方向应力云图 2. ini z add -1 range group one 群one 的所有单元,在z 方向上向下移动1m ;然后合并 命令 gen merge 1e-5 range z 0此命令是接触面单元合并成一个整体,1e-5是容差 3. (基坑开挖步骤):Step 1: create initial model state (建立初始模型)Step 2: excavate trench (开挖隧道) 4. group Top range group Base not 定义(群组Base 以外的为)群组Top 5. plot blo gro 使得各个群组不同颜色显示 6. (两个部分间设置界面;切割法):gen separate Top 使两部分的接触网格分离 为两部分;interface 1 wrap Base Top 在(Base 和Top )这两部分之间添加接触单元;plot create view_int 显示,并创建标题view_int ;plot add surface 显示表面;plot add interface red 界面颜色红色 7. (简单的定义函数及运行函数)new ;def setup 定义函数setup ;numy = 8定义常 量numy 为8;depth = 10.0 定义depth 为10;end 结束对函数的定义;setup 运行函数setup 8. (隧道生成)上部圆形放射性圆柱及下部块体单元体的建立,然后镜像。 9. 模拟模型的材料问题时为什么要去定义某个方向上的初始速度?— 10. 渐变应力施加:apply nstress -1e6 gradient 0,0,1e5 range z 3.464,0 plane dip 60 dd 270 origin .1 0 0;施加法向应力:apply nstress -1e6 range plane dip 60 dd 270 origin .1 0 0 11. d ip dd 确定平面位置使用:(纠结) 12. p rint gp position range id=14647 输出节点坐标 13. a pply sxx -10e6 gradient 0 , 0, 1e5 range z -100 , 0在这个求解方程中,z 为变量,所以xx σ为:65=-1010+10xx z σ?? ;原点(0,0,0) 14. f ree x range x -.1 .1 z 6.9 10.1放松x=0 平面上,z=7,10 这一部分在x 方向的约 束(可以在此处产生破坏) 15. 体积模量K 和剪切模量G 与杨氏模量及泊松比v 之间的转换关系如下: =3(1-2v)E K G=2(1+v) E 16. 一般而言,大多数问题可以采用FLAC 3D 默认的收敛标准(或称相对收敛标准),即当体 系最大不平衡力与典型内力的比率R 小于定值10-5;(也可由用户自定义该值,命令:

Flac3D 常用命令流

FLAC3D 常用命令流 ;模型镜像 gen zone radcylinder size 25 1 25 25 gen zone reflect normal -1 0 0 origin x y z(面上一点);沿X轴镜像,通过对称平面法线向量确定对称面 gen zone reflect normal 0 0 -1 ;沿z轴镜像 ;绘图控制 pl contour szz outline on ;在模型中显示位移-应变曲线 hist gp ydisp 0,0,0 hist zone syy 0,1,0 hist zone syy 1,1,0 pl his -2 -3 vs 1 ;在plot hist m vs n的形式里,m代表y轴,n代表x轴(不管m,n的正负); "-"表示对其值作"mirror" ;对模型进行压缩实验的方法 ;即在模型两侧施加相反方向的速度 ini yvel 1e-7 range y -.1 .1 ini yvel -1e-7 range y 1.9 2.1 ;修改模型的坐标值 ini x add -100 y add -100 z add -100 ;显示云图的同时也显示模型网格轮廓 plot add cont disp outline on ;gradient更精确 ;输入角度、弧度方法pi=π,90°为90.0*degrad def set_vals ptA = 25.0 * sin(pi/2);ptA=25.0 ptB = 25.0 * cos( 60.0*degrad );ptB=12.5 ptC=pi;ptC=3.1415926 end set_vals print ptA ptB ptC ;施加结构单元方法 sel shell id=5 range cylinder end1=(0.0, 0.0,0.0) & end2=(0.0,25.0,0.0) radius=24.5 not plot add sel geom black black cid on scale=0.03 sel node init zpos add -25.0 ;如何显示某一平面 plot create name_plane

FLAC3D快速入门及简单实例

FLAC3D快速入门 及简单实例 李佳宇编 LJY指南针教程

前言 FLAC及FLAC3D是由国际著名学者、英国皇家工程院院士、离散元的 发明人Peter Cundall博士在70年代中期开始研究的,主要面对岩土工程的通 用软件系统,目前已经在全球70多个国家得到广泛应用,在岩土工程学术界 和工业界赢得了广泛的赞誉。前国际岩石力学会主席 C.Fairhurst(1994)对 FLAC程序的评价是:“现在它是国际上广泛应用的可靠程序。” 我从研二(2010年)开始接触FLAC3D,最初的原因是导师要求每一个人至 少学会一个数值计算软件,而他嘴里每天念叨最多的就是FLAC,自己当时对数 值计算一无所知,便答应老师要学会FLAC3D。第一次打开软件界面,我心里 就凉了大半截,面对着一个操作界面跟记事本无异的所谓“功能强大”的岩土工 程专业软件,半点兴趣也提不起来。年底,从项目工地回到学校准备论文开题, 老师对我的开题报告非常不满意,当着全教研室师生的面,劈头盖脸大批一顿, 第二天又找谈话。在巨大的压力和强烈的自尊心驱使下,我硬着头皮开始啃 FLAC3D,一个半月之后,终于有了初步的计算结果,对老师有个交代,我也能 回家过年了。 前面这一段过程可能是大多数FLAC3D初学者的必经阶段,或者是即将 开始软件学习的人惧怕的事情。毫无疑问,FLAC3D极其不友好的界面是阻碍 初学者前进的很大障碍,当然还包括它是一个全英文的软件。但是当你费尽周折 的走进FLAC3D的世界,你就会发现它独特的魅力,比如简洁的界面,快捷的 命令流操作,高效的计算方法,不易报错等等。另外一个拿不上台面的优点就是 它非常小巧,包括Manual在内一共才几十兆大小,而且已经被破解成绿色版, 只要把它和命令流装进U盘,你就可以随便找一个身边功能最强大的电脑开始 计算了,如果你有过ANSYS、ABAQUS等大型软件痛苦的安装经历,你便能 毕业之后,本以为不用再接触数值计算,但工作需要使得我又一次开始与理解“绿色版”的含义,当然还请大家尊重知识产权,支持正版。 FLAC3D进行亲密接触。我的领导给了我很多新思路和很大的支持,如今我的 水平比研究生时有了不少提高。于是,我想把我的经验总结成文,希望对初学 者起

完整word版,FLAC3D 实例命令流1

第1部分命令流按照顺序进行 2-1定义一个FISH函数 new def abc abc = 25 * 3 + 5 End print abc 2-2使用一个变量 new def abc hh = 25 abc = hh * 3 + 5 End Print hh Print abc 2-3对变量和函数的理解 new def abc hh = 25 abc = hh * 3 + 5 End set abc=0 hh=0 print hh print abc print hh new def abc abc = hh * 3 + 5 end set hh=25 print abc set abc=0 hh=0 print hh print abc print hh 2-4获取变量的历史记录 new gen zone brick size 1 2 1 model mohr prop shear=1e8 bulk=2e8 cohes=1e5 tens=1e10

fix x y z range y -0.1 0.1 apply yvel -1e-5 range y 1.9 2.1 plot set rotation 0 0 45 plot block group def get_ad ad1 = gp_near(0,2,0) ad2 = gp_near(1,2,0) ad3 = gp_near(0,2,1) ad4 = gp_near(1,2,1) end get_ad def load load=gp_yfunbal(ad1)+gp_yfunbal(ad2)+gp_yfunbal(ad3)+gp_yfunbal(ad4) end hist load hist gp ydis 0,2,0 step 1000 plot his 1 vs -2 2-5用FISH函数计算体积模量和剪砌模量 new def derive s_mod = y_mod / (2.0 * (1.0 + p_ratio)) b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio)) end set y_mod = 5e8 p_ratio = 0.25 derive print b_mod print s_mod 2-6 在FLAC输入中使用符号变量 New def derive s_mod = y_mod / (2.0 * (1.0 + p_ratio)) b_mod = y_mod / (3.0 * (1.0 - 2.0 * p_ratio)) end set y_mod = 5e8 p_ratio = 0.25 derive gen zone brick size 2,2,2 model elastic prop bulk=b_mod shear=s_mod print zone prop bulk print zone prop shear

命令流共享-FLAC3D模型导入ANSYS中讲解学习

;********************** get nodetotal & zone total ********************* ;; Initialization def initialization IO_READ = 0 IO_WRITE = 1 IO_FISH = 0 IO_ASCII = 1 SuferFile='FLAC-ANSYS_MODEL.dat' array message(1) end initialization def writehead pnt1=zone_head pnt2=gp_head Z_total=0 P_total=0 loop while pnt1 # null Z_total=Z_total+1 pnt1=z_next(pnt1) endloop loop while pnt2 # null P_total=P_total+1 pnt2=gp_next(pnt2) endloop message(1)='!------------------------------------------------------------------ \n' message(1)='!This is a free soft for converting the model from Flac3d to Ansys! \n' message(1)='! Code by haitang 2006-05-30 \n' message(1)='!------------------------------------------------------------------ \n' message(1)='!The Number of Nodes and Elements \n' message(1)=message(1)+'!'+string(P_total)+' , '+string(Z_total)+'\n' message(1)=message(1)+'/prep7 \n' message(1)=message(1)+'ET,1,SOLID45' status = write(message,1) end define GetElements Z1_p=zone_head loop n (1,Z_total) Z1_code=z_code(Z1_p)

FLAC3D实例命令流1

FLAC3D实例命令流1 第1部分命令流按照顺序进行 2-1定义一个FISH函数newdefabc abc=25某3+5Endprintabc 2-2使用一个变量newdefabchh=25 abc=hh某3+5EndPrinthhPrintabc 2-3对变量和函数的理解newdefabchh=25 abc=hh某3+5End etabc=0hh=0printhhprintabcprinthhnewdefabc abc=hh某3+5end ethh=25printabc etabc=0hh=0printhhprintabcprinthh 2-4获取变量的历史记录new genzonebrickize121modelmohr prophear=1e8bulk=2e8cohe=1e5ten=1e10 fi某某yzrangey-0.10.1 applyyvel-1e- 5rangey1.92.1plotetrotation0045plotblockgroupdefget_ad ad1=gp_near(0,2,0)ad2=gp_near(1,2,0)ad3=gp_near(0,2,1)ad4=gp _near(1,2,1)endget_addefload

load=gp_yfunbal(ad1)+gp_yfunbal(ad2)+gp_yfunbal(ad3)+gp_yfun bal(ad4)endhitload hitgpydi0,2,0tep1000 plothi1v-2 2-5用FISH函数计算体积模量和剪砌模量new defderive _mod=y_mod/(2.0某(1.0+p_ratio)) b_mod=y_mod/(3.0某(1.0-2.0某p_ratio))end ety_mod=5e8p_ratio=0.25derive printb_modprint_mod 2-6在FLAC输入中使用符号变量New defderive _mod=y_mod/(2.0某(1.0+p_ratio)) b_mod=y_mod/(3.0某(1.0-2.0某p_ratio))end ety_mod=5e8p_ratio=0.25derive genzonebrickize2,2,2modelelatic propbulk=b_modhear=_modprintzonepropbulkprintzoneprophear 2-7控制循环Newdef某某某 um=0prod=1

相关FLAC3D程序命令流(--单桩承载力)

相关FLAC3D程序命令流(--单桩承载力)

相关FLAC3D程序命令流 单桩竖向受力模拟 new title 竖向加载桩 ;桩长11m,直径0.6m ;建立土体网格 gen zo radcyl p0 (0 0 0) p1 (8 0 0 ) p2 (0 0 -1.5) p3 (0 8 0) & P4 (8 0 -1.5) P5 (0 8 -1.5) P6 (8 8 0) P7 (8 8 -1.5) & P8 (0.3 0 0) P9 (0 0.3 0) P10 (0.3 0 -1.5) P11 (0 0.3 -1.5) & size 3 3 6 15 ratio 1 1 1 1.15 ;粉质粘土2 group fznt2 range z -1.5 0 gen zo radcyl p0 (0 0 -1.5) p1 (8 0 -1.5) p2 (0 0 -3.0) p3 (0 8 -1.5) & p4 (8 0 -3.0) p5 (0 8 -3.0) p6 (8 8 -1.5) p7 (8 8 -3.0) & p8 (0.3 0 -1.5) p9 (0 0.3 -1.5) p10 (0.3 0 -3.0) p11 (0 0.3 -3.0) & size 3 3 6 15 ratio 1 1 1 1.15 ;粉细砂 group fxs range z -3.0 -1.5 gen zo radcyl p0 (0 0 -3.0) p1 (8 0 -3.0) p2 (0 0 -5.5) p3 (0 8 -3.0) & p4 (8 0 -5.5) p5 (0 8 -5.5) p6 (8 8 -3.0) p7 (8 8 -5.5) & p8 (0.3 0 -3.0) p9 (0 0.3 -3.0) p10 (0.3 0 -5.5) p11 (0 0.3 -5.5) & size 3 5 6 15 ratio 1 1 1 1.15 ;粉质粘土4 group fznt4 range z -5.5 -3.0 gen zo radcyl p0 (0 0 -5.5) p1 (8 0 -5.5) p2 (0 0 -7.0) p3 (0 8 -5.5) & p4 (8 0 -7.0) p5 (0 8 -7.0) p6 (8 8 -5.5) p7 (8 8 -7.0) & p8 (0.3 0 -5.5) p9 (0 0.3 -5.5) p10 (0.3 0 -7.0) p11 (0 0.3 -7.0) & size 3 3 6 15 ratio 1 1 1 1.15 ;粉质粘土5 group fznt5 range z -7.0 -5.5 gen zo radcyl p0 (0 0 -7.0) p1 (8 0 -7.0) p2 (0 0 -8.5) p3 (0 8 -7.0) & p4 (8 0 -8.5) p5 (0 8 -8.5) p6 (8 8 -7.0) p7 (8 8 -8.5) & p8 (0.3 0 -7.0) p9 (0 0.3 -7.0) p10 (0.3 0 -8.5) p11 (0 0.3 -8.5) & size 3 3 6 15 ratio 1 1 1 1.15 ;粉土 group ft range z -8.5 -7.0 gen zo radcyl p0 (0 0 -8.5) p1 (8 0 -8.5) p2 (0 0 -11) p3 (0 8 -8.5) & p4 (8 0 -11) p5 (0 8 -11) p6 (8 8 -8.5) p7 (8 8 -11) & p8 (0.3 0 -8.5) p9 (0 0.3 -8.5) p10 (0.3 0 -11) p11 (0 0.3 -11) &

陈育民《FLAC3D基础与工程实例》全部命令流

《FLAC/FLAC 3D基础与工程实例》全部命令流 1 gen zone bri p0 0 0 0 p1 10 0 0 p2 0 10 0 p3 0 0 10 & p4 15 15 0 p5 0 15 15 p6 15 0 10 p7 20 20 20 & size 10 10 10 rat 1。0 0。9 1.1 group brick_1 gen zone bri p0 20 0 0 p1 add 10 0 0 p2 add 0 20 0 p3 add 0 0 15 & size 10 10 10 rat 1。0 0.9 1.1 group brick_2 gen zone bri p0 40 0 0 edge 10 size 10 10 10 rat 1。0 0。9 1.1 group brick_3 plot sur 2—1 new gen zon bri size 3 3 3 model elas prop bulk 3e8 shear 1e8 ini dens 2000 fix z ran z -.1 .1 fix x ran x —.1 .1 fix x ran x 2。9 3。1 fix y ran y —.1 。1 fix y ran y 2.9 3.1 set grav 0 0 -10 solve app nstress -10e4 ran z 3 x 1 2 y 1 2 hist gp vel 0 0 3 hist gp vel 0 3 3 plo hist 1 red plo add hist 2 blue solve 3—1 ;———----——----—--—-—-——-—-—-——--------工程信息 ;Project Record Tree export ;Title:Simple test ;-—-—-—-——-—---———-—---—-——--——--—计算第一步 ;。.. STATE: STATE1 ...。 config grid 10,10 model elastic group ’User:Soil' notnull model elastic notnull group 'User:Soil' prop density=1500.0 bulk=3E6 shear=1E6 notnull group ’User:Soil' fix x y j 1 fix x i 1 fix x i 11 set gravity=9。81 history 999 unbalanced solve save state1。sav ;---—--——-——----—--——————--———---—-计算第二步 ;.。。 STATE: STATE2 。。。。 initial xdisp 0 ydisp 0 initial xvel 0 yvel 0 model null i 4 7 j 8 10 group 'null' i 4 7 j 8 10 group delete ’null’ history 1 xdisp i=4, j=11 solve save state2.sav ;—-—--—————-——---—-——-———-—-———-—绘图命令 ;*** plot commands **** ;plot name: syy plot hold grid syy fill ;plot name: Unbalanced force plot hold history 999 ;plot name: grid plot hold grid magnify 20.0 lred grid displacement ;plot name: Xdis—A plot hold history 1 line 5—1 new ; =============================== ; 定义球体半径和半径方向上单元网格数 ; =============================== def parm rad=10.0 rad_size=5

(完整word版)flac命令流

1、FLAC3D常见命令: 1. FLAC3D是有限元程序吗? 答:不是!是有限差分法。 2。最先需要掌握的命令有哪些? 答:需要掌握gen, ini, app,plo,solve等建模、初始条件、边界条件、后处理和求解的命令. 3。怎样看模型的样子? 答:plo blo gro可以看到不同的group的颜色分布 4. 怎样看模型的边界情况? 答:plo gpfix red 5。怎样看模型的体力分布? 答:plo fap red 6。怎样看模型的云图? 答:位移:plo con dis (xdis,ydis,zdis) 应力:plo con sz (sy, sx, sxy, syz,sxz) 7。怎样看模型的矢量图? 答:plo dis (xdis,ydis, zdis) 8。怎样看模型有多少单元、节点? 答:pri info 9。怎样输出模型的后处理图? 答:File/Print type/Jpg file,然后选择File/Print,将保存格式选择为jpe文件 10. 怎样调用一个文件? 答:File/call或者call命令 10. 如何施加面力? 答:app nstress

11. 如何调整视图的大小、角度? 答:综合使用x, y, z, m,Shift键,配合使用Ctrl+R,Ctrl+Z等快捷键 12. 如何进行边界约束? 答:fix x ran (约束的是速度,在初始情况下约束等效于位移约束) 13. 如何知道每个单元的ID? 答:用鼠标双击单元的表面,可以知道单元的ID和坐标 14。如何进行切片? 答:plo set plane ori (点坐标) norm (法向矢量) plo con sz plane (显示z方向应力的切片) 15. 如何保存计算结果? 答:save +文件名. 16. 如何调用已保存的结果? 答:rest +文件名;或者File / Restore 17。如何暂停计算? 答:Esc 18。如何在程序中进行暂停,并可恢复计算? 答:在命令中加入pause命令,用continue进行继续 19. 如何跳过某个计算步? 答:在计算中按空格键跳过本次计算,自动进入下一步 20。Fish是什么东西? 答:是FLAC3D的内置语言,可以用来进行参数化模型、完成命令本身不能进行的功能21. Fish是否一定要学? 答:可以不用,需要的时候查Mannual获得需要的变量就可以了

FLAC-3D深基坑的开挖与支护的命令流

FLAC D3深基坑的开挖与支护的命令流 一、实例工程 南宁地区地层属于河流阶地二元地层,广泛分布有较厚的圆砾层,国内尚无在类似地层条件下建设地铁基坑的经验,为此,可使用FLAC3D 对基坑开挖的全过程进行三维数值模拟,在对比实测数据的基础上,总结圆砾层中地铁车站深基坑的地下连续墙水平变形及周围地表沉降变形特征。 该基坑位于大学路与明秀路交叉路口处,沿大学东路东西向布置。车站基坑长465m,标准断面宽度为20.7m,为地下两层式结构,底板埋深为15.535m(相对地面),顶板覆土厚度大于3m。本工程主体建筑面积21163.6m2,主要结构形式为双柱三跨框架箱型结构。本工程所处的大学路为南宁市东西向的主要交通枢纽,车流量大,人流密集,地面条件复杂。 基坑施工采用明挖顺作法施工,围护结构为800mm厚地下连续墙+内撑(三道内支撑加一道换撑)的支护体系。第一道支撑采用钢筋混凝土支撑,尺寸为800×900mm,冠梁同时作为第一道钢筋混凝土支撑的围檩。第二、三道支撑及换撑使用钢支撑并施加预加力,直径为609mm,壁厚为t=16mm,斜撑段采用800×1000mm钢筋砼腰梁,其余为2×I45C 钢围檩。 二、模型建立 建模工作由两部分组成,实体模型部分,包括土体和地下连续墙;结构单元部分,包括混凝土支撑和钢支撑。 根据对称性原理,拟选取1/2 的实际工程尺寸进行分析。考虑到实际的基坑长度将近500m,根据以往的经验,选取全部长度的一半虽然能够得到满意的结果,但是由于中间部分的基坑基本处于同样的受力状态,这样会使大部分的计算长度变为重复的计算,降低了计算效率。根据初步计算结果和经验,最终确定的基坑尺寸为,宽度取基坑的最大宽度24m,开挖深度19m,基坑长度36m。根据地勘报告,合并相似土层,模型中共划分了7个土层。 在FLAC3D 中,围护结构可以用衬砌单元(liner)或实体单元模拟。根据Zdravdovi的研究,在二维平面基坑模拟中,分别采用实体单元和梁单元(相当于三维模型中的衬砌单元)计算所产生的墙体变形差别小于4%,而引起地表沉降的主要原因是围护结构变形造成的地层损失,可见上述两种方法计算结果的差别可忽略不计。相比于衬砌单元(liner),实体单元物理模型清晰,参数相对较少且容易确定,因此本例围护结构采用实体单元模拟。 模型的范围必须要足够大以覆盖基坑开挖可能影响的区域,同时又要兼顾计算效率。根据现行国家和个地方规范的要求,参考已有的研究成果,影响范围为4 倍开挖深度,整个模型的尺寸为96×144×39(长×宽×高m)。为减少不必要的网格,网格划分原则为基坑附近密集,远处稀疏,模型共有节点10890个,单元9408个。模型如图所示。

FLAC3D_施加锚杆命令流

new gen zone radtun p0 0,0,0 p1 25,0,0 p2 0,50, 0 p3 0,0,25 size 4 25 4 10 dim 4 4 4 4 ratio 1 1 1 1.1 fill gen zone reflect normal 1 0 0 ori 0 0 0 gen zone reflect normal 0 0 1 ori 0 0 0 mo mohr pro bulk 2.2e9 she 1.3e9 fric 30 coh 1.3e6 ten 1.5e5 ini dens 2000 fix x range x -25.1,-24.9 fix x range x 24.9 25.1 fix y range y 49.9 50.1 fix z range z -25.1 -24.9 fix z range z 24.9 25.1 PLOT SURF PAUSE K sel cable id=1 beg 0, 0, 0 end 0 ,29, 0 nseg 10 sel cable id=1 beg 0,29,0 end 0,35,0 nseg 6 sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 & gr_coh 1 gr_k 1 gr_per 0.0785 range cid 1,10 sel cable id=1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 & gr_coh 10e5 gr_k 2e7 range cid 11,17 sel delete link range id 1 sel link id=100 1 target zone sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id 100 sel cable id=1 pretension 60e3 range cid 1,10 set mech ratio 7.0e-3 solve SAVE 2CALE.SAV plot set plane ori 0 5 0 norm 0 0 1 plot cont syy plane ou on plot add axe plot set rot 90 0 0 plot set mag 2

FLAC3D命令流(挺实用)

FLAC3D命令流(挺实用) 1、怎样查看模型? 答:plot grid 可以查看网格,plot grid num 可以查看节点号。 2、请问在圆柱体四周如何施加约束条件? 答:可以用fix ... ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not,其中r2 fix x range end1 1 0 0 end2 1 4 0 rad 1 end1 1 0 0 end2 1 4 0 rad 1 3、怎么能把一个PLOT的图像数据导出来以便用其他软件绘图?答:用set log on 命令,把数据导出来,转到excel里处理一下,然后用surfer或者什么作图软件绘制就行了。 4、用命令建立模型后,如何显示点的坐标? 答:使用plo blo gro gpnum on命令 5、关于gauss_dev对性质进行高斯正态分布的问题? 答:根据手册上的说明:下面的命令设定一个平均摩擦角为40度,标准方差是±5%。则命令如下:prop friction 40 gauss_dev 2 问题:请问gauss_dev 2中的2是如何计算的?如果把±5%改为±10%,则命令应如何写?40×5%=2 6、reflect问题 问:gen zone radbrick &p0 (0,0,0) p1 (10,0,0) p2 (0,10,0) p3 (0,0,10) &size 3,5,5,7 &ratio 1,1,1,1.5 &dim 1 4 2 fill plot surf gen zone reflect dip 0 dd 90 (对xy面做镜像) gen zone reflect dip 90 dd 90 (对yz面做镜像) (1)dd表示y轴正向顺时针到那条射线的夹角,dip表示对称参照面与xy平面的夹角,对称参照面与xy平面的夹角在xy平面的投影

FLAC3D常用命令

1. apply(缩写:app) 可用来定义边界条件及初始条件: 1)添加应力 格式1:apply szz -0.3395e6 range z -0.1 0.1 格式1:apply szz -0.3395e6 range group pile 格式3:apply szz -0.3395e6 range z -0.1 0.1 group pile 格式4:apply nstress 数值range z 2.9 3.1(或3)x 1 2 y 1 2 2)以一定速度施加位移边界 格式1:apply yvel -1e5 range y -1.9 2.1 ;施加y方向速度-1e5/step 3)添加边界条件 格式:apply szz <数值> grad <梯度> range <范围> 示例1:apply szz -1e9 grad 0 0 8.3e5 range z 0 120 示例2:apply szz -0.6e6 range z 0.05 0.15 group pile 注:<数值>是梯度方向坐标0点的数值,可通过定义坐标范围的上下值与梯度计算得到。 2. range(缩写:ran) 通过range功能,可以使命令作用在一定指定范围的目标上;如果一个命令没有使用range来确定范围,则命令对整个模型有效。 1)利用坐标指定一定的范围 格式1:range z 0 1 格式2:range z 2.9 3.1 x 1 2 y 1 2 2)利用分组来指定范围 格式:range group 1 3)以上两种的复合 格式:range z -0.1 0.1 group pile 4)利用id号来指定一定的范围 格式:range id 0 10 该命令后跟起始id和结束id,这里的id可以是实体单元、网格、结构单元、

FLAC3D命令流(整理版)

精心整理1、怎样查看模型? 答:plotgrid可以查看网格,plotgridnum可以查看节点号。 2、请问在圆柱体四周如何施加约束条件? 答:可以用fix...rancylinderend1end2radiusr1cylinderend1end2radiusr2not,其中r2

flac3D模拟地面曲面方程命令流

new gen zone brick p0 0,167.0072,0 p1 0,0,0 p2 242.6586,167.0072,0 p3 0,167.0072,32.0 size 56,81,11 ;;定义地面曲面方程 DEF eqn537 ARRAY c(45),cx(12),cy(12),v(70) c(1)=30.33465961749452 c(2)=0.3768529280014758 c(3)=-1.411874011262434 c(4)=-0.4466983009380354 c(5)=0.02498771591159658 c(6)=-0.6225933592839830 c(7)=-0.5572409253529970 c(8)=0.2594590432443336 c(9)=-0.2154254299775080 c(10)=0.05611608056180144 c(11)=-0.3573754488804990 c(12)=0.02174566734976303 c(13)=0.4051641051840550 c(14)=0.1744265369646383 c(15)=0.1147014441688443 c(16)=-0.2088509864184064 c(17)=0.1493365036211442 c(18)=0.4311419368791974 c(19)=0.008961362272148653 c(20)=-0.07774920998867369 c(21)=0.1811142002504400 c(22)=-0.1519527001982752 c(23)=0.07915472779033880 c(24)=0.2559732910968667 c(25)=0.1026065146630541 c(26)=-0.2819625969135452 c(27)=-0.09027129904676756 c(28)=0.1130183438471391 c(29)=-0.02169639756233838 c(30)=0.1907291399047935 c(31)=0.1297366637584324 c(32)=0.08806739336824376 c(33)=-0.2260336878126227 c(34)=-0.2115733057775976 c(35)=0.1148217058529207 c(36)=-0.06508949301229160 c(37)=-0.04511427839772801 c(38)=0.1632228872577737

FLAC3D命令流挺实用

1、怎样查瞧模型? 答:plot grid 可以查瞧网格,plot grid num 可以查瞧节点号。 2、请问在圆柱体四周如何施加约束条件? 答:可以用fix 、、、ran cylinder end1 end2 radius r1 cylinder end1 end2 radius r2 not,其中r2

相关文档
相关文档 最新文档