- 积分
- 0
- 注册时间
- 2006-1-19
- 仿真币
-
- 最后登录
- 1970-1-1
|
我用Flac2D的轴对称模型模拟弹性半空间地基在竖向集中简谐力作用下的lamb经典解,场地取20m长,5m高,E=4.25e7,K=2.83e7,G=1.7e7,密度=1700,泊松比=0.25,剪切波速100m/s。简谐力F=1.0×sin(2*pi*freq*dytime),freq=20Hz,作用在中心轴上,监测地表面水平距离0.1m处的竖向位移。我算出的竖向位移峰值是1e-8,但是经典解析解的结果是7.6e-8。我找不出错误出在哪里,我师兄用ANSYS计算,同样用轴对称模型,同样边界条件,结果就很接近经典解析解,而且20m远处边界改为自由端的化结果基本相同,而我的命令流要是进行同样改动就会差别很大。请哪位高手能否给予指正。以下是我写的命令流:
new
config dynamic axi
set dyn off
grid 48 50
gen 0,0 0,5 1,5 1,0 i=1,11 j=1 51
gen 1,0 1,5 20,5 20,0 i=11,49 j=1 51
model elas
property density=1700 bulk=2.83e7 shear=1.7e7
ini syy -85000 var 0 85000
ini sxx -28050 var 0 28050
ini szz -28050 var 0 28050
fix x y j=1
fix x i=49
set grav=10
hist unbal
solve
save D:\FLAC\practise\dyn-p\axi\8-24\a1.sav
set legend on
set plot windows color 0.0 0.0 1.0 1.0
set plot jpg color
set output D:\FLAC\practise\dyn-p\axi\8-24\grid-a.jpg
title
grid
plot pen grid blue fix
set output D:\FLAC\practise\dyn-p\axi\8-24\syy-a.jpg
title
syy
sclin 1 (1,0) (1,5)
plot pen syy blue int=8500 bou brown
set output D:\FLAC\practise\dyn-p\axi\8-24\sxx-a.jpg
title
sxx
sclin 1 (1,0) (1,5)
plot pen sxx blue int=2805 bou brown
set output D:\FLAC\practise\dyn-p\axi\8-24\szz-a.jpg
title
szz
sclin 1 (1,0) (1,5)
plot pen szz blue int=2805 bou brown
;
set dyn on
set large
def wave
wave = sin(2.0*pi*freq*dytime)
end
set freq=20.0
apply yforce 1.0 hist wave i=1 j=51
set dytime 0
hist dytime ;记录2
hist wave ;记录3
hist xdis i=2 j=51 ;记录4
hist ydis i=2 j=51 ;记录5
hist xdis i=1 j=51 ;记录6
hist ydis i=1 j=51 ;记录7
solve dytime=0.2
save D:\FLAC\practise\dyn-p\axi\8-24\a2.sav
hist write 5 vs 2 table 5
print table 5
set output D:\FLAC\practise\dyn-p\axi\8-24\ydisp-a.jpg
title
20Hz-ydisp
plot pen hist 5 vs 2
set output D:\FLAC\practise\dyn-p\axi\8-24\ydisp1-a.jpg
title
20Hz-ydisp
plot pen hist 7 vs 2
[ 本帖最后由 imwnab 于 2007-8-30 16:46 编辑 ] |
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
|