用flac2d做桩动测的接触面问题

Dchenxue

我用flac2d做桩的动测模拟，但是，我无论如何改变接触面的参数，实体桩内的应力和速度曲线都无明显变化，烦请各位大侠指导一下，找出问题所在。
config dynamic
grid 10,20
model elastic
gen (0.0,0.0) (0.0,10.0) (3.0,10.0) (3.0,0.0) ratio 1.0,1.0 i 1 11 j 1 21
model null i=4 j=1 20
model null i=7 j=1 20
ini x add -0.3 y add 0.0 nmregion 5 15
ini x add -0.6 y add 0.0 nmregion 8 15
interface 1 aside from 5,21 to 5,1 bside from 4,21 to 4,1
interface 1 unglued kn=8.0E9 ks=8.0E9 cohesion=7500.0 dilation=0.0 friction=19.0 tbond=0.0 bslip=on
interface 2 aside from 7,21 to 7,1 bside from 8,21 to 8,1
interface 2 unglued kn=8.0E9 ks=8.0E9 cohesion=7500.0 dilation=0.0 friction=19.0 tbond=0.0 bslip=on
group 'pile' region 5 15 notnull
model elastic group 'pile'
prop density=2500.0 bulk=1.39E10 shear=1.04E10 group 'pile'
group 'claya'  region 3 5 notnull
model mohr group 'claya'
prop density=1670 bulk=2.621E7 shear=1.0E6 cohesion=15000.0 friction=15.50 dilation=0.0 tension=0.0 group 'claya'
group 'claya'  region 8 5 notnull
model mohr group 'claya'
prop density=1670 bulk=2.621E7 shear=1.0E6 cohesion=15000.0 friction=15.50 dilation=0.0 tension=0.0 group 'claya'
fix x i=1 j= 1 21
fix x i=11 j= 1 21
fix  y i=1 4  j=1
fix  y i=8 11  j=1
set gravity=10
history 1 yvel i=6, j=21
history 2 yvel i=6, j=20
history 3 yvel i=6, j=1
history 4 syy i=6, j=20
history 5 syy i=6, j=19
history 6 syy i=6, j=1
set dyn off
solve



initial xdisp 0 ydisp 0
initial xvel 0 yvel 0
set dyn on
set dy_dampi local 0.15
set =large
def wave
wave=18.3785e6*exp(-631*dytime)
end
;apply yquiet  from 2,31 to 4,31
apply syy -2 hist wave from 5,21 to 7,21
hist nstep 1
hist wave
hist dytime
hist ydis i=6 j=1
hist ydis i=6 j=20
solve dytime 0.015

happywmj

我也在用flac3d做这个问题，也是困难重重。   是不是地震时间太短了，只有0.015s，另外我觉得应力时程上是不是还有点问题。

yinwen0809

是不是直接用桩单元来算会好些

zhengliminday

yinwen0809 发表于 2013-10-24 15:08
是不是直接用桩单元来算会好些

那个，大哥，我问下那个桩单元在哪里找到呀，或者说桩是怎么建的呀，我是个新手，对这些不太懂，看您能否给个模型什么的，谢谢啦