Подякувати Студентському архіву довільною сумою
Admin
26.02.2023 12:38
Дякуємо, що користуєтесь нашим архівом!
Математичне моделювання електромеханічних пристроїв
Інформація про навчальний заклад
ВУЗ:
Національний університет Львівська політехніка
Інститут:
Не вказано
Факультет:
Не вказано
Кафедра:
Не вказано
Інформація про роботу
Рік:
2004
Тип роботи:
Розрахункова робота
Предмет:
Інші
Частина тексту файла (без зображень, графіків і формул):
Міністерство освіти та науки України
Національний університет
“Львівська політехніка”
кафедра ЕМА
Розрахункова робота
“Математичне моделювання електромеханічних пристроїв”
Математична модель подвоювача частоти
1. Завдання. Скласти математичну модель електромагнітних перехідних процесів для подвоювала частоти, схема якого зображена на рис. 1
рис. 1. Подавач частоти
Характеристика намагнічування сталі задана:
Вебер-амперна характеристика дроселя задана:
Напруги живлення є задані функціями часу:
Вхідні дані:
Сумарні опори обмоток підмагнічування та первинних обмоток:
Сумарні індуктивності розсіяння обмоток підмагнічування та первинних обмоток:
Кількість витків обмоток підмагнічування, первинних та вторинних обмоток:
Довжина силової лінії елементів магнітопровода:
Площа поперечного перетину магнітопроводів:
Розрахункова робота
з предмету "Математичне моделювання електромеханічних перетворювачів |
Виконав: студент групи ЕМА-41 Осадця О.П.
!-------------------- Головний модуль програми ------------------------|
! Підключення модуля даних
Use EMA
! Описова частина
Dimension X_POPR(5),Y(4)
! Читання вхідних даних, початкової умови та вивід їх у файл вихідної
! інформації
CALL DANI
! Час початку інтегрування
T=T0
! Запис у файл вихідної інформації даних сформованих на основі
! початкової умови
CALL REZ(T)
! Формування кореня 0-го кроку інтегрування
CALL POPER(X_POPR)
! Початкове значення лічильника кроку введення результатів
LIC=0
! Мітка початку кроку інтегрування
100 CONTINUE
! Визначення часу на кроці інтегрування
T=T+h
Write(*,*)' Час на кроці інтегрування T= ',T
! Перевірка умови: якщо час на кроці інтегрування більший за час
! припинення інтегрування, то слід припинити роботу програми :
IF(T.GT.Tp) GOTO 110
! Визначення значення напруг на кроці інтегрування
CALL VOLT(T)
! Занулення лічильника ітерацій
NI=0
! Мітка початку ітерації методу Ньютона
200 CONTINUE
! Визначення номеру ітерації
NI=NI+1
! Перевірка умови: якщо кількість ітерацій перевищує 25 (це означає,
! що ітераційна процедура розбігається), слід зробити аварійну
! зупинку програми
IF(NI.GT.25) GOTO 120
! Обчислення нев'язки
CALL NVZ(X_POPR,Y)
! Обчислення похідної нев'язки
CALL MATR_A
! Визначення поправки первинних невідомих
CALL GAUSS(A,Y,4,4,DELTA_XP)
! Визначення первинних невідомих
CALL PR_NV
! Визначення вторинних невідомих
CALL FORMUL
! Прийняття рішення про припинення ітераційного процесу
call NSTOP(XP,DELTA_XP,4,NI,KZ)
! Перевірка умови: якщо KZ=0, то виконуємо наступну ітерацію
IF (KZ.EQ.0) GOTO 200
write(*,*)'ni',ni
WRITE(*,*)'Кількість ітерацій на кроці інтегрування NI= ',NI
! Збільшення значення лічильника кроку виведення результатів
LIC=LIC+1
! Перевірка умови: якщо значення лічильника дорівнює кроку виведення
! результатів, то вони записуються у файл
IF (LIC.EQ.krok) THEN
! Запис результатів обчислень
CALL REZ(T)
! Занулення значення лічильника кроку виведення результатів
LIC=0
END IF
! Формування кореня кроку інтегрування
CALL POPER(X_POPR)
! Повернення на початок наступного кроку інтегрування
GOTO 100
! Мітка аварійної зупинки програми
120 Write(*,*)' Аварійна зупинка програми NI>25 '
! Закриваємо файл результатів
CLOSE(20)
STOP
! Мітка припинення інтегрування
110 Write(*,*)' Розрахунок закінчено ! '
! Закриваємо файл результатів
CLOSE(20)
END
!_____________________________________________________________________|
! Підпрограма для читання вхідних даних, початкової умови та |
! виводу їх у файл вихідної інформації |
!_____________________________________________________________________|
SUBROUTINE DANI
!_____________________________________________________________________|
! Вхідні дані: немає |
! Вихідні дані: |
! XP0 - вектор початкових значень первинних невідомих |
! XV0 - вектор початкових значень вторинних невідомих |
!_____________________________________________________________________|
! Підключення модуля даних
Use EMA
! Описова частина
Character Sim*60,Prop*78
! Відкриваємо файли вхідної та вихідної інформації
open(10,file='dani.dat',status='old')
open(20,file='rez.dat',status='unknown')
open(30,file='rez2.dat',status='unknown')
! Читаємо та виводимо заголовок
Do i=1,4
Read(10,'(A78)')Prop
Write(20,*)Prop
End do
! Читаємо та виводимо вхідні дані
20 Format(A60,E9.4)
Read(10,20)Sim,R0
Write(20,*)Sim,R0
Read(10,20)Sim,R1
Write(20,*)Sim,R1
Read(10,20)Sim,Ls0
Write(20,*)Sim,Ls0
Read(10,20)Sim,Ls1
Write(20,*)Sim,Ls1
Read(10,20)Sim,w0
Write(20,*)Sim,w0
Read(10,20)Sim,w1
Write(20,*)Sim,w1
Read(10,20)Sim,w2
Write(20,*)Sim,w2
Read(10,20)Sim,L
Write(20,*)Sim,L
Read(10,20)Sim,S
Write(20,*)Sim,S
Read(10,20)Sim,Ah
Write(20,*)Sim,Ah
Read(10,20)Sim,Bh
Write(20,*)Sim,Bh
Read(10,20)Sim,Kpsi
Write(20,*)Sim,Kpsi
Read(10,20)Sim,Mpsi
Write(20,*)Sim,Mpsi
Read(10,20)Sim,Ppsi
Write(20,*)Sim,Ppsi
Read(10,20)Sim,U0
Write(20,*)Sim,U0
Read(10,20)Sim,U1
Write(20,*)Sim,U1
Read(10,20)Sim,Tu
Write(20,*)Sim,Tu
Read(10,20)Sim,f
Write(20,*)Sim,f
Read(10,20)Sim,alfa
Write(20,*)Sim,alfa
Read(10,20)Sim,h
Write(20,*)Sim,h
Read(10,20)Sim,To
Write(20,*)Sim,To
Read(10,20)Sim,Tp
Write(20,*)Sim,Tp
30 Format(A60,I3)
Read(10,30)Sim,Krok
Write(20,30)Sim,Krok
! Читаємо та виводимо другий заголовок
Read(10,'(A80)')Prop
Write(20,*)Prop
! Читаємо в масив XP0 та виводимо вектор первинних невідомих
Do i=1,4
Read(10,20)Sim,XP(i)
Write(20,*)Sim,XP(i)
end do
! Читаємо в масив XV0 та виводимо вектор вторинних невідомих
Do i=1,7
Read(10,20)Sim,XV(i)
Write(20,*)Sim,XV(i)
end do
! Закриваємо файл вхідної інформації
Close(10)
! Записуємо у файл вихідної інформації "шапку".
write(20,10)
10 format(/96('-')/'| T,c | I1,A | I0,A | B1,Тл | B2,Тл | Psi11 | Psi12 | Psi01 | Psi02 | PsiD | H1 | H2 |'/96('-'))
Return
END
!____________________________________________________________________|
! Підпрограма для визначення вторинних невідомих |
!____________________________________________________________________|
SUBROUTINE FORMUL
!____________________________________________________________________|
! Вхідні величини : |
! XP - вектор первинних невідомих, обчислений на l-й ітерації. |
! X_POPR - Корінь відомий з попереднього кроку інтегрування. |
! Вихідні величини : |
! XV - вектор вторинних невідомих, обчислений на l-й ітерації. |
!____________________________________________________________________|
! Підключення модуля даних
Use EMA
! Визначення вторинних невідомих
XV(1)=Ls1*XP(1)/2.+w1*XP(3)*S
XV(2)=Ls1*XP(1)/2.-w1*XP(4)*S
XV(3)=Ls0*XP(2)/2.+w0*XP(3)*S
XV(4)=Ls0*XP(2)/2.+w0*XP(4)*S
XV(5)=Kpsi*XP(2)+Mpsi*(1.-exp(-Ppsi*XP(2)))
XV(6)=ah*XP(3)+bh*XP(3)**11
XV(7)=ah*XP(4)+bh*XP(4)**11
Return
End
!-------------------------------------------------
! Підпрограма розв'язування лінійної САР !
! методом Гаусса !
!-------------------------------------------------
SUBROUTINE GAUSS(A,B,N,NM,X)
!-------------------------------------------------
! Вхідні величини: !
! A - матриця коефіцієнтів; !
! B -стовпець вільних членів; !
! N - порядок системи; !
! NM - максимальний порядок системи. !
! Вихідні величини: !
! X - стовпець розв'язків. !
! Матриця коефіцієнтів зберігає свої значення. !
!-------------------------------------------------
DIMENSION A(NM,NM),B(NM),X(NM)
DIMENSION AA(500,500),JP(500),KP(500),AP(500,500), &
BP(500,500),AX(500),AI(500,500),BI(500)
!--Переприсвоєння матриці А-------
DO 110 IR=1,N
DO 110 JR=1,N
AA(IR,JR)=A(IR,JR)
110 CONTINUE
!-----------------------------------------
IS=0
50 AM=0.
IS=IS+1
DO 10 J=1,N
DO 10 K=1,N
IF(AM.GT.ABS(AA(J,K))) GOTO 10
AM=ABS(AA(J,K))
JM=J
KM=K
JP(IS)=JM
KP(IS)=KM
10 CONTINUE
DO 20 K=1,N
20 AP(IS,K)=AA(JM,K)/AA(JM,KM)
BP(JM,IS)=B(JM)/AA(JM,KM)
DO 30 J=1,N
DO 40 K=1,N
40 AI(J,K)=-AA(J,KM)*AP(IS,K)+AA(J,K)
BI(J)=-AA(J,KM)*BP(JM,IS)+B(J)
30 CONTINUE
DO 90 J=1,N
DO 100 K=1,N
100 AA(J,K)=AI(J,K)
B(J)=BI(J)
90 CONTINUE
IF(IS.LT.N) GOTO 50
DO 60 K=1,N
60 X(K)=0.
DO 80 IS=N,1,-1
L=KP(IS)
M=JP(IS)
AX(L)=0.
AP(IS,L)=0.
DO 70 K=1,N
AX(L)=AX(L)+AP(IS,K)*X(K)
70 CONTINUE
X(L)=BP(M,IS)-AX(L)
80 CONTINUE
RETURN
END
!____________________________________________________________________|
! Підпрограма для обчислення похідної вектора нев'язки |
!____________________________________________________________________|
SUBROUTINE MATR_A
!____________________________________________________________________|
! Вхідні величини : |
! XP - вектор первинних невідомих. |
! Вихідні величини : |
! A - похідна вектора нев'язки. |
!____________________________________________________________________|
! Підключення модуля даних
Use EMA
! Обчислення диференційних параметрів
ANU1=ah+11.*bh*XP(3)**10
ANU2=ah+11.*bh*XP(4)**10
ALD=Kpsi-Mpsi*exp(-XP(2))
! Обчислення елементів похідної вектора нев'язки
A(1,1)=-Ls1/h-R1
A(1,2)=0.
A(1,3)=-w1*s/h
A(1,4)=w1*S/h
A(2,1)=0.
A(2,2)=-Ls0/h-ALD/h-R0
A(2,3)=-w0*S/h
A(2,4)=-w0*s/h
A(3,1)=-w1
A(3,2)=-w0
A(3,3)=L*ANU1
A(3,4)=0.
A(4,1)=-w1
A(4,2)=w0
A(4,3)=0.
A(4,4)=-L*ANU2
Return
End
!____________________________________________________________________|
! Підпрограма, яка здійснює припинення ітераційного процесу при |
! розв'язуванні нелінійної системи рівнянь методом Ньютона |
!____________________________________________________________________|
SUBROUTINE N_STOP(KZ)
!____________________________________________________________________|
! Вхідні величини : |
! XP - вектор первинних невідомих, обчислений на l-й ітерації. |
! DELTA_XP - вектор поправок первинних невідомих, |
! обчислений на l-й ітерації. |
! Вихідні величини : |
! KZ - керуюча змінна (KZ=0, то необхідно виконувати наступну |
! ітерацію; KZ=1 - ітераційний процес закінчений). |
!____________________________________________________________________|
! Підключення модуля даних
Use EMA
dimension xpet(4)
! Перевіряємо умову: якщо сота частина однієї з поправок змінює |
! значення вектора первинних невідомих, то ітераційний процес |
! продовжується
do i=1,4
xpet(i)=XP(i)+0.01*DELTA_XP(i)
end do
Do i=1,4
IF (XP(i).NE.XPet(i)) Then
KZ=0
exit
else
KZ=1
end if
end do
return
end
!***********************************************
!* Підпрограма перевірки умови припинення *
!* ітераційної процедури *
!* при розв'язуванні систем нелінійних рівнянь *
!* методом Ньютона *
!***********************************************
SUBROUTINE NSTOP_NEW(XPPOP,XP,N,NI,KZ)
!------------------------------------------------------------
! Вхідні величини:
! XP(N)-вектор первинних невідомих на (l-1)-й ітерації;
! XPPOP(N)-вектор первинних невідомих на l-й ітерації;
! N-кількість первинних невідомих;
! NI-номер ітерації;
! KZ-керуюча змінна:
! якщо KZ=0-ітераційну процедуру слід продовжити;
! якщо KZ=1-ітераційну процедуру слід припинити.
!------------------------------------------------------------
IMPLICIT REAL*8 (A-H,O-T,X,U,Z,Y)
CHARACTER WPPOP*35,WP*35
DIMENSION XP(N),XPPOP(N)
SELECT CASE(NI) ! номер ітерації
CASE(1:3) ! 1-а - 3-я ітерації
NN=12 ! перевірка, чи поправка уточнює чотирнадцяту значущу цифру
CASE(4:6) ! 4-а - 6-а ітерації
NN=11 ! перевірка, чи поправка уточнює тринадцяту значущу цифру
CASE(7:9) ! 7-а - 9-а ітерації
NN=10 ! перевірка, чи поправка уточнює дванадцяту значущу цифру
CASE(10:11) ! 10-а - 11-а ітерації
NN=9 ! перевірка, чи поправка уточнює одинандцяту значущу цифру
CASE(12:13) ! 12-а - 13-а ітерації
NN=8 ! перевірка, чи поправка уточнює десяту значущу цифру
CASE(14:) ! 14-а - 15-а ітерації
NN=7 ! перевірка, чи поправка уточнює дeв'яту значущу цифру
END SELECT
DO I=1,N
IF(DABS(XP(I)).LT.1.D-07) THEN
KZ=1
CYCLE
END IF
WRITE(WPPOP,'(E20.14)')XPPOP(I)
WRITE(WP,'(E20.14)')XP(I)
IF(WP(1:NN+2)//WP(17:20).NE.WPPOP(1:NN+2)//WPPOP(17:20)) THEN
KZ=0
RETURN
ELSE
KZ=1
END IF
END DO
RETURN
END
!-----------------------------------------------------
! Підпрограма перевірки умови припинення !
! ітераційної процедури при розв'язуванні !
! нелінійної системи рівнянь методом Ньютона !
!-----------------------------------------------------
SUBROUTINE NSTOP(X,DX,N,NI,MG)
!-----------------------------------------------------
! Вхідні величини: !
! X(N) - вектор первинних невідомих; !
! DX(N) - вектор поправок первинних невідомих; !
! N - порядок системи рівнянь; !
! NI - номер ітерації; !
! Вихідні величини: !
! MG - керуюча змінна (якщо MG=0, то необхідно вико-!
! нувати наступну ітерацію методу Ньютона, якщо !
! MG=1, то слід припинити ітераційну процедуру. !
!-----------------------------------------------------
! IMPLICIT REAL*8 (A-H,O-T,X,U,Z,Y)
DIMENSION X(N),DX(N)
FAN=0.1
IF(NI.GT.3.AND.NI.LE.6) FAN=0.01
IF(NI.GT.6.AND.NI.LE.9) FAN=0.001
IF(NI.EQ.10.OR.NI.EQ.11) FAN=0.0001
IF(NI.EQ.12.OR.NI.EQ.13) FAN=0.00001
IF(NI.GE.14) FAN=0.000001
DO 10 I=1,N
XP=FAN*DX(I)
A=ABS(X(I))
AA=ABS(XP)
IF(DX(I).EQ.0.) GOTO 10
IF(X(I).EQ.0.) GOTO 5
DO 20 J=0,1000
TEN=10.**J
B=A*TEN
BB=AA*TEN
IF(A.GE.1.) THEN
B=A/TEN
BB=AA/TEN
END IF
IF(B.LT.1.0.AND.B.GE.0.1) GOTO 30
20 CONTINUE
GOTO 30
5 BB=AA
30 MM=BB*10000000
IF(MM.NE.0) GOTO 40
10 CONTINUE
MG=1
GOTO 50
40 MG=0
50 CONTINUE
RETURN
END
!***********************************************
!* Підпрограма перевірки умови припинення *
!* ітераційної процедури *
!* при розв'язуванні систем нелінійних рівнянь *
!* методом Ньютона *
!***********************************************
SUBROUTINE NSTOP_NEW(XPPOP,XP,N,NI,KZ)
!------------------------------------------------------------
! Вхідні величини:
! XP(N)-вектор первинних невідомих на (l-1)-й ітерації;
! XPPOP(N)-вектор первинних невідомих на l-й ітерації;
! N-кількість первинних невідомих;
! NI-номер ітерації;
! KZ-керуюча змінна:
! якщо KZ=0-ітераційну процедуру слід продовжити;
! якщо KZ=1-ітераційну процедуру слід припинити.
!------------------------------------------------------------
IMPLICIT REAL*8 (A-H,O-T,X,U,Z,Y)
CHARACTER WPPOP*35,WP*35
DIMENSION XP(N),XPPOP(N)
SELECT CASE(NI) ! номер ітерації
CASE(1:3) ! 1-а - 3-я ітерації
NN=12 ! перевірка, чи поправка уточнює чотирнадцяту значущу цифру
CASE(4:6) ! 4-а - 6-а ітерації
NN=11 ! перевірка, чи поправка уточнює тринадцяту значущу цифру
CASE(7:9) ! 7-а - 9-а ітерації
NN=10 ! перевірка, чи поправка уточнює дванадцяту значущу цифру
CASE(10:11) ! 10-а - 11-а ітерації
NN=9 ! перевірка, чи поправка уточнює одинандцяту значущу цифру
CASE(12:13) ! 12-а - 13-а ітерації
NN=8 ! перевірка, чи поправка уточнює десяту значущу цифру
CASE(14:) ! 14-а - 15-а ітерації
NN=7 ! перевірка, чи поправка уточнює дeв'яту значущу цифру
END SELECT
DO I=1,N
IF(DABS(XP(I)).LT.1.D-07) THEN
KZ=1
CYCLE
END IF
WRITE(WPPOP,'(E20.14)')XPPOP(I)
WRITE(WP,'(E20.14)')XP(I)
IF(WP(1:NN+2)//WP(17:20).NE.WPPOP(1:NN+2)//WPPOP(17:20)) THEN
KZ=0
RETURN
ELSE
KZ=1
END IF
END DO
RETURN
END
!-----------------------------------------------------
! Підпрограма перевірки умови припинення !
! ітераційної процедури при розв'язуванні !
! нелінійної системи рівнянь методом Ньютона !
!-----------------------------------------------------
SUBROUTINE NSTOP(X,DX,N,NI,MG)
!-----------------------------------------------------
! Вхідні величини: !
! X(N) - вектор первинних невідомих; !
! DX(N) - вектор поправок первинних невідомих; !
! N - порядок системи рівнянь; !
! NI - номер ітерації; !
! Вихідні величини: !
! MG - керуюча змінна (якщо MG=0, то необхідно вико-!
! нувати наступну ітерацію методу Ньютона, якщо !
! MG=1, то слід припинити ітераційну процедуру. !
!-----------------------------------------------------
IMPLICIT REAL*8 (A-H,O-T,X,U,Z,Y)
DIMENSION X(N),DX(N)
FAN=0.1
IF(NI.GT.3.AND.NI.LE.6) FAN=0.01
IF(NI.GT.6.AND.NI.LE.9) FAN=0.001
IF(NI.EQ.10.OR.NI.EQ.11) FAN=0.0001
IF(NI.EQ.12.OR.NI.EQ.13) FAN=0.00001
IF(NI.GE.14) FAN=0.000001
DO 10 I=1,N
XP=FAN*DX(I)
A=DABS(X(I))
AA=DABS(XP)
IF(DX(I).EQ.0.) GOTO 10
IF(X(I).EQ.0.) GOTO 5
DO 20 J=0,1000
TEN=10.**J
B=A*TEN
BB=AA*TEN
IF(A.GE.1.) THEN
B=A/TEN
BB=AA/TEN
END IF
IF(B.LT.1.0.AND.B.GE.0.1) GOTO 30
20 CONTINUE
GOTO 30
5 BB=AA
30 MM=BB*10000000
IF(MM.NE.0) GOTO 40
10 CONTINUE
MG=1
GOTO 50
40 MG=0
50 CONTINUE
RETURN
END
!____________________________________________________________________|
! Підпрограма для обчислення нев'язки |
!____________________________________________________________________|
SUBROUTINE NVZ(X_POPR,Y)
!____________________________________________________________________|
! Вхідні величини : |
! X_POPR- корінь, відомий з попереднього кроку інтегрування. |
! U0_T,U1_T - значення напруг на кроці інтегрування. |
! XP,XV - вектори первинних та вторинних невідомих |
! Вихідні величини : |
! Y - вектор нев'язки. |
!____________________________________________________________________|
! Підключення модуля даних
Use ema
! Описова частина
Dimension X_POPR(5),Y(4)
! Визначення елементів вектора нев'язки
Y(1)=-(XV(1)-X_POPR(1))/h-R1*XP(1)-(XV(2)-X_POPR(2))/h+U1_T
Y(2)=-(XV(3)-X_POPR(3))/h-(XV(5)-X_POPR(5))/h-R0*XP(2)-(XV(4)-X_POPR(4))/h+U0_T
Y(3)=l*XV(6)-w1*XP(1)+w0*XP(2)
Y(4)=-l*XV(7)-w1*XP(1)+w0*XP(2)
! Зміна знаку нев'язки
Do i=1,4
Y(i)=-Y(i)
End do
Return
End
!____________________________________________________________________|
! Підпрограма для збереження кореня біжучого кроку інтегрування |
!____________________________________________________________________|
SUBROUTINE POPER(X_POPR)
!____________________________________________________________________|
! Вхідні величини : |
! XP,XV - вектори первинних та вторинних невідомих |
! Вихідні величини : |
! X_POPR - корінь. |
!____________________________________________________________________|
! Підключення модуля даних
Use EMA
! Описова частина
Dimension X_POPR(5)
! Переприсвоювання значень
Do i=1,5
X_POPR(i)=XV(i)
End do
Return
End
!____________________________________________________________________|
! Підпрограма для визначення первинних невідомих |
!____________________________________________________________________|
SUBROUTINE PR_NV
!____________________________________________________________________|
! Вхідні величини : |
! XP - вектор первинних невідомих, обчислений на (l-1)-й ітерації. |
! DELTA_XP - вектор поправок первинних невідомих, обчислений на l-й |
! ітерації. |
! Вихідні величини : |
! XP - вектор первинних невідомих, обчислений на l-й ітерації. |
!____________________________________________________________________|
! Підключення модуля даних
Use EMA
! Визначення первинних невідомих
DO i=1,4
XP(i)=XP(i)+DELTA_XP(i)
End do
Return
End
Файл вхідної інформації до розрахункової роботи
(ст. гр. ЕМА-41 Осадці О.П.)
1. Опір 1-ї котушки R0 (Ом) .5600E+02
2. Опір 2-ї котушки R1 (Ом) .9000E+03
3. Індуктивність розсіяння Ls0 (Гн) .2000E-03
4. Індуктивність розсіяння Ls1 (Гн) .2000E-03
5. Кількість витків котушки w0 .1000E+03
6. Кількість витків котушки w1 .1500E+04
7. Кількість витків котушки w2 .7500E+03
8. Довжина силової лінії магнітопровода L (м) .3000E+00
9. Площа поперечного перерізу магнітопроводу S (м^2) .3000E-03
10.Коефіцієнт х-ки намагнічування сталі Ah (А/м*Тл) .3050E+02
11.Коефіцієнт х-ки намагнічування сталі Bh (А/м*Тл) .8000E+04
12.Коефіцієнт Веб-амперної х-ки дроселя Kpsi(Вб/А) .8000E+00
13.Коефіцієнт Веб-амперної х-ки дроселя Mpsi(Вб) .1000E+00
14.Коефіцієнт Веб-амперної х-ки дроселя Ppsi(1/А) .1000E+01
15.Амплітуда напруги живлення U0 (В) .3000E+02
16.Амплітуда напруги живлення U1 (В) .1500E+03
17.Стала часу Tu (c) .2000E+03
18.Частота мережі живлення f (Гц) .5500E+02
19.Кут зсуву фаз alfa .5233E+00
20.Величина кроку інтегрування h (с) .1000E-03
21.Час початку інтегрування To (с) .0000E+00
22.Час припинення інтегрування Tp (с) .1000E+00
23.Крок виведення результатів (1,2,...) Krok 1
Нульове наближення невідомих
Первинні невідомі: Cтрум i1 (А) .0000E+00
Cтрум i0 (А) .0000E+00
Індукція B1 (Тл) .0000E+00
Індукція B2 (Тл) .0000E+00
Вторинні невідомі: Потокощеплення Psi11 .0000E+00
Потокощеплення Psi12 .0000E+00
Потокощеплення Psi01 .0000E+00
Потокощеплення Psi02 .0000E+00
Потокощеплення PsiD .0000E+00
Напруженість H1 .0000E+00
Напруженість H2 .0000E+00
| T,c | I1,A | I0,A | B1,Тл | B2,Тл | Psi11 | Psi12 | Psi01 | Psi02 | PsiD | H1 | H2 |
.0000000E+00 .0000000E+00 .0000000E+00 .0000000E+00 .0000000E+00 .0000000E+00 .0000000E+00 .0000000E+00 .0000000E+00 .0000000E+00 .0000000E+00 .0000000E+00
.1000000E-03 .5378648E-04 .1655885E-08 .8817438E-02 -.8817438E-02 .3967853E-02 .3967853E-02 .2645231E-03 -.2645231E-03 .1490296E-08 .2689319E+00 -.2689319E+00
.2000000E-03 .1104865E-03 .4958149E-08 .1811248E-01 -.1811248E-01 .8150627E-02 .8150627E-02 .5433744E-03 -.5433744E-03 .4462334E-08 .5524306E+00 -.5524306E+00
.3000000E-03 .1700304E-03 .9896266E-08 .2787372E-01 -.2787372E-01 .1254319E-01 .1254319E-01 .8362118E-03 -.8362118E-03 .8906640E-08 .8501486E+00 -.8501486E+00
.4000000E-03 .2323454E-03 .1645781E-07 .3808922E-01 -.3808922E-01 .1714017E-01 .1714017E-01 .1142677E-02 -.1142677E-02 .1481203E-07 .1161721E+01 -.1161721E+01
.5000000E-03 .2973552E-03 .2463329E-07 .4874648E-01 -.4874648E-01 .2193595E-01 .2193595E-01 .1462394E-02 -.1462394E-02 .2216996E-07 .1486768E+01 -.1486768E+01
.6000000E-03 .3649803E-03 .3442191E-07 .5983246E-01 -.5983246E-01 .2692465E-01 .2692465E-01 .1794974E-02 -.1794974E-02 .3097972E-07 .1824890E+01 -.1824890E+01
.7000000E-03 .4351384E-03 .4580023E-07 .7133365E-01 -.7133365E-01 .3210019E-01 .3210019E-01 .2140010E-02 -.2140010E-02 .4122021E-07 .2175676E+01 -.2175676E+01
.8000000E-03 .5077436E-03 .5876601E-07 .8323602E-01 -.8323602E-01 .3745626E-01 .3745626E-01 .2497081E-02 -.2497081E-02 .5288940E-07 .2538698E+01 -.2538698E+01
.9000000E-03 .5827076E-03 .7331183E-07 .9552503E-01 -.9552503E-01 .4298633E-01 .4298633E-01 .2865751E-02 -.2865751E-02 .6598065E-07 .2913514E+01 -.2913514E+01
.9999999E-03 .6599390E-03 .8941971E-07 .1081857E+00 -.1081857E+00 .4868365E-01 .4868365E-01 .3245572E-02 -.3245572E-02 .8047775E-07 .3299665E+01 -.3299665E+01
.1100000E-02 .7393439E-03 .1070801E-06 .1212027E+00 -.1212027E+00 .5454130E-01 .5454130E-01 .3636082E-02 -.3636082E-02 .9637206E-07 .3696684E+01 -.3696684E+01
.1200000E-02 .8208257E-03 .1262857E-06 .1345601E+00 -.1345601E+00 .6055215E-01 .6055215E-01 .4036804E-02 -.4036804E-02 .1136571E-06 .4104086E+01 -.4104086E+01
.1300000E-02 .9042857E-03 .1470175E-06 .1482417E+00 -.1482417E+00 .6670888E-01 .6670888E-01 .4447252E-02 -.4447252E-02 .1323157E-06 .4521379E+01 -.4521379E+01
.1400000E-02 .9896230E-03 .1692877E-06 .1622309E+00 -.1622309E+00 .7300401E-01 .7300401E-01 .4866927E-02 -.4866927E-02 .1523589E-06 .4948059E+01 -.4948059E+01
.1500000E-02 .1076736E-02 .1925677E-06 .1765106E+00 -.1765106E+00 .7942989E-01 .7942989E-01 .5295319E-02 -.5295319E-02 .1733110E-06 .5383615E+01 -.5383615E+01
.1600000E-02 .1165522E-02 .2183813E-06 .1910636E+00 -.1910636E+00 .8597872E-01 .8597872E-01 .5731907E-02 -.5731907E-02 .1965431E-06 .5827538E+01 -.5827538E+01
.1700000E-02 .1255881E-02 .2451758E-06 .2058720E+00 -.2058720E+00 .9264253E-01 .9264253E-01 .6176160E-02 -.6176160E-02 .2206582E-06 .6279321E+01 -.6279321E+01
.1800000E-02 .1347716E-02 .2734619E-06 .2209180E+00 -.2209180E+00 .9941325E-01 .9941325E-01 .6627541E-02 -.6627541E-02 .2461157E-06 .6738489E+01 -.6738489E+01
.1899999E-02 .1440943E-02 .3032292E-06 .2361833E+00 -.2361833E+00 .1062826E+00 .1062826E+00 .7085500E-02 -.7085500E-02 .2729062E-06 .7204612E+01 -.7204612E+01
.1999999E-02 .1535494E-02 .3344503E-06 .2516494E+00 -.2516494E+00 .1132424E+00 .1132424E+00 .7549482E-02 -.7549482E-02 .3010053E-06 .7677357E+01 -.7677357E+01
.2099999E-02 .1631335E-02 .3671479E-06 .2672974E+00 -.2672974E+00 .1202840E+00 .1202840E+00 .8018923E-02 -.8018923E-02 .3304331E-06 .8156552E+01 -.8156552E+01
.2200000E-02 .1728486E-02 .4012934E-06 .2831084E+00 -.2831084E+00 .1273990E+00 .1273990E+00 .8493253E-02 -.8493253E-02 .3611640E-06 .8642298E+01 -.8642298E+01
.2300000E-02 .1827053E-02 .4368939E-06 .2990631E+00 -.2990631E+00 .1345786E+00 .1345786E+00 .8971895E-02 -.8971895E-02 .3932045E-06 .9135118E+01 -.9135118E+01
.2400000E-02 .1927271E-02 .4739271E-06 .3151422E+00 -.3151422E+00 .1418142E+00 .1418142E+00 .9454265E-02 -.9454265E-02 .4265343E-06 .9636195E+01 -.9636195E+01
.2500000E-02 .2029573E-02 .5114013E-06 .3313259E+00 -.3313259E+00 .1490968E+00 .1490968E+00 .9939777E-02 -.9939776E-02 .4602611E-06 .1014770E+02 -.1014770E+02
.2600000E-02 .2134682E-02 .5523138E-06 .3475944E+00 -.3475944E+00 .1564177E+00 .1564177E+00 .1042783E-01 -.1042783E-01 .4970825E-06 .1067322E+02 -.1067322E+02
.2700000E-02 .2243727E-02 .5935955E-06 .3639276E+00 -.3639276E+00 .1637677E+00 .1637677E+00 .1091783E-01 -.1091783E-01 .5342359E-06 .1121844E+02 -.1121844E+02
.2800000E-02 .2358415E-02 .6352723E-06 .3803052E+00 -.3803052E+00 .1711376E+00 .1711376E+00 .1140916E-01 -.1140916E-01 .5717450E-06 .1179186E+02 -.1179186E+02
.2900000E-02 .2481231E-02 .6804038E-06 .3967066E+00 -.3967066E+00 .1785182E+00 .1785182E+00 .1190120E-01 -.1190120E-01 .6123634E-06 .1240593E+02 -.1240593E+02
За допомогою файла rez_o.dat результати зображені у вигляді графіків:
Характеристика намагнічування сталі
Вебер-амперна характеристика дроселя
Залежність від часу напруг живлення подавача частоти
Знайдені залежності струмів у часі
Знайдені залежності індукцій у часі
Знайдені залежності потокозчеплень у часі
Знайдена залежність потокозчеплення дроселя у часі
Знайдені залежності напружень у часі
Ділись своїми роботами та отримуй миттєві бонуси!
0%
Оголошення від адміністратора