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

книги из ГПНТБ / Скорецкий, Э. С. Автоматизация отопительных систем

.pdf
Скачиваний:
0
Добавлен:
19.10.2023
Размер:
3.79 Mб
Скачать

М И Н И С Т Е Р С Т В О В Ы С Ш Е Г О И С Р Е Д Н Е Г О С П Е Ц И А Л Ь Н О Г О

О Б Р А З О В А Н И Я СССР

М О С К О В С К И Й О Р Д Е Н А Т Р У Д О В О Г О . К Р А С Н О Г О З Н А М Е Н И И Н Ж Е Н Е Р Н О - С Т Р О И Т Е Л Ь Н Ы Й И Н С Т И Т У Т

имени В. В. К У Й Б Ы Ш Е В А

С.Б . УХОВ

РА С Ч Е Т С О О Р У Ж Е Н И Й

ИО С Н О В А Н И Й

МЕ Т О Д О М К О Н Е Ч Н Ы Х

ЭЛ Е М Е Н Т О В

У ч е б н о е п о с о б и е

М О С К В А — 1973

і

У Д К 627.8:624.15.04(075.8)

r ' 7

 

В

пособии

излагаются

основы

метода

конечных

эле­

ментов

 

применительно

 

к

расчету

 

напряженно-деформи­

рованного

 

состояния

 

в

сложных

 

неоднородных

системах

типа

«сооружение

— основание».

 

Рассматривается

 

связь

метода

 

конечных

элементов с

вариационными

методами

решения

задач

 

теории,

 

упругости.

 

Приводится

 

вывод

основного

 

уравнения

 

для

случая

 

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

зада­

чи

и

процедура

 

решения

методом

конечных

элементов

плоской

задачи

 

для

 

неоднородных

 

изотропных

и транс-

версально-изотропных

 

 

тел. Даются

практические

 

реко­

мендации

для

решения

инженерных

задач,

иллюстрируе­

мые

рядом

примеров.

 

В

приложении

приводится

спра­

вочный

 

материал,

облегчающий

 

понимание

построений

метода

конечных

 

элементов,

и

алгоритм

программы

расчета

на

ЭЦВМ.

 

 

 

 

 

 

 

 

 

 

Учебное

 

пособие

 

 

предназначено

 

для

студентов

и

аспирантов

строительных

вузов.

 

Оно

будет

полезно

инженерам

и

 

научным

работникам,

занимающимся

расчетами

сооружений

 

и

оснований.

 

 

 

 

Сергей Борисович УХОВ

РАСЧЕТ С О О Р У Ж Е Н И Й И О С Н О В А Н И Й М Е Т О Д О М К О Н Е Ч Н Ы Х Э Л Е М Е Н Т О В

(Учебное пособие)

Редактор А. В. Светлова Корректор В. К. Колоткова

Сдано в набор 26/ХІІ 1972 г. Подписано к печати 10/ѴП 1973 г.

Л-101343

Цена 80 коп.

Объем 7,5

Тираж 1000

Заказ 120

 

Московская типография № 10 Соіозполнграфпрома

 

 

при Государственном комитета Совета Министров СССР

 

по делам издательств, полиграфии и книжной торговли.

 

 

Москва. М-114. Шлюзовая

наб., 10.

 

ВВ Е Д Е Н И Е

Внастоящее время во многих научно-исследователь­ ских, проектных институтах и вузах Советского Союза

созданы

и

успешно

работают вычислительные

 

центры

или отделы, обеспечивающие решение сложных

инженер­

ных задач

на

электронных

цифровых

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

машинах

( Э Ц В М ) .

 

 

 

 

 

Применение

Э Ц В М в

инженерных

целях

открыло

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

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

и создания

новых

методов

расчета напряженно - дефор ­

мированного состояния в сложных неоднородных

средах.

Одним из

таких методов, успешно развиваемых

в по­

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

долины или откос сооружения и т. п. могут

проводиться

не . только

на дорогостоящих и трудоемких моделях но

и путем

непосредственных аналитических

расчетов.

В последние годы в отечественной и зарубежной ли­

тературе

публикуется много работ, посвященных теоре­

тическим проблемам метода конечных элементов и конк­ ретным приложениям этого метода к решению широкого класса инженерных задач . Среди наиболее крупных обобщающих трудов следует отметить изданную за ру­

бежом работу Зенкевича и

Чьенга [35]

и

у нас в стра­

н е — работу Л . А.

Розина

[15]. Расчетам

бетонных и

земляных плотин,

взаимодействующих

с

основаниями,

подземных выработок, склонов долин и т. п.

посвящены

работы, выполненные в М И С И , В Н И И Г е ,

В О Д Г Е О ,

Гидропроекте, МГУ и других организациях [|3, 5,

I I , 12,

14,

18,

19—21].

 

 

 

В

М И С І І автором настоящего

пособия,

начиная

с 1970

г, основные положения метода

конечных

элемен­

тов для расчета взаимодействия гидротехнических

сооружений и их оснований

сообщаются

студентам

I V и V

курсов факультета гидротехнического

строительства

в

лекциях по механике

грунтов и

скальных

пород.

В

1971 —1972 г. выпускниками факультета

ГС М И С И вы­

полнялись дипломные работы, включающие расчеты ме­ тодом конечных элементов. В 1972 г. автором были про­ ведены семинары для преподавателей, аспирантов и науч­ ных работников факультета ГС М И С И и других органи­ заций. Составленный автором в 1970 году конспект лек­ ций, явившийся первой редакцией настоящего пособия, использовался дипломниками н аспирантами М И С И и МГУ для изучения основ метода конечных элементов и выполнения конкретных расчетов.

Большой интерес к активному овладению расчетом методом конечных элементов, проявленный преподавате­ лями, научными работниками, аспирантами и студента­

ми,

а т а к ж е почти полное отсутствие учебной

литерату ­

ры

по этому вопросу, побудили автора к

изданию

настоящего пособия. При этом преследовалась цель без

ущерба для строгости максимально упростить

изложе ­

ние теоретических основ н расчетных приемов

метода,

сделав его доступным для студентов и инженеров, не имеющих специальной математической подготовки.

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

В пособии приводятся практические рекомендации

для расчета некоторых

задач, содержащие, в частности,

предложения по реализации расчетов

на Э Ц В М . Теоре­

тические положения и

практические

рекомендации ил-

4

люстрируются примерами. В качестве справочного ма­ териала приводятся основные положения матричной алгебры, необходимые и достаточные для понимания математических преобразований, выполняемых при по­ строении метода конечных элементов, и алгоритм расче­ та одной из рассмотренных задач .

Работа выполнена в лаборатории механики скальных пород кафедры механики грунтов, оснований и фунда­

ментов М И С И под

общим руководством заведующего

кафедрой чл.-корр.

АН С С С Р Н. А. Цытовича и содер­

жит как собственные исследования автора, так и прове­ денное им обобщение материалов, изложенных в опубли­

кованных трудах

отечественных

и

зарубежных

специа­

листов. П а р а г р а ф 4.3 написан при

участии В . В . С е м е н о в а .

Им ж е , под руководством автора, были

разработаны

расчетные схемы, составлены программы и

реализованы

на Э Ц В М расчеты

большинства

примеров,

приведенных

в настоящем пособии.

 

 

 

 

Автор пользуется возможностью выразить искреннюю

признательность

профессорам,

докторам

технических

наук M . М. Гришину, Г. К. Клейну,

С. М. Слисскому и

В. С. Эрнстову, ознакомившимся

с

рукописью

работы

и высказавшим ряд ценных советов и замечаний, учтен­ ных при подготовке пособия к печати. Автор считает своим долгом отметить, что его работе над методом ко­

нечных

элементов в

значительной мере способствовала

научная

стажировка

в 1969 г. в Центральной лаборато ­

рии мостов и дорог

( П а р и ж ) .

Г л а в а I

Т Е О Р Е Т И Ч Е С К ИЕ ОСНОВЫ МЕТОДА КОНЕЧНЫХ ЭЛЕМЕНТОВ

§ 1.1. П О С Т А Н О В К А З А Д А Ч В Т Е О Р И И УПРУГОСТИ

Первым этапом расчета любого объекта является представление его в виде расчетной схемы. П р и этом ин­ женер сталкивается с весьма противоречивыми условия­ ми: во-первых, стремлением отразить в расчетной схеме все существенные особенности реального объекта, услож­ няющие ее, во-вторых, необходимостью разработки до­ статочно простой схемы, позволяющей применять для рас­ чета математический аппарат механики деформируемого тела. От того, насколько успешно удается удовлетворить эти условия, во многом зависит качество расчета и со­ ответствие результатов действительности.

В настоящее время наиболее распространенным д л я решения инженерных задач является аппарат теории упругости. Если задача решается в рамках теории упру­ гости, то обычно расчетная схема представляется в виде некоторого тела, ограниченного начальной (до деформи ­ рования) поверхностью с заданными законами распреде­ ления нагрузок на поверхности и внутри тела. Могут быть и более сложные случаи, например, дополнительное

задание закона

изменения

нагрузок во времени, однако

они в настоящей

работе не

рассматриваются .

Решение задачи в общем виде сводится к определе­ нию в каждой точке тела компонентов перемещений, на-

X пряжений, а иногда и деформаций:

« =

/. (•*• у.

2 ) ;

]

Зд: —/.,

(Х,

У,

Z)\

 

т .ту

\ 1

(Х,

У,

Z)\

(î - i)

s * = / , o ( - k .

 

У - Z

) ;

bx =

L(x,y,z).

 

)

б

В каждой точке тела, таким образом, определяется 15 неизвестных компонентов: три компоненты перемещений

(и, V, ад), шесть компонентов

напряжений

(ст.ѵ, Oy, <xr, хХу,

Xyz, "Г») и шесть компонентов

деформаций

Л-, ви,

ехху,

у У г , угх) • Поэтому для

решения

общей

задачи

нужно

иметь 15 уравнений, которые

можно было

бы применить

к каждой точке внутри

тела,

и

особые уравнения

(гра­

ничные условия), справедливые для любой точки по­ верхности, ограничивающей тело. Совместное решение этих уравнений, в принципе, позволяет определить на­ пряженно-деформированное состояние данного тела при данных нагрузках, т. е. напряженно - деформированное состояние расчетной схемы реального объекта.

Приведем основные уравнения линейной теории упру­

гости д л я решения задач при статическом

воздействии

нагрузок.

 

 

 

 

 

 

 

 

 

 

 

Уравнения

равновесия

(уравнения

статики)

ах

 

'

àу

'

dz

1 г

 

'

 

 

 

дъух

I

да,і

I dz,

 

 

 

 

 

 

дх

 

 

 

 

 

 

 

 

 

 

Геометрические уравнения

(уравнения

Коши)

 

 

du

 

 

du

I

du

I

 

 

s*

 

~дх'

I'v

— 'dJ'T'dT'

I

 

 

 

 

dv

 

 

dv

I

dai

I

 

 

 

 

âw

 

 

dw

I

du

I

 

 

S z

"35"'

^zx~~dT^"dT

 

 

)

 

 

Физические

уравнения

(обобщенный

закон Гука)

о я

=

2Gex

- ) - ЯѲ,

tXy

=

G^xy]

"J

 

oz

=

2Gez

+ 26,

X « = G y „ .

I

 

З д е с ь : G = 2 ( | E

+

^ ; Я = у ^ ; Ѳ = 3 з с р ; s c p = - L ( £ : ( + e y +

+ е г); P-X. ?У,

pZ — проекции

на

 

осп x, y, z объемной

силы, приходящейся

на

единицу

объема

тела .

7

геометрических уравнении:

В зависимости от принятого метода решения задач вместо геометрических уравнений (1—3) могут исполь­ зоваться уравнения неразрывности деформаций, полу­ ченные Сен-Венаном из

ô2l

ду

дх2

дх•ду

 

 

 

 

 

dz2

~І~

 

дуа'~'ду~дг'

 

 

 

 

 

д х 1

'

dz2

'àz-àx

 

 

 

 

(1-5)

 

 

-LEI'*.

 

 

 

дЧ

 

 

 

=2

 

 

dz

дх

dz

 

дх-ду

 

^

ду

 

 

 

_д_

 

 

 

 

-)=2

 

д2е

 

дх

ày

 

dz

дх

 

ду-dz

 

 

 

 

дх

atz,

=

2

д-^у

 

 

 

 

ду

 

 

дг-дх

 

Граничное условия могут быть заданы в виде компо­ нентов сил, приложенных на поверхности тела (статиче­ ские граничные условия):

рА .ѵ =

з х

cos {xv)

-L хху

cos (t/v) - f - ^« cos (zv);

Ï

 

py,

=

tyS

cos {xv)

- f ay cos (yv)

-f- xyz cos (zv);

}

( 1 " 6 )

Р г ѵ =

 

cos (л-ѵ) + Tz y cos (yv) 4 - az cos (zv). J

 

З д е с ь

р д Ѵ

p

,

p z w — проекции

на

оси л*, у,

z внешней

силы,

приходящейся

на

единицу

площади поверхности

тела;

cos(xv),

cos(yv),

c o s ( z v ) — н а п р а в л я ю щ и е

коси­

нусы, т. е. косинусы углов между нормалью к поверх­ ности тела и осями х, у, z в точке приложения силы.

Другим способом определения граничных условий является задание компонентов перемещений точек по­ верхности тела (кинематические граничные условия) . Возможно т а к ж е задание на одних участках поверхно­ сти компонентов перемещений точек, на других —- ком ­ понентов сил (смешанные граничные условия) .

В зависимости от того, какие неизвестные компонен- ' ты прежде всего необходимо определить, применяют различные методы решения задач теории упругости.

Метод перемещений. Принимают за основные. неиз­ вестные перемещения точек упругого тела:

" = Ы*> у, z), v = f2(x, у, z), w = fs{x, у, z) .

Тогда необходимо в физические уравнения (1-4) под­ ставить геометрические уравнения (1-3), т. е. выразить напряжения через перемещения. Затем полученные вы­ ражения подставить в уравнения равновесия (1-2). Получается система трех уравнений с тремя неизвест­ ными компонентами перемещений. Метод перемещений

наиболее часто применяется при решении

задач

мето­

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

 

 

 

 

 

 

Метод сил. Принимают за основные неизвестные на­

пряжения:

 

 

 

 

 

 

 

 

 

ax = L(x,

У,

z),

oy = fs(x,

у,

2 ) ,

az=zft(x,

у,

z),

*ч/ =

М * -

У'

z ) '

"cux =

ft{x,

У,

z),

Tzx = f0(x, у,

z).

Выразив

в

физических

уравнениях

(1-4)

деформации

через напряжения, подставляют их в уравнения нераз­

рывности

(1-5).

З а т е м ,

используя

преобразованные

уравнения

неразрывности

и уравнения

равновесия

(1-2),

получают

шесть

уравнений

с шестью неизвестными

ком­

понентами

напряжений (уравнения БельтрамИ) .

 

Смешанный метод. В этом случае за основные не­ известные принимаются некоторые из компонентов пе­ ремещений и напряжений .

Известно, что далеко не все задачи, представляю­ щие практический интерес, могут быть строго решены

методами классической теории

упругости с доведением

до относительно

простых формул. Во многих случаях

геометрическая

конфигурация

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

вия настолько сложны, что определение напряженно - деформированного состояния в замкнутом виде оказы­ вается крайне трудоемким, а иногда и невозможным. Это относится, в частности, ко многим задачам, с кото­ рыми приходится иметь дело инженеру. Например, имеется строгое решение методом классической теории упругости для плотины треугольного профиля, рассмат­ риваемой как бесконечный клин. Однако, строгое реше­ ние для треугольной контрфорсной плотины сложного сечения, работающей совместно с неоднородным осно­ ванием, в рамках классической теории упругости от­ сутствует.

В сложных случаях обычно прибегают к дальней­ шему упрощению расчетных схем, например, к замене единой системы «сооружение — основание» двумя от­ дельными системами «сооружение» и «основание» с за­ данием на границе между ними некоторого закона рас-

9

Соседние файлы в папке книги из ГПНТБ