考虑轴力的弹性支撑变截面梁的模态分析问题

trueghost

最近做了一个变截面梁的模态分析。
在加轴力之前分析结果和理论结果符合，
但加了轴力之后，总感觉不大对劲。轴力在0到2000N之间时，频率与未加轴力一样， 但轴力位于2000到2650N之间时，频率明显变小，直至屈曲。下面是我的命令流，请各位指教。
finish
/CLEAR,START
pd=0
/prep7
ANTYPE, STATIC
PSTRES, ON
*do,i,1,20,1
et,i,beam54okok.org !定义变截面桥墩单元
mp,dens,i,2700
mp,ex,i,7e10
mp,prxy,i,0.33
mp,gxy,i,2.63e10
R,i,2.755e-4*(1-(i-1)*pd/20),2.072e-9*(1-(i-1)*pd/20)**3,0.00475*(1-(i-1)*pd/20),0.00475*(1-(i-1)*pd/20),
RMODIF,i,9, , ,
RMODIF,i,14, ,
RMODIF,i,5,2.755e-4*(1-i*pd/20),2.072e-9*(1-i*pd/20)**3,0.00475*(1-i*pd/20),0.00475*(1-i*pd/20),
RMODIF,i,11, , ,
RMODIF,i,15, ,
RMODIF,i,13, ,
RMODIF,i,16, , , ,
*enddo
et,21,combin14
keyopt,21,2,1 !纵向弹簧单元
et,22,combin14
keyopt,22,2,!扭转弹簧单元
okok.org
r,21,8.2e3
r,22,8.2e3
r,23,5.63e2

n,1,0,0,0
n,21,0,0.262,0
fill
n,22,0,0,0
n,23,0,0,0
n,24,0,0.262,0

*do,i,1,20,1
type,i
real,i
e,i,i+1
*enddo

type,21
real,21
e,1,22

type,21
real,22
e,21,24

type,22
real,23
e,1,23

f,21,fy,-2600
d,22,all
d,23,all,all
d,24,all
d,1,uy
finish

!（4）进入求解模块, 指定求解选项并求解
/SOLU
OUTPR, BASIC, 1
SOLVE
FINISH
! （5）进入一般后处理模块, 提取频率计算结果
/POST1
ETABLE, STRS, LS, 1
*GET, STRSS, ELEM, 20, ETAB, STRS
*GET, DEF, NODE, 21, U, y
! （6）重新进入求解模块, 用模态求解方法求解
/SOLU
ANTYPE, MODAL
PSTRES, ON
MXPAND, 3
MODOPT, LANB, 3
SOLVE
! （7）提取频率计算结果
*GET, FREQ1, MODE, 1, FREQ
*GET, FREQ2, MODE, 2, FREQ
*GET, FREQ3, MODE, 3, FREQ
*STATUS, PARM
FINISH

[ 本帖最后由 trueghost 于 2006-7-16 13:51 编辑 ]

sysh320

动力学中有这个概念的,可以看看

trueghost

sysh320兄，概念方面我自己觉得已经差不多弄清了
考虑轴力与不考虑轴力的关系为：w^2=w0^2(1-N/Ncr)
但是算出来的结果就是不对
不知为什么？

[ 本帖最后由 trueghost 于 2006-7-16 14:00 编辑 ]

trueghost

另外我用屈曲分析算出来的临界轴力  与一阶频率为0时的轴力 相等

ME!

trueghost 发表于 2006-7-17 09:32

                               
登录/注册后可看大图

另外我用屈曲分析算出来的临界轴力  与一阶频率为0时的轴力 相等

屈曲分析不是应该antype：buckle or 2 屈曲分析吗？modal是模态分析吧！