用flac3d模拟锚杆杆体拉伸试验

jhtdasa

new
sel cable id 1 beg 0 0 0 end 0.05 0 0 nseg 1
sel cable id 1 beg 0.05 0 0 end 2.05 0 0 nseg 40
sel cable id 1 beg 2.05 0 0 end 2.1 0 0 nseg 1
;pl sel geo cid on node on scale 0.02
sel cable id 1 prop emod 20e10 ytension 266e3 xcarea 0.00038
;;;;;;;边界条件;;;;;;;;;;;;;
sel node fix x yr zr range id 3 any id 43 any
sel node ini xvel -1e-5 range id 1 any id 2 any
;;;;;;;;;;测点;;;;;;;;;;
hist sel node xdisp id 1
hist sel cablesel force cid 2
pl his 2 v -1
step 15000
pl sel cable force

实验时，锚杆右端固定，左端以恒定速度1e-5 m/step的速度向左拉伸，运行到15 000步时停止计算，其中位移采用最左端节点位移监测值，轴力采用CID2单元的轴力问题：一运行，锚杆轴力就降到0，请问程序有什么不合理？哪里出问题了

dongxx0221

sel node APPLY xvel -1e-5 range id 1 any id 2 any

jhtdasa

dongxx0221 发表于 2016-3-15 22:04
sel node APPLY xvel -1e-5 range id 1 any id 2 any

您好，我试了一下，apply只能加force或moment，对于xvel不识别

dongxx0221

jhtdasa 发表于 2016-3-16 09:52
您好，我试了一下，apply只能加force或moment，对于xvel不识别

哦  ini绝对是不合适的。可以在两端做实体，两端nod与zone固接，给实体apply xvel. 类似于用锚头锁住拉。

jhtdasa

dongxx0221 发表于 2016-3-16 12:09
哦  ini绝对是不合适的。可以在两端做实体，两端nod与zone固接，给实体apply xvel. 类似于用锚头锁住拉。
...

您好，您说的两端做实体实现了锚杆拉伸，非常感谢!
不知您对sc_length了解吗？
def xl
     whilestepping   
     t=2
     weizhi=s_find(t)              ;2为指定cable单元的CID
     xl=sc_length(weizhi)  ;获取长度         
end
his xl

我想要监测cid2这段在计算过程中长度的变化，为什么我的xl没有循环记录计算过程中的值，一直是0.05