请教关于lamb解的问题

imwnab

我用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 编辑 ]

chensha

解决问题没有？

chensha

上传一个报告。

古襄月

都是老问题