Теория расчета пластин и оболочек. Методические указания
.pdf
M x k |
–0,267 |
|
M y |
|
|
–1 |
|
|
|
k |
|
|
|
||
|
|
|
|
|
|
|
|
–1 |
2,533 |
–1 |
|
–0,084 |
2,169 |
–0,084 |
|
|
–0,267 |
×5,47∙103 |
|
|
|
–1 |
×9,7∙103 |
|
|
|
|
|
|
||
|
|
1 |
|
|
–1 |
|
|
|
H k |
0 |
|
|
×1,53∙103 |
|
|
|
|
|
∆x |
∆y |
|
|
|
|
|
–1 |
1 |
|
|
||
Рис. 5. Операторы для вычисления моментов
Последовательно накладываем эти операторы на узлы сетки. Изгибающие моменты Mx, My посредине защемленных сторон:
MxR = 5,47∙103 (2,53∙wR – 2∙w6) = 5,47(–2∙3,49) = –38,2 кН∙м/м. MyT = 9,7∙103 (2,17∙wT – 2∙w2) = 9,7(–2∙3,27) = –63,4 кН∙м/м.
Крутящий момент H в углу шарнирно-опертых сторон (узел S):
HS = 1,53∙103 (–w7 – w7 – w7) = 1,53∙3(–4,03) = –18,5 кН∙м/м.
Аналогично вычисляются изгибающие и крутящие моменты в остальных узлах сетки. Числовые поля внутренних усилий Mx, My и H представлены на рис.6,7.
0 |
–8,8 |
–9,5 |
–5,9 |
0 |
x |
0 |
–58,7 |
–63,4 |
–39,3 |
0 |
x |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
15,9 |
9,3 |
5,0 |
|
|
–22,1 |
0 |
|
7,3 |
|
8,7 |
6,0 |
|
|
–3,3 |
0 |
34,8 |
19,2 |
9,8 |
|
|
–38,2 |
0 |
|
43,4 |
|
38,8 |
22,3 |
|
|
–5,7 |
0 |
23,1 |
16,6 |
8,7 |
|
|
–30,2 |
0 |
|
27,4 |
|
32,9 |
20,7 |
|
|
–4,5 |
0 |
0 |
0 |
0 |
0 |
|
|
0 |
|
0 |
|
0 |
0 |
0 |
|
|
y |
Mx, кН∙м/м |
|
|
|
|
y |
|
My, кН∙м/м |
|
|
|
|
|||
Рис. 6. Значения изгибающих моментов в узлах сетки
11
4,7 |
0 |
0 |
0 |
3,1 |
x |
17,1 |
8,9 |
–3,2 |
–8,9 |
0 |
|
3,1 |
1,9 |
–0,4 |
–1,9 |
0 |
|
–17,1 |
–8,9 |
3,2 |
8,9 |
0 |
|
–18,5 |
–13,9 |
3,9 |
13,9 |
4,3 |
|
y H, кН∙м/м
Рис. 7. Значения крутящих моментов в узлах сетки
Результаты расчета, а именно эпюры прогибов w и изгибающих моментов Mx, My для средних сечений плиты, эпюра крутящих моментов H в сечении по нижней грани плиты (y = b), а также величины наибольших напряжений в плите приведены на рис.8.
|
|
|
|
|
|
|
|
|
|
|
|
T |
|
|
|
|
|
|
|
|
|
x |
|
|
|
|
|
|
|
9,5 |
63,4 |
|
|||||||||||||||
|
|
|
|
|
|
F |
1 |
|
2 |
3 |
|
|
|
|
b/2 |
|
|
|
|
|
|
|
9,3 |
|
|
8,7 |
|
||||||||||||||||||||
|
4 |
|
|
|
5 |
|
6 |
|
|
|
R |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
7 |
|
8 |
|
|
9 |
|
|
|
|
|
|
b/2 |
|
|
|
|
|
3,3 |
|
|
|
|
|
19,2 |
|
|
|
|
|
|
|
|
|
||||||||||
|
S |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
38,8 |
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
16,6 |
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
5,8 |
|
|
|
|
|
|
|
|
|
|
|
32,9 |
|
||||||||
y |
|
|
|
|
|
a/2 |
|
|
|
|
|
|
|
a/2 |
|
|
|
4,5 |
|
Mx, кН∙м/м |
|
My, кН∙м/м |
|||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
w, мм |
|
|||||||||||||||||||||
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 w, мм |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
3,5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
5,5 |
|
|
5,8 |
38,2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
6 38,2 |
|
|
|
|
|
|
|
|||||||||||
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Mx, кН∙м/м |
x max |
|
|
|
|
5,7МПа |
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
9,8 |
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
34,8 |
|
|
19,2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0,22 103 |
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
5,7 |
My, кН∙м/м |
y max |
|
|
|
6 63,4 |
|
9,5МПа |
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||
|
43,4 |
|
|
38,8 |
22,3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0,22 103 |
|||||||||||||||||||||
|
|
|
|
|
|
|
|
3,9 |
|
13,9 |
|
|
|
|
|
|
4,3 |
|
|
H, кН∙м/м |
max |
|
|
|
6 18,5 |
2,8МПа |
|||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0,22 103 |
|||||||||||||||||||||||||||
18,5 |
13,9 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
Рис. 8. Результаты расчета плиты
12
Расчет плиты МКР на ЭВМ
Выполним расчет плиты на разных сетках узлов с помощью компьютерной программы, реализующей метод конечных разностей. Результаты расчета: прогибы в середине плиты и наибольшие значения моментов в соответствующих узлах приводятся в табл. 3. По приведенным результатам построены графики сходимости решения МКР для прогиба wср и изгибающего момента MyT (рис. 9).
|
|
|
|
|
Таблица 3 |
|
Результаты расчета плиты МКР |
|
|||
|
|
|
|
|
|
Сетка |
Число ненуле- |
wср, мм |
MxR, кН∙м/м |
MyT, кН∙м/м |
HS, кН∙м/м |
узлов |
вых неизвестных |
||||
|
|
|
|
|
|
4×4 |
9 |
5,77 |
–38,2 |
–63,4 |
–18,5 |
|
|
|
|
|
|
4×6 |
15 |
5,31 |
–35,5 |
–71,5 |
–19,2 |
|
|
|
|
|
|
4×8 |
21 |
5,12 |
–34,4 |
–74,5 |
–19,4 |
|
|
|
|
|
|
8×8 |
49 |
5,08 |
–51,3 |
–74,7 |
–22,3 |
|
|
|
|
|
|
8×12 |
77 |
4,95 |
–50,2 |
–77,1 |
–22,5 |
|
|
|
|
|
|
12×12 |
121 |
4,92 |
–54,9 |
–77,9 |
–23,1 |
|
|
|
|
|
|
w, мм
5,8 5,77
5,6
5,4 |
5,31 |
|
|
|
|
|
|
|
|
|
|
||
5,2 |
|
5,12 |
5,08 |
|
||
|
|
|
||||
5,0 |
|
Сетка узлов |
4,95 |
4,92 |
||
|
|
|
|
|||
4×4 |
4×6 |
4×8 |
8×8 |
8×12 |
12×12 |
|
MyT, кН
80 |
77,1 |
77,9 |
|
75 |
74,5 |
74,7 |
70 |
71,5 |
|
|
|
65 63,4
60
Сетка узлов
4×4 4×6 4×8 8×8 8×12 12×12
Рис. 9. Графики сходимости решения МКР
13
Расчет плиты МКЭ на ЭВМ
Выполним расчет плиты на разных сетках узлов с помощью компьютерной программы, реализующей метод конечных элементов. Результаты расчета: прогибы в середине плиты и наибольшие значения моментов в соответствующих узлах приводятся в табл. 4. По приведенным результатам построены графики сходимости решения МКЭ для прогиба wср и момента MyT (рис.10).
|
|
|
|
|
Таблица 4 |
|
Результаты расчета плиты МКЭ |
|
|||
|
|
|
|
|
|
Сетка |
Число ненуле- |
wср, мм |
MxR, кН∙м/м |
MyT, кН∙м/м |
HS, кН∙м/м |
узлов |
вых неизвестных |
||||
|
|
|
|
|
|
4×2 |
13 |
5,45 |
–56,3 |
–67,7 |
–22,6 |
|
|
|
|
|
|
4×4 |
33 |
5,18 |
–50,0 |
–79,3 |
–24,3 |
|
|
|
|
|
|
4×6 |
53 |
5,12 |
–48,9 |
–81,7 |
–25,4 |
|
|
|
|
|
|
8×4 |
73 |
4,96 |
–57,7 |
–75,4 |
–26,5 |
|
|
|
|
|
|
8×6 |
117 |
4,91 |
–56,2 |
–77,8 |
–27,3 |
|
|
|
|
|
|
8×8 |
161 |
4,89 |
–55,7 |
–78,7 |
–27,8 |
|
|
|
|
|
|
w, мм |
|
|
|
|
|
|
|
|
5,8 |
|
|
|
|
|
|
|
|
5,6 |
|
|
|
|
|
|
|
|
5,45 |
|
|
|
|
|
|
|
|
5,4 |
|
|
|
|
|
|
|
|
5,2 |
5,18 |
5,12 |
|
|
|
|
||
|
|
|
|
|
|
|
||
5,0 |
|
|
|
4,96 |
4,91 |
4,89 |
||
|
Сетка узлов |
|
|
|
|
|||
4×2 |
4×4 |
4×6 |
8×4 |
8×6 |
8×8 |
|||
MyT, кН |
|
80 |
79,3 |
75 |
|
81,7 |
78,7 |
77,8 |
75,4 |
70 67,7
65
60
Сетка узлов
4×2 4×4 4×6 8×4 8×6 8×8
Рис. 10. Графики сходимости решения МКЭ
14
Сопоставление решений, полученных МКР и МКЭ
Сравним точность и сходимость решений, полученных МКР и МКЭ. Для этого изобразим графики сходимости решения для прогиба w в середине плиты в зависимости от числа ненулевых узловых неизвестных (рис. 11). При определении числа неизвестных необходимо учитывать, что в МКР каждый узел содержит одно неизвестное (прогиб w), а в МКЭ – три неизвестных (прогиб и два угла поворота).
w, мм |
|
|
|
5,8 |
5,77 |
|
|
МКР |
|
|
|
5,6 |
|
|
|
5,45 |
|
|
|
5,4 |
МКЭ |
|
|
|
|
||
5,2 |
5,31 |
5,18 |
5,12 |
5,12 |
|
||
5,0 |
|
5,08 |
|
|
|
Число ненулевых неизвестных
0 10 20 30 40 50 60
4,96 |
|
|
|
4,92 |
|
|
4,95 |
|
|
|
|
|
|
|
4,91 |
||
|
|
|
|
||
70 |
80 |
90 |
100 |
110 |
120 |
Рис. 11. Графики сходимости решения МКР и МКЭ
Выводы
1.В рассматриваемом случае решения задачи изгиба плиты точность определения основных неизвестных (прогибов плиты) методами конечных разностей и конечных элементов практически одинакова.
2.Сходимость решения по МКЭ более равномерная и стабильная, чем по МКР.
3.Скорость сходимости решения по МКЭ при увеличении числа неизвестных также лучше, чем МКР.
4.Применительно к расчету изгибаемых плит прямоугольного очертания с закрепленными краями МКР менее сложен в реализации, чем МКЭ.
15
Расчетно-графическая работа № 2 РАСЧЕТ ОБОЛОЧКИ МЕТОДОМ КОНЕЧНЫХ ЭЛЕМЕНТОВ
Задание на РГР № 2
Тонкая прямоугольная в плане оболочка размерами 2a×2b и толщиной h шарнирно оперта по всем четырем сторонам и загружена собственным весом (γ=25 кН/м3), а также сосредоточенной силой F (рис. 12). Рассмотреть следующие варианты загружения:
1)сила посередине оболочки (x = 0, y = 0);
2)сила в середине верхней (или нижней) половины (x = 0, y = ±b/2);
3)сила в середине правой (или левой) половины (x = ±a/2, y = 0). Требуется:
1.Рассчитать оболочку методом конечных элементов в ПВК
SCAD или Лира-САПР. Для каждого варианта загружения выписать максимальные значения прогибов и напряжений. Выбрать невыгодный вариант загружения (по наибольшим растягивающим напряжениям).
2.Для невыгодного загружения изобразить эпюры прогибов, изгибающих моментов и нормальных сил по линиям главных кривизн, проходящих через точку приложения силы.
3.Выписать наибольшие значения компонент напряжений,
действующих в опасном сечении оболочки ( вx , нx , вy , нy ).
Построить эпюры нормальных напряжений. Определить вклад нормальных и изгибающих усилий в суммарные напряжения.
Исходные данные принять в соответствии с номером группы (табл. 5) и порядковым номером в списке группы (табл. 6).
|
z |
y |
3 |
1 |
b |
F |
|
f2 |
0 |
x |
|
2 |
|
|
b
f1
a a
Рис. 12. Прямоугольная в плане оболочка
16
|
|
|
|
Таблица 5 |
|
|
Исходные данные по группам |
||||
|
|
|
|
|
|
Номер |
E, кН/м2 |
|
a, м |
|
b, м |
группы |
|
||||
1 |
3,6∙107 |
0,20 |
16 |
|
12 |
|
|
|
|
|
|
2 |
3,7∙107 |
0,20 |
12 |
|
18 |
|
|
|
|
|
|
3 |
3,8∙107 |
0,20 |
18 |
|
14 |
|
|
|
|
|
|
4 |
3,9∙107 |
0,20 |
14 |
|
20 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Таблица 6 |
|||
|
Исходные данные по порядковому номеру в списке группы |
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
№ |
F = 150 кН |
№ |
F = 200 кН |
№ |
F = 250 кН |
|
||||||||
п/п |
|
|
|
п/п |
|
|
|
|
п/п |
|
|
|
|
|
f1, м |
f2, м |
h, м |
f1, м |
f2, м |
h, м |
f1, м |
f2, м |
h, м |
|
|||||
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
5 |
–3 |
0,45 |
11 |
5 |
|
3 |
0, 45 |
21 |
–4 |
6 |
|
0, 50 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
–2 |
5 |
0,40 |
12 |
–2 |
|
5 |
0,50 |
22 |
3 |
5 |
|
0,55 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
4 |
2 |
0,45 |
13 |
–3 |
|
6 |
0, 55 |
23 |
6 |
–3 |
|
0, 60 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
4 |
–3 |
6 |
0,50 |
14 |
–2 |
|
6 |
0,45 |
24 |
5 |
–3 |
|
0,55 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
5 |
5 |
3 |
0,40 |
15 |
6 |
|
–3 |
0,50 |
25 |
4 |
3 |
|
0, 50 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
6 |
6 |
–3 |
0, 50 |
16 |
5 |
|
3 |
0, 55 |
26 |
3 |
4 |
|
0, 60 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
7 |
–3 |
5 |
0,45 |
17 |
3 |
|
6 |
0, 50 |
27 |
–3 |
6 |
|
0,55 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
8 |
2 |
5 |
0,50 |
18 |
5 |
|
–2 |
0, 45 |
28 |
6 |
–4 |
|
0,50 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
9 |
–2 |
5 |
0, 45 |
19 |
3 |
|
5 |
0, 50 |
29 |
4 |
5 |
|
0,55 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
10 |
6 |
–2 |
0,40 |
20 |
6 |
|
–3 |
0,55 |
30 |
5 |
–3 |
|
0, 60 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
17 |
|
|
|
|
|
|
|
|
Пример расчета оболочки
Выполним расчет железобетонной прямоугольной в плане оболочки размерами 2a = 32 м, 2b = 24 м и толщиной h = 0,5 м. Стрелы подъема оболочки по сторонам: f1 = –3 м, f2 = 5 м. Параметры упругости: E = 3,8∙107 кН/м2, ν=0,2. Оболочка шарнирно оперта по контуру и загружена собственным весом γ= 25кН/м3, а также сосредоточенной силой F =200кН, которая может быть приложена в одной из трех точек: 1, 2, 3 (рис. 13).
Требуется:
рассчитать оболочку для трех вариантов загружения;
для невыгодного загружения изобразить эпюры прогибов, изгибающих моментов и нормальных сил по линиям главных кривизн;
для опасного сечения построить эпюры нормальных напряжений.
|
z |
y |
|
3 |
1 |
f2=5м |
b=12м |
|
|
||
2 |
0 |
|
x |
|
|
b=12 |
м |
|
f1=–3м |
a=16м |
a=16м |
Рис. 13. Схема заданной оболочки
Расчет оболочки в ПВК SCAD
Создаем поверхность оболочки согласно формуле
z (x,y) f |
|
x2 |
f |
|
y2 |
f |
|
3 |
x2 |
5 |
y2 |
2 (м), |
|
1 a2 |
2 b2 |
0 |
162 |
122 |
|||||||||
|
|
|
|
|
|
||||||||
где f0 = f1 + f2 = –3 + 5 = 2 м.
Начало системы координат x, y, z размещаем в центре плоскости проекции оболочки. Сетку конечных элементов назначаем 16×16 с разбиением каждой четырехугольной ячейки на четыре трехузловых плоских элемента.
18
Изополя перемещений по оси z (прогибы оболочки) для каждого варианта загружения показаны на рис. 14–16.















200
Рис. 14. Изополя перемещений для загружения № 1












200
Рис. 15. Изополя перемещений для загружения № 2
19








200
Рис. 16. Изополя перемещений для загружения № 3
Максимальные (по модулю) значения прогибов и наибольшие растягивающие напряжения по каждому варианту загружения сведены в табл. 7.
Таблица 7
Результаты расчета оболочки по вариантам загружения
Номер |
wmax, мм |
в |
|
, МПа |
н |
|
, МПа |
в |
|
, МПа |
н |
|
, МПа |
загружения |
|
x max |
|
x max |
|
y max |
|
y max |
|
||||
1 |
0,41 |
|
0,17 |
|
0,58 |
|
0,18 |
|
0,78 |
||||
|
|
|
|
|
|
|
|
|
|
||||
2 |
0,49 |
|
0,22 |
|
0,66 |
|
0,22 |
|
0,86 |
||||
|
|
|
|
|
|
|
|
|
|
||||
3 |
0,39 |
|
0,15 |
|
0,54 |
|
0,16 |
|
0,73 |
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Как видно, наибольшие растягивающие напряжения в оболочке имеют место при приложении сосредоточенной силы в точке 2 (см. рис. 13) (координаты узловой нагрузки: x = 0, y = –6 м). Таким образом, наиболее невыгодным является вариант загружения № 2 (расчетная схема для этого варианта изображена на рис.17).
20
