蠕变例子

new
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
mo mohr
pro bulk 44.4e9 she 14.8e9 fric 26.6 coh 0.3e6
ini dens 2200
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
PLOT SURF
hist id=100 unbal
;
sel cable id=1 beg 0, 0, 0 end 0 ,29, 0 nseg 10
sel cable id=1 beg 0,29,0 end 0,35,0 nseg 6
sel cable id=1 prop emod 3e10 ytension 310e4 xcarea 0.0177 gr_c 1 gr_k 1 gr_per 0.471 range cid 1,10
sel cable id=1 prop emod 3e10 ytension 310e4 xcarea 0.0177 gr_c 10e4 gr_k 2e6 gr_per 0.471 range cid 11,17
sel delete link range id 1 ;1为link编号
sel link id=100 1 target zone ;100为新的link编号 1为节点号 即在节点1上建立编号为100的node-zone连接
sel link attach xdir=rigid ydir=rigid zdir=rigid xrdir=rigid yrdir=rigid zrdir=rigid range id 100
sel cable id=1 pretension 60e3 range cid 1,10 ;cid为结构构件编号
solve
plot set plane ori 0 5 0 norm 0 0 1
plot cont syy plane ou on
plot add axe
plot set rot 90 0 0
plot set mag 2
save 1.sav
;考虑蠕变10天
new
rest 1.sav
config creep
set creep mindt=1e-5 maxdt=0.1
mo cvisc
pro bulk 1.1e10 she 0.83e10 fric 30 coh 5e4 ten 2e6 &
kshear 14.8e9 kvis 45e9 mshear 1.48e9 mvis 600e9
hist id=1000 unbal
;
set creep on
set creep dt auto on
solve age 10
plot set plane ori 0 5 0 norm 0 0 1
plot cont syy plane ou on
plot add axe
plot set rot 90 0 0
plot set mag 2
save 3.sav