手册上的例子!
new
title
axially loaded pile
gen zone brick size 15 15 15 edge=11
mode elastic
prop bulk 5e9 shear 1e9 density=2000
fix z range z -0.01 0.01
fix x range x -0.01 0.01
fix x range x 10.99 11.01
fix y rang y -0.01 0.01
fix y range y 10.99 11.01
ini szz -2.2e5 grad 0 0 2e4
ini sxx -2.2e5 grad 0 0 2e4
ini syy -2.2e5 grad 0 0 2e4
set gravity 0 0 -10
solve
sel pile id=1 begin 5.5 5.5 12.0 end 5.5 5.5 4.0 nseg=32
sel pile prop emod 8e10 nu 0.3 xcarea=0.754 xcj=0 xciy=0 xciz=0 &
per=3.14 cs_sk=1.3e11 cs_scoh=0 cs_sfric=10 cs_nk=1.3e11 cs_ncoh=0 cs_nfric=0 cs_ngap=off
sel delete link range id 29
sel link id=100 2 target zone
sel link attach ydir=free zdir=free range id 100
sel link attach xrdir=free yrdir=free zrdir=free range id 100
sel link attach xdir=nydeform range id 100
sel link constit nydeform 1 area=1 k=5.4e11 ycomp=2.22e5 range id 100
sel node fix x range id=1
sel node ini zvel=-0.5e-8 range id=1
hist id=1 unbal
hist id=2 sel node zdisp id=1
hist id=3 sel pilesel force fx cid=1
plot create history
plot add sel pile coupling stress shear max 5e4
plot add hist -3 vs -2 xmax 12e-5 ymax 10e5
plot set magn 0.8
plot set plane origin 0 5.5 0 normal 0 -1 0
plot add grid plane
plot add axes
plot show
sel set damp combined
cycle 1
cycle 20000
save plie-ax-pile2.sav
请问 sel node ini zvel=-0.5e-8 range id=1中的速度是怎么确定的?我把速度改成-0.5e-5结果大不一样!