DataMuseum.dk

top - download

    Length: 9728 (0x2600)
    Types: TextFile
    Names: »B28«

└─⟦6dbcc9c03⟧ Bits:30000463 DOMUS disk image
    └─ ⟦this⟧ »/B28« 

TextFile

 0010 REM ******** TORBEN GRAVESEN 3U BOX388 *******
 0020 PRINT CHR(12)
 0030 DEF FNC(Ø)=95+Ø-((Ø-1) DIV 32)*64
 0040 DEF FIX(ØØ)=INT(ØØ*10)/10
 0100 PRINT "     DIMMENTIONERING OG STYRING AF EN REAKTOR"
 0110 PRINT 
 0120 PRINT "DETTE POGRAM ER BEREGNET TIL BRUG SAMMEN MED FINN@S ELLER ØLGAARD@S"
 0130 PRINT "NOTER I REAKTORFYSIK  PROGRAMMET GIVER MULIGHED FOR AT"
 0140 PRINT "DIMMENTIONERE OG STYRE EN HOMOGEN REAKTOR,HVOR DER BENYTTES U(235)"
 0150 PRINT "SOM BRÆNDSELS MATERIALE (1)  ELLER DU KAN VÆLGE EN FÆRTIG"
 0160 PRINT "REAKTOR OG BLOT STYRE DEN (2)"
 0170 INPUT " VÆLG 1 ELLER 2 ",VALG
 0200 IF VALG=2 THEN 
 0210   R=2.33; ALFA=2000; S=1.3529; EP=1; P=.9; LM=.05; TAU=.035; AM=.0045
 0215   V=2.47; AF=582; AA=683; AU=2.75; MMM=12; MFM=1.65
 0220   ETA=(V*AF)/(AA+(R*AU))
 0225   B=3*((3.141593/S)^2); V=S*S*S
 0230   PF=EXP(-B*TAU)
 0280   GO TO 0450
 0300 ENDIF 
 0410 EXEC N
 0420 EXEC TU
 0430 EXEC NEUTRON
 0440 EXEC SILON
 0450 EXEC FIREFAKTOR
 0460 EXEC STYR
 0470 GO TO 0450
 0480 END 
 2000 PROC MASSEBEREGN
 2002   IF V=0 THEN 
 2004     PRINT " DET KAN IKKE OPLYSES,DA DU HAR VALGT EN UENDELIG STOR REAKTOR"
 2005     PRINT 
 2006   ENDPROC 
 2007   PRINT 
 2008 ENDIF 
 2010 MILFOR238=R*1.02766
 2020 MILFORMOD=ALFA*MMM*18.65/(MFM*235)
 2030 FORHOLD=MILFOR238+1+MILFORMOD
 2040 VOLUME235=V/FORHOLD
 2050 VOLUME238=MILFOR238*V/FORHOLD
 2060 VOLUMEMOD=MILFORMOD*V/FORHOLD
 2070 REM 
 2080 MASSE235=18.86*VOLUME235
 2090 MASSE238=18.65*VOLUME238
 2100 MASSEMOD=MFM*VOLUMEMOD
 2102 PRINT "FØLGENDE MÆNGDE BRÆNDSELS OG MODERATOR MATERIALE"
 2104 PRINT 
 2110 PRINT "VOLUME AF URAN(235) ER"; VOLUME235; "KUBIKMETER DET VEJER"; 
 2120 PRINT FIX(MASSE235*1000); "KG"
 2130 PRINT "VOLUME AF URAN(238) ER"; VOLUME238; "KUBIKMETER DET VEJER"; 
 2140 PRINT FIX(MASSE238*1000); "KG"
 2150 PRINT "VOLUMEN AF MODERATOR ER"; VOLUMEMOD; "KUBIKMETER DET VEJER"; 
 2160 PRINT FIX(MASSEMOD*1000); "KG"
 2165 PRINT 
 2170 ENDPROC MASSEBEREGN
 5000 PROC N
 5004   PRINT "   START MED AT VÆLGE DIN "
 5006   INPUT "  BERIGNINGSGRAD I % : ",R
 5007   R=(100-R)/R
 5008   PRINT 
 5010   PRINT " ETA ER DEFINERET SOM ANTALLET AF FISSIONS NEUTRONER"
 5020   PRINT " DER UDSENDES PR TERMISK NEUTRON ABSORBERET I URANSTAVENE"
 5025   PRINT 
 5030   V=2.47; AF=582; AA=683; AU=2.75; TAU=0.081
 5040   REM SIGMA(F)(235)=AF   SIGMS(A)(235)=AA   SIGMA(A)(238)=AU
 5050   ETA=(V*AF)/(AA+(R*AU))
 5060   PRINT " ANTAL FISSIONS NEUTRONER (ETA) : "; ETA
 5070 ENDPROC 
 6000 PROC TU
 6004   PRINT 
 6006   PRINT 
 6010   PRINT " F ER DEN TERMISKE UDNYTTELSESFAKTOR "
 6015   PRINT 
 6020   PRINT "    N(235)*SIGMA(A)(235)+N(238)*SIGMA(A)(238)"
 6030   PRINT " F= -------------------------------------------"
 6040   PRINT "    SIGMA(A)(235)+R*SIGMA(A)(238)+ALFA*SIGMA(A)(M)"
 6045   PRINT 
 6048   PRINT "  VÆLG MODERATOR MATERIALE DER ER 5 MULIGHEDER  "
 6050   PRINT "   1  ALM. VAND     H2O"
 6052   PRINT "   2  TUNGT VAND    D2O"
 6054   PRINT "   3  GRAFIT        C"
 6056   PRINT "   4  BERYLLIUM     BE"
 6058   PRINT "   5  ANDET"
 6060   INPUT " VÆLG (1-5): ",VALG
 6062   IF VALG=1 THEN LM=.0288; TAU=.0033; AM=.66; MMM=18; MFM=1
 6064   IF VALG=2 THEN LM=.1; TAU=.012; AM=.0026; MMM=20; MFM=1.10
 6066   IF VALG=3 THEN LM=.05; TAU=.035; AM=.0045; MMM=12; MFM=1.65
 6068   IF VALG=4 THEN LM=.024; TAU=.0098; AM=.009; MMM=9; MFM=1.85
 6070   IF VALG=5 THEN 
 6072     INPUT "  SKRIV TAU , LM , SIGMA(A)(M)  :",TAU,LM,AM
 6074   ENDIF 
 6076   PRINT 
 6079   INPUT "  ALFA=N(M)/N(235)  SKRIV ALFA ! (F.EKS 2000) : ",ALFA
 6080   F=(AA+R*AU)/(AA+R*AU+ALFA*AM)
 6090   PRINT "    F = "; F
 6100   PRINT 
 6110 ENDPROC TU
 7000 PROC SILON
 7010   PRINT 
 7020   PRINT 
 7030   PRINT " FIREFAKTORFORMEL :  K=ETA*EPSILON*P*F*PF*PT "
 7040   PRINT 
 7050   INPUT "  VÆLG EPSILON MELLEM 1 OG 1.04 : ",EP
 7060   INPUT "  OG VÆLG P MELLEM 0.75 OG 0.9 : ",P
 7062 ENDPROC SILON
 7068 PROC FIREFAKTOR
 7070   F=(AA+R*AU)/(AA+R*AU+ALFA*AM); L=LM*LM*(1-F); PT=1/(1+B*L)
 7080   K=ETA*EP*F*P*PF*PT
 7090   PRINT "      K = ",K
 7100   PRINT 
 7110   IF K<1 THEN 
 7120     PRINT " DENNE REAKTOR VIL ALDRIG KOMME I GANG "
 7125     GO TO 7200
 7130   ENDIF 
 7140   IF K>1.1 THEN 
 7150     PRINT " DENNE REAKTOR VIL NEDSMELTE ILØBET AF KORT TID (KINASYNDROMET ! )"
 7155     GO TO 7200
 7160   ENDIF 
 7170   PRINT " DENNE REAKTOR HAR GODE FREMTIDS MULIGHEDER "
 7180   INPUT "   ØNSKER DU AT FORBEDRE REAKTOREN?JA=1 NEJ=2 ",VALG
 7185   IF VALG=2 THEN GO TO 7450
 7190   GO TO 7210
 7200   PRINT "   JEG VIL RÅDE DIG TIL AT ÆNDRE DIN REAKTOR"
 7210   PRINT "   DU KAN VÆLGE AT ÆNDRE PÅ FØLGENDE OMRÅDER"
 7220   PRINT "     1   HELE DIMENTIONERINGEN (STARTE FORFRA)"
 7230   PRINT "     2   BERIGNINGSGRADEN"
 7240   PRINT "     3   DEN TERMISKEUDNYTTELSES FAKTOR F OG MODERATOR MATERIALE"
 7250   PRINT "     4   REAKTORENS STØRRELSE OG FORM"
 7260   PRINT "     5   EPSILON OG P"
 7270   PRINT "     6   INGEN ÆNDRING"
 7280   REM 
 7290   INPUT "  DIT VALG (1-6) : ",VALG
 7300   ON VALG GO TO 7310,7350,7370,7400,7420,7450
 7310   EXEC N
 7320   EXEC TU
 7325   EXEC SILON
 7330   EXEC NEUTRON
 7340   GO TO 7070
 7350   EXEC N
 7360   GO TO 7070
 7370   EXEC TU
 7380   L=LM*LM*(1-F); PT=1/(1+B*L)
 7390   GO TO 7070
 7400   EXEC NEUTRON
 7410   GO TO 7070
 7420   EXEC SILON
 7430   GO TO 7070
 7450 ENDPROC FIREFAKTOR
 8000 PROC NEUTRON
 8010   PRINT 
 8020   PRINT " DENNE DEL AF PROGRAMMET UDREGNER FORSKELLIGE VÆRDIER "
 8030   PRINT " FOR PF OG PT AFHÆNGIGT AF REAKTORENS FORM"
 8035   PRINT " DU STARTER MED AT VÆLGE REAKTORENS FORM, DER ER 4 MULIGHEDER "
 8040   PRINT " 1  SFÆRISK REAKTOR MED RADIUS R "
 8050   PRINT " 2  CYLINDRISK REAKTOR MED HØJDEN H OG RADIUS R"
 8060   PRINT " 3  KUBISK REAKTOR MED KANTLÆNGDEN S"
 8065   PRINT " 4 UENDELIG STOR REAKTOR PF OG PT ER LIG 1"
 8070   INPUT " VÆLG REAKTOR FORM 1,2,3 ELLER 4 : ",VALG
 8080   ON VALG GO TO 8100,8200,8300,8250
 8100   PRINT 
 8110   INPUT "   SKRIV RADIUS R : ",R
 8120   B=(3.141593/R)^2; V=4*3.14159*R*R*R/3
 8130   GO TO 8400
 8200   PRINT 
 8210   INPUT "   SKRIV HØJDEN H ",H,"   OG RADIUS R ",R
 8220   B=(2.405/R)^2+(3.141593/H)^2
 8230   GO TO 8400
 8250   PF,PT=1; B=0; V=0
 8260   GO TO 8570
 8300   PRINT 
 8310   INPUT "   SKRIV KANTLÆNGDEN S : ",S
 8320   B=3*((3.141593/S)^2); V=S*S*S
 8400   PF=EXP(-B*TAU); L=LM*LM*(1-F); PT=1/(1+B*L)
 8430   PRINT 
 8440   PRINT "      PF= "; PF; "   PT= "; PT
 8450   PRINT 
 8570 ENDPROC NEUTRON
 9000 PROC STYR
 9001   PRINT 
 9002   INPUT "  ØNSKER DU AT FÅ OPLYST  BRÆNSELS MÆNGDE? JA=1",VALG
 9004   IF VALG=1 THEN EXEC MASSEBEREGN
 9005   Z=100; TID=0; DT=.2; ALF=.02; TAV=.081
 9010   PRINT " DU HAR NU DIMENTIONERET DIN REAKTOR OG SKAL VÆLGE MATERIALE "
 9020   PRINT " TIL DINE STYRINGS KONTROL STÆNGER ! "
 9030   PRINT 
 9040   PRINT " DU HAR FØLGENDE VALG MULIGHEDER :"
 9050   PRINT " 1  BOR (B)"
 9060   PRINT " 2  CADMIUM (CD)"
 9070   PRINT " 3  EUROPIUM (EU)"
 9080   PRINT " 4  ANDET MATERIALE"
 9090   INPUT " DIT VALG !(1,2,3 ELLER 4) : ",VALG
 9100   IF VALG=1 THEN SIGMAS=755
 9110   IF VALG=2 THEN SIGMAS=2450
 9120   IF VALG=3 THEN SIGMAS=4300
 9130   IF VALG=>4 THEN 
 9140     INPUT " SKRIV ABS.TVÆRSNITTET FOR DIT MATERIALE : ",SIGMAS
 9150   ENDIF 
 9160   PRINT " DU KAN NU STARTE REAKTOREN OP,MEN FØRST SKAL DU OPLYSE"
 9170   PRINT " HVOR MANGE PROCENT AF DINE KOMTROLSTÆNGER DU VIL HAVE "
 9180   INPUT " I REAKTOREN : ",PRC
 9190   ALFAS=PRC*ALF/100
 9300   PRINT " DU KAN NU STYRE DIN REAKTOR VED AT FLYTTE STÆNGERNE OP OG NED "
 9310   PRINT " MENS REAKTOREN KØRE DETTE GØRES FRA TASTATURET "
 9320   PRINT " TASTES  '1' FLYTTES OPAD OG BLIVER VED MED DET INDTIL DU "
 9330   PRINT " TASTER  '2' FOR STOP ELLER '3' FOR NED "
 9340   F=(AA+R*AU)/(AA+R*AU+ALFA*AM+ALFAS*SIGMAS); L=LM*LM*(1-F); PT=1/(1+B*L)
 9345   K=ETA*EP*F*P*PF*PT
 9350   DN=EXP((K-1)/TAV)
 9360   PRINT " NEUTRONERNE FORØGES NU MED FAKTOREN PR SEK : "; DN; " K="; K
 9370   INPUT " VIL DU STARTE ? (JA=1 NEJ=2 OMDIMENTIONERING=3 ) ",VALG
 9375   PRINT 
 9378   IF VALG=3 THEN GO TO 9900
 9380   IF VALG<>1 THEN GO TO 9160
 9390   PRINT CHR(12)
 9400   REPEAT KØRSEL
 9405     TID=TID+DT
 9410     IN 132,OPNED
 9412     IF OPNED=177 AND PRC>0 THEN PRC=PRC-1
 9414     IF OPNED=179 AND PRC<100 THEN PRC=PRC+1
 9416     IF CHR(OPNED)="H" THEN PRC=100
 9418     IF OPNED=177 AND PRC<0 THEN PRC=0
 9420     ALFAS=ALF*PRC/100
 9450     F=(AA+R*AU)/(AA+R*AU+ALFA*AM+ALFAS*SIGMAS); L=LM*LM*(1-F); PT=1/(1+B*L)
 9500     K=ETA*EP*F*P*PF*PT
 9535     Z=Z*EXP((K-1)*DT/TAV)
 9540     EXEC CURSOR(5,10)
 9550     PRINT " DU HAR NU "; INT(Z-100); " PROCENT FLERE NEUTRONER "
 9570     EXEC CURSOR(5,3)
 9580     PRINT "   TID(SEK) : "; INT(TID*10)/10
 9582     PRINT 
 9584     PRINT "  PRCENT  AF KONTROLST. I REAKTOREN : "; PRC
 9586     PRINT 
 9588     PRINT "  K : "; K
 9800   UNTIL CHR(OPNED)="S" OR INT(Z)<=0 OR INT(Z)>200
 9810   PRINT CHR(12)
 9900 ENDPROC STYR
 9910 PROC CURSOR(X,Y)
 9920   PRINT CHR(6); CHR(FNC(X)); CHR(FNC(Y))
 9930 ENDPROC CURSOR
«eof»