Моделирование электротехнических устройств. Учебно-методическое пособие для выполнения контрольных работ по дисциплине «Моделирование электротехническ
.pdf
Уравнения (2) примут вид системы: ax1 + b = y1
ax2 + b = y2
решая ее, находим коэффициенты a и b.
С целью определения отклонений нужно составить таблицу 2:
Таблица 2
X
Y
Y*=aX+b
ε = |Y-Y*|
Аналогичной таблицей завершаются вычисления и при использовании метода средних и метода наименьших квадратов.
Определение параметров эмпирической зависимости по методу средних
Метод состоит в том, что параметры ao, a1,..., am зависимости (1) определяются с использованием условия равенства нулю суммы отклонений во всех точках xi:
n |
n |
yi 0 |
|
i |
xi , a0 , a1 ,..., am |
(3) |
|
i 0 |
i 0 |
|
|
Полученное уравнение |
служит для определения |
параметров |
|
а0,a1,...,am. Ясно, что из одного уравнения нельзя однозначно определить все m+1 параметров. Однако, поскольку других условий нет, равенство (3), путем группировки отклонений i разбивается на систему, состоящую из m+1 уравнений. Например,
.............
n-1 n
Решая эту систему уравнений, можно найти неизвестные параметры.
В процессе группировки необходимо все отклонения располагать по возрастанию номеров. В каждом уравнении должно быть примерно одинаковое количество отклонений.
Пример: Тело, движущееся прямолинейно с неизвестной скоростью Vo , в момент времени t=0 начинает тормозить. Измерялось расстояние Х от начала торможения в следующие моменты времени
t[сек] |
0 |
5 |
10 |
15 |
20 |
25 |
|
|
|
|
|
|
|
X [м] |
0 |
106 |
182 |
234 |
261 |
275 |
|
|
|
|
|
|
|
|
|
|
11 |
|
|
|
Считая движение тела равнозамедленным с постоянным замедлением а, найти приближенные значения Vo и а.
Искомые параметры могут быть найдены из уравнения движения тела, которое представим с помощью эмпирической формулы, используя результаты измерений. Вид эмпирической формулы в данном случае известен из физических соображений - при равнозамедленном движении те-
ла пройденное расстояние является квадратичной функцией времени
X = At2 + Bt + C
Легко установить, что С=0, поскольку Х=0 при t=0. Эмпирическая
формула принимает вид |
|
X = At + Bt |
(4) |
Для определения параметров А и В, нужно получить два уравнения. Воспользуемся методом средних и запишем уравнение (3) для всех точек (кроме начальной):
Запишем вместо этого уравнения путем его расщепления систему из двух уравнений
где = A*t1 + B*t1-Х1 и т.д.
Используя выражение (4) и табличные данные, получаем
(А*52+В*5-106)+(А*102+В*10-182)+(А*152+В*15-234)=0(А*202 + В*20 - 261)+(А*252 + В*25 - 275)=0
или окончательно
175*А + 15*В = 2611025*А + 45*В = 536
Решая эту систему уравнений, находим А = -0.30 и В = 39.07.
Следовательно, эмпирическую формулу (2.4), которая дает приближен-
2
ную связь между пройденным расстоянием и временем, можно записать в виде
X = -0.30*t + 39.07*t
Сравнивая это уравнение с уравнением X = a*t /2 + Vo*t,
2 |
|
получаем оценки для среднего ускорения тела и его начальной скорости: |
|
|
2 |
а = 2*А = -0.60 м/с; |
Vo = B = 39.07 м/c. |
Рассмотренные методы определения параметров эмпирической зависимости (метод выбранных точек и метод средних) являются сравнительно простыми, однако в ряде случаев получаемые с их помощью аппроксимации не обладают достаточной точностью.
12
Определение параметров эмпирической зависимости по методу наименьших квадратов
Запишем сумму квадратов отклонений для всех точек xo, x1,..., xn
n |
n |
2 |
|
xi , a0 , a1 ,...., am yi |
|||
S i2 |
|
||
i |
|
|
|
i 0 |
i 0 |
|
|
Параметры ao,a1,...,am |
эмпирической формулы (1) будем |
||
находить из условия минимума функции S = S( ao,a1,...,am ) В этом состоит метод наименьших квадратов.
С целью упрощения расчетов согласно данному методу рекоменду-
ется составить таблицу, где n - |
количество значений переменной х: |
||||
|
|
|
|
|
|
x0 |
x1 |
|
x2 |
y |
x · y |
1 |
x 1 |
|
. |
y 1 |
. |
. |
x 2 |
|
. |
y 2 |
. |
. |
. |
|
. |
. |
. |
. |
. |
|
. |
. |
. |
1 |
x n |
|
. |
y n |
. |
∑ x0 |
∑ x1 |
|
∑ x2 |
∑ y |
∑ x · y |
По итогам таблицы составляется система уравнений
∑x0· а + ∑ x1· b = ∑ y
∑x1· а + ∑ x2· b = ∑ x · y
решив которую, можно определить наилучшие значения параметров a и b.
Естественно ожидать, что значения найденной функции y=f(x,a,b)в точках х1, х2,...,хn будут отличаться от табличных значений y1, y2,...,yn. Для определения имеющихся отклонений необходимо составить таблицу (см. табл.2)
В соответствии с методом наименьших квадратов сумма квадратов полученных отклонений для заданного вида приближающей функции (и найденных значений параметров a и b) должна быть наименьшей
n |
|
i2 |
(5) |
i 0
Из двух разных приближений одной и той же табличной функции, следуя принципу наименьших квадратов, лучшим нужно считать то, для которого сумма (5) имеет наименьшее значение.
13
Исследовательская задача
Разработка любого электротехнического устройства сопровождается физическим или математическим моделированием. Физическое моделирование связано с большими материальными затратами, поскольку требуется изготовление макетов и их трудоемкое исследование. Часто физическое моделирование просто невозможно из-за чрезвычайной сложности устройства. В этом случае прибегают к математическому моделированию с использованием средств и методов вычислительной техники.
При выполнении рассматриваемого задания контрольных работ в качестве основного инструментального средства используется программа «Electronics Workbench», предназначенная для схемотехнического моделирования аналоговых и цифровых устройств различного назначения, которая отличается простым и легко осваиваемым пользовательским интерфейсом. Как учебная программа EWB обладает весьма важным достоинством, заключающимся в развитии творческого потенциала студентов, которые могут не только выполнять задания преподавателя, но и имеют возможность предложить и апробировать свои технические решения.
Моделирование и расчет электрических цепей постоянного тока с использованием законов Ома и Кирхгофа
Реальные источники электрической энергии можно представить двояко:
1)в виде генератора напряжения, который характеризуется ЭДС Е, численно равной напряжению холостого хода и включенной последовательно с внутренним сопротивлением Ri;
2)в виде генератора тока, который характеризуется током I, численно равным току короткого замыкания реального источника, и параллельно подключенной проводимостью GI.
Переход от генератора напряжения к эквивалентному генератору то-
ка осуществляется по формулам: |
|
|
I=E/Ri, |
Gi= 1/Ri, |
(6) |
а обратный переход от генератора тока к эквивалентному генератору напряжения по формулам:
E=I/Gi, Ri=l/Gi (7)
У идеального генератора напряжения внутреннее сопротивление Ri равно нулю, а у идеального генератора тока равна нулю внутренняя проводимость Gi.
Закон Ома
З-н Ома для участка цепи: ток в проводнике I равен отношению падения
14
напряжения U на участке цепи к ее электрическому сопротивлению R: Закон Ома иллюстрируется схемой на рис. 1, из которой видно, что
на участке цепи с сопротивлением R=100 Ом создается падение напряжения U=9,09В, измеряемое вольтметром. Согласно (8) ток в цепи I=9,09/100=90,9 мА, что и измеряет последовательно включенный в цепь амперметр. Отметим, что в рассматриваемой схеме внутреннее сопротивление амперметра выбрано равным 1012 Ом, т.е. очень малым, а входное сопротивление вольтметра -1012 Ом, т.е. очень большим, чтобы подключение измерительных приборов к цепи не оказывало скольконибудь заметного влияния на результаты измерений.
I=U/R |
(8) |
Рис. 1. Простейшая цепь постоянного тока
Отметим назначение ключей К и Х на рис. 1, управляемых нажатием одноименных клавиш клавиатуры. При размыкании ключа Х в схеме реализуется режим холостого хода, при этом вольтметр U измеряет ЭДС источника Е=10 В, а вольтметр Ui имеет нулевые показания. При замыкании ключа К в схеме реализуется режим короткого замыкания и, согласно (1), ток короткого замыкания Io=E/Ri=l0/10=l А. При этом вольтметр Ui измеряет падение напряжения Ui=Io·Ri=10 В.
З-н Ома для полной цепи: ток в замкнутой электрической цепи равен ЭДС источника Е, деленной на сопротивление всей цепи. Применительно к цепи на рис. 1 ее полное сопротивление равно Ri+R, и на основании закона Ома получаем I=E/(Ri+R)=90,9 мА, что и измеряет амперметр.
Отметим, что на дисплее амперметра отображаются только три старшие разряда измеряемой величины, так что результат измерения округляется.
Обобщенный з-н Ома: ток в замкнутой одноконтурной цепи равен отношению алгебраической суммы всех ЭДС к арифметической сумме всех сопротивлений. Перед расчетом выбирают направление обхода контура и считают это направление за положительное направление тока. При
15
определении алгебраической суммы ЭДС со знаком плюс берут те ЭДС, направления которых совпадают с выбранным положительным направлением тока, и со знаком минус — ЭДС с противоположными направлениями.
В качестве примера рассмотрим изображенную на рис. 2 одноконтурную цепь, состоящую из источников напряжения Е1=120 В, Е2=40 В и резисторов с сопротивлениями Rl=12 Ом и R2=8 Ом. Определим напряжение между точками А и В.
Выберем направление обхода контура по часовой стрелке. В таком случае ЭДС Е1 войдет со знаком "+", поскольку ток от Е1 совпадает с направлением обхода (положительным направлением тока во внешней цепи считается направление от положительного к отрицательному зажиму источника). При обходе же ветви с источником Е2 направление обхода не совпадает с направлением тока, который создается этим источником.
Рис. 2. Одноконтурная цепь с двумя источниками напряжения
Поэтому для схемы на рис. 38 ток в цепи
I=(Е1-E2)/(R1+R2)=80/20=4А.
Так как величина тока получилась положительной, то, следовательно, направление тока совпадает с выбранным. Если бы результат получился отрицательным, то это означает, что действительное направление тока в цепи противоположно выбранному.
Напряжение Uаь между точками А и В определяется с помощью закона Ома для участка цепи. Выберем участок А-Е2-В. Для этого участка закон Ома запишется в следующем виде:
I=(Uab-E2)/R2, откуда Uab=E2+I·R2=40+4·8=72 В.
Законы Кирхгофа Первый закон Кирхгофа: алгебраическая сумма всех токов, втекаю-
щих в любой узел, равна нулю. Токи, втекающие в узел, условно принимаются положительными, а вытекающие из него — отрицательными (или наоборот). Если, например, в узел втекает ток I1, а вытекают токи I2 и I3, то первый закон Кирхгофа может быть записан в виде выражения:
I1-I2-I3=0
16
Второй закон Кирхгофа: алгебраическая сумма ЭДС любого замкнутого контура равна алгебраической сумме падений напряжений на всех участках контура.
При применении второго закона Кирхгофа необходимо учитывать знаки ЭДС и выбранное направление токов на всех участках контура. Направление обхода контура выбирается произвольным; при записи левой части равенства ЭДС, направления которых совпадают с выбранным направлением обхода независимо от направления протекающего через них тока, принимаются положительными, а ЭДС обратного направления принимаются отрицательными. При записи правой части равенства со знаком плюс берутся падения напряжения на тех участках, в которых положительное направление тока совпадает с направлением обхода независимо от направления ЭДС на этих участках, и со знаком минус — на участках, в которых положительное направление тока противоположно направлению обхода.
Общая методика применения законов Кирхгофа для расчета сложных многоконтурных цепей такова. Устанавливается число неизвестных токов, которое равно числу ветвей р. Для каждой ветви задается положительное направление тока. Число независимых уравнений, составляемых по первому закону Кирхгофа, равно числу узлов q (точек соединения не менее чем трех проводников) минус единица, т.е. q-1. Число независимых уравнений, составляемых по второму закону Кирхгофа, равно числу контуров n=p-q+1. Общее число уравнений, составляемых по первому и второму законам Кирхгофа, равно числу неизвестных токов р. Решение этой системы уравнений и дает значения искомых токов.
Для иллюстрации изложенной методики рассмотрим многоконтурную цепь постоянного тока на рис. 3. В этой цепи всего три узла: А, В и С (q=3), следовательно, число независимых уравнений, составляемых по первому закону Кирхгофа, будет на единицу меньше, т.е. два. При числе ветвей цепи р=5 число контуров n=5-З+1=З, следовательно, по второму закону Кирхгофа можно составить три взаимно независимых уравнения. Таким образом, общее число независимых уравнений, составляемых по первому и второму законам Кирхгофа, будет равно числу неизвестных токов в пяти ветвях схемы.
17
Рис. 3. Многоконтурная цепь постоянного тока
Выберем положительные направления токов, которые на схеме обозначены соответствующим включением амперметров. Например, ток I1 течет справа налево и втекает в узел А (положительное направление тока), поскольку отрицательная клемма, отмеченная утолщенной черной линией, находится слева и ток через амперметр будет течь справа налево. Ток I2 вытекает из узла А, поскольку ток через одноименный амперметр будет течь сверху вниз (к отрицательному зажиму, расположенному на нижней грани иконки) и т.д.
Составим систему уравнений Кирхгофа:
для узла А |
I1-I2+I3-I5=0 |
для узла В |
-I1-I3-I14=0 |
для контура ABFA |
E1+E2=I1-R1-I3-R3 |
для контура АВСА |
E3=-I3-R3+I4-R4+I5-R5 |
для контура ADCA |
E2=I2-R2+I5-R5 |
После подстановки в полученные уравнения числовых значений они
приобретают следующий вид: |
I1-I2+I3-I5=0 |
|
I1-I3-I4=0 |
|
6-I1-10·I3=20 |
|
-10-I3+2,5·I4+15·I5=5 |
|
5·I2+15·I5=70 |
Решив полученную систему уравнений, имеем: I1=5А; I2=8А; I3=1А; I4= -6А; I5=2А, что соответствует показаниям приборов. Отрицательный знак для тока I4 означает, что истинное направление этого тока противоположно принятому.
18
Моделирование и расчет сложных электрических цепей постоянного тока методом контурных токов
Метод контурных токов (метод Максвелла) напоминает метод расчета с использованием законов Кирхгофа, однако он несколько проще и поэтому получил большее распространение на практике при расчетах многоконтурных цепей, состоящих из n независимых контуров. Определение токов в ветвях сводится к решению системы n=p-q+1 уравнений для контурных токов I1k, I2k, I3k …; действительный же ток в каждой ветви находится как алгебраическая сумма контурных токов, протекающих через соответствующую ветвь. Выбор направлений контурных токов произволен. Каждая из ветвей сложной электрической цепи должна войти хотя бы в один из анализируемых контуров.
Для иллюстрации рассматриваемого метода расчета обратимся к схеме на рис. 4, на которой выбранные направления контурных токов обозначены стрелками, токи в ветвях контуров — расположением входных зажимов амперметров.
Составим систему уравнений для контуров:
El-E2-E3= I1k (Rl+R2)- I2k R2
E2-E4= I2k (R2+R5+R4)- I3k R4- I1k R2 -E3-E4= I3k (R6+R4)+ I2k R4
После подстановки числовых значений получим:
60=20 I1k –10 I2k
24= -10 I1k +22 I2k +7 I3k -16=7 I2k +22 I3k
Решив эту систему уравнений, найдем контурные токи:
I1k = 5 А, I2k = 4 A, I3k= -2 А.
Теперь найдем истинные токи во всех ветвях. В ветви, где действует ЭДС Е1, истинный ток I1 имеет направление контурного тока I1k и равен ему: I1= I1k= 5А. В ветви с сопротивлением R5 истинный ток I5 имеет направление контурного тока I2k и равен ему: I5 = I2k = 4А. В ветви с сопротивлением R6 истинный ток I6 имеет направление, противоположное контурному току I3k, и равен I6= - I3k = 2 А. В ветви с сопротивлением R2 истинный ток I2 получится наложением контурных токов I1k и I2k и будет иметь направление большего из них, т.е. I2= I1k - I2k =1А. В ветви с сопротивлением R4 истинный ток I4 получится наложением контурных токов I2k и I3k и будет иметь направление большего из них, т.е. I4= I2k + I3k =4+(-2)=2 А. В ветви, где действует ЭДС ЕЗ, истинный ток I3 получится наложением контурных токов I1k и I3k и будет иметь направление тока I1k, т.е. I3= I1k + I3k =5+(-
19
2)=3А. Из сравнения полученных расчетных данных с показаниями приборов на рис. 4 видно, что они полностью совпадают.
.
Рис. 4. Схема сложной трехконтурной цепи
Моделирование и расчет сложных электрических цепей постоянного тока методом наложения
При методе наложения ток в любой ветви схемы рассчитывается как алгебраическая сумма токов, вызываемых в ней каждой ЭДС в отдельности. При этом предполагается, что при анализе воздействия одной ЭДС остальные ЭДС принимаются равными нулю и при расчетах учитываются только их внутренние сопротивления.
При использовании метода наложения потребуются некоторые навыки в преобразовании электрических цепей. Во всех случаях замена одних схем другими, им эквивалентными, не должна приводить к изменению токов и напряжений на участках цепи, не подвергшихся преобразованию. Наиболее распространенные и часто употребляемые преобразования цепей сводятся к следующим.
1. Замена последовательных сопротивлений одним эквивалентным. Эквивалентное сопротивление R цепи, состоящей из последовательно соединенных сопротивлений R1, R2,…,.Rn, равно сумме этих сопротивле-
ний: R= R1+ R2+…+Rn (9) 2. Замена параллельных сопротивлений одним эквивалентным. Со-
противления параллельны, если все они присоединены к одной паре узлов. Эквивалентное сопротивление R цепи, состоящей из ряда параллельно соединенных сопротивлений, определяется из формулы суммирования проводимостей:
l/R=l/ R1+l/ R2+...+l/ Rn |
(10) |
В случае двух параллельно соединенных |
сопротивлений |
R=R1·R2/(R1+R2). |
|
20 |
|
