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

!!!一简支梁长10m，L/2、L/4、边支点处的单元E损伤!!!simplebeamMFfinish/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损伤/prep7b=0.3!!!梁截面参数h=0.5et,1,beam3r,1,b*h,b*h**3/12,h!!!area=b*h,I=b*h^3/12,prxybm=0.2mp,dens,1,robm!!!定义3种弹性模量，未损伤E=3.25E+10mp,prxy,1,prxybmmp,ex,1,3.25E+10mp,dens,2,robm!!!L/2处单元损伤，E=exL2为修正量mp,prxy,2,prxybmmp,ex,2,exL2mp,dens,3,robm!!!L/4处单元损伤，E=exL4为修正量mp,prxy,3,prxybmmp,ex,3,exL4mp,dens,4,robm!!!边支点单元损伤，E=ex20为修正量mp,prxy,4,prxybmmp,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+1latt,1,1,1,lsel,none*enddolsel,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,allsellesize,all,1!!!线的尺寸划分,0.5m一个单元lmesh,all!!!对梁进行网格划分dk,1,uy,,,,uz,rotx,roty!!!增加约束,简支梁dk,21,uy,,,,ux,uz,rotx,rotyfinish/solu!模态求解allselantype,modalacel,,9.8modopt,lanb,8,1,500!频率提取从1Hz到500Hz且阶数不超过8mxpand,8!计算相应频率的振型solvefinish/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,1shape(j,i)=uy(j)*enddoset,next*enddofinish*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