有限元上机实验报告结构数值分析与程序设计上机实验院系: 土木工程与力学学院专业: 土木工程班级: 姓名:学号: 指导教师: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 NTYPE6 4 12 60.1000E+01 0.000 1.0000.0000E+00 1 NODE 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+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 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+00 2、修改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,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,16 17,13,14 17,14,18 18,14,15 18,15,19 19,15,16 19,16,20 0,3 0,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 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+008 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.000000E+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.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.821554E+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.123014E+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+05 3、以例1为对象，研究单元细分对计算结果的影响。

3.1、单元细分如下所示:图1,单元细分图示3.2、输入数据文件修改为:15,16,30,10,1.0E0,0.0,1.0,0.0,1 3,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、输出数据文件修改为:NN NE ND NFIX E ANU T GM NTYPE15 16 30 100.1000E+01 0.000 1.0000.0000E+00 1 NODE 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.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.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-16ELEMENT 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 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.457838E-0115 -0.202577E-01 -0.620125E-01 0.116686E+0016 0.000000E+00 -0.620125E-01 -0.761703E-013.4、与习题1计算结果的对比:表1,单元细分的节点位移对比表节点号 X-DISP Y-DISP 节点号 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-023 -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、用程序完成习题3和4的分析。