【文章內(nèi)容簡(jiǎn)介】 重要的?！　? 結(jié)論　　通過對(duì)機(jī)翼根部固支情況的結(jié)構(gòu)進(jìn)行模態(tài)測(cè)量并與有限元模型的固有振動(dòng)模態(tài)數(shù)值計(jì)算結(jié)果進(jìn)行對(duì)比，利用試驗(yàn)所得模態(tài)參數(shù)對(duì)結(jié)構(gòu)有限元模型作出相應(yīng)調(diào)整，利用調(diào)整后的模型進(jìn)行氣動(dòng)彈性顫振分析，得出如下結(jié)論：　　(1)在模態(tài)試驗(yàn)過程中，不同激振方式對(duì)于模態(tài)測(cè)試結(jié)果有一定的影響，以激振器激勵(lì)所得結(jié)果較有參考價(jià)值。雖然激振桿對(duì)結(jié)構(gòu)的附加剛度的影響使得各階模態(tài)頻率普遍比錘擊法測(cè)得頻率高，但是由于結(jié)構(gòu)存在密集模態(tài)，錘擊法因?yàn)樾枰又笖?shù)窗的原因使模態(tài)識(shí)別結(jié)果的正確性不能保證，而激振器激勵(lì)的方式通過采用合適的激勵(lì)信號(hào)可以很順利地進(jìn)行測(cè)量?！　?2)根據(jù)模態(tài)試驗(yàn)結(jié)果對(duì)計(jì)算模型進(jìn)行修正會(huì)一定程度地影響顫振計(jì)算結(jié)果。本文機(jī)翼模型的顫振形式屬于典型的彎扭耦合顫振：二階彎曲與一階扭轉(zhuǎn)的耦合，且修正前后耦合形式不變?？梢姼鶕?jù)模態(tài)試驗(yàn)結(jié)果對(duì)計(jì)算模型的修正是十分有必要的。（附錄．結(jié)構(gòu)動(dòng)力學(xué)計(jì)算程序）計(jì)算結(jié)構(gòu)力學(xué)編程大作業(yè)時(shí)間： 2007年6月 !!!**************************************************************************** !!! 關(guān)于程序的說明 !!!**************************************************************************** !一、功能： ! 可計(jì)算包括節(jié)點(diǎn)力，一般非節(jié)點(diǎn)力，支座沉降、溫度荷載作用、制造誤差的平 ! 面桁架、梁、剛架及其組合結(jié)構(gòu)的節(jié)點(diǎn)位移與桿端力； ! 可同時(shí)計(jì)算多種工況下的節(jié)點(diǎn)位移與桿端力。 !***************************************************************************** !****************************************************************************** ! ! 二、變量說明： ! NE——單元數(shù)； ! N——結(jié)構(gòu)中自由度數(shù)； ! NJ——節(jié)點(diǎn)數(shù)； ! NS——特殊節(jié)點(diǎn)數(shù)，包括支座節(jié)點(diǎn)、主從節(jié)點(diǎn)（1節(jié)點(diǎn)不做主節(jié)點(diǎn)）、連接桁架的鉸節(jié)點(diǎn)（沒有轉(zhuǎn)角）； ! NAI——結(jié)構(gòu)的單元截面類型數(shù)； ! MT——單元截面類型號(hào)； ! NL——荷載工況數(shù)； ! H——截面高度； ! E——彈性模量； ! JC——單元定位向量數(shù)組； ! X(NJ),Y(NJ)——節(jié)點(diǎn)的X，Y坐標(biāo)值； ! JE(NE,2)——單元兩端節(jié)點(diǎn)碼數(shù)組； ! AI(NAI,2)——按單元類型順序存放A與I，AI（I，1）—第I類單元的截面積，AI（I，2）—第I類單元的 ! 慣性矩； ! MT(NE)——單元所屬單元類型號(hào)； ! JS(NS,4)——特殊節(jié)點(diǎn)信息，JS(I,1)—結(jié)點(diǎn)碼；JS(I,2)，JS(I,3)，JS(I,4)—U，V，CETA約束信息， ! 有約束為1，沒有約束為0；從節(jié)點(diǎn)某位移同主節(jié)點(diǎn)時(shí)位移時(shí)，該位移約束信息填主節(jié)點(diǎn)碼； ! ! PJ(NP,3)——節(jié)點(diǎn)荷載信息數(shù)組；PJ（I，1）—節(jié)點(diǎn)力所在節(jié)點(diǎn)號(hào)；PJ（I，2）—節(jié)點(diǎn)力作用坐標(biāo)方向： ! 坐標(biāo)方向U，V，M分別為1，2，3； PJ（I，3）—節(jié)點(diǎn)力的大小(含正負(fù)號(hào))；U，V方向集中力時(shí)， ! 與坐標(biāo)軸正向同向?yàn)檎?，M按右手法則為正；本程序推導(dǎo)過程取y軸向下為正。 ! ! PF(NF,4)——非節(jié)點(diǎn)荷載數(shù)組，并給出以下類型說明： ! 前6類型數(shù)據(jù)輸法(梯形等可以用疊加法計(jì)算)： ! PF(I,1)單元碼；PF(I,2)類型；PF(I,3)荷載大小；PF(I,4)c值； ! 1——垂直于單元的均布力,大小為q，以坐標(biāo)軸正向?yàn)檎?c為荷載末端距i節(jié)點(diǎn)距離。 ! 2——非節(jié)點(diǎn)集中力P，c為荷載距i節(jié)點(diǎn)距離。 ! 3——非節(jié)點(diǎn)集中力距M，c為荷載距i節(jié)點(diǎn)距離，右手法則判正負(fù)。 ! 4——三角形荷載，c為荷載距i節(jié)點(diǎn)距離，i端為0，距離i端c時(shí)力為q。 ! j端為0的三角形，可按疊加法處理。 ! 5——沿桿軸向均布力,大小為q,c為荷載末端距i節(jié)點(diǎn)距離。 ! 6——沿桿軸向集中力,大小為q,c為荷載末端距i節(jié)點(diǎn)距離。 ! ! 從第7到第9類型（支座沉降）數(shù)據(jù)輸法：PF(I,1)單元碼；PF(I,2)類型；PF(I,3)位移大?。ê?fù)），坐標(biāo)軸正向位為正，轉(zhuǎn)角按右手法則；PF(I,4)沉降所在的單元位移分量，i端為13，j端為46； ! ! 7——沿軸向支座沉降； ! 8——垂直于軸向支座沉降； ! 9——支座轉(zhuǎn)動(dòng)； ! 10——制造誤差,PF(I,1)—制造誤差所在單元，PF(I,2)類型；PF(I,3)誤差大小(含正負(fù)),正負(fù)取決于消除 ! 誤差時(shí)端點(diǎn)的運(yùn)動(dòng)方向，PF(I,4)—誤差所在坐標(biāo)號(hào)； ! 11——溫度荷載，PF(I,1)—荷載所在單元，數(shù)據(jù)形式為： ,如2單元上有溫度荷載，則PF(I,1)=。 ! PF(I,2)—溫度變化值t1, PF(I,3)—溫度變化值t2, PF(I,4)—材料線膨脹系數(shù)。 ! ! TK(NN)——采用一維存儲(chǔ)結(jié)構(gòu)剛度矩陣，上半帶元素（每列第一個(gè)非零元素到對(duì)角元）； ! KD——主元地址數(shù)組，表示結(jié)構(gòu)剛度矩陣的主元在TK中的序號(hào)，KD中最后一個(gè)數(shù)是TK中元素的總個(gè)數(shù)； ! JI——結(jié)構(gòu)剛度矩陣上半帶的非對(duì)角元素在TK中的地址，JI=KD(J)J+I； ! JN(NJ,3)——結(jié)點(diǎn)位移分量編號(hào)數(shù)組，用于存放結(jié)點(diǎn)三個(gè)位移的位移分量號(hào)碼， ! JN(I,1)，JN(I,2)，JN(I,2)分別為結(jié)點(diǎn)I的U，V，CETA分量的位移分量（坐標(biāo)）號(hào)碼； ! ! P(N)——節(jié)點(diǎn)荷載列陣；在回代求位移時(shí)存放位移量； ! F(N)——求得的桿端力列陣； ! FO(6)——等效節(jié)點(diǎn)荷載列陣； ! ! !!!*********************************************************************** !!!********************** 平面結(jié)構(gòu)分析源程序內(nèi)容 ************************* !!!********************************************************************** PROGRAM PFF DIMENSION X(50),Y(50),JE(50,2),MT(30),AI(10,2),JS(20,4),PJ(50,3),PF(50,4),JN(50,3),amp。 KD(150),TK(1000),P(150),F(6),H(50) DOUBLE PRECISION TK,P,F CHARACTER *200 TL OPEN(1,FILE=39。39。,STATUS=39。OLD39。) OPEN(2,FILE=39。39。,STATUS=39。NEW39。) READ(1,70) TL READ(1,70) TL READ(1,*)NE,NJ,NS,NAI,NL,E WRITE(2,10)NE,NJ,NS,NAI,NL,E10 FORMAT(5X,39。PLANE FRAME STRUCTURE ANALYSIS39。/5X,39。**********39。//2X,39。CONTROL PARAMETERS amp。OF STRUCTURE39。/5X,39。NE=39。,I2,8X,39。NJ=39。,I2,8X,39。NS=39。,I2,8X,39。NAI=39。,I2,/5X,39。NL=39。,I2,8X,39。E=39。,) CALL INPUT(NE,NJ,NS,NAI,X,Y,JE,MT,AI,JS,H) !讀入數(shù)據(jù)文件 CALL DJN(NJ,NS,JS,JN,N) !計(jì)算結(jié)構(gòu)自由度數(shù)N，形成結(jié)點(diǎn)位移分量數(shù)組JN CALL ADE(NE,NJ,N,JE,JN,KD,NN) !形成主元地址數(shù)組KD(N) CALL SSM(NE,NJ,NAI,E,N,NN,X,Y,JE,MT,AI,JN,KD,TK) !形成總剛，一維存儲(chǔ)數(shù)組TK(NN) CALL UTDU3(TK,NN,KD,N) !對(duì)總剛進(jìn)行UTDU分解，以用于解方程組 DO 20 LC=1,NL !對(duì)工況循環(huán) READ(1,70)TL READ(1,70)TL READ(1,70)TL READ(1,*)NP,NF !讀入工況信息 WRITE(2,30)LC,NP,NF30 FORMAT(/2X,39。LOAD DATA39。/10X,39。LOAD CASE=39。,I3/10X,39。NP=39。,I3,8X,39。NF=39。,I3) CALL NLV(NE,NJ,NAI,E,N,NP,NF,X,Y,JE,JN,PJ,PF,MT,AI,P,H) !形成總荷載列陣P(N) CALL BACK3(TK,NN,P,N,KD,JN,NJ) !解方程組并輸出結(jié)點(diǎn)位移，存放在數(shù)組P(N)中 WRITE(2,40)40 FORMAT(//4X,39。MEMBEREND FORCES OF ELEMENTS39。/4X,39。ELEMENT39。,13X,39。N39。,17X,39。V39。,17X,39。M39。) DO 60 M=1,NE CALL MQN(M,NE,NJ,NAI,N,NF,E,X,Y,JE,MT,AI,JN,PF,P,F,H) !計(jì)算單元桿端力，存放在數(shù)組F(6)中 WRITE(2,50)M,(F(I),I=1,6) !輸出桿端力50 FORMAT(/1X,I10,3X,39。N1=39。,3X,39。V1=39。,3X,39。M1=39。,39。N2=39。,amp。,3X,39。V2=39。,3X,39。M2=39。,)60 CONTINUE20 CONTINUE70 FORMAT(A) CLOSE(1) CLOSE(2) END SUBROUTINE INPUT(NE,NJ,NS,NAI,X,Y,JE,MT,AI,JS,H) !讀入數(shù)據(jù)文件 DIMENSION X(NJ),Y(NJ),JE(NE,2),MT(NE),AI(NAI,2),JS(NS,4),H(NE) INTEGER NO READ(1,70) TL READ(1,70) TL READ(1,70) TL READ(1,*)(NO,X(I),Y(I),I=1,NJ) READ(1,70) TL READ(1,70) TL READ(1,70) TL READ(1,*)(NO,JE(I,1),JE(I,2),MT(I),H(I),I=1,NE) READ(1,70) TL READ(1,70) TL READ(1,70) TL READ(1,*)(NO,(AI(I,J),J=1,2),I=1,NAI) READ(1,70) TL READ(1,70) TL READ(1,70) TL READ(1,*)((JS(I,J),J=1,4),I=1,NS) WRITE(2,10)(I,X(I),Y(I),I=1,NJ) WRITE(2,20)(I,JE(I,1),JE(I,2),MT(I),I=1,NE) WRITE(2,30)(I,(AI(I,J),J=1,2),I=1,NAI) WRITE(2,40)((JS(I,J),J=1,4),I=1,NS)10 FORMAT(//2X,39。COORDINATES OF JOINTS39。/6X,39。JOINT39。,11X,39。X39。,11X,39。Y39。/(6X,I4,5X,))20 FORMAT(//2X,39。INFORMATION OF ELEMENTS39。/6X,39。ELEMENT39。,amp。 4X,39。JOINTI39。,4X,39。JOINTJ39。,5X,39。TYPE39。/(2X,4I10))30 FORMAT(/7X,39。TYPE39。,10X,39。A39。,12X,39。I39。/(8X,I2,5X,))40 FORMAT(//2X,39。INFORMATION OF SPECIAL JOINTS39。/6X,39。JOINT39。,4X,39。u39。,4x,39。v39。,4x,39。ceta39。/(6X,4I5))