Приложения технологии языка программирования Паскаль в прикладной механике

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

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

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

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

Program Statika;

Var

x, y, a, Pmin, Pmax:Real;

R:Integer;

Begin

Writeln('vvedite ves sterchnya');

Readln(dv);

Writeln('vvedite dliny sterchnya');

Rreadln(dl);

Writeln('vvedite ugol');

Readln(?); {?=60}

Pmin:=(R*(((cos(?)*3.14/180)/cos(?)*3.14/180)+y)); {minimalnaya velichina gruza}

Pmax:=(R*(((cos(a)*3.14/180)/cos(?)*3.14/180)-y)); {maximalnaya velichina gruza}

{pri cos(?)<y sila xmax ne suschestvuet}

Begin

If xmax<0 then x:=xmin+xmax;

If xmax>0 then x:=Pmin-xmax;

End;

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

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

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

Readln;

End.

2.2. ЗАДАЧИ КИНЕМАТИКИ, РЕШАЕМЫЕ ПОСРЕДСТВОМ ЯЗЫКА ПАСКАЛЬ

Скорости точек плоской фигуры могут быть определены аналитическими , графическими или же графоаналитическими методами.

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

Хo₁=f ₁(t) , Yо₁ =f ₂(t) , ? = f₃ (t). (1*)

Проекции скорости точки М на неподвижные оси координат определяется равенствами:

V_x=Vо₁Х-W _z(Y-Yo₁) (2*)

V_y=Vо₁Y + W_z (Х-Хo₁) (3*)

В этих формулах V_x , V_y - искомые проекции скорости точки М на неподвижные оси координат; Vо₁x = Хo₁ , Vо₁y = Yo₁-проекции скорости полюса, начала подвижной системы координат ХY на неподвижные оси координат ; W_z =? -проекция угловой скорости на ось Z , перпендикулярную к плоскости движения ;ХУ координаты точки М в неподвижной системе координат ; Хо₁ ,Yо₁ - координаты полюса О₁ в неподвижной системе осей. Определение координат Х ,Y точки М , по заданным уравнениям движения плоской фигуры (1*) производится по формулам:

Х=Хо₁+Х₁ cos? - Y₁sin ?

Y= Yо₁+ Х₁ sin ? +Y₁ cos ?

Проекции скорости точки М на неподвижные оси координат находятся по формулам:

Vx₁ =Vo₁ x cos ? + Vo₁ у sin ? - Wzy₁ (4*)

Vу ₁= -Vo₁ x sin ? + Vo₁ у сos ? - Wz х₁ (5*)

В этих формулах Vx , Vу- искомые проекции скорости точки М на оси х, у подвижной системы координат , жестко связанной с плоской фигурой ; х у - координаты точки М в подвижной системе осей , остальные величины имеют то же значение , что и в уравнениях (2*) , (3*).

Величины скорости точки М по известным проекциям определяются формулой:

V= vVx ²+ Vу² =v Vx₁² + Vу₁² (6*)

Направляющие косинусы даются равенствами:

cos (V,X)= Vx/V , cos (V,У)= Vy/V (7*)

cos (V,X)= Vx/V , cos (V,У)= Vy/V (8*)

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

Рис. 2. Геометрическая модель плоской фигуры

V=Vo₁+W r , (9*)

в этой формуле

V - искомая скорости точки М;

Vo₁ - скорость полюса О,

W - угловая скорость плоской фигуры;

r - радиус-вектор , проведенный из полюса О в Рис. 3 точку М.

Таким образом, зная скорость какой-либо точки плоской фигуры, выбираем эту точку за полюс. Далее откладываем от точки М, скорость которой подлежит определению, вектор, равный скорости полюса, и вектор W r₁, направлен перпендикулярно к r₁, и равный по величине W r₁.Векторная сумма этих слагаемых и дает искомую скорость точки М, Если скорость точки М известна по направлению, то можно не знать величины вращательной скорости W r₁, так как и эта скорость, и искомая скорость точки М определяется пересечением прямой, совпадающей с направлением скорости точки М, и прямой, перпендикулярной к линии О М, проведенной из конца вектора Vо, отложенного из точки М.

После того, как треугольник или параллелограмм скоростей выражающий равенство (9*), построен, задача может считаться решенной. Величина и направление скорости точки М могут быть найдены по рис.(2) или получены из решения этого треугольника.

Формула распределения скоростей (9*)записывается часто в виде:

V_b=V_a+V_ba (9**)

где V-искомая скорость точки В, V_a-известная скорость точки А , избранной за полюс, а

V_ba= W АВ, V_ba = W AB, (9***)

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

Многие задачи могут быть решены при помощи теоремы о равенстве проекции скоростей концов отрезка плоской фигуры на направление отрезка.

Второй графоаналитический метод определения скоростей плоской фигуры основан на использовании мгновенного центра скоростей этой фигуры. При не поступательном движении плоской фигуры (W=0) в каждый данный момент существует точка тела , скорость которой равна нулю. Эта точка называется мгновенным центром скоростей и обычно обозначается через P. Единственным исключением является случай так называемого мгновенно поступательного движения (W=0), который будет рассмотрен отдельно. Выбирая мгновенный центр за полюс, имеем закон распределения скоростей в плоской фигуре.

V=W r ₁, V=W М Р , (10*)

где V-искомая, скорость произвольной точки М; W-угловая скорость плоской фигуры, r₁-радиус-вектор, проведенный из мгновенного центра скоростей в точку М, называемый мгновенным радиусом. Таким образом, скорости всех точек плоской фигуры являются в данный момент вращательными скоростями вокруг мгновенного центра скоростей. Их величина равна произведению величин угловой скорости на модуль мгновенного радиуса, а направлены они перпендикулярно к мгновенному радиусу. Таким образом, величина скоростей точки плоской фигуры пропорциональны величине мгновенных радиусов:

V_b / ВР=V_a / АР=...=W, (11*)

При графоаналитических методах решения задач рекомендуется такая последовательность действий:

1) Выбрать за полюс ту точку плоской фигуры, направление скорости которой известно;

2) найти другую точку плоской фигуры, направление скорости которой известно;

1) пользуясь формулой распределения скоростей, найти скорость этой точки плоской фигуры;

2) исходя из формулы распределения скоростей определить значение угловой скорости плоской фигуры в данный момент времени;

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

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

1. Выбрать за полюс точку плоской фигуры, скорость которой известна;

2. Пользуясь формулой распределения скоростей, построить скорость другой точки плоской фигуры, у которой известно направление ее скорости;

3. Спроектировать полученный треугольник скоростей на направление прямой, соединяющей обе точки, и найти скорость второй точки;

4. Спроектировать треугольник скоростей на направление, перпендикулярное к прямой, соединяющей обе точки, и найти вращательную скорость второй точки по отношению к полюсу;

5. Разделив вращательную скорость на расстояние от точки до полюса, найти мгновенную угловую скорость плоской фигуры;

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

Если задача решается при помощи мгновенного центра скоростей, рекомендуется такая последовательность:

1) определить положение мгновенного центра плоской фигуры;

2) найти величину мгновенного радиуса той точки плоской фигуры, скорость которой известна и определить угловую скорость плоской фигуры, разделив величину скорости точки на величину мгновенного радиуса;

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

Теперь рассмотрим конкретную задачу по кинематике.

Задача. Для преобразования периодического возвратно-поступательного движения во вращательное движение в двигателях внутреннего сгорания применяют шатунно-кривошипный механизм, состоящий из кривошипа ОА=r, вращающегося вокруг неподвижной точки О в плоскости рисунка, шатуна АВ=l, шарнирно соединенного с кривошипом, и поршня В, двигающегося по горизонтальным направляющим цилиндра, совпадающим с прямой ОВ.

Рис. 3. Геометрическая модель плоского шатунно-кривошипного механизма

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

Program Kinematika;

Var t,f0,k,b1,b2,x1,x2,a,b,l,d,f1,f2,a1,a2,w1,w:real;

Begin

Write('vvedite welichini a b d l w t = ');

Readln(a,b,d,l,w,t);

x1:=(sqr(d+b)-a*a+l*l)/(2*l*(d+b));

f1:=(180/pi)*arctan(sqrt(1-sqr(x1))/x1);

x2:=(sqr(d+b)-a*a-l*l)/(2*l*a);

a1:=(180/pi)*arctan(sqrt(1-sqr(x2))/x2);

x1:=-(l*l+sqr(d-b)-a*a)/(2*l*(d-b));

f2:=(180/pi)*arctan(sqrt(1-sqr(x1))/x1);

x2:=(l*l-sqr(d-b)+a*a)/(2*l*a);

a2:=(180/pi)*arctan(sqrt(1-sqr(x2))/x2);

f0:=w*t;

k:=sqrt(l*l+b*b-2*l*b*cos(f0));

x1:=(k*k+d*d-a*a)/(2*k*d);

b2:=arctan(sqrt(1-sqr(x1))/x1);

x2:=(1/k)*sin(f0);

b1:=arctan(x2/sqrt(1-sqr(x2)));

w1:=(b*k*sin(b1+b2)*w)/(a*a*sin(b2));

Writeln('a1= ', a1:6:3,' a2= ',a2:6:3);

Writeln('f1= ', f1:6:3,' f2= ',f2:6:3);

Writeln('w1= ',w1:6:3);

Readln;

End.

2.3. ЗАДАЧИ ДИНАМИКИ, РЕШАЕМЫЕ В СРЕДЕ СИСТЕМЫ ПАСКАЛЬ

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

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

V_с=(?m_k *V_k)/M, т.е.

x_c=(?m_k*x_k)/M, y_c=(?m_k*y_k)/M, z_c=(?m_k*z_k)/M (1*)

Здесь

V_с=x_c*i+y_c*j+z_c*k

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

W_с=(?m_k*W_k)/M, т.е.

x_c=(?m_k*x_k)/M, y_c=(?m_k*y_k)/M, z_c=(?m_k*z_k)/M, (2*)

Здесь

W_с=x_c*i+y_c*j+z_c*k

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

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

Mx_c=?Fkx, My_c=?Fky, Mz_c=?Fkz, (3*)

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

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

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

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

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

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

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

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

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

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

Mx_c=?m_kx_k, My_c=?m_ky_k, Mz_c=?m_kz_k (4*)

где m_{k -} масса k-ой точки, а x_k, y_k, z_{k -} уравнения ее движения, ввести эти результаты в (3*) и определить неизвестную силу.

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

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

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

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

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

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

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

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

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

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

Mx_c=?F_kx, My_c=?F_ky, Mz_c=?F_kz

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

?F_kx=R_x, ?F_ky=R_y-P₁-P₂, R_x= Mx_c, (1)

R_y= My_c+P₁+P₂ (2)

Для определения сил R_x и R_y остается подсчитать Mx_c и My_c. Вычисление Mx_c и My_c ведется по формулам:

Mx_c=?m_kx_k, My_c=?m_ky_k.

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

Mx_c= m₁ x₁+m₂ x₂ и My_c= m₁ y₁+m₂ y₂, (3).

Где x₁ и y_{1 -} координаты центра тяжести С₁ станины механизма и неподвижного колеса 1, x₂ и y₂ - координаты центра тяжести С₂ подвижного колеса 2.

Как видно из рис., x₁=0, y₁=ОС₁ - постоянная, x₁=C₁ C₂ cosw t=(r₁+r₂) cos w t (угол поворота кривошипа С₁С₂ равен ?=wt, так как по условию w постоянна), y₂=ОС₁+С₁С₂ sinw t=ОС₁+(r₁+r₂) sinw t.

Вычислив вторые производные x₁, y₁, x₂, y₂ по времени t находим x₁=0 y₁=0, x₂=-(r₁+r₂) w² cosw t, y₂=-(r₁+r₂) w² sinw t.

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

Mx_c= -m₂ ( r₁+ r₂ )w² соs wt, (4)

My_c= -m₂( r₁+ r₂ )w² sin wt (5)
После подстановки (4) в (1) и (5) в (2) находим:

R_x = -P₂ /g *( r₁+ r₂ )w² соs wt (6)

R_y= P₁+ P₂ - P₂/g *( r₁+ r₂ )w² sin wt (7)

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

N_y=P₁+ P₂ -P₂ /g *( r₁+r₂ ) w² sin wt

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

N_{y max} = P₁ + P₂+ P₂/g * (r₁+ r₂ ) w²

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

N_{y min} = Р₁ + P₂ - P₂ /g * ( r₁ +r₂ ) w²

В условиях отсутствия болтов механизм может начать подпрыгивать над горизонтальной плоскостью. Это будет иметь место при R_ymin<0, т.е при Р₁ +P₂-P₂/g* (r₁ + r₂) w²<0, откуда следует, что угловая скорость w вращения кривошипа C₁C₂, при которой происходит подпрыгивание механизма, должна удовлетворять неравенству

w > vg*(P₁+P₂) / P₂(r₁+r₂).

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

N_x=P₂/g*(r₁ + r₂)w² coswt.

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

N_xmax=P₂/g*(r₁ + r₂)w²

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

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

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

х₂=х₁ + С₁С₂ cos wt= х₁ + (r₁ + r₂) cos wt.

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

Mx_c =т₁х₁ +т₂ x₂ = (m₁ +m₂)x₁ - m₂ (r₁ + r₂) w² cos wt (8)

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

Мх_с = ?F^e_kx

Так как после среза болтов реакция R_x отсутствует, а внешние силы Р₁ Р₂ и R_у перпендикулярны к оси х, то ?F_kx = 0 и Мх_с = 0. Подставив в это уравнение значение Mx_с из формулы (8), получим

(т₁ +m₂) х₁ -m₂ (r₁ + r₂) w² cos wt = 0,

т. е.

x₁ = Р₂/(Р₁+Р₂ )*(r₁ + r₂) w² cos wt, (9)

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

при t= 0 x₁ =0 и y₁ = 0.

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

x₁= Р₂/Р₁+Р₂ *(r₁ + r₂) w sin wt + D₁

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

x₁= Р₂/Р₁+Р₂ *(r₁ + r₂) w sin wt

Вторично проинтегрировав, находим х₁ = - Р₂/Р₁+Р₂ *(г₁ + r₂) cos wt +D₂. Использовав то, что при t=0, х₁=0, имеем:

D₂ = Р₂/Р₁+Р₂ *(r₁ + r₂)

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

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

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

Mw_e=?F_k+?J_rk ,

где F_k-- внешние силы, a J_rk -- силы инерции в относительном движении.

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

Мх_е =? F_kx^e+ ?J_rkx Му_е = ?F_ky^e + ?J_rky,

k=1

Мz_е = ?F_kz^e + ?J_rkz

k=1

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

Изобразив все внешние силы системы Р₁, Р_2, R_x и R_y (см. рис. в), добавляем центробежную силу инерции в относительном движении

J_rn = -Р₂ /g*w_rn. Так как точка С₂ в относительном движении описывает окружность с центром С₁ радиуса С₁С₂ = r₁+ r₂, то, центро-стремительное ускорение w_rn, направлено от С₂ к С₁ и, следовательно, центробежная сила инерции в относительном движении J_rn направлена противоположно. По модулю

J_rn = -Р₂ /g*w_rn= Р₂ /g*(r₁+ r₂)w²

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

Мх_е =? F_kx^e+ ?J_rkx , Му_е = ?F_ky^e + ?J_rky,

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

получим

Mx_e =R_x+J_rn coswt, My_e =R_e -- P₁-- P₂+J_rn sinwt,

Так как х_e = х₁ ,y_e=y₁ , J_rn =P₂/g*(r₁+r₂) w², то

Мх₁=R_x+P₂/g(r₁+r₂)w²coswt , (10)

My₁=R_y-P₁- Р₂ +P₂/g (r₁ + r₂) w² sinwt. (11)

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

R_x+P₂/g(r₁+r₂)w²coswt =0, (12)

R_y- -P₁- Р_г +P₂/g (r₁ + r₂) w² sinwt , (13)

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

R_y = P₁ - Р_г +P₂ /g (r₁ + r₂) w² sinwt. (14)

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

R_x= P₂ / g (r₁+r₂ )w²coswt. (15)

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

R_ymin = P₁ + Р_г - P₂ /g *(r₁ + r₂) w², а R_{y min}<0 , то

P₁ +Р₂ -P₂ /g *(r₁ + r₂) w²<0

откуда w>v g*(P₁+P₂)/(P₂(r₁+r₂ ))

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

Мх₁ = P₂/g*(r₁ + r₂) w² coswt ,

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

x₁=P₂ /(P₁+ P₂ )*(r₁ + r₂ ) w²cos 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;

End.

ЗАКЛЮЧЕНИЕ

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

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

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.: Справочное руководство пользователя. Техническая документация.