Пиротехническая химия
Главная Начинающим пиротехникам Статьи Добавить статью Добавить материалы на сайт Поиск по сайту Карта книг Карта сайта
Книги в помощь
Военная история Изготовление и применение ВВ Пиротехника в военном деле Разное по пиротехнике Физика в пиротехнике Химия ВВ и составов
Новые книги
Суворов С. "Бронированная машина пехоты БМП -3 часть 1" (Военное дело)

Яковлев Г.П. "122 мм самоходная пушка образца 1944 г." (Военное дело)

Суарес Г. "Тактическое преимущество " (Военное дело)

Стодеревский И.Ю. "Автобиография записки офицера спецназа ГРУ " (Военное дело)

Семиколенков Н.П. "стрельба из танковых пулеметов " (Военное дело)
Применение ЭВМ для термодинамических расчетов металлургических процессов - Синярев Г.Б.
Синярев Г.Б., Ватолин Н.А., Трусов Б.Г., Моисеев Г.К. Применение ЭВМ для термодинамических расчетов металлургических процессов . Под редакцией Щепкин А.А. — М.: Наука, 1982. — 267 c.
Скачать (прямая ссылка): primenenevm1982.djvu
Предыдущая << 1 .. 82 83 84 85 86 87 < 88 > 89 90 91 92 93 94 .. 103 >> Следующая

Q ************* ************ *** **
N1=N+1
DO 1 1=1,N
1 M(I)=I NN=N-I
Q*«**
DO 8 K=I,NN КК=К+1
II (К•GT•NCH) GO ТО 7 W=O.
DO 3 J=K1NCH DO 3 I=K1NCH Z=DABS(A(I1J)) II (W-Z) 2,3,3
2 W=Z К1=1 К 2=J
3 CONTINUE
IF (К.EQ.К1) GO ТО 5 DO 4 J=K,N1 S=A(K,J) A(K,J)=A(K1 ,J)
4 A(KI1J)=S
5 IF (K.EQ.K2) GO ТО 7 DO 6 1=1,N S=A(I1K)
A(I,K)=A(I,К2)
6 A(I,K2)=S КЗ=М(К)
М(К)=М(К2) М(К2)=КЗ
7 DO 8 I=KK1N Z=A(I, К)/А(К,К) DO 8 J=K,N1
8 A(I,J)=A(I,J)-Z*A(K,J)
Q**«*
DO 11 J=I1N NN=J-I I=NW S=O.
II (N.EQ.I) GO ТО 10 DO 9 K=I9NN KI=NI-K
9 S=A(KI9N1)*A(I,K1)+S
10 КЗ=М(1)
A(I,N1) = (A(I,N1)-S)/A(I,( * I)
11 T(K3)=A(I,N1) RETURN
END
SUBROUTINE RECODE (E9R) COMMON /SYMB/K(94),A(44) DIMENSION R(9)
Q******************************
С ФОРМАТНОЕ ПРЕОБРАЗОВАНИЕ
С ДЕЙСТВИТЕЛЬНОГО ЧИСЛА
с* *****************************
R(D=AO»
IF (E.LT.O.O) R(1)=A(27)
F=ABS(E)
R(7)=A(26)
IF (I•LT.0.1) R(7)=A(27) N=O
IF (F.EQ.O) GO TO 5
1 IF (F-O.0999996) 2,3,3
2 I=F*10.0 N=N+T
GO TO 1
3 IF (F-O.999996) 5,4,4
4 F=F*0.1
212
с****
N=N+1
GO ТО З
N1=N/10
B(8)=A(N1+35)
N=N-N1*10
B(9)=A(N+35)
F=F+0.5E-5
DO 6 I=I,5
B=F*10.0
N=INT(B)
F=B-FLOAT(N)
B(I+1)=A(N+35)
RETURN
E(ID
SUBROUTINE HEABT (U9UH) DOUBLE PRECISION AC1RZ9T1
* A9STR9X9EX9S9H9CP9Y9EK9
* BO9U9W9XI9EPS9AI9A29XL9
* W
COUUOH /GEN/KDB(11),NC(15),
* AC(10)/ID/DN(40)9K(7)9
* P1(10)9P2(10)9EL(20)9
* NO(20)9MX(80)/BIBL/
* F(35099) 9NM(350,4)/BEZ/
* BZ(100)9T(100)//
* A(1009101)9STB(101)9
* X(350)9EX(350)9NJI(21)
* /PR/S(350)9H(350)9CP(350) DIMENSION PP(6)9W(100)
Q*******************************
С ИТЕРАЦИОННЫЙ ПРОЦЕСС
С ОПРЕДЕЛЕНИЯ РАВНОВЕСНОГО
С СОСТАВА ДЛЯ ЗАДАННОЙ ПАРЫ
С ТЕРМОДИНАМИЧЕСКИХ ПАРАМЕТРОВ
с*******************************
DATA PP/2HP=92HT=,2HV=9
* 2HS=92HI=92HU=/ В0=1.98719 А1=-1000. А2=1000.
NC(15)=0 IT=O XL=-50. ITEB=O
IP1=NC(10)/1000
IP2=NC(10)-IP1*1000
IP3=NC(11)/1000
IP4=NC(11)-IP3*1000
IP5=NC(12)/1000
IP6=NC(12)-IP5*1000
L1=K(5)+NC(4)
L2=NC(3)
L3=K(5)
IM1=L1+NC(7)+7
IF (NC(5).EQ.0) IM1=DI1-1
IM2=IM1+1
1 IF(K(I).EQ.2) BZ(2)=AC(5) IF(K(2).EQ.2) BZ(2)=AC(6) CALL PBOPS (BZ(2))
2 EPS=O.
DO 3 I=I9M DO 3 J=I9MM
3 A(I9J)=O. DO 6 J=I,2 L4=K(J) A(J9L4)=1.0
IF (L4-3) 49594
4 A(J9MM)=AC(J+4) GO TO 6
5 A(J92)=-1.0/BZ(2) A(J9MM)=DLOG(B0*BZ(2)/
* AC(J+4))-1.0
6 CONTINUE A(3,1)=-DEXP(-BZ(3)) DO 7 1=19L3
7 A(I+3,MM)=EL(I) DO 8 1=4,6 IK=I+L1
8 A(UC,1)=-1.
IF (NC(5)) 10,10,9
9 A(IMI ,IMI )=1 ./Al A(IMI ,MM) = I .0
213
10 IF (NC(6)) 12,12,11
11 A(IM2,IM2)=1 ./A1 A(IM2,MM)=KO
12 UC=LI +7
13 XI=-1D+6
DO 43 1=1 ,L2 JI=NM(I,3)
CALL STRING(I,EX( I),RZ( 2)) IF ( Jl .EQ. 1.OR.Jl .GT.4)
* GO TO 16 JA=3
14 A(JA,2)=STR(1)+A(JA,2)
A(JA,3)=STR(2)+2.0*EX(I)+
* A(JA,3)
A( JA,MM)=A( JA,MM)-STR( 3 )+
* 2.0*RZ(3)*EX(I) DO 15 J=I,L1
15 A(JA,J+6)=STR(J+3)+
* A(JA,J+6)
16 DO 22 J2=1,L3
IF(NJKJ2).BQ.O) GO TO 22
JA=J2+3
NJ=NJI(J2)
GO TO (17,21,21,21,18,20),
* Jl
17 EK=DEXP(RZ(IK)) A(JA, IK) =NJ*SC A(JA9MM)=(RZ(JJC)-I. )*NJ*
* EK+A(JA,MM) GO TO 22
18 EK=DEXP(RZ(IMI)) Y=NJ*EK
A(JA,IMI )=Y*EX(I)+
* A(JA1IMI) A(JA,MM)=Y*EX(I)*RZ(IM1 ) +
* A(JA1MM)
19 A(JA,2)=STR(1)»Y+A(JA 9 2) A(JA9 3)=STR(2)*Y+A(JA9 3) A(JA,MM)=A(JA,MM)-STR(3)*Y DO 80 J= 1 ,Ll
80 A ( JA, J+6 ) =STR( J+3) * Y+
* A(JA,J+6) GO TO 22
20 EK=DEXP(RZ(IM2)) Y=NJ*EK
A(JA,IM2)=Y*EX(I)+
* A(JA,IM2)
A( JA,MM) =Y*EX( I) *RZ( IM2 )+
* A(JA1MM)
GO TOI9
21 Y=NJ
GO TO 19
22 CONTINUE
?***«
JA=LI+4
OO TO (23,24,24,24,26,27),
* J1
23 A(JA, 2)=CP(I)*EH/RZ( 2)+
* A(JA, 2) A(JA,IK)=S(I)*EK A(JA,MM)=( (RZ(JJC)-I . )* -
* S(I)+CP(I))*EK+A(JA,MM) GO TO 28
24 A(JA,2)=EX(I)*CP(I)/RZ(2)
* +A(JA,2) A(JA,3)=A(JA,3)-R0*EX(I) A(JA,MM)=(CP(I)-R0*(RZ(3)
* +1.))*EX(I)+A( JA1MM) Y=S(I)-R0*(1.+X(I)+RZ(3)) GO TO 25
26 Y=EK*EX(I) A(JA,2)=CP(I)*Y/RZ(2)+
* A(JA,2)
A(JA9IMI)=(S(I)-R0*X(I))*
* Y+A(JA,IM1)
A(JA,MM)=(S(I)*RZ(IM1)+
* (1.0-RZ(JJlI )*X(I))*R0+
* CP(I))*Y+A(JA,MM) Y=((1.0-X(I))*R0+S(I))*EK GO TO 25
27 Y=EK*EX(I) A(JA,2)=CP(I)*Y/RZ(2)+
* A(JA, 2)
A(JA,IM2)=(S(I)-R0*X(I))*
214
* Y+A(JA9IM2) A(JA9MM)=(S(I)*RZ(IM2) +
* (1.0-RZ(IM2)*X(I))*R0+
* CP(I ))*Y+A (JA9MM) Y=((1.0-X(I))*R0+S(I))*EK
25 A(JA,2)=STR(1)*Y+A( J A 92) A(JA93)=STR(2)*Y+A(JA93) A(JA9MM)=A(JA9MM)-
* STR(3)*Y DO 81 J=19L1
81 A(JA,J+6)=STR(J+3)*Y+
* A(JA9J+6)
o***«
28 U=H(I)
DO 35 J2=1,2 JA=LI+4+J2
GO TO (29,30,30,30,32,
* 33), Л
29 A(JA,IK)=U*EK A(JA92)=CP(I)*EK+A(JA,2) A(JA9MM)=((RZ(IK)-I.)*U+
* CP(I)*RZ(2))*EK+A(JA9MM) GO TO 34
30 A(JA92)=CP(I)*EX(I)+
* A(JA,2) A(JA9MM)=CP(I)*EX(I)*
* RZ(2)+A(JA,MM) Y=U
GO TO 31
32 A(JA,2)=CP(I)*EK*EX(I)+
* A(JA,2) Y=U*EK
A(JA1IMI)=Y*EX(I)+
* A(JA1IMI)
A(JA9MM)=(U*RZ(IM1)+
* CP(I)*RZ(2))*EK*EX(I)+
* A(JA1MM) GO TO 31
33 A(JA,^)=CP(I)*EK*EX(I) +
* A(JA,2) Y=U *EK
A(JA9IM2)=Y*EX(D +
* A(JA1IM2)
A(JA9MM) = (U*RZ(IM2) +
* CP(I)*RZ(2))*EK*EX(I)+
* A(JA9MM)
31 A(JA92)=STR(1)*Y+A(JA,2) A(JA93)=STR(2)*Y+A(JA93) A(JA9MM)=A(JA9MM)-
* STR(3)*Y DO 82 J=I,Ll
82 A(JA9 J+6) =STR( J+3 )*Y+
* A(JA9J+6)
34 U=H(I)-R0*RZ(2)
Предыдущая << 1 .. 82 83 84 85 86 87 < 88 > 89 90 91 92 93 94 .. 103 >> Следующая
Реклама
 
 
Авторские права © 2010 PiroChem. Все права защищены.