Рефетека.ру / Коммуникации и связь

Курсовая работа: Расчет жесткого стержня

Содержание


1. Задание

2. Схема нагруженного стержня

3. Исходные данные

4. Построение системы линейных алгебраических

5. Вывод формул проверки, достоверности вычисления опорных реакций

6. Вывод рабочих формул определение внутренних усилий стержня

7. Численный метод решения СЛАУ - метод Гаусса

8. Обоснование применения метода Гаусса

9. Блок - схема алгоритма

10. Программа

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

Литература

1. Задание


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

Вариант - 82-4г. Схема - 2.

Численный метод решения СЛАУ - метод Гаусса.


2. Схема нагруженного стержня


Расчет жесткого стержня

Расчет жесткого стержняP1, P2-сосредоточенная сила, Н

q4 - интенсивность распределенной нагрузки, H/м

C1, C2 - отрезок балки, м

L1, L2 - пролет балки, м

М1, M2 - круговой момент, HРасчет жесткого стержням

3. Исходные данные


P1=15kH P2=30kH L1=6м L2=12м

M1=10kHм M2=35kHм С1=3м C2=2м

L1=6м L2=12м q4=10kH

Расчет жесткого стержня

4. Построение системы линейных алгебраических


уравнений для определения опорных реакций.

Преобразуем исходную систему:

отбросим опорные стержни и заменим их опорными

реакциями (R1; R2; R3)

интенсивность распределённой нагрузки заменим эквивалентной

силой (F4 = q4Расчет жесткого стержняc2)

зададим систему координат.


Расчет жесткого стержняРасчет жесткого стержня


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


Расчет жесткого стержня:

Расчет жесткого стержня

Расчет жесткого стержня Расчет жесткого стержня

Расчет жесткого стержня

Расчет жесткого стержня


Представил уравнения равновесия балки в форме системы линейных алгебраических уравнений (СЛАУ).


Расчет жесткого стержня


Матричная форма записи СЛАУ вычисление опорных реакций балки

AРасчет жесткого стержняR=B

А - матрица коэффициентов при неизвестных

R - матрица неизвестных

В - матрица свободных членов


Расчет жесткого стержня Расчет жесткого стержня

Расчет жесткого стержня


5. Вывод формул проверки, достоверности вычисления опорных реакций


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


Y=R1-P1+R2=0

X=R3-P2-F4=0


6. Вывод рабочих формул определение внутренних усилий стержня


На рассматриваемом стержне выделим четыре участка длиной S (длина отрезка от начала до точки сечения стержня), для которых составим формулы для вычисления внутренних усилий: поперечной силы Q и изгибающего момента М.

s - отрезок от начала до точки сечения балки

I cечение Расчет жесткого стержня


Расчет жесткого стержня


II cечение Расчет жесткого стержня


Расчет жесткого стержня


Расчет жесткого стержня

III cечение Расчет жесткого стержня


Расчет жесткого стержня


Расчет жесткого стержня


IV cечение Расчет жесткого стержня


Расчет жесткого стержня


Расчет жесткого стержняРасчет жесткого стержня


В точках границ Расчет жесткого стержня, Расчет жесткого стержня,Расчет жесткого стержняорганизуем вычисления поперечной силы Q слева (и QQ справа), изгибающего момента М слева (и MМ справа) от рассматриваемых точек.

1 точка границ: Расчет жесткого стержня


Расчет жесткого стержня


Расчет жесткого стержня

Расчет жесткого стержня


2 точка границ: Расчет жесткого стержня


Расчет жесткого стержня


Расчет жесткого стержня

Расчет жесткого стержня


3 точка границ: Расчет жесткого стержня


Расчет жесткого стержня


Расчет жесткого стержня

Расчет жесткого стержня


7. Численный метод решения СЛАУ - метод Гаусса


Численный метод Гаусса относится к точным методам решения системы линейных алгебраических уравнений. Он основан на приведении матрицы коэффициентов Расчет жесткого стержняк треугольному виду. Процесс поиска решения системы линейных алгебраических уравнений выполняется в два хода: прямой ход и обратный ход.

Прямой ход исключения переменных выполняется путём преобразования коэффициентов СЛАУ, коэффициенты при неизвестных обращаются в нуль, начиная со второго по формулам:


Расчет жесткого стержня; Расчет жесткого стержня; Расчет жесткого стержня, где

Расчет жесткого стержня; Расчет жесткого стержня; Расчет жесткого стержня


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

Обратный ход (последовательное нахождение неизвестных

Расчет жесткого стержня) выполняется по формулам:


Расчет жесткого стержня; Расчет жесткого стержня; Расчет жесткого стержня; Расчет жесткого стержня, где

Расчет жесткого стержня; Расчет жесткого стержня


В результате формируется матрица неизвестных: Расчет жесткого стержня Метод Гаусса для решения СЛАУ применим при условии, что все диагональные элементы матрицы Расчет жесткого стержня отличны от нуля, т.е. Расчет жесткого стержня, где Расчет жесткого стержня.


8. Обоснование применения метода Гаусса


Исходная СЛАУ имеет на главной диагонали элементы равные нулю:


Расчет жесткого стержня Расчет жесткого стержня Расчет жесткого стержня


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

Для того чтобы использовать численный метод Гаусса для решения данной СЛАУ необходимо её преобразовать. Для этого необходимо применить к исходной СЛАУ схему выбора главных элементов. В исходной СЛАУ переставим уравнения местами: первое уравнение поставим на второе место, второе уравнение поставим на третье место, третье уравнение поставим на первое место.

В результате на главной диагонали матрицы А отсутствуют члены равные нулю.


Расчет жесткого стержня


Для повышения точности получаемого решения СЛАУ матрица А должна быть диагонально преобладающей:


Расчет жесткого стержня, Расчет жесткого стержня Расчет жесткого стержня


Преобразованная СЛАУ имеет вид:


Расчет жесткого стержня

Расчет жесткого стержня

Расчет жесткого стержня


Условия применения метода Гаусса выполняются, следовательно, метод Гаусса можно использовать для решения преобразованной СЛАУ.


9. Блок - схема алгоритма


Расчет жесткого стержня


Расчет жесткого стержня

10. Программа


CLS

SCREEN 12

WINDOW (20, 20) - (-20, - 20)

N = 3

PRINT "Программу составил студент гр.320851 Клычников А.В."

50 PRINT " Расчет жесткого стержня "

PRINT " Исходные данные"

INPUT "Интенсивность распределения нагрузки q4 (кH/м) ="; q4

INPUT "Отрезок балки С1 (м) ="; C1

INPUT "Пролет балки L1 (м) ="; L1

INPUT "Отрезок балки C2 (м) ="; c2

INPUT "Пролет балки L2 (м) ="; L2

INPUT "Круговой момент M1 (кH*м) ="; M1

INPUT "Круговой момент M2 (кH*м) ="; M2

INPUT "Сосредоточенная сила P1 (кH) ="; P1

INPUT "Сосредоточенная сила P2 (кH) ="; P2

PRINT " "

IF C1 > 0 THEN GOTO 10 ELSE GOTO 40

10 IF c2 > 0 THEN GOTO 20 ELSE GOTO 40

20 IF L1 > C1 THEN GOTO 30 ELSE GOTO 40

30 IF L2 > c2 THEN GOTO 60 ELSE GOTO 40

40 PRINT "Ошибка ввода": GOTO 50

60 F = q4 * c2

DIM A (N, N), R (N), B (N)

A (1,1) = - (L1 - C1): A (1,2) = 0: A (1,3) = 0

A (2,1) = 0: A (2,2) = L1 - C1: A (2,3) = L2

A (3,1) = - (L1 - C1): A (3,2) = 0: A (3,3) = L2

B (1) = P1 * (L1 - C1) - M1 - F * (C1/2) - M2 - P2 * c2

B (2) = F * (L2 - c2/2) - M1 + P2 * (L2 - c2) - M2

B (3) = - P1 * (L1 - C1) - M1 + F * (L2 - c2/2) - M2 + P2 * (L2 - c2)

FOR I = 1 TO N - 1

FOR J = I + 1 TO N

A (J, I) = - A (J, I) / A (I, I)

FOR K = I + 1 TO N

A (J, K) = A (J, K) + A (J, I) * A (I, K): NEXT K

B (J) = B (J) + A (J, I) * B (I): NEXT J

NEXT I

R (N) = B (N) / A (N, N)

FOR I = N - 1 TO 1 STEP - 1: H = B (I)

FOR J = I + 1 TO N: H = H - R (J) * A (I, J): NEXT J

R (I) = H / A (I, I)

NEXT I

R1 = R (1): R2 = R (2): R3 = R (3)

X = R1 - P1 + R2

Y = R3 - P2 - F

PRINT " Результаты "

PRINT "Опорная реакция в точке 1 R1="; R (1); "kН"

PRINT "Опорная реакция в точке 2 R2="; R (2); "kН"

PRINT "Опорная реакция в точке 3 R3="; R (3); "kН"

PRINT "Y="; Y; " X="; X

PRINT

PRINT " Таблица ординат эпюр Q и M "

PRINT " S Q M QQ MM"

FOR s = 0 TO L1 + L2

IF s >= 0 AND s < C1 THEN

Q = 0

M = 0

GOTO 70

END IF

IF s > C1 AND s < L1 THEN

Q = R1 - P1

M = P1 * (s - C1) - R1 * (s - C1) + M1

GOTO 70

END IF

IF s > L1 AND s < L1 + L2 - c2 THEN

Q = 0

M = P1 * (L1 - C1) - R1 * (L1 - C1) + M1

GOTO 70

END IF

IF s > L1 + L2 - c2 AND s <= L1 + L2 THEN

Q = - P2 - q4 * (s - L1 - L2 + c2)

M = P1 * (L1 - C1) - R1 * (L1 - C1) + M1 + M2 + P2 * (s - L1 - L2 + c2) + q4 * (s - L1 - L2 + c2) * (s - L1 - L2 + c2) / 2

GOTO 70

END IF

IF s = C1 THEN

Q = R1 - P1

M = M1

QQ = R2 - P1 + R1

MM = - M1 - R2 * (L1 - s) + P2 * (L2 - c2) - M2 - R3 * L2 + F * (L2 - c2/2)

GOTO 80

END IF

IF s = L1 THEN

Q = R1 - P1 + R2

M = P1 * (s - C1) - R1 * (s - C1) + M1

QQ = R2

MM = P2 * (L2 - c2) - M2 - R3 * L2 + F * (L2 - c2/2)

GOTO 80

END IF

IF s = L1 + L2 - c2 THEN

Q = - P2

M = M2 + P1 * (L1 - C1) - R1 * (L1 - C1) + M1 + F * (L1 - C1) / 2 - 30

QQ = R3 - P2 - F

MM = - M2 - R3 * c2 + F

GOTO 80

END IF

70 PRINT USING "##. ## ####. #### ####. ####"; s; Q; M: GOTO 90

80 PRINT USING "##. ## ####. #### ####. #### ####. #### ####. ####"; s; Q; M; QQ; MM

90 NEXT s

A$ = INPUT$ (1)

LINE (10,8) - (18,8), 8

LINE (10,3) - (10, 20), 8

FOR Z = 10 TO 18 STEP.5

LINE (Z, 7.9) - (Z, 8.1), 8

FOR W = 3 TO 20 STEP.5

LINE (9.9, W) - (10.1, W), 8

NEXT W

NEXT Z

LINE (10, - 3) - (18, - 3), 8

LINE (10, 0) - (10, - 18), 8

FOR Z = 10 TO 18 STEP.5

LINE (Z, - 2.9) - (Z, - 3.1), 8

FOR W = - 18 TO 0 STEP.5

LINE (9.9, W) - (10.1, W), 8

NEXT W

NEXT ZFOR T = 0 TO L1 + L2 STEP.001

IF T >= 0 AND T < C1 THEN

Q = 0

M = 0

V1 = Q

U1 = M

GOTO 100

END IF

IF T > C1 AND T < L1 THEN

Q = R1 - P1

M = P1 * (T - C1) - R1 * (T - C1) + M1

V2 = Q

U2 = M

GOTO 100

END IF

IF T > L1 AND T < L1 + L2 - c2 THEN

Q = 0

M = P1 * (L1 - C1) - R1 * (L1 - C1) + M1

V3 = Q

U3 = M

GOTO 100

END IF

IF T > L1 + L2 - c2 AND T <= L1 + L2 THEN

Q = - P2 - q4 * (T - L1 - L2 + c2)

M = P1 * (L1 - C1) - R1 * (L1 - C1) + M1 + M2 + P2 * (T - L1 - L2 + c2) + q4* * (T - L1 - L2 + c2) * (T - L1 - L2 + c2) / 2

GOTO 100

END IF

100 PSET (T / 3 + 10, Q / 3 + 8), 4

PSET (T / 3 + 10, M / 3 - 3), 5

NEXT T

T = C1: GOTO 110

110 Q = R1 - P1

M = M1

PSET (T / 3 + 10, Q / 3 + 8), 4

PSET (T / 3 + 10, M / 3 - 3), 5

LINE (T / 3 + 10, V1/3 + 8) - (T / 3 + 10, Q / 3 + 8), 4

LINE (T / 3 + 10, U1/3 - 3) - (T / 3 + 10, M / 3 - 3), 5

T = L1: GOTO 120

120 Q = R1 - P1 + R2

M = P1 * (T - C1) - R1 * (T - C1) + M1

PSET (T / 3 + 10, Q / 3 + 8), 4

PSET (T / 3 + 10, M / 3 - 3), 5

LINE (T / 3 + 10, V2/3 + 8) - (T / 3 + 10, Q / 3 + 8), 4

LINE (T / 3 + 10, U2/3 - 3) - (T / 3 + 10, M / 3 - 3), 5

T = L1 + L2 - c2: GOTO 130

130 Q = - P2

M = M2 + P1 * (L1 - C1) - R1 * (L1 - C1) + M1 + F * (L1 - C1) / 2

PSET (T / 3 + 10, Q / 3 + 8), 4

PSET (T / 3 + 10, M / 3 - 3), 5

LINE (T / 3 + 10, V3/3 + 8) - (T / 3 + 10, Q / 3 + 8), 4

LINE (T / 3 + 10, U3/3 - 3) - (T / 3 + 10, M / 3 - 3), 5

END


11. Форма ввода - вывода информации


Программу составил студент гр.320851 Клычников А.В."

Расчет жесткого стержня

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

Интенсивность распределения нагрузки q4 (кH/м) = 10

Отрезок балки c1 (м) = 3

Пролет балки L1 (м) = 6

Отрезок балки c2 (м) = 2

Пролет балки L2 (м) = 12

Круговой момент M1 (кH*м) = 10

Круговой момент M2 (кH*м) = 35

Сосредоточенная сила P1 (кH) = 15

Сосредоточенная сила P2 (кH) = 30

Результаты

Опорная реакция в точке 1 R1=56.6668kН

Опорная реакция в точке 2 R2=-41.6667kН

Опорная реакция в точке 3 R3=50kН

Y=0 X=Расчет жесткого стержня

PRINT " Таблица ординат эпюр Q и M "


x Q M QQ MM

0.0000 0.0000 0.0000

1.0000 0.0000 0.0000

2.0000 0.0000 0.0000

3.0000 41.6667 10.0000 0.0000 0.0000

4.0000 41.6667 -31.6667

5.0000 41.6667 -73.3334

6.0000 0.0000 -115.0000 -41.6667 -115.0000

7.0000 0.0000 -115.0000

8.0000 0.0000 -115.0000

9.0000 0.0000 -115.0000

10.0000 0.0000 -115.0000

12.0000 0.0000 -115.0000

13.0000 0.0000 -115.0000

14.0000 0.0000 -115.0000

15.0000 0.0000 -115.0000

16.0000 -30.0000 -80.0000 0.0000 -115.0000

17.0000 -40.0000 -45.0000

18.0000 -50.000 0.0000

Проверка по оси X =0


Расчет жесткого стержняРасчет жесткого стержня


Программу составил студент Лазарев В.А. гр.320851


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


Эпюры поперечной силы Q и изгибающего момента М.

Q (kH) M (kHм)

Расчет жесткого стержня


Анализ результатов показал, что наиболее напряженное сечение стержня находится в точке с координатой S=14м, Q=-40 kH, M=-80kHм.


Литература


1. Данилина Н.И. Численные методы. - М.: Выш. шк. 1976г. - 368 с.

2. Дъяков В.П. Справочник по алгоритмам и программам на языке Бейсик для ПЭВМ. - М.: Наука, 1987г. - 240с.


Рефетека ру refoteka@gmail.com