有限元上机实验报告

《有限元上机实验报告》由会员分享，可在线阅读，更多相关《有限元上机实验报告（20页珍藏版）》请在装配图网上搜索。

1、结构数值分析与程序设计 上机实验院系：土木工程与力学学院专业：土木工程班级：姓名：学号：指导教师：1、调试教材P26-30程序FEM1。1.1、输入数据文件为：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、输出数据文件为：NN NE ND NFIX E ANU T GM NTYPE64 1260.1000E+01 0.000 1.0000.0000E+00 1NODEX-LOADY-LOAD10.000000E+00-0.100

2、000E+0120.000000E+000.000000E+0030.000000E+000.000000E+0040.000000E+000.000000E+0050.000000E+000.000000E+0060.000000E+000.000000E+00NODEX-DISPY-DISP1 -0.879121E-15-0.325275E+0120.879121E-16-0.125275E+013 -0.879121E-01-0.373626E+0040.117216E-15-0.835165E-1550.175824E+00-0.293040E-1560.175824E+000.263

3、736E-15ELEMENT X-STR Y-STR XY-STR1-0.879121E-01-0.200000E+010.439560E+0020.175824E+00-0.125275E+010.256410E-153-0.879121E-01-0.373626E+000.307692E+0040.000000E+00-0.373626E+00-0.131868E+002、修改FEM1，计算P31例2-2。2.1、输入数据文件为:20.0,24.0,40.0,5.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,1

4、010,6,710,7,1111,7,811,8,1213,9,1013,10,1414,10,1114,11,1515,11,1215,12,1617,13,1417,14,1818,14,1518,15,1919,15,1619,16,200,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,402.2、输出数据文件为：NN NE ND NFIX E ANU T GM NTYPE20 24 4050.2000E+11 0.167 1.0000.0

5、000E+00 1NODEX-LOADY-LOAD10.000000E+00-0.112500E+0520.000000E+000.000000E+0030.000000E+000.000000E+0040.000000E+000.000000E+0050.000000E+00-0.225000E+0560.000000E+000.000000E+0070.000000E+000.000000E+0080.000000E+000.000000E+0090.000000E+00-0.225000E+05100.000000E+000.000000E+00110.000000E+000.00000

6、0E+00120.000000E+000.000000E+00130.000000E+00-0.225000E+05140.000000E+000.000000E+00150.000000E+000.000000E+00160.000000E+000.000000E+00170.000000E+00-0.112500E+05180.000000E+000.000000E+00190.000000E+000.000000E+00200.000000E+000.000000E+00ELEMENT X-STRY-STRXY-STR1-0.135626E+060.821075E+040.409751E

7、+052-0.471063E+05-0.229544E+04-0.412344E+053-0.450136E+050.102356E+050.457906E+0540.432811E+050.278738E+04-0.362994E+0550.455430E+050.163319E+050.500168E+0560.138924E+060.659433E+04-0.367475E+057-0.118919E+06-0.142881E+050.486599E+058-0.379873E+05-0.749746E+04-0.269849E+059-0.385520E+05-0.108789E+05

8、0.506495E+05100.395343E+05-0.183501E+04-0.224095E+05110.381311E+05-0.102382E+050.461386E+05120.117794E+060.418256E+04-0.285486E+0513-0.821554E+05-0.148737E+050.433242E+0514-0.273330E+05-0.125220E+05-0.785096E+0415-0.274155E+05-0.130158E+050.484341E+05160.262559E+05-0.710079E+04-0.177296E+04170.25568

9、1E+05-0.112192E+050.430881E+05180.850804E+05-0.199044E+03-0.127141E+0519-0.308284E+05-0.131056E+050.306629E+0520-0.123014E+05-0.236279E+050.137015E+0521-0.105684E+05-0.132504E+050.387534E+05220.646544E+04-0.469412E+050.238659E+05230.122400E+05-0.123633E+050.349753E+05240.349927E+05-0.730929E+050.155

10、524E+053、以例1为对象，研究单元细分对计算结果的影响。3.1、单元细分如下所示:15图1:单元细分图示3.2、输入数据文件修改为：15,16,30,10,1.0E0,0.0,1.0,0.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

11、,3,7,13,21,22,24,26,28,303.3、输出数据文件修改为:NN NE ND NFIX E ANU T GM NTYPE15 16 30 100.1000E+01 0.000 1.0000.0000E+00 1NODEX-LOADY-LOAD0.000000E+00 -0.100000E+010.000000E+000.000000E+000.000000E+000.000000E+000.000000E+000.000000E+000.000000E+000.000000E+000.000000E+000.000000E+000.000000E+000.000000E+00

12、0.000000E+000.000000E+000.000000E+000.000000E+00100.000000E+000.000000E+00110.000000E+000.000000E+00120.000000E+000.000000E+00130.000000E+000.000000E+00140.000000E+000.000000E+00150.000000E+000.000000E+00NODEX-DISPY-DISP1-0.971690E-15-0.440129E+012-0.273631E-16-0.240129E+013-0.236860E+00-0.142960E+0

13、140.193766E-17-0.121245E+0150.841994E-01-0.102402E+016-0.257230E-01-0.406152E+0070.117892E-15-0.712244E+0080.783953E-01-0.389010E+0090.159794E+00-0.223862E+00100.149665E+00-0.310062E-01110.990200E-16-0.474829E-15120.148530E+00-0.271029E-15130.205578E+00-0.145183E-15140.225835E+00-0.257352E-16150.225

14、835E+000.761703E-16ELEMENT X-STRY-STRXY-STR1-0.473719E+00-0.400000E+010.971690E+0020.168399E+00-0.237768E+010.188429E+003-0.473719E+00-0.811160E+000.650631E+004-0.219845E+00-0.811160E+000.296808E+0050.156791E+00-0.100041E+010.323234E+0060.168399E+00-0.100041E+010.188429E+0070.162797E+00-0.127002E+01

15、0.170951E+008-0.219845E+00-0.364580E+000.432351E+009-0.202577E-01-0.364580E+000.733912E-02100.297060E+00-0.142449E+010.222672E-15110.156791E+00-0.778019E+000.253099E+00120.114095E+00-0.778019E+00-0.701347E-01130.162797E+00-0.447725E+000.119364E+00140.405153E-01-0.447725E+00-0.457838E-0115-0.202577E-

16、01-0.620125E-010.116686E+00160.000000E+00-0.620125E-01-0.761703E-013.4、与习题1计算结果的对比:表1 :单元细分的节点位移对比表节点号 X-DISP Y-DISP 节点号X-DISP Y-DISPX-误差Y-误差1-8.79E-16 -3.25E+001-9.72E-16 -4.40E+00-9.53E-02-2.61E-0128.79E-17-1.25E+0041.94E-18-1.21E+004.44E+013.32E-023-8.79E-02-3.74E-016-2.57E-02 7.62E-172.42E+00-4.

17、91E+1541.76E-012.64E-16119.90E-17-4.75E-161.78E+15-1.56E+0051.76E-01-2.93E-16132.26E-01-1.45E-02-2.21E-01-1.00E+004、用程序完成习题3和4的分析。课本第二章习题3,单元划分如下图所示:图2 :习题2-2图输入数据文件为：8,6,16,0,2E11,0.3,0.1,0,16,3,25,6,25,2,14,5,18,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.0 0输出数据文件为：NN NE ND NFI

18、X E ANU T GM NTYPENODEX-LOADY-LOAD10.000000E+000.000000E+0020.000000E+000.000000E+0030.707000E+03-0.707000E+0340.000000E+000.000000E+0050.000000E+000.000000E+0060.000000E+000.000000E+007-0.707000E+030.707000E+0380.000000E+000.000000E+00NODEX-DISPY-DISP1-0.177257E+000.493387E-012-0.177257E+00-0.25812

19、7E-013-0.177257E+00-0.100964E+004-0.102106E+000.493387E-015-0.102106E+00-0.258127E-016-0.102106E+00-0.100964E+0070.481974E-010.493387E-0180.481974E-01-0.258128E-01ELEMENT X-STRY-STRXY-STR1-0.590520E+04-0.800872E+040.668666E+042-0.776742E+04-0.729983E+040.755878E+0430.146578E+04-0.497822E+04-0.259785

20、E+044-0.709961E+030.530443E+04-0.460698E+045-0.405635E+04-0.960508E+040.794023E+046-0.218427E+04-0.531371E+040.443640E+045、采用SAP2000分析上述算例。(1) 上面习题1,输出结果如下表所示:表2 :上机题1spa2000计算的单元应力结果TABLE: Element Stresses - Area PlanesAreaAreaElemJointOutputCaseCaseTypeS11S22S33S12SMaxTextTextTextTextTextN/mm2N/mm

21、2N/mm2N/mm2N/mm2116DEADLinStatic-0.00008791-0.00200.00043960.000008296111DEADLinStatic-0.00008791-0.00200.00043960.000008296114DEADLinStatic-0.00008791-0.00200.00043960.000008296321DEADLinStatic0.0001758-0.001253000.0001758323DEADLinStatic0.0001758-0.001253000.0001758322DEADLinStatic0.0001758-0.0012

22、53000.0001758431DEADLinStatic-0.00008791-0.00037360-0.00030770.0001085434DEADLinStatic-0.00008791-0.00037360-0.00030770.0001085432DEADLinStatic-0.00008791-0.00037360-0.00030770.0001085544DEADLinStatic6.746E-20-0.00037360-0.00013190.00004185542DEADLinStatic0-0.00037360-0.00013190.00004185545DEADLinSt

23、atic-6.746E-20-0.00037360-0.00013190.00004185表3 :上机题1spa2000计算的节点位移结果TABLE: Joint ReactionsJointOutputCaseCaseTypeF1F2F3M1M2M3TextTextTextNNNN-mmN-mmN-mm1DEADLinStatic-131.87000002DEADLinStatic00439.560003DEADLinStatic-87.910626.370005DEADLinStatic00-65.930006DEADLinStatic219.7800000对比前面fortran编程计算得

24、出的结果知，两者的差距不大，有限单元 法的计算精度较高。(2) 输出结果如下表所示：表4 :习题2-2spa2000计算考察点输出结果TABLE: Element Stresses - Area PlanesAreaS11S22S33S12SMaxSMinSAngleSVMTextN/mm2N/mm2N/mm2N/mm2N/mm2N/mm2DegreesN/mm210.2440.0086920-0.0190.2460-4.6590.24220.1580.01300.0310.164011.6740.161270.1570.0093360-0.020.160-7.6240.157280.0770

25、.00546300.0290.087-0.0045819.3840.089290.0770.0076210-0.0210.0830-15.8860.08230-0.002050.00073600.0270.026-0.02746.4880.04631-0.001320.0050680-0.0230.025-0.021-48.9520.0432-0.08-0.0017600.0250.00553-0.08773.7130.0933-0.0790.0023070-0.0250.009415-0.086-74.1880.09134-0.158-0.0016700.0220.001468-0.1618

26、2.0120.16235-0.158-0.000450-0.0280.004322-0.162-80.2560.16536-0.2350.00196900.0180.00332-0.23785.6960.238(3) 三角形单元细分。图3 :单元细分的spa模型图示spa2000输出单元细分计算结果如下：表5 : spa2000计算单元细分各节点位移输出值TABLE: Joint DisplacementsJointTextU1mmU2mmU3mmR1RadiansR2RadiansR3Radians100-4.43779000200-2.437790003-0.2530470-1.38592

27、8000400-1.338626000500-0.601012000600000070.1438120000080.2280790000090.24515300000100.24515300000110.1658620-0.024095000120.1743990-0.216827000130.1321370-0.458571000140.0490040-0.9355100015-0.0353910-0.371081000表6 : spa2000计算单元细分各单元应力输出值TABLE: Element Stresses - Area PlanesAreaAreaElemJointS11S22S

28、12SMaxSMinSAngleSVMTextTextTextN/mm2N/mm2N/mm2N/mm2N/mm2DegreesN/mm2111-0.00051-0.0040.0010520-000.0041892249.8E-05-0.0022-0.0004030.0001667-0-9.6730.002355332-0.00051-0.0009-0.000750.00007189-0-37.6260.0015164415-0.00017-0.00090.00026240-000.00094585540.000264-0.00148-0.0001420.0002759-0-4.650.0016

29、4266139.8E-05-0.00095-0.000320.0001877-0-15.6580.00114977148.45E-05-0.00095-0.0001590.0001082-0-8.4940.0010368812-0.00017-0.00031-0.0003550.0001228-0-39.4280.000679970.000288-0.001200.0002876-000.001369101050.000264-0.00092-0.0001310.0002786-0-6.2410.001098111170.000169-0.000921.167E-050.0001687-00.

30、6160.0010121212138.45E-05-0.00043-0.0001880.0001456-0-17.9870.00058131313123.42E-05-0.000435.368E-050.00004023-06.4630.0004612141412-1.7E-05-4.8E-05-0.0001130.00008187-0-41.0950.00020115159-4.1E-20-4.8E-057.929E-050.00005878-036.5480.0001455161615-1.7E-05-0.00031-1.71E-050-000.0003018(4) 、课本第二章习题3的s

31、pa2000分析单元部分应力结果如下表所示：表7 : spa2000计算单元细分各单元应力输出值TABLE: Element Stresses - Area PlanesAreaElem Joint S11Text Text N/mm2S22S33N/mm2 N/mm226.84727.2936.84704.7410S12N/mm27.2936.847SMaxSMinSAngleSVMN/mm2N/mm214.14-0.447Degrees-45N/mm214.36912.982-0.94839.7213.481.887E-152.5533.6E-152.553902.5532.553-4.4

32、E-15-2E-152.553-4E-152.553节点位移计算结果如下表所示：表8 : spa2000计算单元细分各单元应力输出值TABLE: Joint DisplacementsJointU1U2U3R1R2R31mmmmmmRadiansRadiansRadians2-0.00011100.0000370003-0.00011500.0001260004-0.00009100.0003150005-0.00003800.0000410006-0.00002600.0001130007-0.00006600.0000170008-0.00004200.00008400013.769E-0

33、600.000291000注：部分输出列省略6、TRUSS：计算一个不少于5杆的桁架6.1、计算桁架如下图所示：10KN泗 / 泅图4 :桁架计算简图6.2、输入数据文件为：5,7,10,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、输出数据文件为:NN NE ND NFIX E57103 0.2000E+08ELEMENTNODE1NODE2 AREA1530.1000E+012540.1000E+013340.1000E+014

34、420.1000E+015320.1000E+016210.1000E+017310.1000E+01NODEX-COORDY-COORD10.1000E+03 0.5000E+0220.5000E+02 0.5000E+0230.5000E+02 0.0000E+0040.0000E+00 0.5000E+0250.0000E+00 0.0000E+00IFIX= 79 10NODEX-LOADY-LOAD10.0000E+00-0.1000E+0520.0000E+000.0000E+0030.0000E+000.0000E+0040.0000E+000.0000E+0050.0000E+

35、000.0000E+00NODEX-DISPY-DISP10.5000E-01-0.3164E+0020.2500E-01-0.1457E+003-0.5000E-01-0.1457E+0040.3694E-16-0.2500E-015-0.5000E-16-0.2500E-16ELEMENTFORCE1 -0.2000E+052 -0.1000E+053 0.1414E+054 0.1000E+055 0.0000E+006 0.1000E+057 -0.1414E+057、FRAME:计算一个不少于5杆的刚架7.1、刚架的计算图示如下：7.2、输入数据文件:6,6,18,6,2000000

36、00.1,2,0.010,0.00023,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、输出数据文件NN NE ND NFIX E66186 0.2000E+09ELEMENTNODE1NODE2 AREAFINT112 0.1000E-010.2000E-03234 0.1000E-010.2000E-03331 0.1000E-010.2000E-03442

37、 0.1000E-010.2000E-03553 0.1000E-010.2000E-03664 0.1000E-010.2000E-03NODEX-COORDY-COORD10.0000E+000.6600E+0120.6000E+010.6600E+0130.0000E+000.3600E+0140.6000E+010.3600E+0150.0000E+000.0000E+0060.6000E+010.0000E+00IFIX= 13 14 15 16 17 18nodalforces from applied loadsNODEX-LOADY-LOADMOMENT10.000000E+0

38、0-0.900000E+02-0.900000E+0220.000000E+00-0.900000E+020.900000E+0230.000000E+000.000000E+000.000000E+0040.000000E+000.000000E+000.000000E+0050.000000E+000.000000E+000.000000E+0060.000000E+000.000000E+000.000000E+00NODEX-DISPY-DISPTHITA10.444916E-04-0.297000E-03-0.153025E-022-0.444910E-04-0.297000E-03

39、0.153025E-023-0.534624E-04-0.162000E-030.352660E-0340.534624E-04-0.162000E-03-0.352660E-0350.581326E-18-0.162000E-18-0.154054E-186-0.581326E-18-0.162000E-180.154054E-18reaction nodal forces from the equationsNODEX-LOADY-LOADMOMENT1-0.410133E-04-0.899997E+02-0.899985E+0220.398815E-04-0.899997E+020.89

40、9984E+0230.538512E-040.167285E-03-0.345905E-034-0.497890E-040.163391E-030.356474E-035-0.598073E+010.900000E+020.684686E+0160.598072E+010.900000E+02-0.684685E+01ELEMENT AXI-FSHEAR-QMOMENT-M10.296609E+020.900000E+020.695967E+02-0.296609E+020.900000E+02-0.695967E+022-0.356416E+020.278123E-060.470214E+0

41、10.356416E+02-0.278123E-06-0.470214E+0130.900000E+02-0.296609E+02-0.193859E+02-0.900000E+020.296609E+02-0.695967E+0240.900000E+020.296609E+020.193859E+02-0.900000E+02-0.296609E+020.695967E+0250.900000E+020.598073E+010.684686E+01-0.900000E+02-0.598073E+010.146838E+0260.900000E+02-0.598072E+01-0.684685E+01-0.900000E+020.598072E+01-0.146838E+02