МІНІСТЕРСТВО ОСВІТИ І НАУКИ УКРАЇНИ
НАЦІОНАЛЬНИЙ УНІВЕРСИТЕТ “Львівська політехніка”
АНАЛІЗ ЛІНІЙНИХ ЕЛЕКТРИЧНИХ КІЛ
ЗАВДАННЯ ТА МЕТОДИЧНІ НАСТАНОВИ
ДО КОМПЛЕКСНОЇ РОЗРАХУНКОВОЇ РОБОТИ
для студентів базового напрямку
6.1601 „Інформаційна безпека”
Затверджено
на засіданні кафедри теоретичної
та загальної електротехніки
протокол № від . .2007 р.
Львів – 2007
Завдання та методичні настанови з дисципліни “Теорія електричних кіл ” для студентів базового напряму „Інформаційна безпека”/ Укл. Гамола О.Є., Соколовський М.О. – Львів: Видавництво Національного університету “Львівська політехніка” 2007, – __ с.
Укладачі Гамола О.Є., канд.техн.наук, доц.,
Соколовський М.О. канд.техн.наук, доц.
Відповідальний за випуск Стахів П.Г., д-р техн. наук, проф.
Рецензенти Пеленський Р.А., д-р, техн.наук, проф.
АНАЛІЗ ЛІНІЙНИХ ЕЛЕКТРИЧНИХ КІЛ
Мета: Набуття практичних навиків розрахунку лінійних електричних кіл та аналізу перехідних процесів в них.
В завданні до комплексної розрахункової роботи (КРР) подано варіантів. Варіант завдання видається викладачем. Результати виконання КРР оформляють у вигляді пояснювальної записки, що має титульний лист, де вказують назви міністерства, університету, кафедри, тему курсової роботи, назву дисципліни, хто виконав роботу: спеціальність, група, ПІБ студента, хто прийняв: посада та ПІБ викладача. Пояснювальну записку та титульний лист оформляють на аркушах паперу формату 297х210 (А4).
Матеріал пояснювальної записки рекомендується викладати послідовно згідно з виданим завданням. В тексті записки формули подаються в загальному вигляді і лише потім в них підставляють числові значення. Графічний матеріал записки (схеми, графіки тощо), рекомендується виконувати на міліметровому папері формату А4. Його розташовують у місцях відповідного текстового матеріалу. До рисунків роблять підписи з тематичною назвою. В кінці пояснювальної записки приводять список використаної літератури.
Задача 1. Для заданого лінійного електричного кола (табл. 1, рис. 1–33) з джерелами змінної напруги EMBED Equation.3 та струму EMBED Equation.3 частотою f=50 Гц:
1. Записати систему рівнянь за методом контурних струмів та обчислити коефіцієнти цієї системи лінійних алгебричних рівнянь.
2. Записати систему рівнянь за методом вузлових напруг та обчислити коефіцієнти цієї системи лінійних алгебричних рівнянь.
3. Сформувати відповідні матриці та записати матричне рівняння методу контурних струмів.
4. Розв’язати матричне рівняння методу контурних струмів в середовищі MATLAB та визначити струми віток електричного кола.
5. Сформувати відповідні матриці та записати матричне рівняння методу вузлових напруг.
6. Розв’язати матричне рівняння методу вузлових напруг в середовищі MATLAB та визначити струми віток електричного кола.
7. Скласти рівняння балансу потужностей.
Задача 2. В заданому електричному колі (табл. 2, рис. 34–93) відбувається комутація. Параметри кола подано в табл. 2. Потрібно:
1. Визначити перехідну величину (струм або напругу), вказану в табл. 2, після комутації класичним методом при дії джерела постійної напруги або джерела постійного струму.
2. За отриманим аналітичним виразом побудувати графік перехідної величини в інтервалі часу від t=0 до t=5/||. На графіках показати окремо вільну та вимушену складові перехідної величини. Ліворуч від осі ординат зобразити шукану величину до комутації.
6. За допомогою чисельного методу розрахувати і побудувати шукану перехідну величину для електричного кола у випадку генерування джерелом енергії імпульсного сигналу (табл.3, рис. 1–8).
Примітка:
1. Тривалість імпульсу прийняти tім=2/| |.
2. Амплітуда імпульсу Um або Jm дорівнюють відповідним значенням Е та J з табл. 2.
МЕТОДИЧНІ НАСТАНОВИ
Аналіз усталеного режиму лінійних електричних кіл здійснюється переважно методом контурних струмів або методом вузлових напруг.
Система рівнянь за методом контурних струмів має такий вигляд:
EMBED Equation.3 ;
EMBED Equation.3 ;
; (1)
EMBED Equation.3 ;
де N=p-(q-1) – кількість незалежних (головних) контурів; p – кількість віток електричного кола (без віток з джерелами струму); q – кількість вузлів електричного кола; Z11, ZNN – власні опори контурів; Z12, ZN-1N – спільні опори контурів N-1 і N; ІK1, IKN – комплексні діючі значення контурних струмів; ЕK1, ЕKN – комплексні діючі значення контурних ЕРС.
За методом вузлових напруг записується така система рівнянь:
EMBED Equation.3 ;
EMBED Equation.3 ;
; (2)
EMBED Equation.3 ;
де Q=q-1 – кількість незалежних вузлів; Y11, YQQ – власні вузлові провідності; Y12, YQ-1Q – спільні вузлові провідності вузлів Q-1 і Q; U10, UQ0 – комплексні діючі значення вузлових напруг (відносно базового); I11, IQQ – комплексні діючі значення вузлових струмів.
Матричне рівняння методу контурних струмів
EMBED Equation.3 ; (3)
де EMBED Equation.3 – контурна матриця; EMBED Equation.3 – матриця повних опорів віток електричного кола; EMBED Equation.3 –матриця-стовпець контурних струмів; EMBED Equation.3 –матриця-стовпець ЕРС джерел напруг; EMBED Equation.3 –матриця-стовпець струмів джерел струмів.
Рівняння методу вузлових напруг у матричній формі має вигляд:
EMBED Equation.3 ; (4)
де EMBED Equation.3 – вузлова матриця (матриця інциденцій); EMBED Equation.3 – матриця повних провідностей віток електричного кола; EMBED Equation.3 –матриця-стовпець вузлових напруг (відносно базового); EMBED Equation.3 –матриця-стовпець ЕРС джерел напруг; EMBED Equation.3 –матриця-стовпець струмів джерел струмів.
Розглянемо як записуються рівняння за цими методами і як формуються відповідні матриці на такому прикладі.
Приклад 1. На рис.1 подано схему електричного кола з такими параметрами: R1=10 Ом; R2=5 Ом; XL1=10 Ом; XL2=5 Ом; XC=10 Ом; e1=502sin(314t+45) В; j5=22sin(314t) А.
EMBED Visio.Drawing.11
Рис.1 Схема електричного кола
Проведемо аналіз електричного кола методом контурних струмів.
Визначаємо: кількість віток (без джерел струму ДС) p=4; кількість вузлів q=3; кількість незалежних контурів N=p-(q-1)=2. Струм контура, утвореного віткою з ДС є відомим EMBED Equation.3 . Цей контур показує шлях замикання струму ДС.
Запишемо систему рівнянь за методом контурних струмів:
EMBED Visio.Drawing.11
EMBED Visio.Drawing.11
Рис. 2. Схема електричного кола для методу контурних струмів
EMBED Equation.3 ; (5)
EMBED Equation.3 ;
де EMBED Equation.3 ; EMBED Equation.3 ; EMBED Equation.3 ; EMBED Equation.3 Ом; EMBED Equation.3 Ом;
EMBED Equation.3 Ом; EMBED Equation.3 Ом; EMBED Equation.3 Ом; EMBED Equation.3 Ом
EMBED Equation.3 Ом; EMBED Equation.3 В; EMBED Equation.3 В.
Сформуємо необхідні матриці для запису матричного рівняння контурних струмів на підставі графа заданого електричного кола (рис. 5). Під час формування контурної матриці EMBED Equation.3 дотримуємося правила: +1 записуємо в клітинках, утворених перетином рядка k і стовця p, якщо напрям обходу контура k і напрям вітки p збігаються; –1 – якщо напрям обходу контура k і напрям вітки p не збігаються; 0 – якщо вітка p не належить контуру k. В матриці-стовпці EMBED Equation.3 ЕРС джерел напруг записуємо +Еp, якщо напрям вітки і напрям ЕРС в цій вітці збігаються; – Еp, якщо напрям вітки і напрям ЕРС в цій вітці не збігаються; 0 – ЕРС в даній вітці відсутня. Аналогічно формується матриця-стовпець EMBED Equation.3 струмів ДС: +Ji – записуємо в рядку p, якщо вітка p входить в контур утворений ДС і напрям обходу цього контура й напрям вітки p збігаються; – Ji – записуємо в рядку p, якщо вітка p входить в контур утворений ДС і напрям обходу цього контура й напрям вітки p не збігаються; 0 – в рядку p, якщо вітка p не входить в контур ДС.
EMBED Visio.Drawing.11
Рис. 5. Граф електричного кола для формування матриць EMBED Equation.3 , EMBED Equation.3 ; EMBED Equation.3 .
На підставі цього записуємо відповідні матриці з врахуванням того, що контур ДС замикається через вітки 3 і 4
EMBED Equation.3 ; EMBED Equation.3 ; (6)
EMBED Equation.3 ; EMBED Equation.3 ; EMBED Equation.3
Далі формуємо m-файл для розв’язку матричного рівняння методу контурних струмів в середовищі MATLAB
loop_currenst_method.m
G=[1 1 0 1; 0 -1 1 0] % контурна матриця
Z1=10+i*10 % елемент матриці Z
Z2=i*5 % – // –
Z3=-i*10 % – // –
Z4=5 % – // –
E1=50+i*50 % елемент матриці Е
J5=2 % елемент матриці J
Z=[Z1 0. 0. 0.; 0. Z2 0. 0.; 0. 0. Z3 0.; 0. 0. 0. Z4] % матриця опорів віток
E=[E1; 0.; 0.; 0.] % матриця ЕРС ДН
J=[0.; 0.; J5; J5] % матриця струмів ДС
ZG=Z*G’
ZV=G*ZG % матриця контурних опорів
ZJ=Z*J
EZJ=E-ZJ
B=G*EZJ
X=ZV\B % матриця контурних струмів
IV=G’*X+J % матриця струмів віток
Під час розрахунку отримано такі результати
EMBED Equation.3 ; EMBED Equation.3 ; EMBED Equation.3
EMBED Equation.3
Проведемо аналіз електричного кола методом вузлових напруг:
Кількість незалежних рівнянь за цим методом Q=q–1=2. За базовий вибираємо вузол “0”. Запишемо систему рівнянь за методом вузлових напруг
EMBED Equation.3 ; (7)
EMBED Equation.3 ;
де EMBED Equation.3 ; EMBED Equation.3 ; EMBED Equation.3 ; EMBED Equation.3 См; EMBED Equation.3 ; EMBED Equation.3 ; EMBED Equation.3 ; EMBED Equation.3 ; EMBED Equation.3 ; EMBED Equation.3 ;
EMBED Equation.3 ; EMBED Equation.3 .
EMBED Visio.Drawing.11
EMBED Visio.Drawing.11
Рис. 6. Схема електричного кола для формування матриць EMBED Equation.3 , EMBED Equation.3 ; EMBED Equation.3 .
Сформуємо необхідні матриці для запису матричного рівняння вузлових напруг на підставі графа заданого електричного кола (рис. 6). Під час формування вузлової матриці EMBED Equation.3 дотримуємося правила: +1 записуємо в клітинках, утворених перетином рядка q і стовця p, якщо вітка p спрямована від вузла q; –1 – якщо вітка p спрямована до вузла q; 0 – якщо вітка p не належить до вузла q. Матриці-стовпці ЕРС EMBED Equation.3 та ДС EMBED Equation.3 формуються аналогічно як у методі контурних струмів.
EMBED Visio.Drawing.11
Рис. 7. Граф електричного кола для методу вузлових напруг
На підставі цього записуємо відповідні матриці методу вузлових напруг
EMBED Equation.3 ; EMBED Equation.3 ; (8)
EMBED Equation.3 ; EMBED Equation.3 ; EMBED Equation.3
Далі формуємо m-файл для розв’язку матричного рівняння методу вузлових напруг в середовищі MATLAB
nodal_voltage_method.m
P=[-1 1 1 0; 0 -1 -1 1] % вузлова матриця
Y1=0.05-i*0.05 % елемент матриці Y
Y2=-i*0.2 % – // –
Y3=i*0.1 % – // –
Y4=0.2 % – // –
E1=50+i*50 % елемент матриці Е
J5=2. % елемент матриці J
Y=[Y1 0. 0. 0.; 0. Y2 0. 0.; 0. 0. Y3 0.; 0. 0. 0. Y4] % матриця провідностей віток
V=[E1; 0.; 0.; 0.] % матриця ЕРС ДН
C=[0.; 0.; J5; J5] % матриця струмів ДС
YP=Y*P'
A=P*YP % матриця вузлових провідностей
YV=Y*V
CYV=C-YV
B=P*CYV
X=A\B % матриця вузлових напруг
IV=Y*(P'*X+V) % матриця струмів віток
Під час розрахунку отримано такі результати
EMBED Equation.3 ; EMBED Equation.3 ; EMBED Equation.3
EMBED Equation.3
Аналіз перехідного процесу лінійного електричного кола відповідно з умовою задачі 2 проведемо класичним методом.
Класичний метод аналізу перехідного процесу лінійного електричного кола полягає в знаходженні розв’язку лінійного неоднорідного диференційного рівняння EMBED Equation.3, як суми часткового розв’язку лінійного неоднорідного диференційного рівняння – вимушеної складової перехідної величини EMBED Equation.3 і загального розв’язку лінійного однорідного диференційного рівняння – вільної складової перехідної величини EMBED Equation.3
EMBED Equation.3. (9)
Перехідними величинами, щодо яких формуються диференційні рівняння, є струм il у вітці з котушкою індуктивності та напруга uc на конденсаторі
EMBED Equation.3, (10)
для яких початкові умови визначаються за 1-им та 2-им законами комутації відповідно
EMBED Equation.3. (11)
Розрахунок перехідного процесу класичним методом проводимо за алгоритмом, поданим блок схемою.
EMBED Visio.Drawing.11
Чисельні методи аналізу перехідного процесу в лінійному електричному колі передбачають інтегрування диференційних рівнянь, складених методом змінних стану, у формі Коші EMBED Equation.3. Для отримання єдиного розв’язку систему диференційних рівнянь необхідно доповнити початковими умовами EMBED Equation.3. Суть чисельних методів аналізу полягає у послідовному обчисленні наближених значень вектор-функції EMBED Equation.3 на множині точок аргументу EMBED Equation.3 в інтервалі [a, b] визначення вектора-функції EMBED Equation.3 при EMBED Equation.3 з кроком EMBED Equation.3. Відомі різноманітні методи чисельного розв’язку рівнянь у формі Коші. Наприклад, рекурентна формула для однокрокового методу Рунге-Кутта другого порядку
EMBED Equation.3, (12)
де EMBED Equation.3; EMBED Equation.3,
чи однокрокового неявного методу Ейлера
EMBED Equation.3, (13)
для реалізації якого застосовують метод Ньютона на (l+1)-му кроці
EMBED Equation.3; (14)
EMBED Equation.3, (15)
де
EMBED Equation.3; EMBED Equation.3. (16)
Для здійснення ітерації за (16) на кожному кроці необхідно мати нульове наближення EMBED Equation.3, яке може бути обчислено на основі явної формули Ейлера
EMBED Equation.3. (17)
Розглянемо приклад розрахунку перехідного процесу електричного кола.
Приклад 2. Для заданого електричного кола (рис. 8) визначити струми у вітках після комутації. Е=50 В; R1=5 Ом; R2=8 Ом; R3=10 Ом; L=100 мГн; С=200 мкФ.
EMBED Visio.Drawing.11
Рис. 8. Схема електричного кола
Розрахуємо задачу класичним методом.
1. Сформуємо диференційне рівняння щодо струму iL заданого електричного кола на підставі законів Кірхгофа та рівнянь елементів:
Система рівнянь за законами Кірхгофа
i1–i2–il=0,
ur1+ ur2=E (18)
ur2–ul–ur3=0
Рівняння елементів
ur1=R1i1, ur2=R2i2, ur3= R3il, ul=Ldil/dt, (19)
Після відповідних перетворень отримаємо диференційне рівняння щодо перехідного струму il
Ldil/dt+(R1R2+R1R3+R2R3)/(R1+R2)il=ER2/(R1+R2). (20)
Визначаємо незалежну початкову умову il(0). Для цього необхідно розрахувати усталений режим кола до комутації з врахування того, що опір котушки в колі постійного струму в усталеному режимі дорівнює нулю (рис.9).
Оскільки вітка з R2 замкнута накоротко, то i1(0-)=iL(0-)=E/R1=50/5=10А. На підставі 1-го закону комутації отримуємо
iL(0)= iL(0-)=10 A.
EMBED Visio.Drawing.11
Рис. 9. Схема електричного кола для усталеного режиму до комутації
2. Визначаємо вимушену складову струму ilвим. Для цього необхідно розрахувати усталений режим кола після комутації з урахуванням того, що опір котушки в колі постійного струму в усталеному режимі дорівнює нулю (рис. 10). Вимушену складову струму ilвим можна також отримати з диференційного рівняння (20) стану електричного кола після комутації, враховуючи, що в усталеному режимі dilвим/dt=0.
ilвим=ER2/(R1R2+R1R3+R2R3)=508/(58+510+810)=2,35 A,
EMBED Visio.Drawing.11
Рис. 9. Схема електричного кола для усталеного режиму після комутації
3. Визначаємо вільну складову струму ilвіл Для цього спочатку сформуємо характеристичне рівняння.
3.1. Характеристичне рівняння може бути сформовано декількома способами. Для заданого кола розглянемо два способи: шляхом запису вхідного опору кола після комутації щодо будь-якої розімкненої вітки змінному струму (рис. 10); шляхом алгебризації ( EMBED Equation.3 ) диференційного однорідного рівняння.
Визначаємо корінь характеристичного рівняння
EMBED Equation.3,
EMBED Equation.3 с-1.
EMBED Visio.Drawing.11
Рис. 10. Схема для визначення вхідного опору кола
На підставі кореня характеристичного рівняння вільну складову перехідної величини iLвіл шукатимемо у такому вигляді
EMBED Equation.3, (22)
3.3. Визначаємо сталу інтегрування. Для цього запишемо рівняння для вільної складової перехідної величини iLвіл на момент часу t=0:
EMBED Equation.3. (23)
Значення величини iLвіл(0) знаходимо з рівняння
EMBED Equation.3А.
Звідки А1=7,65 А.
Отже, перехідна величина il описується таким виразом:
EMBED Equation.3. (24)
Решту перехідних величин знаходимо з рівнянь елементів та рівнянь за законами Кірхгофа:
i1–i2–il=0,
R2i2– Ldil/dt – R3il=0,
звідки EMBED Equation.3
EMBED Equation.3. (20)
Рис. 11. Часова діаграма перехідного струму il
Для розрахунку перехідного процесу чисельним методом запишемо диференційне рівняння щодо змінної стану iL (20) у формі Коші:
EMBED Equation.3
Обчислимо перехідний струм EMBED Equation.3 чисельним однокроковим методом Рунге-Кутта другого порядку за його рекурентною формулою
EMBED Equation.3, (12)
де EMBED Equation.3; EMBED Equation.3; EMBED Equation.3 ; EMBED Equation.3
Далі формуємо m-файл для розрахунку перехідного процесу чисельним методом в середовищі MATLAB
dtk=1e-4
t(1)=0.
t2=0.
R1=5
R2=8
R3=10
L=0.1
E=50
R123=R1*R2+R1*R3+R2*R3
R12=R1+R2
LR=1/(L*R12)
ilk=10.
ilk1(1)=10.
for i=2:400
k1k=dtk*LR*(-R123*ilk+E*R2)
k2k=dtk*LR*(-R123*(ilk+k1k)+E*R2)
ilk=ilk+(k1k+k2k)/2
ilk1(i)=ilk
t2=t2+dtk
t(i)=t(i-1)+dtk
i=i+1
end
plot(t,ilk1,'-k','LineWidth',2),'FontSize',14
xlabel('time t','FontSize',14);
ylabel('current il','FontSize',14);
grid on
Рис. 12. Часова діаграма перехідного струму iL , отримана чисельним методом
Таблиця 1
Таблиця 2
НАВЧАЛЬНЕ ВИДАННЯ
КОМПЛЕКСНА РОЗРАХУНКОВА РОБОТА
ЗАВДАННЯ ТА МЕТОДИЧНІ НАСТАНОВИ
до комплексної розрахункової роботи
з дисципліни “Теорія електричних кіл”
для студентів базового напрямку
6.1601 „Інформаційна безпека”
Укладачі Гамола Орест Євгенович Соколовський Михайло Олександрович
Редактор
Комп’ютерне складання