Programa DINÂMICA

1 CLS
2 BORDER 0: PAPER 0: INK 7: CLS
3 PRINT "INSTITUTO POLITÉCNICO DE LISBOA": PRINT
4 PRINT "INSTITUTO"
5 PRINT "SUPERIOR"
6 PRINT "DE"
7 PRINT "ENGENHARIA"
8 PRINT "DE"
9 PRINT "LISBOA"
10 PRINT : PRINT
11 PRINT "PROGRAMA: DINÂMICA"
12 PRINT
13 PRINT "Método de Holzer"
14 PRINT
15 PRINT "AUTORES: Copyright (c) 1993": PRINT
16 PRINT "Maria Helena S. J. Cunha Lopes"
17 PRINT "Eduardo F. P. Cunha Lopes J."
18 PRINT : PRINT
19 PRINT "Centro de Estudos de Engenharia Civil/Centro de Calculo - ISEL"
20 PAUSE O: CLS
30 PRINT AT 10,28;"t(s)"
40 PRINT AT 3,1;"f(Hz)"
50 PLOT 50,50: DRAW 0,100
60 PLOT 50,100: DRAW 200,0
70 PLOT 50,100: DRAW 5,-40: DRAW 5,80: DRAW 5,-90: DRAW 5,100: DRAW 5,-80: DRAW 5,95:DRAW 5,-70: DRAW 5,40: DRAW 5,-50: DRAW 5,20: DRAW 5,-4O: DRAW 5,70: DRAW 5,-65: DRAW 5,60: DRAW 5,-55: DRAW 5,5O DRAW 5,-45: DRAW 5,4O: DRAW 5,-35: DRAW 5,30: DRAW 5,-20;; DRAW 5,10: DRAW 5,-5: DRAW 5,0
80 INPUT "Numero de pisos "; n
90 INPUT "Modos de vibração pretendidos "; m
100 DIM m(n): DIM K(n): DIM V(n): DIM X(n): DIM F(n): DIM D(n): DIM W(m): DIM S(m): DIM T(m): DIM c(m)
105 INPUT "Direcção da análise (XX/YY)";D$
110 FOR i=1 TO n
115 PRINT AT 18,0; INVERSE 1; "PISO ";i
120 INPUT "Massa do piso (daKN*seg-2/cm) ",m(i)
130 INPUT "Rigidez do piso (daKN/cm) ",K(i)
140 NEXT i
145 PRINT AT 18,0;"
150 FOR j=1 TO m
151 INPUT "Arbitre a frequência W2 (seq-2) ";S(j)
152 GO TO 154
153 INPUT "Arbitre novamente W2 (seq-2) ",S(j)
154 CLS
155 PRINT TAB 8;"X";TAB 14;"DX";TAB 20;"V";TAB 26;"F": PRINT AT 2,O;"PISO 1"
160 LET X (1) = 1: PRINT AT 2,8;INT ( X (1)*10/10) / 10: LET D(1) = 1: PRINT AT 2,14;INT(D(1)*10 )/10
170 LET V(1)=K(1)*D(1): PRINT AT 2,20; INT (V(1)*10/10
180 LET F(1)=m(1)*S(j)*X(1) : PRINT AT 2,26; INT ( F(1) *10)/10
190 FOR i=2 TO n
200 LET V(i)=V (i -1)-F ( i -1)
210 LET D ( i) =V( i) /K (i)
220 LET X ( i)=X ( i -1)+D ( i)
230 LET F ( i ) =m( i ) *S ( j ) *X ( i )
235 PR I NT " PISO " ; i ; TAB 8 ; I NT ( X ( i ) * 10) /10; TAB 14; I NT ( D ( i) * l O) / l O ; TAB 20 ; I NT V( i)*10)/10; TAB 26;INT ( F( i) *10) / 10
240 NEXT i
250 LET r=V(n)-F(n)
260 PRINT : PRINT "erro= ";r: PAUSE O
270 IF ABS (r)>0.5 THEN GO TO 153
280 LET W(j)-SQR (S(j))
290 LET T ( j ) =2*PI /W(j)
300 LET c(j)=1/T(j)
310 PRINT : PRINT "Frequência f= ";INT ( c ( j)*100)/100;" Hz"
320 NEXT j
340 INPUT "DESEJA IMPRIMIR ? (S/N) ";I$
350 IF I$="n" OR I$="N" THEN GO TO 360
355 INPUT "DESEJA TÍTULO ? (S/N) "; T$
356 IF T$="n" OR T$="N" THEN GO TO 3500
357 IF T$="s" OR T$="S" THEN GO TO 400
360 INPUT "DESEJA NOVO CALCULO ? (S/N) ";N$
370 IF N$="s" OR N$="S" THEN GO TO 1
380 INPUT "CONFIRMA A OPÇÃO ? ; S/N) ";0$
330 IF 0$="s" OR 0$="S" THEN NEW
400 LPRINT "*********************************************************" : LPRINT
500 FORMAT LPRINT "U"
600 LPRINT , CHR$ 14;
700 LPRINT " S E I S M I C A N A L Y S I S "
800 FORMAT LPRINT "E"
900 LPRINT , CHR$ 10;
1OOO LPRINT : LPRINT : LPRINT "**************************************************": LPRINT: LPRINT
1500 FORMAT LPRINT "U"
2000 LPRINT ,CHR$ 14;
2500 LPRINT "DYNAMICS - HOLZER METHOD"
3000 FORMAT LPRINT "E"
3100 LPRINT ,CHR$10;
3200 LPRINT : LPRINT : LPRINT
3500 LPRINT "-------------------"
4000 LPRINT "DIRECTION - "; D$
4500 LPRINT "-------------------"
5000 LPRINT
5100 LPRINT TAB 10;"Mass";TAB 20;"Stiffness": LPRINT
5200 FOR i=1 TO n
5300 LPRINT "LEVEL ";i;TAB l O;m (i);TAB 20 ; K(i)
5400 NEXT i
5500 LPRINT : LPRINT "MODES OF VIBRATION FREQUENCIES"
6000 LPRINT
6500 FOR j=1 TO m
7000 LPRINT "MODE ";j;" : f= ";INT (c(j)*100)/100;" Hz"
7500 LPRINT
8000 NEXT j
9000 RUN
9900 GO TO 1
9999 SAVE "DINÂMICA" LINE 1

Copyright (c) 1993 Eduardo Cunha Lopes