ansys模态分析,提取竖向位移

一简支梁长10m，L/2、L/4、边支点处的单元E损伤

simplebeamMF

finish

/clear

/config,nproc,3

----------------参数定义---------------

robm = 2500梁的密度

exL2 = 3.25E+10*(1-0.0) L/2处的一个单元损伤

exL4 = 3.25E+10*(1-0.0) L/4处的一个单元损伤

ex20 = 3.25E+10*(1-0.0)边支点单元20损伤

/prep7

b = 0.3梁截面参数

h = 0.5

et,1,beam3

r ,1,b*h, b*h**3/12, h area=b*h, I=b*h^3/12,

prxybm=0.2

mp,dens,1,robm定义3种弹性模量，未损伤E=3.25E+10 mp,prxy,1,prxybm

mp,ex ,1,3.25E+10

mp,dens,2,robm L/2处单元损伤，E=exL2为修正量mp,prxy,2,prxybm

mp,ex ,2,exL2

mp,dens,3,robm L/4处单元损伤，E=exL4为修正量mp,prxy,3,prxybm

mp,ex ,3,exL4

mp,dens,4,robm边支点单元损伤，E=ex20为修正量mp,prxy,4,prxybm

mp,ex ,4,ex20

*do,i,1,21,1立21个关键点,0.5m一个单元建

k,i,(i-1)*0.5

*enddo

*do,i,1,20,1建立20个单元，全部用第1种材料

l,i,i+1

latt,1,1,1,

lsel,none

*enddo

lsel,s,,,5修改第5单元(L/4处)的材料常数，第3种材料latt,3,1,1,

lsel,s,,,10修改第10单元(L/2处)的材料常数，第2种材料latt,2,1,1,

lsel,s,,,20修改第20单元(支点处)的材料常数，第4种材料latt,4,1,1,

allsel

lesize,all,1线的尺寸划分,0.5m一个单元

lmesh,all对梁进行网格划分

dk, 1,uy,,,,uz,rotx,roty增加约束,简支梁

dk,21,uy,,,,ux,uz,rotx,roty

finish

/solu!模态求解

allsel

antype,modal

acel,,9.8

modopt,lanb,8,1,500!频率提取从1Hz到500Hz且阶数不超过8 mxpand,8!计算相应频率的振型

solve

finish

/post1

set,list

set,first

pldisp,1

set,next

pldisp,1

set,1,5

pldisp,1

*dim,freq ,array,6提取前6阶竖弯频率和振型

*dim,shape,array,21,8

*do,i,1,3,1

*get,freq(i),mode,i,freq提取前3阶频率

*enddo

*get,freq(4),mode,5,freq提取5，6阶频率，4阶为平动振型，不提取*get,freq(5),mode,6,freq

*get,freq(6),mode,8,freq第7阶也有平动

set,first

*do,i,1,8,1提取前8阶振型,4、7阶为平动，需去掉

*do,j,1,21,1

shape(j,i)=uy(j)

*enddo

set,next

*enddo

finish

*cfopen,outfreq,txt,写出频率

*vwrite,freq(1)

(f12.6)

*cfclos

*cfopen,outshape,txt,写出振型,第4阶为平动振型，不输出

*vwrite,shape(1,1),shape(1,2),shape(1,3),shape(1,5),shape(1,6),shape(1,8)

(6f12.6)

*cfclos