книги из ГПНТБ / Скорецкий, Э. С. Автоматизация отопительных систем
.pdfМ И Н И С Т Е Р С Т В О В Ы С Ш Е Г О И С Р Е Д Н Е Г О С П Е Ц И А Л Ь Н О Г О
О Б Р А З О В А Н И Я СССР
М О С К О В С К И Й О Р Д Е Н А Т Р У Д О В О Г О . К Р А С Н О Г О З Н А М Е Н И И Н Ж Е Н Е Р Н О - С Т Р О И Т Е Л Ь Н Ы Й И Н С Т И Т У Т
имени В. В. К У Й Б Ы Ш Е В А
С.Б . УХОВ
РА С Ч Е Т С О О Р У Ж Е Н И Й
ИО С Н О В А Н И Й
МЕ Т О Д О М К О Н Е Ч Н Ы Х
ЭЛ Е М Е Н Т О В
У ч е б н о е п о с о б и е
М О С К В А — 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