..>ANSYS動力分析(18)-隨機振動分析-實例(1)2010-09-2607:41:23|分類：ANSYS動力分析|標簽：隨機振動實例模型飛機機翼psd|舉報|字號訂閱PSD實例：模型飛機機翼的隨機振動說明：確定由于施加在機翼根部的Y向加速度PSD，在模型飛機機翼中造成的位移和應力。假設機翼在Z=0處固支。操作指南1.去除數據庫并讀入文件wing.inp以創立幾何模型和網格。2.定義材料屬性：3.施加邊界條件。提示：選擇在areas上施加位移約束，拾取Z=0處所有的Areas，約束所有自由度。4.定義新分析為Model，使用BlockLanczos方法，抽取和擴展前15個自然模態。然后求解CurrentLS。5.查看模態形狀，如圖為前4階振型。6.使用所顯示的PSD譜，執行PSDSpectrum分析。首先定義分析類型為Spectrum分析類型為PSD，使用全部模態，計算單元應力：注意激活"Calculateelemstresses〞選項。7.在根底上施加指定的PSD譜(注意：確保PSD的單位是G2/Hz)。施加Y向鼓勵(方法是：在根底節點上施加單位Y向位移)。設置常阻尼比0.02：設置有關參數–重力加速度值注意：響應譜類型選擇Accel(g**2/Hz)，否則后面的PSD譜應該輸入實際加速度值：定義PSD譜表格：繪制PSD譜表格曲線以檢查輸入值：激活和設置模態組合參數：選擇所有模態。設置計算內容：計算各模態的參與因子:在輸出窗口中可以看到參與因子的計算結果：求解：10.在一般后處理中，查看相對位移和應力(載荷步3)。首先查看ResultsSummary：載荷步3為set17查看載荷步3的相對位移和應力：在ReadResults中讀取載荷步3的結果：使用BysetNumber，輸入DataSetNumber為17。查看載荷步3的相對位移：查看載荷步3的SZ應力：查看載荷步4的相對位移–實際上是速度：查看載荷步5的相對位移–實際是加速度：小結：在POST1中，查看載荷步3的相對位移和應力。–-可以直接使用應力云圖，比方SZ，與屈服應力比較嗎？–-對載荷步1怎么樣"–-由1sigma分量應力導出的等效/主應力(equivalent/principalstressesderivedfrom1sigmacomponentstresses)有效嗎"–-載荷步4為速度結果–-載荷步5為加速度結果11.在TimeHistoryPostprocessor中，對翼尖處的一個節點創立Uy的PSD響應。繪制對數-對數曲線圖。–進入時間歷程后處理時，首先出現如下對話框，在其中選擇：Createresponsepower…然后執行命令‘StoreData’并承受默認值。這是計算ResponsePSD所必須的。在定義變量對話框中，點擊add添加變量。選擇翼尖處的一個節點創立Uy變量。在列表中選擇Uy：變量定義完畢：計算變量2的PSD功率譜，存為變量3：顯示變量3的曲線：將圖形改為對數–對數形式，操作如下：設置有關參數：*,y軸的標記、*,y軸的類型(線性或對數)按更改后的參數重新繪制變量圖形：顯示變量2：附錄：輸入文件wing.inp的內容!Modalanalysise*erciseforDynamicsSeminar

!Modelairplanewing

!finish

/clear

/udoc,1,date,offcomplete=0!如果將此句改為complete=1，將會執行完整的分析，否則只是生成模型/filnam,wing

/title,Modalanalysisofamodelairplanewing!Elementtypesandmaterials

/prep7

et,1,82

et,2,95*if,complete,eq,1,then

mp,e*,1,38000

mp,dens,1,1.033e-3/12

mp,pr*y,1,0.3

*endif!Createandmeshthearea

k,1

k,2,2

k,3,2.3,0.2

k,4,1.9,0.45

k,5,1,0.25

l,1,2

l,1,5

bspline,2,3,4,5,,,-1,0,0,-1,-0.25,0

al,all

esize,0.13

amesh,all!E*trudevolumemeshfromthearea

type,2

e*topt,esize,10

e*topt,aclear,1

ve*t,all,,,0,0,10

/view,,1,1,1

eplot

save!Analysistypeandoptions

*if,complete,eq,1,then

/solu

antype,modal

modopt,lanb,15

m*pand,15,,,yes!Constraints

asel,,loc,z,0

da,all,all

asel,all!*solve

finish/solution

/title,RandomVibrationsanalysis(PSD)ofanairplanewing

antype,spectrum

spopt,psd,15,yes

psdunit,1,accg,386.4

psdfrq

psdfrq,1,,20,100,400,600

psdval,1,0.075,0.1,0.1,0.025

asel,s,loc,z,0

da,all,uy,1!applye*citationinUYdirection

asel,all

outre,nsol,all

!computeparticipationfactors

pfact,1,basepsdout,disp,rel

psdout,velo,abs

psdout,acel,abs

psdcom

solve

finish/post1

set,3,1

ntip=node(0,0,10)

uytip=uy(ntip)plns,u,y

finish*vwrite

(/,5*,'UYdisplacementatnodelocated