隨機(jī)振動(dòng)分析實(shí)例

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

!Modelairplanewing

!finish

/clear

/udoc,1,date,offplete=0!如果將此句改為plete=1，將會(huì)執(zhí)行完整的分析，否則只是生成模型/filnam,wing

/title,Modalanalysisofamodelairplanewing!Elementtypesandmaterials

/prep7

et,1,82

et,2,95*if,plete,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,plete,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

!puteparticipationfactors

pfact,1,basepsdout,disp,rel

psdout,velo,abs

psdout,acel,abs

psd

solve

finish/post1

set,3,1

ntip=node(0,0,10)

uytip=uy(ntip)plns,u,y

finish*vwrite

(/,5*,'UYdisplacementatnode