Задачи динамики, решаемые в среде системы Паскаль 


Мы поможем в написании ваших работ!



ЗНАЕТЕ ЛИ ВЫ?

Задачи динамики, решаемые в среде системы Паскаль



 

Теорема о движении центра масс системы материальных точек.

В случае сохранения скорости центра масс. Теорема о движении центра масс системы материальных точек. Зависимость между скоростью центра масс и скоростями точек материальной системы имеет вид:

Vс=(∑ mk *Vk)/M, т.е.

x c =(∑mk*xk)/M, yc=(∑mk*yk)/M, zc=(∑mk*zk)/M (1*)

Здесь

V с =xc*i+yc*j+zc*k

Зависимость между ускорением центра масс и с ускорениями точек материальной системы выражается соотношением

W с =(∑mk*Wk)/M, т.е.

xc=(∑mk*xk)/M, yc=(∑mk*yk)/M, zc=(∑mk*zk)/M, (2*)

Здесь

W с =xc*i+yc*j+zc*k

Напомним формулировку теоремы о движении центра масс: центр масс системы материальных точек движется как материальная точка, масса которой равна массе материальной системы и к которой приложены все внешние силы, действующие на систему: M*Wc=∑Fk.

Та же теорема, записанная в проекциях на оси декартовых координат имеет вид:

Mxc=∑Fkx, Myc=∑Fky, Mzc=∑Fkz, (3*)

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

Задачи динамики поступательного движения твердого тела решаются посредством теоремы о движении центра масс системы материальных точек.

Действительно, применив эту теорему, мы определим уравнение траектории, скорость и ускорение центра тяжести твердого тела. При поступательном движении твердого тела траектории всех точек одинаковы, а скорости и ускорения их соответственно равны.

С помощью теоремы о движении центра масс можно решать прямые и обратные задачи динамики. Последовательность решения задач:

1. Изобразить на рисунке все внешние силы системы;

2. Выбрать систему осей координат;

3. Записать теорему о движении центра масс (3*) в проекциях на декартовы оси координат;

4. Вычислить суммы проекций всех внешних сил системы на оси декартовых координат и подставить их в (3*);

5. В зависимости от условия решать прямую, либо обратную задачи динамики.

В некоторых прямых задачах бывают заданы все внешние силы, корме одной, массы всех материальных точек системы и законы их движения. Тогда после выполнения первых четырех пунктов для вычисления левых частей уравнения (3*) надо воспользоваться вспомогательными формулами

Mxc=∑mkxk, Myc=∑mkyk, Mzc=∑mkzk (4*)

где mk – масса k-ой точки, а xk, yk, zk – уравнения ее движения, ввести эти результаты в (3*) и определить неизвестную силу.

В некоторых обратных задачах бывают заданы все внешние силы, массы всех точек системы и законы движения всех точек, кроме одной и требуется определить движение этой точки. Тогда, после выполнения первых четырех пунктов также следует воспользоваться формулами (4*), полученные результаты ввести в левые части уравнений (3*) и затем найти искомый закон движения точки.

Если в состав системы входят тела с непрерывным распределением масс, то следует записать координаты xk, yk, zk центров тяжести этих тел и затем воспользоваться формулами (4*).

Задача. Эпициклический механизм, расположенный в вертикальной плоскости, установлен на горизонтальной идеально гладкой плоскости и прикреплен к ней болтами K и L. Зубчатое колесо 1 радиуса r1 неподвижно. С2 – центр тяжести зубчатого колеса 2 весом Р2 и радиусом r2. С1 – центр тяжести станины А и колеса 1, общий вес которых равен Р1. Массой кривошипа С1С2, вращающегося с постоянной угловой скоростью w, пренебречь. В начальный момент кривошип занимал правое горизонтальное положение. Определить:

1) нормальное давление механизма на плоскость,

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

3) наибольшее горизонтальное усилие, действующее на болты,

движение центра тяжести С1 станины механизма после среза болтов K и L.

Рис. 4. Геометрическая модель наклонной плоскости.

 

Решение. Материальная система состоит из двух масс: неподвижного колеса 1 со станиной и подвижного колеса 2. Изобразим внешние силы этой системы: Р1 – вес станины и неподвижного колеса 1, Р2 – вес подвижного колеса 2, R y – суммарная нормальная реакция плоскости, R x – суммарная тангенциальная реакция болтов K и L. Направим ось Oy по вертикали через точку С1, ось x – вдоль горизонтальной плоскости направо.

Запишем теорему о движении центра масс системы в проекциях на оси x и y:

Mxc=∑Fkx, Myc=∑Fky, Mzc=∑Fkz

В данной задаче

∑ Fkx = Rx, ∑ Fky = Ry - P 1 - P 2, Rx = Mxc, (1)

Ry= Myc+P1+P2 (2)

Для определения сил Rx и Ry остается подсчитать Mxc и Myc. Вычисление Mxc и Myc ведется по формулам:

Mxc=∑mkxk, Myc=∑mkyk.

В данном случае

Mxc= m1 x1+m2 x2 и Myc= m1 y1+m2 y2, (3).

Где x1 и y1 координаты центра тяжести С1 станины механизма и неподвижного колеса 1, x2 и y2 – координаты центра тяжести С2 подвижного колеса 2.

Как видно из рис., x1 =0, y1 =ОС1 – постоянная, x1=C1 C2 cosw t=(r1+r2) cos w t (угол поворота кривошипа С1С2 равен φ=wt, так как по условию w постоянна), y2=ОС11С2 sinw t=ОС1+(r1+r2) sinw t.

Вычислив вторые производные x1, y1, x2, y2 по времени t находим x1 =0 y1 =0, x2=-(r1+r2) w2 cosw t, y2=-(r1+r2) w2 sinw t.

Внеся эти значения в формулы (3), получим:

Mxc = -m2 (r1+ r2 )w 2соs wt,             (4)

Myc = -m2(r1+ r2 )w 2 sin wt                (5)
После подстановки (4) в (1) и (5) в (2) находим:

Rx = -P2 /g *(r1+ r2 )w2 со s wt           (6)

Ry= P1+ P2 - P2/g *(r1+ r2 )w 2 sin wt (7)

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

N y = P1+ P2 -P2 /g *(r1+r2 ) w 2 sin wt

Наибольшее давление:

Ny max = P 1 + P 2 + P 2 / g * (r 1 + r 2 ) w 2

Наименьшее давление:

Ny min = Р1 + P 2 - P 2 / g * (r 1 + r 2 ) w 2

В условиях отсутствия болтов механизм может начать подпрыгивать над горизонтальной плоскостью. Это будет иметь место при Rymin <0, т.е при Р1 +P2-P2/g* (r1 + r2) w2 <0, откуда следует, что угловая скорость w вращения кривошипа C 1 C 2, при которой происходит подпрыгивание механизма, должна удовлетворять неравенству

w > √ g*(P1+P2) / P2(r1+r2).

Горизонтальное давление, действующее на болты, направлено противоположно Rх (см. формулу (6)), причем

Nx=P2/g*(r1 + r2)w2 coswt.

Наибольшее давление равно

Nxmax=P2/g*(r1 + r2)w2

Допустим, что под действием, силы Nx произошел срез болтов.
Тогда весь механизм начнет двигаться по идеально гладкой горизонтальной плоскости.

На рис. б изображен механизм в положении, когда точка С1 сместилась с оси у направо на х 1. Так как станина механизма находится в движении относительно оси х, то х 1является функцией времени t.

Из чертежа видно, что в данном случае

х 21 + С1С 2 cos wt = х1 + (r 1 + r 2) cos wt.

Следовательно,

M xc1 х 12 x 2 = (m 1 +m 2) x 1 – m 2 (r 1 + r 2) w 2 cos wt (8)

Теорема о движении центра масс системы материальных точек в проекции на ось х имеет вид

Мхс = ∑ Fe kx

Так как после среза болтов реакция Rx отсутствует, а внешние силы Р1 Р 2 и Rу перпендикулярны к оси х, то ∑ F kx = 0 и Мхс = 0. Подставив в это уравнение значение Mx с из формулы (8), получим

1 + m 2) х1 - m 2 (r 1 + r 2) w 2 cos wt = 0,

т. е.                                                 

x 1 = Р 2 /(Р 1 + Р 2)* (r1 + r2) w2 cos wt, (9)

Это - дифференциальное уравнение движения центра тяжести С1 станины механизма по идеально гладкой горизонтальной плоскости при отсутствии болтов. Для интегрирования уравнения (9) должны быть известны начальные условия движения точки С1. Так как в момент среза болтов точка C1 находилась на оси у и была в покое, то начальные условия движения записываются в виде:

при t = 0 x1 =0 и y1 = 0.

Проинтегрировав дифференциальное уравнение (9), получим:

x 1 = Р 2 / Р 1 + Р 2 * (r1 + r2) w sin wt + D1

После подстановки начального условия движения t = 0 и x 1 = 0 имеет D 1 = 0, т. е

x 1 = Р 2 / Р 1 + Р 2 * (r1 + r2) w sin wt

Вторично проинтегрировав, находим х1 = - Р212 *(г1 + r2) cos wt +D2. Использовав то, что при t=0, х1=0, имеем:

D 2 = Р 212 * (r1 + r2)

т.е. x 1 = Р 2 / Р 12 * (r 1 + r 2 )(1- cos wt).

Итак, центр тяжести С1 станины механизма в случае отсутствия болтов совершает гармонические колебания с амплитудой Р 212 * (r 1 + r 2) и круговой частотой, равной угловой скорости w вращения кривошипа С1С2.

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

Mwe =∑ Fk+∑Jrk,

где F k внешние силы, a Jrk силы инерции в относительном движении.

В проекциях на оси декартовых координат имеем:

Мхе =∑ Fkxe + ∑ Jrkx Муе = ∑ Fkye + ∑ Jrky ,

k =1

Мzе = ∑ Fkze + ∑ Jrkz

k =1

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

Изобразив все внешние силы системы Р 1, Р2, Rx и R y (см. рис. в), добавляем центробежную силу инерции в относительном движении

Jrn = - Р2 /g* w rn. Так как точка С2 в относительном движении описывает окружность с центром С1 радиуса С1С2 = r 1 + r2, то, центро­стремительное ускорение wrn, направлено от С2 к С1 и, следовательно, центробежная сила инерции в относительном движении Jr n направлена противоположно. По модулю

Jr n = - Р2 /g*wrn= Р2 /g*(r 1 + r2)w2

Вращательная сила инерции в относительном движении J = -Р2 /g*wr τ равна нулю, так как кривошип вращается равномерно. Применив дифференциальные уравнения переносного поступательного движения материальной системы в проекциях на оси х и у:

Мхе =∑ Fkxe + ∑ Jrkx, Муе = ∑ Fkye + ∑ Jrky ,

k =1 k =1 k =1 k =1

получим

Mx e = Rx + Jrn coswt, Mye = Re — P 1 — P 2 + Jrn sinwt,

Так как х e = х1 ,ye=y1 , Jrn =P2/g*(r1+r2) w 2, то

Мх 1 = Rx+P2/g(r1+r2)w2coswt, (10)

My1=Ry-P1- Р 2 +P2/g (r1 + r2) w2 sinwt.                   (11)

В случае механизма, закрепленного болтами, центр тяжести С1 колеса 1 и станины неподвижен, т. е. х11= 0, и дифференциальные уравнения принимают вид

Rx+P2/g(r1+r2)w2coswt =0, (12)

Ry- -P1- Рг +P2/g (r1 + r2) w2 sinwt, (13)

откуда вытекает, что проекция нормальной реакции плоскости равна

Ry = P1 - Рг +P2 /g (r1 + r2) w2 sinwt.                       (14)

 Проекция на ось х горизонтальной силы реакции болтов равна

Rx= P2 / g (r1+r2 )w2coswt. (15)

Условие подпрыгивания определяем из (14), считая R у min отрицательным. Так как

Rymin = P 1 + Рг - P 2 / g *(r 1 + r 2) w 2, а Ry min <0, то

P 1 2 - P 2 / g *(r 1 + r 2) w 2 < 0

откуда w >√ g *(P 1 + P 2)/(P 2 (r 1 + r 2 ))

Для определения закона движения центра тяжести CL колеса 1 истанины механизма после среза болтов надо в формуле (10) положить Rx = 0. Тогда

Мх1 = P 2 / g *(r 1 + r 2) w 2 coswt,

Т.е. приходим к уравнению (9):

x1=P2 /(P1+ P2 )*(r1 + r2 ) w2cos wt,

решение которого было получено выше.

На основе разработанного алгоритма решения задачи по кинематике составим Паскаль – программу.

Program DINAMIKA;

Var

w,r1,r2,P1,P2,t,NxMax,Ny,x1:Real;

Const

g=9.8;

Begin

Writeln ('vvedite radius r1');

Readln (r1);

Writeln ('vvedite radius r2');

Readln (r2);

Writeln ('vvedite ves P1');

Readln (P1);

Writeln ('vvedite ves P2');

Readln (P2);

Writeln ('vvedite vremya');

Readln (t);

w:=sqrt((g*(P1+P2))/(P2*(r1+r2)));

Ny:=P1+P2-(P2/g)*(r1+r2)*w*w*cos(w)*t;

NxMax:=P2/g*(r1+r2)*w*w;

x1:=P2/P1+P2*(r1+r2)*(1-cos(w)*t);

Writeln ('w:=',w);

Writeln ('Ny:=',Ny:8:6);

Writeln ('NxMax:=',NxMax:8:6);

Writeln ('x1:=',x1:8:6);

Readln;

E nd.


ЗАКЛЮЧЕНИЕ

 

Целью курсовой работы являлась изучение полного спектра функциональных возможностей языка программирования Паскаль для решения задач прикладной механики.

Задачами данной работы являлись:

1. Освоение полного спектра функциональных возможностей языка программирования Паскаль;

2.  Постановка и решение задач прикладной механики традиционным способом;

3.  Решение задач механики в среде языка программирования Паскаль.

Методами работы при выполнении поставленных задач:

1. Теоретический анализ научно-технической литературы по языку программирования Паскаль;

2. Математическое моделирование задач прикладной механики;

3. Компьютерное решение задач прикладной механики.

На основе проведенного курсового исследования на тему «Приложения технологии языка программирования паскаль в прикладной механике» можно сформулировать следующие выводы:

1. Язык программирования высокого уровня Паскаль обладает широким спектром логических конструкций и функций, необходимых для успешного решения задач прикладной механики.

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


ЛИТЕРАТУРА

1. Бать М.И., Джанелидзе Г., Кельзон А.С. Теоретическая механика в примерах и задачах. Т.1. М.: Просвещение, 2000.

2. Бать М.И., Джанелидзе Г., Кельзон А.С. Теоретическая механика в примерах и задачах. Т.2. М.: Просвещение, 2000.

3. Бочкин А. И. Методика преподавания информатики. - Минск: Высшая школа, 1998.

4. Блашкин И.И., Буров А.А. Новые возможности Turbo-Pascal 6.0. — Спб.: Изд-во «Макет», 1992.

5. Бородич Ю.С. и др. Паскаль для персональных компьютеров: Справ. пособие/ Ю.С.Бородич, А.Н.Вальвачев, А.И.Кузьмич. — Мн.: Выш. шк.: БФ ГИТМП «НИКА», 1991.

6. Васильев П.П. Турбо Паскаль — мой друг: М.: Компьютер, ЮНИТИ, 1995.

7. Великов В.П., Новая информатика в школе // Информатика и образование. – 1986. - №1.

8. Вычислительная техника и программирование. Под редакцией А. В. Петрова М., Высшая школа, 1990.

9. Голубева О.В. Теоретическая механика. Изд-во «Высшая школа». М.: 1968.

10. Донцов Д.А. Самые нужные программы для Windows. Популярный самоучитель.- Спб.: Питер, 2006.

11. Джордейн Р. Справочник программиста персональных компьютеров типа IBM PC, XT, AT: Пер. с англ./ Предисл. Н.В.Гайского. — М.: Финансы и статистика, 1991.

12. Зозуля Ю. Компьютер на 100 % - Спб.: Питер, 2006.

13. Зуев Е.А. Язык программирования Turbo Pascal 6.0. — М.: Унитех, 1992.

14. Информатика. Базовый курс: Учеб. пособ. для студентов технических вузов / С.В. Симонович, Г. Евсеев, В. И. Мухаровский и др.; под ред. Симоновича – Спб.: Питер, 2005.

15.  Информатика: Учеб. пособ. для пед. спец. вузов / А.Р. Есаян, В.И. Ефимов, Л.П. Липецкая и др. - М.: Просвещение, 1991.

16. Лапчик М. П. Методика преподавания информатики. М.: Посвещение, 2001.

17. Левин А. Самоучитель полезных программ 3-е изд.- Спб.: Питер, 2003.Турбо Паскаль 7.0 - К.: Издательская группа BHV, 1998.

18. Марченко А. И., Марченко Л. И. Программирование в среде Turbo-Pascal 7.0-М., Бином Универсал, К.: Юниор, 1997.

19. Мизрохи А.М. Turbo Pascal и объектно-ориентированное программирование. — М.: Финансы и статистика, 1992.

20. Немнюгин С.А. Turbo Pascal. Программирование на языке высокого уровня. Учебник для вузов. 2-е изд.- Спб.: Питер, 2005.

21. Рывкин К.А. Справочник школьника по информатике. 7-11 кл. - М.: ООО Изд. дом «Оникс 21 век», 2005.

22. Справочник по процедурам и функциям Borland Pascal with Objects 7.0. — Киев: «Диалектика», 1993.

23. Фарафонов В.В. Турбо Паскаль 7.0. Начальный курс: учеб. пособие. - М.: Кнорус, 2006.

24. Фёдоров А. Особенности программирования на Borland Pascal. — Киев: «Диалектика», 1994.

25. Хершель Р. Турбо Паскаль/ 2-е изд., перераб. — Вологда: МП «МИК», 1991.

26. POWER TOOLS PLUS. Процедуры поддержки для Turbo Pascal 4.0.: Справочное руководство пользователя. Техническая документация.



Поделиться:


Последнее изменение этой страницы: 2020-03-14; просмотров: 216; Нарушение авторского права страницы; Мы поможем в написании вашей работы!

infopedia.su Все материалы представленные на сайте исключительно с целью ознакомления читателями и не преследуют коммерческих целей или нарушение авторских прав. Обратная связь - 44.197.114.92 (0.127 с.)