土木华科大有限元大作业Word文件下载.docx
《土木华科大有限元大作业Word文件下载.docx》由会员分享,可在线阅读,更多相关《土木华科大有限元大作业Word文件下载.docx(31页珍藏版)》请在冰豆网上搜索。
3,1,2
5,2,4
2,5,3
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
数据文件输出:
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-1
5-0.835165E-1550.175824E+00-0.293040E-15
60.175824E+000.263736E-15
ELEMENTX-STRY-STRXY-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+00调试成功,与sap2000结果相同2.修改FEM1,计算P31页例2.2.
例2-2.如图所示的简支梁,梁高w=3m,长2l=18m,厚度h=1m,承受均布载荷
q=10kN/m2。
已知102E=210/m,0.167Nμ⨯=。
按平面应力问题进行求解。
受均布荷载作用的简支梁
由于对称性,取梁右半边进行有限元分析,如下图所示。
划分144个单元,91个节点。
简支梁的有限元网格划分(1、建立数据文件E2.DAT。
E2.DAT
91,144,182,8,2.0E10,0.167,1.0,0.0,12,8,19,8,23,9,210,9,34,10,311,10,45,11,412,11,56,12,513,12,6
7,13,614,13,7
9,15,8
16,15,910,16,917,16,1011,17,10
18,17,1112,18,1119,18,1213,19,1220,19,1314,20,1321,20,1416,22,1523,22,1617,23,1624,23,1718,24,1725,24,1819,25,1826,25,1920,26,1927,26,1021,27,2028,27,2123,29,2230,29,2324,30,2331,30,2425,31,2432,31,2526,32,2533,32,2627,33,2634,33,2728,34,2735,34,2830,36,2937,36,3031,37,3038,37,3132,38,3139,38,3233,39,3240,39,3334,40,3341,40,3435,41,3442,41,3537,43,3644,43,37
38,44,37
45,44,38
39,45,38
46,45,39
40,46,39
47,46,40
41,47,40
48,47,41
42,48,41
49,48,42
44,50,43
51,50,44
45,51,44
52,51,45
46,52,45
53,52,46
47,53,46
54,53,47
48,54,47
55,54,48
49,55,48
56,55,49
51,57,50
58,57,51
52,58,51
59,58,52
53,59,52
60,59,53
54,60,53
61,60,54
55,61,54
62,61,55
56,62,55
63,62,56
58,64,57
65,64,58
59,65,58
66,65,59
60,66,59
67,66,60
61,67,60
68,67,61
62,68,61
69,68,62
63,69,62
70,69,63
65,71,64
72,71,65
66,72,65
73,72,66
67,73,66
74,73,67
68,74,67
75,74,68
69,75,68
76,75,69
70,76,69
77,76,70
72,78,71
79,78,72
73,79,72
80,79,73
74,80,73
81,80,74
75,81,74
82,81,75
76,82,75
83,82,76
77,83,76
84,83,77
79,85,78
86,85,79
80,86,79
87,86,80
81,87,80
88,87,81
82,88,81
89,88,82
83,89,82
90,89,83
84,90,83
91,90,84
0.0,3.0
0.0,2.5
0.0,1.5
0.0,0.5
0.75,3.0
0.75,2.5
0.75,2.0
0.75,1.5
0.75,1.0
0.75,0.5
0.75,0.0
1.5,3.0
1.5,
2.5
2.0
1.5,1.5
1.5,1.0
1.5,0.5
1.5,0.0
2.25,
3.0
2.25,2.5
2.25,2.0
2.25,1.5
2.25,1.0
2.25,0.5
2.25,0.0
3.0,3.0
3.0,2.5
3.0,2.0
3.0,1.5
3.0,1.0
3.0,0.5
3.0,0.0
3.75,3.0
3.75,2.5
3.75,2.0
3.75,1.5
3.75,1.0
3.75,0.5
3.75,0.0
4.5,3.0
4.5,2.5
4.5,2.0
4.5,1.5
4.5,1.0
4.5,0.5
4.5,0.0
5.25,3.05.25,2.55.25,2.05.25,1.55.25,1.05.25,0.5
5.25,0.0
6.0,3.06.0,2.56.0,2.06.0,1.56.0,1.0
6.0,0.5
6.0,0.0
6.75,3.0
6.75,2.5
6.75,2.0
6.75,1.5
6.75,1.0
6.75,0.5
6.75,0.0
7.5,3.0
7.5,2.5
7.5,2.0
7.5,1.5
7.5,1.0
7.5,0.5
7.5,0.0
8.25,3.0
8.25,2.5
8.25,2.0
8.25,1.5
8.25,1.0
8.25,0.5
8.25,0.0
9.0,3.0
9.0,2.5
9.0,2.0
9.0,1.5
9.0,1.0
9.0,0.5
9.0,0.0
1,3,5,7,9,1
1,13,182
(2、扩大源程序主程序段中有关数组的维数,将分布荷载转化为等效节点荷载,在源程序的主程序中直接输入,从而形成新的源程序FEM2。
(3、编辑,连接,形成FEM2.EXE。
运行,形成计算结果文件OUT2.TXT。
911441828.20000E+11.167001.00000.00000E+001
1.000000E+00-.375000E+09
2.000000E+00.000000E+00
3.000000E+00.000000E+00
………..
………..(略去
85.000000E+00-.375000E+09
86.000000E+00.000000E+00
87.000000E+00.000000E+00
88.000000E+00.000000E+00
89.000000E+00.000000E+00
90.000000E+00.000000E+00
91.000000E+00.000000E+00
1-.494173600E-25-.296178800E+02
2-.302584000E-25-.297209800E+02
3-.156707600E-25-.297803200E+02
4-.248758100E-27-.298067200E+02
5.147426300E-25-.297991000E+02
6.287587200E-25-.297723200E+02
85-.692795800E+01-.148261300E+01
86-.455980100E+01-.145083100E+01
87-.225408900E+01-.137430300E+01
88.364165100E-01-.121739900E+01
89.238245500E+01-.953019100E+00
90.489660900E+01-.558089600E+00
91.781694000E+01-.491989900E-25
ELEMENTX-STRY-STRXY-STR
1-.235965000E+11.183573300E+09.177544800E+10
2-.157376100E+11.865337500E+07-.286918400E+10
3-.157828200E+11-.262075100E+09.220028800E+10
4-.813478600E+10.397898800E+09-.246771900E+10
5-.825516400E+10-.322907500E+09.237660400E+10
6.134268100E+09-.873846200E+0