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

книги / Метод конечных элементов

..pdf
Скачиваний:
4
Добавлен:
12.11.2023
Размер:
8.24 Mб
Скачать

в общую систему координат и рассылается в массив канонических уравнений основной части в соответствии с номерами перемещений суперузлов. При формировании столбцов свободных членов основ­ ной части учитывается наличие супернагрузок, которые уже приве­ дены к узловым при обработке алгоритмических классов. Узловая нагрузка для определенного суперэлемента выбирается из соответ­ ствующего суперкласса, переводится в общую систему координат и рассылается по номерам суперузлов в соответствующий столбец свободных членов основной части. Затем производится полный рас­ чет основной части с выпечатыванием перемещений узлов и усилий в КЗ, непосредственно входящих у эту часть.

Перемещения и усилия К З входящих в суперэлементы, опреде­ ляются следующим образом: для рассматриваемого суперэлемента строится матрица направляющих косинусов и из соответствующего* суперкласса выбирается обработанная неполным прямым ходом си­ стема канонических уравнений. Бронируется q (по количеству ос­ новных нагружений) одномерных массивов Bt [1 niQi, в которых должны' разместиться столбцы свободных членов, а затем и значе­ ния перемещений. Если в некоторых нагружениях элемент имеет супернагрузку, то соответствующие ей столбцы (также обработан­ ные по неполному прямому ходу) выбираются из нужного алгорит­ мического класса (рис. 81 — столбцы В а , B i2, B is) и размещаются в соответствующие массивы. Затем выбираются перемещения, ко­ торые относятся к суперузлам и получены из расчета основной части, переводятся в местную систему координат и записываются в конце одномерных массивов. После этого, начиная с уравнения (nt — /210), производится обратный ход, в результате чего в масси­ вах B tj вычисляются искомые перемещения для всех узлов рассмат­ риваемого суперэлемента. Вычисление усилий и напряжений для К З еходящих в рассматриваемый суперэлемент, производится ра­ нее описанным способом. Перемещения и усилия выдаются нащечать

и происходит

переход к рассмотрению следующего суперэлемента..

Результаты

счета — перемещения, усилия, напряжения — вы­

даются на печать в табличной форме. Так, для перемещения указы­ ваются: номер узла, вид перемещения, номер нагружения. Для уси­ лий и напряжений указываются номер К З номер процедуры, покоторой вычисляется матрица жесткости этого К З номера первых, двух узйов. Затем для каждого нагружения выпечатывается вели­ чина напряжений в центральной точке К З а при необходимости (сб этом указывается в исходных данных в массиве КЭ) и величины узловых усилий. Индексация перемещений, напряжений и усилий соответствует общепринятым обозначениям теории упругости.

Сравнительные характеристики вычислительного комплекса «Мираж.», аналогичных зарубежных разработок «Сезам», «АСКА», «Титус» и отечественной разработки «МКЭ> приведены в табл. 39. Как видим, возможности вычислительногокомплекса «Мираж» мало в чем уступают возможностям вычислительных комплексов «Сезам», «АСКА», «Титус», несмотря на то что последние реализованы на

Возможности вычислительного комплекса

«Сезам»

«Титус» «АСКА»

«мкэ»

«Мираж»

Объект изучения и харак-

 

 

 

теристики вычислительного

созданного на ЭВМ

 

комплекса

 

 

 

«УНИВАК-1108»

«Минск-22»

(язык — рабочие

(язык — ФОРТРАН)

 

коды)

Конечный элемент стер- ж ень:

пространственный

с

жесткими

встав­

ками

 

с

учетом сдвига

 

на упругом основании

криволинейный

 

Плоское напряженное со-

стояние:

 

прямоугольник

 

треугольник

четы­

произвольный

рехугольник прямоугольник и тре­

угольник

с

промежу­

точными

узлами по

сторонам

 

 

 

Тонкая

плита:

 

 

прямоугольник

 

треугольник

 

четы­

произвольный

рехугольник

 

осно­

Плита

на упругом

вании:

 

 

 

 

прямоугольник

 

треугольник

упругого

элементы

 

основания,

работаю­

щие

за

 

пределами

плиты

 

 

 

Тонкая

оболочка:

 

прямоугольник

 

треугольник

 

четы­

произвольный

рехугольник

 

 

Трехмерное

напряженное

состояние:

 

 

 

тетраэдр

 

 

 

параллелепипед

призма произвольный выпук­ лый 8-угольник

ЕС Т Ь

ЕС Т Ь

ЕС Т Ь

Е С Т Ь

ЕС Т Ь

Е С Т Ь

ЕС Т Ь Н Е Т

ЕС Т Ь

ЕС Т Ь

Е С Т Ь

Е С Т Ь

ЕС Т Ь

Е С Т Ь

ЕС Т Ь

Н Е Т

Н ЕТ

Н Е Т

ЕС Т Ь

ЕС Т Ь

Е С Т Ь

Н Е Т

Е С Т ЬН Е Т

ЕС Т Ь

Е С Т Ь

Е С Т Ь

Е С Т Ь

ЕС Т Ь

ЕС Т Ь

ЕС Т Ь

Е С Т Ь ЕС Т Ь

ЕС Т Ь

ЕС Т Ь

Е С Т Ь

Е С Т Ь

ЕС Т Ь

•НЕТ

ЕС Т Ь

Н ЕТ

ЕС Т Ь

ЕС Т Ь

Н ЕТ

Н Е Т

ЕС Т Ь

ЕС Т Ь

ЕС Т Ь

Е С Т Ь ЕС Т Ь

Н Е Т

Е С Т Ь

ЕС Т Ь Н ЕТ

ЕС Т Ь

Н Е Т

- ЕС Т Ь

Н ЕТ

 

Н ЕТН Е Т

Н Е Т

ЕС Т Ь

Н ЕТ

Н ЕТ

ЕС Т Ь

Н Е Т

Н Е Т Н ЕТ

Н Е Т

ЕС Т Ь

Н Е Т

Н Е Т

Н ЕТ

Н Е Т

ЕС Т Ь

ЕС Т Ь

Е С Т Ь

Н ЕТ

Н ЕТ

ЕС ТЬ

ЕС Т Ь

Е С Т Ь

ЕС Т Ь

Н ЕТ

ЕС Т Ь

ЕС Т Ь

ЕС Т Ь

Н ЕТ

Н Е Т

Н Е Т

Н Е Т

ЕС Т Ь

ЕНСЕТТЬ

ЕС Т Ь

Н Е Т ЕС Т Ь

ЕС Т Ь

Н ЕТ

ЕС Т Ь

 

Н ЕТЕС ТЬ

ЕС Т Ь

Н ЕТ

ЕС Т Ь

ЕС Т Ь

Е С Т Ь

ЕС ТЬ

Н ЕТ

Н ЕТ

 

 

 

 

 

Возможности вычислительного комплекса

 

 

 

 

 

«Сезат»

«Титус»

сАСКА»

♦мкэ»

«Мираж»

Объект

изучения и харак­

 

 

 

 

 

теристики вычислительного

 

озданного на ЭВМ

 

 

комплекса

 

 

 

 

 

 

 

 

 

 

 

 

 

« УНИВАК-1108»

«Минск-22»

 

 

 

 

 

(язык — рабочие

 

 

 

 

 

(язык — ФОРТРАН)

коды)

 

 

 

 

 

 

 

 

Ограничения:

 

 

 

 

 

 

 

 

количество рангов су ­

 

 

 

 

 

перэлементной рекур­

6

0

0

0

1

сии

 

 

 

 

возможность решения

ЕС Т Ь

ЕС ТЬ

ЕС ТЬ

Н ЕТ

ЕСТЬ

композитных

 

систем

количество элементов

8 00

не ука­

6900

1000

1024

количество

узлов

со

 

зано

 

 

 

200

То же

465

1000

512

связями

загруже-

количество

6

 

720

3640

64

ний

 

 

суперуз­

 

количество

 

200

_

_

_

64

лов

в суперэлементе

количество

неизвест­

 

 

 

 

 

ных

для

 

основной

 

 

 

 

 

системы или

для

су ­

3 0 0 0

5000

14000

3000

3500

перэлемента

 

 

 

ширина ленты

 

1000

5000

14000

60

1700

Сервис:

 

 

попол­

 

 

 

 

 

возможность

 

Е С Т Ь

ЕС Т Ь

Е С Т Ь

ЕС ТЬ

ЕС ТЬ

нения библиотеки КЭ

блок поиска формаль­

ЕС Т Ь

ЕС Т Ь

ЕС Т Ь

Н ЕТ

ЕС ТЬ

ных

ошибок

информа­

сокращение

 

 

 

 

 

ции

при наличии

ре­

ЕС Т Ь

ЕС Т Ь

ЕС Т Ь

Н ЕТ

ЕС Т Ь

гулярностей

 

 

 

табличная выдача ре­

ЕС ТЬ

ЕСТЬ

ЕС Т Ь

ЕС Т Ь

ЕС ТЬ

зультатов счета

 

возможность

 

задания

 

 

 

 

 

исходной информации

 

 

 

 

 

в различных

систе­

Н Е Т

Н ЕТ

Н ЕТ

Н ЕТ

ЕС ТЬ

мах

координат

 

более мощной ЭВМ «УНИВАК-1Ю8» (по параметрам она сравнима с ЭВМ «ЕС-1065»).

Опыт эксплуатации программного комплекса «Мираж» лег в ос­ нову разработки вычислительного комплекса «Супер» на ЭВМ «Минск-32». Использование более мощной ЭВМ дало возможность значительно расширить сервисную часть и библиотеку конечных элементов, увеличить быстродействие, подключить новые блоки, которые автоматизируют расчеты, связанные со строительным про­ ектированием и основанные на положениях СНиПов. Так, кроме

расчета системы, состоящей из произвольного набора КЭ, на ста­ тические воздействия в виде местной или узловой нагрузки, задан­ ных перемещений и температурных воздействий, по вычислитель­ ному комплексу «Супер» можно: выполнить расчеты на сейсмиче­ ские воздействия и на ветровую нагрузку с учетов пульсации ветра; выбрать расчетные сочетания усилий для стержневых систем, а так­ же подобрать необходимую площадь арматуры для всех случаев напряженного состояния железобетонного сечения -или определить напряжения в заданных точках сечения для КЭ Из изотропного-ма­ териала. По вычислительному комплексу «Супер» можно решить, например, такие задачи: рассчитать конструкцию высотного здания, о учетом совместной пространственной работы диафрагм жестко­ сти, каркасных включений, плит перекрытий и фундаментной плиты на упругом основании на статические нагрузки, сейсмические воз­ действия в различных направлениях, ветровую Нагрузку с учетом, пульсации ветра; выбрать расчетные сочетания усилий; подобрать сечение арматуры для железобетонных конструкций; рассчитать сложную машиностроительную конструкцию на динамические воз­ действия в виде гармонических колебаний, удара или импульса

иопределить напряжения в сечениях ее элементов.

Впрограммах «Мираж» и «Супер» решен важный вопрос баланса между универсальностью комплекса и специализированностью про­ грамм. Так, по информации, содержащейся в исходных данных, управляющий блок производит настройку комплекса на решаемую задачу. При этом, с одной стороны, сохраняется возможность реше­ ния широкого класса задач, а с другой стороны, состав исходных, данных, выходные результаты, время решения конкретной частной задачи остаются характерными для специализированных программ. Это позволяет по современным представлениям классифицировать такие программы, как пакеты сложной структуры со специализи­ рованным монитором.

Внастоящее время появилось много вычислительных комплек­

сов, выполненных на уровне промышленных разработок. К наибо­ лее важным из них следует отнести вычислительные комплексы «Прочность» (Киевский инженерно-строительный институт, ЭВМ БЭСМ-6), «Рассудок» (организация «Киевпроект», ЭВМ. М-222), «Гамма» (Киевский зональный НИИ экспериментального проекти­ рования, ЭВМ БЭСМ-6), «Каскад» (Ленинградский кораблестрои­ тельныйинститут, ЭВМ «Минск-32»), «Корпус» (Ленинградский промстройпроект, ЭВМ М-222), «Марс-ЕС» (Центральный НИИ проектирования и автоматизации строительства, ЭВМ ЕС), «Марк» (фирма IBM, ЭВМ IBM-370), «Настран» (Ведомство морского и воздушного флота США, ЭВМ IBM-370).

Коллектив вычислительных комплексов «Супер» и «Рассудок» разработал комплекс «Лира» на ЕС ЭВМ, который во многом близок к вычислительному комплексу «Супер», но значительно превосхо­ дит’ его в сервисной части, быстродействии и удобстве расшире­ ния,

А р г и р и с Д

ж . Энергетические теоремы

и расчет

 

конструкций.—

В

кн .:

Современные

методы

расчета

 

слож ных

статически

 

неопределимых

систем

/П од ред.

А. П . Фи

л

и н а .— Л .: Судпромгиз,

1961.

 

 

конструкции

А р г и р и с Д ж . Современные

достижения

в

методах

расчета

с применением

матриц

/П од

ред.

А.

ФС. м и р н о ва .— М .:

Стройиздат,

 

1968.

 

х о в Н . И .,

Л

у ж и н О. В . Приложение

методов

теории

упругости

Б е з у

и пластичности-к

решению инженерных зад ач .— М .: Высш ая школ а ,">1974.

 

Д е к л у

Ж . Метод конечных

элементов.—

М .: Мир, 1976.

 

 

 

 

 

 

 

Е р ж а н о в Ж .

С. ,

К а р

и м б а е в Т Метод.Д. конечных элементов в зада­

 

чах

механики

горных пород.—

А лм а-А та: Н аука,

1975.

 

1975.

 

 

 

З е н к е в и чО. М етод элементов

в техн и к е .—

М .: М .:

Мир,

 

 

 

М

а с л е н н и к о вА .

Н . Метод

конечных

элементов.— В

кн.: Справочник

по

теории

упругости

/П од

ред. П .

МВ. а

р

в

а

к

а

 

и

А.

Ф.

Р я б оа .в—

 

К иев: Б удивэльник,

1 97 К

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

М а с л е н н и к о вА.

Н . Расчет строительных конструкций методом ко­

нечных

элем ентов.—

Л .: И зд -во Л ен и н !р . инж .-строит, ин-та,

1977.

 

 

 

М етод

конечных элементов в механике сплошной среды /П од

ред. П. М. В а р-

в а к а .—

Киев: И зд-во

Киев,

автодор. ин-та,

1976.

 

 

с

использованием

Методы

расчета

стерж невы х

 

систем, пластин

и

оболочек

ЭВМ : В 2-х ч. /П од ред.А . Ф . С м и р н о ва .—

М .: Стройиздат,

1976.

 

 

М .:

М и х л

и н

С.

 

Г .

 

Вариационные

методы

 

в

математической

физике.—

Н ау к а,

1970.

 

 

 

 

 

 

 

 

методы расчета

судовы х

конструкций.—

Л .:

П о с т н о в В .А. Численные

Судостроение,

1977.

 

 

Х а р х у р иИм.

Я .

 

Метод

конечных

 

Элементов

в

П о с т н о в

 

В.

А. ,

 

 

,

расчетах

судовы х

конструкц и й .— Л .: Судостроение,

1974.

 

 

 

конеч­

Р

о з и н Л . А. Расчет'гидротехнических

сооружений на ЭЦВМ : Метод

ных

элем ентов.—

Л .:

Энергия,

1971.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Л .:

Р

о з

и н Л . А . Стержневые

системы как системы конечных элементов.—

И зд -во

Л енингр. гос.

ун -та, 1976.

 

 

 

в

приложении

к

упругим

систе­

Р о з и н Л .

А . М етод

конечных

элементов

м ам .—

М .: Стройиздат,

1977.

 

 

 

 

задач

для

упруги х

си стем .—

Л .:

Р

о з

ин

Л . А. Вариационная постановка

И зд-во

Ленингр. го с.-у н -та,

1978.

 

Г .

А.

М етод

ск1нченних

елемент1в :—

 

С

а в

у

л

а Я . Г ., Ш

и

н

к

а

р е н к о

 

Льв1в: Вид-во Льв1в. держ’. ун-ту,, 1976.

 

конечных элементов.—

М .: Мир,

С t р е н г Э .,

Ф

и к с

Д ж .

 

Теория метода

1977.

 

 

 

 

 

 

 

С.

 

П. ,

 

Г у д ь еД р.

Теория

уп ругости .—

М .:

Н аука,

Т и м о ш е н к о

 

 

 

1975.

 

 

 

 

 

 

 

 

 

 

в

статике стерж невы х

си стем .—

Л .:

Стройиздат,

Ф и л и н А . П . Матрицы

B

a t

h e

K .- J .,W i

l s

o

n

E .

N um erical

M ethods

in

F in ite

Elem ent A naly­

s is .—

New Jerse y :

P re n tice -H all,

1975.

 

 

A nalysis: F u n d am en tals.—

New

G a l l a g h e r R i c h a r

d

HFin. ite

Elem en t

Jersey :

P re n tice -H all,

 

1976.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

К

о l a

г

V ., K r a

t

o

c hv i l

J .

,

L e

i t

n

e

г

F .,

Z

e n

T s e

k

A.

V ypocet

plosnysh

a

prostorovych

k o n stru k d

m etodou

 

konelm ych

p rv k u .—

P rah

S N T L ;

1972.

 

 

L

a

r

r y AJ pplied.

F in ite

Elem ent

A n alysis.— London:

Jo h n

S

e g

e

r

l i

n d

 

W iley ,

1975.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Предисловие

• .

 

 

 

 

.

 

 

.

 

 

. .

 

 

3

 

Глава 1. Фундаментальные вариационные принципы и методы

 

 

 

 

§

1.

Вариационные

принципы

 

.

 

 

.

.

 

 

 

 

4

 

§

2.

Метод

Рэлея —

Ритца

 

. .

 

 

 

 

17

 

§

3 .

Метод

Бубнова— Галеркина

 

 

 

 

 

 

25

 

Глава

 

2.

Метод

конечных

элементов

(МКЭ)

 

 

 

 

 

 

 

§

4.

Основные

положения

М КЭ

 

 

 

 

 

 

 

 

 

28

 

§ 5. Схема решения задач по МКЭ

 

 

 

 

 

 

 

 

38

 

Глава

 

3.

Стержневые системы

 

 

 

 

 

 

 

 

 

 

 

 

§ 6.

Основные

гипотезы.

 

Функционал

полной потенциальной

энергии41

 

§

7 .

 

Конечный

элемент

растянутого

(сж атого)

стерж ня

 

. .

42

 

§

8.

Конечный

элемент

изгибаемого

стержня . . . . . .

47.

 

§

9.

 

Конечный элемент стержня в условиях пространственной

нагрузки52

 

§

10.

 

 

Примеры

 

расчета

стержневых

систем

 

 

. .

 

 

55. .

 

Глава 4.

Плоская

задача

теории

упругости

 

 

 

 

энергии61

 

§ 1 1 .

 

Основные гипотезы. Функционал полной потенциальной

.

§

12.

 

 

Прямоугольный

КЭ

 

* . .

 

 

 

 

 

62

§

13.

 

 

Треугольны й

КЭ

 

. .

 

 

 

 

 

70

.

§ 14. Примеры расчета балок-стенок

 

 

 

 

 

 

 

75 .

 

Глава 5.

Пластины

 

 

 

Функционал

полной

потенциальнойэнергии

80

 

§

15.

 

 

Основные

 

гипотезы.

 

§

16.

 

 

Прямоугольный

КЭ

с двенадцатью степенями

свободы

 

.

82»

 

§

17*

 

 

Прямоугольный КЭ с шестнадцатью степенями свободы

 

, .

91.

 

§

18.

 

 

Треугольны й

КЭ

плиты

пластин

 

 

 

 

 

 

93

 

§

19.

 

 

Примеры расчета

тонких

 

 

 

 

 

 

 

96

 

Г лава

 

6.

Многослойные

балки

и

пластины

 

 

 

 

 

 

 

§

20 .

Основные гипотезы и функционал полной потенциальной

энергии

 

§

21 .

для

балки

 

элемент

 

 

.....................

 

»

 

 

 

100

 

 

 

Конечный

 

многослойной

 

балки

 

 

 

 

106

 

§

22 .

Основные гипотезы и функционал полной потенциальной

энергии

 

§ 2 3 .

для

пластины

 

 

.

 

.

 

. . . . .

 

 

ПО.

 

 

 

Прямоугольный

КЭ

многослойной' пластины

.

 

 

114.

 

§

24 .

 

 

Примеры

 

расчета

многослойных

 

балок

и пластин

 

 

117

 

Глава

 

7. Массивы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

§

2 5 .

 

 

Основные гипотезы. Функционал потенциальной энергии

 

121

 

§

26 .

 

 

Конечный

элемент

в форме тетраэдра

.

*

 

 

 

122

 

§

2 7 .

 

Конечный элемент в форме параллелепипеда

 

 

 

125

 

Глава

 

8.

Оболочки

 

 

 

 

 

 

 

 

 

 

 

 

 

 

127

 

§

28 .

 

Основные гипотезы. Функционал потенциальной энергии

 

 

 

§

2 9 .

 

Прямоугольный КЭ оболочки

двоякой

кривизны

 

 

128

 

§

30.

 

 

И спользование плоских

КЭ для

расчета

оболочки

 

 

131

 

Глава

 

9.

Комбинированные

системы

 

 

 

 

 

 

 

 

132

 

§

31.

 

 

П ластина,

подкрепленная

ребрами

 

 

 

 

 

 

§

32.

 

 

Пластина

 

на

упругом

основании

 

 

 

 

 

134

 

§

33.

 

 

Балка

на

упругом

 

основании

 

 

 

 

 

 

 

 

141

 

Глава

 

10.

Численная

реализация

МКЭ

 

 

 

 

 

 

 

 

 

§

34 .

 

Основные

 

вопросы

реализации МКЭ на современных ЭВМ . . .

 

§

35.

 

Алгоритм

 

построения матриц

жесткости

 

.

 

 

 

1 4 6

 

 

 

 

 

 

 

147 .

 

§

36.

 

Составление систем

 

канонических

 

уравнений

 

 

 

157

 

§

3 7 .

Реш ение

систем

 

уравнений

высоких-

порядков

КЭ

 

1 6 0

 

§

38 .

 

Определение

напряжений

и

деформаций

по

области

 

162

 

§

39 .

 

Применение

суперэлементов

для

расчета

слож ны х систем

1 6 2

 

§

40 .

Сущ ествую щ ие

вычислительные

комплексы

 

 

 

165

 

Список литературы

 

 

 

 

 

 

 

 

 

 

 

 

 

 

175

 

Соседние файлы в папке книги