陈育民书中模拟预应力锚杆方式一：通过删除，建立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
sel link id 100 1 target zone
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
solve
其中
sel delete link range cid 1
sel link id 100 1 target zone
sel link attach xdir=rigid ydir=rigid xrdir=rigid zrdir=rigid range id 100
为删除链接1，建立链接100，并把新链接100约束以模拟预应力锚墩。个人觉得这个方法有点繁琐，可以改为
sel link attach xdir=rigid ydir=rigid xrdir=rigid zrdir=rigid range id 1
这样程序可以简化为
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 link attach xdir=rigid ydir=rigid xrdir=rigid zrdir=rigid range id 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
这样程序会简化一些，经过验证，运算结果与书中例子相同。
