动力分析总是出错，网格变形很大（附简单源程序）
大家好，我建立了很简单的一个模型来分析自由场动力问题，静力分析是顺利进行的，但是施加了finn后，运行程序就会网格变形很大，只要是网格两侧顶部，不知道是什么原因，请大家指点，谢谢！config dyn gw ats
def wave
wave = 0.5 * (1.0 - cos(2.0*pi*freq * dytime))
end
set freq=2.0
grid 16 10
mod m
prop dens 1500 shear 2.34e8 bulk 5.09e8
prop fric 30
set grav 10
set flow off
prop poros 0.5
water dens 1000 bulk 2e9 ;tens 1e10
fix x i=1
fix x i=17
fix y x j=1
;
set dyn off
hist unbal
hist ydis i 5 j 5
hist ydis i 5 j 11
call ininv.fis
set wth 10 k0x 0.428 k0z 0.428
ininv
solve elastic
m finn
prop dens 1500 shear 2.34e8 bulk 5.09e8
prop fric 30
prop poros 0.5
prop ff_latency=100
prop ff_switch = 1
prop ff_c1=0.08 ff_c2=5
set dyn on
set large
apply xacc 1 hist wave yvel=0j 1
ini xvel=0 yvel=0 xdis=0 ydis=0
apply ff
;
set dytime 0
hist reset
hist dytime
hist xvel i 1 j 11
hist xvel i 8 j 11
hist xvel i 17 j 11
hist xvel i 1 j 1
hist xvel i 8 j 1
hist xvel i 17 j 1
hist xacc i 1 j 11
hist xacc i 8 j 11
hist xacc i 17 j 11
hist xacc i 1 j 1
hist xacc i 8 j 1
hist xacc i 17 j 1
hist pp i 1 j 5
hist pp i 8 j 5
hist pp i 16 j 5
hist pp i 1 j 1
hist pp i 8 j 1
hist pp i 16 j 1
set dy_damp=rayleigh 0.05 4
solve dytime 5 本来觉得是较为简单的一种分析情况，没想到也有这么多的问题. 希望你能够指导下我做动力分析啊！QQ：312169218 先把静力部分算对。你的文件静力计算都有提示错误信息的。
config dyn gw ats
def wave
wave = 0.5 * (1.0 - cos(2.0*pi*freq * dytime))
end
set freq=2.0
g 16 10
mod m
prop dens 1500 shear 2.34e8 bulk 5.09e8
prop fric 30
set grav 10
prop poros 0.5
water dens 1000 ;tens 1e10
fix x i=1
fix x i=17
fix y x j=1
set dyn off
hist unba
hist ydis i 5 j 5
hist ydis i 5 j 11
call ininv.fis
set flow off
set wth=10 k0x=0.48 k0z=0.48
ininv
solve elastic
页:
[1]