Исследование удлинения стержня конической формы

Исследование удлинения стержня конической формы

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

Стальной cтержень конической формы круглого поперечного сечения (рис. 1) наглухо закреплен в конце O и подвергается действию продольной силы N, приложенной к концу стержня на расстоянии L от места закрепления. Наименьший диаметр равен 2r₁, наибольший - 2r₂.

Исследовать удлинение стального и медного стержней на различных расстояниях z и для различных радиусов r₁ и r₂. Построить графики зависимости ΔL(z).

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

Длина стержня L=5 м

Модуль упругости для стали E=2000000 MПa

Сила N=50 H

Радиус r₁=0,1 м

Радиус r₂=0,2 м, r₂ =0,3 м, r₂ =0,4 м, r₂ =1 м.

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

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

Радиус сечения стержня на расстоянии z от левого конца равен

Площадь сечения стержней на расстоянии z равна

Удлинение бруса находится по формуле

, где

Е - модуль упругости материала

Вычисление интеграла методом трапеций.

Интеграл оценивается вычислением суммы площадей элементарных трапеций со сторонами, равными значениям f(x) в начале и конце элементарного отрезка. Это приближение равносильно замене функции отрезком прямой, соединяющей значения f(x) в начальной и конечной точках отрезка (рис. 2).

Рис. 2. Метод трапеции

Площадь каждого элементарного сегмента разбиения считается по формуле

где

Тогда площадь искомой фигуры будем искать по формуле:

Следовательно, формула трапеций для численного интегрирования имеет вид:

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

1. Вводим исходные данные

l, N, m, r₁, r₂, r₂₂, r₂₃, r₂₄, E;

2. Выводим исходные данные

l, N, m, r₁, r₂, r₂₂, r₂₃, r₂₄, E;

Вычисляем удлинение стержня при r₂=0,2 м. с использованием процедуры TRAP

3.      TRAP (m, 0, l, E, r₂, Dl);

Вычисляем удлинение стержня при r₂=0,3 м. с использованием процедуры TRAP

4.      TRAP (m, 0, l, E, r₂₂, Dl);

Вычисляем удлинение стержня при r₂=0,4 м. с использованием процедуры TRAP

5.      TRAP (m, 0, l, E, r₂₃, Dl);

Вычисляем удлинение стержня при r₂=1 м. с использованием процедуры TRAP

6.      TRAP (m, 0, l, E, r₂₄, Dl);

Алгоритм функции f

1.

Алгоритм процедуры TRAP

.

. Для i=1, m+1

2.1

.2

3. Для i=2, m+1

;

. Для i=1, m+1

Вывод(Int_i).

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

Схема головной программы

Функция f

Процедура TRAP

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

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

kyrs16;crt;Vect=array [1..100] of real;=function (z, r1, r, l:real):real;

{$F+} function f (z, r1, r, l:real):real;:=pi*sqr (r1+(r-r1)*z/l);

{$F-}m:integer; n, E, r1, r2, r22, r23, r24, l:real; f1, f2:text; dl:vect;TRAP (m:integer; xn, xk, E:real; r:real; var int:vect);i:integer; h:real; x, y:vect;:=(xk-xn)/m;i:=1 to m+1 do begin[i]:=xn+(i-1)*h;[i]:=1/f (x[i], r1, r, l);;i:=2 to m+1 do[i]:=int [i-1]+h/2*N/e*(y[i]+y [i-1]);i:=1 to m+1 do writeln (f2, x[i]:5:2,' ', int[i]:10:8);;;(f1,'dan16.txt');(f1);(f2,'res16.res');(f2);(f1, l, n, m, r1, r2, r22, r23, r24, E);(f2,' Kyrsovoi proekt');(f2,' Issledovanie ydlineniya sterznei koni4eskoi formi');(f2,' Isxodnie dannie');(f2,'l=', l:2:0,' N=', N:3:0,' m=', m:2,' r1=', r1:3:1,' r2=', r2:3:1);(f2,' r22=', r22:3:1,' r23=', r23:3:1,' r24=', r24:3:1,' E=', E:6:0);(f2,'Rezultati raboti');(f2,'Udlinenie sterznya pri r2=0.2 m.');(f2,'x, m dl, mm');(m, 0, l, E, r2, dl);(f2,'Udlinenie sterznya pri r2=0.3 m.');(f2,'x, m dl, mm');(m, 0, l, E, r22, dl);(f2,'Udlinenie sterznya pri r2=0.4 m.');(f2,'x, m dl, mm');(m, 0, l, E, r23, dl);(f2,'Udlinenie sterznya pri r2=1 m.');(f2,'x, m dl, mm');(m, 0, l, E, r24, dl);('PA6OTA 3ABEPIIIEHA');(f1); close(f2); until keypressed;

end.

7. Результаты работы программы

Kyrsovoi proektydlineniya sterznei koni4eskoi formidannie= 5 N= 50 m=20 r1=0.1 r2=0.2=0.3 r23=0.4 r24=1.0 E=2000000rabotisterznya pri r2=0.2m., m dl, mm

.00 0.00000000

.25 0.00018970

.50 0.00036213

.75 0.00051955

.00 0.00066384

.25 0.00079658

.50 0.00091910

.75 0.00103254

.00 0.00113787

.25 0.00123594

.50 0.00132746

.75 0.00141307

.00 0.00149333

.25 0.00156872

.50 0.00163968

.75 0.00170658

.00 0.00176976

.25 0.00182953

.50 0.00188614

.75 0.00193986

.00 0.00199089sterznya pri r2=0.3m., m dl, mm

.00 0.00000000

.25 0.00018168

.50 0.00033297

.75 0.00046090

.00 0.00057051

.25 0.00066547

.50 0.00074854

.75 0.00082181

.00 0.00088694

.25 0.00094519

.50 0.00099761

.75 0.00104504

.00 0.00108815

.25 0.00112750

.50 0.00116357

.75 0.00119676

.00 0.00122739

.25 0.00125575

.50 0.00128208

.75 0.00130660

.00 0.00132948sterznya pri r2=0.4m., m dl, mm

.00 0.00000000

.25 0.00017469

.50 0.00030876

.75 0.00041493

.00 0.00050110

.25 0.00057244

.50 0.00063247

.75 0.00068369

.00 0.00072792

.25 0.00076648

.50 0.00080041

.75 0.00083049

.00 0.00085734

.25 0.00088146

.50 0.00090324

.00 0.00094103

.25 0.00095753

.50 0.00097269

.75 0.00098666

.00 0.00099959sterznya pri r2=1m., m dl, mm

.00 0.00000000

.25 0.00014678

.50 0.00022165

.75 0.00026722

.00 0.00029792

.25 0.00032002

.50 0.00033670

.75 0.00034975

.00 0.00036022

.25 0.00036882

.50 0.00037601

.75 0.00038211

.00 0.00038735

.25 0.00039190

.50 0.00039588

.75 0.00039941

.00 0.00040254

.25 0.00040535

.50 0.00040788

.75 0.00041017

.00 0.00041226

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

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

Удлинение стержня напрямую зависит от его радиуса. В пример этому служат результаты вычислений удлинения стержня в зависимости от радиуса, описанные выше.

Литература

1.   Рапаков Г.Г., РжеуцкаяС.Ю. Тurbo Pascal для студентов и школьников. - СПБ.: БХВ - Петербург, 2004. - 352 с.:ил.

2.      Анципорович П.П., Алейникова О.И., Булгак Т.И., Луцко Н.Я. Информатика. Учебно-метод. Пособие к лабораторным работам для студ. машиностроит. спец. В 4 ч. - Мн.: БНТУ, 2009.

.        П.И. Рудаков, М.А. Федотов. Основы языка Pascal учебный курс 2-ое издание исправленное. Москва. «Радио и связь», «Горячая линия - Телеком». 2000.

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