安徽建筑工程學院計算結(jié)構(gòu)力學ppt課件

1、整理課件計算結(jié)構(gòu)力學計算結(jié)構(gòu)力學第八章第八章 結(jié)構(gòu)分析結(jié)構(gòu)分析程序設(shè)計程序設(shè)計首頁首頁上頁上頁返回返回下頁下頁整理課件8-1 概述概述: :程序設(shè)計的基本概念與要點程序設(shè)計的基本概念與要點 至此至此, ,我們已完成了結(jié)構(gòu)矩陣分析桿系有我們已完成了結(jié)構(gòu)矩陣分析桿系有限元法基本原理的教學內(nèi)容限元法基本原理的教學內(nèi)容, ,本章主要介紹計本章主要介紹計算機實現(xiàn)過程。算機實現(xiàn)過程。 程序設(shè)計程序設(shè)計: :當人們需要計算機完成科學計當人們需要計算機完成科學計算算, ,數(shù)據(jù)處理等計算工作時數(shù)據(jù)處理等計算工作時, ,必須事先恰當?shù)乇仨毷孪惹‘數(shù)匕才藕糜嬎悴襟E安排好計算步驟, ,每一步的計算內(nèi)容以及在什每一步的

2、計算內(nèi)容以及在什么條件下進行這一部分計算么條件下進行這一部分計算, ,這就是所謂編制這就是所謂編制計算機程序。計算機程序。首頁首頁上頁上頁返回返回下頁下頁整理課件程序設(shè)計應(yīng)注意以下幾點：程序設(shè)計應(yīng)注意以下幾點：=保證程序的正確性，可通過考題校驗保證程序的正確性，可通過考題校驗=使程序具有高效率使程序具有高效率, ,并盡可能降低費用并盡可能降低費用, ,求解求解方程組最費機時方程組最費機時(80%(80%左右左右),),應(yīng)設(shè)計再解功能應(yīng)設(shè)計再解功能=使程序便于調(diào)試、修改、擴充和完善使程序便于調(diào)試、修改、擴充和完善, ,既要既要有通用性有通用性, ,又要留下可擴充修改的窗口又要留下可擴充修改的窗口

3、結(jié)構(gòu)矩陣分析原理結(jié)構(gòu)矩陣分析原理結(jié)構(gòu)分析基礎(chǔ)結(jié)構(gòu)分析基礎(chǔ)數(shù)值計算方法數(shù)值計算方法程序設(shè)計程序設(shè)計算法語言算法語言計算結(jié)果分析計算結(jié)果分析上機調(diào)試上機調(diào)試從結(jié)構(gòu)矩陣分析原理到最終計算機實現(xiàn)解決具從結(jié)構(gòu)矩陣分析原理到最終計算機實現(xiàn)解決具體問題體問題, ,主要有以下過程模塊主要有以下過程模塊, ,用框圖可表示為：用框圖可表示為：首頁首頁上頁上頁返回返回下頁下頁整理課件8-2 結(jié)構(gòu)分析程序設(shè)計的框圖設(shè)計結(jié)構(gòu)分析程序設(shè)計的框圖設(shè)計根據(jù)結(jié)構(gòu)方程組的建立與求解來進行程根據(jù)結(jié)構(gòu)方程組的建立與求解來進行程序設(shè)計序設(shè)計結(jié)構(gòu)方程組的建立與求解是結(jié)構(gòu)分析的核心問結(jié)構(gòu)方程組的建立與求解是結(jié)構(gòu)分析的核心問題題, ,如建立

4、結(jié)構(gòu)方程組的有限元方法如建立結(jié)構(gòu)方程組的有限元方法, ,解線性代解線性代數(shù)方程組的消元分解法數(shù)方程組的消元分解法, ,這兩個方法一經(jīng)確定這兩個方法一經(jīng)確定, ,程序設(shè)計的大致思路與過程也就基本確定了。程序設(shè)計的大致思路與過程也就基本確定了。1. 1. 編程要點編程要點首頁首頁上頁上頁返回返回下頁下頁整理課件由單元定位向量組織整體流程圖的運行由單元定位向量組織整體流程圖的運行實施實施 從形成結(jié)構(gòu)方程組從形成結(jié)構(gòu)方程組K=PK=P過程中過程中KK與與PP的形成的形成, ,到計算結(jié)構(gòu)內(nèi)力和反力的過程到計算結(jié)構(gòu)內(nèi)力和反力的過程, ,都都離不開離不開MWMW的組織。的組織。應(yīng)具有良好的通用性應(yīng)具有良好的

5、通用性 本程序的編制特點是利用特殊結(jié)點的約束信本程序的編制特點是利用特殊結(jié)點的約束信息，可模擬六種類型的桿系結(jié)構(gòu)，故具有廣泛息，可模擬六種類型的桿系結(jié)構(gòu)，故具有廣泛的適用性。的適用性。首頁首頁上頁上頁返回返回下頁下頁整理課件2. 設(shè)計流程圖設(shè)計流程圖主程序主程序,數(shù)據(jù)輸入數(shù)據(jù)輸入形成總剛形成總剛K對對K消元分解消元分解形成荷載列陣形成荷載列陣P對對P正消，回代正消，回代求出結(jié)點位移求出結(jié)點位移求內(nèi)力與反力求內(nèi)力與反力 形成形成MW 形成單剛形成單剛K1234567對對各各組組數(shù)數(shù)據(jù)據(jù)循循環(huán)環(huán)對對各各組組荷荷載載循循環(huán)環(huán)首頁首頁上頁上頁返回返回下頁下頁整理課件主、數(shù)據(jù)輸入主、數(shù)據(jù)輸入KJX:形成

6、形成KNXFJ:消元分解消元分解YDX:形成形成 PNXFJ:回代求解回代求解YWY:打印位移打印位移QDL:求求F與與RQJW21222324252627 結(jié)束結(jié)束 開始開始 QJZHQCHQKDQXSDKXPGPDXL1112131415161718ZERO1ZERO2JZZQMW010203043.框框圖圖設(shè)設(shè)計計首頁首頁上頁上頁返回返回下頁下頁整理課件4. 4. 本程序設(shè)計的模塊功能介紹本程序設(shè)計的模塊功能介紹模塊模塊：具有單一的獨立的功能塊：具有單一的獨立的功能塊, ,由子由子程序或自定義函數(shù)所組成。程序或自定義函數(shù)所組成。模塊可根據(jù)其功能進一步劃分模塊可根據(jù)其功能進一步劃分, ,依

7、次分解依次分解成較低級的模塊成較低級的模塊, ,模塊之間通過調(diào)用而組模塊之間通過調(diào)用而組成一個協(xié)同的程序；成一個協(xié)同的程序；這種可通過自上而這種可通過自上而下進行分解下進行分解, ,并可通過自上而下的調(diào)用并可通過自上而下的調(diào)用, ,一級一級地組成程序是程序設(shè)計的重要一級一級地組成程序是程序設(shè)計的重要方法。方法。首頁首頁上頁上頁返回返回下頁下頁整理課件各個模塊的功能在很大程度上是獨立的各個模塊的功能在很大程度上是獨立的, ,因而不同的模塊可以由不同的人來完成。因而不同的模塊可以由不同的人來完成。例如例如, ,不太了解矩陣力學的人也可以設(shè)計不太了解矩陣力學的人也可以設(shè)計消元分解及正消回代的子程序消

8、元分解及正消回代的子程序SUB. NXFJSUB. NXFJ。模塊的相互獨立性不僅方便程序設(shè)計模塊的相互獨立性不僅方便程序設(shè)計, ,也也方便程序調(diào)試。方便程序調(diào)試。調(diào)試時自下而上一塊一塊調(diào)試時自下而上一塊一塊地進行。這時主要著重調(diào)試模塊接口和上地進行。這時主要著重調(diào)試模塊接口和上一級模塊一級模塊, ,而不必調(diào)試已通過的模塊。而不必調(diào)試已通過的模塊。亦可以根據(jù)模塊功能組成其它程序。亦可以根據(jù)模塊功能組成其它程序。首頁首頁上頁上頁返回返回下頁下頁整理課件本程序的模塊設(shè)有三個級別本程序的模塊設(shè)有三個級別:0:0、1 1、2 2=0 0級表示各個程序均可調(diào)用級表示各個程序均可調(diào)用=1 1級僅供級僅供2

9、 2級調(diào)用。級調(diào)用。首頁首頁上頁上頁返回返回下頁下頁整理課件5. 5. 本程序的靜力計算功能本程序的靜力計算功能(1)(1)結(jié)構(gòu)形式結(jié)構(gòu)形式: :可對六種類型的桿系進行計可對六種類型的桿系進行計算分析算分析連續(xù)梁桁架排架剛架框排架連續(xù)梁桁架排架剛架框排架梁、桁組合結(jié)構(gòu)?？捎闪?、桁組合結(jié)構(gòu)?？捎蒍TX(4,NJT)JTX(4,NJT)進行分進行分類。類。(2)(2)材料材料: :各向同性各向同性, ,按按EI(EA)EI(EA)進行分組進行分組, ,分分組數(shù)為組數(shù)為NAINAI首頁首頁上頁上頁返回返回下頁下頁整理課件(3)(3)荷載類型荷載類型結(jié)點荷載信息結(jié)點荷載信息NPJ:NPJ:受載結(jié)點數(shù)受

10、載結(jié)點數(shù) 需輸入結(jié)點荷載信息數(shù)組需輸入結(jié)點荷載信息數(shù)組PJZ(NPJ,2),NPJPJZ(NPJ,2),NPJ行行,2,2列列: :第一列輸?shù)谝涣休擩D.x,JD.x,第二列輸大小第二列輸大小( (與坐標一致為正，與坐標一致為正，反之為負反之為負) )。JDJD是結(jié)點號，是結(jié)點號，x x是方向是方向(1(1為為X X向向, 2, 2為為Y Y向向, 3, 3為為z z向向) )。首頁首頁上頁上頁返回返回下頁下頁整理課件單元荷載信息單元荷載信息單元荷載信息數(shù)組單元荷載信息數(shù)組PMZ(NPM,3):PMZ(NPM,3):NPM:NPM:受載單元數(shù)受載單元數(shù)NPMNPM行行,3,3列列第一列輸?shù)谝涣?/p>

11、輸M.L:M.L:第第M M單元，第單元，第L L類荷載；類荷載；第二列輸荷載位置第二列輸荷載位置( (距始端距離距始端距離x)x)第三列輸荷載大小第三列輸荷載大小( (與坐標一致為正，反之為負與坐標一致為正，反之為負) )L L類：共六類類：共六類, ,見講義附表見講義附表(P63)(P63)或參考教材表或參考教材表5-1(P131),5-1(P131),程序程序PSTDYPSTDY的子程序的子程序SUB.DXLSUB.DXL中留有中留有用戶入口用戶入口, ,可接入這六種以外的荷載?？山尤脒@六種以外的荷載。首頁首頁上頁上頁返回返回下頁下頁整理課件8-3 單元定位向量的主線作用單元定位向量的主

12、線作用 從程序設(shè)計框圖可以看出從程序設(shè)計框圖可以看出: :程序設(shè)計程序設(shè)計的每一個環(huán)節(jié)都離不開單元定位向量的每一個環(huán)節(jié)都離不開單元定位向量, ,故故它在程序設(shè)計中起到組織者的作用它在程序設(shè)計中起到組織者的作用。我們。我們稱之為主線作用。稱之為主線作用。 回顧回顧: :單元定位向量單元定位向量是按單元結(jié)點編是按單元結(jié)點編號順序由結(jié)點的未知量編號所組成的一個號順序由結(jié)點的未知量編號所組成的一個列向量?？捎闪邢蛄??？捎蒍W(3,NJ)JW(3,NJ)直接生成直接生成, ,其作用其作用主要有：主要有：首頁首頁上頁上頁返回返回下頁下頁整理課件1.1.裝備結(jié)構(gòu)剛度矩陣裝備結(jié)構(gòu)剛度矩陣（1 1）按方陣存貯）

13、按方陣存貯SUB.KJX1SUB.KJX1DO 10 I=1,6DO 10 I=1,6L=MW(I)L=MW(I)IF (L.LE.0) GOTO 10IF (L.LE.0) GOTO 10DO 20 J=1,6DO 20 J=1,6K=MW(J)K=MW(J)IF (K.LE.0) GOTO 20IF (K.LE.0) GOTO 20ZK(L,K)=ZK(L,K)+DK(I,J)ZK(L,K)=ZK(L,K)+DK(I,J)2020CONTINUECONTINUE1010CONTINUECONTINUE首頁首頁上頁上頁返回返回下頁下頁整理課件思考題：若形成上三角陣如何改動思考題：若形成上三角

14、陣如何改動（2 2）一維變帶寬上三角按行存貯）一維變帶寬上三角按行存貯SUB.KJXSUB.KJXDO 10 I=1,6DO 10 I=1,6L=MW(I)L=MW(I)IF (L.LE.0) GOTO 10IF (L.LE.0) GOTO 10II=KD(L)II=KD(L)DO 20 J=1,6DO 20 J=1,6K=MW(J)K=MW(J)IF (IF (K.LE.LK.LE.L) GOTO 20) GOTO 20IJ=II-L+KIJ=II-L+KZK(IJ)=ZK(IJ)+DK(I,J)ZK(IJ)=ZK(IJ)+DK(I,J)2020CONTINUECONTINUE1010CON

15、TINUECONTINUE首頁首頁上頁上頁返回返回下頁下頁整理課件2.2.形成形成PP見見PSTDYPSTDY中中SUB.YDXSUB.YDX3.3.單元位移的形成單元位移的形成見見PSTDYPSTDY中中SUB.YWYSUB.YWY由由(存在存在PP中，中，存在存在D(NE,6)D(NE,6)中中4.F4.Fj jRRj j亦由亦由MWMW的第的第j j個分量是否為零來個分量是否為零來判斷判斷, ,見見PSTDYPSTDY中中SUB.QDLSUB.QDL首頁首頁上頁上頁返回返回下頁下頁整理課件5.5.在在KDKD數(shù)組中求帶寬亦用到數(shù)組中求帶寬亦用到MWMW見見PSTDYPSTDY中中SUB.

16、QKDSUB.QKD帶寬公式：帶寬公式：NDK=NDK=單元兩端未知量編號最大差值單元兩端未知量編號最大差值+1+16.6.更為重要的是單元定位向量還體現(xiàn)了單元間更為重要的是單元定位向量還體現(xiàn)了單元間的相互聯(lián)接，以及對結(jié)構(gòu)邊界條件的處理，如的相互聯(lián)接，以及對結(jié)構(gòu)邊界條件的處理，如主從關(guān)系、無效未知量的處理等，可通過對特主從關(guān)系、無效未知量的處理等，可通過對特殊結(jié)點的約束信息數(shù)組殊結(jié)點的約束信息數(shù)組JTX(4,NJT)JTX(4,NJT)來模擬實際來模擬實際結(jié)構(gòu)中的復雜關(guān)系。結(jié)構(gòu)中的復雜關(guān)系。換句話說，我們在程序的換句話說，我們在程序的功能中所提到的六種類型的桿系結(jié)構(gòu)，程序設(shè)功能中所提到的六種類

17、型的桿系結(jié)構(gòu)，程序設(shè)計中最后的區(qū)別形式就是計中最后的區(qū)別形式就是MWMW，這點尤為重要這點尤為重要。首頁首頁上頁上頁返回返回下頁下頁整理課件對各種類型的結(jié)構(gòu)，我們設(shè)計了統(tǒng)一的單元剛對各種類型的結(jié)構(gòu)，我們設(shè)計了統(tǒng)一的單元剛度矩陣形式，由度矩陣形式，由MWMW來直接裝配總剛！來直接裝配總剛！六種類型桿系結(jié)構(gòu)的六種類型桿系結(jié)構(gòu)的JTXJTX數(shù)組數(shù)組1.1.連續(xù)梁連續(xù)梁1234510001 10001 10001 10001100011101010000JTX 12345首頁首頁上頁上頁返回返回下頁下頁整理課件2.2.桁架桁架1234500100010110001 10001 10001 100011

18、0001JTX03.3.排架排架256789100144111040111000111JTX考考慮慮軸軸向向變變形形 12543124356987EI EI EA EA 首頁首頁上頁上頁返回返回下頁下頁整理課件4.4.框架框架1234560100103111000110001100011000111000011JTX不不考考慮慮軸軸向向變變形形 56111111JTX考考慮慮軸軸向向變變形形123465首頁首頁上頁上頁返回返回下頁下頁整理課件5.5.框排架框排架12345678910010011001100105511110001 10001 10001 10001 10001 10001 1

19、0001 10001 10001 100010000000111JTX 不不考考慮慮軸軸向向變變形形348910221112211100111JTX考考慮慮軸軸向向變變形形12346571098首頁首頁上頁上頁返回返回下頁下頁整理課件6.6.框排架框排架124567110013006110013016010001 10001 10001010001JTX2367451,EI EA首頁首頁上頁上頁返回返回下頁下頁整理課件根據(jù)各程序模塊的功能，本程序設(shè)計了下列控根據(jù)各程序模塊的功能，本程序設(shè)計了下列控制變量與循環(huán)變量：制變量與循環(huán)變量： NENE單元總數(shù)單元總數(shù) NJNJ結(jié)點總數(shù)結(jié)點總數(shù) NJTN

20、JT特殊結(jié)點數(shù)特殊結(jié)點數(shù) NJZNJZ支座結(jié)點數(shù)支座結(jié)點數(shù) NAIEANAIEA或或EIEI的分組數(shù)的分組數(shù) NPJNPJ結(jié)點荷載數(shù)結(jié)點荷載數(shù) NPMNPM單元荷載數(shù)單元荷載數(shù)8-4 變量與數(shù)組設(shè)計變量與數(shù)組設(shè)計首頁首頁上頁上頁返回返回下頁下頁整理課件 N N未知量總數(shù)未知量總數(shù) MM單元序號單元序號 LOADLOAD荷載分組序號荷載分組序號 I I、JJ單元剛度矩陣的行列號單元剛度矩陣的行列號 L L，KK結(jié)構(gòu)剛度矩陣的行列號結(jié)構(gòu)剛度矩陣的行列號 IIII，LLLL對角線元素地址對角線元素地址 N1-KDN1-KD數(shù)組的個數(shù)，數(shù)組的個數(shù)， NZYNZY結(jié)構(gòu)剛度矩陣的元素數(shù)結(jié)構(gòu)剛度矩陣的元素數(shù)

21、首頁首頁上頁上頁返回返回下頁下頁整理課件整型數(shù)組有整型數(shù)組有 JH(2JH(2，NE)NE)單元兩端的結(jié)點號單元兩端的結(jié)點號 MW(6)MW(6)單元定位向量單元定位向量 JW(3JW(3，NJ)NJ)結(jié)構(gòu)結(jié)點未知量編號結(jié)構(gòu)結(jié)點未知量編號 JTX(4JTX(4，NJT)NJT)特殊結(jié)點信息，特殊結(jié)點信息， JZH(NJZ)JZH(NJZ)支座維點號支座維點號 NLX(2NLX(2，NL)NL)每組荷載信息，即每組荷載信息，即NPJNPJ數(shù)和數(shù)和NPMNPM數(shù)數(shù) NAI(1NAI(1，NE)NE)截面特性分組號截面特性分組號 KD(NI)KD(NI)結(jié)構(gòu)剛陣一維存貯時主元素地址；結(jié)構(gòu)剛陣一維存貯

22、時主元素地址；首頁首頁上頁上頁返回返回下頁下頁整理課件雙精度實型數(shù)組雙精度實型數(shù)組 X(NJ)X(NJ)結(jié)點坐標值結(jié)點坐標值 Y(NJ)Y(NJ)結(jié)點坐標值結(jié)點坐標值 SL(NE)SL(NE)單元長度單元長度 CX(NE)CX(NE)單元的單元的coscos SY(NE) SY(NE)單元的單元的sinsin EA(NAI) EA(NAI)單元的單元的EAEA EI(NAI) EI(NAI)單元的單元的EIEI XS(7) XS(7)整體坐標系下各單元的整體坐標系下各單元的7 7個常數(shù)，個常數(shù)， XSA(NEXSA(NE，7)7)整體坐標系下各單元的整體坐標系下各單元的7 7個常數(shù)個常數(shù)首頁首

23、頁上頁上頁返回返回下頁下頁整理課件 DK(6 DK(6，6)6)單元剛度矩陣，單元剛度矩陣， T(6,6)T(6,6)坐標變換矩陣，坐標變換矩陣， ZK(NZY)ZK(NZY)結(jié)構(gòu)剛度矩陣，一維存貯結(jié)構(gòu)剛度矩陣，一維存貯 P(N)P(N)荷載列陣荷載列陣 D(NED(NE，6)6)單元兩端的結(jié)點位移，后存結(jié)點力單元兩端的結(jié)點位移，后存結(jié)點力 F(6)F(6)單元結(jié)點力單元結(jié)點力 FE(6)FE(6)等效結(jié)點力等效結(jié)點力 FG(6)FG(6)整體坐標系下的單元結(jié)點力整體坐標系下的單元結(jié)點力 DG(6)DG(6)整體坐標系下單元結(jié)點位移整體坐標系下單元結(jié)點位移 FLZ(NJZFLZ(NJZ，3)3

24、)支座反力支座反力 PMZ(NPMPMZ(NPM，3)3)單元荷載信息單元荷載信息 PJZ(NPJPJZ(NPJ，2)2)結(jié)點荷載信息結(jié)點荷載信息 首頁首頁上頁上頁返回返回下頁下頁整理課件 在采用在采用FORTRANFORTRAN語言編制程序時，一定語言編制程序時，一定要摘清楚數(shù)據(jù)如何傳遞。鑒于要摘清楚數(shù)據(jù)如何傳遞。鑒于FORTRANFORTRAN語言語言的模塊化性質(zhì)，各程序的數(shù)據(jù)一般可通過下的模塊化性質(zhì)，各程序的數(shù)據(jù)一般可通過下列三種方式進行傳遞：列三種方式進行傳遞：啞實結(jié)合啞實結(jié)合 COMMONCOMMON塊塊 數(shù)據(jù)文件數(shù)據(jù)文件方法方法還可運用在機器設(shè)備與外部設(shè)備還可運用在機器設(shè)備與外部設(shè)

25、備( (終終端端) )的數(shù)據(jù)傳遞，的數(shù)據(jù)傳遞，如本程序設(shè)計就采用如本程序設(shè)計就采用2 2個個OPENOPEN語句，建立了輸入、輸出數(shù)據(jù)文件。語句，建立了輸入、輸出數(shù)據(jù)文件。我們在這里介紹的程序設(shè)計我們在這里介紹的程序設(shè)計各子程序間各子程序間的數(shù)的數(shù)據(jù)傳遞均采用方法據(jù)傳遞均采用方法 。8-5 數(shù)據(jù)傳遞與動態(tài)數(shù)組設(shè)計數(shù)據(jù)傳遞與動態(tài)數(shù)組設(shè)計首頁首頁上頁上頁返回返回下頁下頁整理課件子程序的一般形式為子程序的一般形式為SUBROUTINE QJW(NJSUBROUTINE QJW(NJ，NJTNJT，JTXJTX，JWJW，N)N)RETURNRETURNENDEND其中括號內(nèi)的其中括號內(nèi)的NJNJ，N

26、 N即為形式參數(shù)或稱即為形式參數(shù)或稱為為虛擬變量虛擬變量，或稱，或稱啞元啞元。它可以是變量名字、。它可以是變量名字、數(shù)組名字或數(shù)組元素。數(shù)組名字或數(shù)組元素。1 1、啞實結(jié)合的數(shù)據(jù)傳遞方式、啞實結(jié)合的數(shù)據(jù)傳遞方式首頁首頁上頁上頁返回返回下頁下頁整理課件子程序中的形式參數(shù)沒有確切的數(shù)值，這就子程序中的形式參數(shù)沒有確切的數(shù)值，這就是虛擬變量的由來。只有在調(diào)用該子程序時，是虛擬變量的由來。只有在調(diào)用該子程序時，才對形式參數(shù)賦值，或賦予其實在的存貯空才對形式參數(shù)賦值，或賦予其實在的存貯空間。間。如：如：CALL QJW(NJCALL QJW(NJ，NJTNJT，JTXJTX，JWJW，N)N)這時這時N

27、JNJ，N N稱為稱為實在參數(shù)實在參數(shù)( (實元實元) )。注意：。注意：形式參數(shù)和實在參數(shù)的類型應(yīng)一致，個數(shù)應(yīng)形式參數(shù)和實在參數(shù)的類型應(yīng)一致，個數(shù)應(yīng)相等，但參數(shù)名可以不相同相等，但參數(shù)名可以不相同。首頁首頁上頁上頁返回返回下頁下頁整理課件形式參數(shù)形式參數(shù)的作用可分為兩種：的作用可分為兩種：=一種是從主程序或其它子程序來賦值的，也一種是從主程序或其它子程序來賦值的，也就是通過啞實結(jié)合接收從外面輸入的數(shù)據(jù)，作就是通過啞實結(jié)合接收從外面輸入的數(shù)據(jù)，作為本子程序計算的依據(jù)，好象是加工廠的原料為本子程序計算的依據(jù)，好象是加工廠的原料一樣一樣=另一種是本程序模塊計算的結(jié)果，通過啞實另一種是本程序模塊計算

28、的結(jié)果，通過啞實結(jié)合傳遞到調(diào)用處，這是向外傳遞的數(shù)據(jù)，好結(jié)合傳遞到調(diào)用處，這是向外傳遞的數(shù)據(jù)，好象是加工廠的產(chǎn)品一樣。如：象是加工廠的產(chǎn)品一樣。如：子程序：子程序：SUB. QJW(NJSUB. QJW(NJ，NJTNJT，JTXJTX，JWJW，N)N)主程序：主程序：CALL QJW(NJCALL QJW(NJ，NJTNJT，JTXJTX，JWJW，N)N)首頁首頁上頁上頁返回返回下頁下頁整理課件動態(tài)數(shù)組又稱為可調(diào)數(shù)組動態(tài)數(shù)組又稱為可調(diào)數(shù)組，如：，如：JH(2,NE)JH(2,NE)、JW(3,NJ)JTX(4,NJT)JW(3,NJ)JTX(4,NJT)、ZK(NZY)ZK(NZY)等，

29、由于這里等，由于這里NENE、NJNJ、NJTNJT，NZYNZY等均為數(shù)組變量，沒有確切等均為數(shù)組變量，沒有確切的數(shù)字，即數(shù)組的大小未能得到確切的定義，的數(shù)字，即數(shù)組的大小未能得到確切的定義，因而可調(diào)數(shù)組在主程序中是不允許出現(xiàn)的。因而可調(diào)數(shù)組在主程序中是不允許出現(xiàn)的。在主程序中只能出現(xiàn)確切定義的數(shù)組，如我們在主程序中只能出現(xiàn)確切定義的數(shù)組，如我們在前二個大作業(yè)里所介紹的數(shù)組如在前二個大作業(yè)里所介紹的數(shù)組如JH(2,20)JH(2,20)、ZK(50,50)ZK(50,50)等，這樣才能在在等，這樣才能在在DIMENSIONDIMENSION語句中語句中予以確切定義。予以確切定義。2 2、主程

30、序中動態(tài)數(shù)組的設(shè)計、主程序中動態(tài)數(shù)組的設(shè)計首頁首頁上頁上頁返回返回下頁下頁整理課件根據(jù)子程序的形式參數(shù)根據(jù)子程序的形式參數(shù)( (啞元啞元) )的定義，的定義，在子在子程序中可以出現(xiàn)動態(tài)數(shù)組，程序中可以出現(xiàn)動態(tài)數(shù)組，但必須在調(diào)用但必須在調(diào)用時進行啞實結(jié)合，才能進行運算。時進行啞實結(jié)合，才能進行運算。由于在實際計算中大多數(shù)數(shù)組的大小是隨具體由于在實際計算中大多數(shù)數(shù)組的大小是隨具體問題的不同而變化的，但在主程序中又不允許問題的不同而變化的，但在主程序中又不允許出現(xiàn)動態(tài)數(shù)組，這就給我們在出現(xiàn)動態(tài)數(shù)組，這就給我們在DIMENSIONDIMENSION語句語句中如何定義數(shù)組帶來不少困難：既受到計算機中如何

31、定義數(shù)組帶來不少困難：既受到計算機內(nèi)存的限制又應(yīng)使數(shù)組有足夠的存貯空間，這內(nèi)存的限制又應(yīng)使數(shù)組有足夠的存貯空間，這對于在微機上解決大型工程問題顯得更重要。對于在微機上解決大型工程問題顯得更重要。首頁首頁上頁上頁返回返回下頁下頁整理課件我們知道，我們知道，二維或高維數(shù)組在計算機內(nèi)部都是二維或高維數(shù)組在計算機內(nèi)部都是按列存放的，按列存放的，即在計算機內(nèi)部都是按一維數(shù)組即在計算機內(nèi)部都是按一維數(shù)組 的方式來存貯的，這表示數(shù)組變量之間存在一的方式來存貯的，這表示數(shù)組變量之間存在一定的關(guān)系，如對于：定的關(guān)系，如對于：JHJH2 22020數(shù)組：數(shù)組：JH(2JH(2，3)JH(6)3)JH(6)，JH(

32、1JH(1，4)JH(7)4)JH(7)，這說明數(shù)組一經(jīng)定義，這種關(guān)系便確定，計算這說明數(shù)組一經(jīng)定義，這種關(guān)系便確定，計算機立即機立即“了解了解”。此外，一個數(shù)組變量。此外，一個數(shù)組變量( (下標下標變量變量) )的下標值加的下標值加1 1，就是緊跟在它后面的下標，就是緊跟在它后面的下標值，值，這叫做數(shù)組變量這叫做數(shù)組變量( (下標變量下標變量) )的后繼函數(shù)的后繼函數(shù)，計算機的處理功能使程序會自動按后繼函數(shù)找計算機的處理功能使程序會自動按后繼函數(shù)找到下一個元素，到下一個元素，這叫做下標的自動后繼性質(zhì)這叫做下標的自動后繼性質(zhì)。首頁首頁上頁上頁返回返回下頁下頁整理課件如在主程序如在主程序DIME

33、NSIONDIMENSION語句中定義說明了語句中定義說明了JH(2JH(2，20)20)，則在調(diào)用時實元用，則在調(diào)用時實元用JHJH或或JH(1JH(1，1)1)或或JH(1)JH(1)調(diào)用的效果是相同的，調(diào)用的效果是相同的，如在主程序中用一維定如在主程序中用一維定義說明了義說明了JH(40)JH(40)，則用，則用JHJH或或JH(1)JH(1)調(diào)用效果亦調(diào)用效果亦是相同。至于調(diào)用的是相同。至于調(diào)用的JHJH具體形式，則由子程序具體形式，則由子程序?qū)烧{(diào)數(shù)組對可調(diào)數(shù)組JH(2JH(2，NE)NE)進行說明確定，換句話進行說明確定，換句話說，說，利用下標的自動后繼性質(zhì)，采用首元素調(diào)利用下標的

34、自動后繼性質(zhì)，采用首元素調(diào)用，即可在主程序中出現(xiàn)動態(tài)數(shù)組用，即可在主程序中出現(xiàn)動態(tài)數(shù)組。根據(jù)算例的需要，我們在主程序中僅開設(shè)了整、根據(jù)算例的需要，我們在主程序中僅開設(shè)了整、實二個大數(shù)組：實二個大數(shù)組： 整：整：IA(1000)IA(1000)，暫定其大小，暫定其大小NIA=1000NIA=1000 實：實：A(10000)A(10000)，暫定其大小，暫定其大小NA=10000NA=10000首頁首頁上頁上頁返回返回下頁下頁整理課件然后按照程序中所要出現(xiàn)數(shù)組名字順序，定出然后按照程序中所要出現(xiàn)數(shù)組名字順序，定出每個數(shù)組的第每個數(shù)組的第個元素地址，個元素地址，由所輸入的變量由所輸入的變量或巳確定

35、的變量確定這一點并不難，或巳確定的變量確定這一點并不難，參考下圖：參考下圖：整：整：實：實：首頁首頁上頁上頁返回返回下頁下頁整理課件JH(2JH(2，NE)NE)的第一個元素地址為的第一個元素地址為IA(1)IA(1)NLX(2NLX(2，NL)NL)的第一個元素地址為的第一個元素地址為IA(K1)IA(K1)，這里，這里(K1)=1+2(K1)=1+2* *NENEJW(3JW(3，NJ)NJ)的第一個元素地址為的第一個元素地址為IA(K2)IA(K2)，這里，這里(K2)=K1+2(K2)=K1+2* *NLNL對于各雙精度實型數(shù)組在對于各雙精度實型數(shù)組在A A中第一個元素地址中第一個元素

36、地址亦可由圖方便推出。如將某數(shù)組的第一個元素亦可由圖方便推出。如將某數(shù)組的第一個元素作為實元調(diào)用，作為實元調(diào)用，通過啞實結(jié)合過程中的下標自通過啞實結(jié)合過程中的下標自動后繼，當子程序啞元表中相應(yīng)的啞元為動態(tài)動后繼，當子程序啞元表中相應(yīng)的啞元為動態(tài)數(shù)組時，便可完全按照子程序中該數(shù)組的數(shù)組時，便可完全按照子程序中該數(shù)組的DIMENSIONDIMENSION語句的定義在主程序相應(yīng)數(shù)組中得語句的定義在主程序相應(yīng)數(shù)組中得到反映，即在主程序中實現(xiàn)了動態(tài)數(shù)組。到反映，即在主程序中實現(xiàn)了動態(tài)數(shù)組。首頁首頁上頁上頁返回返回下頁下頁整理課件由此可知，只要找出控制各數(shù)組大小的一些變由此可知，只要找出控制各數(shù)組大小的一

37、些變量，即可確定各數(shù)組的量，即可確定各數(shù)組的第一個元素第一個元素。仔細研究。仔細研究這些變量，發(fā)現(xiàn)其中某些變量這些變量，發(fā)現(xiàn)其中某些變量( (如如N N，N1N1，NZY)NZY)等可通過另一些變量由程序計算確定，我們將等可通過另一些變量由程序計算確定，我們將后者稱為后者稱為主控變量主控變量，須按其出現(xiàn)的次序在程序，須按其出現(xiàn)的次序在程序中首先輸入。本程序的主控變量為：中首先輸入。本程序的主控變量為：NENE、NJNJ、NJTNJT、NJZNJZ、NLNL、NAINAI。為促進同學們今后的工。為促進同學們今后的工程應(yīng)用能力，我們這里所提供的程序程應(yīng)用能力，我們這里所提供的程序PSTDYPSTD

38、Y并并有動力分析內(nèi)容，故主控變是還增加一個有動力分析內(nèi)容，故主控變是還增加一個MJMJ。首頁首頁上頁上頁返回返回下頁下頁整理課件MJMJ：擬求振型數(shù)。：擬求振型數(shù)。 在在PSTDYPSTDY中，我們將中，我們將MJMJ采用屏幕輸入，并采用屏幕輸入，并兼作靜、動力分析的程序運行控制，在屏幕的兼作靜、動力分析的程序運行控制，在屏幕的提示下，若鍵人提示下，若鍵人 MJ=0MJ=0僅作靜力計算僅作靜力計算 MJ=3MJ=3做完靜力計算后，再求前三階頻做完靜力計算后，再求前三階頻率和振型，具體請參閱算例率和振型，具體請參閱算例1 1，不另贅述。，不另贅述。 另外，在主程序中還備有這兩大類型數(shù)組另外，在主

39、程序中還備有這兩大類型數(shù)組實際使用的元素個數(shù)顯示，如實際使用的元素個數(shù)顯示，如發(fā)生發(fā)生IAIA或或A A的溢的溢出現(xiàn)象，只需在內(nèi)存允許的范圍內(nèi)修改這四條出現(xiàn)象，只需在內(nèi)存允許的范圍內(nèi)修改這四條語句，見主程序的語句，見主程序的001001、002002、004004、005005句。句。首頁首頁上頁上頁返回返回下頁下頁整理課件1 1、根據(jù)設(shè)計框圖，各子程序及其功能匯總?cè)缦拢?、根?jù)設(shè)計框圖，各子程序及其功能匯總?cè)缦拢?2 KJX22 KJX：形成結(jié)構(gòu)剛陣：形成結(jié)構(gòu)剛陣ZK(NZY)ZK(NZY)2323、25 25 NXFJNXFJ：消元分解法解線代方程組：消元分解法解線代方程組2323：剛陣：剛

40、陣ZKZK的消元的消元2525：對：對PP的正消回代的正消回代24 YDX24 YDX：形成：形成PP26 YWY26 YWY：打印桿端位移：打印桿端位移27 QDL27 QDL：計算：計算FF及及RR8-6 源程序設(shè)計源程序設(shè)計首頁首頁上頁上頁返回返回下頁下頁整理課件28(28(補充補充) MJX) MJX：形成結(jié)構(gòu)質(zhì)量的矩陣：形成結(jié)構(gòu)質(zhì)量的矩陣29(29(補充補充) DYNA) DYNA：逆迭代法計算結(jié)構(gòu)前：逆迭代法計算結(jié)構(gòu)前MJMJ階振型階振型和頻率和頻率11 QJW11 QJW：結(jié)點未知量編號數(shù)組：結(jié)點未知量編號數(shù)組12 QJZH12 QJZH：形成支座結(jié)點號：形成支座結(jié)點號13 DC

41、H13 DCH：形成單元常數(shù)：形成單元常數(shù)14 QKD14 QKD：形成主元地址數(shù)組：形成主元地址數(shù)組15 QXS15 QXS：形成單剛系數(shù)數(shù)組：形成單剛系數(shù)數(shù)組16 DKX16 DKX：形成單剛矩陣：形成單剛矩陣17 PGP17 PGP：疊加荷載，形成：疊加荷載，形成PP的具體計算，的具體計算，首頁首頁上頁上頁返回返回下頁下頁整理課件18 DXL18 DXL：形成：形成FFE E 01 ZER0101 ZER01：對向量充零：對向量充零02 ZER0202 ZER02：對二維數(shù)組充零：對二維數(shù)組充零03 JZZ03 JZZ：矩陣轉(zhuǎn)置：矩陣轉(zhuǎn)置04 QMW04 QMW：形成單元定位向量。：形成

42、單元定位向量。這里的順序由二位數(shù)表示，十位數(shù)上有這里的順序由二位數(shù)表示，十位數(shù)上有0 0，1 1，2 2級，其中高級別的可調(diào)用低級別的模塊，個位數(shù)級，其中高級別的可調(diào)用低級別的模塊，個位數(shù)則表示被調(diào)用的順序。則表示被調(diào)用的順序。首頁首頁上頁上頁返回返回下頁下頁整理課件由由 OPEN(1OPEN(1，F(xiàn)ILE=QAZ.TXT)FILE=QAZ.TXT)建立數(shù)據(jù)文件建立數(shù)據(jù)文件 READ(1READ(1，* *) NE) NE，NJNJ，NJT,NJZNJT,NJZ，NLNL，NAINAI( (見主程序見主程序語句標號語句標號003) 003) READ(1READ(1，* *) JH) JH，N

43、LXNLX，JMHJMH，JTXJTX( (見主程序語句標號見主程序語句標號006)006)READ(1READ(1，* *) EA) EA，EIEI，X X，Y Y( (見程序語句標號見程序語句標號007)007)READ(1READ(1，* *) PJZ(NPJ) PJZ(NPJ，2)2)( (見見SUBSUBYDXYDX語句標號語句標號008) 008) READ(1READ(1，* *) PMZ(NPM) PMZ(NPM，3)3)( (見見SUBSUBYDXYDX浯句標號浯句標號009)009)2、數(shù)據(jù)文件的形成、數(shù)據(jù)文件的形成首頁首頁上頁上頁返回返回下頁下頁整理課件在在PSTDYPS

44、TDY程序中，我們用程序中，我們用OPENOPEN語句建立了數(shù)據(jù)輸語句建立了數(shù)據(jù)輸入文件入文件QAEQAETXTTXT及結(jié)果輸出文件及結(jié)果輸出文件 FCADFCADTXTTXT。全。全部的原始數(shù)據(jù)，除擬求振型數(shù)部的原始數(shù)據(jù)，除擬求振型數(shù)MJMJ兼作靜、動力計兼作靜、動力計算的運行控制從屏幕輸入外，均可在算的運行控制從屏幕輸入外，均可在QAEQAETXTTXT文文件上從終端讀人。件上從終端讀人。全部的數(shù)據(jù)均采用自由格式，由全部的數(shù)據(jù)均采用自由格式，由READ(1READ(1，* *) )語句語句讀人，這里讀人，這里“1 1”表示通道號，表示通道號，“* *”表示自由格表示自由格式。自由格式要求兩

45、個數(shù)之間用式。自由格式要求兩個數(shù)之間用“，分隔，分隔，且每個數(shù)所占的位數(shù)不限，但要與對應(yīng)變量的隱且每個數(shù)所占的位數(shù)不限，但要與對應(yīng)變量的隱含數(shù)型一致。含數(shù)型一致。3. 算例算例首頁首頁上頁上頁返回返回下頁下頁整理課件1 1數(shù)據(jù)文件可在每行開始直接輸入數(shù)據(jù)，不用數(shù)據(jù)文件可在每行開始直接輸入數(shù)據(jù)，不用說明符號說明符號“C”C”，也沒有語句標號區(qū)及續(xù)行區(qū)，也沒有語句標號區(qū)及續(xù)行區(qū)，但每行不應(yīng)超過但每行不應(yīng)超過7272列列2 2每個輸入語句的結(jié)尾不要加任何標點符號，每個輸入語句的結(jié)尾不要加任何標點符號，但實型數(shù)據(jù)的小數(shù)點除外但實型數(shù)據(jù)的小數(shù)點除外3 3下一個語句不要與上一個輸入語句接著輸入，下一個語句

46、不要與上一個輸入語句接著輸入，而應(yīng)別起一行而應(yīng)別起一行建立數(shù)據(jù)文件時還要注意以下幾點建立數(shù)據(jù)文件時還要注意以下幾點首頁首頁上頁上頁返回返回下頁下頁整理課件算例算例1 1作圖示剛架的靜力計算，并求出其基本作圖示剛架的靜力計算，并求出其基本頻率與振型。頻率與振型。任取任取l=4m,m=3.0=4m,m=3.010102 2kg/mkg/mP=2KN,E=2.1P=2KN,E=2.110106 6KN/mKN/m2 2I=4I=410105 5cmcm4 4=4=41010-3-3m m4 4EI=8.4EI=8.410103 3KNmKNm2 2則則各結(jié)點標號與單元劃分如圖示。各結(jié)點標號與單元劃分

47、如圖示。1243h2l1.5m4EIEIEImmP首頁首頁上頁上頁返回返回下頁下頁整理課件1 1第一個輸入語句在主程序第一個輸入語句在主程序003003句，要求輸入控句，要求輸入控制變量制變量6 6個：個： NENE， NJNJ， NJTNJT， NJZNJZ， NLNL， NAINAI 3 3， 4 4， 4 4， 2 2， 1 1， 2 2注意，這六個整形數(shù)成為一行，句尾不加任何符注意，這六個整形數(shù)成為一行，句尾不加任何符號。號。首頁首頁上頁上頁返回返回下頁下頁整理課件2 2第二個輸入語句在主程序第二個輸入語句在主程序006006句，要求輸入整句，要求輸入整型數(shù)組型數(shù)組IAIA，實際上輸入

48、，實際上輸入JHJH，NLXNLX，JMHJMH，JTXJTX，其，其中中JHJH有有2 2* *NE=6NE=6個數(shù)，個數(shù)，NLXNLX有有2 2* *NL=2NL=2個數(shù)，個數(shù)，JMHEJMHE有有NE=3NE=3個數(shù)，個數(shù)，JTXJTX有有4 4* *NJT=16NJT=16個數(shù)，總共個數(shù)，總共2929個整型個整型數(shù)，可分為兩行，但這兩行之間要用逗號數(shù)，可分為兩行，但這兩行之間要用逗號“，”隔開。隔開。3,1,4,2,1,2,1,0,1,1,2,1,0,10001,0,2,1001,13,1,4,2,1,2,1,0,1,1,2,1,0,10001,0,2,1001,10001,0,3,1

49、,1,1,4,1,1,1 0001,0,3,1,1,1,4,1,1,1 注意：這里特殊節(jié)點約束信息注意：這里特殊節(jié)點約束信息JTX(4JTX(4，NJT)NJT)的填的填寫是一項非常細致的工作，需根據(jù)具體力學模型，寫是一項非常細致的工作，需根據(jù)具體力學模型，參照參照1717中的有關(guān)說明，認真填寫。中的有關(guān)說明，認真填寫。首頁首頁上頁上頁返回返回下頁下頁整理課件3 3第第3 3個輸入語句在主程序個輸入語句在主程序007007句，要求輸入雙句，要求輸入雙精度型數(shù)組精度型數(shù)組A A，實際上輸入，實際上輸入EAEA，EIEI，X X，y y。亦采。亦采用自由格式，用自由格式，EAEA有有2 2個數(shù)，個

50、數(shù)，EIEI也有也有2 2個數(shù)，個數(shù)，X X有有4 4個個數(shù)，數(shù)，Y Y也有也有4 4個數(shù)，總共有個數(shù)，總共有1212個數(shù)，可一行輸入：個數(shù)，可一行輸入：0,0,8400,33600,0,8,0,8,4,4,0,00,0,8400,33600,0,8,0,8,4,4,0,04 4在主程序?qū)奢d的循環(huán)中，要輸入每一組荷在主程序?qū)奢d的循環(huán)中，要輸入每一組荷載的數(shù)據(jù)，載的數(shù)據(jù)，PSTDYPSTDY程序仍采用通常的荷載信息集程序仍采用通常的荷載信息集約方法，在集中荷載作用時，要給出節(jié)點荷載的約方法，在集中荷載作用時，要給出節(jié)點荷載的個數(shù)，荷載作用在哪號節(jié)點的哪個位移方向以及個數(shù)，荷載作用在哪號節(jié)點的

51、哪個位移方向以及荷載的大小，這些信息可以定義一個數(shù)組荷載的大小，這些信息可以定義一個數(shù)組PJZ(NPJPJZ(NPJ，2)2)來存放。來存放。首頁首頁上頁上頁返回返回下頁下頁整理課件PJZ(IPJZ(I，1)1)可填實型數(shù)可填實型數(shù)JD.xJD.x其中：其中：JDJD一荷載作用的節(jié)點號；一荷載作用的節(jié)點號； x x一荷載的作用方向，可在總體坐標系下一荷載的作用方向，可在總體坐標系下考慮?？紤]。x=1x=1，荷載作用沿，荷載作用沿X X方向；方向；x=2x=2，荷載作用，荷載作用沿沿Y Y方向方向,x=3,x=3，繞，繞Z Z軸正向作用的力矩。軸正向作用的力矩。 PJZ(1PJZ(1，2)2)填

52、荷載的大小，與坐標方向一致填荷載的大小，與坐標方向一致者為正，者為正，I I表示節(jié)點荷載序號。表示節(jié)點荷載序號。 在單元荷載作用下，在單元荷載作用下，PSTDYPSTDY的靜力計算可解的靜力計算可解決六種類型的單元荷載問題，我們只要給出如下決六種類型的單元荷載問題，我們只要給出如下信息：信息： 首頁首頁上頁上頁返回返回下頁下頁整理課件(1)(1)在哪個單元在哪個單元M M上有什么類型上有什么類型L L的荷載作用；的荷載作用；(2)(2)在右手系下荷載離單元左節(jié)點的距離；在右手系下荷載離單元左節(jié)點的距離；(3)(3)荷載大小荷載大小Q Q，其符號規(guī)定與坐標方向一致為正。，其符號規(guī)定與坐標方向一致

53、為正。于是也可定義一個數(shù)組于是也可定義一個數(shù)組PMZ(NPMPMZ(NPM，3)3)來存放上述來存放上述信息。信息。 PMZ(IPMZ(I，1)1)填填M.LM.L； PMZ(IPMZ(I，2)2)填填x x； PMZ(IPMZ(I，3)3)填填Q Q；首頁首頁上頁上頁返回返回下頁下頁整理課件其中，其中，L=1L=1，表示左端有局部均布荷載，表示左端有局部均布荷載Q Q作用；作用；L=2L=2，表示離左端，表示離左端x x處有集中荷載處有集中荷載Q Q作用；作用；L=3L=3，表示離左端，表示離左端x x處有集中力矩處有集中力矩Q Q作用；作用；L=4L=4，表示左端有局部三角形荷載作用，其，

54、表示左端有局部三角形荷載作用，其x x處為處為 Q QL=5L=5，表示左端有均布軸力荷載，表示左端有均布軸力荷載Q Q作用；作用；L=6L=6，表示離左端，表示離左端x x處有集中軸力處有集中軸力Q Q作用；作用；首頁首頁上頁上頁返回返回下頁下頁整理課件在子程序在子程序YDXYDX中，有兩個輸入語句，應(yīng)分別輸入中，有兩個輸入語句，應(yīng)分別輸入PJZPJZ與與PMZPMZ。本例只有一組結(jié)點荷載。本例只有一組結(jié)點荷載NL=1NL=1，沒有單，沒有單元荷載元荷載NPM=0NPM=0，結(jié)點荷載數(shù)，結(jié)點荷載數(shù)NPJ=1NPJ=1，所以第四個輸，所以第四個輸入語句在子程序入語句在子程序SUBSUBYDXY

55、DX中標號為中標號為008008句，輸入句，輸入結(jié)點荷載信息結(jié)點荷載信息PJZPJZ，共，共2 2* *NPJ=2NPJ=2個數(shù)，也占個數(shù)，也占行：行： 1.11.1，20002000注意，由于注意，由于NPM=0NPM=0，單元荷載輸入信息，單元荷載輸入信息PME(NPMPME(NPM，3)3)的輸入語句的輸入語句009009，程序自動跳轉(zhuǎn)，該語句可以，程序自動跳轉(zhuǎn)，該語句可以不輸入。不輸入。首頁首頁上頁上頁返回返回下頁下頁整理課件5 5至此，靜力計算的所有信息輸入完畢，如不至此，靜力計算的所有信息輸入完畢，如不做動力計算，可直接從屏幕上輸做動力計算，可直接從屏幕上輸“0”0”，程序運，程序

56、運行結(jié)束，并輸出結(jié)果。如需要進行動力計算時，行結(jié)束，并輸出結(jié)果。如需要進行動力計算時，應(yīng)輸入所求振型數(shù)應(yīng)輸入所求振型數(shù)MJMJ，從屏幕上輸入。本例，從屏幕上輸入。本例MJ=1MJ=1，可直接輸入可直接輸入“1”1”。當。當MJ0MJ0時，還應(yīng)在時，還應(yīng)在QAZQAZTXTTXT文件上繼續(xù)輸入動力計算的有關(guān)數(shù)據(jù)，這樣，第文件上繼續(xù)輸入動力計算的有關(guān)數(shù)據(jù)，這樣，第五個輸入語句在子程序五個輸入語句在子程序SUBSUBMJXMJX的的010010句，輸入句，輸入ERM(NAI)ERM(NAI)桿件質(zhì)量線密度數(shù)據(jù)，本題有桿件質(zhì)量線密度數(shù)據(jù)，本題有NAI=2NAI=2個個數(shù)，也占一行：數(shù)，也占一行： 300300，450450首頁首頁上頁上頁返回返回下頁下頁整理課件解解 根據(jù)圖根據(jù)圖5 55 5，按程序中輸入語句的順序依次，按程序中輸入語句的順序依次填寫數(shù)據(jù)。由算例，可知：填寫數(shù)據(jù)。由算例，可知： 1 1第一個輸入語句標號第一個輸入語句標號003003，要求輸人控制變量，要求輸人控制變量6 6個：個：NENE，NJNJ，NJTNJT，NJZNJZ，NLNL，NAINAI 9 9， 8 8， 2 2， 2 2， 1 1， 2 2這這6 6個數(shù)恰好一行，在文件中叫做一個記錄。個數(shù)恰好一行，