![](/user_photo/_userpic.png)
книги из ГПНТБ / Гениев Г.А. Теория пластичности бетона и железобетона
.pdfо
Z -const
U -const
О
Рис. 82. |
Рис. 83. |
Рис. 84. |
. В
Рис. 86. |
Рис. 87. |
161
1 1 — 1 0 1 8
(рис. 84). Построение решения различно в зависимости от величины угла (3 в точке О. Если этот угол равен углу Р на отрезке ОА в точке О, то поле характеристик имеет вид, показанный на рис. 84. Если же характеристика z=
= const, проведенная из точки О при движении |
к ней |
вдоль ОВ, лежит внутри области АОВ (рис. 85), |
то эта |
область разбивается на две части — ВОА' и А'ОА. Пер вая из них будет находиться в условиях предыдущего слу чая при задании значений t и р на характеристике ОА'. Эти значения определяются решением для области А'ОА — вырожденного случая начальной характеристи ческой задачи.
Численные методы решения исходной системы урав нений (5.19). Наиболее простыми практическими метода ми решения конкретных задач являются приближенные методы построения полей характеристик, связанные с пе реходом к конечно-разностным соотношениям. При этом используются различные рассмотренные выше свойства характеристик. Эти методы нашли широкое применение в статике сыпучей среды, газовой динамике и т. п. В тео рии идеальной пластичности эти методы развиты в рабо тах В. В. Соколовского [100], Р. Хилла [113], В. Праге ра и П. Ходжа.
Ниже приведены некоторые приближенные решения исходной системы уравнений (5.19).
В дальнейшем удобно принять за неизвестные величи ны параметры %и (3. Значения искомых величин %, (3, х, у в точках пересечения т-й характеристики £=const и 1-й характеристики Tj= const обозначим через %т ,г, рто,г, хт,и
Ут,1• ЕСЛИ ВеЛИЧИНЫ %m—1 ,г> %т,1—Ь Рт—1,1, Рт,1—1 >
... известны, то для вычисления значений %т,г, рт ,г
справедливы формулы |
|
|
|
|
%т,1 |
' *1/’ |
(5.71) |
|
Рт,/ = 5 т ~ Т1г |
||
|
|
||
Для приближенного вычисления значений xm,i и ym,i |
|||
заменим дифференциальные |
уравнения |
характеристик |
|
(5.21) и (5.22) разностными уравнениями вида |
|||
9 т , I |
У т — 1 , 1 ~ { Х т , 1 |
Х т — \ ,/) |
X |
|
X tg ( Pm—U |
|
(5.72) |
|
|
|
|
У т , 1 |
У m , l - \ — { X m , l |
X m , l - 1) X |
X ^(Р*,Л_1 + Vm . l - 1)'
162
Решая эти уравнения относительно хт,и Ут,и получим следующие рекуррентные формулы:
|
Ут—1,1~~ УтЛ—1~Ь xm,l—1 |
( РтЛ—1 ~Ь Ут,1—l) ~ |
> |
||
Xm,l = |
|
|
|
|
|
|
|
хт—1,1 ( Pm—\,l |
|
(5.73) |
|
|
|
'*§ ( Pm—1,1 ~ |
Vm—l,l) |
||
|
|
|
|||
^ , 1 = |
Ут-1,1 + {Xm,l |
|
|
<5'74) |
|
Рассмотрим более подробно |
приближенные решения в |
||||
применении к краевым задачам. |
|
|
|||
1. |
Задача Коши. Разделим дугу АВ на малые отрез |
||||
ки (0, 0), (1, 1) ... |
(рис. 86). |
Значения %и р в узлах, |
бли |
||
жайших к дуге АВ, например в узле(т, m + 1), находим |
|||||
по условиям на характеристиках |
z.—const (5.21) и « = |
||||
= const |
(5.22): |
|
|
|
|
|
1 |
Pm,m+ 1 |
%m,m |
Pm.m ’ |
(5.75) |
|
|
|
|
|
Координаты узлов определяются по рекуррентным формулам (5.73) и (5.74). В остальных узлах сетка ли ний скольжения строится по схеме начальной характери стической задачи.
2. Начальная характеристическая задача. Отрезки линий скольжения ОА и ОВ делим на малые части (/, 0), (2, 0), ..., (m, 0)\ (0, 1), (0, 2), ..., (О, п). На основании первого свойства характеристик находим значения функ ций х и р в узле (m, п) (рис. 87):
%m,n %т,0 ' Ч п |
%0,0» |
|
(5.76) |
Рт.п = Рт,О “Ь Ро,п |
Ро,0‘ |
Координаты узловых точек |
находятся по формулам |
(5.73), (5.74).
В вырожденном случае известен угол р в вершине О. Разделим этот угол на малые части Ро,о, Ро,ь •••» Ро,п, при этом Ро,п будет углом между характеристиками О—п и ОА (рис. 88). Значение рТО)Т1 в узле ш, п находим по вто рому соотношению (5.76). Параметр %разрывен в точке О, и его нельзя непосредственно вычислить по первой формуле (5.76). Определим сначала %\,п в точках (/, п)
и * |
163 |
(см. рис. 88). Для точки (1, 0) значения хко и Pi,0 заданы; тогда известен параметр
|
( 5 '7 7 ) |
постоянный вдоль линии |
и —const, проходящей через |
(1, 0). Отсюда |
|
%1,п = |
2Д + Pi,я • |
В дальнейшем можно пользоваться формулой (5.76), заменяя значение хо,о через xi.n-
3. Смешанная задача. В общем случае смешанной за дачи (рис. 85) решение в области ОАА' строится как в вырожденном случае начальной характеристической за дачи. Рассмотрим область ОА'В (см. рис. 89). Разделим ОА ' на малые отрезки (/, 0) (2, 0) ... Методом последо вательных приближений находим значения искомых функций х и р в точке (1, 1). Точки (2, 1), (3, 1)... вычис ляются как в начальной характеристической задаче, точ ка (2, 2) опять находится путем последовательных при ближений.
Разрывные решения. Соотношения на линии разрыва напряжений. В некоторых случаях невозможно построить решения с непрерывным распределением напряжений. При этом существуют решения с разрывными напряже ниями, удовлетворяющие граничным условиям (такие разрывы называются сильными).
Метод С. А. Христиановича [116], примененный ра нее к определению уравнений характеристических линий, позволяет получить и уравнение линии разрыва. Для это го достаточно в приведенных выше соотношениях (5.20) приравнять нулю знаменатель каждой дроби. Выполнив это, найдем уравнение линии разрыва в виде
^ = t g ( p ± Y). |
(5.78) |
ах
Касательная к линии разрыва образует с осью х либо угол p-j-y, либо угол р—у. Это означает, что линия раз рыва может быть характеристической или огибающей характеристик.
Вдоль линии разрыва производные от t и р по коорди натам х, у, а также и производные от компонент напря жения ох, оу, %Ху по координатам обращаются в бесконеч ность. Иными словами, компоненты напряжений претер-
164
певают на линиях разрыва конечные скачки. Наличие линий разрыва в решении указывает на невозможность су ществования чисто пластического напряженного состоя ния с непрерывным распределением напряжений. Поми мо пластических областей будут существовать упругие зоны, вырождающиеся в пределе в линии разрыва.
Необходимым и достаточным условием существования линии разрыва является равенство нулю вдоль этой ли нии функционального определителя преобразования:
Д = D (*»У) |
________ sin2Y________ |
= q / 5 ym |
D (£ ,rj) |
cos (p + y) cos (P — y) |
5 | 5ц |
Для доказательства высказанного положения выра зим производные от компонент напряжения ох, ау, гху по координатам х, у через производные от х, у по перемен ным ц, например:
дах __ дах |
5 | |
. дах |
5ц |
|
дх |
5£ |
дх |
5ц |
дх |
Используя соотношения |
(5.18) |
и принимая во внима |
|||
ние зависимости (5.50) и (5.52), получим |
|
|
|||
дах = 2 |
(И -Г 0cos 2|3)- |
|
t |
dt |
|
дх |
|
sin 2у |
5ц |
дх |
|
(t + |
Т 0cos 2|3) |
|
+ t |
(5.80) |
|
sin 2у |
дх |
и аналогичные выражения для остальных производных. Внося затем в (5.80) частные производные
5 | = |
J _ |
jty |
и _3ц ^ |
____1 |
ду |
дх |
Д |
5ц |
дх |
А |
5£ ’ |
найдем |
|
|
|
|
fox |
{t + Т0cos 2(5) |
- t dg |
|
|
дх |
sin 2y |
5ц |
|
|
(t + |
Tocos 2P)- |
dy_ |
(5.81) |
|
5S |
||||
|
sin 2y |
|
Все члены правых частей этих соотношений (исключая
по условию ограничены. Следовательно, если одна
из производных обращается в бесконечность, то функцио нальный определитель преобразования Д(х, у) =0. Таким
165
образом, доказана необходимость этого условия для су-
„ |
3£ |
Зц |
3 | |
Зц |
ществования линии разрыва. Равенство------ ------------ - — |
||||
|
дх |
ду |
ду |
дх |
1
= — показывает, что это условие является также и дос
таточным.
Можно показать, что огибающая характеристик в плоскости х, у является линией разрыва. Рассмотрим, к примеру, семейство характеристик | = const, параметри
ческие уравнения |
которого |
даются |
функциями х = |
||||||
= х(1, ц), |
у= у(1, т]) |
при |
go=const. Для нахождения |
||||||
уравнения |
огибающей |
продифференцируем функции |
|||||||
х = х(1, г|), |
у = у{\, л) |
по £, считая ц функцией от §, и ре |
|||||||
зультат дифференцирования приравняем нулю: |
|
||||||||
|
|
дх |
. |
дх |
Зг| |
_ |
|
|
|
|
|
"ai |
+ |
~дц |
зГ “ |
’ |
|
(5.82) |
|
|
|
ду |
. |
ду |
Зц |
= Q |
|
|
|
|
|
|
|
|
|||||
|
|
3| |
^ |
3П Э6 |
|
|
|
|
|
Исключая |
, получим соотношение |
|
|||||||
|
|
А (х,у) |
дх ду |
|
дх |
ду |
= q |
(5.83) |
|
|
|
|
Зг] |
Зц |
3£ |
|
|||
|
|
|
|
|
|
которое показывает, что вдоль огибающей характеристик функциональный определитель преобразования А(х, у) = = 0, т. е. огибающая является линией разрыва. Аналогич ные результаты для обобщенного условия пластичности
получены В. В. Соколовским [100]. |
|
|
|
Возвращаясь к исследованию интегралов уравнений |
|
пластичности, рассмотрим частное |
решение системы |
|
уравнений (5.54), соответствующее случаям: |
||
1 ) |
^^const, т] = г(о = const; |
|
2 ) |
£ = £o = const, r]=^const. |
|
|
В первом случае одно семейство характеристик в плос |
|
кости х, у состоит из прямых z = const |
(5.59), второе на |
ходится путем интегрирования уравнения — = tg (p —у).
йх
Проведенное исследование показывает, что это реше ние справедливо лишь до линии разрыва, представляющей собой огибающую прямолинейных характеристик и опре деляемой уравнением*
* (l + - J ) + cos2 (р + у) Ф'(р) = 0. |
(5.84) |
1 6 6
Это уравнение получено по интегралу (5.59) уравнений пластичности.
Аналогичное уравнение легко может быть найдено и во втором случае, при g = £0= const.
Рис. 90.
При разрывах в напряжениях (так же как и в случае идеальной пластичности) выполняются простые соотно шения, вытекающие из уравнений равновесия и условия пластичности.
На рис. 90 представлены отрезок линии разрыва N и компоненты напряжений, приложенные к бесконечно малому элементу, лежащему на линии N. Значения ком
167
понент напряжений по разные стороны линии разрыва отличаются индексами «-f-» и «—». Имеем
o f — 07 |
I |
(5.85) |
т+ = х~. |
||
П П |
|
|
Разрыв возможен лишь в тангенциальной составляю щей 0(. Условие пластичности (5.10), справедливое по обе стороны линии разрыва, разрешим относительно ot; имеем
Поскольку N — линия разрыва, в формуле (5.86) для
значений a t и аТ нужно взять верхний и нижний знаки соответственно. Величина скачка в at будет равна:
= |
27V r„-T=„+(7-j + - f S ’.) (5.87) |
или, используя более удобное обозначение,
[ о,] = 4 У 2Т0о„ - r l + (73 + ф S * ). (5.88)
В силу условий непрерывности напряжений ап и хп на линии разрыва на последней должны выполняться сле дующие соотношения:
/2+ ^
Г /+ cos2 (p+ — ф) = ^ г- + *“ соз2 ( Г - <р); (5.89)
t+ sin 2 ( р+ — <р) = t sin 2 ( (J |
—■ф) |
или |
|
~ + /co s2 (p — ф) J = |
0; |
(5.90)
[t sin 2 ф — ср)] = 0.
Согласно этим соотношениям, линия разрыва N в каждой точке не является биссектрисой угла, образуемо го одноименными линиями скольжения (рис. 91) [39].
168
Переходя к исследованию исходной системы уравне ний (5.11), (5.12), отметим, что условие (5.11) отличает ся от условия (5.9) коэффициентом в фигурных скобках. Изменение этого коэффициента в зависимости от вида напряженного состояния и его аппроксимация прямой
линией |
представлены на |
рис. 92. |
На рис. 93 показаны |
кривая прочного сопротивления для плоской деформации при условии пластичности (5.9) — II и кривая прочного сопротивления (5.11), полученная в результате аппрок симации уточняющего коэффициента — /. Максимальное отклонение аппроксимирующей прямой
/*(о) |
tg а 1 |
(5.91) |
|
RcRv |
RcRP} Rc — Rv |
от истинной кривой в области смешанного напряженно го состояния растяжения-сжатия составляет приблизи тельно 10%, а в области сжатия не превышает 4% (см. рис. 92). Из полученных ниже решений было уста новлено, что проведенная аппроксимация в большом диа пазоне напряженных состояний обеспечивает получение результатов, близких к истинным.
Подставляя (5.91) в (5.11), получим окончательное
выражение условия |
пластичности при плоском |
дефор |
|||||
мированном состоянии (см. рис. 93): |
|
|
|||||
1 + Зта = |
( r f |
+ |
— ) |
(5 + |
»)т), |
(5.92) |
|
|
|
T l j |
|
|
|
||
|
|
|
|
|
|
||
6т1 |
■ i; = |
tg« |
/ |
371 |
(5.93) |
||
RcRp |
\ |
ЯсЯ, |
|||||
|
169
а — угол наклона аппроксимирующей прямой к оси т;
; Р |
O'! + < Т 2 |
(5.94) |
|
||
Т0 ’ |
|
|
Т0 = |
|
(5.95) |
|
|
|
Не останавливаясь на подробном анализе исходной |
||
системы уравнений (5.12) и (5.92), |
как это было сделано |
для системы уравнений (5.9) и (5.12), отметим только основные моменты.
Разрешающие уравнения для плоской деформации.
На основании (5.92) и (5.94) |
имеем |
|
|
|
|||
|
|
<*i = Р + |
1 |
|
|
(5.96) |
|
|
|
о2 = р - |
■Т0х |
|
|
||
|
|
) |
|
|
|
||
или |
|
|
|
|
|
|
|
|
\ |
_ Т0(1 + Зт2) |
_ |
S1 |
+ |
Т0т. |
(5.97) |
|
<*2 J |
6 (Н- цт) |
|
6 г 0 |
|||
|
|
|
|
|
|||
Напряжения |
на произвольной площадке аж, °V> |
||||||
в соответствии с (5.97) определяются соотношениями |
|||||||
о* |
J |
Г0(1+Зт») |
si |
|
, То т cos 2р; |
(5.98) |
|
tfy |
J |
6 (1 + ЦТ) |
6 т0 |
|
|
|
|
|
|
t ху ~ То Tsin2p, |
|
|
(5.99) |
где р — угол между положительным направлением оси х и направлением максимального главного нормального напряжения.
Подставляя (5.98) и (5.99) в систему дифференциаль ных уравнений равновесия (5.12), получим систему двух разрешающих дифференциальных уравнений относитель но неизвестных функций т и р :
(£ + cos 2р) — |
sin 2р — — 2т fsin 2р — — |
|||
дх |
ду |
|
\ |
дх |
— cos2p-^') = |
0; |
|
||
|
д\ |
+ |
2т |
(5Л0°) |
s i n 2 p ^ + ( £ — cos 2р) — |
cos2p — + |
|||
ОХ |
ду |
|
\ |
дх |
+ |
sin2p ^ ) |
= |
О, |
|
|
ду ! |
|
|
|
170