有限元08上机作业
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
调试书本26到30页程序
开列数组维数:
DIMENSION LOC(4,3),CX(6),CY(6),IFIX(6),F(12),
1GK(12,12),STRES(4,3),BAK(4,3,6)
结点集中力输入:
DO 10 I=1,ND
10 F(I)=0.0
F(2)=-1.0
数据文件输入:
6,4,12,6,1.0E0,0.0,1.0,0.0,1
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
数据文件输出:
NN NE ND NFIX E ANU T GM NTYPE
6 4 12 6 0.1000E+01 0.000 1.000 0.0000E+00 1
NODE X-LOAD Y-LOAD
1 0.000000E+00 -0.100000E+01
2 0.000000E+00 0.000000E+00
3 0.000000E+00 0.000000E+00
4 0.000000E+00 0.000000E+00
5 0.000000E+00 0.000000E+00
6 0.000000E+00 0.000000E+00
NODE X-DISP Y-DISP
1 -0.879121E-15 -0.325275E+01
2 0.879121E-16 -0.125275E+01
3 -0.879121E-01 -0.373626E+00
4 0.117216E-1
5 -0.835165E-15
5 0.175824E+00 -0.293040E-15
6 0.175824E+00 0.263736E-15
ELEMENT X-STR Y-STR XY-STR
1 -0.879121E-01 -0.200000E+01 0.439560E+00
2 0.175824E+00 -0.125275E+01 0.256410E-15
3 -0.879121E-01 -0.373626E+00 0.307692E+00
4 0.000000E+00 -0.373626E+00 -0.131868E+00 通过Sap2000进行结果检验,结果符合.
计算书本31页例题2-2
开列数组维数:
DIMENSION LOC(144,3),CX(91),CY(91),IFIX(8),F(182),
1GK(182,182),STRES(144,3),BAK(144,3,6)
数据读入程序:
OPEN(5,FILE='E2.DAT',STATUS='OLD')
OPEN(6,FILE='OUT2.')
READ(5,*)NN,NE,ND,NFIX,E,ANU,T,GM,NTYPE
WRITE(6,105)NN,NE,ND,NFIX,E,ANU,T,GM,NTYPE
105 FORMAT(2X,'NN NE ND NFIX E ANU T GM
1 NTYPE'/4I4,2X,E10.4,2F7.3,2X,E10.4,I3)
READ(5,*)(IFIX(K),K=1,NFIX)
单元编码程序:
LOC(1,1)=2
LOC(1,2)=8
LOC(1,3)=1
LOC(2,1)=8
LOC(2,2)=2
LOC(2,3)=9
DO 25 I=1,10
LOC(I+2,1)=LOC(I,1)+1
LOC(I+2,2)=LOC(I,2)+1
25 LOC(I+2,3)=LOC(I,3)+1
DO 35 I=1,132
LOC(I+12,1)=LOC(I,1)+7
LOC(I+12,2)=LOC(I,2)+7
35 LOC(I+12,3)=LOC(I,3)+7
节点坐标程序:
DO 45 J=1,13
DO 45 I=1,7
CX(I+7*J-7)=0.75*J-0.75
45 CY(I+7*J-7)=3.5-0.5*I
结点集中力输入:
DO 10 I=1,ND
10 F(I)=0.0
F(2)=-3.75
F(170)=-3.75
DO 15 I=8,78,7
15 F(2*I)=-7.5
数据文件输入:
91,144,182,8,2.0E10,0.167,1.0,0.0,1
1,3,5,7,9,11,13,182
数据文件输出:
NN NE ND NFIX E ANU T GM NTYPE
91 144 182 8 0.2000E+11 0.167 1.000 0.0000E+00 1
(结点荷载太多,再此不显示:)
NODE X-DISP Y-DISP(只显示最大最小位移结点)
91 0.766665E-07 -0.492018E-23(x方向位移最大)
85 -0.694272E-07 -0.148263E-07(x方向位移最小)
90 0.476889E-07 -0.558120E-08(y方向位移最小)
4 -0.570359E-2
5 -0.294638E-06(y方向位移最大)
5 0.141856E-23 -0.294665E-06(y方向位移最大)
ELEMENT X-STR Y-STR XY-STR
12 0.244167E+03 0.868761E+01 -0.193180E+02(x方向应力最大)
13 -0.236666E+03 -0.136781E+02 0.267470E+02(x方向应力最小)
11 0.157851E+03 0.128788E+02 0.312907E+02(y方向应力最大)
144 0.601269E+02 -0.213207E+03 0.400848E+02(y方向应力最小)142 -0.141647E+01 -0.158214E+03 0.533252E+02(xy剪应力最大)
2 -0.157505E+0
3 -0.457531E+00 -0.277912E+02(xy剪应力最小)