有限元上機(jī)實(shí)驗(yàn)報(bào)告.doc

1、有限元上機(jī)實(shí)驗(yàn)報(bào)告結(jié)構(gòu)數(shù)值分析與程序設(shè)計(jì)上機(jī)實(shí)驗(yàn)院系 :土木工程與力學(xué)學(xué)院專業(yè) :土木工程班級(jí):姓名:學(xué)號(hào) :指導(dǎo)教師 :1、調(diào)試教材 P26-30 程序 FEM1。1.1 、輸入數(shù)據(jù)文件為 :6,4,12,6,1.0E0,0.0,1.0,0.0,13,1,25,2,43,2,56,3,50.0,2.00.0,1.01.0,1.00.0,0.01.0,0.02.0,0.01,3,7,8,10,121.2 、輸出數(shù)據(jù)文件為 :NN NE ND NFIX E ANU T GM NTYPE6 4 12 60.1000E+01 0.000 1.0000.0000E+00 1NODE X-LOAD Y-

2、LOAD1 0.000000E+00 -0.100000E+012 0.000000E+00 0.000000E+003 0.000000E+00 0.000000E+004 0.000000E+00 0.000000E+005 0.000000E+00 0.000000E+006 0.000000E+00 0.000000E+00NODE X-DISP Y-DISP1 -0.879121E-15 -0.325275E+012 0.879121E-16 -0.125275E+013 -0.879121E-01 -0.373626E+004 0.117216E-15 -0.835165E-155

3、 0.175824E+00 -0.293040E-156 0.175824E+00 0.263736E-15ELEMENT X-STR Y-STR XY-STR1 -0.879121E-01 -0.200000E+01 0.439560E+002 0.175824E+00 -0.125275E+01 0.256410E-153 -0.879121E-01 -0.373626E+00 0.307692E+004 0.000000E+00 -0.373626E+00 -0.131868E+002、修改 FEM1,計(jì)算 P31例 2-2 。2.1 、輸入數(shù)據(jù)文件為 :20.0,24.0,40.0,5

4、.0,2E10,0.167,1.0,0.0,1.05,1,26,5,26,2,37,6,37,3,47,4,89,5,69,6,1010,6,710,7,11 11,7,811,8,12 13,9,10 13,10,14 14,10,11 14,11,15 15,11,12 15,12,1617,13,14 17,14,18 18,14,15 18,15,19 19,15,16 19,16,20 0,30,20,10,02.25,32.25,22.25,12.25,04.5,34.5,24.5,14.5,06.75,36.75,26.75,16.75,09,39,29,19,01,3,5,7,

5、402.2 、輸出數(shù)據(jù)文件為 :NN NE ND NFIX E ANU T GM NTYPE20 24 40 50.2000E+11 0.167 1.0000.0000E+00 1 NODE X-LOAD Y-LOAD1 0.000000E+00 -0.112500E+052 0.000000E+00 0.000000E+003 0.000000E+00 0.000000E+004 0.000000E+00 0.000000E+005 0.000000E+00 -0.225000E+056 0.000000E+00 0.000000E+007 0.000000E+00 0.000000E+00

6、8 0.000000E+00 0.000000E+009 0.000000E+00 -0.225000E+0510 0.000000E+00 0.000000E+0011 0.000000E+00 0.000000E+0012 0.000000E+00 0.000000E+0013 0.000000E+00 -0.225000E+0514 0.000000E+00 0.000000E+0015 0.000000E+00 0.000000E+0016 0.000000E+00 0.000000E+0017 0.000000E+00 -0.112500E+0518 0.000000E+00 0.0

7、00000E+0019 0.000000E+00 0.000000E+0020 0.000000E+00 0.000000E+00ELEMENT X-STR Y-STR XY-STR1 -0.135626E+06 0.821075E+04 0.409751E+052 -0.471063E+05 -0.229544E+04 -0.412344E+053 -0.450136E+05 0.102356E+05 0.457906E+054 0.432811E+05 0.278738E+04 -0.362994E+055 0.455430E+05 0.163319E+05 0.500168E+056 0

8、.138924E+06 0.659433E+04 -0.367475E+057 -0.118919E+06 -0.142881E+05 0.486599E+058 -0.379873E+05 -0.749746E+04 -0.269849E+059 -0.385520E+05 -0.108789E+05 0.506495E+0510 0.395343E+05 -0.183501E+04 -0.224095E+0511 0.381311E+05 -0.102382E+05 0.461386E+0512 0.117794E+06 0.418256E+04 -0.285486E+0513 -0.82

9、1554E+05 -0.148737E+05 0.433242E+0514 -0.273330E+05 -0.125220E+05 -0.785096E+0415 -0.274155E+05 -0.130158E+05 0.484341E+0516 0.262559E+05 -0.710079E+04 -0.177296E+0417 0.255681E+05 -0.112192E+05 0.430881E+0518 0.850804E+05 -0.199044E+03 -0.127141E+0519 -0.308284E+05 -0.131056E+05 0.306629E+0520 -0.1

10、23014E+05 -0.236279E+05 0.137015E+0521 -0.105684E+05 -0.132504E+05 0.387534E+0522 0.646544E+04 -0.469412E+05 0.238659E+0523 0.122400E+05 -0.123633E+05 0.349753E+0524 0.349927E+05 -0.730929E+05 0.155524E+053、以例 1 為對(duì)象，研究單元細(xì)分對(duì)計(jì)算結(jié)果的影響。3.1 、單元細(xì)分如下所示 :圖 1, 單元細(xì)分圖示3.2 、輸入數(shù)據(jù)文件修改為 :15,16,30,10,1.0E0,0.0,1.0,0

11、.0,13,1,25,2,43,2,56,3,58,4,75,4,79,5,86,5,910,6,912,7,118,7,1213,8,129,8,1314,9,1310,9,1415,10,140.0,2.00.0,1.50.5,1.50.0,1.00.5,1.01.0,1.00.0,0.50.5,0.51.0,0.51.5,0.50.0,0.00.5,0.01.0,0.01.5,0.02.0,0.01,3,7,13,21,22,24,26,28,303.3 、輸出數(shù)據(jù)文件修改為 :NN NE ND NFIX E ANU T GM NTYPE15 16 30 100.1000E+01 0.0

12、00 1.0000.0000E+00 1NODE X-LOAD Y-LOAD1 0.000000E+00 -0.100000E+012 0.000000E+00 0.000000E+003 0.000000E+00 0.000000E+004 0.000000E+00 0.000000E+005 0.000000E+00 0.000000E+006 0.000000E+00 0.000000E+007 0.000000E+00 0.000000E+008 0.000000E+00 0.000000E+009 0.000000E+00 0.000000E+0010 0.000000E+00 0.

13、000000E+0011 0.000000E+00 0.000000E+0012 0.000000E+00 0.000000E+0013 0.000000E+00 0.000000E+0014 0.000000E+00 0.000000E+0015 0.000000E+00 0.000000E+00NODE X-DISP Y-DISP1 -0.971690E-15 -0.440129E+012 -0.273631E-16 -0.240129E+013 -0.236860E+00 -0.142960E+014 0.193766E-17 -0.121245E+015 0.841994E-01 -0

14、.102402E+016 -0.257230E-01 -0.406152E+007 0.117892E-15 -0.712244E+008 0.783953E-01 -0.389010E+009 0.159794E+00 -0.223862E+0010 0.149665E+00 -0.310062E-0111 0.990200E-16 -0.474829E-1512 0.148530E+00 -0.271029E-1513 0.205578E+00 -0.145183E-1514 0.225835E+00 -0.257352E-1615 0.225835E+00 0.761703E-16ELE

15、MENT X-STR Y-STR XY-STR1 -0.473719E+00 -0.400000E+01 0.971690E+002 0.168399E+00 -0.237768E+01 0.188429E+003 -0.473719E+00 -0.811160E+00 0.650631E+004 -0.219845E+00 -0.811160E+00 0.296808E+005 0.156791E+00 -0.100041E+01 0.323234E+006 0.168399E+00 -0.100041E+01 0.188429E+007 0.162797E+00 -0.127002E+01

16、 0.170951E+008 -0.219845E+00 -0.364580E+00 0.432351E+009 -0.202577E-01 -0.364580E+00 0.733912E-0210 0.297060E+00 -0.142449E+01 0.222672E-1511 0.156791E+00 -0.778019E+00 0.253099E+0012 0.114095E+00 -0.778019E+00 -0.701347E-0113 0.162797E+00 -0.447725E+00 0.119364E+0014 0.405153E-01 -0.447725E+00 -0.4

17、57838E-0115 -0.202577E-01 -0.620125E-01 0.116686E+0016 0.000000E+00 -0.620125E-01 -0.761703E-013.4 、與習(xí)題 1 計(jì)算結(jié)果的對(duì)比 :表 1, 單元細(xì)分的節(jié)點(diǎn)位移對(duì)比表節(jié)點(diǎn)號(hào) X-DISP Y-DISP節(jié)點(diǎn)號(hào) X-DISP Y-DISP X- 誤差 Y- 誤差1 -8.79E-16 -3.25E+00 1 -9.72E-16 -4.40E+00 -9.53E-02 -2.61E-012 8.79E-17 -1.25E+00 4 1.94E-18 -1.21E+00 4.44E+01 3.32E-02

18、3 -8.79E-02 -3.74E-01 6 -2.57E-02 7.62E-17 2.42E+00 -4.91E+154 1.76E-01 2.64E-16 11 9.90E-17 -4.75E-16 1.78E+15 -1.56E+005 1.76E-01 -2.93E-16 13 2.26E-01 -1.45E-02 -2.21E-01 -1.00E+006 1.17E-16 -8.35E-16 15 2.06E-01 -4.06E-01 -1.00E+00 -1.00E+004、用程序完成習(xí)題3 和 4 的分析。課本第二章習(xí)題3，單元?jiǎng)澐秩缦聢D所示:圖 2, 習(xí)題 2-2 圖 輸入數(shù)

19、據(jù)文件為 : 8,6,16,0,2E11,0.3,0.1,0,1 6,3,2 5,6,2 5,2,1 4,5,1 8,5,48,4,70.0,0.01.0,0.02.0,0.00.0,1.01.0,1.02.0,1.00.0,3.01.0,3.00輸出數(shù)據(jù)文件為 :NN NE ND NFIX E ANU T GM NTYPE8 6 16 00.2000E+12 0.300 0.1000.0000E+00 1 NODE X-LOAD Y-LOAD1 0.000000E+00 0.000000E+002 0.000000E+00 0.000000E+003 0.707000E+03 -0.7070

20、00E+034 0.000000E+00 0.000000E+005 0.000000E+00 0.000000E+006 0.000000E+00 0.000000E+007 -0.707000E+03 0.707000E+038 0.000000E+00 0.000000E+00NODE X-DISP Y-DISP1 -0.177257E+00 0.493387E-012 -0.177257E+00 -0.258127E-013 -0.177257E+00 -0.100964E+004 -0.102106E+00 0.493387E-015 -0.102106E+00 -0.258127E

21、-016 -0.102106E+00 -0.100964E+007 0.481974E-01 0.493387E-018 0.481974E-01 -0.258128E-01ELEMENT X-STR Y-STR XY-STR1 -0.590520E+04 -0.800872E+04 0.668666E+042 -0.776742E+04 -0.729983E+04 0.755878E+043 0.146578E+04 -0.497822E+04 -0.259785E+044 -0.709961E+03 0.530443E+04 -0.460698E+045 -0.405635E+04 -0.

22、960508E+04 0.794023E+046 -0.218427E+04 -0.531371E+04 0.443640E+045、采用 SAP2000分析上述算例。(1) 上面習(xí)題 1，輸出結(jié)果如下表所示 :表 2, 上機(jī)題 1spa2000 計(jì)算的單元應(yīng)力結(jié)果TABLE: Element Stresses - Area Planes Area AreaElem Joint OutputCase CaseType S11 S22 S33 S12 SMax SMin SAngle SVM Text Text Text Text Text N/mm2 N/mm2 N/mm2 N/mm2 N/m

23、m2 N/mm2 Degrees N/mm2 1 1 6 DEAD LinStatic - 0.00008791 -0.002 0 0.0004396 0.000008296 -0.002096 12.346 0.0021 1 1 1DEAD LinStatic -0.00008791 -0.002 0 0.0004396 0.000008296 -0.002096 12.346 0.0021 1 1 4 DEAD LinStatic -0.00008791 -0.002 0 0.00043960.000008296 -0.002096 12.346 0.0021 3 2 1 DEAD Lin

24、Static 0.0001758 -0.001253 0 0 0.0001758 -0.001253 0 0.001349 3 2 3 DEAD LinStatic0.0001758 -0.001253 0 0 0.0001758 -0.001253 0 0.001349 3 2 2 DEAD LinStatic 0.0001758 -0.001253 0 0 0.0001758 -0.001253 0 0.001349 4 3 1 DEAD LinStatic -0.00008791 -0.0003736 0 -0.0003077 0.0001085 -0.00057 - 32.548 0.

25、00063134 3 4 DEAD LinStatic -0.00008791 -0.0003736 0 -0.0003077 0.0001085 -0.00057 -32.548 0.00063134 3 2 DEAD LinStatic -0.00008791 -0.0003736 0 -0.0003077 0.0001085 -0.00057 -32.548 0.00063135 4 4 DEAD LinStatic 6.746E-20 -0.0003736 0 -0.0001319 0.00004185 -0.0004155 -17.609 0.00043795 4 2 DEAD Li

26、nStatic 0 -0.0003736 0 -0.0001319 0.00004185 -0.0004155-17.609 0.00043795 4 5 DEAD LinStatic -6.746E-20 -0.0003736 0 -0.0001319 0.00004185 -0.0004155 -17.609 0.0004379表 3, 上機(jī)題 1spa2000 計(jì)算的節(jié)點(diǎn)位移結(jié)果TABLE: Joint ReactionsJoint OutputCase CaseType F1 F2 F3 M1 M2 M3 Text Text Text N N N N-mm N-mm N-mm 1 DE

27、AD LinStatic -131.87 0 0 0 0 0 2 DEAD LinStatic 0 0439.56 0 0 0 3 DEAD LinStatic -87.91 0 626.37 0 0 0 5 DEAD LinStatic 0 0-65.93 0 0 0 6 DEAD LinStatic 219.78 0 0 0 0 0對(duì)比前面 fortran編程計(jì)算得出的結(jié)果知，兩者的差距不大，有限單元法的計(jì)算精度較高。(2) 輸出結(jié)果如下表所示 :表 4, 習(xí)題 2-2spa2000 計(jì)算考察點(diǎn)輸出結(jié)果TABLE: Element Stresses - Area PlanesArea S1

28、1 S22 S33 S12 SMax SMin SAngle SVM Text N/mm2 N/mm2 N/mm2N/mm2 N/mm2 N/mm2 Degrees N/mm2 1 0.244 0.008692 0 -0.019 0.246 0 -4.6590.242 2 0.158 0.013 0 0.031 0.164 0 11.674 0.161 27 0.157 0.009336 0 -0.02 0.16 0 -7.624 0.157 28 0.077 0.005463 0 0.029 0.087 -0.00458 19.3840.089 29 0.077 0.007621 0 -0.

29、021 0.083 0 -15.886 0.082 30 -0.002050.000736 0 0.027 0.026 -0.027 46.488 0.046 31 -0.00132 0.005068 0 -0.0230.025 -0.021 -48.952 0.04 32 -0.08 -0.00176 0 0.025 0.00553 -0.08773.713 0.09 33 -0.079 0.002307 0 -0.025 0.009415 -0.086 -74.188 0.091 34 -0.158 -0.00167 0 0.022 0.001468 -0.161 82.012 0.162

30、 35 -0.158 -0.000450 -0.028 0.004322 -0.162 -80.256 0.165 36 -0.235 0.001969 0 0.0180.00332 -0.237 85.696 0.238(3) 三角形單元細(xì)分。圖 3, 單元細(xì)分的 spa 模型圖示spa2000 輸出單元細(xì)分計(jì)算結(jié)果如下:表 5,spa2000 計(jì)算單元細(xì)分各節(jié)點(diǎn)位移輸出值TABLE: Joint DisplacementsJoint U1 U2 U3 R1 R2 R3 Text mm mm mm Radians Radians Radians 1 0 0 -4.43779 0 0 0 2

31、0 0 -2.43779 0 0 0 3 -0.253047 0 -1.385928 0 0 0 4 0 0 -1.338626000500-0.6010120006000000708 0.228079 0 0 0 0 0 9 0.245153 0 0 0 0 0 10 0.245153 0 0 0 0 0 110.165862 0 -0.024095 0 0 0 12 0.174399 0 -0.216827 0 0 0 13 0.132137 0 -0.458571 0 0 0 14 0.049004 0 -0.93551 0 0 0 15 -0.035391 0

32、 -0.371081 0 00表 6,spa2000 計(jì)算單元細(xì)分各單元應(yīng)力輸出值TABLE: Element Stresses - Area PlanesArea AreaElem Joint S11 S22 S12 SMax SMin SAngle SVM Text Text TextN/mm2 N/mm2 N/mm2 N/mm2 N/mm2 Degrees N/mm2 1 1 1 -0.00051 -0.0040.001052 0 -0 0 0.004189 2 2 4 9.8E-05 -0.0022 -0.000403 0.0001667 -0 -9.673 0.002355 3 3

33、2 -0.00051 -0.0009 -0.00075 0.00007189 -0 -37.6260.001516 4 4 15 -0.00017 -0.0009 0.0002624 0 -0 0 0.0009458 5 5 40.000264 -0.00148 -0.000142 0.0002759 -0 -4.65 0.001642 6 6 13 9.8E-05 -0.00095 -0.00032 0.0001877 -0 -15.658 0.001149 7 7 14 8.45E-05 -0.00095 -0.000159 0.0001082 -0 -8.494 0.001036 8 8

34、 12 -0.00017 -0.00031 -0.000355 0.0001228 -0 -39.428 0.00067 9 9 7 0.000288 -0.0012 0 0.0002876 -0 0 0.001369 10 10 5 0.000264 -0.00092 -0.000131 0.0002786 -0 -6.2410.001098 11 11 7 0.000169 -0.00092 1.167E-05 0.0001687 -0 0.616 0.00101212 12 13 8.45E-05 -0.00043 -0.000188 0.0001456 -0 -17.987 0.000

35、5813 1313 12 3.42E-05 -0.00043 5.368E-05 0.00004023 -0 6.463 0.0004612 14 14 12 -1.7E-05 -4.8E-05 -0.000113 0.00008187 -0 -41.095 0.000201 15 15 9 -4.1E-20 -4.8E-05 7.929E-05 0.00005878 -0 36.548 0.0001455 16 16 15 -1.7E-05 -0.00031 -1.71E-05 0 -0 0 0.0003018(4) 、課本第二章習(xí)題 3 的 spa2000 分析單元部分應(yīng)力結(jié)果如下表所示

36、:表 7,spa2000 計(jì)算單元細(xì)分各單元應(yīng)力輸出值TABLE: Element Stresses - Area Planes AreaElem Joint S11 S22 S33 S12 SMax SMin SAngle SVMText Text N/mm2 N/mm2 N/mm2 N/mm2 N/mm2 N/mm2 Degrees N/mm21 2 6.847 6.847 0 -7.293 14.14 -0.447 -45 14.3692 2 7.293 4.741 0 6.847 12.982 -0.948 39.72 13.483 1 1.887E-15 2.553 0 3.6E-1

37、5 2.553 0 90 2.5534 1 2.553 -4.4E-15 0 -2E-15 2.553 -4E-15 0 2.5535 4 6.847 6.847 0 -7.293 14.14 -0.447 -45 14.3696 4 4.741 7.293 0 -6.847 12.982 -0.948 -50.28 13.48節(jié)點(diǎn)位移計(jì)算結(jié)果如下表所示:表 8,spa2000 計(jì)算單元細(xì)分各單元應(yīng)力輸出值TABLE: Joint DisplacementsJoint U1 U2 U3 R1 R2 R3 1 mm mm mm Radians Radians Radians 2 - 0.0001

38、11 0 0.000037 0 0 0 3 -0.000115 0 0.000126 0 0 0 4 -0.000091 00.000315 0 0 0 5 -0.000038 0 0.000041 0 0 0 6 -0.000026 0 0.000113 0 0 07 -0.000066 0 0.000017 0 0 0 8 -0.000042 0 0.000084 0 0 0 1 3.769E-06 00.000291 0 0 0注 , 部分輸出列省略6 、TRUSS:計(jì)算一個(gè)不少于5 桿的桁架6.1 、計(jì)算桁架如下圖所示 :圖 4, 桁架計(jì)算簡(jiǎn)圖6.2 、輸入數(shù)據(jù)文件為 :5,7,10,

39、3,200000005,3,1.05,4,1.03,4,1.04,2,1.03,2,1.02,1,1.03,1,1.0100.0,50.050.0,50.050.0,0.00,50.00,07,9,106.3 、輸出數(shù)據(jù)文件為 :NN NE ND NFIX E 5 7 10 3 0.2000E+08ELEMENT NODE1 NODE2 AREA 1 5 3 0.1000E+012 5 4 0.1000E+013 3 4 0.1000E+014 4 2 0.1000E+015 3 2 0.1000E+016 2 1 0.1000E+017 3 1 0.1000E+01NODE X-COORD

40、Y-COORD1 0.1000E+03 0.5000E+022 0.5000E+02 0.5000E+023 0.5000E+02 0.0000E+004 0.0000E+00 0.5000E+025 0.0000E+00 0.0000E+00IFIX= 7 9 10NODE X-LOAD Y-LOAD1 0.0000E+00 -0.1000E+052 0.0000E+00 0.0000E+003 0.0000E+00 0.0000E+004 0.0000E+00 0.0000E+005 0.0000E+00 0.0000E+00NODE X-DISP Y-DISP1 0.5000E-01 -

41、0.3164E+002 0.2500E-01 -0.1457E+003 -0.5000E-01 -0.1457E+004 0.3694E-16 -0.2500E-015 -0.5000E-16 -0.2500E-16ELEMENT FORCE1 -0.2000E+052 -0.1000E+053 0.1414E+054 0.1000E+055 0.0000E+006 0.1000E+057 -0.1414E+057、 FRAME:計(jì)算一個(gè)不少于 5 桿的剛架7.1 、剛架的計(jì)算圖示如下 :圖 5, 剛架計(jì)算簡(jiǎn)圖 7.2 、輸入數(shù)據(jù)文件 : 6,6,18,6,200000000. 1,2,0.0

42、10,0.0002 3,4,0.010,0.00023,1,0.010,0.00024,2,0.010,0.00025,3,0.010,0.00026,4,0.010,0.00020,6.66,6.60,3.66,3.60,06,013,14,15,16,17,180,21,-302,-30.7.3 、輸出數(shù)據(jù)文件NN NE ND NFIX E 6 6 18 6 0.2000E+09ELEMENT NODE1 NODE2 AREA FINT 1 1 2 0.1000E-01 0.2000E-032 3 4 0.1000E-01 0.2000E-033 3 1 0.1000E-01 0.2000E-034 4 2 0.1000E-01 0.2000E-035 5 3 0.1000E-01 0.2000E-036 6 4 0.1000E-01 0.2000E-03NODE X-COORD Y-COORD 1 0.0000E+00 0.6600E+01 2 0.6000E+01 0.6600E+01 3 0.0000E+00 0.3600E+01 4 0.6000E+01 0.3600E+015 0.0000E+00 0.0000E+006 0.6000E+01 0.0000E+0