我的硕士论文《预应力锚索加固基坑的三维数值研究》，2007年被评为省级优硕论文

lmm

;此例说明需要删掉哪个铰结，只需要将其对应的节点的号码cid附上即可
;而附上id则删掉了这根锚索上所有的的铰结
;另外，在图中显示的被删掉的铰结为红颜色，未删的则为黑色
;还说明在.log文件中没有被注释的语句也被程序所示别了。
set log on
set logfile 999.log
gen zone brick size 10 10 10
sel cable id 1 begin 3.5 5 5 end 6 5 5 nseg 4
sel cable id 2 begin 6 5 5 end 10 5 5 nseg 6
sel delete link range cid 1
sel delete link range cid 3
sel delete link range cid 7
sel delete link range cid 8
sel delete link range cid 9
sel delete link range cid 6
sel link 2 target node tgt_num 6
plot
add sel geom link on
show

lmm

set log on
set logfile 2.1.log
model null range group exc1
;第一层锚索
;锚固段
sel cable id=1 begin 0.2825 0.05 1.20 end 0.3825 0.05 1.20 nseg=3
sel cable id=1 begin 0.2825 0.15 1.20 end 0.3825 0.15 1.20 nseg=3
sel cable id=1 begin 0.2825 0.25 1.20 end 0.3825 0.25 1.20 nseg=3
sel cable id=1 begin 0.2825 0.35 1.20 end 0.3825 0.35 1.20 nseg=3
sel cable id=1 begin 1.63 0.05 1.20 end 1.73 0.05 1.20 nseg=3
sel cable id=1 begin 1.63 0.15 1.20 end 1.73 0.15 1.20 nseg=3
sel cable id=1 begin 1.63 0.25 1.20 end 1.73 0.25 1.20 nseg=3
sel cable id=1 begin 1.63 0.35 1.20 end 1.73 0.35 1.20 nseg=3
;自由段
sel cable id=1 begin 0.3825 0.05 1.20 end 0.7125 0.05 1.20 nseg=7
sel cable id=1 begin 0.3825 0.15 1.20 end 0.7125 0.15 1.20 nseg=7
sel cable id=1 begin 0.3825 0.25 1.20 end 0.7125 0.25 1.20 nseg=7
sel cable id=1 begin 0.3825 0.35 1.20 end 0.7125 0.35 1.20 nseg=7
sel cable id=1 begin 1.30 0.05 1.20 end 1.63 0.05 1.20 nseg=7
sel cable id=1 begin 1.30 0.15 1.20 end 1.63 0.15 1.20 nseg=7
sel cable id=1 begin 1.30 0.25 1.20 end 1.63 0.25 1.20 nseg=7
sel cable id=1 begin 1.30 0.35 1.20 end 1.63 0.35 1.20 nseg=7
;sel delete link range cid 33
;sel delete link range cid 41
;sel delete link range cid 49
;sel delete link range cid 57
;sel delete link range cid 66
;sel delete link range cid 74
;sel delete link range cid 82
;sel delete link range cid 90

;sel link 2 target node tgt_num 33
;sel link attach xd rigid yd rigid zd rigid range &
;xr rigid yr rigid zr rigid cid 2

;sel link 6 target node tgt_num 41
;sel link attach xd rigid yd rigid zd rigid &
;xr rigid yr rigid zr rigid range cid 6

;sel link 10 target node tgt_num 49
;sel link attach xd rigid yd rigid zd rigid &
;xr rigid yr rigid zr rigid range cid 10
  
;sel link 14 target node tgt_num 57
;sel link attach xd rigid yd rigid zd rigid &
;xr rigid yr rigid zr rigid range cid 14

;sel link 17 target node tgt_num 66
;sel link attach xd rigid yd rigid zd rigid &
;xr rigid yr rigid zr rigid range cid 17

;sel link 21 target node tgt_num 74
;sel link attach xd rigid yd rigid zd rigid &
;xr rigid yr rigid zr rigid range cid 21

;sel link 25 target node tgt_num 82
;sel link attach xd rigid yd rigid zd rigid &
;xr rigid yr rigid zr rigid range cid 25

;sel link 29 target node tgt_num 90
;sel link attach xd rigid yd rigid zd rigid &
;xr rigid yr rigid zr rigid range cid 29

macro maosuo 'emod 3.6e9 xcarea 2.83e-5 ytens 3.82e7 gr_per 2.51e-2'
macro jiang 'gr_coh 2e5 gr_k 4.31e9 gr_fric 38'
macro ziyou 'gr_coh 0 gr_k 0 gr_fric 0'
macro maosuo1 'maosuo jiang' ;内锚固段锚索参数
macro maosuo2 'maosuo ziyou' ;自由段锚索参数

sel node apply xforce=80 range x 0.7115 0.7135
sel node apply xforce=-80 range x 1.299 1.301

sel cable id 1 prop maosuo1 range x 0.2825 0.3825
sel cable id 1 prop maosuo1 range x 1.63 1.73
sel cable id 1 prop maosuo2 range x 0.3825 0.7125 z 1.10 1.30
sel cable id 1 prop maosuo2 range x 1.30 1.63 z 1.10 1.30

由此例可以看出作者在结构单元部分关于link一部分的一些了解，原来作者是把自由段和锚固段分别定义为锚索1和锚索2，这样出现了很多问题，就是关于删除多余的铰结，并重新生成node－node铰结，并且铰结还要分别赋予各个方向的连接条件，并分别赋予铰结的属性，比如刚度，粘结力，摩擦角等。但是只要把锚固段和自由段定义为一根锚索，则使得这些问题问题化为乌有，问题得到大大的简化。

lmm

gen zone brick p0 0 0 0 p1 8 0 0 p2 0 8 0 p3 0 0 8 size 10 10 10
model mohr
prop bulk=6.7e8 shear=3.1e8 fric=41 coh=4.5e3 tens=3.5e3
fix x range x -.1 .1
fix y range y -.1 .1
fix y range y 9.9 10.1
fix z range z -.1 .1
apply szz=-1.8e5 range x 4 6 z 7.9 8.1
sel cable id 1 begin 4.5 1.6 4.5 end 8 1.6 4.5 nseg 5
sel cable id 1 begin 4.5 3.2 4.5 end 8 3.2 4.5 nseg 5
sel cable id 1 begin 4.5 4.8 4.5 end 8 4.8 4.5 nseg 5
sel cable id 1 begin 4.5 6.4 4.5 end 8 6.4 4.5 nseg 5

macro maosuo1 'xcarea=8.5e-3 emod=200e9 ytens=1e10 gr_k=0.0 gr_coh=0.0 gr_per=0.314 gr_fric=0.0'
macro maosuo2 'xcarea=8.5e-2 emod=200e9 ytens=1e10 gr_k=7e6 gr_coh=1e4 gr_per=0.314 gr_fric=20'
sel delete link range x 7.9 8.1
sel link net range x 7.9 8.1
sel link attach xd rig yd rig zd rig xr rig yr rig zr rig range x 7.9 8.1
sel cable id=1 prop maosuo2 range x 4.5 5.9
sel cable id=1 prop maosuo1 range x 5.9 8
sel node apply force 3.2e5 0 0 range x 7.9 8.1
solve

lmm

gen zone brick p0 0 0 0 p1 8 0 0 p2 0 8 0 p3 0 0 8 size 10 10 10
model mohr
prop bulk=6.7e8 shear=3.1e8 fric=41 coh=4.5e3 tens=3.5e3
fix x range x -.1 .1
fix y range y -.1 .1
fix y range y 9.9 10.1
fix z range z -.1 .1
apply szz=-1.8e5 range x 4 6 z 7.9 8.1
sel cable id 1 begin 4.5 1.6 4.5 end 8 1.6 4.5 nseg 5
sel cable id 1 begin 4.5 3.2 4.5 end 8 3.2 4.5 nseg 5
sel cable id 1 begin 4.5 4.8 4.5 end 8 4.8 4.5 nseg 5
sel cable id 1 begin 4.5 6.4 4.5 end 8 6.4 4.5 nseg 5

macro maosuo1 'xcarea=8.5e-3 emod=200e9 ytens=1e10 gr_k=0.0 gr_coh=0.0 gr_per=0.314 gr_fric=0.0'
macro maosuo2 'xcarea=8.5e-2 emod=200e9 ytens=1e10 gr_k=7e6 gr_coh=1e4 gr_per=0.314 gr_fric=20'
sel delete link range x 7.9 8.1
sel link net range x 7.9 8.1
sel link attach xd rig yd rig zd rig xr rig yr rig zr rig range x 7.9 8.1
sel cable id=1 prop maosuo2 range x 4.5 5.9
sel cable id=1 prop maosuo1 range x 5.9 8
sel node apply xforce=3.2e5 range x 7.9 8.1
solve

lmm

gen zon brick p0 0.725 0 1.1 p1 1.3 0 1.1 p2 0.725 0.4 1.1 &
p3 0.7 0 1.3 p4 1.3 0.4 1.1 p5 0.7 0.4 1.3 p6 1.3 0 1.3 &
p7 1.3 0.4 1.3 size 25 7 4
group exc1
gen zon brick p0 0.75 0 0.9 p1 1.30 0 0.90 p2 0.75 0.40 0.90 p3 0.725 0 1.10 &
p4 1.30 0.40 0.90 p5 0.725 0.40 1.10 p6 1.30 0 1.10 p7 1.30 0.40 1.10 &
size 25 7 4
group exc2 range group exc1 not
gen zon brick p0 0.775 0 0.70 p1 1.30 0 0.70 p2 0.775 0.40 0.70 &
p3 0.75 0 0.90 p4 1.30 0.40 0.70 p5 0.75 0.40 0.90 p6 1.30 0 0.90 p7 1.30 0.40 0.90 &
size 25 7 4
group exc3 range group exc1 not group exc2 not
gen zon brick p0 0.80 0 0.50 p1 1.30 0 0.50 p2 0.80 0.40 0.50 p3 0.775 0 0.70 &
p4 1.30 0.40 0.50 p5 0.775 0.40 0.70 p6 1.30 0 0.70 p7 1.30 0.40 0.70 size 25 7 4
group exc4 range group exc1 not group exc2 not group exc3 not

gen zone brick p0 0 0 0.5 p1 0.3 0 0.5 p2 0 0.4 0.5 p3 0 0 1.3 size 9 7 16 ratio 0.87 1 1
gen zone brick p0 0.3 0 0.5 p1 0.8 0 0.5 p2 0.3 0.4 0.5 p3 0.3 0 1.3 p4 0.8 0.4 0.5 &
p5 0.3 0.4 1.3 p6 0.7 0 1.3 p7 0.7 0.4 1.3 size 25 7 16
gen zone brick p0 1.3 0 0.5 p1 1.7 0 0.5 p2 1.3 0.4 0.5 p3 1.3 0 1.3 size 20 7 16
gen zone brick p0 1.7 0 0.5 p1 2 0 0.5 p2 1.7 0.4 0.5 p3 1.7 0 1.3 size 9 7 16 ratio 1.15 1 1

gen zone brick p0 0 0 0 p1 0.3 0 0 p2 0 0.4 0 p3 0 0 0.5 size 9 7 7 ratio 0.87 1 0.91
gen zone brick p0 0.3 0 0 p1 0.8 0 0 p2 0.3 0.4 0 p3 0.3 0 0.5 size 25 7 7 ratio 1 1 0.91
gen zone brick p0 0.8 0 0 p1 1.7 0 0 p2 0.8 0.4 0 p3 0.8 0 0.5 size 45 7 7 ratio 1 1 0.91
gen zone brick p0 1.7 0 0 p1 2 0 0 p2 1.7 0.4 0 p3 1.7 0 0.5 size 9 7 7 ratio 1.15 1 0.91
model null range group exc1
model null range group exc2
model null range group exc3
model null range group exc4

;sel shell id=1 range x 0.7 0.8 z 0.5 1.3 direction 0.99228 0 0.12403
;sel shell id=1 range direction 0.99228 0 0.12403 x 0.7 0.8 ;和上句效果一样
sel shell id=1 range plane normal 0.99228 0 0.12403 origin     dist     
plot show

lmm

gen zone brick p0 0 0 0 p1 8 0 0 p2 0 8 0 p3 0 0 8 size 10 10 10
model mohr
prop bulk=6.7e9 shear=3.1e9 fric=41 coh=4.5e3 tens=3.5e3
fix x range x -.1 .1
fix x range x 9.9 10.1
fix y range y -.1 .1
fix y range y 9.9 10.1
fix z range z -.1 .1
apply szz=-1.8e5 range x 4 6 z 7.9 8.1
sel cable id 1 begin 4.5 1.6 4.5 end 8 1.6 4.5 nseg 5
sel cable id 1 begin 4.5 3.2 4.5 end 8 3.2 4.5 nseg 5
sel cable id 1 begin 4.5 4.8 4.5 end 8 4.8 4.5 nseg 5
sel cable id 1 begin 4.5 6.4 4.5 end 8 6.4 4.5 nseg 5
;sel node apply xforce=3.2e5 range cid=2;第一根锚索的自由端
;sel node apply xforce=3.2e5 range cid=10；第二根锚索的锚固段近端
sel node apply xforce=3.2e5 range x 7.9 8.1；四根锚索的自由端
macro maosuo1 'xcarea=8.5e-3 emod=200e9 ytens=1e10 gr_k=0.0 gr_coh=0.0 gr_per=0.314 gr_fric=0.0'
macro maosuo2 'xcarea=8.5e-2 emod=200e9 ytens=1e10 gr_k=7e6 gr_coh=1e4 gr_per=0.314 gr_fric=20'
sel cable id=1 prop maosuo2 range x 4.5 5.9
sel cable id=1 prop maosuo1 range x 5.9 8
step 100
save 1.sav



gen zone brick p0 0 0 0 p1 8 0 0 p2 0 8 0 p3 0 0 8 size 10 10 10
model mohr
prop bulk=6.7e8 shear=3.1e8 fric=41 coh=4.5e3 tens=3.5e3
fix x range x -.1 .1
fix x range x 9.9 10.1
fix y range y -.1 .1
fix y range y 9.9 10.1
fix z range z -.1 .1
apply szz=-1.8e5 range x 4 6 z 7.9 8.1
sel cable id 1 begin 4.5 1.6 4.5 end 8 1.6 4.5 nseg 5
sel cable id 1 begin 4.5 3.2 4.5 end 8 3.2 4.5 nseg 5
sel cable id 1 begin 4.5 4.8 4.5 end 8 4.8 4.5 nseg 5
sel cable id 1 begin 4.5 6.4 4.5 end 8 6.4 4.5 nseg 5
sel cable id 1 pretension 3.2e5；此处是把预应力施加于整个锚索上，包括锚固段和自由段
macro maosuo1 'xcarea=8.5e-3 emod=200e9 ytens=1e10 gr_k=0.0 gr_coh=0.0 gr_per=0.314 gr_fric=0.0'
macro maosuo2 'xcarea=8.5e-2 emod=200e9 ytens=1e10 gr_k=7e6 gr_coh=1e4 gr_per=0.314 gr_fric=20'
sel cable id=1 prop maosuo2 range x 4.5 5.9
sel cable id=1 prop maosuo1 range x 5.9 8
step 100
save 1.sav
此处为两种关于预应力的施加方法
其中，通过节点施加预应力的方法必须通过计算，其轴力、剪应力等才会在图形中显示出来；而通过直接施加预应力的方法在没有经过计算的情况下也可以在图形中看到轴力、剪应力等。其原因可能是通过sel node apply，即节点施加预应力的方法相当于给定模型一个力边界条件，边界条件只有通过计算才能看到锚索体的轴力、剪应力等，即锚索体的轴力、剪应力是这个边界条件所产生的效果；而通过sel cable pretension命令施加的预应力则本身就是锚索单元的自身属性，可以在图形中看到其施加的效果，即相当于我们可以在图形中看到所施加的边界条件。

lmm

上面的命令流是我在2004年下半年学习FLAC时做的一些算例，因为当时没有任何关于锚索方面可以参考的命令流文件，包括这个论坛上也没有，所以当时对锚索的摸索花了很多的精力，大家可以从我上面的一些算例上看出当时对一些命令的运用以及该命令所产生的效果都是自己通过简单的算例来试验得出来的，所以我想给大家提供一个经验：如果你对哪一条命令不解的话，你可以针对这个命令做一个小的算例，看看命令产生的效果是不是和想象的一样，这样你对该命令就会有更深刻的理解

gongsunde

用flac做基坑模拟，不错！

276774950

不错，谢谢楼主

江苏萧肖

谢谢您的上传哦  我的QQ是389439159  还有东西要向你请教呢？

csu_lijun

很好的帖子，建议置顶

csu_lijun

macro maosuo2 'xcarea=8.5e-2 emod=200e9 ytens=1e10 gr_k=7e6 gr_coh=1e4 gr_per=0.314 gr_fric=20'
xcarea=8.5e-2
gr_per=0.314
这两个如何取值，是锚索直径还是注浆体直径或者分开取值
望楼主解释

lmm

xcarea是锚索的断面积
gr_per是锚索加灌浆体的总周长

arnoldray

什么时候俺也能像楼主如此厉害啊

secondlone

下来学习学习先

nndd521

谢谢分享

hxbhwj

谢谢分享

hnchinabbs

谢谢楼主共享

cumt_whs

感谢楼主共享，麻烦楼主能否把命令流文件（txt格式的）共享一下，非常谢谢。

zdpjll

鼓励，支持+学习