flac3d工程实例

合集下载

FLAC3D基本原理及简单实例

FLAC3D基础知识
• 其中，体积模量K和剪切模量G与杨氏模量E和泊松比v有以下关系：
E 3(1 2 ) E G 2(1 ) K
9 KG 3K G 3K 2G G 2(3K G ) E
或
摩尔-库伦塑性模型需要材料参数有：（1）密度（2）体积模量（3）剪切模量（4）内摩擦角（5）粘聚力（6）抗拉强度如果不指定这些材料参数，其值将会自动默认为零。
3D
生成网格
执行变更
定义材料本构关系和性质定义边界、初始条件
计算结果保存及调用
图形绘制及结果输出
FLAC3D基础知识
指定材料模型
• 一旦完成了网格的生成，就必须给模型中的所有单元指定一种或者更多的材料模型及相应的性质。这可以用两个命令MODEL和 PROPERTY来完成。FLAC中有十种内置的材料模型，一般只用三种模型：MODEL null，MODEL elastic和MODEL mohr。 • MODEL null指的是从模型中去除的或开挖的材料； MODEL elastic 指的是各向同性弹性材料行为； MODEL mohr指的是摩尔-库伦塑性行为。 • MODEL elastic和MODEL mohr需要通过PROPERTY命令指定材料的性质，弹性模型需要的材料参数有： • （1）密度 • （2）体积模量 • （3）剪切模量
f t 3 t
式中，是摩擦角，C是粘聚力， t 是张拉强度，且有：
N
3
张拉强度不超过值，最大值由下式给定：
1 sin 1 sin
t max
c tan
2.2 FLAC3D常用材料本构模型
Mohr-Coulomb模型
流动法则

FLAC及FLAC3D基础与工程实例51-2

FLAC及FLAC3D基础与工程实例51-2FLAC及FLAC3D基础与工程实例51第1章FLAC、FLAC3D的功能与特性自R.W；数值模拟技术的优势在于有效延伸和扩展了分析人员的；本章重点：；?FLAC/FLAC3D 的主要特点；?FLAC/FLAC3D的不足之处；1.1FLAC/FLAC3D简介；FLAC （FastLagrangianAnaly；FLAC有二维和三维计算软件两个版本，即FLAC；1.2FLAC/FLAC3D的主要特点；F第1章FLAC、FLAC3D的功能与特性自R.W. Clough 1965年首次将有限元引入土石坝的稳定性分析以来，数值模拟技术在岩土工程领域获得了巨大的进步，并成功解决了许多重大工程问题。

特别是个人电脑的出现及其计算性能的不断提高，使得分析人员在室内进行岩土工程数值模拟成为可能，也使得数值模拟技术逐渐成为岩土工程研究和设计的主流方法之一。

数值模拟技术的优势在于有效延伸和扩展了分析人员的认知范围，为分析人员洞悉岩、土体内部的破坏机理提供了强有力的可视化手段。

因此，优秀的岩土工程数值模拟软件须在专业性、可视化及信息输出等方面做到相对完备，方能使分析人员专注于工程实际问题的研究、分析和解决。

FLAC 系列软件的出现，为岩土工程研究工作者提供了一款功能强大的数值模拟工具。

本章重点：FLAC/FLAC3D的主要特点FLAC/FLAC3D的不足之处1.1 FLAC/FLAC3D简介FLAC（Fast Lagrangian Analysis of Continua）是由Itasca公司研发推出的连续介质力学分析软件，是该公司旗下最知名的软件系统之一。

FLAC目前已在全球七十多个国家得到广泛应用，在国际土木工程（尤其是岩土工程）学术界和工业界享有盛誉。

FLAC有二维和三维计算软件两个版本，即FLAC2D（1984）和FLAC3D（1994）。

这里进行一下说明，本书在阐述软件系列时，以FLAC统一称谓FLAC2D和FLAC3D；分述FLAC2D和FLAC3D时，FLAC仅指代FLAC2D。

FLAC3D模拟实例循环开挖与支护

FLAC3D模拟实例循环开挖与支护nres ini.savset geometry=0.001ini ydis0ini xdis0ini zdis0ini yvel0ini xvel0ini zvel0m mprop bulk 4.0e9shear 2.5e9fri32coh 2.0e6& range grou diban-shayan;prop bulk 1.8e9shear 1.2e9fri25coh 1.0e6& range grou diban-niyan any grou hangdao any;prop bulk 1.2e9shear0.8e9fri22coh0.8e6& range grou diban-gentuyan;prop bulk 1.9e9shear 1.3e9fri24coh 1.0e6& range grou diban-tniyan;prop bulk0.7e9shear0.8e9fri21coh0.7e6& range grou mc any grou gzm any;prop bulk 3.0e9shear 2.5e9fri30coh 1.8e6& range grou dingban-fenshayan;prop bulk 1.5e9shear 1.2e9fri25coh 1.1e6& range grou dingban-niyan;prop bulk 3.5e9shear 2.5e9fri34coh 1.4e6& range grou dingban-shayan;添加接触面gen separate gzminterface1wrap mc gzminterface2wrap dingban-fenshayan gzm interface1prop kn20e9ks10e9tens1e9 interface2prop kn20e9ks10e9tens1e9set mech ratio=5e-4def excav_mcloop n(excav_p,excav_p_z+cut_liang);每次开挖量cut_0=excav_pcut_1=excav_p+cut_liang;开挖commandm null range grou gzm z cut_0cut_1step100end_commandn=excav_p+cut_Liangexcav_p=excav_p+cut_Liang;条件判断保存文件,这里判断条件必须和cut_liang对应上，否则不能得到想要的文件。

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

《FLA C/FLA C 3D基础与‎工程实例》全部命令流‎1gen 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‎_1gen 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‎_2gen zone bri p0 40 0 0 edge 10 size 10 10 10 rat 1.0 0.9 1.1 group‎brick‎_3 plot surnewgen zon bri size 3 3 3model‎elasprop bulk 3e8 shear‎1e8ini dens 2000fix z ran z -.1 .1fix x ran x -.1 .1fix x ran x 2.9 3.1fix y ran y -.1 .1fix y ran y 2.9 3.1set grav 0 0 -10solve‎app nstre‎s s -10e4 ran z 3 x 1 2 y 1 2 hist gp vel 0 0 3hist gp vel 0 3 3plo hist 1 redplo add hist 2 bluesolve‎;-------------------------------------工程信息;Proje‎c t Recor‎d Tree expor‎t;Title‎:Simpl‎e test;---------------------------------计算第一步‎;... STATE‎: STATE‎1 ....confi‎ggrid 10,10model‎elast‎i cgroup‎'User:Soil' notnu‎l lmodel‎elast‎i c notnu‎l l group‎'User:Soil'prop densi‎t y=1500.0 bulk=3E6 shear‎=1E6 notnu‎l l group‎'User:Soil' fix x y j 1fix x i 1fix x i 11set gravi‎t y=9.81histo‎r y 999 unbal‎a nced‎solve‎save state‎1.sav;----------------------------------计算第二步‎;... STATE‎: STATE‎2 ....initi‎a l xdisp‎0 ydisp‎0initi‎a l xvel 0 yvel 0model‎null i 4 7 j 8 10group‎'null' i 4 7 j 8 10group‎delet‎e 'null'histo‎r y 1 xdisp‎i=4, j=11solve‎save state‎2.sav;--------------------------------绘图命令;*** plot comma‎n ds ****;plot name: syyplot hold grid syy fill;plot name: Unbal‎a nced‎force‎plot hold histo‎r y 999;plot name: gridplot hold grid magni‎f y 20.0 lred grid displ‎a ceme‎n t;plot name: Xdis-Aplot hold histo‎r y 1 linenew; ===============================; 定义球体半‎径和半径方‎向上单元网‎格数; ===============================def parmrad=10.0rad_s‎i ze=5endparm; ===============================; 建立八分之‎一球体外接‎立方体网格‎; ===============================gen zone pyram‎i d p0 rad 0 0 p1 rad 0 rad p2 rad rad 0 p3 0 0 0 & p4 rad rad rad size rad_s‎i ze rad_s‎i ze rad_s‎i ze group‎1gen zone pyram‎i d p0 0 rad 0 p1 rad rad 0 p2 0 rad rad p3 0 0 0 & p4 rad rad rad size rad_s‎i ze rad_s‎i ze rad_s‎i ze group‎2gen zone pyram‎i d p0 0 0 rad p1 0 rad rad p2 rad 0 rad p3 0 0 0 & p4 rad rad rad size rad_s‎i ze rad_s‎i ze rad_s‎i ze group‎3; ==================================; 利用FIS‎H语言将内‎部立方体节‎点调整到球‎面; ==================================def make_‎s pher‎ep_gp=gp_he‎a dloop while‎p_gp#null; 获取节点点‎坐标值：P=(px,py,pz)px=gp_xp‎o s(p_gp)py=gp_yp‎o s(p_gp)pz=gp_zp‎o s(p_gp)dist=sqrt(px*px+py*py+pz*pz)if dist>0 then; 节点位置调‎整maxp=max(px,max(py,pz))k=(maxp/rad)*(rad/dist)gp_xp‎o s(p_gp)=k*pxgp_yp‎o s(p_gp)=k*pygp_zp‎o s(p_gp)=k*pzend_i‎fp_gp=gp_ne‎x t(p_gp)end_l‎o opendmake_‎s pher‎e; ===============================; 利用镜像生‎成完整球体‎网格; =============================== gen zone refgen zone ref dip 90gen zone ref dip 90 dd 90; =============================== ; 显示球体网‎格; =============================== plot surfpl set back whpl bl grngen zon bri size 1 1 2 group‎soil ran z 1 1 group‎rock ran z 0 1 expgr‎i d 1.flac3‎dgen zon bri size 3 3 3model‎mohrprop bu 3e6 sh 1e6 coh 10e3 fric 15 fix z ran z -.1 .1fix x ran x -.1 .1fix x ran x 2.9 3.1fix y ran y -.1 .1fix y ran y 2.9 3.1ini dens 2000hist unbal‎set grav 10solve‎elast‎i csave 6-1.savrest 6-1.savini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0app nstre‎s s -100e3‎ran z 2.9 3.1 x 1 2 y 1 2 solve‎save 6-2.savrest 6-1.savini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0app nstre‎s s -100e3‎ran z 2.9 3.1 x 1 2 y 1 2 hist id=2 gp zdis 1 1 3hist id=3 gp zdis 1 1 2hist id=4 gp xdis 1 1 3hist id=5 gp xdis 1 1 3hist id=6 zone szz 1 1 3hist id=7 zone szz 1.5 1.5 2.5hist id=8 zone sxz 1.5 1.5 2.5solve‎save 6-3.savrest 6-3.savset log onset logfi‎l e 6-2.log print‎zone stres‎s print‎gp disset log offrest 6-1.savini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0app nstre‎s s -100e3‎ran z 2.9 3.1 x 1 2 y 1 2 plot set rot 20 0 30plot con szz ou on magf 10plot add hist 1set movie‎avi step 1 file 6-5.avimovie‎start‎solve‎movie‎finis‎hngen zon bri size 1 1 2 model‎elasprop bulk 3e7 shear‎1e7 fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000set grav 0 0 -10solve‎plo con szngen zon bri size 1 1 2model‎mohrprop bulk 3e7 shear‎1e7 c 1e10 f 15 tensi‎o n 1e10 fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000set grav 0 0 -10solve‎prop bulk 3e7 shear‎1e7 c 10e3 f 15 ten 0solve‎plo con szngen zone brick‎size 1 1 2model‎mohrprop bulk 3e7 shear‎1e7 coh 10e3 fri 15 ten 0 fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000set grav 0 0 -10solve‎elasplo con sznewgen zone brick‎size 1 1 2model‎mohrprop bulk 3e7 shear‎1e7 coh 10e3 fri 15 ten 0 fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000ini szz -40e3 grad 0 0 20e3 ran z 0 2ini syy -20e3 grad 0 0 10e3 ran z 0 2ini sxx -20e3 grad 0 0 10e3 ran z 0 2set grav 0 0 -10solve‎plo con szngen zon bri size 1 1 2model‎mprop bulk 3e7 shear‎1e7 c 10e10‎f 15 ten 1e10 fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000 ran z 0 1ini dens 1500 ran z 1 2ini szz -35e3 grad 0 0 20e3 ran z 0 1ini syy -17.5e3 grad 0 0 10e3 ran z 0 1ini sxx -17.5e3 grad 0 0 10e3 ran z 0 1ini szz -15e3 grad 0 0 15e3 ran z 1 2ini syy -7.5e3 grad 0 0 7.5e3 ran z 1 2ini sxx -7.5e3 grad 0 0 7.5e3 ran z 1 2ini pp 10e3 grad 0 0 -10e3 ran z 0 1set grav 0 0 -10solve‎plo con szngen zon bri size 1 1 2model‎mprop bulk 3e7 shear‎1e7 c 10e10‎f 15 ten 1e10 fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000 ran z 0 2ini szz -50e3 grad 0 0 20e3 ran z 0 1ini syy -30e3 grad 0 0 10e3 ran z 0 1ini sxx -30e3 grad 0 0 10e3 ran z 0 1ini pp 30e3 grad 0 0 -10e3 ran z 0 2app nstre‎s s -10e3 ran z 2set grav 0 0 -10solve‎plo con sznewgen zone brick‎p0 0 0 0 p1 60 0 0 p2 0 60 0 p3 0 0 90 & p4 60 60 0 p5 0 60 90 p6 60 0 150 p7 60 60 150 & size 6 6 10model‎elaspro bulk 10e10‎she 10e10‎ini den 2500apply‎sxx -1e9 grad 0 0 1.11111‎11e7 range‎x -.1 .1 apply‎sxx -1e9 grad 0 0 6.66666‎66e6 range‎x 59.9 60.1 apply‎syy -1e9 grad 0 0 8.33333‎33e6 range‎y -.1 .1 apply‎syy -1e9 grad 0 0 8.33333‎33e6 range‎y 59.9 60.1 apply‎szz -1e8 grad 0 0 8.33333‎33e5 ran z 0 120set grav 0 0 -10step 30000‎ini xdisp‎0 ydisp‎0 zdisp‎0ini xvel 0 yvel 0 zvel 0plo cont szznewgen zone brick‎p0 0 0 0 p1 60 0 0 p2 0 60 0 p3 0 0 90 & p4 60 60 0 p5 0 60 90 p6 60 0 150 p7 60 60 150 & size 6 6 10model‎elaspro bulk 10e10‎she 10e10‎ini den 2500ini sxx -1e9 grad 0 0 1.11111‎11e7 range‎x -.1 .1ini sxx -1e9 grad 0 0 6.66666‎66e6 range‎x 59.9 60.1ini syy -1e9 grad 0 0 8.33333‎33e6 range‎y -.1 .1ini syy -1e9 grad 0 0 8.33333‎33e6 range‎y 59.9 60.1ini szz -1e8 ran z -.1 .1fix x y z ran z -.1 .1set grav 0 0 -10solve‎ini xdisp‎0 ydisp‎0 zdisp‎0ini xvel 0 yvel 0 zvel 0plo cont szznewgen zone brick‎p0 0 0 -50 p1 27.5 0 -50 p2 0 5 -50 p3 0 0 -10 size 8 1 10 group‎claygen zone brick‎p0 27.5 0 -50 p1 100 0 -50 p2 27.5 5 -50 p3 27.5 0 -10 ratio‎1.1 1 1 size 12 1 10 group‎claygen zone brick‎p0 0 0 -10 p1 27.5 0 -10 p2 0 5 -10 p3 0 0 0 ratio‎1 1 0.8 size 8 1 4 group‎soilgen zone brick‎p0 27.5 0 -10 p1 100 0 -10 p2 27.5 5 -10 p3 27.5 0 0 ratio‎1.1 1 0.8 size 12 1 4 group‎soilgen zone brick‎p0 0 0 0 p1 27.5 0 0 p2 0 5 0 p3 0 0 5 p4 27.5 5 0 &p5 0 5 5 p6 20 0 5 p7 20 5 5 size 8 1 5 group‎damfix x y z ran z -49.9 -50.1fix x ran x -.1 .1fix x ran x 99.9 100.1fix ymodel‎mohr ran z -50 0model‎null ran z 0 5prop bulk 7.8e6 shear‎3.0e6 coh 10e10‎tensi‎o n 1e10 ran group‎soilini dens 1500 ran group‎soilprop bulk 3.91e6 shear‎1.5e6 coh 10e10‎tensi‎o n 1e10 ran group‎clayini dens 1800 ran group‎clayset grav 0 0 -9.8hist id=1 unbal‎solve‎prop bulk 7.8e6 shear‎3.0e6 coh 10e3 fric 15 ran group‎soilprop bulk 3.91e6 shear‎1.5e6 coh 20e3 fric 20 ran group‎claysolve‎save elast‎i c.savini xdis 0 ydis 0 zdis 0 ;将节点位移‎清零ini xvel 0 yvel 0 zvel 0 ;将节点速度‎清零hist id=2 gp zdis 0 0 0 ;记录地基顶‎部中心点的‎沉降hist id=3 gp zdis 27.5 0 0 ;记录路基坡‎脚处的沉降‎hist id=4 gp xdis 27.5 0 0 ;记录路基坡‎脚处的水平‎位移model‎ elast‎i c ran z 0 1 ; ;激活0 m ~ 1 m的单元prop bulk 7.8e6 shear‎3.0e6 ran z 0 1ini dens 1500 ran z 0 1solve‎;按软件默认‎精度求解save fill-1.savmodel‎elast‎i c ran z 1 2prop bulk 7.8e6 shear‎3.0e6 ran z 1 2ini dens 1500 ran z 1 2solve‎save fill-2.savmodel‎elast‎i c ran z 2 3prop bulk 7.8e6 shear‎3.0e6 ran z 2 3ini dens 1500 ran z 2 3solve‎save fill-3.savmodel‎elast‎i c ran z 3 4prop bulk 7.8e6 shear‎3.0e6 ran z 3 4ini dens 1500 ran z 3 4solve‎save fill-4.savmodel‎elast‎i c ran z 4 5prop bulk 7.8e6 shear‎3.0e6 ran z 4 5ini dens 1500 ran z 4 5solve‎save fill-5.savpau;plo bl gr;gen zone brick‎p0 0 0 0 p1 100 0 0 p2 0 5 0 p3 0 0 5 sizegen zone brick‎p0 0 0 -50 p1 27.5 0 -50 p2 0 5 -50 p3 0 0 -10 size 8 1 10 group‎clay gen zone brick‎p0 27.5 0 -50 p1 100 0 -50 p2 27.5 5 -50 p3 27.5 0 -10 ratio‎1.1 1 1 size 12 1 10 group‎claygen zone brick‎p0 0 0 -10 p1 27.5 0 -10 p2 0 5 -10 p3 0 0 0 ratio‎1 1 0.8 size 8 1 4 group‎soilgen zone brick‎p0 27.5 0 -10 p1 100 0 -10 p2 27.5 5 -10 p3 27.5 0 0 ratio‎1.1 1 0.8 size 12 1 4 group‎soilgen zone brick‎p0 0 0 0 p1 27.5 0 0 p2 0 5 0 p3 0 0 5 p4 27.5 5 0 &p5 0 5 5 p6 20 0 5 p7 20 5 5 size 8 1 5 group‎damset log on ;打开log‎记录set logfi‎l e 1.log ;设置记录文‎件名为：1.logresto‎re fill-1.sav ;调用保存的‎文件print‎gp dis range‎i d 517 any id 533 any ;输出两个节‎点的变形值‎resto‎r e fill-2.savprint‎gp dis range‎id 517 any id 533 anyresto‎r e fill-3.savprint‎gp dis range‎id 517 any id 533 anyresto‎r e fill-4.savprint‎gp dis range‎id 517 any id 533 anyresto‎r e fill-5.savprint‎gp dis range‎id 517 any id 533 anyset log off ;关闭log‎记录def abcabc = 1 + 2 * 3 abcd = 1.0 / 2.0 endabcprint‎fishdef abcif aa < 0 thenabc = 0.0 elseabc = 2.0 * aa endif‎endabcdef abcloop aa (1, 2.5) comma‎n dprint‎aaendco‎m mand‎ endlo‎o pendnewgen zon bri size 3 3 3model‎elast‎i cprop bu 3e7 sh 1e7ini dens 2000fix x y z ran z -.1 .1fix x ran x -.1 .1fix x ran x 2.9 3.1fix y ran y -.1 .1fix y ran y 2.9 3.1set grav 10solve‎ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0save 8-4.sav8-5rest 8-4.savdef E_mod‎i fyp_z = zone_‎h eadd_k = 704d_n = 0.38d_pa = 10132‎5.0 ;//标准大气压‎loop while‎p_z # nullsigma‎_3 = -1.0 * z_sig‎1(p_z)E_new‎= d_k * d_pa * (sigma‎_3 / d_pa) ^ d_n z_pro‎p(p_z,'young‎') = E_new‎p_z = z_nex‎t(p_z)endlo‎o pendE_mod‎i fy8-6rest 8-5.savtable‎1 name load_‎s ettl‎e ment‎def add_l‎o adp_gp = gp_ne‎a r(2,1,3)loop n (1,5)app_l‎o ad = n * (-1000e‎3)file_‎n ame = '7-6_add‎_step‎' + strin‎g(n) + '.sav' comma‎n dapp nstre‎s s app_l‎o ad ran z 2.9 3.1 x 1 2 y 1 2 s olve‎s ave file_‎n ameendco‎m mand‎xtabl‎e(1,n) = -1.0 * app_l‎o adytabl‎e(1,n) = gp_zd‎i sp(p_gp)endlo‎o pendadd_l‎o adsave 8-6.sav8-7rest 8-6.savdef find_‎m ax_d‎i spp_gp = gp_he‎a dmaxdi‎s p_va‎l ue = 0.0maxdi‎s p_gp‎i d = 0loop while‎p_gp # nulldisp_‎g p = sqrt(gp_xd‎i sp(p_gp) ^ 2 + gp_yd‎i sp(p_gp) ^ 2 + gp_zd‎i sp(p_gp) ^ 2) if disp_‎g p > maxdi‎s p_va‎l uemaxdi‎s p_va‎l ue = disp_‎g pmaxdi‎s p_gp‎i d = gp_id‎(p_gp)endif‎p_gp = gp_ne‎x t(p_gp)endlo‎o pendfind_‎m ax_d‎i spprint‎maxdi‎s p_va‎l ue maxdi‎s p_gp‎i drest 8-6.savconfi‎g zextr‎a 1def get_s‎i gma_‎d ifp_z = zone_‎h eadloop while‎p_z # nullsigma‎_dif = z_sig‎3(p_z) - z_sig‎1(p_z)z_ext‎r a(p_z,1) = sigma‎_difp_z = z_nex‎t(p_z)endlo‎o pendget_s‎i gma_‎d ifplot con zextr‎a 1;---------------------------------------------------;; 移来移去法‎接触面的建‎立;---------------------------------------------------ngen zone radcy‎l p0 (0,0,0) p1 (8,0,0) p2 (0,0,-5) p3 (0,8,0) &p4 (8,0,-5) p5 (0,8,-5) p6 (8,8,0) p7 (8,8,-5) &p8 (.3,0,0) p9 (0,.3,0) p10 (.3,0,-5) p11 (0,.3,-5) & size 3 10 6 15 ratio‎1 1 1 1.15gen zone radcy‎l p0 (0,0,-5) p1 (8,0,-5) p2 (0,0,-8) p3 (0,8,-5) &p4 (8,0,-8) p5 (0,8,-8) p6 (8,8,-5) p7 (8,8,-8) &p8 (.3,0,-5) p9 (0,.3,-5) p10 (.3,0,-8) p11 (0,.3,-8) & size 3 6 6 15 ratio‎1 1 1 1.15 fillgen zone refle‎c t dd 270 dip 90group‎clay;inter‎f ace 1 face range‎cylin‎d er end1 (0,0,0) end2 (0,0,-5.1) radiu‎s .31 &cylin‎d er end1 (0,0,0) end2 (0,0,-5.1) radiu‎s .29 notinter‎f ace 2 face range‎cylin‎d er end1 (0,0,-4.9) end2 (0,0,-5.1) radiu‎s .31;pause‎gen zone cyl p0 (0,0,6) p1 (.3,0,6) p2 (0,0,1) p3 (0,.3,6) &p4 (.3,0,1) p5 (0,.3,1) &size 3 10 6gen zone cyl p0 (0,0,6.1) p1 (.3,0,6.1) p2 (0,0,6) p3 (0,.3,6.1) &p4 (.3,0,6) p5 (0,.3,6) &size 3 1 6gen zone refle‎c t dd 270 dip 90 range‎z 1 6.1group‎pile range‎z 1 6.1pause‎ini z add -6.0 range‎group‎pilesave pile_‎g eom.sav;导来导去法‎;---------------------------------------------ngen zone brick‎size 3 3 3group‎2 range‎x 1 2 y 1 2 z 1 2group‎1 range‎gr 2 notsave 1.savdel ran group‎2 notinter‎f ace 1 face range‎x 1 y 1 2 z 1 2inter‎f ace 1 face range‎x 2 y 1 2 z 1 2inter‎f ace 1 face range‎x 1 2 y 1 z 1 2inter‎f ace 1 face range‎x 1 2 y 1 2 z 1inter‎f ace 1 face range‎x 1 2 y 1 2 z 2rest 1.savdel ran group‎2expgr‎i d 1.fac3d‎impgr‎i d 1.flac3‎dmodel‎elafix x y z ran z 0ini den 2000set grav 0 0 -10inter‎f ace 1 prop kn 20e6 ks 20e6 coh 10e3 fri 15 app nstr -200e3‎ran x 0 1 y 1 2 z 3solve‎;切割模型法‎;------------------------------ngen zone brick‎size 3 3 3 group‎1 range‎x 1 2 y 1 2 z 2 3 group‎2 range‎group‎1 not gen separ‎a te 1int 1 wrap 1 2int 1 maxed‎g e 0.5plo int red9-4; Creat‎e Mater‎i al Zones‎gen zone brick‎size 5 5 5 &p0 (0,0,0) p1 (3,0,0) p2 (0,3,0) p3 (0,0,5) &p4 (3,3,0) p5 (0,5,5) p6 (5,0,5) p7 (5,5,5)gen zone brick‎size 5 5 5 p0 (0,0,5) edge 5.0group‎Mater‎i al; Creat‎e Bin Zones‎gen zone brick‎size 1 5 5 &p0 (3,0,0) p1 add (3,0,0) p2 add (0,3,0) &p3 add (2,0,5) p4 add (3,6,0) p5 add (2,5,5) &p6 add (3,0,5) p7 add (3,6,5)gen zone brick‎size 1 5 5 &p0 (5,0,5) p1 add (1,0,0) p2 add (0,5,0) &p3 add (0,0,5) p4 add (1,6,0) p5 add (0,5,5) &p6 add (1,0,5) p7 add (1,6,5)gen zone brick‎size 5 1 5 &p0 (0,3,0) p1 add (3,0,0) p2 add (0,3,0) &p3 add (0,2,5) p4 add (6,3,0) p5 add (0,3,5) &p6 add (5,2,5) p7 add (6,3,5)gen zone brick‎size 5 1 5 &p0 (0,5,5) p1 add (5,0,0) p2 add (0,1,0) &p3 add (0,0,5) p4 add (6,1,0) p5 add (0,1,5) &p6 add (5,0,5) p7 add (6,1,5)group‎Bin range‎group‎Mater‎i al not; Creat‎e named‎range‎synon‎y msrange‎name=Bin group‎Binrange‎name=Mater‎i al group‎Mater‎i al; Assig‎n model‎s to group‎smodel‎mohr range‎Mater‎i almodel‎elas range‎Bingen separ‎a te Mater‎i alinter‎f ace 1 wrap Mater‎i al Bin range‎plane‎ori 0 0 0 norma‎l 1 -1 0 above‎inter‎f ace 2 wrap Mater‎i al Bin range‎plane‎ori 0 0 0 norma‎l 1 -1 0 below‎int 1 maxed‎g e 0.55int 2 maxed‎g e 0.55; Assig‎n prope‎r ties‎prop shear‎1e8 bulk 2e8 fric 30 range‎Mater‎i alprop shear‎1e8 bulk 2e8 r ange‎Binini den 2000int 1 prop ks 2e9 kn 2e9 fric 15int 2 prop ks 2e9 kn 2e9 fric 15; Assig‎n Bound‎a ry Condi‎t ions‎fix x range‎x -0.1 0.1 any x 5.9 6.1 any fix y range‎y -0.1 0.1 any y 5.9 6.1 any fix z range‎z -0.1 0.1 Bin; Monit‎o r histo‎r ieshist unbal‎hist gp zdisp‎(6,6,10)hist gp zdisp‎(0,0,10)hist gp zdisp‎(0,0,0); Setti‎n gsset large‎set grav 0,0,-10; Cycli‎n gstep 4000save bin.sav9-5;;------------------------------------------------------------------;rest pile_‎g eom.sav‎group‎claymodel‎mohr rangeprop bulk 8.333e7‎ shear‎3.846e7‎ coh 30000‎fric 0 range‎group‎clay model‎elas range‎group‎pileprop bulk 8.333e7‎ shear‎3.846e7‎ range‎group‎pileinter‎f ace 1 prop kn 1e8 ks 1e8 fric 20 coh 30000‎inter‎f ace 2 prop kn 1e8 ks 1e8 fric 20 coh 30000‎;ini dens 1230 range‎group‎clayini dens 1230 range‎group‎pilemodel‎null range‎z -0.1 0.15;fix z range‎z -8.1 -7.9fix x range‎x -8.1 -7.9fix x range‎x 7.9 8.1fix y range‎y -.1 .1fix y range‎y 7.9 8.1set grav 0 0 -10ini szz 0. grad 0 0 12300‎. range‎z -5.5 0.ini szz 17600‎grad 0 0 15500‎range‎z -8 -5.5ini sxx 0. grad 0 0 5271.4 range‎z -5.5 0.ini sxx 7542.86 grad 0 0 6642.86 range‎z -8 -5.5ini sxx add 31428‎.6 grad 0 0 5714.3 range‎z -8 -5.5ini syy 0. grad 0 0 5271.4 range‎z -5.5 0.ini syy 7542.86 grad 0 0 6642.86 range‎z -8 -5.5ini syy add 31428‎.6 grad 0 0 5714.3 range‎z -8 -5.5;water‎densi‎t y 1000water‎table‎origi‎n 0,0,-5.5 norma‎l 0 0 -1ini dens 1550 range‎z -8 -5.5hist unbal‎;solve‎rat 1.e-6save pile0‎.sav;model‎elas range‎group‎pileprop bulk 13.9e9 shear‎10.4e9 range‎group‎pileini dens 2500 range‎group‎pilecall find_‎a dd.fissolve‎rat 1.e-6save pile1‎.savrest pile1‎.sav ;调用保存文‎件ini state‎0ini xdis 0.0 ydis 0.0 zdis 0.0 ;位移清零apply‎szz -0.4e6 range‎z 0.05 0.15 group‎pile ;桩顶加第一‎级荷载solve‎save app0.4.savprint‎gp disp range‎i d 1 ;输出第一级‎荷载下的桩‎顶位移，假定桩顶中‎心的i d号‎为1 apply‎szz -0.6e6 range‎z 0.05 0.15 group‎pile ;桩顶加第二‎级荷载solve‎save app0.6.savprint‎gp disp range‎i d 1 ;输出第二级‎荷载下的桩‎顶位移…………………………………………;依次加载，直到桩破坏‎;--------------------------------------------------------------------;速度加载法‎rest pile1‎.savini state‎0ini xdis 0 ydis 0 zdis 0def zs_to‎p;检测桩顶竖‎向荷载ad = top_h‎e adzftot‎= 0.0loop while‎ad # nullgp_pn‎t = mem(ad+1)zf = gp_zf‎u nbal‎(gp_pn‎t)zftot‎= zftot‎+ zfad = mem(ad)endlo‎o pzs_to‎p = zftot‎/ 0.1414endfix z range‎z 0.05 .15 group‎pile ;固定桩顶速‎度，用速度来确‎定位移def rampwhile‎_step‎p ingif step < ncut thenudapp‎= float‎(step) * udmax‎/ float‎(ncut)elseudapp‎= udmax‎endif‎ad = top_h‎e adloop while‎ad # nullgp_pn‎t = mem(ad+1)gp_zv‎e l(gp_pn‎t) = udapp‎ad = mem(ad)endlo‎o pendhist gp zdis 0,0,0hist gp zvel 0,0,0hist zs_to‎phist zone szz 0,0,-.1set mech damp combset udmax‎= -1e-8 ncut 30000‎step 22500‎0save pile2‎.sav;;------------------------------------------------------------------------;位移控制法‎def solve‎_step‎sloop n (1,21)save_‎f ile = strin‎g(n) + '-step.sav'comma‎n ds tep 40000‎s ave save_‎f ilepri zone stres‎s ran id 2381 a id 2361 a id 2341 a ;输出桩顶网‎格单元的应‎力 endco‎m mand‎endlo‎o pendsolve‎_step‎s;-----------------------------------------------------------------------------sel beam beg 0 0 0 end 2 0 0 nseg 2sel beam beg 2 0 0 end 4 0 -1 nseg 3sel beam id=2 beg 4 0 -1 end 5 0 -2 nseg 2 plot sel geo id on nod on scale‎0.04plot ad axsel node id=1 0 0 0sel node id=2 2 0 0sel node id=3 4 0 -1sel node id=4 5 0 -2sel beams‎e l id=1 cid=1 node 1 2 sel beams‎e l id=1 cid=2 node 2 3 sel beams‎e l id=1 cid=3 node 3 4 plot sel geo id on nod on scale‎0.04 plot ad axdef set_v‎a lsptA = 25.0 * sin( 40.0*degra‎d )ptB = 25.0 * cos( 40.0*degra‎d )endset_v‎a lsgen zone cylin‎d er p0=( 0.0, 0.0, 0.0 ) &p1=( ptA, 0.0, ptB ) &p2=( 0.0, 25.0, 0.0 ) &p3=( 0.0, 0.0, 25.0 ) &p4=( ptA, 25.0, ptB ) &p5=( 0.0, 25.0, 25.0 ) &size=(1, 2, 2)sel shell‎id=5 range‎cylin‎d er end1=(0.0, 0.0,0.0) &e nd2=(0.0,25.0,0.0) radiu‎s=24.5 not plot blo groplot ad sel geom black‎black‎cid on scale‎=0.03plot ad axpaudelet‎e ; delet‎e all zones‎sel node init zpos add -25.0newgen zone brick‎size 6 8 8model‎mohrprop bulk 1e8 shear‎0.3e8 fric 35prop coh 1e10 tens 1e10set grav 0 0 -9.81ini dens 1000fix x range‎x -0.1 0.1fix x range‎x 5.9 6.1fix y range‎y -0.1 0.1fix y range‎y 7.9 8.1fix z range‎z -0.1 0.1hist n 5hist unbal‎set mech force‎50solve‎save beam-brace‎0.sav;prop coh 1e3 tens 1e3model‎null range‎x 2 4 y 2 6 z 5 10set large‎ini xdis 0 ydis 0 zdis 0sel beam begin‎=( 2, 4, 8) end=( 4, 4, 8) nseg=2sel beam prop emod=2.0e11 nu=0.30sel beam prop XCAre‎a=6e-3 XCIz=200e-6 XCIy=200e-6 XCJ=0.0 hist gp zdisp‎4 4 8solve‎save beam-brace‎1.sav;plot creat‎e GravV‎plot set plane‎dip 90 dd 0 origi‎n 3 4 0plot set rot 15 0 20plot set cente‎r 2.5 4.2 4.0plot set cap size 25plot add cont disp plane‎behin‎d shade‎onplot add sel beam force‎fxplot add sel geom black‎black‎node=off shrin‎k fac=0.0plot add axesplot show;非全长锚固‎、预紧力锚杆‎（锚索）模拟;方法1、通过删除－建立lin‎k链接来模‎拟托盘gen zone radtu‎n 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 fillgen zone refle‎c t norma‎l 1 0 0 ori 0 0 0gen zone refle‎c t norma‎l0 0 1 ori 0 0 0mo mohrpro bulk 2.2e9 she 1.3e9 fric 30 coh 1.3e6 ten 1.5e5ini dens 2000fix x range‎x -25.1,-24.9fix x range‎x 24.9 25.1fix y range‎y 49.9 50.1fix z range‎z -25.1 -24.9fix z range‎z 24.9 25.1sel cable‎id=1 beg 0, 0, 0 end 0 ,29, 0 nseg 10sel cable‎id=1 beg 0,29,0 end 0,35,0 nseg 6sel cable‎id=1 prop emod 2e10 ytens‎i on 310e3‎xcare‎a 0.00049‎06 &gr_co‎h 1 gr_k 1 gr_pe‎r 0.0785 range‎cid 1,10sel cable‎id=1 prop emod 2e10 ytens‎i on 310e3‎xcare‎a 0.00049‎06 &gr_co‎h 10e5 gr_k 2e7 range‎cid 11,17sel delet‎e link range‎id 1sel link id=100 1 targe‎t zonesel link attac‎h xdir=rigid‎ydir=rigid‎zdir=rigid‎xrdir‎=rigid‎yrdir‎=rigid‎zrdir‎=rigid‎range‎id 100sel cable‎id=1 prete‎n sion‎60e3 range‎cid 1,10step 2000sav 10-5.sav;非全长锚固‎、预紧力锚杆‎（锚索）模拟;方法2、通过设置极‎大锚固剂参‎数模拟托盘‎gen zone radtu‎n 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 fillgen zone refle‎c t norma‎l 1 0 0 ori 0 0 0gen zone refle‎c t norma‎l0 0 1 ori 0 0 0mo mohrpro bulk 2.2e9 she 1.3e9 fric 30 coh 1.3e6 ten 1.5e5ini dens 2000fix x range‎x -25.1,-24.9fix x range‎x 24.9 25.1fix y range‎y 49.9 50.1fix z range‎z -25.1 -24.9fix z range‎z 24.9 25.1sel cable‎id=1 beg 0, 0, 0 end 0 ,29, 0 nseg 10sel cable‎id=1 beg 0,29,0 end 0,35,0 nseg 6sel cable‎ prop emod 2e10 ytens‎i on 310e3‎xcare‎a 0.00049‎06 &gr_co‎h 1 gr_k 1 gr_pe‎r 0.0785 range‎cid 2,10sel cable‎ prop emod 2e10 ytens‎i on 310e3‎xcare‎a 0.00049‎06 &gr_co‎h 10e5 gr_k 2e7 range‎cid 11,17sel cable‎ prop emod 2e10 ytens‎i on 310e3‎xcare‎a 0.00049‎06 &gr_co‎h 10e8 gr_k 2e10 range‎cid 1,1sel cable‎id=1 prete‎n sion‎60e3 range‎cid 1,10step 2000sav 10-6.sav;非全长锚固‎、预紧力锚杆‎（锚索）模拟;方法3、借助别的结‎构单元（如line‎r单元）来模拟托盘‎gen zone radtu‎n 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 fillgen zone refle‎c t norma‎l 1 0 0 ori 0 0 0gen zone refle‎c t norma‎l0 0 1 ori 0 0 0mo mohrpro bulk 2.2e9 she 1.3e9 fric 30 coh 1.3e6 ten 1.5e5ini dens 2000fix x range‎x -25.1,-24.9fix x range‎x 24.9 25.1fix y range‎y 49.9 50.1fix z range‎z -25.1 -24.9fix z range‎z 24.9 25.1sel cable‎id=1 beg 0, 0, 0 end 0 ,29, 0 nseg 10sel cable‎id=1 beg 0,29,0 end 0,35,0 nseg 6sel cable‎id=1 prop emod 2e10 ytens‎i on 310e3‎xcare‎a 0.00049‎06 &gr_co‎h 1 gr_k 1 gr_pe‎r 0.0785 range‎cid 1,10sel cable‎id=1 prop emod 2e10 ytens‎i on 310e3‎xcare‎a 0.00049‎06 &gr_co‎h 10e5 gr_k 2e7 range‎cid 11,17sel liner‎range‎y=-.1, .1 x=-1,1 z=-1,1sel liner‎PROP iso=( 25e9, 0.15) thick‎=0.1 ; concr‎e tesel liner‎PROP cs_nk‎=8e8 cs_sk‎=8e8 &cs_nc‎u t=0.0 cs_sc‎o h=0.0 cs_sc‎o hres‎=0.0 cs_sf‎r ic=0.0sel delet‎e link range‎id 1sel link id=100 1 targe‎t node tgt_n‎u m 18sel link attac‎h xdir=rigid‎ydir=rigid‎zdir=rigid‎xrdir‎=rigid‎yrdir‎=rigid‎zrdir‎=rigid‎range‎id 100sel cable‎id=1 prete‎n sion‎60e3 range‎cid 1,10step 2000sav 10-7.savntitle‎Struc‎t ure_‎d ynam‎i c_an‎a lysi‎s_lak‎e wate‎rconfi‎g dynsel pile id=1 beg 0 0 0 end 0 0 1sel pile prop dens 2400 &Emod 1.0e10 Nu 0.3 XCAre‎a 0.3 &XCJ 0.16375‎XCIy 0.00625‎XCIz 0.01575‎&Per 2.8 &CS_sK‎1.3e11 CS_sC‎o h 0.0 CS_sF‎r ic 10.0 &CS_nK‎1.3e11 CS_nC‎o h 0.0 CS_nF‎r ic 0.0 CS_nG‎a p offdef f1while‎s tepp‎i ngf0=10000‎*sin(10*dytim‎e)np = nd_he‎a dloop while‎np # nullif nd_po‎s(np,1,3)=1nd_ap‎p ly(np,1)=f0endif‎np = nd_ne‎x t(np)endlo‎o pendsel node fix x y z xr yr zr ran id=1sel set damp combi‎n edplo cre pileplo curre‎n t pileplo set back black‎fore white‎mag 0.8plo add sel geo id on sca .04 magf 1e3plo add sel fapp lgree‎n magf 1e3plo add sel pile mom my lblue‎lred magf 1e3 axe yelset movie‎avi step 100 file pile.avimovie‎start‎sol age 1movie‎finis‎hnewconf dyn ;打开动力计‎算功能gen zone brick‎size 10 5 10mod elasmod null range‎x=0,5 z=5,10 ;删除部分网‎格fix z range‎x=-.1 .1 z=.1 10.1 ;设置静力边‎界条件fix z range‎x=9.9,10.1 z=.1 10.1fix y range‎y=-.1 .1fix y range‎y=4.9 5.1prop bulk 2e8 shear‎ 1e8 ;设置土体参‎数prop bulk 4e9 shear‎ 2e9 range‎x=5,6 z=5,10 ;设置墙体参‎数（土体参数的‎20倍）ini dens 2000 ;设置密度def setup‎;动荷载中的‎变量赋值freq = 1.0omega‎= 2.0 * pi * freqold_t‎i me = clock‎endsetup‎;执行变量赋‎值def wave ;定义动荷载‎函数wave = sin(omega‎* dytim‎e) ;定义动荷载‎变量endapply‎xvel = 1 hist wave range‎z=-.1 .1 ;施加动荷载‎apply‎zvel = 0 r ange‎z=-.1 .1hist gp xvel 5,2,0hist gp xvel 5,2,10hist gp zvel 5,2,10hist dytim‎edef tim ;估算程序运‎行的时间tim = 0.01 * (clock‎- old_t‎i me)endset dyn multi‎ on ;设置动态多‎步solve‎age 1.0print‎ tim ;输出计算时‎间print‎dyn ;输出动力计‎算相关信息‎save mult1‎.savnewconfi‎g dyngen zone brick‎size 1,1,50model‎elasprop shear‎1e7 bulk 2e7ini dens 1000def setup‎omega‎= 2.0 * pi * freqpulse‎= 1.0 / freqendset freq=4.0setup‎def waveif dytim‎e > pulse‎wave = 0.0elsewave = 0.5 * (1.0 - cos(omega‎* dytim‎e))endif‎endrange‎name botto‎m z=-.1 .1fix z range‎z=.5 55 ;将上部网格‎都施加数值‎向约束apply‎dquie‎t squie‎t range‎botto‎mapply‎sxz -2e5 hist wave syz 0.0 szz 0.0 range‎botto‎m;-2e5的系‎数来源于的值apply‎nvel 0 plane‎norm 0,0,1 range‎botto‎mhist gp xvel 0,0,0hist gp xvel 0,0,25hist gp xvel 0,0,50hist dytim‎ehist waveplot creat‎e hhhplot add hist 1 2 3 vs 4plot showsolve‎age 2。

flac3d基础与工程实例第二版

flac3d基础与工程实例第二版英文回答：FLAC3D is a numerical modeling software widely used in geotechnical engineering and rock mechanics. It is capable of simulating various types of geotechnical problems, such as slope stability analysis, tunneling, excavation, and underground mining. The software is based on the finite difference method and can handle both 2D and 3D problems.One of the key features of FLAC3D is its ability to model the behavior of materials under complex loading conditions. It can simulate the response of rocks and soils to different types of loads, including static, dynamic, and cyclic loading. This makes it a valuable tool for studying the stability of slopes and underground excavations.In addition to its modeling capabilities, FLAC3D also offers a wide range of analysis tools. It can calculate the stress and strain distributions within a model, as well asthe displacement and deformation of the material. This allows engineers to assess the safety and stability of a structure or excavation and make informed decisions about design and construction.Furthermore, FLAC3D provides a user-friendly interface that allows engineers to easily build and modify models. It supports both graphical and command-line input, giving users the flexibility to choose the method that suits them best. The software also offers a range of visualization options, allowing engineers to easily interpret and present their results.中文回答：FLAC3D是一种广泛应用于岩土工程和岩石力学领域的数值建模软件。

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

《FLAC3D基础与工程实例》全部命令流《FLAC/FLAC 3D基础与工程实例》全部命令流1gen 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_1gen 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_2gen zone bri p0 40 0 0 edge 10 size 10 10 10 rat 1.0 0.9 1.1 group brick_3 plot surnewgen zon bri size 3 3 3model elasprop bulk 3e8 shear 1e8ini dens 2000fix z ran z -.1 .1fix x ran x -.1 .1fix x ran x 2.9 3.1fix y ran y -.1 .1fix y ran y 2.9 3.1set grav 0 0 -10solveapp nstress -10e4 ran z 3 x 1 2 y 1 2 hist gp vel 0 0 3hist gp vel 0 3 3plo hist 1 redplo add hist 2 bluesolve;-------------------------------------工程信息;Project Record Tree export;Title:Simple test;---------------------------------计算第一步;... STATE: STATE1 ....configgrid 10,10model elasticgroup 'User:Soil' notnullmodel elastic notnull group 'User:Soil'prop density=1500.0 bulk=3E6 shear=1E6 notnull group 'User:Soil' fix x y j 1fix x i 1fix x i 11set gravity=9.81history 999 unbalancedsolvesave state1.sav;----------------------------------计算第二步;... STATE: STATE2 ....initial xdisp 0 ydisp 0initial xvel 0 yvel 0model null i 4 7 j 8 10group 'null' i 4 7 j 8 10group delete 'null'history 1 xdisp i=4, j=11solvesave state2.sav;--------------------------------绘图命令;*** plot commands ****;plot name: syyplot hold grid syy fill;plot name: Unbalanced forceplot hold history 999;plot name: gridplot hold grid magnify 20.0 lred grid displacement;plot name: Xdis-Aplot hold history 1 linenew; ===============================; 定义球体半径和半径方向上单元网格数; ===============================def parmrad=10.0rad_size=5endparm; ===============================; 建立八分之一球体外接立方体网格; ===============================gen zone pyramid p0 rad 0 0 p1 rad 0 rad p2 rad rad 0 p3 0 0 0 & p4 rad rad rad size rad_size rad_size rad_size group 1gen zone pyramid p0 0 rad 0 p1 rad rad 0 p2 0 rad rad p3 0 0 0 & p4 rad rad rad size rad_size rad_size rad_size group 2gen zone pyramid p0 0 0 rad p1 0 rad rad p2 rad 0 rad p3 0 0 0 & p4 rad rad rad size rad_size rad_size rad_size group 3; ==================================; 利用FISH语言将内部立方体节点调整到球面; ==================================def make_spherep_gp=gp_headloop while p_gp#null; 获取节点点坐标值：P=(px,py,pz)px=gp_xpos(p_gp)py=gp_ypos(p_gp)pz=gp_zpos(p_gp)dist=sqrt(px*px+py*py+pz*pz)if dist>0 then; 节点位置调整maxp=max(px,max(py,pz))k=(maxp/rad)*(rad/dist)gp_xpos(p_gp)=k*pxgp_ypos(p_gp)=k*pygp_zpos(p_gp)=k*pzend_ifp_gp=gp_next(p_gp)end_loopendmake_sphere; ===============================; 利用镜像生成完整球体网格; =============================== gen zone refgen zone ref dip 90gen zone ref dip 90 dd 90; =============================== ; 显示球体网格; =============================== plot surf pl set back whpl bl grngen zon bri size 1 1 2 group soil ran z 1 1 group rock ran z 01 expgrid 1.flac3dgen zon bri size 3 3 3model mohrprop bu 3e6 sh 1e6 coh 10e3 fric 15 fix z ran z -.1 .1fix x ran x -.1 .1fix x ran x 2.9 3.1fix y ran y -.1 .1fix y ran y 2.9 3.1ini dens 2000hist unbalset grav 10solve elasticsave 6-1.savrest 6-1.savini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0app nstress -100e3 ran z 2.9 3.1 x 1 2 y 1 2 solvesave 6-2.savrest 6-1.savini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0app nstress -100e3 ran z 2.9 3.1 x 1 2 y 1 2 hist id=2 gp zdis 1 1 3hist id=3 gp zdis 1 1 2hist id=4 gp xdis 1 1 3hist id=5 gp xdis 1 1 3hist id=6 zone szz 1 1 3hist id=7 zone szz 1.5 1.5 2.5hist id=8 zone sxz 1.5 1.5 2.5solvesave 6-3.savrest 6-3.savset log onset logfile 6-2.log print zone stress print gp disset log offrest 6-1.savini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0app nstress -100e3 ran z 2.9 3.1 x 1 2 y 1 2 plot set rot 20 0 30plot con szz ou on magf 10plot add hist 1set movie avi step 1 file 6-5.avimovie startsolvemovie finishngen zon bri size 1 1 2 model elasprop bulk 3e7 shear 1e7 fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000set grav 0 0 -10solveplo con szngen zon bri size 1 1 2model mohrprop bulk 3e7 shear 1e7 c 1e10 f 15 tension 1e10 fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000set grav 0 0 -10solveprop bulk 3e7 shear 1e7 c 10e3 f 15 ten 0solveplo con szngen zone brick size 1 1 2model mohrprop bulk 3e7 shear 1e7 coh 10e3 fri 15 ten 0 fix z ran z 0 fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000set grav 0 0 -10solve elasplo con sznewgen zone brick size 1 1 2model mohrprop bulk 3e7 shear 1e7 coh 10e3 fri 15 ten 0 fix z ran z 0 fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000ini szz -40e3 grad 0 0 20e3 ran z 0 2ini syy -20e3 grad 0 0 10e3 ran z 0 2ini sxx -20e3 grad 0 0 10e3 ran z 0 2set grav 0 0 -10solveplo con szngen zon bri size 1 1 2model mprop bulk 3e7 shear 1e7 c 10e10 f 15 ten 1e10 fix z ran z 0 fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000 ran z 0 1ini dens 1500 ran z 1 2ini szz -35e3 grad 0 0 20e3 ran z 0 1ini syy -17.5e3 grad 0 0 10e3 ran z 0 1ini sxx -17.5e3 grad 0 0 10e3 ran z 0 1ini szz -15e3 grad 0 0 15e3 ran z 1 2ini syy -7.5e3 grad 0 0 7.5e3 ran z 1 2ini sxx -7.5e3 grad 0 0 7.5e3 ran z 1 2ini pp 10e3 grad 0 0 -10e3 ran z 0 1set grav 0 0 -10solveplo con szngen zon bri size 1 1 2model mprop bulk 3e7 shear 1e7 c 10e10 f 15 ten 1e10 fix z ran z 0fix x ran x 0fix x ran x 1fix y ran y 0fix y ran y 1ini dens 2000 ran z 0 2ini szz -50e3 grad 0 0 20e3 ran z 0 1ini syy -30e3 grad 0 0 10e3 ran z 0 1ini sxx -30e3 grad 0 0 10e3 ran z 0 1ini pp 30e3 grad 0 0 -10e3 ran z 0 2app nstress -10e3 ran z 2set grav 0 0 -10solveplo con sznewgen zone brick p0 0 0 0 p1 60 0 0 p2 0 60 0 p3 0 0 90 & p4 60 60 0 p5 0 60 90 p6 60 0 150 p7 60 60 150 & size 6 6 10model elaspro bulk 10e10 she 10e10ini den 2500apply sxx -1e9 grad 0 0 1.1111111e7 range x -.1 .1 apply sxx -1e9 grad 0 0 6.6666666e6 range x 59.9 60.1 apply syy -1e9 grad 0 0 8.3333333e6 range y -.1 .1 apply syy -1e9 grad 0 0 8.3333333e6 range y 59.9 60.1 apply szz -1e8 grad 0 0 8.3333333e5 ran z 0 120set grav 0 0 -10step 30000ini xdisp 0 ydisp 0 zdisp 0ini xvel 0 yvel 0 zvel 0plo cont szznewgen zone brick p0 0 0 0 p1 60 0 0 p2 0 60 0 p3 0 0 90 & p4 60 60 0 p5 0 60 90 p6 60 0 150 p7 60 60 150 & size 6 6 10model elaspro bulk 10e10 she 10e10ini den 2500ini sxx -1e9 grad 0 0 1.1111111e7 range x -.1 .1ini sxx -1e9 grad 0 0 6.6666666e6 range x 59.9 60.1ini syy -1e9 grad 0 0 8.3333333e6 range y -.1 .1ini syy -1e9 grad 0 0 8.3333333e6 range y 59.9 60.1ini szz -1e8 ran z -.1 .1fix x y z ran z -.1 .1set grav 0 0 -10solveini xdisp 0 ydisp 0 zdisp 0ini xvel 0 yvel 0 zvel 0plo cont szz。

FLAC3D的实例应用分析

FLAC3D的实例应用分析首先是岩土工程领域。

FLAC3D可以用于模拟岩土体的力学行为，预测在不同荷载作用下的岩土体变形和破坏，为设计和施工提供依据。

例如在基岩边坡稳定性分析中，FLAC3D可以模拟边坡在自然的和工程加载下的变形和破坏，评估边坡的稳定性，并优化边坡设计。

另外，FLAC3D还可以用于模拟土体动力响应，预测地震荷载下土体的动力特性和地震响应，为抗震设计提供参考。

其次是矿产资源开发领域。

FLAC3D可以模拟矿山开采过程中岩体的破坏和变形，评估开采对周围环境的影响，提供合理的采矿方案。

比如在隧道开挖中，FLAC3D可以模拟隧道的开挖和支护过程，评估围岩的稳定性，指导隧道支护设计和施工。

此外，FLAC3D还可以用于矿山坍塌、局部塌陷和裂隙水压力分布等现象的模拟与分析。

第三是地下空间开发领域。

FLAC3D可以模拟地下空间的开挖、支护和使用过程，预测开挖对周围建筑物的影响，评估地下空间的稳定性和安全性。

例如在地铁隧道施工中，FLAC3D可以模拟盾构掘进和地面沉降过程，评估地下水位、水压及地表沉降对周围土体的影响，指导施工方案的调整与优化。

最后是地质灾害研究领域。

FLAC3D可以模拟地质灾害的发生过程，了解其机理和演化规律，评估灾害对人类和环境的影响，提出相应的防灾措施。

例如在滑坡研究中，FLAC3D可以模拟土体的滑动过程，预测滑坡位置、速度和影响范围，为滑坡防治提供科学依据。

此外，FLAC3D还可以用于模拟地震、火山喷发和地下水位变化等灾害事件的发生和演化。

综上所述，FLAC3D在岩土工程、矿产资源开发、地下空间开发和地质灾害研究等领域有着广泛的应用。

它的模拟能力和计算精度使其成为解决实际问题的重要工具，为工程设计和决策提供准确、可靠的技术支持。

FLAC3D工程实例之MILESFIND_ADDCONT_EXCAVPROCAV

FLAC3D工程实例之MILESFIND_ADDCONT_EXCAVPROCAV1、FLAC3D工程实例之MILES.dat（calculate and store Von Mises stress in zone extension1）config zextra 1def mises; --- calculate and store Von Mises stress in zone extension 1 ---p_z = zone_headmax_mises = 0.0loop while p_z # nullmstr = (sxx(p_z) + syy(p_z) + szz(p_z)) / 3.dsxx = sxx(p_z) - mstrdsyy = syy(p_z) - mstrdszz = szz(p_z) - mstrdsxy = sxy(p_z)dsxz = sxz(p_z)dsyz = syz(p_z)vmstr2 = 1.5 * (dsxx*dsxx + dsyy*dsyy + dszz*dszz)vmstr2 = vmstr2 + 3. * (dsxy*dsxy + dsxz*dsxz + dsyz*dsyz) if vmstr2 > 0.0 thenz_extra(p_z,1) = sqrt(vmstr2)elsez_extra(p_z,1) = 0.0endifmax_mises = max(max_mises,z_extra(p_z,1))p_z = z_next(p_z)end_loopendmisesplot cont zextra 1 alias 'Von Mises Stress' averageprint max_mises2、FLAC3D工程实例之FIND_ADD.dat（find gridpoints at pile cap）; find gridpoints at pile cap;def find_addtop_head = nullgp_pnt = gp_headloop while gp_pnt # nullif gp_zpos(gp_pnt) > 0.05 thennew = get_mem(2)mem(new) = top_headmem(new+1) = gp_pnttop_head = newendifgp_pnt = gp_next(gp_pnt)endloopendfind_add3、FLAC3D工程实例之CONT_EX.dat（完成当前开挖继续以下开挖）; finishes stepping for the current excavation step; and continues following excavation steps;def cont_excavrest_step = 3000*(step/3000+1)-stepif rest_step = 3000 thenrest_step = 0end_ifcommandstep rest_stepend_command;store displacements in tablesxtable(1,cut) = 3.0 * cutytable(1,cut) = zdisp(ipt_surf)xtable(2,cut) = 3.0 * cutytable(2,cut) = zdisp(ipt_crown)xtable(3,cut) = 3.0 * cutytable(3,cut) = zdisp(ipt_spring)commandsave m1.savend_commandif cut=10 thencommandsave m1_30.savend_commandend_ifexcavend4、FLAC3D工程实例之cav.pro（剑桥模型介质的圆柱形空腔膨胀排水结果）;-----------------------------------------------------------; cav.pro; results from undrained cylindrical cavity expansion; in Cam-Clay medium;-----------------------------------------------------------res cav.sav; table 10 : pore pressure; table 11 : radial-component of effective stress (xx-stress); table 12 : axial-component of effective stress (yy-stress); table 13 : tangential-component of effective stress (zz-stress) ; rval : ln(r/a0)def b_tablepnt = zone_headloop while pnt # nullxval = xcen(pnt)yval = ycen(pnt)rval = 0.5*ln(xval*xval+yval*yval)c_pp = z_pp(pnt)table(10,rval) = c_pptable(11,rval) = -(sxx(pnt)*cc+syy(pnt)*ss+sxy(pnt)*sc) - c_pp table(12,rval) = -(sxx(pnt)*ss+syy(pnt)*cc-sxy(pnt)*sc) - c_pp table(13,rval) = -szz(pnt) - c_pppnt = z_next(pnt)end_loopendb_table; results ---plot his 9 vs 4 ; Fig 9.3plot his -2 3 vs 4 ; Fig 9.4plot table 10 line 11 line 12 line 13 line ; Fig 9.5ret5、FLAC3D工程实例之CAV（剑桥模型介质的圆柱形空腔膨胀排水）;----------------------------------------------------------- ; cav.dat; Undrained cylindrical cavity expansion in Cam-Clay medium;----------------------------------------------------------- title Undrained cylindrical cavity expansion in Cam-Clay medium config fluid; --- model geometry ---gen zone brick p0 1.0 0.0 1.0 p1 200.0 0.0 1.0 &p2 1.0 0.0 0.0 p4 200.0 0.0 0.0 &p3 0.9877 0.1564 1.0 p6 197.5377 31.2869 1.0 &p5 0.9877 0.1564 0.0 p7 197.5377 31.2869 0.0 &size 31 1 1 ratio 1.1 1 1def find_idpnt1 = gp_near(1.,0.,1.)pnt2 = gp_near(1.,0.,0.)pnt3 = gp_near(0.9877,0.1564,1.)pnt4 = gp_near(0.9877,0.1564,0.)id1 = gp_id(pnt1)id2 = gp_id(pnt2)id3 = gp_id(pnt3)id4 = gp_id(pnt4)pntz = z_near(1.,0.07,0.5)cc = cos(4.5*degrad)ss = sin(4.5*degrad)sc = 2.*cc*sscc = cc*ccss = ss*ssb_vel= 1.e-5b_vx = 0.9877*b_velb_vy = 0.1564*b_velendfind_idrange name out cyl end1 0 0 -10 end2 0 0 10 rad 197.5377not range name in cyl end1 0 0 -10 end2 0 0 10 rad 1.001 ; --- model properties ---mo cam-clayprop shear 74.prop mm 1.2 lambda 0.15 kappa 0.03 mp1 1.0 mv_l 2.3mo fl_iso; --- boundary conditions ---fix y z range y -.001 .001apply nvel 0.0 plane dip 90.0 dd -9.0 origin 0.0 0.0 0.0 &range in not y 0.07 200.0apply dvel 0.0 plane dip 90.0 dd -9.0 origin 0.0 0.0 0.0 &range in not y 0.07 200.0apply nstress -1.65 range outfix x y z range inini xvel b_vel range id id1 id1ini xvel b_vel range id id2 id2ini xvel b_vx yvel b_vy range id id3 id3ini xvel b_vx yvel b_vy range id id4 id4ini sxx -1.65 syy -1.65 szz -3.; model settings ---set fluid biot onset fluid offset large; --- fish functions ---; ... initial specific volume, tangent bulk modulus, porosity ... def set_proppnt = zone_heads1 = -z_sxx(pnt)s2 = -z_syy(pnt)s3 = -z_szz(pnt)p0 = (s1 + s2 + s3) / 3.q0 = sqrt(((s1-s2)*(s1-s2)+(s2-s3)*(s2-s3)+(s3-s1)*(s3-s1))*0.5) e0 = q0 / (z_prop(pnt,'mm') * p0)pc = p0 * (1. + e0 * e0)p1 = z_prop(pnt,'mp1')vc = z_prop(pnt,'mv_l') - z_prop(pnt,'lambda') * ln(pc/p1) v0 = vc - z_prop(pnt,'kappa') * ln(p0/pc)b0 = v0 * p0 / z_prop(pnt,'kappa')n0 = (v0 - 1.) / v0bb = 10. * b0bi = 100. * b0commandprop mpc pc bulk_bound bbini biot_mod biend_commandend; ... numerical variables ...def pathsrad = z_sxx(pntz)*cc+z_syy(pntz)*ss+z_sxy(pntz)*scp_fl = z_pp(pntz)c_tp = -(z_sxx(pntz)+z_syy(pntz)+z_szz(pntz)) / 3.c_ep = z_prop(pntz,'cp')c_q = z_prop(pntz,'cq')c_cq = c_ep * z_prop(pntz,'mm')c_v = z_prop(pntz,'cv')c_a = 1. + gp_xdisp(pnt2)c_b = z_prop(pntz,'bulk')c_pc = z_prop(pntz,'mpc')end; --- histories ---hist nstep 500 hist pathhist sradhist p_flhist c_ahist c_ephist c_cqhist c_tphist c_vhist c_qhist c_bhist c_pchist unbal; --- test --- set_prop step 100000 save cav.sav。

FLACFLAC3D基础与工程实例_记录

《FLACFLAC3D基础与工程实例》阅读札记目录一、FLACFLAC3D软件概述 (2)1. 软件背景与简介 (3)1.1 FLACFLAC3D的发展历程 (4)1.2 软件的应用领域及特点 (5)2. 软件安装与运行环境 (6)2.1 系统要求 (7)2.2 安装步骤 (8)2.3 运行环境配置 (10)二、FLACFLAC3D基础知识 (11)1. 基本概念与术语 (13)1.1 有限元分析原理 (14)1.2 离散元法简介 (14)1.3 FLACFLAC3D中的相关术语解释 (15)2. 软件操作界面及功能模块 (17)2.1 操作界面介绍 (18)2.2 主要功能模块说明 (20)2.3 菜单功能详解 (20)三、工程实例分析 (22)1. 地质工程实例 (23)1.1 工程背景及问题定义 (25)1.2 模型建立与参数设置 (26)1.3 结果分析与讨论 (27)2. 土木工程实例 (29)2.1 工程概况与建模目的 (30)2.2 建模过程及计算步骤 (31)2.3 结果展示与工程应用 (32)四、FLACFLAC3D应用技巧与注意事项 (33)1. 建模技巧与优化方法 (34)1.1 建模策略及优化思路 (35)1.2 网格划分与模型简化技巧 (36)1.3 参数设置与模型验证方法 (38)2. 数据分析与处理方法 (40)2.1 数据采集与整理方法 (41)2.2 结果分析与图表展示技巧 (42)一、FLACFLAC3D软件概述3D是一种广泛使用的岩土力学与有限元分析软件。

它是一套专门用来分析连续介质中的物理力学现象的强大工具，主要应用于土木、矿山、隧道等领域，能针对各种复杂的工程问题进行数值建模和模拟分析。

3D以其高效、灵活的数值分析能力，为工程师提供了强大的技术支持。

其主要特点包括：多功能：3D能够模拟多种物理过程，包括应力分析、稳定性分析、流体流动分析等，适用于多种工程场景。

利用FLAC3D对山西某电厂单桩静载试验模拟

利用FLAC3D对山西某电厂单桩静载试验模拟近年来，随着我国经济的不息进步和工业化进程的推行，电厂建设变得越来越重要。

山西某电厂作为该地区的重要能源供应单位，对电厂的建设质量和安全性有着严格的要求，尤其是电厂的地基承载能力。

为了准确评估电厂基础的承载能力，单桩静载试验成为不行或缺的一环。

本文将介绍利用FLAC3D软件对山西某电厂单桩静载试验进行模拟的探究。

一、探究背景山西某电厂的单桩静载试验旨在评估桩基的承载能力和变形状况，为后续的设计和施工提供依据。

传统的现场试验方法费时费劲，而且存在一定的风险，因此利用数值模拟方法进行试验模拟成为一种重要的探究手段。

二、试验方案设计山西某电厂的试验方案设计如下：1. 试验桩径为1m，桩长为30m，桩身接受混凝土浇筑；2. 试验载荷沿桩长方向逐渐增加，并保持静载；3. 在试验过程中对桩体的应力、变形等参数进行监测和记录。

三、FLAC3D软件建模基于FLAC3D软件的数值模拟方法可以对试验的力学行为进行准确的猜测和分析。

以下是FLAC3D建模的步骤：1. 导入电厂的工程地质资料，包括土层厚度、强度参数等；2. 对土体力学参数进行划分，将地层分为不同的单元；3. 设定桩的几何外形和材料特性，并与土体单元互相毗连；4. 在桩顶施加试验载荷，并监测各关键位置的应力和变形。

四、数值模拟结果分析通过FLAC3D软件的数值模拟，可以得出以下结果：1. 桩身变形特征：随着试验载荷的增加，桩身会发生不同程度的弯曲和沉降，呈现出非线性的变形特征；2. 桩身应力分布：试验载荷导致桩身上下部分产生不同的应力分布，载荷较大的位置应力集中；3. 土体响应变化：土体受到桩身载荷的作用，会发生一定的变形和应力分布，但相对桩身而言，土体的变形较小。

五、结果验证和谈论与现场试验结果对比显示，FLAC3D软件模拟的单桩静载试验具有较高的模拟精度。

桩身的变形和应力分布与实际试验数据吻合较好，验证了数值模拟的可行性。

FLAC3D案例之混凝土桩的轴向和侧向荷载

FLAC3D案例之混凝土桩的轴向和侧向荷载一、问题描述：计算单个混凝土桩基础在轴向和横向荷载下的荷载-挠度响应。

首先对桩施加100kN 的轴向载荷，然后桩的顶部最大水平位移为4cm。

确定桩体轴向荷载与极限承载力的关系，计算桩身侧向荷载-位移曲线。

桩身直径0.6米，长度5米，嵌在均匀的粘土层中。

地下水位为5.5米。

混凝土桩和粘土的性能如表1所示。

二、建模过程：通过桩轴的垂直平面是本模型的对称平面。

FLAC3D模型的坐标轴位于桩顶原点处，z轴沿桩轴向上。

模型网格如图1所示。

模型的顶部z=0处是一个自由曲面，模型底部z=-8m处为z向固定，模型两侧|x|=8m处为滚动边界，y=8m处为滚动边界。

桩的轴向承载力是桩身表面摩阻力和桩端承载力的函数。

通过在桩墙和粘土之间设置一个界面来模拟桩墙的表面摩擦阻力。

界面的摩擦和粘聚特性代表了混凝土与粘土之间的摩擦阻力。

对于这个示例，摩擦20°角和凝聚力30 kPa被假定为接口的属性。

在桩顶和粘土之间设置第二个界面。

Figure 1: FLAC3D grid for vertical and lateral loading of a concrete pile in clay.为模拟更加真实的工程环境，采用了两种不同的界面：一种在桩壁上，另一种在桩基础上。

在前面的命令中，区域面是分开的，这样两个通用的网格点也会分开。

在桩身安装前，首先将模型引入重力荷载作用下的平衡应力状态。

在z=-5.5m处创建一个水平的地下水位，粘土的湿密度被分配到地下水位以下的区域。

在下一步的分析中，模型在桩身安装完成后进入平衡状态。

通过将桩带的性质由代表粘土材料的性质改为代表混凝土桩材料的性质来模拟安装。

平衡状态下竖向应力分布，包括桩身重量，如图2所示。

Figure 2: Contours of vertical stress at the initial stress-state, including the weight of the pile.通过在桩顶施加竖向速度计算桩的轴向极限承载力。

FLAC3D实例分析教程(2024)

温度场建模
支持多种温度场建模方法，如有限元法、有限差分法等，可自定义温度边界条件和初始条件。
热传导过程模拟
通过数值计算方法求解热传导方程，得到温度场中的温度分布和热流密度等参数。
热应力计算
基于热弹性力学理论，计算由温度梯度引起的热应力和变形，以及温度对岩土体力学性质的影响。
29
07
总结与展望
地下工程
软件可模拟地下洞室、巷道、地铁等地下工程的开挖和支护过程，分析围岩稳定性和支护结构受力情况。
基坑工程
软件可用于分析基坑开挖过程中的土压力分布、支护结构变形和稳定性等问题。
2024/1/30
隧道工程
FLAC3D可应用于隧道工程的开挖、支护和衬砌设计，评估隧道施工对周围岩体的影响。
5
操作界面及基本设置
边界条件
设定模型的底部和侧面为固定边界，顶部为自由边界；考虑地下水的影响，设置相应的孔隙水压力和渗透系数。
2024/1/30
18
求解过程与结果分析
求解过程
采用FLAC3D内置的求解器进行计算，包括初始地应力平衡和后续加载过程。
结果分析
提取边坡的位移、应力、应变等计算结果，分析边坡的变形和破坏模式。
2024/1/30
操作界面
FLAC3D的操作界面包括菜单栏、工具栏、模型视图、属性视图和输出窗口等部分，方便用户进行建模、分析和后处理等操作。
基本设置
在使用FLAC3D前，需要进行一些基本设置，如选择单位制、设置材料参数、定义边界条件等。这些设置将直接影响模拟结果的准确性和可靠性。
建模流程
FLAC3D的建模流程包括建立几何模型、划分网格、定义材料属性、施加边界条件和初始条件等步骤。用户可根据实际需求选择合适的建模方法和工具。

FLAC3D的实例应用分析

矿震震级 11.5
(M L )
矿震数目 5721
(次 )
1.62.0 1675
2.12.5 2.63.0 3.13.6
298
63
19
合计 7776
矿震震级、频度与时间特征
1100m
F25
z y x
计算模型
F1-A F1
F26
3800m
F16 F18
F7-1 3000m
1999-2000
beam
cable
2. 锚索(cable)单元
pile
3. 桩(pile)单元
4. 壳(shell)单元 5. 格栅(geogrid)单元
shell
土工织物；土工格栅
6. 衬砌(liner)单元
geogrid
liner
1.5﹑FLAC3D接触单元
接触单元原理
1. 三角形单元(无厚度!) 2. 三种工作模式
• FLAC3D适用模拟计算岩土材料力学行为，特别适合模拟大变形和扭曲，包括材料的高度非线性(应变硬化/软化)、不可逆剪切破坏和压密、粘弹(蠕变)、孔隙介质的应力—渗流耦
合、热—力耦合以及动力学问题等。
1﹑FLAC3D软件简介
4. 采用“混合离散法”用以精确模拟塑性坍塌和塑性流动。这种方法比有限元法中采用的渐进迭代更为有效。
-20
X0
Y
10
20
30
40
50
60
70
80
提
纲
一﹑FLAC3D软件简介二﹑FLAC3D应用实例三﹑FLAC3D软件应用四﹑FLAC3D模拟技巧
3.1、建模与模拟过程
基本步骤 1. 根据研究目的对实际模型进行构思与概化，计算模型

FLAC3D实例分析教程

FLAC3D实例分析教程假设我们要分析一个简单的边坡稳定性问题。

下面是具体的步骤：1.建立几何模型：首先，我们需要建立一个几何模型，包括边坡的形状和岩土层的属性。

在FLAC3D中，我们可以通过在网格上定义顶点和连线来创建边坡的形状。

然后，我们可以设置每个区域的岩土层属性，如密度、强度和摩擦角等。

确保模型的几何和岩土层属性与实际情况相符。

2.设定边界条件：接下来，我们需要设定边界条件，即模拟中的约束和加载条件。

在边坡稳定性问题中，我们可以设定边坡底部的约束条件，如水平位移和垂直位移。

此外，我们还可以为边坡施加水平和垂直方向的荷载，模拟边坡于不同加载条件下的行为。

3.运行模拟：在完成模型和边界条件的设置后，我们可以开始运行模拟。

FLAC3D使用多线程计算，能够利用多核处理器的能力来进行快速计算。

我们可以选择设置时间步长和计算精度等参数。

模拟运行完毕后，FLAC3D将输出边坡在不同加载条件下的应力、位移和变形等结果。

4.结果分析：最后，我们需要对模拟结果进行分析和解释。

FLAC3D 提供了丰富的结果显示和分析功能。

我们可以通过绘制曲线图、生成动画和查看计算网格等方式来可视化和分析结果。

根据模拟结果，我们可以评估边坡的稳定性，并提出针对性的建议和改进方案。

在实际应用中，我们还可以使用FLAC3D的其他高级功能来进一步分析和优化边坡设计。

例如，我们可以引入土体的非线性行为模型，模拟地下水流和渗流等复杂的工程问题。

此外，FLAC3D还支持参数化建模和优化分析，可以帮助工程师迅速评估不同方案的可行性和性能。

总结起来，FLAC3D是一个强大的岩土工程分析软件，可以用于解决各种实际问题。

通过学习和应用FLAC3D的基本使用方法和分析技巧，工程师可以更好地理解和评估岩土工程问题，为工程设计和施工提供有力支持。

FLAC3D实例分析教程2

FLAC3D实例分析教程2FLAC3D实例分析教程2实例说明：在一个地下隧道工程中，我们需要分析围岩的稳定性。

隧道的尺寸为10mx6mx20m，围岩由砂岩组成，其物理特性如下：- 密度：2.5g/cm³-弹性模量：40GPa-泊松比：0.25我们将在FLAC3D中建立一个三维模型，并进行围岩的稳定性分析。

步骤1：建立模型在FLAC3D中，首先需要创建一个新的项目文件。

点击“File”-“New”-“Project”来创建一个新的项目文件，并保存为适当的文件名。

然后，点击“Grid”-“Generate”来生成一个新的网格。

在对话框中，输入隧道的尺寸，并选择合适的网格密度。

点击“Apply”来生成网格。

步骤2：定义围岩属性在FLAC3D中，可以通过定义不同的材料属性来模拟不同的岩石类型。

点击“Model”-“Material”-“New”来定义一个新的材料，并设置其物理特性。

在对话框中，输入材料的密度、弹性模量和泊松比。

点击“OK”来保存材料属性。

步骤3：生成围岩在FLAC3D中，可以通过定义不同的围岩属性来模拟围岩中的不同部分。

点击“Grid”-“Approval”来选择需要定义材料的单元，并在对话框中选择刚刚定义的材料。

点击“Apply”来应用材料属性。

步骤4：施加边界条件在FLAC3D中，可以通过定义不同的边界条件来模拟不同的荷载情况。

在本实例中，我们将施加一个围压荷载，并固定隧道的底部。

点击“Model”-“Boundary Condition”来定义边界条件。

在对话框中，选择围压荷载并输入荷载大小。

点击“Apply”来应用边界条件。

步骤5：运行模拟在FLAC3D中，可以通过点击“Model”-“Run”来运行模拟。

在运行模拟之前，可以选择运行的时间步长、计算方法和收敛准则。

点击“OK”来开始运行模拟。

步骤6：结果分析在FLAC3D中，可以通过查看不同的结果图来分析模拟结果。

点击“Post”-“Plot”来选择需要查看的结果图，并选择合适的结果类型。

FLAC及FLAC3D基础与工程实例51-2

FLAC及FLAC3D基础与工程实例51第1章FLAC、FLAC3D的功能与特性自R.W；数值模拟技术的优势在于有效延伸和扩展了分析人员的；本章重点：；?FLAC/FLAC3D的主要特点；?FLAC/FLAC3D的不足之处；1.1FLAC/FLAC3D简介；FLAC （FastLagrangianAnaly；FLAC有二维和三维计算软件两个版本，即FLAC；1.2FLAC/FLAC3D的主要特点；F第1章 FLAC、FLAC3D的功能与特性自R.W. Clough 1965年首次将有限元引入土石坝的稳定性分析以来，数值模拟技术在岩土工程领域获得了巨大的进步，并成功解决了许多重大工程问题。

数值模拟技术的优势在于有效延伸和扩展了分析人员的认知范围，为分析人员洞悉岩、土体内部的破坏机理提供了强有力的可视化手段。

因此，优秀的岩土工程数值模拟软件须在专业性、可视化及信息输出等方面做到相对完备，方能使分析人员专注于工程实际问题的研究、分析和解决。

FLAC 系列软件的出现，为岩土工程研究工作者提供了一款功能强大的数值模拟工具。

本章重点：? FLAC/FLAC3D的主要特点? FLAC/FLAC3D的不足之处1.1 FLAC/FLAC3D简介FLAC（Fast Lagrangian Analysis of Continua）是由Itasca公司研发推出的连续介质力学分析软件，是该公司旗下最知名的软件系统之一。

FLAC目前已在全球七十多个国家得到广泛应用，在国际土木工程（尤其是岩土工程）学术界和工业界享有盛誉。

FLAC有二维和三维计算软件两个版本，即FLAC2D（1984）和FLAC3D（1994）。

这里进行一下说明，本书在阐述软件系列时，以FLAC统一称谓FLAC2D和FLAC3D；分述FLAC2D和FLAC3D时，FLAC仅指代FLAC2D。

FLAC3D隧道耦合实例

作一个盾构隧道开挖，地下水位设为z=12.6m，（隧道中心线设为z=0m）。

分了7层土，给定初始条件的代码如下：config fluid（建模过程省）model fl_isotropicprop perm 1e-12 poro 0.48 range z 9.1 12.6prop perm 1e-12 poro 0.55 range z -0.9 9.1prop perm 1e-12 poro 0.5 range z -7.4 -0.9prop perm 1e-12 poro 0.5 range z -11.4 -7.4prop perm 1e-12 poro 0.45 range z -18.4 -11.4prop perm 1e-12 poro 0.45 range z -30.4 -18.4ini fdens 1e3ini fmod 8.5e7ini sat 0 range z 12.6 15.1ini sat 1 range z -30.4 12.6ini pp 0 grad 0 0 -1e4 range z -30.4 12.6fix pp range x -.1 .1fix pp range x 39.9 40.1fix pp range y -.1 .1fix pp range y 119.9 120.1fix pp range z 12.5 12.7fix pp range z -30.5 -30.3;material mechanic parametersmodel mohrdef derives_mod1=E_mod1/(2.0*(1.0+p_ratio1))b_mod1=E_mod1/(3.0*(1.0-2.0*p_ratio1)) s_mod2=E_mod2/(2.0*(1.0+p_ratio2))b_mod2=E_mod2/(3.0*(1.0-2.0*p_ratio2))s_mod3=E_mod3/(2.0*(1.0+p_ratio3))b_mod3=E_mod3/(3.0*(1.0-2.0*p_ratio3))s_mod4=E_mod4/(2.0*(1.0+p_ratio4))b_mod4=E_mod4/(3.0*(1.0-2.0*p_ratio4))s_mod5=E_mod5/(2.0*(1.0+p_ratio5))b_mod5=E_mod5/(3.0*(1.0-2.0*p_ratio5))s_mod6=E_mod6/(2.0*(1.0+p_ratio6))b_mod6=E_mod6/(3.0*(1.0-2.0*p_ratio6))s_mod7=E_mod7/(2.0*(1.0+p_ratio7))b_mod7=E_mod7/(3.0*(1.0-2.0*p_ratio7))endset E_mod1=2.40e6 p_ratio1=0.25 E_mod2=5.0e6 p_ratio2=0.30 E_mod3=11.0e6 p_ratio3=0.25 &E_mod4=8.5e6 p_ratio4=0.25 E_mod5=11.5e6 p_ratio5=0.27 E_mod6=9.25e6 p_ratio6=0.22 &E_mod7=12.40e6 p_ratio7=0.2deriveprop bulk b_mod5 shear s_mod5 cohe 18.0e3 fric 18.0 ten 55.398e3range z 12.6 15.1prop bulk b_mod3 shear s_mod3 cohe 6.3e3 fric 21.0 ten 16.412e3 range z 9.1 12.6prop bulk b_mod1 shear s_mod1 cohe 13.2e3 fric 10.0 ten 74.861e3 range z 0 9.1prop bulk b_mod1 shear s_mod1 cohe 13.2e3 fric 10.0 ten 74.861e3 range z -0.9 0prop bulk b_mod2 shear s_mod2 cohe 15.3e3 fric 10.0 ten 86.771e3 range z -7.4 -0.9prop bulk b_mod4 shear s_mod4 cohe 22.0e3 fric 20.0 ten 60.445e3 range z -11.4 -7.4prop bulk b_mod6 shear s_mod6 cohe 3.0e3 fric 25.0 ten 6.434e3 range z -18.4 -11.4prop bulk b_mod7 shear s_mod7 cohe 3.0e3 fric 25.0 ten 6.434e3 range z -30.4 -18.4;boundary conditionsfix x range x -0.1 0.1fix x range x 39.9 40.1fix y range y -0.1 0.1fix y range y 119.9 120.1fix x range z -30.5 -30.3fix y range z -30.5 -30.3fix z range z -30.5 -30.3interface 1 prop kn 3e9 ks 1e9 fric 20 coh 3e5 interface 2 prop kn 3e9 ks 1e9 fric 20 coh 3e5;stress conditionsset grav 0 0 -10ini dens 1.87e3 range z 12.6 15.1ini dens 1.87e3 range z 9.1 12.6ini dens 1.76e3 range z 0 9.1ini dens 1.76e3 range z -0.9 0ini dens 1.84e3 range z -7.4 -0.9ini dens 2.0e3 range z -11.4 -7.4ini dens 1.89e3 range z -18.4 -11.4ini dens 1.93e3 range z -30.4 -18.4ini szz -28.237e4 grad 0 0 1.87e4 range z 12.6 15.1 ini szz -28.237e4 grad 0 0 1.87e4 range z 9.1 12.6ini szz -27.236e4 grad 0 0 1.76e4 range z 0 9.1ini szz -27.236e4 grad 0 0 1.76e4 range z -0.9 0ini szz -27.164e4 grad 0 0 1.84e4 range z -7.4 -0.9 ini szz -25.98e4 grad 0 0 2.0e4 range z -11.4 -7.4ini szz -27.234e4 grad 0 0 1.89e4 range z -18.4 -11.4 ini szz -26.498e4 grad 0 0 1.93e4 range z -30.4 -18.4ini sxx -28.237e4 grad 0 0 1.87e4 range z 12.6 15.1 ini sxx -28.237e4 grad 0 0 1.87e4 range z 9.1 12.6 ini sxx -27.236e4 grad 0 0 1.76e4 range z 0 9.1ini sxx -27.236e4 grad 0 0 1.76e4 range z -0.9 0ini sxx -27.164e4 grad 0 0 1.84e4 range z -7.4 -0.9 ini sxx -25.98e4 grad 0 0 2.0e4 range z -11.4 -7.4ini sxx -27.234e4 grad 0 0 1.89e4 range z -18.4 -11.4 ini sxx -26.498e4 grad 0 0 1.93e4 range z -30.4 -18.4ini syy -14.905e4 grad 0 0 9.871e3 range z 12.6 15.1 ini syy -13.597e4 grad 0 0 8.833e3 range z 9.1 12.6ini syy -16.834e4 grad 0 0 1.239e4 range z 0 9.1ini syy -16.834e4 grad 0 0 1.239e4 range z -0.9 0ini syy -16.783e4 grad 0 0 1.296e4 range z -7.4 -0.9ini syy -19.117e4 grad 0 0 9.806e3 range z -11.4 -7.4ini syy -21.550e4 grad 0 0 7.672e3 range z -18.4 -11.4ini syy -21.252e4 grad 0 0 7.834e3 range z -30.4 -18.4（取控制点省）solvesave iniconditions.sav初始平衡后的PP如下图。

1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。

地下建筑结构
4.分步开挖支护计算模型
对基坑放坡开挖BC段,可以采取对开挖土体和基坑围岩土体分别建模,建模中用到基本单元的八节点定义方式。然后用null定义开挖土体部分的模块来模拟,如图4。
本工程实例中基坑坡率较大,近乎直立开挖,也可直接采取每步开挖支护中,再多次分步开挖得到阶梯状开挖面的方式近似模拟坡面开挖。经模拟计算验证,该建模方式对基坑开挖支护数值模拟计算结果影响很小,可忽略不计。
本工程沉降水准观测,对基坑周边边坡土体和建筑物沉降监测共布设10+8个沉降观测点。其中,在基坑周边土体共布设10个水平位移和沉降观测点,东面两栋距基坑较近的被测建筑物的首层柱上共为3个沉降观测点。其中对周边建筑物进行13次沉降观测。观测周期与工程进度密切联系,基坑开挖期间,土体扰动对周围环境影响较大,沉降速率较大,故保证每隔2天一次观测,时刻注意环境动态,后期施工过程中, 根据实际情况相对减小监测频率。
2、边界条件
边界条件分为位移边界条件和受力边界条件利用FLAC3D中fish语言程序命令:fix(位移方向)range(位移约束面),可在网格模型中设定位移边界条件。模型除x=0面及基坑顶面为自由面不设置位移边界条件外,其余面均采用法向约束。坑底边界静止不动,假设为固定铰支,限制三个方向的位移;模型x轴及y轴方向两侧避均施加边界约束条件,限制该临空方向的位移;基坑顶面为设为自由面,不加约束。地下水己进行处理,可不考虑。利用FLAC3D中fish语言程序命令:apply szz(荷载大小)range(荷载范围),可在网格模型中定义受力边界条件。本工程CD段基坑周边己有建筑物荷载,荷载大小为90KN。
图19第四步开挖支护x方向位移等值线图
图20第五步开挖支护x方向位移等值线图
图21第五步完成后CD段x方向位移典型剖面图
地下建筑结构
图中左侧表示x方向的位移。下面的数值为位移大小的变化范围,从上到下由负值到正值,负值代表向基坑内侧移动,正值代表背离基坑开挖面。对应于右侧图中同颜色区域的位移值。
图4 坡面开挖建模示意图
地下建筑结构
整个基坑开挖支护模拟严格按照工程施工顺序,分层分段开挖与支护。如表2所示,CD段分五步开挖,BC段分四步开挖,因CD段地面标高较高,故先开挖CD段第一层土体,然后开挖BC段第一层土体, 完成第一步开挖支护后,进行模拟计算,保存计算成果;然后依次进行CD段和BC第二步开挖支护,再次计算;依次类推完成整个基坑开挖支护过程。
地下建筑结构
5.3分步开挖水平位移分析对基坑分步开挖进行模拟分析。由于基坑位移控制方向为x和z方向,故水平位移只列出x方向的等值线图。分步开挖水平位移等值线图及典型剖面图见图16到21。
图16第一步开挖支护x方向位移等值线图
图17第二步开挖支护x方向位移等值线图
地下建筑结构
图18第三步开挖支护x方向位移等值线图
地下建筑结构
4、从基坑位移分布图结合不平衡力模拟计算结果知,基坑开挖隐患位于CD段,其主要危险因素为距基坑5.7m远的90KN 局部建筑荷载。因支护工艺采用小型机械施工,为充分反映周边建筑物对基坑稳定的相互影响,模拟过程中忽略了基坑
顶面的施工荷载以及雨水等不确定因素的影响,可知在基坑
BC段不平衡力很小,直接开挖产生的位移也可忽略不计。
初始地应力场生成的水平位移(x方向)和垂直位移(z方向)结果如图10、图11。
地下建筑结构
图10土体自然固结x方向位移
图11土体自然固结z方向位移
地下建筑结构
由结果可以看出:
1、x方向初始地应力固结的位移图中,前后两面因添加位移约束条件,故使得模拟值与固结位移值稍有偏差,主要表现在模型局部荷载区域靠近后侧x方向位移约束处,但总体而言,可以代表场地初始地应力场的水平固结位移分布。 2、x方向位移图表明:土体因承受局部建筑荷载而产生水平位移, 最大水平位移出现在荷载下方的土体表面,最大值仅为1.6mm。一般而言,在地面无局部荷载(或荷载沿整个场区表面均匀分布)且各土层为成平面各向异性体的理想情况下,土体自然固结不会产生水平方向的位移。故在基坑设计计算以及数值模拟分析中均不能忽略周边建筑物的荷载作用。 3、水平位移向外影响范围约为荷载宽度的一倍,远离荷载处水平位移非常小,可不计。模型前后两个边界因定义位移边界条件,故位移值为零。
级,地基等级为二级。
基坑大致呈矩形,东西向宽25.50m;南北向长45.00m,基坑开挖面积约1150.oom2,基坑周长约140m。地下建筑结构
1.2工程地质及水文地质条件
1、工程地质条件本基坑位于三级阶地。开挖后,侧壁中上部为填土和一般粘性土,部分坡段下部为老粘性土,部分坡段侧壁下部仍为一般粘性土, 坑底接近老粘性土,总体而言,侧壁中上部土质条件相对较差,坑底土质良好,对整体稳定有利。场地内地基岩土层划分为四个层组共六个亚层。 2、水文地质条件场地经人工整平后相对平坦。场地地下水分为上层滞水及砂土层中空隙承压水。（1）层杂填土,结构松散,属中等透水层,含少量上层滞水,但不能形成统一自由水位,受大气降水及地表径流补给; (2) 二层三层属弱透水层,一般起隔水作用。4、5层为粉砂夹粉土、粉质粘土,粉砂,属含水层,为空隙承压水。地下建筑结构
地下建筑结构
4、z方向固结沉降位移图表明:最大沉降值出现在局部荷载下方土体表面,最大固结位移约为3.47cm,且荷载对沉降位移的影响范围不大,边界土体初始地应力沉降固结位移约为Zcm。 5、因受土体自重作用,同一垂直平面内,土体固结沉降位移最大值
出现在土体顶面,因存在沉降位移的累积,随着深度增大位移值逐
M5累计水平位移在2.5mm左右,累计沉降值和水平位移发展趋
势均在容许范围内,且几乎没有差异沉降,不会对基坑边坡和
周边环境产生危害。满足基坑稳定和周边建筑物安全要求。地下建筑结构
3.模型和参数设置
l、计算范围
通过分析该深基坑工程的地质条件,周边环境及基坑要求,结合收集的资料,不考虑地下水的影响,对整个基坑的东面进行三维数值模拟。基坑东面分BC段和CD段两段:BC段长20m,基坑顶面标高-2.8m,基坑底面标高-9.4m,开挖深度6.6m,采用预应力锚杆支护结构,放坡坡率1:0.3;CD段长25m,基坑顶面标高-1.4m,基坑底面标高-6.85m,开挖深度5.45m,采用超前钢管桩复合预应力锚杆支护结构,垂直开挖。本计算模型在基坑宽度方向(x方向)取25m(开挖面后影响区取 15m,开挖基坑内影响区取10m),长度方向(y方向)取65m(两侧各加宽 10m影响区),深度方向(z方向)取20m。在建立FLAC3D计算模型时,
图4微型钢管桩超前支护图5第一步开挖与支护模型图图6第二步开挖与支护模型图图7第三步开挖与支护模型
地下建筑结构
图8第四步开挖与支护模型
图9第五步开挖完成后模型
地下建筑结构
5.数值模拟计算
5.1初始地应力场的生成
基坑开挖前,先进行土体的自然固结计算,即初始地应力场的生成,然后消除固结过程产生的位移。初始地应力场的生成选用分阶段弹塑性求解法进行生成。基本步骤为:生成网格模型--定义模型为摩尔--库伦弹塑性本构模型--设置土体强度参数、在模型上添加所有的物理边界条件 (包括位移边界条件和受力边界条件)--并设置土体密度和重力加速度--输入sofve elastic命令--软件按照默认设置自动分阶段求解---保存计算结果至initial.sav文件---自然固结位移清零。
渐变小,在模型底面边界,因定义边界约束条开挖位移分析通过零模型(null)来模拟基坑开挖,逐层开挖工况可以通过逐层设置开挖区材料为零模型来模拟。但是在无支护开挖时,分不开挖模拟和一次性开挖模拟结果差别不大,故在对
无支护开挖的基坑模拟计算直接通过一次性开挖到位后进行
地下建筑结构
表1周边建筑物水平位移及沉降监测成果汇总表
图1典型沉降观测点沉降与时间关系曲线图
地下建筑结构
图中: ①-⑤分别表示第一到五次分步开挖工况。从上图表可以看出:
1 、东侧已有建筑物位移随基坑开挖处于不断增大的趋
势,但随着基坑开挖支护完成后一段时间,位移趋于稳定。
2 、监测点 M7 、 M8 累计沉降量在 2.5mm 左右 , 基坑边 M4 、
河南城建学院
交通运输工程学院
地下建筑结构 L/O/G/O
1.工程概况及地质情况
1.1工程概况
湖北省高级人民法院位于武昌区首义路新建湖北省高级人民法院二号楼办公楼。办公楼高9层,总建筑面7922m2(不含地下室),采用框架结构,高38.50m,基坑开挖深度约3.85~7.05m。抗震设防分类为丙类,该工程重要性等级为二级,场地等级为三
地下建筑结构
2.2沉降观测成果该工程监测工作从2007年11月27日开始,至2008年3月24
日结束。分析中将第一次,累计观测时间未118天。将沉降监
测结果整理见表1。为了能更直观看出周围建筑物沉降随时间的变化规律, 对表中数据进行比较后,选取具有代表性的已有建筑物上M7、 M8点和基坑顶面坑边M4、M8点(见图1),将各点的沉降观测数据,绘制成随时间的沉降曲线图如图1。
地下建筑结构
2.位移监测成果
基坑监测由多方面内容组成,检测项目的选择应根据基坑工程的安全等级、周围环境的复杂程度和实际工作条件等因素而定。周边建筑物沉降观测是基坑周围环境监测其中的一个内容。在深基坑开挖过程中,为了掌握临近建筑物的沉降情况,确保周边环境的安全,需进行沉降观测。
2.1基坑监测方案
运算(solve)。基坑无支护开挖模拟水平、垂直位移及位移速率计算值见图12、13。
地下建筑结构
图12无支护开挖模拟水平位移计算
图13无支护开挖模拟水平位移速率计算
图14无支护开挖模拟沉降计算
图15无支护开挖模拟沉降速率计算