（急）请教Ansys模拟拉索振动的问题！

knighthit

本人用ansys模拟斜拉索在简谐激励荷载作用下拉索的振动问题，开始的时候定义了线性阻尼，后来把阻尼去掉之后在13阶荷载的激励下希望拉索发生13阶振动，以此测试斜拉索的对数衰减率。结果发现去掉阻尼后拉索发生的振动中含有13阶成分，将13阶模态提取出来后，发现13阶模态自身的对数衰减率已经达到了0.01，这与理想的无阻尼衰减相矛盾，希望高人给予指点。
有好心人能帮助一下请联系我：QQ：39727517.
下面是命令流：

/TITLE,DAMPER Analysis
/CONFIG,NRES,1000000
/PREP7
Pi=4*ATAN(1)
!****************输入常数*********
EX=1.95E11
DENSY=8367.922962
A=0.009275
T=4611257.934         !索的张力
ISTRN=T/A/EX
cita=26.7305/180*Pi
L=448.04
!*******************阻尼器常数**************
DK1=0          !大阻尼器刚度
DC1=0          !大阻尼器阻尼系数（2阶优化）
Dk2=0          !小阻尼器刚度
DC2=0          !小阻尼器阻尼系数
!*******************
fm=DENSY*A*9.81*L**2/8/T
ET,1,LINK10
ET,2,COMBIN14
R,1,A,ISTRN
R,2,DK1,DC1,,
R,3,DK2,DC2,,
MP,EX,1,EX   
MP,PRXY,1,0.3                 
MP,DENS,1,DENSY   
!************定义节点坐标****
*DIM,XD,,150
*DIM,YD,,150
*DO,I,1,4,1
  XD(I)=10.66*COS(cita)/3*(I-1)
  YD(I)=(4*fm/L**2/COS(cita)/COS(cita)*XD(I)+(TAN(cita)-4*fm/L/COS(cita)))*XD(I)
*ENDDO
*DO,I,5,101,1
XD(I)=(L-11)*COS(cita)/97*(I-4)+XD(4)
YD(I)=(4*fm/L**2/COS(cita)/COS(cita)*XD(I)+(TAN(cita)-4*fm/L/COS(cita)))*XD(I)
*ENDDO
!**********生成单元***********
*DO,I,1,101,1
N,I,XD(I),YD(I),0 bn
*ENDDO
N,102,XD(4)/COS(cita)**2,0,0
N,103,XD(2)/COS(cita)**2,0,0
*del, XD,
*del, YD,
TYPE,1
REAL,1
MAT,1
*DO,I,1,100,1
EN,I,I,I+1
*ENDDO
TYPE,2
REAL,2
EN,101,4,102
TYPE,2
REAL,3
EN,102,2,103
D,1  , , , , , ,UX,UY,UZ, , ,
D,101, , , , , ,UX,UY,UZ, , ,
D,102, , , , , ,UX,UY,UZ, , ,
D,103, , , , , ,UX,UY,UZ, , ,
sstif,on
FINISH
!********静力分析***************
/SOLU
ANTYPE,STATIC
ACEL,0,9.81,0
SOLVE
OUTRES,ALL,ALL   
FINISH

!********模态分析***************
/SOLU
ANTYPE,MODAL
MODOPT,LANB,30,0.1,30
pstres,on                    !声明预应力
mxpand,30                    !输出模态阶数
solve
finish
!********瞬态分析***************
/solu
antype,trans                !瞬态分析
trnopt,full                 !指定瞬态分析选项(求解方法)
NLGEOM,ON                   !大变形
sstif,on                    !在非线性分析中激活预应力效果
timint,on
autots,on                   !使用自动时间步长
nsubst,1                    !每一荷载步的子步数  如果使用autots时间步长，则为定义第一荷载步的长度   
kbc,1                       !制定载荷为阶跃载荷(递增载荷,0)
IC,y,all,0,0         
*do,i,1,2800,1
    time,0.04*i                !荷载步结束时间
    *if,i,le,300,then              
        f,44,fy,-30000*sin(2*Pi*3.5671*i*0.04)*cos(cita)
        f,44,fx,30000*sin(2*Pi*3.5671*i*0.04)*sin(cita)
    *else
        f,44,fy,0         
        f,44,fx,0
    *endif
solve
lswrite,i                   !将荷载步的结果写入文件
*enddo
outres,all,all              !规定写入数据库的求解信息(所有求解项，每一步都写)
finish
下面是15号节点的振动信号
提取出的13阶模态
file:///C:/Documents%20and%20Settings/YDJ/桌面/三阶振动.JPG

[ 本帖最后由 16443 于 2008-4-28 09:19 编辑 ]

knighthit

提取的13阶振动
file:///C:/Documents%20and%20Settings/YDJ/桌面/三阶振动.JPG

aircraft

楼主,你的建模与划分网格过程,
没仔细看.
看了下你的求解过程:
提几点意见:
1.模态分析打开了PSTRES开关,
为啥静力分析阶段却没有看到这个开关.
2.模态分析用的是PSTRES开关
而瞬态分析用的是SSTIF开关
也就是说:模态分析考虑的线性问题
而瞬态分析考虑的非线性问题
这里我有点不明白.
还望楼主指教!

knighthit

我实际主要关心的是后面的瞬态分析，前面模态分析只是为了计算频率，为了后面瞬态分析的激励频率提供依据。