Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

Теория расчета пластин и оболочек. Методические указания

.pdf
Скачиваний:
0
Добавлен:
15.08.2026
Размер:
969 Кб
Скачать

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,6107

0,20

16

 

12

 

 

 

 

 

 

2

3,7107

0,20

12

 

18

 

 

 

 

 

 

3

3,8107

0,20

18

 

14

 

 

 

 

 

 

4

3,9107

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=

b=12м

 

 

2

0

 

x

 

 

b=12

м

 

f1=–

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

Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]