Матэматычная мадэль сінхроннага генератара (Лабараторная работа № 10), страница 9

*Sn,Un,Ufn,Ifxx,Ifn,Jrt,Rs,Rf,Xd11,Xd1,Xd,X2,Td0,Td11,              81

*Snagr,CosFI,GAM                                                    82

CLOSE(UNIT=1)                                                      90

WRITE(5,*)'SINGEN: ИСХОДНЫЕ ДАННЫЕ СЧИТАНЫ С ДИСКА'                100

C        ОБРАЗОВАНИЕ ИМЕНИ ФАЙЛА РЕЗУЛЬТАТОВ РАСЧЕТА В ВИДЕ AAAA.REZ, ГДЕ

C        AAAA - ЭТО ИМЯ ФАЙЛА ИСХОДНЫХ ДАННЫХ, А ТАКЖЕ ИМЕНИ ФАЙЛА AAAA.GRF,

C        В КОТОРЫЙ ПОМЕЩАЕТСЯ ВЫХОДНАЯ ИНФОРМАЦИЯ В ВИДЕ, УДОБНОМ ДЛЯ ВЫВОДА

C        ОСЦИЛЛОГРАММ НА ЭКРАН ДИСПЛЕЯ В ГРАФИЧЕСКОМ РЕЖИМЕ

Z2(1:4)=Z1(1:4)                                                    110

Z3(1:4)=Z1(1:4)                                                    111

Z2(5:8)='.REZ'                                                     112

Z3(5:8)='.GRA'

C        ОТКРЫТИЕ ФАЙЛОВ РЕЗУЛЬТАТОВ НА ДИСКЕ

OPEN(UNIT=3,FILE=Z2,STATUS='NEW')                                  120

OPEN(UNIT=4,FILE=Z3,STATUS='NEW')                                  120

C        ЗАПИСЬ ИСХОДНЫХ ДАННЫХ В ВЫХОДНОЙ ФАЙЛ С ИМЕНЕМ AAAA.REZ

WRITE(3,*)'SINGEN - РАСЧЕТ ТРЕХФАЗНОГО К.З.В СИНХРОННОМ ГЕН-РЕ'    130

WRITE(3,*)' ИСХОДНЫЕ ДАННЫЕ ДЛЯ СИСТЕМЫ ДИФ.УРАВНЕНИЙ'             131

WRITE(3,*)NUR,INTWR,T0,TKON,H                                      132

WRITE(3,*)' ИСХОДНЫЕ ДАННЫЕ ДЛЯ ГЕНЕРАТОРА'                        133

WRITE(3,*)Sn,Un,Ufn,Ifxx,Ifn,Jrt,Rs,Rf,Xd11,Xd1,Xd,                134

*X2,Td0,Td11,Snagr,CosFI,GAM                                        135

C        РАСЧЕТ ПАРАМЕТРОВ ГЕНЕРАТОРА, ОТСУТСТВУЮЩИХ В КАТАЛОГЕ

OMEGA=314.15928                                                    140

Xq=0.9*Xd                                                          141

Xasig=Xd11-0.025                                                   142

Xad=Xd-Xasig                                                       143

Xaq=Xq-Xasig                                                       144

Xfsig=Xad*(Xd1-Xasig)/(Xd-Xd1)                                     145

Xf=Xfsig+Xad                                                       146

Xusigd=1./(1./(Xd11-Xasig)-1./Xfsig-1./Xad)                        147

Xud=Xusigd+Xad                                                     148

Xq11=X2*X2/Xd11                                                    149

Xusigq=Xaq*(Xq11-Xasig)/(Xq-Xq11)                                  150

Xuq=Xaq+Xusigq                                                     151

Ibaz=SQRT(2./3.)*Sn/Un                                             152

Ubaz=SQRT(2./3.)*Un                                                153

Sbaz=1.5*Ubaz*Ibaz                                                 154

Zbaz=Ubaz/Ibaz                                                     155

Rso=Rs/Zbaz                                                        156

Rfo=Xf/(OMEGA*Td0)                                                 157

Xud11=Xusigd+1./(1./Xasig+1./Xad+1./Xfsig)                         158

Rudo=Xud11/(OMEGA*Td11)                                            159

Xuq11=Xusigq+1./(1./Xasig+1./Xaq)                                  160

Ruqo=Xuq11/(OMEGA*Td11)                                            161

Ifbaz=Xad*Ifxx/1.05                                                162

Hj=Jrt*OMEGA*OMEGA/Sbaz                                            163

C         PACЧЕТ ПАРАМЕТРОВ РЕЖИМА ГЕНЕРАТОРА ДЛЯ МОМЕНТА ВРЕМЕНИ Т=0.

SinFI=SQRT(1.-CosFI*CosFI)                                         170

FI=ACOS(CosFI)                                                     171

Inagr=SQRT(2.)*Snagr/(SQRT(3.)*Un*Ibaz)                            172

TETA=ATAN(Inagr*(Xq*CosFI-Rso*SinFI)/(1.+Inagr*SinFI*(Rso+Xq)))    173

Id=-Inagr*SIN(TETA+FI)                                             174

Iq=Inagr*COS(TETA+FI)                                              175

Ud=-SIN(TETA)                                                      176

Uq=COS(TETA)                                                       177

Eq=Uq-Id*Xd+Rso*Inagr*COS(TETA+FI)                                 178

If=Eq/Xad                                                          179

Uf=If*Rfo                                                          180

PSId=Xd*Id+Xad*If                                                  181

PSIq=Xq*Iq                                                         182