本帖最后由 niej2008 于 2012-8-19 16:08 编辑
这是我的命令流,麻烦老师帮我看下,梁结构很简单,沿z轴长度为1m、半径为0.02m的圆截面直梁,网格划分100等份,作用空间载荷,载荷从右往左移动速度为0.01m/s。只运行模态分析时一阶固频为125,运行瞬态分析后再查看固频就变为0.25000E-10 ,以及速度、加速度曲线异常,是我程序哪里错了吗
fini
/cle
/config,nres,2000
*set,long,1
*set,ne,100
*set,dl,long/ne
*set,nn,ne+1
*set,ffx,75.08
*set,ffy,170.31
*set,ffz,75.20
*set,vv,0.01
*set,dt,dl/vv
/prep7
et,1,beam44
mp,ex,1,2.07E11
mp,nuxy,1,0.28
mp,dens,1,7850
sectype,1,beam,csolid
secdata,0.02,30,6
*do,i,1,nn
n,i,,,(i-1)*dl
*enddo
*do,i,1,ne
e,i,i+1
*enddo
fini
/sol
antype,modal
modopt,lanb,5
mxpand,5
d,1,all
d,nn,ux,,,,,uy
allsel,all
solve
save
fini
/post1
set,list
fini
!若不加下面语句,一阶固有频率为125,加上后一阶频率为0.25000E-10 ,相差得很大,不知道是不是程序哪编错了
/sol
antype,4
trnopt,full
timint,off
outres,basic,all
kbc,1
time,1e-10
nsub,4,8
sstif,on
acel,0,9.8
alls
solve
timint,on
*do,ii,1,nn,1
outres,all,all
time,ii*dt
nsubst,5
fdele,all,all
nsel,s,loc,z,1-(ii-1)*dl
f,all,fx,ffx
f,all,fx,-ffy
f,all,fx,-ffz
allsel
solve
*enddo
save
fini
/post1
set,ne/2+1,last
/eshape,1
plnsol,u,x,0,1
plnsol,s,eqv,0,1
allsel
fini
/post26
NUMVAR,200
nsol,2,ne/2+1,u,x,ux_m
xvar,1
plvar,2
deriv,3,2,1,,,,,1
plvar,3
deriv,4,3,1,,,,,1
plvar,4 |