flac3d模拟锚杆三种方法

n
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
m m
pro bul 2.2e9 she 1.3e9 fric 30 coh 1.3e6 ten 1.5e5
ini den 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

;方式一：通过删除，建立link模拟
sel cable id 1 beg 0 0 0 end 0 29 0 nse 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 16
sel delete link range cid 1 ;或id＝1，link编号从锚索头到尾依次编号1……18，无11（11是连接点处）
sel link id 100 1 target zone ;以前link（1）改名为100，id可以用cid
sel link attach xdir=rigid ydir=rigid xrdir=rigid zrdir=rigid range id 100 ;id可以用cid
sel cable id 1 pretension 60e3 range cid 1,10

pl set rot 90 0 0
pl sel cabl force
pl add sel geom cid on
solve


;方式二：通过赋三段属性模拟，即将锚索的端头、自由段、锚固段赋不同的属性模拟预应力锚索。
sel cable id 1 beg 0 0 0 end 0 29 0 nse 10
sel cable id 1 beg 0 29 0 end 0 35 0 nseg 6
sel cable prop emod 2e10 ytension 310e3 xcarea 0.0004906 &
gr_coh 1 gr_k 1 gr_per 0.0785 range cid 2 10
sel cable id 1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 &
gr_coh 10e5 gr_k 2e7 range cid 11 16
sel cable id 1 prop emod 2e10 ytension 310e3 xcarea 0.0004906 &
gr_coh 10e8 gr_k 2e10 range cid 1 1
sel cable id 1 pretension 60e3 range cid 1,10

pl set rot 90 0 0
pl sel cabl force
pl add sel geom cid on
solve


;方式三：借助别的结构单元（如liner单元）来模拟托盘，删除锚索端头的link，然后建立新的link，&
;新link的target为liner上的node
sel cable id 1 beg 0 0 0 end 0 29 0 nse 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 16
sel liner range y -.1 .1 x -1 1 z -1 1
sel liner prop iso=(25e9,0.15) thick 0.1 ;concrete
sel liner prop cs_nk=8e8 cs_sk=8e8 &
cs_ncut 0.0 cs_scoh 0.0 cs_scohres 0.0 cs_sfric 0.0
sel delete link range id 1
sel link id 100 1 target node tgt_num 18
sel link attach xdir=rigid ydir=rigid xrdir=rigid zrdir=rigid range id 100
sel cable id 1 pretension 60e3 range cid 1,10

pl set rot 90 0 0
pl sel cabl force
pl add sel geom cid on
step 2000
save sy.sav


Flac3d输出：

1 Edit copy to clipboard 粘贴到word

2 set log on

Set logfile A.txt

Print sel cab

force

Set logfile B.txt

Pri sel cab str

Set log off

3 视频记录

Config dyn ；命令开头

Set movie avi step 1000 file name.avi

Movie start

Solve

Movie finish

4 solve

hist write 1 vs step file aaaaaa.xls
his write 1 vs step begin 1 end 2000 file aaaaaaa.his


本人是FLAC3D的初学者，现遇到一个工程，就是对一个公路边坡衡重式挡墙进行加固，采用的是预应力锚索抗滑桩方式加固该挡墙（因为是做这方面的优化设计研究，所以采用的预应力锚索抗滑桩）。
为了后处理方便，桩采用pile单元来模拟，预应力锚索采用cable单元。
现遇到以下几个问题：
1.pile和cable单元在连接处是否共接点
2.施加的预应力是否作用到了桩上
' a7 M) @" h: }' p7 p; } 3.pile单元是一条线，由于挡墙墙背倾斜很小，建模时我忽略了这个倾斜，认为墙背垂直，这样桩和挡墙间素填的混凝土也就忽略了。加pile单元的时候，pile的线就和挡墙的墙背所在的面重合了，不知道这样是否合理。0 r6 U0 O9 u9 c y& |1 |4 M3 B8 H
以下是关于cable的建立以及处理预应力的情况的部分命令流，我按先加桩和先加锚索分别模拟了：
1.先加锚索：
sel cable id=1 begin=(67.8398 45.8598 3) end=(95.0291 33.1813 3) nseg=15
sel cable id=1 begin=(67.8398 45.8598 9) end=(95.0291 33.1813 9) nseg=15
sel cable id =1 pro xcarea=1.77e-4 emod=195e9 ytens=1.32e9 &
gr_k=1 gr_coh=1 gr_per=0.3454 range cid=1,10 any cid=16,25 any
sel cable id =1 pro xcarea=1.77e-4 emod=195e9 ytens=1.32e9 &
gr_fric=60 gr_k=11.5e6 gr_coh=10e5 gr_per=0.3454 &
range cid=11,15 any cid=26,30 any
sel delete link range id 10
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 delete link range id 17
sel link id=200 17 target zone
sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id 200
sel cable id=1 pretension=0.1e6 range cid=1,10 any cid=16,25 any
sel cable pro slide=on
plot add sel cable force yellow red
plot add sel cable grout slip
sel pile id=1 begin=(67.8398 46.3598 3) end=(67.8398 45.8598 3)（为了考虑桩与锚索共结点，特别在桩和锚索的连接点处将桩分两次建立）
sel pile id=1 begin=(67.8398 45.8598 3) end=(67.8398 28.5598 3) nseg=10
sel pile id=1 begin=(67.8398 46.3598 9) end=(67.8398 45.8598 9)
sel pile id=1 begin=(67.8398 45.8598 9) end=(67.8398 28.5598 9) nseg=10
sel pile id=1 pro emod=30e9 nu=0.3 xcarea=3.6 xciy=1.728 xciz=0.675 xcj=2.403 &
per=7.8 cs_sk=1.3e11 cs_scoh=1.0e10 cs_sfric=0 &
cs_nk=1.3e9 cs_ncoh=1.0e04 cs_nfric=0 cs_ngap=off
sel pile pro slide=on