地震波加载

machi

在动力计算的时候出现这个错误，显示command那个地方有错误？
Illegal geometry in zone 1

GZhang2468

我的经验是自行设定的动力时间步过大导致的

dongxx0221

可能是zone的参数赋值问题

dingdinggou110

应该是模型的赋值没有赋上，你找一下zone的位置。

machi

dingdinggou110 发表于 2017-3-20 15:27
应该是模型的赋值没有赋上，你找一下zone的位置。

请问5.0的模型建立对网格有什么要求吗？我的模型分成7组，网格简化成size1 1 1，很快可以算完，网格很多时就有问题，我应该在哪方面修改呢？会不会是网格的问题

machi

dongxx0221 发表于 2017-3-19 13:11
可能是zone的参数赋值问题

请问5.0的模型建立对网格有什么要求吗？我的模型分成7组，网格简化成size1 1 1，很快可以算完，网格很多时就有问题，我应该在哪方面修改呢？会不会是网格的问题

machi

是变形过大的问题，这是边界条件的原因吗？参数我看了一下，没问题

AirTure

是不是底部边界锁定了

machi

是模型网格不均匀的问题，请问Flac在动力计算中对网格有什么要求吗？怎么改变

machi

machi 发表于 2017-3-22 11:17
是变形过大的问题，这是边界条件的原因吗？参数我看了一下，没问题

;;;;;;;jinglibufen;;;;;;;
;;;;;;;moxingwanggejianli;;;;;;;;
gen zone uwedge  p0(-8 20 0)  p1(-488 20 0) p2(-8 -20 0) p3(-8 260 240)  p4(-488 -20 0)&
p5(-8 -260 240)  size 20 20 24 group 1;;1-1  

gen zone brick p0(-8 -20 0) p1(-6 -20 0) p2(-8 20 0) p3(-8 -260 240) p4(-6 20 0)&
p5(-8 260 240) p6(-6 -260 240) p7(-6 260 240) size 1 20 24  group 2;;2-1

gen zone brick p0(-6 -20 0) p1(-5 -20 0) p2(-6 20 0) p3(-6 -260 240) p4(-5 20 0) size 1 20 24  group 3;;3-1
;&
;p5(-5 260 240) p6(-5 -260 240) p7(-5 260 240)


gen zone brick p0(-5 -20 0) p1(5 -20 0) p2(-5 20 0) p3(-5 -260 240) p4(5 20 0)&
p5(-5 260 240) p6(5 -260 240) p7(5 260 240) size 5 20 24  group 4;;4

gen zone brick p0(5 -20 0) p1(6 -20 0) p2(5 20 0) p3(5 -260 240) p4(6 20 0)&
p5(5 260 240) p6(6 -260 240) p7(6 260 240) size 1 20 24  group 3;;3-2

gen zone brick p0(6 -20 0) p1(8 -20 0) p2(6 20 0) p3(6 -260 240) p4(8 20 0)&
p5(6 260 240) p6(8 -260 240) p7(8 260 240) size 1 20 24  group 2;;2-2

gen zone uwedge  p0(8 -20 0)  p1(488 -20 0) p2(8 20 0) p3(8 -260 240)  p4(488 20 0)&
p5(8 260 240)  size 20 20 24 group 1;;1-2

;;;;;;;;;;;;;;;
range name 1 group 1
range name 2 group 2
range name 3 group 3
range name 4 group 4


model mohr
prop cohesion = 10000.0 friction = 42.3 dilation = 12.0 tension = 0.0 range 1
prop cohesion = 30000.0 friction = 40.4 dilation = 12.0 tension = 0.0 range 2
prop cohesion = 37000.0 friction = 36.3 dilation = 10.0 tension = 0.0 range 3
prop cohesion = 30000.0 friction = 25 dilation = 14.0 tension = 0.0 range 4
ini den 2260 range 1
ini den 2210 range 2
ini den 2320 range 3
ini den 2350 range 4


;model fl_iso ran gro 1 gro 2 gro 3 gro 4 gro 6 a gro 7 a
plot zone
;plot block group
;pause

;;;;;;;????????????????????boundary conditions;;;;;;;;
fix x y z range plane norm 0,1,-1.0 origin -8 260 240 dist 0.1
fix x y z range z -0.01 0.01
fix x y z range plane norm 0,-1,-1.0 origin -8 -260 240 dist 0.1

his id 1 unbal
set grav=9.81 ;large
;set large
set dyn=off
;pause
;solve

set dynamic on
ini  xdis=0 ydis=0 zdis=0

table 100 read  shuiping.acc  ;;;;;;;;;;;read acc history from located directory .acc file


;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;apply xacc history on bottom of model
free x
apply xacc 0.1 his table 100 dynamic range  z -0.01 0.01
apply xacc 0.1 his table 100 dynamic range  plane norm 0,1,-1.2 origin -8 394 295 dist 0.1     z=0,295                                
apply xacc 0.1 his table 100 dynamic range  plane norm 0,-1,-2 origin -8 -335 295 dist 0.1    z=0,295
ini damp hyst default -3.156 1.904
set dyn damp rayleigh=0.005 1.5 stiffness                                          
set large
set dynamic time=0.0
hist nstep 100
set dynamic dt 0.0005

def solve_ages                                    
  loop n(1,15)                                    
    sol_n = n * 2.0    ;分15步保存，每2秒保存一个结果
    sol_file = 'Dyn-x' + string(n) + 'x2s.sav'
    command
      sol age sol_n
      save sol_file            
    endcommand                                    
    endloop                                       
end                                               
solve_ages                                       

dongxx0221

set dynamic dt 0.0005  可能dt值太大；
加 set dyn multi on

dongxx0221

no valid gp in the range

machi

dongxx0221 发表于 2017-4-8 14:33
no valid gp in the range

怎么发文件给您，这个命令流复制有问题

machi

dongxx0221 发表于 2017-4-8 14:33
no valid gp in the range

没剪切模量，table100在下面

hanzp

GZhang2468 发表于 2020-11-29 21:20
我的经验是自行设定的动力时间步过大导致的

可以的,我试了一下，果然你说的对哦