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

Математическое обозначение I0 M0 Mc n Δφp φ t φp φt
Идентификатор I0 M0 Mc n dfp fi t fip fit
Математическое обозначение Δφt i ε ωcp int Tp Tt Md1 Md C
Идентификатор dft i b wcp int Tp Tt Md1 Md C

6. Текст программы

program kurs; {Курсовая работа студента Лабоцкого Д.В.}

{Исследование вращательного движения вала Вариант 13}

uses crt;

type Big=array[1..30] of real;

var Md,fi,w,t,b,int:Big;

n,i :integer;

fe:text;

C,Mc,I0,Wcp,fip,dfp,fit,dft,Tp,Tt,M0:real;

begin clrscr;

assign(fe,'kurs-13v.rez');rewrite(fe);writeln(fe);

writeln(fe,' ':15,'Определение параметров вращательного движения',

' тела');

writeln(fe);

writeln(fe, ' ':40,' Лабоцкий Д.В.');

writeln(fe);

writeln(fe, ' ':30,'Вариант 13');

writeln('Введите исходные данные');

write('Момент инерции тела равен I0= ');readln(I0);

write('Коэффициент для движущего момента равен М0= ');readln(M0);

write('Момент сопротивления равен Мc= ');readln(Mc);

write('Угол разгона fip= ');readln(fip);

write('Количество интервалов разбиения n= ');readln(n);

writeln('Исходные данные занесены в файл результатов');

writeln(fe);

writeln(fe,' ':25,'Исходные данные');

writeln(fe);

writeln(fe,' ':10,'Момент инерции тела равен I0= ',I0:5:2,' кг/м2');

writeln(fe,' ':10,'Коэффициент движущего момента М0= ',M0:5:2,' нм');

writeln(fe,' ':10,'Момент сопротивления Мc= ',Mc:5:2,' нм');

writeln(fe,' ':10,'Угол разгона fip= ',fip:5:2,' рад');

writeln(fe,' ':10,'Количество интервалов разбиения n= ',n:2);

dfp:=fip/n;

fi[1]:=0;W[1]:=0;t[1]:=0;

Md[1]:=M0+ln(fi[1]+1)+sqrt(fi[1]);

b[1]:=(Md[1]-Mc)/I0;

for i:=2 to (n+1) do begin

fi[i]:=fi[i-1]+dfp;

Md[i]:=M0+ln(fi[i]+1)+sqrt(fi[i]);

int[i]:=(((Md[i]-Mc)+(Md[i-1]-Mc))*dfp)/2;

W[i]:=sqrt((2/I0)*(I0*sqr(W[i-1])/2+int[i]));

Wcp:=(W[i]+W[i-1])/2;

t[i]:=t[i-1]+(fi[i]-fi[i-1])/Wcp;

b[i]:=(W[i]-W[i-1])/(t[i]-t[i-1]) end;

Tp:=t[n+1];

write(fe,' ':10);

for i:=1 to 57 do

write(fe,'_');writeln(fe);

writeln(fe,' ':10,'I',' ':10,'I',' ':10,'I',' ':16,'I',' ':16,'I');

writeln(fe,' ':8,' I fi[i] I W[i]  I b[i] I ',

' t[i] I ');

writeln(fe,' ':10,'I',' ':10,'I',' ':10,'I',' ':16,'I',' ':16,'I');

write(fe,' ':10);

for i:=1 to 57 do write(fe,'-');writeln(fe);

for i:=1 to n+1 do

writeln(fe,' ':7,i:2,' I',' ',fi[i]:7:3,' I',' ',W[i]:7:3,' I',' ',

b[i]:7:3,' I',' ',t[i]:7:3,' I');

writeln('Произведен расчет параметров разгона');

fit:=(I0*sqr(W[n+1]))/(2*Mc);

b[n+1]:=-Mc/I0;

dft:=fit/n;

for i:=n+2 to (2*n+1) do begin

fi[i]:=fi[i-1]+dft;

W[i]:=sqrt((2/I0)*((I0*sqr(W[i-1])/2)-(Mc*(fi[i]-fi[i-1]))));

Wcp:=(W[i]+W[i-1])/2;

t[i]:=t[i-1]+(fi[i]-fi[i-1])/Wcp;

b[i]:=(W[i]-W[i-1])/(t[i]-t[i-1]) end;

for i:=n+1 to (2*n+1) do

writeln(fe,' ':7,i:2,' I',' ',fi[i]:7:3,' I',' ',W[i]:7:3,' I',' ',

b[i]:7:3,' I',' ',t[i]:7:3,' I');

writeln('Произведен расчет параметров торможения');

write(fe,' ');

for i:=1 to 60 do

write(fe,'_');

writeln(fe);

writeln(fe);

Tt:=t[2*n+1]-t[n+1];

writeln(fe,' Быстродействие для угла разгона равно Tp= ',

Tp:7:3,' сек');

writeln(fe,' Быстродействие для угла торможения равно Tt= ',

Tt:7:3,' сек');

close(fe);

writeln('Результаты вычислений занесены в файл kurs-13v.rez');

repeat until keypressed

end.


7. Распечатка результатов.

Определение параметров вращательного движения тела

Лабоцкий Д.В.

Вариант 13

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

Момент инерции тела равен I0= 2.50 кг/м2

Коэффициент движущего момента М0= 15.50 нм

Момент сопротивления Мc= 10.00 нм

Угол разгона fip= 0.20 рад

Количество интервалов разбиения n= 10

I I I I I

I fi[i] I W[i] I b[i] I t[i] I

I I I I I

---------------------------------------------------------

1 I 0.000 I 0.000 I 2.200 I 0.000 I

2 I 0.020 I 0.299 I 2.232 I 0.134 I

3 I 0.040 I 0.425 I 2.280 I 0.189 I

4 I 0.060 I 0.522 I 2.308 I 0.231 I

5 I 0.080 I 0.605 I 2.333 I 0.267 I

6 I 0.100 I 0.678 I 2.354 I 0.298 I

7 I 0.120 I 0.745 I 2.374 I 0.326 I

8 I  0.140 I 0.807 I 2.393 I 0.352 I

9 I 0.160 I 0.865 I 2.411 I 0.376 I

10 I 0.180 I 0.919 I 2.428 I 0.398 I

11 I 0.200 I 0.971 I 2.444 I 0.419 I

11 I 0.200 I 0.971 I -4.000 I 0.419 I

12 I 0.212 I 0.921 I -4.000 I 0.432 I

13 I 0.224 I 0.868 I -4.000 I 0.445 I

14 I 0.235 I 0.812 I -4.000  I 0.459 I

15 I 0.247 I 0.752 I -4.000 I 0.474 I

16 I 0.259 I 0.686 I -4.000 I 0.490 I

17 I 0.271 I 0.614 I -4.000 I 0.509 I

18 I 0.282 I 0.532 I -4.000 I 0.529 I

19 I 0.294 I 0.434 I -4.000 I 0.554 I

20 I 0.306 I 0.307 I -4.000 I 0.585 I

21 I 0.318 I 0.000 I -4.000 I  0.662 I

Быстродействие для угла разгона равно Tp= 0.419 сек

Быстродействие для угла торможения равно Tt= 0.243 сек

8. Графическое представление результатов



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

Анализ результатов показывает:

а) В начальный момент времени скорость равна нулю, тело начинает двигаться с начальным ускорением β=2,2 м/с2

б) При перемещении тела его скорость увеличивается, так как значение силы, действующей на тело, увеличивается при перемещении, а ускорение уменьшается.

в) После того, как движущая сила перестала действовать, тело начало двигаться по инерции с постоянным ускорением, а скорость за счет силы трения уменьшается до остановки тела.

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


Литература

1.      Офицеров Д.В., Старых В.А. Программирование в интегрированной среде Турбо-Паскаль: Справ. пособие.—Мн.: Беларусь, 1992.

2.      Петров А.В. и др. Вычислительная техника и программирование: Курсовая работа/ А.В. Петров, М.А. Титов, П.Н. Шкатов; Под ред. А.В. Петрова.—М.: Высш. школа, 1992.

3.      Поляков Д.Б., Круглов И.Ю. Программирование в среде Турбо-Паскаль: Версия 5.5.—М.: Изд-во МИА, А/О Росвузнаука, 1992.

4.      Фигурнов В.Э. IBM PC для пользователя: Краткий курс.—Сокращенная версия 7-го издания.—М.:ИНФРА, 1999.

5.      Н.Я. Луцко.,П.П. Анципорович., Информатика Контрольные работы и курсовое проектирование: Учебно-методическое пособие для студентов-заочников машиностроительных специальностей


Информация о работе «Моделирование движения на плоскости»
Раздел: Информатика, программирование
Количество знаков с пробелами: 12594
Количество таблиц: 2
Количество изображений: 9

Похожие работы

Скачать
29980
3
2

... порыв ветра, интенсивность горения пороха, влажность и т.д. Из всех случайных величин рассмотрим только сопротивление воздуха. Учитывая то, что сопротивление воздуха направленно в противоположную сторону движения тела [4], то силу сопротивления можно представить в виде двух составляющих: вертикальной и горизонтальной, при этом, чем выше скорость, тем больше сопротивление и наоборот. Выделим ...

Скачать
74075
20
1

... проведения коррекции траектории МКА, моделирования процесса, и расчет потребного топлива для проведения коррекции траектории. 3) Исследование динамики системы коррекции траектории при стабилизации углового положения в процессе проведения коррек­ции траектории МКА. 2.4. МОДЕЛИРОВАНИЕ ДВИЖЕНИЯ ЦЕНТРА МАСС МКА 2.4.1.УРАВНЕНИЕ ДВИЖЕНИЯ КА Рассмотрим невозмущенное движение материальных точек М и m в ...

Скачать
156021
0
0

... или снега, в следах пробуксовки отсутствует поверхностный слой, покрывающий дорогу. О направлении движения буксующего автомобиля в этих случаях можно судить по направлению движения частиц дорожного полотна, выбрасываемых буксующим колесом. Наиболее ярко это явление наблюдается при буксовке колеса на снегу, песке, на сырой грунтовой дороге. Выбрасываемые буксующим колесом частицы располагаются по ...

Скачать
36871
3
34

... диаметрах критического сечения представлены на рисунке 2.24 Рисунок 2.24 - Зависимость оптимальной высоты поднятия фурмы от давления при различных диаметрах критического сечения сопла Лаваля 3. Численное исследование движения жидкости Приведены уравнения Навье - Стокса установившегося осесимметричного движения несжимаемой вязкой жидкости в переменных функция тока - вихрь. Проведено ...

0 комментариев


Наверх