有限元上机实验报告.docx
- 文档编号:12698436
- 上传时间:2023-04-21
- 格式:DOCX
- 页数:27
- 大小:81.39KB
有限元上机实验报告.docx
《有限元上机实验报告.docx》由会员分享,可在线阅读,更多相关《有限元上机实验报告.docx(27页珍藏版)》请在冰豆网上搜索。
有限元上机实验报告
结构数值分析与程序设计上机实验
院系:
土木工程与力学学院
专业:
土木工程
班级:
姓名:
学号:
指导教师:
1、调试教材P26-30程序FEM1。
1.1、输入数据文件为:
6,4,12,6,1.0E0,0.0,1.0,0.0,1
3,1,2
5,2,4
3,2,5
6,3,5
0.0,2.0
0.0,1.0
1.0,1.0
0.0,0.0
1.0,0.0
2.0,0.0
1,3,7,8,10,12
1.2、输出数据文件为:
NNNENDNFIXEANUTGMNTYPE
641260.1000E+010.0001.0000.0000E+001
NODEX-LOADY-LOAD
10.000000E+00-0.100000E+01
20.000000E+000.000000E+00
30.000000E+000.000000E+00
40.000000E+000.000000E+00
50.000000E+000.000000E+00
60.000000E+000.000000E+00
NODEX-DISPY-DISP
1-0.879121E-15-0.325275E+01
20.879121E-16-0.125275E+01
3-0.879121E-01-0.373626E+00
40.117216E-15-0.835165E-15
50.175824E+00-0.293040E-15
60.175824E+000.263736E-15
ELEMENTX-STRY-STRXY-STR
1-0.879121E-01-0.200000E+010.439560E+00
20.175824E+00-0.125275E+010.256410E-15
3-0.879121E-01-0.373626E+000.307692E+00
40.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.0
5,1,2
6,5,2
6,2,3
7,6,3
7,3,4
7,4,8
9,5,6
9,6,10
10,6,7
10,7,11
11,7,8
11,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,2
0,1
0,0
2.25,3
2.25,2
2.25,1
2.25,0
4.5,3
4.5,2
4.5,1
4.5,0
6.75,3
6.75,2
6.75,1
6.75,0
9,3
9,2
9,1
9,0
1,3,5,7,40
2.2、输出数据文件为:
NNNENDNFIXEANUTGMNTYPE
20244050.2000E+110.1671.0000.0000E+001
NODEX-LOADY-LOAD
10.000000E+00-0.112500E+05
20.000000E+000.000000E+00
30.000000E+000.000000E+00
40.000000E+000.000000E+00
50.000000E+00-0.225000E+05
60.000000E+000.000000E+00
70.000000E+000.000000E+00
80.000000E+000.000000E+00
90.000000E+00-0.225000E+05
100.000000E+000.000000E+00
110.000000E+000.000000E+00
120.000000E+000.000000E+00
130.000000E+00-0.225000E+05
140.000000E+000.000000E+00
150.000000E+000.000000E+00
160.000000E+000.000000E+00
170.000000E+00-0.112500E+05
180.000000E+000.000000E+00
190.000000E+000.000000E+00
200.000000E+000.000000E+00
ELEMENTX-STRY-STRXY-STR
1-0.135626E+060.821075E+040.409751E+05
2-0.471063E+05-0.229544E+04-0.412344E+05
3-0.450136E+050.102356E+050.457906E+05
40.432811E+050.278738E+04-0.362994E+05
50.455430E+050.163319E+050.500168E+05
60.138924E+060.659433E+04-0.367475E+05
7-0.118919E+06-0.142881E+050.486599E+05
8-0.379873E+05-0.749746E+04-0.269849E+05
9-0.385520E+05-0.108789E+050.506495E+05
100.395343E+05-0.183501E+04-0.224095E+05
110.381311E+05-0.102382E+050.461386E+05
120.117794E+060.418256E+04-0.285486E+05
13-0.821554E+05-0.148737E+050.433242E+05
14-0.273330E+05-0.125220E+05-0.785096E+04
15-0.274155E+05-0.130158E+050.484341E+05
160.262559E+05-0.710079E+04-0.177296E+04
170.255681E+05-0.112192E+050.430881E+05
180.850804E+05-0.199044E+03-0.127141E+05
19-0.308284E+05-0.131056E+050.306629E+05
20-0.123014E+05-0.236279E+050.137015E+05
21-0.105684E+05-0.132504E+050.387534E+05
220.646544E+04-0.469412E+050.238659E+05
230.122400E+05-0.123633E+050.349753E+05
240.349927E+05-0.730929E+050.155524E+05
3、以例1为对象,研究单元细分对计算结果的影响。
3.1、单元细分如下所示:
图1:
单元细分图示
3.2、输入数据文件修改为:
15,16,30,10,1.0E0,0.0,1.0,0.0,1
3,1,2
5,2,4
3,2,5
6,3,5
8,4,7
5,4,7
9,5,8
6,5,9
10,6,9
12,7,11
8,7,12
13,8,12
9,8,13
14,9,13
10,9,14
15,10,14
0.0,2.0
0.0,1.5
0.5,1.5
0.0,1.0
0.5,1.0
1.0,1.0
0.0,0.5
0.5,0.5
1.0,0.5
1.5,0.5
0.0,0.0
0.5,0.0
1.0,0.0
1.5,0.0
2.0,0.0
1,3,7,13,21,22,24,26,28,30
3.3、输出数据文件修改为:
NNNENDNFIXEANUTGMNTYPE
151630100.1000E+010.0001.0000.0000E+001
NODEX-LOADY-LOAD
10.000000E+00-0.100000E+01
20.000000E+000.000000E+00
30.000000E+000.000000E+00
40.000000E+000.000000E+00
50.000000E+000.000000E+00
60.000000E+000.000000E+00
70.000000E+000.000000E+00
80.000000E+000.000000E+00
90.000000E+000.000000E+00
100.000000E+000.000000E+00
110.000000E+000.000000E+00
120.000000E+000.000000E+00
130.000000E+000.000000E+00
140.000000E+000.000000E+00
150.000000E+000.000000E+00
NODEX-DISPY-DISP
1-0.971690E-15-0.440129E+01
2-0.273631E-16-0.240129E+01
3-0.236860E+00-0.142960E+01
40.193766E-17-0.121245E+01
50.841994E-01-0.102402E+01
6-0.257230E-01-0.406152E+00
70.117892E-15-0.712244E+00
80.783953E-01-0.389010E+00
90.159794E+00-0.223862E+00
100.149665E+00-0.310062E-01
110.990200E-16-0.474829E-15
120.148530E+00-0.271029E-15
130.205578E+00-0.145183E-15
140.225835E+00-0.257352E-16
150.225835E+000.761703E-16
ELEMENTX-STRY-STRXY-STR
1-0.473719E+00-0.400000E+010.971690E+00
20.168399E+00-0.237768E+010.188429E+00
3-0.473719E+00-0.811160E+000.650631E+00
4-0.219845E+00-0.811160E+000.296808E+00
50.156791E+00-0.100041E+010.323234E+00
60.168399E+00-0.100041E+010.188429E+00
70.162797E+00-0.127002E+010.170951E+00
8-0.219845E+00-0.364580E+000.432351E+00
9-0.202577E-01-0.364580E+000.733912E-02
100.297060E+00-0.142449E+010.222672E-15
110.156791E+00-0.778019E+000.253099E+00
120.114095E+00-0.778019E+00-0.701347E-01
130.162797E+00-0.447725E+000.119364E+00
140.405153E-01-0.447725E+00-0.457838E-01
15-0.202577E-01-0.620125E-010.116686E+00
160.000000E+00-0.620125E-01-0.761703E-01
3.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-01
2
8.79E-17
-1.25E+00
4
1.94E-18
-1.21E+00
4.44E+01
3.32E-02
3
-8.79E-02
-3.74E-01
6
-2.57E-02
7.62E-17
2.42E+00
-4.91E+15
4
1.76E-01
2.64E-16
11
9.90E-17
-4.75E-16
1.78E+15
-1.56E+00
5
1.76E-01
-2.93E-16
13
2.26E-01
-1.45E-02
-2.21E-01
-1.00E+00
6
1.17E-16
-8.35E-16
15
2.06E-01
-4.06E-01
-1.00E+00
-1.00E+00
4、用程序完成习题3和4的分析。
课本第二章习题3,单元划分如下图所示:
图2:
习题2-2图
输入数据文件为:
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,4
8,4,7
0.0,0.0
1.0,0.0
2.0,0.0
0.0,1.0
1.0,1.0
2.0,1.0
0.0,3.0
1.0,3.0
0
输出数据文件为:
NNNENDNFIXEANUTGMNTYPE
861600.2000E+120.3000.1000.0000E+001
NODEX-LOADY-LOAD
10.000000E+000.000000E+00
20.000000E+000.000000E+00
30.707000E+03-0.707000E+03
40.000000E+000.000000E+00
50.000000E+000.000000E+00
60.000000E+000.000000E+00
7-0.707000E+030.707000E+03
80.000000E+000.000000E+00
NODEX-DISPY-DISP
1-0.177257E+000.493387E-01
2-0.177257E+00-0.258127E-01
3-0.177257E+00-0.100964E+00
4-0.102106E+000.493387E-01
5-0.102106E+00-0.258127E-01
6-0.102106E+00-0.100964E+00
70.481974E-010.493387E-01
80.481974E-01-0.258128E-01
[ELEMENTX-STRY-STRXY-STR
1-0.590520E+04-0.800872E+040.668666E+04
2-0.776742E+04-0.729983E+040.755878E+04
30.146578E+04-0.497822E+04-0.259785E+04
4-0.709961E+030.530443E+04-0.460698E+04
5-0.405635E+04-0.960508E+040.794023E+04
6-0.218427E+04-0.531371E+040.443640E+04
5、采用SAP2000分析上述算例。
(1)上面习题1,输出结果如下表所示:
表2:
上机题1spa2000计算的单元应力结果
TABLE:
ElementStresses-AreaPlanes
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/mm2
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
1
DEAD
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.0004396
0.000008296
-0.002096
12.346
0.0021
3
2
1
DEAD
LinStatic
0.0001758
-0.001253
0
0
0.0001758
-0.001253
0
0.001349
3
2
3
DEAD
LinStatic
0.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.0006313
4
3
4
DEAD
LinStatic
-0.00008791
-0.0003736
0
-0.0003077
0.0001085
-0.00057
-32.548
0.0006313
4
3
2
DEAD
LinStatic
-0.00008791
-0.0003736
0
-0.0003077
0.0001085
-0.00057
-32.548
0.0006313
5
4
4
DEAD
LinStatic
6.746E-20
-0.0003736
0
-0.0001319
0.00004185
-0.0004155
-17.609
0.0004379
5
4
2
DEAD
LinStatic
0
-0.0003736
0
-0.0001319
0.00004185
-0.0004155
-17.609
0.0004379
5
4
5
DEAD
LinStatic
-6.746E-20
-0.0003736
0
-0.0001319
0.00004185
-0.0004155
-17.609
0.0004379
表3:
上机题1spa2000计算的节点位移结果
TABLE:
JointReactions
Joint
OutputCase
CaseType
F1
F2
F3
M1
M2
M3
Text
Text
Text
N
N
N
N-mm
N-mm
N-mm
1
DEAD
LinStatic
-131.87
0
0
0
0
0
2
DEAD
LinStatic
0
0
439.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
对比前面fortran编程计算得出的结果知,两者的差距不大,有限单元法的计算精度较高。
(2)输出结果如下表所示:
表4:
习题2-2spa2000计算考察点输出结果
TABLE:
ElementStresses-AreaPlanes
Area
S11
S22
S33
S12
SMax
SMin
SAngle
SVM
Text
N/mm2
N/mm2
N/mm2
N/mm2
N/mm2
N/mm2
Degrees
N/mm2
1
0.244
0.008692
0
-0.019
0.246
0
-4.659
0.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.384
0.089
29
0.077
0.007621
0
-0.021
0.083
0
-15.886
0.082
30
-0.00205
0.000736
0
0.027
0.026
-0.027
46.488
0.046
31
-0.00132
0.005068
0
-0.023
0.025
-0.021
-48.952
0.04
32
-0.08
-0.00176
0
0.025
0.00553
-0.087
73.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
35
-0.158
-0.00045
0
-0.028
0.004322
-0.162
-80.256
0.165
36
-0.235
0.001969
0
0.018
0.00332
-0.237
85.696
0.238
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 有限元 上机 实验 报告