Курсовая работа: Электродинамический преобразователь энергии с тиристорной схемой питания
6.2 Программа (написана на языке TURBO
BASIC)
LET
h = .00001
Bo
= 1
La
= .00235
m1
= 100
D
= .8
w
= 40
C
= .00015
ksi
= 2 * 10 ^ 7
nu
= 4000
E
= 10
R
= .00105
Ra
= .05
Pi
= 3.141592654#
z
= 1
t1
= .00007
t2
= .00621:
t56
= t1
integral
= 0
integral2
= 0
LET
schet = 1
INPUT
"параметры выводить на экран? n-нет";
q1$
IF
q1$ = "n" OR q1$ = "N" THEN q = 0 ELSE q = 1
SCREEN
12
Uc
= 0
LOCATE
1, 45: PRINT "время
t(мс)"
LOCATE
2, 45: PRINT "белая линия
I(А)"
LOCATE
3, 45: PRINT "синяя линия
Uc(В)"
LOCATE
4, 45: PRINT "фиолетовая линия X(мм)"
LOCATE
5, 45: PRINT "красная линия V(мм/с)"
LOCATE
6, 45: PRINT "зеленая линия
a(m/S)"
LOCATE
7, 45: PRINT "коричневая линия Ua(В)"
0
LET i1 = I
LET Uc1 = Uc
LET x1 = X
LET V1 = V
LET xc = 0
1
LET k1i = E / La - (R / La) * i1 - Uc / La - (Bo * Pi * D * w) * V / La
LET i1 = I + h * .5 * k1i
LET k2i = E / La - (R / La) * i1 - Uc / La - (Bo * Pi * D * w) * V / La
LET i1 = I + k2i * h * .5
LET k3i = E / La - (R / La) * i1 - Uc / La - (Bo * Pi * D * w) * V / La
LET i1 = I + k3i * h
LET k4i = E / La - (R / La) * i1 - Uc / La - (Bo * Pi * D * w) * V / La
LET di = h * (2 * k2i + k1i + 2 * k3i + k4i) / 6
2
LET k1Uc = I / C
LET Uc1 = Uc + h * .5 * k1Uc
LET k2Uc = I / C
LET Uc1 = Uc + h * .5 * k2Uc
LET k3Uc = I / C
LET Uc1 = Uc + h * k3Uc
LET k4Uc = I / C
LET dUc = h * (2 * k2Uc + k1Uc + 2 * k3Uc + k4Uc) / 6
3
LET k1x = V
LET x1 = X + k1x * h * .5
LET k2x = V
LET x1 = X + k2x * h * .5
LET k3x = V
LET x1 = X + k3x * h
LET k4x = V
LET dx = (k1x + 2 * k2x + 2 * k3x + k4x) * h / 6
4
LET k1V = (Bo * Pi * D * w * I) / m1 - (nu / m1) * V1 - (ksi / m1) * X
LET V1 = V + k1V * h * .5
LET k2V = (Bo * Pi * D * w * I) / m1 - (nu / m1) * V1 - (ksi / m1) * X
LET V1 = V + k2V * h * .5
LET k3V = (Bo * Pi * D * w * I) / m1 - (nu / m1) * V1 - (ksi / m1) * X
LET V1 = V + k3V * h
LET k4V = (Bo * Pi * D * w * I) / m1 - (nu / m1) * V1 - (ksi / m1) * X
LET dv = h * (2 * k2V + k1V + 2 * k3V + k4V) / 6
LET I = I + di
LET Uc = Uc + dUc
LET X = X + dx
LET V = V + dv
LET integral2 = integral2 + ABS(I) * h
LET integral = integral + h * (ABS(I)) ^ 2
LET a = (Bo * Pi * D * w * I - nu * V - ksi * X) / m1
LET Ua = La * di / h + Bo * Pi * D * w * V + I * Ra
LET tall = tall + h
IF q = 1 THEN
LOCATE 1, 1: PRINT "t="; tall * 1000, " "
LOCATE 2, 1: PRINT "I="; I, " "
LOCATE 3, 1: PRINT "Uc="; Uc, " "
LOCATE 4, 1: PRINT "X="; X * 1000, " "
LOCATE 5, 1: PRINT "V="; V * 1000, " "
LOCATE 6, 1: PRINT "a="; a, " "
LOCATE 7, 1: PRINT "Ua="; Ua, " "
END IF
PSET (t * 20000 * .01 / t2, 250 - I * .5)
PSET (t * 20000 * .01 / t2, 250 - Uc * .1), 3
PSET (t * 20000 * .01 / t2, 250)
PSET (t * 20000 * .01 / t2, 250 - V * 100), 4
PSET (t * 20000 * .01 / t2, 250 - X * 100000), 5
PSET (t * 20000 * .01 / t2, 250 - a / 3), 2
PSET (t * 20000 * .01 / t2, 250 - Ua * .1), 6
LET t = t + h
IF t > t1 THEN
IF (((I - ikontr) / I) < .001) AND (param = 0) THEN
LOCATE 14, 45
PRINT "Тпп(мс)=";
tall * 1000; : INPUT zxc
LET param = 1
END IF
LET ikontr = I
GOTO 5
END IF
GOTO
0
5
i2 = I
V2 = V
x2 = X
7
LET k1i = E / La - (R / La) * i2 - (Bo * Pi * D * w * V / La)
LET i2 = I + h * .5 * k1i
LET k2i = E / La - (R / La) * i2 - (Bo * Pi * D * w * V / La)
LET i2 = I + k2i * h * .5
LET k3i = E / La - (R / La) * i2 - (Bo * Pi * D * w * V / La)
LET i2 = I + k3i * h
LET k4i = E / La - (R / La) * i2 - (Bo * Pi * D * w * V / La)
LET di = h * (2 * k2i + k1i + 2 * k3i + k4i) / 6
8
LET k1x = V
LET x2 = X + k1x * h * .5
LET k2x = V
LET x2 = X + k2x * h * .5
LET k3x = V
LET x2 = X + k3x * h
LET k4x = V
LET dx = (k1x + 2 * k2x + 2 * k3x + k4x) * h / 6
9
LET k1V = (Bo * Pi * D * w * I) / m1 - (nu / m1) * V2 - (ksi / m1) * X
LET V2 = V + k1V * h * .5
LET k2V = (Bo * Pi * D * w * I) / m1 - (nu / m1) * V2 - (ksi / m1) * X
LET V2 = V + k2V * h * .5
LET k3V = (Bo * Pi * D * w * I) / m1 - (nu / m1) * V2 - (ksi / m1) * X
LET V2 = V + k3V * h
LET k4V = (Bo * Pi * D * w * I) / m1 - (nu / m1) * V2 - (ksi / m1) * X
LET dv = h * (2 * k2V + k1V + 2 * k3V + k4V) / 6
10
LET I = I + di
LET X = X + dx
LET V = V + dv
LET integral2 = integral2 + ABS(I) * h
LET integral = integral + h * (ABS(I)) ^ 2
LET a = (Bo * Pi * D * w * I - nu * V - ksi * X) / m1
LET Ua = La * di / h + Bo * Pi * D * w * V + I * Ra
LET tall = tall + h
11
IF q = 1 THEN
LOCATE 1, 1: PRINT "t="; tall * 1000, " "
LOCATE 2, 1: PRINT "I="; I, " "
LOCATE 3, 1: PRINT "Uc="; Uc, " "
LOCATE 4, 1: PRINT "X="; X * 1000, " "
LOCATE 5, 1: PRINT "V="; V * 1000, " "
LOCATE 6, 1: PRINT "a="; a, " "
LOCATE 7, 1: PRINT "Ua="; Ua, " "
END IF
PSET (t * 20000 * .01 / t2, 250 - I * .5)
PSET (t * 20000 * .01 / t2, 250 - Uc * .1), 3
PSET (t * 20000 * .01 / t2, 250)
PSET (t * 20000 * .01 / t2, 250 - V * 100), 4
PSET (t * 20000 * .01 / t2, 250 - X * 100000), 5
PSET (t * 20000 * .01 / t2, 250 - a / 3), 2
PSET (t * 20000 * .01 / t2, 250 - Ua * .1), 6
LET t = t + h
IF t > t1 + t2 THEN GOTO 12
GOTO 5
12
LET Uc = -Uc:
LET tkontr = t:
121
LET Uc3 = Uc
LET i3 = I
LET x3 = X
LET V3 = V
13
LET k1i = E / La - (R / La) * i3 - (Uc / La) - (Bo * Pi * D * w / La) * V
LET i3 = I + h * .5 * k1i
LET k2i = E / La - (R / La) * i3 - (Uc / La) - (Bo * Pi * D * w / La) * V
LET i3 = I + k2i * h * .5
LET k3i = E / La - (R / La) * i3 - (Uc / La) - (Bo * Pi * D * w / La) * V
LET i3 = I + k3i * h
LET k4i = E / La - (R / La) * i3 - (Uc / La) - (Bo * Pi * D * w / La) * V
LET di = h * (2 * k2i + k1i + 2 * k3i + k4i) / 6
14
LET k1Uc = I / C
LET Uc3 = Uc + h * .5 * k1Uc
LET k2Uc = I / C
LET Uc3 = Uc + h * .5 * k2Uc
LET k3Uc = I / C
LET Uc3 = Uc + h * k3Uc
LET k4Uc = I / C
LET dUc = h * (2 * k2Uc + k1Uc + 2 * k3Uc + k4Uc) / 6
15
LET k1x = V
LET x3 = X + k1x * h * .5
LET k2x = V
LET x3 = X + k2x * h * .5
LET k3x = V
LET x3 = X + k3x * h * .5
LET k4x = V
LET dx = (k1x + 2 * k2x + 2 * k3x + k4x) * h / 6
16
LET k1V = (Bo * Pi * D * w * I) / m1 - (nu / m1) * V3 - (ksi / m1) * X
LET V3 = V + k1V * h * .5
LET k2V = (Bo * Pi * D * w * I) / m1 - (nu / m1) * V3 - (ksi / m1) * X
LET V3 = V + k2V * h * .5
LET k3V = (Bo * Pi * D * w * I) / m1 - (nu / m1) * V3 - (ksi / m1) * X
LET V3 = V + k3V * h
LET k4V = (Bo * Pi * D * w * I) / m1 - (nu / m1) * V3 - (ksi / m1) * X
LET dv = h * (2 * k2V + k1V + 2 * k3V + k4V) / 6
17
LET I = I + di
LET Uc = Uc + dUc
LET X = X + dx
LET V = V + dv
LET integral2 = integral2 + ABS(I) * h
LET integral = integral + h * (ABS(I)) ^ 2
LET a = (Bo * Pi * D * w * I - nu * V - ksi * X) / m1
LET Ua = La * di / h + Bo * Pi * D * w * V + I * Ra
LET tall = tall + h
IF Uc > 0 AND xc = 0 THEN
LOCATE 8, 45
LET tvost = (t - tkontr) * 10 ^ 6
PRINT "tvost(¬Є‘)=";
(t - tkontr) * 10 ^ 6
LET xc = 1
END IF
Страницы: 1, 2, 3, 4 |
