在陈文化老师写的那本ansys参考书上看到关于车辆过桥的问题,我用FLAC3D写了简单程序,也想试试,可问题多多,现把命令流贴上,敬请大家讨论,指点迷津!
sel beam id=1 begin=( 1, 0, 0) end=( 9, 0, 0) nseg=9
sel beam id=1 prop emod=3e10 nu=0.20 density=2500 xcarea=0.1 xcj=0 xciy=0.000083 xciz=0.0083
sel node fix x y z xr yr zr range id=1
sel node fix x y z xr yr zr range id=2 ;为简单见,就直接全部约束了两端
set grav 0 0 -10
hist unbal
solve
save static.sav
res static.sav
config dynamic
sel node init xdisp 0
sel node init ydisp 0
sel node init zdisp 0
sel node init xrdisp 0
sel node init yrdisp 0
sel node init zrdisp 0
sel node init xvel 0
sel node init yvel 0
sel node init zvel 0
sel node init xrvel 0
sel node init yrvel 0
sel node init zrvel 0
set dyn on
;---------定义荷载-----------
def load
loop i (1,8)
t1=0.02*i
t2=0.02*(i-1)
if dytime<=t1 then
if dytime>t2 then
command
sel node apply force=(0,0,-1e6) range id=i+2
load=force
endcommand
endif
endif
endloop
end
hist reset
hist unbal
hist dytime
hist sel node zvel id=5 ;好像不能监测加速度!
solve age 0.3
save dynamic.sav
现在的问题是:荷载施加不上去,用print fish查看,或在FLAC3D中通过如下途径看:
plot--base show
plotitems---add-----structural element---applied nodal force,显示节点荷载为零
因此问题的关键本人认为在于这几句:
command
sel node apply force=(0,0,-1e6) range id=i+2
load=force
endcommand
请教下该怎么处理?