Galīgo elementu metode

2.7. GEM UZDEVUMA RISINĀŠANAS ALGORITMS 37
1020 DEFINT I-N: CLS
1030 PRINT ”=================================”
1040 PRINT ” Program written by ”
1050 PRINT ” T.R.Chandrupatla and A.D.Belegundu ”
1060 PRINT ”=================================”
1070 INPUT ” Number of Elements =”; NE
1080 INPUT ” Number of Loads =”; NL
1090 INPUT ”Num of Constr Nodes =”; ND
1100 NN = NE + 1
1110 NBW = 2
1120 DIM X(NN), A(NE), NU(ND), U(ND), DT(NE),AL(NL),ILD(NL)
1130 DIM S(NN, NBW), F(NN), E(NE),
ALP(NE),STRESS(NE),REACT(ND)
1140 PRINT ”Elem#, Area, E, Alpha, Temp Rise”
1150 FOR I = 1 TO NE
1160 INPUT N, A(N), E(N), ALP(N), DT(N): NEXT I
1170 PRINT ”Node#, Coordinate”
1180 FOR I = 1 TO NN
1190 INPUT N, X(N): NEXT I
1200 FOR I = 1 TO NN
1210 F(I) = 0
1220 FOR J = 1 TO NBW
1230 S(I, J) = 0: NEXT J: NEXT I
1240 PRINT ”Node#, Displacement”
1250 FOR I = 1 TO ND
1260 INPUT NU(I), U(I): NEXT I
1270 IF NL = 0 THEN GOTO 1330
1280 PRINT ”Node#, Applied Load”
1290 FOR I = 1 TO NL
1300 INPUT N, F(N)
1310 ILD(I)=N : AL(I)=F(N) : NEXT I
1320 REM *** STIFFNESS MATRIX ***
1330 FOR I = 1 TO NE
1340 I1 = I: I2 = I + 1
1350 X21 = X(I2) - X(I1): EL = ABS(X21)
1360 EAL = E(I) * A(I) / EL
1370 TL = E(I) * ALP(I) * DT(I) * A(I) * EL / X21
1380 REM *** TEMPERATURE LOADS ***
1390 F(I1) = F(I1) - TL
1400 F(I2) = F(I2) + TL
1410 REM *** ELEMENT STIFFNESSES ***