边界设置了约束，但是加桩后，边界处桩的位移和里面的...

龙一乜

边界设置了约束，但是加桩后，边界处桩的位移和里面的桩位移不一样。[img=0,1]file:///C:\Users\h\AppData\Roaming\Tencent\Users\344557819\QQ\WinTemp\RichOle\_E%J9$C98TT[`N[YZ}T]14G.png[/img]计算好长时间也不收敛。哪位大神帮忙看看是哪里出了错。。。。


new
gen zone b p0 0 0 0 p1 30 0 0 p2 0 15 0 p3 0 0 26 p4 30 15 0 size 30 15 26 ratio 0.95 1.00 0.95
gen zone b p0 0 0 26 p1 30 0 26 p2 0 15 26 p3 0 0 34 p4 30 15 26 size 30 15 16 ratio 0.95 1.00 1
gen zone b p0 30 0 0 p1 70 0 0 p2 30 15 0 p3 30 0 26 p4 70 15 0 size 40 15 26 ratio 1.05 1.00 0.95
gen zone b p0 30 0 26 p1 70 0 26 p2 30 15 26 p3 30 0 34 p4 70 15 26 size 40 15 16 ratio 1.05 1.00 1
gen zone b p0 30 0 34 p1 70 0 34 p2 30 15 34 p3 40 0 44 p4 70 15 34 p5 40 15 44 p6 70 0 44 p7 70 15 44 size 40 15 10 ratio 1.05 1.00 1.05

;pause
gen merge 1e-6
attach face
group tt range z 0 44

model mohr range group tt
prop bulk 21e6 shear 9.62e6 fric 25 coh 1e10 dilation 0 tension 1e10 poisson 0.3 range group tt
ini dens 1700 range group tt

fix x y z range z -0.1 0.1
fix y range y -0.1 0.1
fix y range y 14.9 15.1
fix x range x -0.1 0.1
fix x range x 69.9 70.1

set gravity 0 0 -9.8
hist n 5
hist unbal


hist gp zdisp 30 0 34
hist gp zdisp 30 10 34
hist gp zdisp 40 0 44
hist gp zdisp 40 10 44

solve

prop bulk 21e6 shear 9.62e6 fric 25 coh 20000 poisson 0.3 range group tt

solve
save tt1-1.sav

ini xdisp 0 ydisp 0 zdisp 0
ini xvel 0 yvel 0 zvel 0

sel pile id 1 begin 30 0 8 end 30 0 34 nseg=52
sel pile id 2 begin 30 5 8 end 30 5 34 nseg=52
sel pile id 3 begin 30 10 8 end 30 10 34 nseg=52
sel pile id 4 begin 30 15 8 end 30 15 34 nseg=52

sel pile pro d 2500 emod 3.0e10 nu 0.25 xcarea 6 xciy 4.5 xciz 2 xcj 6.5 per 10 cs_scoh 30000 cs_sfric 25 cs_sk 4.3e9 cs_ncoh 30000 cs_nfric 25 cs_nk 4.3e9 cs_ng=on slide=on

prop bulk 21e6 shear 9.62e6 fric 25 coh 20000 poisson 0.3 range group tt


hist reset
hist n 5
hist unbal
hist sel pilesel moment my 30 0 34
HIST SEL NODE xDISP ID=1
HIST SEL PILESEL FORCE FX CID=1

solve
save tt0-2.sav