Ирина Эланс

Автор который поможет с любыми образовательными и учебными заданиями

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

БЕЛОРУССКИЙ НАЦИОНАЛЬНЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ

Кафедра «Теория механизмов и машин»

Определение параметров поступательного

движения тела на плоскости

КУРСОВАЯ РАБОТА

по курсу «Информатика»

Исполнитель Вихновский А.С.

Группа 303117

Руководитель Луцко Н.Я.

2009 г.

СОДЕРЖАНИЕ

Введение

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

1. Постановка задачи

Тело массой m, на которое действуют движущая сила F_Д=F_Д(S) и сила сопротивления F_С, разгоняется на участке пути S_Р. После этого действие движущей силы прекращается ( сила F_С продолжает действовать), начинается торможение, в процессе которого тело пройдет до остановки расстояние S_Т за счет накопленной при разгоне кинетической энергии. F_Д=F₀+c^s+1, m = 3.0 кг, F_С=15 Н, F₀=20 H, с=0,4, S_Р=1.0 м, N=10.

Требуется:

определить зависимости от пути S скорости v(s), ускорения a(s), времени t(s);
установить время T_р прохождения телом участка S_P и время Т_т прохождения участка S_Т;
по полученным данным построить графики v(s), a(s), t(s) для интервала перемещения [0, S_P+ S_Т].

Расчетная схема приведена на рис. 1.1.

Рис. 1.1. Расчетная схема для определения параметров

движения при поступательном движении.

2. Математическая модель движения

Анализ поступательного движения тела показывает, что исходными данными для определения параметров движения (перемещения, скорости, ускорения, времени) являются масса m, движущая сила F_Д, сила F_С, а также начальные значения параметров движения. При использовании дискретной модели задачи весь путь разбивается на некоторое количество элементарных участков длиной ΔS=S_i-S_i-1 (рис. 2.1).

Рис. 2.1. Расчетная схема

На каждом интервале связь кинематических параметров описывается теоремой об изменении кинетической энергии:

(2.1)

Откуда можно выразить скорость движения:

(2.2)

При определении времени Δt прохождения участка ΔS будем считать скорость движения постоянной, равной средней скорости в пределах участка:

откуда

(2.3)

Аналогично, предполагая, что ускорение a_i на участке ΔS постоянно, имеем

(2.4)

Применим построенную математическую модель к расчету параметров поступательного движения тела на участке разгона [0, S_p] и на участке торможения [S_p, S_p+S_T] (рис. 2.2).

Рис. 2.2.

Разобьем каждый из участков движения на n равных элементарных участков длиной и соответственно. Полученные промежуточные положения тела пронумеруем от 1 до 2n+1. Переменная i определяет номер промежуточного положения тела. К участку разгона относятся положения с номерами от 1 до n+1.

Начальные параметры движения в положении i=1 считаются известными и равными S₁=0, v₁=0, t₁=0. Начальное ускорение а₁ определяется из закона Ньютона , который в нашем случае при i=1 примет вид

где F_Д(S₁)=F₀+c^S_¹⁺¹.

Для остальных положений тела при i=2,...,n+1 параметры движения определяются в соответствии с математической моделью по формулам

(2.5)

(2.6)

(2.7)

(2.8)

(2.9)

Интеграл в формуле (2.6) содержит аналитически заданную подынтегральную функцию f(s) = F_Д(S)-F_c с переменной интегрирования S. Он может быть вычислен:

точно — с использованием первообразной по формуле Ньютона-Лейбница;
приближенно - по методу трапеций.

Для нахождения интеграла используем метод трапеций:

(2.10)

Расчет параметров движения на участке торможения требует предварительного определения его длины S_Т. При этом исходим из условия, что вся накопленная при разгоне кинетическая энергия расходуется на преодоление силы сопротивления Fc, совершающей работу А_с=F_c*S_T, т.е.

откуда

(2.11)

Начальные параметры для участка торможения, соответствующие положению i=n+1, частично являются известными. Так, из процесса разгона получены S_n+1, v_n+1, t_n+1. При переходе к торможению имеет место разрыв функции ускорения. Новое значение ускорения, соответствующее началу участка торможения, равно .

Параметры движения в промежуточных положениях участка торможения при i=n+2,…,2n+1 определяются следующим образом:

;

Быстродействие на участке разгона будет равно T_P=t_n+1, а на участке торможения – T_t=t_2n+1-t_n+1.

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

Исходные данные (ввод): m, F₀, F_C, C, S_P, n.
Для первого положения S₁=0, v₁=0, t₁=0, F_Д(S₁)=F₀+c^S_¹⁺¹, a₁=(F_д(S₁)- F_c)/m
Для остальных положений при i=2,…,n+1

S_i=S_i-1+ΔS_p

int вычисляется по формуле трапеций

4.3.

4.4. ;

4.5. ;

4.6. .

5. Вывод параметров движения для участка разгона при i=1,…,n+1

5.1. Вывод i, S_i, v_i, a_i, t_i

6. Вывод быстродействия для участка разгона T_p=t_n+1

Для участка торможения

9. .

10. Для i=n+2, …,2n+1:

10.1. ;

10.2. ;

10.3. ;

10.4. ;

10.5. .

11. Вывод параметров движения для торможения при i=n+1, …, 2n+1

11.1. Вывод i, S_i, v_i, a_i, t_i

12. Вывод быстродействия для участка торможения Т_Т= t_2n+1-t_n+1.

4.Схема алгоритма решения задачи

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

Математическое

обозначение

F_с

S_р

F₀

ΔS_р

ΔS_Т

Т_р

Т_т

Идентификатор

dSp

dSt

Математическое

обозначение

S_Т

int

V_ср

Идентификатор

int

vsr

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

Program kurs_17 {Вихновский А.С., группа 303117}

{Определение параметров поступательного движения}

{тела на плоскости}

{Вариант 27}

uses crt;

type mas=array[1..300] of real;

Var

v,t,S,a,fdS:mas;

m,v0,dt,Sp,St,dSp,dSt,c:real;

int,tt,tp,fc,F0,vsr,k:real;

i,n:integer;

fu:text;

begin

clrscr;

assign(fu,'kurs_17.rez');

rewrite(fu);

writeln(fu,' ':20, 'Определение параметров поступательного');

writeln(fu,' ':25, 'движения тела на плоскости');

writeln(fu);

writeln(fu,' ':25,'‚Вихновский А.С., группа 303117');

writeln(fu);

writeln(fu,' ':34,'Вариант 27');

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

write('Масса m=');

readln(m);

write('Начальный момент сил f0=');

readln(f0);

write('Сила сопротивления fc=');

readln(fc);

write('Путь участка разгона Sp=');

readln(Sp);

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

readln(n);

write('С=');

readln(c);

writeln(fu);

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

writeln(fu,' ':20,'Масса m=',m:5:2,' кг');

writeln(fu,' ':20,'Начальный момент сил f0=',f0:5:2,' Н');

writeln(fu,' ':20,'Сила сопротивления fc=',fc:5:2,' Н');

writeln(fu,' ':20,'Путь участка разгона Sp=',Sp:5:2,' м');

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

writeln(fu,' ':20,'C=',C:5:2);

clrscr;

dSp:=Sp/n;

S[1]:=0;

v[1]:=0;

t[1]:=0;

fds[1]:=f0+exp(ln(c)*(s[1]+1));

a[1]:=(fds[1]-fc)/m;

for i:=2 to n+1 do

begin

S[i]:=S[i-1]+dSp;

fds[i]:=f0+exp(ln(c)*(s[i]+1));

int:=(fds[i]+fds[i-1]-2*fc)*dSp/2;

v[i]:=sqrt((2/m)*((m*v[i-1]*v[i-1]/2)+int));

vsr:=(v[i]+v[i-1])/2;

t[i]:=t[i-1]+(S[i]-S[i-1])/vsr;

a[i]:=(v[i]-v[i-1])/(t[i]-t[i-1]);

end;

for i:=1 to 80 do write (fu,'-');

writeln(fu);

writeln(fu,' ':23,'I',' ':9,'I',' ':9,'I',' ':9,'I');

writeln(fu,' ':21,'i I S I v I a I t ');

writeln(fu,' ':23,'I',' ':9,'I',' ':9,'I',' ':9,'I');

for i:=1 to 80 do write (fu,'-');

writeln(fu);

for i:=1 to n+1 do

writeln(fu,' ':20,i:2,' I ',S[i]:7:5,' I ',v[i]:7:5,' I ',a[i]:7:5,' I ',t[i]:7:3);

tp:=t[n+1];

writeln(fu);

writeln(fu,' ':20,'Быстродействие на участке разгона Tp=',tp:5:2);

writeln(fu);

St:=m*sqr(v[n+1])/(2*fc);

a[n+1]:=-fc/m;

dSt:=St/n;

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

begin

S[i]:=S[i-1]+dSt;

v[i]:=sqrt(2*(m*sqr(v[i-1])/2-fc*dSt)/m);

vsr:=(v[i]+v[i-1])/2;

t[i]:=t[i-1]+dSt/vsr;

a[i]:=(v[i]-v[i-1])/(t[i]-t[i-1]);

end;

for i:=1 to 80 do write (fu,'-');

writeln(fu);

writeln(fu,' ':23,'I',' ':9,'I',' ':9,'I',' ':10,'I');

writeln(fu,' ':21,'i I S I v I a I t ');

writeln(fu,' ':23,'I',' ':9,'I',' ':9,'I',' ':10,'I');

for i:=1 to 80 do write (fu,'-');

writeln(fu);

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

writeln(fu,' ':20,i:2,' I ',S[i]:7:5,' I ',v[i]:7:5,' I ',a[i]:7:5,' I ',t[i]:7:3);

writeln(fu);

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

writeln(fu,' ':20,'Быстродействие на участке торможения Тт =',tt:5:2);

close(fu);

writeln('Rabota zakonchena');

repeat until keypressed

end.

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

Определение параметров поступательного

движения тела на плоскости

Вихновский А.С., группа 303117

Вариант 27

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

Масса m= 3.00 кг

Начальный момент сил f0=20.00 Н

Сила сопротивления fc=15.00 Н

Путь участка разгона Sp= 1.00 м

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

C= 0.40

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

I I I I

i I S I v I a I t

I I I I

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

1 I 0.00000 I 0.00000 I 1.80000 I 0.000

2 I 0.10000 I 0.59903 I 1.79416 I 0.334

3 I 0.20000 I 0.84583 I 1.78300 I 0.472

4 I 0.30000 I 1.03441 I 1.77281 I 0.579

5 I 0.40000 I 1.19277 I 1.76352 I 0.668

6 I 0.50000 I 1.33181 I 1.75504 I 0.748

7 I 0.60000 I 1.45711 I 1.74730 I 0.819

8 I 0.70000 I 1.57201 I 1.74024 I 0.885

9 I 0.80000 I 1.67868 I 1.73380 I 0.947

10 I 0.90000 I 1.77864 I 1.72792 I 1.005

11 I 1.00000 I 1.87298 I 1.72256 I 1.060

Быстродействие на участке разгона Tp= 1.06

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

I I I I

i I S I v I a I t

I I I I

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

11 I 1.00000 I 1.87298 I -5.00000 I 1.060

12 I 1.03508 I 1.77687 I -5.00000 I 1.079

13 I 1.07016 I 1.67525 I -5.00000 I 1.099

14 I 1.10524 I 1.56705 I -5.00000 I 1.121

15 I 1.14032 I 1.45081 I -5.00000 I 1.144

16 I 1.17540 I 1.32440 I -5.00000 I 1.169

17 I 1.21048 I 1.18458 I -5.00000 I 1.197

18 I 1.24557 I 1.02588 I -5.00000 I 1.229

19 I 1.28065 I 0.83762 I -5.00000 I 1.267

20 I 1.31573 I 0.59229 I -5.00000 I 1.316

21 I 1.35081 I 0.00000 I -5.00000 I 1.434

Быстродействие на участке торможения Tт= 0.37

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

Рис. 8.1. Зависимость скорости от пути

Рис. 8.2. Зависимость времени от пути

Рис. 8.3. Зависимость ускорения от пути

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

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

а) тело, движение которого мы исследовали, остановилось при S≈1,35 м, значение скорости в этой точке равно 0;

б) в точке перехода от ускоренного движения к замедленному S=1 м, скорость в этой точке максимальна v≈1, 87 м/с (рис.8.1.);

в) ускорение для пути разгона положительно, а для пути торможения – постоянно и отрицательно (рис. 8.3);

г) функция t(s) – монотонно возрастающая;

д) быстродействие для пути разгона равно 1,06 с, а для угла торможения – 0,37 с.

10. Литература

Офицеров Д. В., Старых В. А. Программирование в интегрированной среде Турбо-Паскаль: Справ. Пособие. – Мн.: Беларусь, 1992.
Луцко Н. Я. Информатика. Контрольные работы и курсовое проектирование: Учебно-методическое пособие для студентов-заочников машиностроительных специальностей. – Мн.: БГПА, 2001.
Поляков Д.Б., Круглов И.Ю. Программирование в среде Турбо-Паскаль: Версия 5.5. – М.: Издательство МИА, АО «Росвузнаука», 1992г.
Быкадоров Ю.А., Кузнецов А.Т., Шербаф А.И. «Информатика и вычислительная математика 10-11»: Учебное пособие. – Мн., Беларусь, 2000г.
Фаронов В. В. Турбо Паскаль 7.0. Начальный курс. Учебное пособие. – М.: «Нолидж», 1998.

Определение параметров поступательного движения тела на плоскости 📙 Курсовая → 🆔 149709