- 积分
- 0
- 注册时间
- 2007-11-29
- 仿真币
-
- 最后登录
- 1970-1-1
|
本帖最后由 cloud328 于 2010-6-23 19:40 编辑
flac3d 手册上动力分析的原例子,怎么运行时提示“unrecognized parameter 4(not)”???盼高手指点迷津!!!不胜感激!急急急!!!
config dynamic
title
Verification of FLAC-3D with SHAKE-Linear Elastic Case
;--- Grid generation and model properties ---
gen zone brick size 1,1,16
ini x mul 3.048 y mul 3.048 z mul 3.048
model elastic
prop bulk 150e6 shear 150e6
prop bulk 300e6 shear 300e6 range z=24.384,36.576
ini dens 1800
ini dens 2000 range z=24.384,36.576
range name=bottom z=-0.5,0.5
;--- Pre-compute some constants ---
def setup
omega=6.0*pi
alfa=2.2
beta=3.75e-1
gamma=8.0
n12=gp_near(0.0,0.0,33.528)
n13=gp_near(0.0,0.0,36.576)
n14=gp_near(0.0,0.0,39.624)
end
setup
;--- FISH function to calculate shear strain in layers 12 and 13 ---
def shrstr12
shrstr12=(gp_xdisp(n13)-gp_xdisp(n12))/(gp_zpos(n13)-gp_zpos(n12))
shrstr13=(gp_xdisp(n14)-gp_xdisp(n13))/(gp_zpos(n14)-gp_zpos(n13))
end
;--- Histories ---
hist unbal
his dytime
his zone sxz 0,0,37
his zone sxz 0,0,34
his gp xacc 0,0,48.768
his gp xacc 0,0,0
his shrstr12
his shrstr13
;--- Boundary Conditions ---
fix y z range=(not bottom)
;--- FISH function to generate input wave ---
def acc_p
acc_p=sqrt(beta*exp(-alfa*dytime)*dytime^gamma)*sin(omega*dytime)
end
;--- Application of acceleration ---
apply xacc 1.0 his acc10 yacc=0 zacc=0 range=bottom
set dyn=on, dyn damp=rayleigh 0.1 3.0
solve age 14.0
save shake.sav |
本帖子中包含更多资源
您需要 登录 才可以下载或查看,没有账号?注册
×
|