- •Общий ход расчета методом перемещений
- •Учет односторонней связи с основанием
- •Общий алгоритм расчета
- •2. Решение
- •Прямоугольные плиты на упругом основании
- •Общие положения и составление системы разрешающих уравнений
- •8.4.2. Матрица жесткости прямоугольного элемента плиты
- •Матрица жесткости прямоугольного кэ плиты
- •Матрица жесткости прямоугольного кэ плиты
- •Матрица жесткости прямоугольного кэ плиты
- •Матрица жесткости прямоугольного кэ плиты
- •Матрица усилий прямоугольного кэ плиты
- •Учет односторонней связи с основанием
- •Матрица жесткости прямоугольного кэ плиты к примеру 8.4
- •Матрица жесткости прямоугольного кэ плиты к примеру 8.4
- •Матрица жесткости прямоугольного кэ плиты к примеру 8.4
- •Матрица погонных усилий Ng прямоугольного кэ плиты к примеру 8.4
Матрица усилий прямоугольного кэ плиты
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Zi |
|
Dxμ |
2Dxa |
2Dμb |
–3Dxa/a |
Dxa |
0 |
0 |
0 |
0 |
–3Dμb/b |
0 |
Dμb |
|
||
|
–3Dxa/a |
–Dxa |
0 |
Dxμ |
–2Dxa |
2Dμb |
–3Dμb/b |
0 |
Dμb |
0 |
0 |
0 |
|
||
|
0 |
0 |
0 |
–3Dμb/b |
0 |
–Dμb |
Dxμ |
–2Dxa |
–2Dμb |
–3Dxa/a |
–Dxa |
0 |
|
Zj |
|
|
–3Dμb/b |
0 |
–Dμb |
0 |
0 |
0 |
–3Dxa/a |
Dxa |
0 |
Dxμ |
2Dxa |
–2Dμb |
|
||
= |
Dyμ |
2Dμa |
2Dyb |
–3Dμa/a |
Dμa |
0 |
0 |
0 |
0 |
–3Dyb/b |
0 |
Dyb |
|
||
|
–3Dμa/a |
–Dμa |
0 |
Dyμ |
–2Dμa |
2Dyb |
–3Dyb/b |
0 |
Dyb |
0 |
0 |
0 |
|
Zk |
|
|
0 |
0 |
0 |
–3Dyb/b |
0 |
–Dyb |
Dyμ |
–2Dμa |
–2Dyb |
–3Dμa/a |
–Dμa |
0 |
|
||
|
–3Dyb/b |
0 |
–Dyb |
0 |
0 |
0 |
–3Dμa/a |
Dμa |
0 |
Dyμ |
2Dμa |
–2Dyb |
|
||
|
–8Dka/b |
–Dkb |
–Dka |
8Dka/b |
–Dkb |
Dka |
–8Dka/b |
Dkb |
Dka |
8Dka/b |
Dkb |
–Dka |
|
Zl |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
В таблице 8.6: |
|
|
|
|
|
|
|
|
|
|
|
-
Учет односторонней связи с основанием
Из принятых положительных направлений усилий и перемещений (см. рис. 8.13) и поверхности прогиба (8.66) видно, что контакт нижней поверхности элемента с поверхностью упругого основания будет осуществляться при Zi > 0; Zj > 0; Zl > 0; Zk > 0; > 0; > 0; > 0; > 0; < 0; < 0; < 0 и < 0. В противном случае при допущении односторонней связи будет происходить отлипание нижней поверхности элемента от поверхности основания.
На основании этого диагональную матрицу единичных функций (8.21) можно записать в виде:
hg = [hi hi hi hi hj hi hi hj hj hi hi hj] (8.72)
Пример 8.4. Требуется рассчитать плиту, свободно стоящую на упругом основании (рис. 8.15) с коэффициентом постели k0 = 12,7 Н/см3 при нагрузке силами F1 = 23 кН и F2 = 46 кН. Материал плиты – бетон (жесткость E = 3,1∙106 Н/см2, μ = 0,15). Расчет произвести без учета односторонней связи с основанием.
Решение. 1. Определим жесткостные характеристики плиты (8.69):
μx = μy = 0,15; ν = 1 – μx μy =0,9775; G = 3,1∙106/2(1+ 0,15) = 1,349∙ 106 Н/см2;
D = Dx = Dy = 3,1∙106∙403/(12∙0,9775) = 169,139∙ 108 Н∙см;
Dk = 1,349∙106∙403/6 = 0,851D.
2. Составим дискретную схему плиты, разбив ее на конечные элементы. Рассматриваемая плита (см. рис. 8.15) при заданном загружении имеет четыре оси симметрии. В силу этого при расчете можно рассмотреть 1/4 часть плиты, а уравнения метода перемещений, используя диагональную ось симметрии, составить для 1/8 части. На рис. 8.16 показана сетка разбивки четверти плиты на конечные элементы со сторонами а = в = 100 см. В силу симметрии для 1/8 части плиты можно сформулировать следующие граничные условия: в узлах 1, 2, 3 әw/әy = 0; в узлах 4, 7, 9 әw/әy = әw/әx.
Таким образом, разрешающая система уравнений буде 21 – го порядка.
3. Вычислим матрицы жесткости конечных элементов, приняв за общий множитель величину D/a : матрицу - по таблице 8.2; матрицу - по таблице 8.3. Так размеры всех элементов одинаковы, то и матрицы жесткости будут одинаковы. Матрицы , и их сумма приведены, соответственно в таблицах 8.7 – 8.9. Так как общие оси координат совпадают о местными, то блоки матриц жесткости
4. В соответствие с заданной расчетной схемой (см. рис. 8.15) и схемой разбивки на конечные элементы (см. рис. 8.16) матрицы свободных членов системы уравнений для каждого узла будут:
5. Составим канонические уравнения в общем виде на основании (8.64) из условий равновесия каждого узла:
Узел 1. riiZ1 + rijZ2 + rilZ4 + rikZ5 + R1 =0;
Узел 2. rjiZ1 + (rjj+rii)Z2 + rijZ3 + rjlZ4 + (rjk+ril)Z5 + rikZ6 + R2 =0;
Узел 3. rjiZ2 + rjjZ3 + rjlZ5 + rjkZ6 + R3 =0;
Узел 4. 0,5(rli+rji)Z1 + 0,5(rlj+rjl)Z2 +0,5(rll + rjj + rii)Z4 +
+0,5(rlk+ rij + rjk + ril)Z5 + 0,5rikZ7 + R4 =0;
Узел 5. rkiZ1 + (rkj+rli)Z2 + rljZ3 +(rkl + rji)Z4 +(rkk+ rll + rjj + rii)Z5+
+(rlk+ rij)Z6 +(rjk+ rli)Z7 + rikZ8 + R5 =0;
Узел 6. rkiZ2 + rkjZ3 +(rkl+rji)Z5+(rkk+rjj)Z6 + rjlZ7 + rjkZ8 + R6 =0;
Узел 7. 0,5rkiZ4 +0,5(rkj+ rli + rkl + rji)Z5+ 0,5(rlj+rji)Z6 +
+0,5(rkk+ rll + rjj + rii)Z7+0,5(rlk+ rij + ril + rjk)Z8+0,5 rikZ9 + R7 =0;
Узел 8. rkiZ5 +rkjZ6 +(rkl + rji)Z7 +(rkk+rjj + 0,5rjl)Z8+rikZ8 + R8 =0;
Узел 9. 0,5rkiZ7 +0,5(rkj + rkl)Z8 +0,5rkkZ9 + R9 =0.
С учетом граничных условий в полученные выражения подставим матрицы неизвестных, пронумеровав их по порядку, начиная с узла 1:
Таблица 8.7