資源描述：
《ansys模態(tài)分析-提取豎向位移.docx》由會員上傳分享，免費(fèi)在線閱讀，更多相關(guān)內(nèi)容在教育資源-天天文庫。
1、!!!一簡支梁長10m，L/2、L/4、邊支點(diǎn)處的單元E損傷!!!simplebeamMFfinish/clear/config,nproc,3!!!!!!!!----------------參數(shù)定義---------------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)!!!邊支點(diǎn)單元20損傷/prep7b=0.3!!!梁截面參數(shù)h=0.5et,1,beam3r,1,b*h,b*h**3/12,h!!!area=
2、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!!!邊支點(diǎn)單元損傷，E=ex20為修正量mp,prxy,4,prxybmmp,ex,4,ex20*do,i,1,
3、21,1!!!立21個關(guān)鍵點(diǎn),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處)的材料常數(shù)，第3種材料latt,3,1,1,lsel,s,,,10!!!修改第10單元(L/2處)的材料常數(shù)，第2種材料latt,2,1,1,lsel,s,,,20!!!修改第20單元(支點(diǎn)處)的材料常數(shù)，第4種材料latt,4,1,1,allsellesize,all,1!!!線的尺寸劃分,0.5m一個單元lmesh,al
4、l!!!對梁進(jìn)行網(wǎng)格劃分dk,1,uy,,,,uz,rotx,roty!!!增加約束,簡支梁dk,21,uy,,,,ux,uz,rotx,rotyfinish/solu!模態(tài)求解allselantype,modalacel,,9.8modopt,lanb,8,1,500!頻率提取從1Hz到500Hz且階數(shù)不超過8mxpand,8!計算相應(yīng)頻率的振型solvefinish/post1!!!set,list!!!set,first!!!pldisp,1!!!!!!set,next!!!pldisp,1!!!!!!set,1,5!!!pldisp,1!!!*dim,freq,array,6
5、!!!提取前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*cfo
6、pen,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
