Двигатель из состояния покоя приводит в движение вал рабочей машины. Механическая характеристика двигателя состоит из двух прямолинейных участков: М_Д1=a₁+b₁щ при и М_Д2=a₂+b₂щ при . Приведенный к валу двигателя момент инерции вращающихся частей машины и двигателя J, момент сопротивления на том же валу равен М_С.

Определить зависимость скорости вала двигателя от времени щ(t). Вычислить скорость щ_уст и время начала t_уст установившегося движения вала. Построить графики зависимостей щ(t) и М_Д(t).

Исходные данные:

a1=96 Н*м

a₂=2100 Н*м

b₁=0,56 Н*м*с

b₁=20 Н*м*с

М_С=60 Н*м

J=12 кг*м²

n=20

t=0,1 с

е=0,1

2. Математическая модель решения задачи

В качестве математической модели задачи используем дифференциальное уравнение движения вала двигателя

Тогда на прямолинейном участке при задача Коши примет вид:

(1)

Значение щ_опр определим (см. рис. 1) из условия М_Д(щ_опр)=М_Д(щ_опр) или a₁+b₁щ_опр=a₂-b₂щ_опр.

Откуда

.

Решив задачу Коши на интервале аргумента методом Эйлера и Рунге-Кутта четвертого порядка точности, получим таблично заданную зависимость t(щ). Построенную таблицу можно рассматривать в качестве зависимости щ(t) и использовать для определения М_Д(t).

На прямолинейном участке при строим задачу Коши в виде

(2)

Для каждого значения аргумента t_i=t_i-1+t, начиная с i=n+2 и используя t_опр=t_n+1, щ_опр=щ_n+1, решаем задачу Коши (2) выбранным методом и получаем последовательность значений функции щ_n+2, щ_n+3, щ_n+4 и т.д. Процесс вычислений прекращаем при выполнении условия . Тогда в качестве щ_уст берем последнее вычисленное значение щ_i и t_уст=t_i.

Решение дифференциальных уравнений методом Эйлера.

Пусть задано дифференциальное уравнение первого порядка или y'=F (x, y)

На интервале [x₀, x_n] разобьём на n частей и получим X_0,…,X_n.

X_i=X₀+(i-1) h или X_i=X_i-1+h₁, где .

Соответствующее значения y₁=y*(x_i), где y*(x_i) - приближенное значение дифференциального уравнения.

Для получения численного решения дифференциального уравнения уравнение заменяется уравнениями относительно значений функции y*(x). Эти уравнения называют разностными. Простейшие разностные уравнения для заданного дифференциального уравнения имеют вид.

y_i+1-y_i=h*F (x, y)

y_i+1=y_i+h*F (x, y) - формула Эйлера.

Алгоритм метода Эйлера.

1) Ввод исходных данных(x₀, x_n,n, y₀).

2)

3) x[0]=x₀; y[0]=y₀.

4) Для i=1, n

4.1.) x[i]=x [i-1]+h;

4.2) y[i]=y [i-1]+h*f (x[i-1], y [i-1]);

5) Для i=0, n

5.1) Вывод x[i], y[i].

3. Алгоритм решения задачи

1. Вводим исходные данные

a1, a2, b1, b2, Mc, J, n, dt, epsilon;

2. wopr:=(a2-a1)/(b1+b2); tn:=0; wn:=0; wk:=wopr;

Решаем дифференциальное уравнение с использованием процедуры EILER

3. EILER (n, wn, wk, tn, fun1, w, t)

Записываем в файл RESULTS.RES полученную табличную зависимость w(t) на промежутке 0<=w<=wopr;

4. Для i=1…n+1

4.1. writeln (t[i], w[i], Md[i]);

Записываем в файл RESULTS.RES topr и wopr

5. writeln (topr, wopr);

6. i:=n+2;

Определяем функцию w(t) на промежутке w>=wopr с использованием цикла с постусловием

7. repeat

7.1. t[i]:=t [i-1]+dt;

7.2. w[i]:=w [i-1]+dt*fun2 (w[i-1]);

7.3. m:=i; i:=i+1;

until abs (w[i-1] - w[i])<=epsilon;

Записываем в файл RESULTS.RES полученную табличную зависимость w(t) на промежутке w>=wopr

8. Для i=n+2…m

writeln (t[i], w[i], Md[i]).

Записываем в файл результатов tyst, wyst, Mdyst.

9. writeln (tyst, wyst, Mdyst).

Алгоритм процедуры EILER

1. t[1]:=tn; t [n+1]:=tk; v[1]:=vn; h:=(t [n+1] - t[1])/n;

2. Для i=2…n+1

2.1 t[i]:=t[1] +(i-1)*h;

2.2 v[i]:=v [i-1]+h*f (t[i-1]);

Алгоритм функции fun1

1. fun1:=12/(2100+0.56*w-60)

Алгоритм функции fun2

1. fun2:=(2100-20*w-60)/12

4. Схема алгоритма

Функция fun1

Размещено на http://www.allbest.ru/

w:real

Функция fun2

Размещено на http://www.allbest.ru/

w:real

Процедура EILER

Размещено на http://www.allbest.ru/

n:integer; tn, tk, vn:real; f:fun;

Var t, v:vect

5. Таблица идентификаторов

вал двигатель скорость программа

3avisimost w(t) pri 0<=w<=wopr

t w Md

0.00 0.00 96.00

0.03 4.87 98.73

0.06 9.75 101.46

0.09 14.62 104.19

0.11 19.49 106.92

0.14 24.37 109.65

0.17 29.24 112.38

0.20 34.11 115.10

0.23 38.99 117.83

0.26 43.86 120.56

0.28 48.74 123.29

0.31 53.61 126.02

0.34 58.48 128.75

0.37 63.36 131.48

0.40 68.23 134.21

0.43 73.10 136.94

0.45 77.98 139.67

0.48 82.85 142.40

0.51 87.72 145.13

0.54 92.60 147.85

0.57 97.47 150.58

topr= 0.57s wopr=97.47rad/s Mdopr=150.58H*s

3avisimost w(t) pri w>wopr

t w Md

0.67 98.23 135.49

0.77 98.85 122.91

0.87 99.38 112.42

0.97 99.82 103.68

1.07 100.18 96.40

1.17 100.48 90.34

1.27 100.74 85.28

1.37 100.95 81.07

1.47 101.12 77.56

1.57 101.27 74.63

1.67 101.39 72.19

1.77 101.49 70.16

1.87 101.58 68.47

tyst=101.58s wyst= 1.87rad/s Mdyst=68.47H*s

8. Графическая часть

9. Анализ результатов

В результате работы программы была получена зависимость скорости вала двигателя от времени щ(t) и установившиеся время t_уст и скорость щ_уст.

За время t двигатель приобрел максимальную угловую скорость щ, которая в дальнейшей работе двигателя изменялась только в пределах заданной относительной погрешности измерений е.

Литература

1. Рапаков Г.Г., Ржеуцкая С.Ю. Тurbo Pascal для студентов и школьников. - СПБ.: БХВ - Петербург. 2004. - 352 с.

2. Анципорович П.П., Алейникова О.И., Булгак Т.И., Луцко Н.Я. Информатика. Учебно-метод. Пособие к лабораторным работам для студ. машиностроит. спец. В 4 ч. - Мн.: БНТУ, 2009.

Размещено на Allbest.ru