有限元程序課程設(shè)計(jì)

1、王順金：有限元程序課程設(shè)計(jì) 重慶大學(xué)本科學(xué)生課程設(shè)計(jì)任務(wù)書課程設(shè)計(jì)題目有限元程序設(shè)計(jì)學(xué)院資源及環(huán)境科學(xué)學(xué)院專業(yè)工程力學(xué)年級(jí)2010級(jí)已知參數(shù)和設(shè)計(jì)要求：1.獨(dú)立完成有限元程序設(shè)計(jì)。2.獨(dú)立選擇計(jì)算算例，并能通過(guò)算例判斷程序的正確性。3.獨(dú)立完成程序設(shè)計(jì)報(bào)告，報(bào)告內(nèi)容包括理論公式、程序框圖、程序本體、計(jì)算算例，算例結(jié)果分析、結(jié)論等。學(xué)生應(yīng)完成的工作：1. 復(fù)習(xí)掌握有限單元法的基本原理。2. 掌握彈性力學(xué)平面問(wèn)題3節(jié)點(diǎn)三角形單元或4節(jié)點(diǎn)等參單元有限元方法的計(jì)算流程，以及單元?jiǎng)偠染仃?、等效?jié)點(diǎn)載荷、節(jié)點(diǎn)應(yīng)變、節(jié)點(diǎn)應(yīng)力和高斯積分等的計(jì)算公式。3. 用Fortran語(yǔ)言編寫彈性力學(xué)平面問(wèn)題3節(jié)點(diǎn)三角形單

2、元或4節(jié)點(diǎn)等參單元的有限元程序。4. 在Visual Fortran 程序集成開發(fā)環(huán)境中完成有限元程序的編輯和調(diào)試工作。5. 利用編寫的有限元程序，計(jì)算算例，分析計(jì)算結(jié)果。6. 撰寫課程設(shè)計(jì)報(bào)告。目前資料收集情況（含指定參考資料）：1.王勖成，有限單元法，北京：高等教育出版社，2002。2.O.C. Zienkiewicz, R. L. Taylor, Finite Element Method, 5th Eition, McGraw-Hall Book Company Limited, 2000。3.張汝清，董明，結(jié)構(gòu)計(jì)算程序設(shè)計(jì)，重慶：重慶大學(xué)出版社，1988。課程設(shè)計(jì)的工作計(jì)劃：1. 第

3、1周星期一上午：教師講解程序設(shè)計(jì)方法，程序設(shè)計(jì)要求和任務(wù)安排。2. 第1周星期一至星期二完成程序框圖設(shè)計(jì)。3. 第1周星期三至第2周星期四完成程序設(shè)計(jì)。4. 第2周星期五完成課程設(shè)計(jì)報(bào)告。任務(wù)下達(dá)日期 2013 年 6 月 6 日完成日期 2013 年 07 月 03 日指導(dǎo)教師 （簽名）學(xué) 生 （簽名）一、前言有限單元法是在當(dāng)今工程分析中獲得最廣泛應(yīng)用的數(shù)值計(jì)算方法，由于其通用性和有效性，受到工程技術(shù)界的高度重視。伴隨著計(jì)算機(jī)科學(xué)和技術(shù)的快速發(fā)展，現(xiàn)已成為計(jì)算機(jī)輔助設(shè)計(jì)和計(jì)算機(jī)輔助制造的重要組成部分。由于有限元法是通過(guò)計(jì)算機(jī)實(shí)現(xiàn)的，因此有限元程序的編制及相關(guān)軟件的研發(fā)就變得尤為重要。從二十世

4、紀(jì)五十年代以來(lái)，有限元軟件的發(fā)展按目的和用途可分為專用軟件和大型通用商業(yè)軟件，而且軟件往往集成了網(wǎng)格自動(dòng)劃分，結(jié)果分析和顯示等前后處理功能，而且隨著時(shí)間的發(fā)展，大型通用商業(yè)軟件的功能由線性擴(kuò)展到非線性，由結(jié)構(gòu)擴(kuò)展到非結(jié)構(gòu)等等，這一系列強(qiáng)大功能的實(shí)現(xiàn)與運(yùn)用都要求我們對(duì)有限元法的基礎(chǔ)理論知識(shí)有較為清楚的認(rèn)識(shí)以及對(duì)程序編寫的基本能力有較好掌握。1.元分析程序的理論基礎(chǔ)作為彈性力學(xué)微分方程的等效積分形式，虛位移原理和虛應(yīng)原理分別是平衡方程與力的邊界條件和幾何方程與位移邊界條件的等效積分形式。將物理方程引入虛位移原理和虛應(yīng)力原理可以分別導(dǎo)出最小位能原理和最小余能原理，它們本質(zhì)上和等效積分的伽遼金“弱”形

5、式相一致，這是建立彈性力學(xué)有限原方程一般表達(dá)格式的理論基礎(chǔ)。對(duì)于通過(guò)彈性力學(xué)變分原理建立的彈性力學(xué)問(wèn)題有限元方法，其未知場(chǎng)變量是位移，以節(jié)點(diǎn)位移為基本未知量，并以最小位能原理為基礎(chǔ)建立的有限單元為位移元。彈性力學(xué)平面問(wèn)題有限元分析表達(dá)格式的建立步驟一般為：a的模型進(jìn)行單元離散，一般采用三角形單元或四邊形單元，不過(guò)由于四邊形單元具有更高精度，應(yīng)用更為普遍，一般而言一次或二次單元已足以滿足精度要求。離散后對(duì)單元進(jìn)行編號(hào)，并且給定單元各節(jié)點(diǎn)的整體編碼以及局部坐標(biāo)系下按逆時(shí)針局部編碼。b坐標(biāo)系下對(duì)各單元構(gòu)造形函得到由單元各節(jié)點(diǎn)位移表示的單元位移形式，進(jìn)而得到單元?jiǎng)偠染仃?，利用等參元性質(zhì)和雅閣比矩陣進(jìn)行

6、組集，建立整體剛度矩陣。c等效結(jié)點(diǎn)載荷列陣并組集成結(jié)構(gòu)節(jié)點(diǎn)等效載荷列陣得到單元各節(jié)點(diǎn)位移與結(jié)構(gòu)節(jié)點(diǎn)等效載荷列陣的線代方程組，求解得到位移，進(jìn)而可得應(yīng)力應(yīng)變。 根據(jù)以上一般步驟，編制相應(yīng)的計(jì)算機(jī)程序并采用數(shù)值積分方法處理有關(guān)數(shù)學(xué)計(jì)算就可以得到完整的彈性力學(xué)問(wèn)題有限元分析程序。2. 問(wèn)題有限元分析教學(xué)程序本程序可對(duì)二維彈性力學(xué)問(wèn)題行有限元分析計(jì)算。計(jì)算采用的單元形式為四邊形四節(jié)點(diǎn)單元或者四邊形八節(jié)點(diǎn)單元，對(duì)于對(duì)稱和非對(duì)稱矩陣采用變帶寬存儲(chǔ)方法，最終輸出結(jié)果為單元各節(jié)點(diǎn)的位移。二、平面4、8節(jié)點(diǎn)有限元公式及計(jì)算原理2.1基本公式（1） 有限元平衡方程（2） 單元?jiǎng)偠染仃嚕?）單元等效結(jié)點(diǎn)載荷2.2單

7、元位移插值及坐標(biāo)變換（1） 通過(guò)Serendipity四邊形單元格式構(gòu)造插值函數(shù)。對(duì)于4節(jié)點(diǎn)單元，插值函數(shù)為： 對(duì)于8節(jié)點(diǎn)單元，插值函數(shù)為： 其中： （i=1,2,3,4）（2） 位移插值 （3） 坐標(biāo)變換 2.3單元?jiǎng)偠鹊挠?jì)算（1）彈性矩陣平面應(yīng)力問(wèn)題：；平面應(yīng)變問(wèn)題：；式中E和v分別為材料Young氏模量和Poisson比。（2） 應(yīng)變矩陣形狀函數(shù)對(duì)局部坐標(biāo)偏導(dǎo)：，形狀函數(shù)對(duì)整體坐標(biāo)偏導(dǎo)：Jacobi矩陣：Jacobi矩陣的逆：其中：（3） 單元?jiǎng)偠染仃囉?jì)算單元?jiǎng)偠染仃嘒auss積分：?jiǎn)卧獎(jiǎng)偠染仃囎泳仃嚕?（i=1，2）2.4等效結(jié)點(diǎn)載荷（1） 集中載荷直接施加在結(jié)點(diǎn)上（2） 體積力產(chǎn)生的

8、等效結(jié)點(diǎn)載荷向量體積力產(chǎn)生的等效結(jié)點(diǎn)載荷Gauss積分：體積力產(chǎn)生的等效結(jié)點(diǎn)載荷子向量：其中：若體積力作用方向于整體坐標(biāo)系的y軸正方向角度，則：2.5單元應(yīng)變和應(yīng)力計(jì)算（1） 單元應(yīng)變單元中任意點(diǎn)的應(yīng)變：?jiǎn)卧腥我恻c(diǎn)的應(yīng)變：其中可作為Gauss積分點(diǎn)或單元結(jié)點(diǎn)的整體坐標(biāo)，其于局部坐標(biāo)間的關(guān)系為：,(2) 單元應(yīng)力三、程序結(jié)構(gòu)（1）有限元程序系統(tǒng)的組成及分析過(guò)程我們本次設(shè)計(jì)主要針對(duì)二維彈性力學(xué)靜力學(xué)問(wèn)題。有限元程序系統(tǒng)通常包括處理、有限元程序本體和后處理三部分。有限元軟件系統(tǒng)的組成及分析過(guò)程如下圖所示：CAD模型、材料參數(shù)、邊界條件、分析類型定義、網(wǎng)格劃分等。有限元分析本體程序前 處 理后 處

9、理核心：各種計(jì)算方法的實(shí)現(xiàn)。以圖形、曲線和表格等方式表達(dá)、分析計(jì)算結(jié)果。（2）有限元計(jì)算程序框圖 輸入離散模型數(shù)據(jù)型數(shù)據(jù)計(jì)算單元?jiǎng)偠汝嚱M集結(jié)構(gòu)剛度矩陣計(jì)算單元等效結(jié)點(diǎn)載荷組集結(jié)構(gòu)結(jié)點(diǎn)載荷列陣引入位移邊界條件求解線性方程組其它輔助計(jì)算輸出結(jié)果結(jié) 束單元循環(huán)形成 K形成 P消除 K 的奇異求解 Ka=P得結(jié)點(diǎn)位移 a計(jì)算應(yīng)力、應(yīng)變等（3）有限元程序主要編程技巧a.結(jié)構(gòu)模塊化主程序：組織調(diào)用子程序，過(guò)程清晰，不涉及復(fù)雜計(jì)算。子模塊：實(shí)現(xiàn)特定功能，可用多個(gè)子程序組成。子程序：實(shí)現(xiàn)簡(jiǎn)單功能，便于調(diào)用和擴(kuò)展。b.覆蓋技術(shù)在主程序中定義一個(gè)整數(shù)組和一個(gè)實(shí)型數(shù)組，分別存放整型和實(shí)型邊界數(shù)組。設(shè)計(jì)動(dòng)態(tài)數(shù)組表。動(dòng)

10、態(tài)數(shù)組覆蓋技術(shù)。（4）有限元程序主程序打開輸入數(shù)據(jù)，調(diào)用PMESH，調(diào)用PROFIL計(jì)算K的對(duì)角元地址和方程號(hào)，分配數(shù)組大小，調(diào)用ASSEM組集整體K，調(diào)用PLOAD組集載荷列向量，調(diào)用DATRI對(duì)剛度矩陣進(jìn)行三角分解，調(diào)用DASOL解方程，調(diào)用PRTDIS輸出結(jié)點(diǎn)位移，調(diào)用TRESSS計(jì)算GAUSS點(diǎn)和結(jié)點(diǎn)應(yīng)變和應(yīng)力。（5）變量說(shuō)明NUMNP：節(jié)點(diǎn)總數(shù)NUMEL：?jiǎn)卧倲?shù)NUMMAT：材料種類NLOAD：載荷數(shù)NDM：每個(gè)節(jié)點(diǎn)坐標(biāo)方向數(shù)NDF：節(jié)點(diǎn)自由度NEN：每個(gè)單元節(jié)點(diǎn)數(shù)。ID（K，J）：節(jié)點(diǎn)J的第K個(gè)自由度約束編號(hào)X（K，J）：節(jié)點(diǎn)J的第K個(gè)坐標(biāo)IX（K，J）：J單元中K節(jié)點(diǎn)的全局節(jié)點(diǎn)

11、編號(hào)F（K，J）：結(jié)點(diǎn)J第K方向上的集中載荷F（NDF，NLOAD）：集中載荷X（NDM，NUMNP）：結(jié)點(diǎn)坐標(biāo)IX（NEN，NUMEL）：?jiǎn)卧xF(K,J)：節(jié)點(diǎn)J在K方向上集中載荷大小ALFL=TRUE：非對(duì)稱矩陣集合ALFL=FALSE：對(duì)稱矩陣集合EE：楊氏彈性模量XNU：泊松比ITYPE：?jiǎn)栴}類型：?jiǎn)卧獎(jiǎng)偠染仃嚕狠d荷和內(nèi)力向量AD：對(duì)角元素AU：上三角元素AL：下三角元素JP：上（下）三角元素每行（列）最后元素LD：一個(gè)單元中各自由度方程數(shù)JD：列高ID：邊界條件，=0：無(wú)約束；=1：固定約束。（6）數(shù)據(jù)輸入格式整數(shù)：I5格式小數(shù)：F10.4格式（7）程序本體（見附錄）4、 程序驗(yàn)

12、證算例算例：平面板受軸向拉伸長(zhǎng)3m寬0.4m的平面板左端鉸接固定，右端受到水平方向向右、大小為10N的均布載荷，模型如下圖所示：F=10NF=10NF=10N1284519其中結(jié)構(gòu)參數(shù)：F=10N,L=3m,h=0.4m 材料參數(shù)：E=200GPa,=0.31、建立坐標(biāo)系，將板劃分成6個(gè)如上圖所示的四邊形四節(jié)點(diǎn)單元，進(jìn)行單元編號(hào)，均布載荷處理為右端三節(jié)點(diǎn)每個(gè)節(jié)點(diǎn)受水平向右大小10N的力，豎直方向不受力。1.1程序input文件中輸入：example(NEN=4) 12 6 1 3 2 2 4 1 1 1 0.0000 0.0000 2 0 0 1.0000 0.0000 3 0 0 2.000

13、0 0.0000 4 0 0 3.0000 0.0000 5 1 1 0.0000 0.2000 6 0 0 1.0000 0.2000 7 0 0 2.0000 0.2000 8 0 0 3.0000 0.2000 9 1 1 0.0000 0.4000 10 0 0 1.0000 0.4000 11 0 0 2.0000 0.4000 12 0 0 3.0000 0.4000 1 1 2 6 5 2 2 3 7 6 3 3 4 8 7 4 5 6 10 9 5 6 7 11 10 6 7 8 12 11 200.0E+9 0.3000 1 4 10.0000 0.0000 8 10.000

14、0 0.0000 12 10.0000 0.00001.2調(diào)試程序得到結(jié)果為：nodal displacements 1 0.0000E+00 0.0000E+00 2 0.3638E-09 0.2813E-10 3 0.7369E-09 0.2030E-10 4 0.1126E-08 0.2437E-10 5 0.0000E+00 0.0000E+00 6 0.3609E-09 -0.2981E-23 7 0.7429E-09 -0.1095E-22 8 0.1104E-08 -0.2132E-22 9 0.0000E+00 0.0000E+00 10 0.3638E-09 -0.2813E

15、-10 11 0.7369E-09 -0.2030E-10 12 0.1126E-08 -0.2437E-102、建立坐標(biāo)系，將板劃分成6個(gè)如上圖所示的四邊形八節(jié)點(diǎn)單元，進(jìn)行單元編號(hào)，均布載荷處理為右端三節(jié)點(diǎn)每個(gè)節(jié)點(diǎn)受水平向右大小10N的力，豎直方向不受力。其中結(jié)構(gòu)參數(shù)：F=10N,L=3m,h=0.4m 材料參數(shù)：E=200GPa,=0.32.1程序input文件中輸入：example(NEN=8) 29 6 1 3 2 2 8 1 1 1 0.0000 0.0000 2 0 0 1.0000 0.0000 3 0 0 2.0000 0.0000 4 0 0 3.0000 0.0000 5

16、1 1 0.0000 0.2000 6 0 0 1.0000 0.2000 7 0 0 2.0000 0.2000 8 0 0 3.0000 0.2000 9 1 1 0.0000 0.4000 10 0 0 1.0000 0.4000 11 0 0 2.0000 0.4000 12 0 0 3.0000 0.4000 13 0 0 0.5000 0.0000 14 0 0 1.5000 0.0000 15 0 0 2.5000 0.0000 16 1 1 0.0000 0.1000 17 0 0 1.0000 0.1000 18 0 0 2.0000 0.1000 19 0 0 3.0000

17、 0.1000 20 0 0 0.5000 0.2000 21 0 0 1.5000 0.2000 22 0 0 2.5000 0.2000 23 1 1 0.0000 0.3000 24 0 0 1.0000 0.3000 25 0 0 2.0000 0.3000 26 0 0 3.0000 0.3000 27 0 0 0.5000 0.4000 28 0 0 1.5000 0.4000 29 0 0 2.5000 0.4000 1 1 2 6 5 13 17 20 16 2 2 3 7 6 14 18 21 17 3 3 4 8 7 15 19 22 18 4 5 6 10 9 20 24

18、 27 23 5 6 7 11 10 21 25 28 24 6 7 8 12 11 22 26 29 25 200.0E+9 0.3000 1 4 10.0000 0.0000 8 10.0000 0.0000 12 10.0000 0.00002.2調(diào)試程序得到結(jié)果為：nodal displacements 1 0.0000E+00 0.0000E+00 2 0.3734E-09 0.2492E-10 3 0.7549E-09 0.4768E-10 4 0.1185E-08 0.6785E-10 5 0.0000E+00 0.0000E+00 6 0.3770E-09 -0.1562E-2

19、2 7 0.7371E-09 -0.3614E-22 8 0.1086E-08 -0.7600E-22 9 0.0000E+00 0.0000E+00 10 0.3734E-09 -0.2492E-10 11 0.7549E-09 -0.4768E-10 12 0.1185E-08 -0.6785E-10 13 0.1849E-09 0.2673E-10 14 0.5567E-09 0.1616E-10 15 0.9277E-09 0.1568E-11 16 0.0000E+00 0.0000E+00 17 0.3737E-09 0.1232E-10 18 0.7517E-09 0.2430E

20、-10 19 0.1114E-08 0.3177E-10 20 0.1810E-09 -0.5215E-23 21 0.5559E-09 -0.2565E-22 22 0.9462E-09 -0.5238E-22 23 0.0000E+00 0.0000E+00 24 0.3737E-09 -0.1232E-10 25 0.7517E-09 -0.2430E-10 26 0.1114E-08 -0.3177E-10 27 0.1849E-09 -0.2673E-10 28 0.5567E-09 -0.1616E-10 29 0.9277E-09 -0.1568E-113.數(shù)據(jù)分析其理論計(jì)算過(guò)程

21、如下，對(duì)于平面應(yīng)力問(wèn)題有： 其中 ： 代入計(jì)算得 任一距原點(diǎn)水平距離為x的點(diǎn)其位移為：理論解與數(shù)值解對(duì)比圖如下：可見接近自由端處兩曲線很接近，當(dāng)X=3m時(shí)，u=15x10E-9m，對(duì)應(yīng)輸出文件中的節(jié)點(diǎn)3、7、11的輸出值和3、7、11點(diǎn)有限元解較為接近，誤差為很小。同時(shí)可以發(fā)現(xiàn)靠近固定端處理論解與數(shù)值解誤差變大。分別對(duì)x=2處的理論解和有限元解進(jìn)行比較發(fā)現(xiàn)誤差較小。考慮到實(shí)際情況固定端的鉸支約束條件對(duì)實(shí)際位移的影響以及離散單元數(shù)目較少等可能造成的誤差，可以認(rèn)為該程序是正確的。另外數(shù)值由可以看出同的對(duì)稱性，符合基本力學(xué)理論結(jié)果。程序可用于平面問(wèn)題的求解，可獲得較好的位移分布。程序適用于平面應(yīng)力、

22、平面應(yīng)變狀態(tài)結(jié)構(gòu)4、8節(jié)點(diǎn)單元有限元分析，可求出結(jié)構(gòu)的位移、應(yīng)力以及應(yīng)變。五、算例1、算例：選擇超靜定梁中點(diǎn)受集中載荷作用的結(jié)構(gòu)，如下圖所示：其中結(jié)構(gòu)參數(shù)：F=10kN,L=4m,h=0.4m 材料參數(shù)：E=200GPa,=0.32、單元網(wǎng)格劃分及規(guī)定對(duì)上述結(jié)構(gòu)采用 101,102三種單元網(wǎng)格劃分形式，分別采用四結(jié)點(diǎn)和八結(jié)點(diǎn)四邊形單元。 規(guī)定單元編號(hào)先沿 y 正方向，再沿 x 正方向。為了方便比較，所有編號(hào)位置固定，例如第12結(jié)點(diǎn)，不管是101單元四結(jié)點(diǎn)，還是101單元八結(jié)點(diǎn)，其位置不變。其不同之處在于八結(jié)點(diǎn)增加的結(jié)點(diǎn)編號(hào)在四結(jié)點(diǎn)的基礎(chǔ)上插入，而不是替換原有的結(jié)點(diǎn)編號(hào)。特殊點(diǎn)如下所示：1161

23、1217221233116281712322(1).10X1單元四結(jié)點(diǎn)input文件中輸入：WSJ_10X1(NEN=4) 22 10 1 1 2 2 4 1 1 1 0.0000 0.0000 2 0 0 0.4000 0.0000 3 0 0 0.8000 0.0000 4 0 0 1.2000 0.0000 5 0 0 1.6000 0.0000 6 0 0 2.0000 0.0000 7 0 0 2.4000 0.0000 8 0 0 2.8000 0.0000 9 0 0 3.2000 0.0000 10 0 0 3.6000 0.0000 11 1 1 4.0000 0.0000

24、12 1 1 4.0000 0.4000 13 0 0 3.6000 0.4000 14 0 0 3.2000 0.4000 15 0 0 2.8000 0.4000 16 0 0 2.4000 0.4000 17 0 0 2.0000 0.4000 18 0 0 1.6000 0.4000 19 0 0 1.2000 0.4000 20 0 0 0.8000 0.4000 21 0 0 0.4000 0.4000 22 1 1 0.0000 0.4000 1 1 2 21 22 2 2 3 20 21 3 3 4 19 20 4 4 5 18 19 5 5 6 17 18 6 6 7 16

25、17 7 7 8 15 16 8 8 9 14 15 9 9 10 13 14 10 10 11 12 13 200.0E+9 0.3000 1 17 0.0000 -10.0E+310X1單元四結(jié)點(diǎn)input文件中輸出：nodal displacements 1 0.0000E+00 0.0000E+00 2 -0.2029E-06 -0.2673E-06 3 -0.3048E-06 -0.8380E-06 4 -0.3055E-06 -0.1509E-05 5 -0.2059E-06 -0.2083E-05 6 0.2994E-21 -0.2330E-05 7 0.2059E-06 -0.

26、2083E-05 8 0.3055E-06 -0.1509E-05 9 0.3048E-06 -0.8380E-06 10 0.2029E-06 -0.2673E-06 11 0.0000E+00 0.0000E+00 12 0.0000E+00 0.0000E+00 13 -0.2015E-06 -0.2671E-06 14 -0.3018E-06 -0.8376E-06 15 -0.3012E-06 -0.1510E-05 16 -0.1986E-06 -0.2077E-05 17 -0.1853E-21 -0.2365E-05 18 0.1986E-06 -0.2077E-05 19 0

27、.3012E-06 -0.1510E-05 20 0.3018E-06 -0.8376E-06 21 0.2015E-06 -0.2671E-06 22 0.0000E+00 0.0000E+00(2).10X1單元八結(jié)點(diǎn)input文件中輸入：WSJ_10X1(NEN=8) 53 10 1 1 2 2 8 1 1 1 0.0000 0.0000 2 0 0 0.4000 0.0000 3 0 0 0.8000 0.0000 4 0 0 1.2000 0.0000 5 0 0 1.6000 0.0000 6 0 0 2.0000 0.0000 7 0 0 2.4000 0.0000 8 0 0

28、2.8000 0.0000 9 0 0 3.2000 0.0000 10 0 0 3.6000 0.0000 11 1 1 4.0000 0.0000 12 1 1 4.0000 0.4000 13 0 0 3.6000 0.4000 14 0 0 3.2000 0.4000 15 0 0 2.8000 0.4000 16 0 0 2.4000 0.4000 17 0 0 2.0000 0.4000 18 0 0 1.6000 0.4000 19 0 0 1.2000 0.4000 20 0 0 0.8000 0.4000 21 0 0 0.4000 0.4000 22 1 1 0.0000

29、0.4000 23 0 0 0.2000 0.0000 24 0 0 0.6000 0.0000 25 0 0 1.0000 0.0000 26 0 0 1.4000 0.0000 27 0 0 1.8000 0.0000 28 0 0 2.2000 0.0000 29 0 0 2.6000 0.0000 30 0 0 3.0000 0.0000 31 0 0 3.4000 0.0000 32 0 0 3.8000 0.0000 33 1 1 4.0000 0.2000 34 0 0 3.6000 0.2000 35 0 0 3.2000 0.2000 36 0 0 2.8000 0.2000

30、 37 0 0 2.4000 0.2000 38 0 0 2.0000 0.2000 39 0 0 1.6000 0.2000 40 0 0 1.2000 0.2000 41 0 0 0.8000 0.2000 42 0 0 0.4000 0.2000 43 1 1 0.0000 0.2000 44 0 0 0.2000 0.4000 45 0 0 0.6000 0.4000 46 0 0 1.0000 0.4000 47 0 0 1.4000 0.4000 48 0 0 1.8000 0.4000 49 0 0 2.2000 0.4000 50 0 0 2.6000 0.4000 51 0

31、0 3.0000 0.4000 52 0 0 3.4000 0.4000 53 0 0 3.8000 0.4000 1 1 2 21 22 23 42 44 43 2 2 3 20 21 24 41 45 42 3 3 4 19 20 25 40 46 41 4 4 5 18 19 26 39 47 40 5 5 6 17 18 27 38 48 39 6 6 7 16 17 28 37 49 38 7 7 8 15 16 29 36 50 37 8 8 9 14 15 30 35 51 36 9 9 10 13 14 31 34 52 35 10 10 11 12 13 32 33 53 3

32、4 200.0E+9 0.3000 1 17 0.0000 -10.0E+310X1單元八結(jié)點(diǎn)input文件中輸出：nodal displacements 1 0.0000E+00 0.0000E+00 2 -0.2896E-06 -0.3882E-06 3 -0.4446E-06 -0.1198E-05 4 -0.4464E-06 -0.2169E-05 5 -0.2954E-06 -0.2984E-05 6 -0.1067E-19 -0.3336E-05 7 0.2954E-06 -0.2984E-05 8 0.4464E-06 -0.2169E-05 9 0.4446E-06 -0.11

33、98E-05 10 0.2896E-06 -0.3882E-06 11 0.0000E+00 0.0000E+00 12 0.0000E+00 0.0000E+00 13 -0.2881E-06 -0.3880E-06 14 -0.4419E-06 -0.1198E-05 15 -0.4415E-06 -0.2170E-05 16 -0.2928E-06 -0.2992E-05 17 0.1160E-19 -0.3406E-05 18 0.2928E-06 -0.2992E-05 19 0.4415E-06 -0.2170E-05 20 0.4419E-06 -0.1198E-05 21 0.

34、2881E-06 -0.3880E-06 22 0.0000E+00 0.0000E+00 23 -0.1582E-06 -0.1219E-06 24 -0.3874E-06 -0.7546E-06 25 -0.4638E-06 -0.1684E-05 26 -0.3906E-06 -0.2617E-05 27 -0.1666E-06 -0.3254E-05 28 0.1666E-06 -0.3254E-05 29 0.3906E-06 -0.2617E-05 30 0.4638E-06 -0.1684E-05 31 0.3874E-06 -0.7546E-06 32 0.1582E-06 -

35、0.1219E-06 33 0.0000E+00 0.0000E+00 34 0.7195E-09 -0.3635E-06 35 0.1556E-08 -0.1195E-05 36 0.2019E-08 -0.2174E-05 37 0.3347E-08 -0.3008E-05 38 0.7962E-21 -0.3388E-05 39 -0.3347E-08 -0.3008E-05 40 -0.2019E-08 -0.2174E-05 41 -0.1556E-08 -0.1195E-05 42 -0.7195E-09 -0.3635E-06 43 0.0000E+00 0.0000E+00 4

36、4 0.1574E-06 -0.1217E-06 45 0.3852E-06 -0.7543E-06 46 0.4599E-06 -0.1684E-05 47 0.3863E-06 -0.2615E-05 48 0.1608E-06 -0.3252E-05 49 -0.1608E-06 -0.3252E-05 50 -0.3863E-06 -0.2615E-05 51 -0.4599E-06 -0.1684E-05 52 -0.3852E-06 -0.7543E-06 53 -0.1574E-06 -0.1217E-06(3).10X2單元四結(jié)點(diǎn)input文件中輸入：WSJ_10X2(NEN=

37、4) 33 20 1 1 2 2 4 1 1 1 0.0000 0.0000 2 0 0 0.4000 0.0000 3 0 0 0.8000 0.0000 4 0 0 1.2000 0.0000 5 0 0 1.6000 0.0000 6 0 0 2.0000 0.0000 7 0 0 2.4000 0.0000 8 0 0 2.8000 0.0000 9 0 0 3.2000 0.0000 10 0 0 3.6000 0.0000 11 1 1 4.0000 0.0000 12 1 1 4.0000 0.4000 13 0 0 3.6000 0.4000 14 0 0 3.2000 0.4

38、000 15 0 0 2.8000 0.4000 16 0 0 2.4000 0.4000 17 0 0 2.0000 0.4000 18 0 0 1.6000 0.4000 19 0 0 1.2000 0.4000 20 0 0 0.8000 0.4000 21 0 0 0.4000 0.4000 22 1 1 0.0000 0.4000 23 1 1 0.0000 0.2000 24 0 0 0.4000 0.2000 25 0 0 0.8000 0.2000 26 0 0 1.2000 0.2000 27 0 0 1.6000 0.2000 28 0 0 2.0000 0.2000 29 0 0 2.4000 0.2000 30 0 0 2.8000 0.2000 31 0 0 3.2000 0.2000 32 0 0 3.6000 0.2000 33 1 1 4.0000 0.2000 1 1 2 24 23 2 2 3 25 24 3 3 4 2