- 积分
- 0
- 注册时间
- 2011-9-30
- 仿真币
-
- 最后登录
- 1970-1-1
|
我用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 |
|