Исследование движения тела, брошенного под углом к горизонту, с учетом горизонтального сопротивления

Исследование движения тела, брошенного под углом к горизонту, с учетом горизонтального сопротивления

Содержание

. Введение

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

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

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

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

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

. Текст программы на языке Pascal

. Результаты решения задачи

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

. Приложение №1 (Microsoft Exel XP)

. Приложение №2 (MathCAD 2001)

. Литература

1. Введение

На практике часто встречается задача, когда тело брошено под углом к горизонту с начальной скоростью и с учётом сопротивления. Поэтому довольно актуально исследовать характер движения такого тела на одном типичном примере. В соответствие с математической моделью задачи для исследования характера движения тела необходимо найти решение задачи Коши для системы дифференциальных уравнений первого и второго порядка, решить её, используя метод Эйлера, и построить необходимые графики.

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

Тело массой брошено под углом к горизонту с начальной скоростью  (рис.1). На тело действует сила сопротивления , направленная горизонтально влево.

Исследовать характер движения тела, построив графики зависимостей υx(t), x(t), υy(t), y(t) и υ(t).

Значения исходных данных для тестирования: , , , , , .

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

Математическую модель задачи построим, используя второй закон Ньютона. Дифференциальные уравнения движения в проекциях на оси ОХ и ОY запишутся в виде:

 и .

Начальные условия движения тела известны из постановки задачи , , , . Таким образом, для исследования характера движения тела необходимо найти решение задачи Коши (сист. 1):

 сист. 1 сист. 2.

Преобразуем её к системе дифференциальных уравнений первого порядка с начальными условиями (сист. 2).

Для решения системы используем метод Эйлера. В первом положении в соответствии с постановкой задачи определим , , ,

, .

Обозначим Δt шаг изменения аргумента t. Тогда для i = 2, 3, …

, , ,

, , .

Процесс вычислений продолжаем до тех пор, пока тело не упадёт на землю, что соответствует выполнению условия . В результате получим таблично заданные зависимости υx(t), x(t), υy(t), y(t) и υ(t).

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

1.      начальные данные:

1.1 υ0=5 м/c

.2 Δt=0,025 с

.3 m=0,5 кг

.4 k=2,8 кг/с

.5 g=9,8 м/с2

.6 α=30o

2.      данные в первой точке пути:

2.1 i=1

.2 t1=0

.3 υx1=υ0 cosα

2.4 x1=0

2.5 υy1=υ0 sinα

.6 y1=0

3.      пока  повторять:

3.1 i=i+1

3.2

3.3

3.4

3.5

3.6

3.7

. n = i-1

. i =1..n:

.1. вывод ti, Vxi, xi, Vyi, yi, Vi:

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

Max

Vybor

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

. Текст программы на языке Pascal

движение тело программа коши

uses graph,crt;v0=5; dt=0.025; alpha=pi/6; k=2.8; g=9.8; m=0.5;tab=array[1..40,1..6] of real;rez:text;,i:byte;:tab;

{====================================================}MyGraphInit;graphdriver,GraphMode,ErrorCode:integer;:=detect;(graphdriver,GraphMode,'');:=graphresult;Errorcode<>grOk then('proizoshla oshibka grafiki: ',grapherrormsg(errorCode));('rabota programmy prervana');;;;

{====================================================}osi_koordinat(A:tab; n,z,k:byte);i:byte;,dy,x:real;:string;:=A[n,1]/10;:=A[n-k,z]/10;(0,1,1);(100,20,100,320);(100,320,600,320);i:=0 to 10 do(99,20+30*i,101,20+30*i);z=4 then x:=A[n-k,z]-i*dy*2x:=A[n-k,z]-i*dy;(x:4:3,s);(50,15+30*i);(s);(100+50*i,319,100+50*i,321);:=i*dt;(x:4:3,s);(90+50*i,324,s);;I:=1 to 10 do;(1,1,1);(8);(103,20+30*(i-1), 600,20+30*(i-1));(100+50*i,318,100+50*i,20);(0,1,1);;;

{====================================================}MAX(const alpha,dt,k,m,g,v0:real; var A:tab; var n:byte);i:byte;:=1;[i,1]:=0;[i,2]:=v0*cos(alpha);[i,3]:=0;[i,4]:=v0*sin(alpha);[i,5]:=0;[i,6]:=v0;A[i,5]>=0 do:=i+1;[i,1]:=A[i-1,1]+dt;[i,2]:=A[i-1,2]-dt*k/m*A[i-1,2];[i,3]:=A[i-1,3]+dt*A[i-1,2];[i,4]:=A[i-1,4]-dt*g;[i,5]:=A[i-1,5]+dt*A[i-1,4];[i,6]:=sqrt(sqr(A[i,2])+sqr(A[i,4]));;:=i-1;;

{=====================================================}graph_Vx(A:tab; n:byte);i,z,k:byte;;:=2;:=n-1;_koordinat(a,n,z,k);(red);i:=1 to n do(round(100+900*A[i,1]),round(320-69*A[i,2]),round(100+900*A[i+1,1]),round(320-69*A[i+1,2]));(100,400,'gfrafik zavisimosti Vx(t)');;;

{=====================================================}graph_X(A:tab; n:byte);i,z,k:byte;;:=3;:=0;_koordinat(a,n,z,k);(100,320);(red);i:=1 to n do(round(100+A[i,1]*500/A[n,1]),round(320-A[i,3]*300/A[n,3]));(100,400,'gfrafik zavisimosti X(t)');;;

{====================================================}graph_Vy(A:tab; n:byte);i,z,k:byte;;:=4;:=n-1;_koordinat(a,n,z,k);(red);(round(100+A[i,1]*500/A[n,1]),round(320-(A[i,4]-A[n,4])*300/(A[1,4]-A[n,1])));i:=1 to n do(round(100+A[i,1]*500/A[n,1]),round(320-(A[i,4]-A[n,4])*300/(A[1,4]-A[n,4])));(100,400,'grafik zavisimosti Vy(t)');;;

{=====================================================}graph_V(A:tab; n:byte);i,z,k:byte;;:=6;:=n-1;_koordinat(a,n,z,k);(100,20);(red);i:=1 to n do(round(100+A[i,1]*500/A[n,1]),round(20+(A[i,6]-A[1,6])*140/(A[n,6]-A[1,6])));(100,400,'grafik zavisimosti V(t)');;;

{=====================================================}graph_Y(A:tab; n:byte);i,z,k:byte;;:=5;:=0;_koordinat(a,n,z,k);(100,320);(red);i:=1 to n do(round(100+A[i,1]*500/A[n,1]),round(320-A[i,5]*22.5/A[n,5]));(100,400,'grafik zavisimosti Y(t) ');;;

{=====================================================}vybor(A:tab; n:byte);z:byte;:=0;;;('1-grafik Vx(t)');('2-grafik X(t)');('3-grafik Vy(t)');('4-grafik Y(t)');('5-grafik V(t)');('sdelaite vybor nazhatiem');('sootvetstvyuschei knopki');(z);;(z<>1) and (z<>2) and (z<>3) and (z<>4) and (z<>5) thencase z of

: graph_Vx(A,n);

: graph_x(A,n);

: graph_Vy(A,n);

: graph_y(A,n);

: graph_V(A,n);;;;;

{=====================================================}(alpha,dt,k,m,g,v0,b,n);(b,n);(rez,'rez.txt');(rez);(rez,' raspe4atka rezul''tatov');(rez,'issledovanie dvizheniya tela brozhenogo pod uglom k gorizontu');(rez,'s u4etom gorizontal''nogo soprotivleniya ');(rez);(rez,'KOZEL A.N. gr.103814 ');(rez);(rez,'variant_30');(rez);(rez,'ishodnye dannye:');(rez,'massa tela m=',m:4:2,' kg');(rez,'na4al''naya skorost'' tela Vo=',v0:2,' m/s');(rez,'koeficient soprotivleniya sredy K=',k:3:2,' kg/s');(rez,'uskorenie Sv.Padeniya g=',g:3:2,' m/s^2');(rez,'promezhutok vremeni dt=',dt:6:4,' s');(rez,'ugol broska alpha=',alpha*180/pi:2:0,' grad');(rez);writeln(rez);(rez,'__________________________________________________');(rez,'| i t Vx x Vy y V|');(rez,'|_________________________________________________|');i:=1 to n do(rez,'|',i:3,' | ',B[i,1]:5:3,' | ',B[i,2]:5:3,' | ',B[i,3]:5:3,' | ',B[i,4]:7:3,' | ',B[i,5]:5:3,' | ',[i,6]:4:3,'|');(rez,'------------------------------------------------------');(rez);

end.

. Результаты решения задачи

atka rezul'tatovdvizheniya tela brozhenogo pod uglom k gorizontu

s u4etom gorizontal'nogo soprotivleniyaA.N. gr.103814_30dannye:tela m=0.50 kgal'naya skorost' tela Vo= 5 m/ssoprotivleniya sredy K=2.80 kg/sSv.Padeniya g=9.80 m/s^2vremeni dt=0.0250 sbroska alpha=30 grad

______________________________________________________

| I        t      Vx      x          Vy        y         V|

|_____________________________________________________|

| 1 | 0.000 | 4.330 | 0.000 | 2.500 | 0.000 | 5.000|

| 2 | 0.025 | 3.724 | 0.108 | 2.255 | 0.063 | 4.353|

| 3 | 0.050 | 3.203 | 0.201 | 2.010 | 0.119 | 3.781|

| 4 | 0.075 | 2.754 | 0.281 | 1.765 | 0.169 | 3.271|

| 5 | 0.100 | 2.369 | 0.350 | 1.520 | 0.213 | 2.814|

| 6 | 0.125 | 2.037 | 0.409 | 1.275 | 0.251 | 2.403|

| 7 | 0.150 | 1.752 | 0.460 | 1.030 | 0.283 | 2.032|

| 8 | 0.175 | 1.507 | 0.504 | 0.785 | 0.309 | 1.699|

| 9 | 0.200 | 1.296 | 0.542 | 0.540 | 0.328 | 1.404|

| 10 | 0.225 | 1.114 | 0.574 | 0.295 | 0.342 | 1.153|

| 11 | 0.250 | 0.958 | 0.602 | 0.050 | 0.349 | 0.960|

| 12 | 0.275 | 0.824 | 0.626 | -0.195 | 0.351 | 0.847|

| 13 | 0.300 | 0.709 | 0.647 | -0.440 | 0.346 | 0.834|

| 14 | 0.325 | 0.610 | 0.664 | -0.685 | 0.335 | 0.917|

| 15 | 0.350 | 0.524 | 0.680 | -0.930 | 0.318 | 1.068|

| 16 | 0.375 | 0.451 | 0.693 | -1.175 | 0.294 | 1.259|

| 17 | 0.400 | 0.388 | 0.704 | -1.420 | 0.265 | 1.472|

| 18 | 0.425 | 0.333 | 0.714 | -1.665 | 0.229 | 1.698|

| 19 | 0.450 | 0.287 | 0.722 | -1.910 | 0.188 | 1.931|

| 20 | 0.475 | 0.247 | 0.729 | -2.155 | 0.140 | 2.169|

| 21 | 0.500 | 0.212 | 0.735 | -2.400 | 0.086 | 2.409|

| 22 | 0.525 | 0.182 | 0.741 | -2.645 | 0.026 | 2.651|

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

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

Тело пролетело расстояние 0,741 метра за время 0,525 секунды. Скорость тела уменьшалась, пока оно не достигло наивысшей точки своей траектории (0,351 метра). Пройдя данную точку, проекция скорости на ось y (υy) поменяла своё направление(значения стали отрицательными), и тело начало разгоняться. Разгон продолжался до соприкосновения с землей.

Следовательно, исходя из результатов, можно сделать вывод:

чтобы увеличить дальность и время полёта тела, необходимо уменьшить массу тела, подобрать оптимальный угол броска (45о) и место проведения эксперимента для уменьшения величины вертикального сопротивления (ветер, дождь и тд.), увеличить начальную скорость броска.

В данной задаче мы проводили расчеты через равные промежутки времени Δt=0,025 с. При этом в последней точке, в которой тело находилось максимально близко к земле, координата y составляла 0,026 метра. Следовательно, это даёт существенную погрешность. Чтобы уменьшить величину погрешности эксперимента, необходимо уменьшить значение Δt, например, до Δt =0,01 с.

10. Приложение №1 (Microsoft Exel XP)

11. Приложение №2 (MathCAD 2001)

Исходные данные для тестирования

Данные в первой точке пути

Графики

График зависимости Vx(t)

График зависимости x(t)

График зависимости Vy(t)

График зависимости y(t)

График зависимости V(t)

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

1) П.П.Анципорович, О.И.Алейникова, Т.И.Булгак, Н.Я.Луцко

Учебно-методическое пособие к лабораторным работам для студентов машиностроительных специальностей. Часть 1. Минск: БНТУ, 2009.

) Конспект лекций.