Смекни!
smekni.com

Розрахунок трьохфазного мостового випрямляча (стр. 5 из 5)

;

.

Отримане рівняння визначає стійкість у "малому" початкової нелінійної системи в околі установившогося режиму

.

Застосуємо дискретне перетворення Лапласа до отриманого рівняння:

,

де

– дискретне зображення Лапласа
;
– початкове значення вектору
:
. Перетворивши отримане рівняння отримаємо:

.

В результаті відкриття визначника

, знаходимо характеристичний многочлен:

;

.

Зробимо заміну змінних

;

;

,
.

Корені даного рівняння знаходяться в середині кола одиничного радіуса і тому можна зробити висновок, що дана схема являється стійкою.

6. Висновки

У даній курсовій роботі був виконаний розрахунок перехідного процесу, що виникає при включенні пристрою, і стійкості в "малому" сталому режимі, на прикладі трьохфазного мостового випрямляча із ШІМ першого роду. При виконанні курсової роботи була спроектована принципова схема пристрою. Схему такого випрямляча можна віднести до схем зі постійною структурою й аналізувати неї на інтервалах роботи, де змінюється напруга.. Для розрахунку перехідного процесу були складені диференціальні рівняння й по них побудовані графіки, елементи силової частини при цьому замінялися їхніми ідеалізованими моделями. Також був проведений розрахунок стійкості системи, що показав, що дана система є стійкою.


7. Література

1. Г.С. Найвельт и др. "Источники электропитания радиоэлектронной аппаратуры. Справочник". М.: Радио и связь, 1985

2. Перетворювальна техніка. Частина 1: Підручник./ В. С. Руденко, В. Я. Ромашко, В. Г. Морозов. – К.: ІСДО, 1996.– 262 с.

3. Замкнутые системы преобразования электрической энергии. Под ред. Д-ра техн. наук В.Я. Жуйкова. – Киев, "Техника", 1989.

4. Методические указания к курсовой работе по курсу "Преобразовательная техника" для студентов специальности "Промышленная электроника" всех форм обучения. Сост.: В. С. Руденко, В. Г. Морозов, В. Я. Ромашко. – Киев: КПИ, 1984. – 56с.


Додаток 1.


Додаток 2.

function main;

Rv1=37*10^(-3);

Rv2=Rv1;

Rdr=80*10^(-3);

Rd=0.648;

Ld=160*10^(-6);

C=39*10^(-3);

E=20*sqrt(6);

U0=0;

I0=0;

t=0;

K=-3500;

U_start=0;

Kfb=20;

Pi=3.1415;

L=Pi/7;

s=0;

Z=0;

Y=[0, 0];

t=0;

PhaseShift=2*Pi/3;

Ufb=0;

Uref=5;

D=0;

Dx=0;

Dt=0;

Dtx=0;

Ufb=0;

O=0;

P=0;

w=2*Pi*300;

for f=1:100

s=0;

Z=0;

Dtx=0;

Dx=0;

U_start=(Y(end,2)*5/36-Uref)*Kfb;

if (U_start<0)

U_start=0;

end

Dtx(1)=0;

Dtx(2)=-U_start/K;

if (Dtx(2)>((Pi/3)/w))

Dtx(2)=(Pi/3)/w;

end

Dx(1)=U_start;

Dx(2)=K*Dtx(2)+U_start;

Y0=[Y(end,1),Y(end,2)];

[s,Z]=ode113(@F_first, [0 (Dtx(2)+10^-8)], Y0);

s=s+t(length(t));

Y=[Y;Z];

t=[t;s];

n=0;

Ox=0;

Px=0;

for x=0:0.000001:Dtx(2)+10^-8

n=n+1;

Ox(n)=E*sin(w*x+PhaseShift);

Px(n)=x;

end

Px=Px+P(length(P));

P=[P, Px];

O=[O, Ox];

PhaseShift=PhaseShift-Pi/3;

PhaseShift=PhaseShift+(Dtx(2)+10^-8)*w;

Dtx(3)=(Pi/3)/w;

Dx(3)=0;

s=0;

Z=0;

Y0=[Y(end,1),Y(end,2)];

[s,Z]=ode113(@F_first, [0 (Dtx(3)-Dtx(2))+10^-8], Y0);

s=s+t(length(t));

Y=[Y;Z];

t=[t;s];

n=0;

Ox=0;

Px=0;

for x=0:0.000001:(Dtx(3)-Dtx(2))

n=n+1;

Ox(n)=E*sin(w*x+PhaseShift);

Px(n)=x;

end

Px=Px+P(length(P));

P=[P, Px];

O=[O, Ox];

PhaseShift=PhaseShift+(Dtx(3)-Dtx(2))*w;

Dtx=Dtx+Dt(length(Dt));

D=[D, Dx];

Dt=[Dt, Dtx];

end

for x=1:length(t)

Ufb(x)=(Y(x, 2)*5/36-Uref)*Kfb;

if (Ufb(x)<0)

Ufb(x)=0;

end

end

I=Y;

I(:, 2)=[];

subplot(4,1,1);

plot(t,I);

xlabel('t(s)');

ylabel('I(A)');

title('Inductor Current (A)');

hold on;

I=Y;

I(:, 1)=[];

subplot(4,1,2);

plot(t,I);

title('Output Voltage (V)');

xlabel('t(s)');

ylabel('U(V)');

subplot(4,1,3);

plot(P,O);

xlabel('t(s)');

ylabel('U(V)');

title('Secondary Winding Voltage (V)');

subplot(4,1,4);

plot(t,Ufb);

hold on;

plot(Dt, D);

title('Ufb (V)');

xlabel('t(s)');

ylabel('U(V)');

function G = F_first(t, y)

Ud=E*sin(w*t+PhaseShift);

G=[-1/Ld*(-Ud+y(1)*(Rv1+Rv2+Rdr)+y(2)); 1/C*(y(1)-y(2)/Rd)];

end

end