求助！（急）关于pile单元和cable的连接和施加预应力问题

hy1984

本人是FLAC3D的初学者，现遇到一个工程，就是对一个公路边坡衡重式挡墙进行加固，采用的是预应力锚索抗滑桩方式加固该挡墙（因为是做这方面的优化设计研究，所以采用的预应力锚索抗滑桩）。
   为了后处理方便，桩采用pile单元来模拟，预应力锚索采用cable单元。
   现遇到以下几个问题：
   1.pile和cable单元在连接处是否共接点
   2.施加的预应力是否作用到了桩上
   3.pile单元是一条线，由于挡墙墙背倾斜很小，建模时我忽略了这个倾斜，认为墙背垂直，这样桩和挡墙间素填的混凝土也就忽略了。加pile单元的时候，pile的线就和挡墙的墙背所在的面重合了，不知道这样是否合理。
   以下是关于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 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 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

2。先加桩
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

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 dens=3000 &
gr_k=1 gr_coh=1 gr_per=0.3454 range cid=23,32 any cid=38,47 any
sel cable id =1 pro xcarea=1.77e-4 emod=195e9 ytens=1.32e9 dens=3000 &
gr_fric=60 gr_k=11.5e6 gr_coh=10e5 gr_per=0.3454 &
range cid=33,37 any cid=48,52 any

sel delete link range id 27
sel delete link range id 3（有点不明白为什么这种方式会在桩与锚索的连接处会出现两个link，恳请高手和前辈给予解答）
sel link id=100 2 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 43
sel delete link range id 16
sel link id=200 14 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=23,32 any cid=38,47 any
sel cable pro slide=on
   
   这两种方式处理时都没有出现错误，也能算的收敛，但计算的结果有些令人不满意（桩的弯距太小，与手算相差太多（相差两个数量级），而剪力较大，位移还算可以，手算是在不考虑挡墙作用下按现有设计方法计算桩的内力），我怀疑是不是共接点和预应力的施加出现问题，请高手和前辈们帮我解答下，采用这两种方式是否达到了以下的目的：
  1.pile和cable单元在桩与锚梭连接处（例如在点67.8398  45.8598  3处）共接点
   2.施加的预应力是作用到了桩上
   3.pile单元是一条线，由于挡墙墙背倾斜很小，建模时我忽略了这个倾斜，认为墙背垂直，这样桩和挡墙间素填的混凝土也就忽略了。加pile单元的时候，pile的线就和挡墙的墙背所在的面重合了，不知道这样是否合理。

在这个问题上我困惑了很久，望前辈们尽快给予答复，十分感谢！！或者能给出关于pile单元和cable单元共节点的处理并施加预应力的完整命令流。

hy1984

本人邮箱为huang158374476@yahoo.com.cn
希望能与前辈们多多交流
再请教下3D边坡中滑面是怎样考虑和模拟的

hy1984

请各位前辈和高手们给予指导啊

weekly007

同样想知道答案

csu_lijun

支持一下
希望有高手来解答

csu_lijun

请教楼主
读了读你的程序，怎么锚索好像是分段的，锚固段和自由段是间隔的，你的有四段。
这个是怎么个考虑

streamdream

希望有高手来解答

xiaoxue730

不明白pile为什么要建两次啊？

小珍鹏

同样想知道解答，楼主问题解答了吗？可否共享下。