Решение системы дифференциальных уравнений (8) осуществляется методом Рунге-Кутта с использованием среды MATLAB через М-файл.
Соответствие машинных идентификаторов М-файла среды MATLAB и физических величин показано в таблице 9.3.
Таблица 3 – Соответствие машинных идентификаторов и физических величин
Физические величины | Машинные идентификаторы | Примечание |
t | V (1) | Начальные условия |
Р СТ | Р(С) | Потери в стали |
Р МПО, Р МЛО | Р(7), Р(8) | Потери в меди паза и лобовых частей соответственно |
Jхаi, = 2…4 | Vi(2),….. | Температура хладагента |
b | b | Температурный коэффициент меди |
G | G | Тепловая проводимость |
m | m | Масса |
Сi, = 1…3 | С i(1),…3 | Удельная теплоемкость |
Произведя соответствие физических величин с машинным идентификатором М-файла, записываем полученные уравнения в
М-файл, показанный на рисунке 9.4.
Рисунок 9.4 – М-файл расчета тепловых процессов в среде MATLAB
После записи М-файла, необходимо его запустить.
Для этого в главном окне программы MATLAB делаем ссылку на выполнение расчета тепловых потерь активных частей якоря генератора. Ссылка записывается следующим образом:
|
|
[t,y]=ode45(@anchor,[0,20],[20,20,20];
plot(t,y),grid,title('anchor'),legend('y'),
где [t,y] – переменные функции;
ode45 – методом Рунге-Кутта;
@anchor – название функции, которое совпадает с М-файлом;
[0,20] – начальные условия интервала решения систем уравнений;
[20,20,20] – начальные условия окружающей среды.