- •Введение
- •Лекция 2
- •Точное решение
- •Лекция 3
- •Лекция 4
- •Лекция 5
- •2. Приближенные методы
- •2.1 Метод конечных разностей
- •Перейдем к нормальному усилию в рамках мкр:
- •Лекция 3
- •2.2 Метод Бубнова - Галеркина
- •Лекция 4
- •2.3 Метод Ритца – Тимошенко
- •Лекция 5
- •2.4 Метод наименьших квадратов
- •Лекция 6
- •2.5 Метод конечных элементов
Лекция 2
Точное решение
Уравнение (5) представляет собой дифференциальное уравнение второго порядка, решаемое в квадратурах, т. е. оно решается последовательным интегрированием. Проинтегрировав его дважды
(6)
мы получим решение с точностью до констант интегрирования С1 и С2.
Эти константы определяются из граничных условий (из условий закрепления торцов стержня). Граничные условия могут быть геометрическими (ГГУ), когда на торцах заданы перемещения, статическими (СГУ), когда на торцах заданы силы, или смешанными, когда на одном торце заданы ГГУ, а на другом – СГУ.
Рассмотрим примеры решения при различных граничных условиях.
Пусть, например, нагрузка изменяется по линейному закону
. (7)
Дифференциальное уравнение (5) примет вид (ЕА=const):
. (8)
Дважды интегрируем уравнение (7):
(9)
Рассмотрим геометрические ГУ (рис. 5 для варианта А).
Рисунок 5. Геометрические граничные условия (вариант А)
Подчиним решение (9) геометрическим граничным условиям
(10)
откуда
(11)
Внося полученные значения констант во второе соотношение (9), получаем выражения для перемещения
. (12)
Здесь введена безразмерная координата
. (13)
На основании соотношения (4) получаем выражение для продольного усилия
. (14)
Произведем ручной счет по формулам (13) и (15) в пяти равноотстоящих точках. Результаты представим в виде таблицы 1 и графиков (рис. 6 и 7).
Таблица 1
|
0 |
0.25 |
0.5 |
0.75 |
1.0 |
|
0 |
15/64 (0.234) |
3/8 (0.375) |
21/64 (0.328) |
0 |
|
1.0 |
13/16 (0.813) |
1/4 (0.25) |
-11/16 (-0.688) |
-2 |
Рис 6. Изменение продольного перемещения по длине стержня (вариант А)
Рис. 7. Изменение продольного усилия по длине стержня (вариант А)
Запишем программу, численно реализующую выражения (13), (15), на алгоритмическом языке «Паскаль».
PROGRAM STTR;
uses crt;
const
qo=1.0; l=1.0; EA=1.0; m=4;
var i: integer;
qoln,qolu,dx,x,x3: real;
u,N: array[1..m+1] of real;
BEGIN
clrscr;
qoln:=qo*l/6;
qolu:=qoln*l/EA;
dx:=l/m;
writeln;
writeln(‘ Результат решения’);
writeln;
writeln(‘Координата Перемещение Усилие’);
for i:=1 to m+1 do begin
x:=dx*(i-1)/l;
x2:=sqr(x);
u[i]:={qolu*}x*(1-x2);
N[i]:={qoln*}(1-3*x2);
writeln;
writeln(‘ x=’,x:5:3,
‘ u=’,u[i]*6:7:4,’ N=’,N[i]*6:7:4);
end;
readln;
END.