Расчет балки

advertisement
Расчет балки –стенки методом конечных разностей (МКР).
1. Постановка задачи
На рис.1 представлена балка – стенка ( L, H  )
Рис.1Балка-стенка
Размеры балки - стенки: длина L  16 м , высота H  12 м , толщина
  1м , длина пилонов a  2 м .
Упругие постоянные материала: модуль упругости E  30000МПа ,
коэффициент Пуассона   0,12 .
Интенсивность поверхностной нагрузки q  3МН / м2 , длина участка
линейного изменения нагрузки Lq  4 м
Поверхностные силы по толщине балки-стенки распределены
равномерно, поэтому все слои балки будут деформироваться одинаково и
можно рассматривать НДС в окрестности точек срединной плоскости.
Закон изменения опорных реакций по длине пилона принимаем
постоянным (рис.2). Тогда из уравнения равновесия определяем
интенсивность опорных реакций q0
Рис. 2 К определению реакций опор.
2
 Fy  0 , (2 *3  4 / 2  3*8)  2 q 0 2  0 , q0  9МН / м
Уравнение совместности деформаций для функции напряжений
 ( x, y) в операторной форме имеет вид
2
2
 2 (x, y)
 2 ( x, y)
(1)

)(

)  0,
x 2 y 2
x 2
y 2
или в традиционной форме
(2)
 4 (x, y)
 4 ( x, y)  4 ( x, y)
.

2


0
x 4
x 2y 2
y 4
Уравнение (2) необходимо решать с учетом условий на контуре
балки-стенки. Для задания этих условий
используют «рамную
аналогию». Рассмотрим стержневую систему, очертание которой
совпадает с контуром балки-стенки (рис.3).
Интенсивность
поверхностных сил по толщине балки-стенки не меняется, поэтому в
стержневой системе переходим к погонным нагрузкам, собирая их по
толщине торцов (   1м ). Стержневая система находится в равновесии
под действием заданной нагрузки и реактивных сил (рис.3). На рис. 7
показана эпюра изгибающих моментов, а на рис.8 эпюра продольных сил.
(
Рис. 3 Рама с нагрузкой
Рис. 5 К построению эпюры N
Рис.4 Эпюра М, верхний ригель
Рис.6 Эпюра М, нижний ригель
Рис.7 Эпюра М изгибающих
Рис. 8 Эпюра N продольных сил в
моментов в раме
раме
В соответствии с рамной аналогией функция напряжений в точках
контура равна изгибающему моменту в раме, а производная от функции
по нормали к контуру равна продольной силе.
Первая группа условий имеет вид:
 (0, y)  M 01(y)  0 ,  (L, y)  M 3 2 (y)  0 ,
(3)
 (x,0)  M 03 (x) ,  (x,H)  M1 2 (x) .
Вторая группа условий запишется так:


 (0, y)  N01(y)  9 ,
 (L, y)  N3 2 (y)  9 ,
x
x
(4)


 (0, x)  N03 (x)  0 ,
 (H, x)  N12 (x)  0 .
y
y
Таким образом, необходимо найти функцию напряжений  ( x, y)
удовлетворяющую дифференциальному уравнению (2) и условиям на
контуре (3),(4).
2. Метод конечных разностей (МКР)
Приведем формулы численного дифференцирования для вычисления
первой и второй производных функции f ( x) в точке xi :
(5)
f( x  h)  f( xi  h)
f ( xi )  i
 0(h 2 ) ,
2h
f( x  h)  2  f( xi )  f( xi  h)
f ( xi )  i
 0(h 2 ) .
(6)
2
h
Формулы (5), (6) являются приближенными и имеют второй порядок
точности по h . Здесь h радиус окрестности точки xi .
Покажем, как практически пользоваться данными формулами.
Рассмотрим функцию f ( x)  x3 и вычислим 1-ю и 2-ю производные в
точке xi =1, используя формулы (5),(6).
Отступим от точки xi влево и вправо на величину h  0.1
(рис.9)
Рис. 9 Шаблон для вычисления производных
Вычислим первую производную
f( xi  h)  f( xi  h) 1,13  0,93
f ( xi ) 

 3,01 .
Точное
2h
2  0,1
производной
f (x)  3  x 2 , f (1)  3 ,
следовательно,
значение
погрешность
формулы (5) порядка h 2  0,01 .
Вычислим вторую производную
f( xi  h)  2  f( xi )  f( xi  h) 0,93  2 13  1,13
Точное
f ( xi ) 

6.
h2
0,12
значение второй производной f (x)  6  x , f (1)  6 , т.е. результат
вычислений по формуле (6) совпал с точным значением.
Покажем, как использовать формулы (5),(6) при вычислении
частных производных функции двух переменных в точке ( xi , y j ) .
 ( xi  h x , y j )   ( xi  h x , y j )

.
 ( xi , y j ) 
x
2h x
 ( xi , y j  h y )   ( xi , y j  h y )

 ( xi , y j ) 
.
y
2h y
2
 ( xi  h x , y j )  2   ( xi , y j )   ( xi  h x , y j )
x
h x2
 ( xi , y j  h y )  2   ( xi , y j )   ( xi , y j  h y )
 ( xi , y j ) 
2
2
 ( xi , y j ) 
2
(7)
(8)
.
.
y
h y2
Формулы (7)-(10) записаны на шаблоне точек (рис.10).
В формулах (7)-(10) введем обозначения :
xi 1  xi  h x , xi 1  xi  h x , y j 1  y j  h y , y j 1  y j  h y
(9)
(10)
(11)
Рис.10 Шаблон для вычисления частных производных
Сучетом обозначений (11) формулы (7)-(10) запишутся так:
 ( xi 1, y j )   ( xi 1, y j )

.
 ( xi , y j ) 
x
2h x
 ( xi , y j 1)   ( xi , y j 1)

 ( xi , y j ) 
.
y
2h y
2
 ( xi 1, y j )  2   ( xi , y j )   ( xi 1, y j )
x
h x2
 ( xi , y j 1)  2   ( xi , y j )   ( xi , y j 1)
 ( xi , y j ) 
2
2
(12)
.
 ( xi , y j ) 
y 2
h y2
Опустим аргументы, оставим только индексы,
 ( xi , y j )  i , j ,  ( xi 1, y j )  i 1, j ,  ( xi 1, y j )  i 1, j ,
 ( xi , y j 1)  i , j 1 ,  ( xi , y j 1)  i , j 1
(13)
в результате формулы можно записать в более компактном виде:
i 1, j i 1, j
i , j 1 i , j 1


i , j 
,
i, j 
y
2h y
x
2h x
2
i 1, j 2  i , j i 1, j
x
h x2
i , j 1 2  i , j i , j 1
, 
2 i j
2
,
i , j 
y 2
h y2
Формулы (14) записаны в соответствии с шаблоном (рис.12)
(14)
Рис.12 Индексный шаблон
3. Конечно-разностная аппроксимация уравнения совместности
деформаций и условий на контуре.
На прямоугольной области наведем сетку (рис.13). Вдоль оси x
L 16
область разобьем на 8 отрезков с шагом hx    2 м , а вдоль оси y на
8 8
H 12
6 с шагом hy    2 м . Таким образом, получим равномерную сетку,
6
6
так как hx  h y  h  2 м (рис.13).
Дифференциальное уравнение (1) запишем в произвольной
внутренней точке области ( xi , y j )
(
2

2
)(
 2i , j

 2i , j
)  0,(i  1..7; j  1..5)
x
y
x
y
Вторые частные производные заменим по формулам (14)
(
2
x 2

2
2
2
1
y
h2
)[
2
2
2
(i 1, j  i 1, j  4  i 1, j  i, j 1  i, j 1)]  0,
(15)
(16)
(i  1..7; j  1..5)
В квадратных скобках в формуле (16) записан конечно-разностный
оператор Лапласа. Если к каждому из 5 слагаемых применить конечноразностный оператор Лапласа (сдвигая индексы), то в результате получим
Неизвестные, входящие в уравнение (17) показаны на шаблоне (рис.14)
Рис.13 Дискретизация срединной плоскости балки – стенки
Рис. 14 Шаблон конечно-разностного уравнения в точке ( xi , y j )
Все 35 уравнений во внутренних точках области представлены ниже:
[ ( 2 0, 08 0, 12 0, 28 1, 020 1, 18 1, 22 2, 08 2, 12 2, 21, 33, 1
1, -1-1, 1 )
h 40, ( 2 1, 08 1, 12 1, 28 2, 020 2, 18 2, 22 3, 0
8 3, 12 3, 22, 34, 12, -10, 1 )
h 40, ( 2 2, 08 2, 12 2, 2
8 3, 020 3, 18 3, 22 4, 08 4, 12 4, 23, 35, 13, -11, 1 )
h4
0, ( 2 3, 08 3, 12 3, 28 4, 020 4, 18 4, 22 5, 08 5, 12 5, 2
4, 36, 14, -12, 1 )
h 40, ( 2 4, 08 4, 12 4, 28 5, 020 5, 1
8 5, 22 6, 08 6, 12 6, 25, 37, 15, -13, 1 )
h 40, ( 2 5, 0
8 5, 12 5, 28 6, 020 6, 18 6, 22 7, 08 7, 12 7, 26, 38, 1
6, -14, 1 )
h 40, ( 2 6, 08 6, 12 6, 28 7, 020 7, 18 7, 22 8, 0
8 8, 12 8, 27, 39, 17, -15, 1 )
h 40, ( 2 0, 18 0, 22 0, 3
8 1, 120 1, 28 1, 32 2, 18 2, 22 2, 31, 43, 21, 0-1, 2 )
h4
0, ( 2 1, 18 1, 22 1, 38 2, 120 2, 28 2, 32 3, 18 3, 22 3, 3
2, 44, 22, 00, 2 )
h 40, ( 2 2, 18 2, 22 2, 38 3, 120 3, 2
8 3, 32 4, 18 4, 22 4, 33, 45, 23, 01, 2 )
h 40, ( 2 3, 1
8 3, 22 3, 38 4, 120 4, 28 4, 32 5, 18 5, 22 5, 34, 46, 2
4, 02, 2 )
h 40, ( 2 4, 18 4, 22 4, 38 5, 120 5, 28 5, 32 6, 1
8 6, 22 6, 35, 47, 25, 03, 2 )
h 40, ( 2 5, 18 5, 22 5, 3
8 6, 120 6, 28 6, 32 7, 18 7, 22 7, 36, 48, 26, 04, 2 ) h 4
0, ( 2 6, 18 6, 22 6, 38 7, 120 7, 28 7, 32 8, 18 8, 22 8, 3
7, 49, 27, 05, 2 )
h 40, ( 2 0, 28 0, 32 0, 48 1, 220 1, 3
8 1, 42 2, 28 2, 32 2, 41, 53, 31, 1-1, 3 )
h 40, ( 2 1, 2
8 1, 32 1, 48 2, 220 2, 38 2, 42 3, 28 3, 32 3, 42, 54, 3
2, 10, 3 )
h 40, ( 2 2, 28 2, 32 2, 48 3, 220 3, 38 3, 42 4, 2
8 4, 32 4, 43, 55, 33, 11, 3 )
h 40, ( 2 3, 28 3, 32 3, 4
8 4, 220 4, 38 4, 42 5, 28 5, 32 5, 44, 56, 34, 12, 3 ) h 4
0, ( 2 4, 28 4, 32 4, 48 5, 220 5, 38 5, 42 6, 28 6, 32 6, 4
5, 57, 35, 13, 3 )
h 40, ( 2 5, 28 5, 32 5, 48 6, 220 6, 3
8 6, 42 7, 28 7, 32 7, 46, 58, 36, 14, 3 )
h 40, ( 2 6, 2
8 6, 32 6, 48 7, 220 7, 38 7, 42 8, 28 8, 32 8, 47, 59, 3
7, 15, 3 )
h 40, ( 2 0, 38 0, 42 0, 58 1, 320 1, 48 1, 52 2, 3
8 2, 42 2, 51, 63, 41, 2-1, 4 )
h 40, ( 2 1, 38 1, 42 1, 5
8 2, 320 2, 48 2, 52 3, 38 3, 42 3, 52, 64, 42, 20, 4 )
h4
0, ( 2 2, 38 2, 42 2, 58 3, 320 3, 48 3, 52 4, 38 4, 42 4, 5
3, 65, 43, 21, 4 )
h 40, ( 2 3, 38 3, 42 3, 58 4, 320 4, 4
8 4, 52 5, 38 5, 42 5, 54, 66, 44, 22, 4 )
h 40, ( 2 4, 3
8 4, 42 4, 58 5, 320 5, 48 5, 52 6, 38 6, 42 6, 55, 67, 4
5, 23, 4 )
h 40, ( 2 5, 38 5, 42 5, 58 6, 320 6, 48 6, 52 7, 3
8 7, 42 7, 56, 68, 46, 24, 4 )
h 40, ( 2 6, 38 6, 42 6, 5
8 7, 320 7, 48 7, 52 8, 38 8, 42 8, 57, 69, 47, 25, 4 ) h 4
0, ( 2 0, 48 0, 52 0, 68 1, 420 1, 58 1, 62 2, 48 2, 52 2, 6
1, 73, 51, 3-1, 5 )
h 40, ( 2 1, 48 1, 52 1, 68 2, 420 2, 5
8 2, 62 3, 48 3, 52 3, 62, 74, 52, 30, 5 )
h 40, ( 2 2, 4
8 2, 52 2, 68 3, 420 3, 58 3, 62 4, 48 4, 52 4, 63, 75, 5
3, 31, 5 )
h 40, ( 2 3, 48 3, 52 3, 68 4, 420 4, 58 4, 62 5, 4
8 5, 52 5, 64, 76, 54, 32, 5 )
h 40, ( 2 4, 48 4, 52 4, 6
8 5, 420 5, 58 5, 62 6, 48 6, 52 6, 65, 77, 55, 33, 5 ) h 4
0, ( 2 5, 48 5, 52 5, 68 6, 420 6, 58 6, 62 7, 48 7, 52 7, 6
6, 78, 56, 34, 5 )
h 40, ( 2 6, 48 6, 52 6, 68 7, 420 7, 5
8 7, 62 8, 48 8, 52 8, 67, 79, 57, 35, 5 )
h 40 ]
Заметим, что уравнения, записанные в точках сетки вблизи контура,
содержат неизвестные в законтурных точках. Например, в точке ( x1, y1)
( 2 0, 08 0, 12 0, 28 1, 020 1, 18 1, 22 2, 08 2, 12 2, 21, 33, 1
1, -1-1, 1 )
h 40
Неизвестные в законтурных точках исключаем в уравнениях,
используя условия на контуре (см.(4))


 (0, y)  N01(y)  9 ,
 (L, y)  N3 2 (y)  9 ,
x
x


 (0, x)  N03 (x)  0 ,
 (H, x)  N12 (x)  0 .
y
y
Например, на левой кромке контура

1
откуда
находим
0, j  (1, j  1, j )  N01(y j )  9 ,
x
2h
1, j  2h(1, j  9) , т.е. значения функции в законтурных точках
выражаются через значения функции во внутренних точках.
Значения функции  в точках контура задаются с учетом условий на
контуре (см.(3))
 (0, y)  M 01(y)  0 ,  (L, y)  M 3 2 (y)  0 ,
 (x,0)  M 03 (x) ,  (x,H)  M1 2 (x) .
Например, на кромке 0-1 и на кромке 1-2
0, j  M 0 1 (y j )  0 , i,6  M12 (xi ) .
Таким образом, выполнен переход от дифференциального уравнения в
частных производных (2) с условиями на контуре (3),(4) к системе
линейных алгебраических уравнений (СЛАУ) относительно неизвестных
значений функции во внутренних узлах сетки i , j . Получаемое решение
СЛАУ является приближенным и при уменьшении шага сетки h  0
сходится к точному (при этом число уравнений возрастает ). На рис.15
показан портрет матрицы СЛАУ.
Рис. 15 Портрет матрицы системы уравнений
Матрица системы симметричная и имеет разряженную ленточную
структуру с числом неизвестных в типовом уравнении равным 13.
На рис.16 показана функция напряжений, полученная в результате
решения СЛАУ.
На рис.17,18,19 представлены эпюры нормальных и касательных
 2
 2
 2
напряжений:  x  2 , y  2 , xy  
, вычисленные с помощью
xy
y
x
соответствующих конечно-разностных формул (14).
Рис. 16 Функция напряжений. Лицевой и торцевой контуры есть
очертание эпюр моментов в ригелях рамы. Углы наклона левой и правой
кромок равны продольной силе в стойках
Рис. 17 Напряжения  x
Рис.18 Напряжения  y
Рис.19 Напряжения  x y
Рис.20 Главное напряжение  max
Рис. 21 Главное напряжение  min
Рис. 22 Эквивалентное (приведенное) напряжение  i   max   min
На рис. 20, 21 представлены эпюры главных напряжений  max ,  min .
Для балки-стенки из хрупкого материала можно проверить прочность по
первой теории прочности
1   max  R ,  3   min  R  .
На рис.22 представлены эпюры эквивалентных напряжений по 3-ей
теории прочности (материал пластичный).
На рис.23 для балки стенки, загруженной сосредоточенной силой
показано векторное поле главных напряжений (половина балки)
Рис.23 Векторное поле главных напряжений
По направлению максимальных главных напряжений (растягивающих)
выполняется армирование железобетонной балки (рис.24)
Рис.4.13
Рис.24 Траектории главных напряжений при изгибе и схема армирования
железобетонной балки.
Download