model new
zone import 'a.f3grid'
zone face skin
model history name '1' mechanical ratio-average
model history name '2' mechanical unbalanced-maximum
;约束边界条件
zone gridpoint fix velocity-y
zone gridpoint fix velocity-z range position-z -0.1 0.1
zone gridpoint fix velocity-x range position-z -0.1 0.1
zone gridpoint fix velocity-x range position-x -0.1 0.1
zone gridpoint fix velocity-x range position-x 109.9 110.1
zone cmodel assign mohr-coulomb
zone property bulk 0.196e9 she 0.0752e9 fric 25 cohesion 10e10 ten 1e10 density 2000 range group 'RRJC'
zone property bulk 1e9 she 0.6e9 fric 38 cohesion 10e10 ten 1e10 density 2500 range group 'YT'
model gravity (0, 0, -10)
model solve
zone property bulk 0.196e9 she 0.0752e9 fric 25 cohesion 0.08e6 ten 0.02e6 density 2000 range group 'RRJC'
zone property bulk 1e9 she 0.6e9 fric 38 cohesion 0.6e6 ten 0.5e6 density 2500 range group 'YT'
model solve
model save '01-ini'
以上是生成初始地应力
zone gridpoint initialize displacement 0 0 0
zone gridpoint initialize velocity 0 0 0
zone initialize state 0
model configure dynamic
model dynamic active on
zone dynamic multi-step on
model largestrain on
zone history name 'N11-acc-x' acc-x position (25,0.5,30)
zone history name 'M21-acc-x' acc-x position (80,0.5,30)
model history dynamic time-total
zone gridpoint free velocity-x range position-z -0.1 0.1
zone dynamic damping local 0.157
fish define setup
local freq = 1.5
global omega = 2.0 * math.pi * freq
end
@setup
fish define wave
wave = math.sin(omega * dynamic.time.total)
end
fish history @wave
zone face apply acceleration-x 1 fish @wave range p-z -1 1
zone dynamic free-field on
model solve time-total 4
model save '03-dyn-1'
以上是输入正弦波进行动力计算
请问动力计算过程有什么错误？为什么不收敛呢？
探究动力响应规律需要对软弱夹层和岩体直接运用interface界面？
