找回密码
 注册
Simdroid-非首页
查看: 87|回复: 2

[动力分析] 关于“参数数组索引超出范围”的错误提示如何解决,命.....

[复制链接]
发表于 2016-7-21 16:08:39 | 显示全部楼层 |阅读模式 来自 浙江舟山
new
config fluid dynamic
set dyn off
gen zone brick size 2 10 2 p0 (0,0,0)  p1 (.3,0 ,0) p2 (0,0,-10) p3 (0,.3,0) group soil1

gen zone brick size 2 10 5 p0 (0 ,0.3,0) p1 (0.3,0.3,0) p2 (0,0.3,-10) p3 (0,6,0) p4 (0.3,0.3,-10) &
                           p5 (0,6,-10) p6 (2, 6, 0)  p7 (2 ,6 ,-10) group soil2

gen zone brick size 2 10 5 p0 (0.3,0.3,0) p1 (0.3,0,0) p2 (0.3,0.3,-10) p3 (6,2,0) p4 (0.3,0,-10) &
                           p5 (6,2,-10) p6 (6, 0, 0)  p7 (6 ,0 ,-10) group soil2
                           
gen zone brick size 5 10 5 p0 (0.3,0.3,0) p1 (6,2,0) p2 (0.3,0.3,-10) p3 (2,6,0) p4 (6,2,-10) &
                           p5 (2,6,-10) p6 (8, 8, 0)  p7 (8 ,8 , -10) group soil2
                           
         
gen zone brick size 2 10 2 p0 (0,6,0) p1 (2,6,0) p2 (0,6,-10) p3 (0,10,0) p4 (2,6,-10) &
                           p5 (0,10,-10) p6 (2, 10, 0)  p7 (2 ,10 , -10) group soil2                  
                                                                                   
gen zone brick size 2 10 6 p0 (6,2,0) p1 (6,0,0) p2 (6,2,-10) p3 (20,2,0) p4 (6,0,-10) &
                           p5 (20,2,-10) p6 (20, 0, 0)  p7 (20 ,0 , -10) group soil2     
                                      
                    
gen zone brick size 5 10 2 p0 (2,6,0) p1 (8,8,0) p2 (2,6,-10) p3 (2,10,0) p4 (8,8,-10) &
                           p5 (2,10,-10) p6 (8, 10, 0)  p7 (8,10 , -10) group soil2              
                  
gen zone brick size 5 10 6 p0 (8,8,0) p1 (6,2,0) p2 (8,8,-10) p3 (20,8,0) p4 (6,2,-10) &
                           p5 (20,8,-10) p6 (20, 2, 0)  p7 (20 ,2 , -10) group soil2      
                           
                          
gen zone brick size 6 10 2 p0 (8,8,0) p1 (20,8,0) p2 (8,8,-10) p3 (8,10,0) p4 (20,8,-10) &
                           p5 (8,10,-10) p6 (20, 10, 0)  p7 (20 ,10 , -10) group soil2         
                                      
gen zone reflec dip 90 dd 0
gen zone reflec dip 90 dd 90   


;边界条件
fix x y z ran z -9.9 -10.1
fix x ran x -19.9 -20.1
fix x ran x  19.9  20.1
fix y ran y  -9.9 -10.1
fix y ran y   9.9  10.1

; 属性
model mo range z 0 -10.1
model fl_iso ran z 0 -10.1
ini fdens 1000
set grav 10.0
ini pp 0 grad 0 0 -10e3  ran z 0 -10.1
hist unbal

set mech on fluid off
ini fmod 0

prop dens 1040 por 0.66 bu 14.71e6 sh 5.640e6  f 25 c 10e3  perm 1e-13 range z  0 -1      
prop dens 1440 por 0.46 bu 2.941e7 sh 1.128e7  f 30 c  0    perm 1e-8  range z -1 -9      
prop dens 1100 por 0.60 bu 14.71e6 sh 5.640e6  f 25 c 15e3  perm 1e-13 range z -9 -10.1   

ini szz 0 grad 0 0 17e3   ran z 0 -1  ;soila     
ini syy 0 grad 0 0 13.5e3 ran z 0 -1            
ini sxx 0 grad 0 0 13.5e3 ran z 0 -1            
                                                
ini szz 2e3 grad 0 0 19e3   ran z -1 -9  ;soilb   
ini syy 1e3 grad 0 0 14.5e3 ran z -1 -9            
ini sxx 1e3 grad 0 0 14.5e3 ran z -1 -9            
                                                
ini szz -16e3 grad 0 0 17e3   ran z -9 -10.1  ;soilc
ini syy -8e3 grad 0 0 13.5e3  ran z -9 -10.1         
ini sxx -8e3 grad 0 0 13.5e3  ran z -9 -10.1        

solve

save 初始平衡.sav               


rest 初始平衡.sav
ini xd 0 yd 0 zd 0 xv 0 yv 0 zv 0
set mech on fluid off
del range group soil1 z -10 0

;生成桩
gen zone brick size 4 13 4 p0 (-.3,-.3,15)  p1 (.3,-.3 ,15) p2 (-.3,-.3,2) p3 (-.3,.3,15) group pile2


interface 1 face range x -0.29 -.31 y -.29 -.15 z -10 0
interface 1 face range x -0.29 -.31 y  .15  .29 z -10 0

interface 1 face range x  0.29  .31 y -.29 -.15 z -10 0
interface 1 face range x  0.29  .31 y  .15  .29 z -10 0

interface 1 face range y -0.29 -.31 x -.29 .29 z -10 0
interface 1 face range y  0.29  .31 x -.29 .29 z -10 0

interface 2 face range x -0.29 -.31 y  -.15  .15 z -10 0
interface 2 face range x  0.29  .31 y  -.15  .15 z -10 0

;interface 1 face range x  -.29  .29 y -.29 .29 z -10.1 -9.9

int 1 prop kn 1e8 ks 1e8 f 20 c 0    range z -1 -9
int 1 prop kn 1e8 ks 1e8 f 15 c 7e3 range z 0 -1
int 1 prop kn 1e8 ks 1e8 f 15 c 10e3 range z -9 -10.1

int 2 prop kn 1e8 ks 1e8 f 20 c 0    range z -1 -9
int 2 prop kn 1e8 ks 1e8 f 15 c 7e3 range z 0 -1
int 2 prop kn 1e8 ks 1e8 f 15 c 10e3 range z -9 -10.1

interface 1 perm = off
interface 2 perm = off  
interface 1 effe = off
interface 2 effe = off


ini z add -12 ran gro pile2

fix z ran z -10.1 -9.9
;分组
group soil range  group pile2 not
group soila range group soil z 0 -1
group soilb range group soil z -1 -9
group soilc range group soil z -9 -10.1

;桩属性
group pile1 range group pile2 z 2 3
group pile  range group pile2 z 2 -10
model elastic ran gro pile
prop bu 1.67e10 sh 7.69e9  ran gro pile
ini dens 2400 ran gro pile ;z 0 2

;桩顶重物

;属性
model elastic ran gro pile1
prop bu 1.67e10 sh 7.69e9  ran gro pile1
ini dens 2400 ran gro pile1


;土属性
model mo range group soila any group soilb any group soilc any
prop dens 1040 por 0.66 bu 14.71e6 sh 5.640e6  f 25 c 10e3   range group soila
prop dens 1440 por 0.46 bu 2.941e7 sh 1.128e7  f 30 c  0     range group soilb
prop dens 1100 por 0.60 bu 14.71e6 sh 5.640e6  f 25 c 15e3   range group soilc



;pau
;水属性
model fl_iso range group  soila any group soilb any group soilc any ;group toushui any
model fl_nul range group pile   a group pile1 a

prop   perm 1e-13  range group soila
prop   perm 1e-8   range group soilb
prop   perm 1e-13  range group soilc

ini fdens 1000 ran gro  soila any group soilb any group soilc any

ini pp 0 grad 0 0 -10e3  ran z 0 -10.1 gro  soila any group soilb any group soilc any

set mech on fluid off
ini fmod 0  ran gro  soila any group soilb any group soilc any  
ini sat = 0 range group pile any group pile1 any
;pau
def add_top
  loop n (1,10)
    ini_add = n * 2400
    save_file = string(n) + 'G.sav'
    command
      ini dens ini_add ran gro pile1
      set mech rat 1e-5
      solve
      save save_file
    endcommand
  endloop
end
add_top

set mech rat 1e-6
solve
save 加桩后初始应力.sav


rest 加桩后初始应力.sav

ini state 0
set dyn on mech on fluid on
set dyn mult on
set fluid time 0
set mech  rat 1e-100
set fluid pcut on ;土体不产生负孔压
model finn ran group soilb
prop dens 1440 por 0.46 bu 2.941e7 sh 1.128e7  f 30 c 0  perm 1e-8 ran group soilb  
prop ff_latency=50 ran group soilb                                                  
                                    
prop ff_switch = 0 ran group soilb                                                  
prop ff_c1=0.8   ff_c2=0.79 ran group soilb                                         
prop ff_c3=0.45  ff_c4=0.73 ran group soilb                                         


ini fmod 2e8 ran gro soila any group soilb any group soilc any

;阻尼

ini  damp local 0.1571 range group pile   a group pile1 a
ini  damp local 0.3142 range group soila a group soilc a
ini  damp local 0.4713 range group soilb

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

free x range z -10.1 -9.9
;底部速度施加
def setup
   freq=5.0
   ampl=2
   omega = 2.0 * pi * freq
end
setup
def sine_wave
  vv = ampl/omega*cos(omega*dytime)
  if dytime < 2.0
    sine_wave = dytime / 2.0 * vv
    else
    if dytime < 15.0
      sine_wave = vv  
      else
      if dytime <= 17.0
          sine_wave = (17.0 - dytime) / 2.0 * vv
      endif
    endif
  endif
  if dytime > 17.0
     sine_wave = 0.0
  endif
end

apply  xvel=1.0 hist sine_wave range z -9.9  -10.1


apply ff
call ppr.dat


;历史变量监测
his reset
hist id 1 unbal
hist id 2 dytime
hist id 3 sine_wave
           
            

;加速度(底部输入加速度时程)
hist id 111 gp xacc    0   0   -10
hist id 112 gp xacc   0.3  0   -10
hist id 113 gp xacc  1.4   0   -10
hist id 114 gp xacc  2.5   0   -10
hist id 115 gp xacc  3.6   0   -10
hist id 116 gp xacc  4.7   0   -10
hist id 117 gp xacc  5.9   0   -10
hist id 118 gp xacc    8   0   -10
hist id 119 gp xacc   12   0   -10
hist id 120 gp xacc   16   0   -10
           
           
def get_ppr_zone
  whilestepping
   
   p_z43 = z_near(1 , 0.1, -5.5)
  ppr_431 = z_extra(p_z43,2)  
  ppr_432 = z_extra(p_z43,3)
  ppr_433 = z_extra(p_z43,6)  
  ppr_434 = z_extra(p_z43,7)
   
   p_z44 = z_near(2, 0.1, -5.5)
  ppr_441 = z_extra(p_z44,2)  
  ppr_442 = z_extra(p_z44,3)  
  ppr_443 = z_extra(p_z44,6)  
  ppr_444 = z_extra(p_z44,7)
                                                        
   p_z46 = z_near(4, 0.1, -5.5)
  ppr_461 = z_extra(p_z46,2)  
  ppr_462 = z_extra(p_z46,3)  
  ppr_463 = z_extra(p_z46,6)  
  ppr_464 = z_extra(p_z46,7)
                             
   p_z49 = z_near(12, 0.1, -5.5)
  ppr_491 = z_extra(p_z49,2)  
  ppr_492 = z_extra(p_z49,3)  
  ppr_493 = z_extra(p_z49,6)  
  ppr_494 = z_extra(p_z49,7)
  
end

;监测 超孔压比、超孔压、有效应力

his id 639 ppr_431          ;(1  0.1 -5.5)
his id 640 ppr_432
his id 641 ppr_433
his id 642 ppr_434

his id 643 ppr_441          ;(2 0.1 -5.5)
his id 644 ppr_442
his id 645 ppr_443
his id 646 ppr_444

his id 651 ppr_461          ;(4 0.1 -5.5)
his id 652 ppr_462
his id 653 ppr_463
his id 654 ppr_464

his id 663 ppr_491          ;(12 0.1 -5.5)
his id 664 ppr_492   
his id 665 ppr_493   
his id 666 ppr_494   

set hist_rep 300                    
set dyn dt 4e-5                     
def solve_ages                     
   loop n (1,20)  
     time_n = n   
     save_file = 'Dyn-' + string(n) + 'x1s.sav'
     command                        
       sol age time_n                     
       save save_file               
     endcommand                     
   endloop                           
end                                 
solve_ages                    
发表于 2016-7-21 16:45:36 | 显示全部楼层 来自 湖南长沙
Simdroid开发平台
你好哥们 你的实体单元建立的桩弯矩是怎么提取的?可以交流一下吗/
回复 不支持

使用道具 举报

 楼主| 发表于 2016-7-21 18:11:37 | 显示全部楼层 来自 浙江舟山
你好,兄弟,我现在还是新手,这是《FLAC3D基础与工程实例》一书中第十七章的命令流文件,具体的我还要消化消化
回复 不支持

使用道具 举报

您需要登录后才可以回帖 登录 | 注册

本版积分规则

Archiver|小黑屋|联系我们|仿真互动网 ( 京ICP备15048925号-7 )

GMT+8, 2024-3-28 17:43 , Processed in 0.028818 second(s), 11 queries , Gzip On, MemCache On.

Powered by Discuz! X3.5 Licensed

© 2001-2024 Discuz! Team.

快速回复 返回顶部 返回列表