книги / Проектирование стальных мостов с учетом пластических деформаций
..pdfоткрытого контура, описываемая уравнением стесненного круче ния:
Ф1У (* )- С /« Ф' (*) = - т (2). |
(2.4) |
а также теория тонкостенных стержней замкнутого контура, опи сываемая уравнением:
Ф1У (г)—к- Ф" (г) |
(г) |
т " (г) |
, при к — |
|
|
ОУк |
|
(2.5)
где Jc — направленный полярный момент инерции сечения.
Между двумя типами стержней открытого и замкнутого контура лежит широкий класс конструкций смешанного типа, например, открытый контур, усиленный связями, или система коробчатых се чений, связанных между собой плитой в одном (верхнем) уровне. Универсальная теория таких стержней отсутствует, видимо, по при чине сложного описания депланации при кручении стержней со смешанным сечением. Практически для таких типов сечений можно использовать одно из уравнений (2.4) или (2.5), при этом погреш ность решения остается неустановленной (в пределах различия пере численных теорий).
Теория тонкостенных стержней позволяет эффективно и доста точно просто рассчитать сложные системы мостов (неразрезные, ван тово-балочные) и получать при этом удовлетворительные результаты. Однако следует отметить невозможность по этому методу учета де формаций контура поперечного сечения и неравномерности распре деления напряжений по ширине поперечного сечения при симметрич ной нагрузке. В сравнении с классической теорией в сечениях тонко стенного стержня учитывают дополнительные силовые факторы для открытого (замкнутого) профиля:
бимомент
Я /в |
|
Я /в |
Ф " (г) или В{ю |
Ф " Ю+ |
т; |
изгибно-крутящий момент |
я /в |
я /в |
|
------я / а Ф '" (г) или Мв = - - ^ Ф |
' ' ' + - ^ т \ |
Дополнительный (седьмой) кинематический фактор сечения — депланация 0 = Ф' (г) для открытых типов сечений и мера деплана ции для замкнутых р = Ф' (z)/p —- М„/(рОУс).
Таким образом, учет пространственной работы по теории тонко стенных стержней приводит к дополнительным нормальным и каса тельным напряжениям, действующим в сечении. Это результат стес ненного (изгибного) кручения стержней.
Теория п л и т н о - б а л о ч н ы х конструкций (складчатых оболочек) базируется на решениях прикладной теории упругости, в частности пластинок. Теория позволяет достаточно полно выявить работу всех элементов пролетного строения — главных балок, плит, стенок балок и др. Однако эта теория не лишена недостатков — трудности учета переменного сечения по длине, граничных условий в сложных системах.
Метод к о н е ч н ы х элементов получил достаточно широкое распространение в последние 10—15 лет благодаря своей универ сальности. Это численный метод, основанный на условном разделе нии упругого тела (конструкции) на ряд элементов, решение для ко торых имеет простую форму, а в ряде случаев элементарно. Объеди няются элементы в узлах, в которых полностью удовлетворяются условия неразрывности перемещений и условия равновесия.
Рассмотрим конечный элемент прямоугольной формы при плос ком напряженном состоянии и изгибе для расчета системы, состоя щей из множества таких элементов [60].
Перемещения точек прямоугольного элемента в е г о п л о с к о с т и при плоском напряженном состоянии ддя ортотропной плас тины может быть задано в таком виде:
и (*. = |
(2/4* + 6/бху)- |
Рх |
(З/i х2+ 2/о x)+f7y - |
|
|
Е-х |
|
|
|
|
—~ТГ“ h 1/2+/а; |
( 2. 6) |
||
|
1 |
|
Иг/ |
|
|
|
(2/4 y+ 3/e х/а)— /7 *— |
||
ü (*» «/) = —^— (6/i xy + 2f2 у) —— |
||||
|
Ьу |
|
t x |
|
|
11 |
ftx |
-«•+/.. |
|
|
G |
Эти перемещения получены из решения плоской задачи теории
упругости, когда |
функция напряжений задана в виде <р (х, у) = |
= /х*3 + /2*2 + f |
3*y+hy2 + hy3> причем совместность деформаций |
обеспечивается при любых коэффициентах внутри элемента. По стоянные коэффициенты /0, f7 и /8 в выражении (2.6) .возникли при получении перемещений. При решении плоской задачи методом пере мещений в каждом узле вводят два неизвестных линейных переме щения U и Vi. Постоянные коэффициенты fu ..., /8 могут быть лёгко выражены через перемещения четырех узлов элемента, если соста вить соответствующую систему уравнений (для каждого узла по два уравнения). Напряжения внутри элемента проще всего найти через функцию напряжений ср (х, у).
Важным моментом является получение м а т р и ц ы |
ж е с т |
к о с т и элемента, которая необходима для составления |
систем |
канонических уравнений метода перемещений для рассчитываемой конструкции.
Порядок составления матрицы жесткости для прямоугольного элемента толщиной Л, шириной а и высотой b следующий.
В уравнения (2.6) для и и и подставляют координаты узлов эле мента, что дает матрицу А восьмого порядка. Получают обратную матрицу А”1, из которой берут пять строк в соответствии с числом постоянных fi в функции напряжений, т. е. рассматривают в даль нейшем матрицу Ag"1. Вводят матрицу физикомеханических пара метров элемента, в данном случае из ортотропного материала
с |
|
Ex |
(Ру Еу) |
0 |
”| |
Y |
(p.x Ex) |
Ey |
О |
I |
|
|
0 |
0 |
vG J |
||
|
|
||||
при р х Е а = Ру Ey, v = |
l — |
|
|
|
|
Матрицу жесткости рассматриваемого прямоугольного элемента
получают из выражения [60]: |
|
|
|
|
||
[ |
0.5а |
0.5b |
-| |
А5- ' . |
|
|
|
f |
f |
В'СВ dxdy |
(2.7) |
||
Матрица В имеет вид: |
—0,5а 0,*5b |
J |
|
|
||
|
|
|
|
|
||
o ü *L |
о |
2 |
6У |
|
||
СУ |
Еу |
|
Ех |
Ех |
|
|
6х |
2 |
0 |
J l h |
N |
« |
|
ЕУ |
ЕУ |
~ Ех |
||||
1 |
— Ъ |
* |
||||
0 |
0 |
0 |
0 |
|
||
G |
|
|||||
|
“ |
|
|
|
Заметим,что матрицы С и В входят в следующие уравнения связи:
Ох |
|
|
г /11 |
|
|
|
|
fl |
|
|
|
|
|
|
|
|
|
Т.\*у |
|
|
/з |
|
|
|
|
П |
|
|
|
|
|
|
L и _ |
|
|
Перемещение точек прямоугольного элемента из |
е г о |
п л о с |
|||
к о с т и при изгибе |
может быть задано в виде: |
|
|
||
™ (*» « /)= /i+ /a |
x + h У+ fi x * + h |
У2"h/ G XIJ+ |
|
||
H - h x 2y + |
f s x y 2 + |
f 9 л 3 + f i о y 2 + f n |
-v3 y - |- / i 2 |
.Vf/3 . |
( 2 . 8 ) |
Данный полином удовлетворяет однородному дифференциаль ному уравнению изгибаемой пластины. На каждый узел (угол) плас тинки накладываются связи — одна линейная (прогиб) w и две уг
ловые (углы |
поворота) ~ и |
Для получения м а т р и ц ы |
ж е с т к о с т и |
необходимо предварительно получить три матрицы |
А, В, и С. Матрица А получается в результате подстановки коорди нат узлов элемента в уравнение (2.8) и выражения ^ и ^ . Для каж
дого узла имеем 3 уравнения и порядок матрицы А будет равен 12. Матрица В входит в преобразование
^ д2 w ~ дх2
д2 w ду2
д2 w |
= [В] [fql при В= |
|
|
дхду |
|
d2w
дудх
“0 0 0 2 0 0 |
(2у) |
0 |
(6*) |
0 |
(6ху) |
0 |
|
|||
— |
0 0 0 0 2 0 |
0 |
(2*) |
0 |
(%) |
0 |
(6ху) |
|||
0 0 0 0 0 1 |
(2а) |
(2у) |
0 |
0 |
(За-2) |
(3ÿ2) |
||||
|
0 0 0 0 0 1 |
(2*) |
(2У) |
0 |
0 |
(За-2) |
(3</2) |
|||
Матрица С входит в преобразование |
|
|
|
|||||||
|
|
|
д2 w |
|
|
|
|
|
|
|
|
|
|
дх* |
|
|
|
|
|
|
|
Мх |
' |
|
д3 w |
|
~ |
Dx |
(Pît &х) |
0 |
0 - |
|
Му |
= |
—ICI |
ду2 |
при С = |
(Ру Dy) |
Dy |
0 |
0 |
||
Мху |
d2w |
0 |
0 |
Як |
0 |
|||||
|
|
|
|
|||||||
Мух |
|
|
дхду |
|
|
0 |
0 |
0 |
я« |
|
|
|
|
d2 w |
|
|
|
|
|
|
|
|
|
|
дудх |
|
|
|
|
|
|
|
|
Ах ~ |
Ex h3 |
|
Еьг& |
|
Gk3 |
|
|
|
|
|
12v У |
и |
12v 1 |
к == |
6 ;и * |
|
|
|
После введения этих выражений матрица жесткости изгибаемого прямоугольного элемента:
г —h (А |
В'СВ dxdy ]А- 1 . |
(2.9) |
Используя полученные выражения (2.7) и (2.9), можно дать раз вернутые формулы для матриц жесткости при плоском напряженном состоянии и изгибе пластины.
2.4. ОСОБЕННОСТИ ПРИМЕНЕНИЯ МЕТОДОВ РАСЧЕТА ПРОСТРАНСТВЕННЫХ СИСТЕМ К СТАЛЬНЫМ МОСТАМ
Существующие теории и методы расчета пространственных сис тем могут бытьэффективио использованы для расчета стальных мос тов. В зависимости от вида конструкции необходимо учитывать осо бенности применения того или иного метода расчета.
Рнс. 2.9. Основная система для кинематически неопределимого блока:
1 — внешние кинематические воздействия
Методы с т е р ж н е в ы х систем в канонической форме широ ко используются для расчета при проектировании мостов. В част ности, метод перемещений применен автором для пространствен ного расчета сквозных пролетных строений под железную дорогу. В них конструкция проезжей части (балочная клетка) посредством специальных диафрагм включена в совместную работу с главными фермами, что значительно повысило степень кинематической неоп ределимости системы.
При реализации метода расчета на ЭВМ целесообразно выделять из системы к и н е м а т и ч е с к и н е о п р е д е л и м ы е блоки, в качестве которых, в данном случае, можно принимать кон струкцию проезжей части, заключенную между двумя соседними по перечными балками.
Такой подход позволяет, с одной стороны, сократить число не известных, определяемых на одном этапе расчета, а с другой — стандартизировать расчет балочной клетки проезжей части, для ко торой вычисление коэффициентов уравнений имеет ряд тонкостей (учет эксцентриситетов, схем диафрагм).
Пример для кинематически неопределимого блока проезжей час ти показан на рис. 2.9. Выделены внешние кинематические воздей ствия и внешние силы по отношению к блоку. В качестве внешних кинематических воздействий выступают линейные и угловые пере мещения, которые являются неизвестными для пролетного строения в целом и должны быть определены. В расчете блока их величины принимают равными единице при определении реакций. Кинемати чески неопределимый блок имеет внутреннюю неопределимость (см. рис. 2.9), которая раскрывается методом перемещений.
После решения глобальной системы уравнений к блоку обраща ются повторно с целью определения усилий в элементах блока. В ре ализации метода на ЭВМ «вход» и «выход» блока оформляют так, что бы можно было вычислять как реакции блока на внешние воздейст вия, так и усилия во всех его элементах.
Рассмотренный подход к расчету сквозных пролетных строений реализован автором в программе СК ПО]. Опыт показал целесооб разность разработки специализированных программ для определен
ного класса конструкций при |
полной автоматизации расчета — |
|
получения расчетных усилий в элементах. |
|
|
Методы т о н к о с т е н н ы х |
с т е р ж н е й |
открытого и замк |
нутого профиля можно использовать для расчета |
сплошностенчатых |
конструкций пролетных строений, в том числе вантовых, а также для решения широкого круга вопросов пространственной работы соору жений. Предположение о недеформируемости контура поперечного сечения позволяет свести к минимуму число неизвестных интеграль ных факторов.
Вобычной постановке задачи, в методе перемещений неизвест ны в каждом узле три угловых, три линейных перемещения и депланация. Возможность учета переменности поперечного сечения по длине пролета, а также любого типа упругих связей в узлах делает метод тонкостенных стержней достаточно эффективным, особенно в случае рассмотрения моста как единой системы. Известна программ ма автора [47], реализующая данный метод на ЭВМ.
Врамках гипотезы о недеформируемости контура поперечного сечения можно решать вопросы неравномерного распределения на пряжений по ширине пролетного строения (случай симметричного загружения поперечного сечения). Для этого нужно разложить про
дольные деформации в ряд по координатным функциям и учесть сдви ги. Такой подход позволяет с единых позиций подойти к расчету как при изгибе, так и кручении любых типов сечений — открытых, замк нутых, комбинированных. Соответствующая программа для ЭВМ
разработана H. М. Митропольским |
под |
руководством автора Л |
Методы п л и т н о - б а л о ч н ы х |
конструкций существуют в |
|
различных -модификациях и используются |
как в форме метода сил, |
так и метода перемещений. В общем случае предпочтение, видимо*, нужно отдать методу перемещений, хотя имеются весьма положив тельные примеры использования и метода сил. Данный метод точен для широкого класса конструкций в рамках прикладной теории уп ругости и поэтому часто используется для расчета мостовых кон струкций. Пластины, входящие в состав расчетной схемы, могут быть изотропными или анизотропными, что важно для стальных конструкций с ребристыми плитами. Наиболее целесообразной мо дификацией данного метода для стальных мостов — вариант расчет ной схемы по типу «пластина — стержень—пластина», позволяющей1
1 Более детально с реализацией теории тонкостенных стержней на ЭВМ можно ознакомиться в литературе [47].
рассчитывать точно как коробчатые системы, так и пролетные строе ния открытого профиля.
Поскольку методом плитно-балочных конструкций решаются практически все вопросы пространственного расчета данного клас са конструкций, целесообразно распространить его и на сложные системы — неразрезные, с дискретными поперечными связями, слож ного очертания в плане. В последнем случае можно использовать метод расширения заданной системы [3].
Опыт практического применения метода плитно-балочных конструкций1, реализованного в ряде программ (МП-4 автора, Спика 10. М. Егорушкина и др.), оказался весьма положительным. В то же время есть возможности совершенствования рассматриваемо го метода и его возможности нельзя считать исчерпанными. В частно сти, можно учитывать переменность высоты балок по длине пролета путем компенсации этой переменности дополнительными нагрузка ми. Вообще, принимая разрезное пролетное строение в качестве основной системы можно рассчитывать весьма многообразные типы мостов.
Метод к о н е ч н ы х э л е м е н т о в существует в различных модификациях и нашел широкое распространение в расчетах самых разнообразных конструкций. Глобальная система уравнений мето да может быть получена как методом Бубнова-Галеркина, так и с использованием энергетических представлений. Потенциальная энергия системы для метода перемещений может быть записана в виде:
п |
т |
|
|
Я = 2 |
P t Z t - ' S , |
f WjdV, |
(2.10) |
*= 1 |
/=1 |
V |
|
где Р — внешние силы, число которых п\ Zt — искомые узловые |
перемет |
щения; т — число элементов, на которые разделена рассчитываемая система; Wj—внутренняя удельная потенциальная энергия в точке элемента объемом
Для определения неизвестных перемещений Z, имеем п линей ных уравнений
дП
— “ = 0 при i = l, 2, ... * л. |
(2. 11) |
Важный этап расчета — построение матрицы жесткостей эле мента в методе перемещений. Для случая анизотропной пластинки соответствующие формулы приведены в п. 2.3.
-При реализации метода часто оказывается необходимым исполь зовать так называемые суперэлементы (кинематически неопредели мые блоки), что позволяет уменьшить порядок глобальной системы уравнений. Для стальных пролетных строений в качестве суперэле мента можно принять блок на всю ширину поперечного сечения и длиной, определяемой особенностями конструкции, например, длн-
-1 В литературе даются практические примеры реализации метода [29, 42, 44].
7м
Рис. 2.10. Расчетная схема перемещении и усилий в опоре под пилон:
/ -i-210 — перемещения узлов. Перемещения в метрах |
даны увеличенными |
в Е раз» где |
|
Е — модуль упругости |
бетона в мегапаскалях; усилия Nrop и ААпер — в |
мегапыотонах |
|
ной, равной шагу |
поперечных связей. |
Тогда блоки «стыкуются» |
только по контуру поперечного сечения и матрица глобальной систе мы уравнений имеет ленточную структуру. При расчете протяжен ных систем, каковыми являются мосты, данный подход будет ес тественным.
Методы стержневых систем могут быть интерпретированы с по зиций метода конечных элементов и в этом случае он будет точным методом в рамках принятых кинематических гипотез. Более того, возможно представление сплошного тела в виде решетчатой (стерж невой) модели. Это расширяет область применения классических методов расчета стержневых систем. В качестве примера рассчитана опора под пилон вантового моста (рис. 2.10) с использованием про граммы ПС-1 [47]. Решена задача о плоском напряженном состоя нии пластины с переменной по высоте толщиной. Замена сплошного тела перекрестной стержневой системой в виде комплекса квадрат ных элементов позволила определить напряженное и деформиро ванное состояние опоры моста г.
Для реализации метода конечных элементов разработаны про граммы, из которых можно отметить вычислительные комплексы «СПРИНТ», «ЛИРА», «СУПЕР».1
1 Детали определения жесткостей элементов и учета коэффициента Пуас сона можно найти в литературе [47].
3. к о м б и н и р о в а н н ы й м е т о д п р и к л а д н о й
ДЕФОРМАЦИОННОЙ ТЕОРИИ ПЛАСТИЧНОСТИ
3.1.ОСОБЕННОСТИ ТЕОРИИ ПЛАСТИЧНОСТИ
ИМЕТОДА УПРУГИХ РЕШЕНИЙ
Наиболее развита в настоящее время теория малых деформаций упругопластических тел, что соответствует расчету конструкций с учетом ф и з и ч е с к о й н е л и н е й и о с т и. Полагают [21], что единственный дефект современной математической теории плас тичности — отсутствие доказательства теоремы существования ре шения основной краевой задачи, что, однако, не служит серьезным препятствием для широкого использования теории пластичности в технических приложениях.
Первая математическая теория пластического течения металлов была создана в XIX в. Б. Сен-Венаном и основывалась на гипотезе о пропорциональности девиаторов напряжений и скоростей дефор маций при условии текучести Треска (наибольших касательных на пряжений). За прошедший период были созданы различные вариантты теории пластичности, но как отмечает А. А. Ильюшин [18], все же не было создано по степени соответствия опыту сопоставимой с теорией упругости Коши-Навье. Для одних типов нагружения тел более правильными оказывались результаты расчетов по теории те чения, для других — по теории деформаций и т. д.
Основы общей математической теории пластичности разработаны А. А. Ильюшиным [18], но практическое использование затруднено; однако можно устанавливать границы обоснованного применения ее вариантов, а также найти правильные соотношения для ранее неисследованных процессов деформирования.
Анализом основных теорий пластичности — теории вязкоплас тических течений, малых упругопластических деформаций, идеаль ной пластичности [18] с позиций общей математической теории — показано, что теория малых упругопластических деформаций стро го справедлива для процессов простого нагружения и дает удовлет ворительные результаты для процессов, заметно отклоняющихся от простых нагружений. Главная п р о б л е м а т е о р и и плас тичности — установление соотношений между напряжениями и де формациями, для чего используется ряд гипотез и предположений. Наиболее общей является теория пластичности, основанная на те ории течения и рассматривающая скорости пластических деформа ций и напряжений. Она базируется на фундаментальном квазитермодинамическом постулате Друккера, гласящем [27], что, если элемент тела находится в каком-либо напряженном состоянии и под влиянием внешнего воздействия в нем возникают дополнительные напряжения, то как в процессе приложения дополнительных напря жений, так и за полный цикл приложения и снятия дополнительных напряжений, внешнее воздействие совершает положительную рабо
ту, т. е. из элемента, имеющего исходное напряженное состояние, нельзя выделить полезную энергию в процессе приложёния и сня тия дополнительных напряжений.
Для поликристаллических тел используется п р и н ц и п м а к с и м а л ь н о г о сопротивления пластическим деформированиям [27], по которому каждому приращению пластической деформации при данном напряженном состоянии аи соответствует приращение работы, равное или большее значения работы, совершаемой прира
щением деформации dz\f при любом другом напряженном состоя нии off, соответствующем области, ограниченной поверхностью те кучести или на самой поверхности текучести, т. е.
33
ДД {ач ~ а*1) <Чр) > °-
Принцип максимального сопротивления пластическому деформи рованию имеет два важных следствия: 1) поверхность течения в про странстве напряжений или обобщенных сил никогда не может быть вогнутой; 2) поверхности течения должны соответствовать уравне ния, связывающие напряжение и деформацию, т. е. вектор прира щения деформации в данной точке поверхности течения перпенди кулярен к поверхности в той же точке и пропорционален вектору внешней нормали (напряжению).
Кроме этого, в теориях пластичности используют предположе ние о том, что действительные деформации всегда можно предста вить в виде суммы упругих ге и пластических sp деформаций, при чем упругие подчиняются закону Гука.
Помимо теории течения, как уже было указано, существует так же деформационная теория пластичности Генки-Илыошина, ко торая не рассматривает скорости пластических деформаций и напря жений и в связи с этим имеет серьезные дефекты [19,21]. Однако не смотря на это, признается возможность ее использования в прак тических приложениях, если пути нагружения всех элементов кон струкции не слишком сильно отклоняются от пропорционального нагружения. При этом ограничения нагружений более слабые для случая материалов с сингулярными поверхностями текучести, т. е. когда на границе между упругой и пластической областями имеются ребра или угловые точки. При полностью произвольных путях на гружения деформационная теория, естественно, неприменима; одна ко следует учитывать, что для каждого элемента конструкции мож но заранее определить режим нагружения, и ограничение отпадет.
В. В. Новожиловым, Л. М. Качановым, В. Т. Койтером и другими исследователями показано, что при активном процессе нагружения уравнения деформационной теории (упругопластических деформат ций) совпадают с уравнениями нелинейно-упругого тела, что, с од ной стороны, ограничивает ее применение при сложных зигзагооб разных путях нагружения, а с другой — открывает возможности ис пользования теорем и методов нелинейной теории упругости .(имеет ся в виду физическая нелинейность).