【正文】 T OF INERTIA 1 DATA OF ELEMENT 1 2 1 2 3 1 4 3 1 DATA OF SUPPORT 1 0 0 0 4 0 0 0 NODAL COORDINATE N0. NODE XX ZZ 1 2 3 4 TOTAL DEGREES OF FREEDOMNX= 6 NUMBER [K] ELEMENT NH= 21 WORKING CASE= 1 NUMBERS OF NODAL LOAD= 1 NODAL NUMBER OF LOAD= 2 VALUE OF NODAL LOAD XFORCE ZFORCE YMOMEMT NUMBER OF ELEMENTAL LOAD= 2 ELEMENTAL NUMBER OF LOAD 1 2 VAIUE OF FIXING FORCE ELEMEMT INX IQZ IMM JNX JQZ J=MM 1 2 NODAL EQUIVALENT LOAD NODAL DISPLACEMENT ELEMENT XX ZZ ANGLE 1 2 3 4 MEMBER END FORCE AXIAL FORCE SHEARING FORCE BENDING MOMEMT NE= 1 FMJ= FMK= NE= 2 FMJ= FMK= NE= 3 FMJ= FMK= SUPPOT REACTION NUMBER OF NODE XFORCE ZFORCE BENDING MOMEMT 1 4 *************************************************************************本例題在輸入數(shù)據(jù)結(jié)束后，形成的的Fortran執(zhí)行文件源代碼如下：STRUCTURE DEFINE FILE(BY BUILDER )NODE,41,0,02,0,43,44,0ELEMENT,31,1,2,12,2,3,13,4,3,1RESTRAINT,21,1,1,14,1,1,1MATERIAL,11,2000000,.05,.005,.01NODE_CENTRE_LOAD,12,20,0,0ELEMENT_CENTRE_LOAD,12,0,40,0ELEMENT_DISTRIBUTE_LOAD,11,0,4,30,0FIXTURE_FORCE,0，點擊計算。學會矩陣位移法的原理和程序化的實現(xiàn)，主要內(nèi)容有：基本信息的介紹、剛度矩陣到總剛度矩陣、邊界條件的處理、非節(jié)點荷載的等效轉(zhuǎn)換、建立剛度方程、計算桿端力的實現(xiàn)。 通過用對所開發(fā)的基于VB ，檢驗了程序的可用性和良好的可視化性能。[6] Basic 編程及實例分析教程[M].北京：清華大學出版社，2007.[7] 徐春林，趙明誠，.[8] 龍馭球，（H）[M].北京：高等教育出版社，2001。只是對簡單的平面鋼架做了分析，對于復雜的結(jié)構(gòu)還需要結(jié)合有限元的程序進行開發(fā)。學習矩陣位移法Fortran編程語言，與VB ，明白了怎樣把矩陣位移法Fortran程序編輯成VB 。從介紹矩陣位移法的原理入手，介紹了矩陣位移法程序化的實現(xiàn)方法，同時提供了比較典型的例題分析。（4） 例題一平面剛架如圖11（ａ）所示，剛架各桿尺寸及所受荷載如圖上所示，各桿的彈性模量＝橫截面面積＝、慣性矩＝，求內(nèi)力及支座反力。該系統(tǒng)通過上述方法，不僅解決了數(shù)據(jù)可視化的現(xiàn)實問題，也實現(xiàn)了Fortran和VB之間無縫連接，為系統(tǒng)的后處理創(chuàng)造了條件。Fortran源文件的生成采用文件操作結(jié)合print過程和文本框來實現(xiàn)。聲明Openprocess和C1oseH andle或Openprocess和GetExitCodeProcess（2） 黑屏問題的解決。W in32API提供的進程操作函數(shù)可解決上述外殼程序的異步運行問題。輸出的數(shù)據(jù)進行處理，使之形成自動形成表格。VB ，交換海量數(shù)據(jù)非常經(jīng)濟，效率也非常高。 基于VB 混合兩種語言Fortran和VB開發(fā)科學計算軟件，一般有兩種方法：一是用VB界面設(shè)計和控制程序， Fortran程序被編譯為動態(tài)鏈接庫，由VB程序調(diào)用它。但是它們還是有一定區(qū)別的, 具體比較如表 2 所示。兩個文本框的作用是為了顯示返回的結(jié)果, 以檢驗程序是否正確執(zhí)行。Fortran 子程序代碼:SUBROUT INE ARRAYT EST(n,arr)! ＄AT T RIBUTES DLLEXPORT: : ARRAYTEST INT EGER n,i,j REAL( 4) arr( 3, 3) DO i= 1, 3 DO j= 1, 3 arr(i,j) = 10. 1* i+ j END DO END DO n= n+ 10 END SUBROUT INE編譯生成 exam. DLL 文件, 將文件拷到 C:WINDOWSSYST EM 子 目錄下。aliasname語法及組成部分如下:subnamen其中, n 為棧的大小, 它是 4 的倍數(shù)。具體的語法是:　 Public Declare Sub subname lib “l(fā)ibname”Alias “aliasname”( arglist) Declare語句的各部分說明如下: Public: 用于聲明在所有模塊中對所有其他過程均合法的過程。對于多維數(shù)組, VB和Fortran都是按列傳遞的。 而且同時要注意 Fortran 和 VB , 具體見表 1。: 缺省約定、 DCALL調(diào)用約定。還有注意的是,必須明確外部例程的實現(xiàn)機制、調(diào)用約定（調(diào)用約定就是明確程序如何調(diào)用過程, 傳遞過程以及命名規(guī)范）、例程名和參數(shù)列表。（2）調(diào)用 Fortran 動態(tài)鏈接庫 此方法是將Fortran的源代碼分成若干個子程序或函數(shù), 文件。Shell 函數(shù)語法如下:Shell(pathname, [window style] )pathname 是可執(zhí)行文件的程序名, 含有所必需的目錄、文件夾、驅(qū)動器。在計算結(jié)束后, 將計算數(shù)據(jù)以數(shù)字或者圖形的形式顯示在 VB 的界面上。把已知條件按設(shè)定的格式輸入就能運算出問題的解。這時方程的階數(shù)變?yōu)閞= (3) 通過這方程解出未知位移,再加上已知位移就得到了所有結(jié)點位移。為坐標變換矩陣的轉(zhuǎn)置矩陣。對于析架、剛架等，各單元軸線方向各異。因為拉(壓)與扭轉(zhuǎn)的單元剛度矩陣形式相同，這樣就只有兩種單元剛度矩陣。為結(jié)點位移列陣。彎扭組合s=3。n= ,s為結(jié)點位移的自由度數(shù)。公式拓展為：+這一步利用矩陣乘法進行運算。 建立求解剛度方程各結(jié)點的位移是剛度方程中的未知量。這樣處理后，不論桿件兩端的約束情況如何，統(tǒng)一看作固定端，用位移法計算桿端支座反力，再將反力反向加在桿端，作為等效結(jié)點荷載。本設(shè)計在處理界面時更簡化，直接在VB ，分為三列：約束號、結(jié)點號、約束類型（包括固定支座、鉸支座、水平輥軸支座、豎直輥軸支座、水平滑移制作、豎直滑移支座） 非結(jié)點荷載的等效前邊在我們已經(jīng)提到的數(shù)據(jù)輸入，結(jié)構(gòu)的荷載因為作用位置的不同，分為結(jié)點荷載和非結(jié)點荷載。 邊界條件的處理具體方法為單獨開辟4個列寬的與結(jié)構(gòu)矩陣存儲有關(guān)結(jié)點的約束數(shù)據(jù)。 總剛度矩陣 總剛度矩陣反映了整個結(jié)構(gòu)的基礎(chǔ)上引入邊界條件的剛度而形成的整體剛度矩陣法。單元的局部坐標系中可將沿i、j 的方向