книги / Несущая способность конструкций при повторных нагружениях
..pdfПредположим, |
что |
на |
первом |
||
этапе с ростом температуры |
общего |
||||
нагрева диска |
t0 (т) перепад |
темпе |
|||
ратур |
по радиусу |
tx (р) быстро воз |
|||
растает от нуля до максимума (ко |
|||||
торый |
достигается |
одновременно |
|||
с Q , |
после чего |
начинает убывать, |
|||
и температурное поле во всем диске |
|||||
выравнивается |
по уровню t*Q. |
Затем |
|||
следует стационарный |
режим |
опре |
|||
деленной длительности, который за |
|||||
канчивается охлаждением (для упро |
щения примем, что тепловые напря жения на последнем этапе пренебре жимо малы).
При анализе условий возникновения прогрессирующего раз рушения примем для простоты, что сопротивление деформирова нию материала диска на этапах нагрева (определяемое диаграм мой кратковременного деформирования) не зависит от темпера туры. При работе на стационарных режимах сопротивление де формированию характеризуется изохронной кривой [43], отвеча ющей суммарному времени прошедших циклов (рис. 5.10). Ее параметры — предел ползучести, предел длительной прочности — убывают с увеличением времени выдержки и также для упрощения считаются не зависящими от температуры.
Для расчета на приспособляемость (согласно методике, рас смотренной в § 22) диаграммы кратковременного и длительного деформирования аппроксимируются диаграммами идеального упругопластического тела. Для удобства предельное напряжение при стационарных режимах представим в виде произведения kost где (7S — предел текучести аппроксимированной кривой кратко временного деформирования; k — отношение предела ползучести к пределу текучести, зависящее здесь только от суммарной дли тельности стационарных режимов (TS = 7W, где N — число цик лов); в общем случае k может учитывать также влияние темпера туры на сопротивление кратковременному и длительному де формированию (пределы текучести и ползучести).
Введенные упрощения позволяют более отчетливо выявить особенности влияния длительности нагружения в задачах при способляемости. В рамках рассматриваемого метода при наличии необходимой информации можно легко от них отка заться.
Для расчета условий возникновения прогрессирующего раз рушения воспользуемся соотношением (5.30), учитывая, что для плоских дисков без лопаточной (контурной) нагрузки при темпе ратурном поле (5.36) наилучшую верхнюю оценку параметров предельного цикла дает механизм полного разрушения. Эта оценка обычно совпадает с точным решением [7].
132
Учитывая, что для плоского диска температурам (5.36) отве чают окружные термоупругие напряжения, равные [44]
Сфт = |
aEti (т) |
, |
(5.37) |
преобразуем условие (5.30) |
к виду |
|
|
|
|
- |
I min [4 - |
(I-Зр*) ]ф. |
(5.38) |
||
ГДе |
п |
СОо : |
=3gas |
1 |
/0 |
40’s |
|
о |
|
||||||
|
Ро |
|
|
yb2 * |
tl = Ч Ё ' |
|
|
Здесь при переходных режимах /е = |
k0 = 1, а при стационар |
||||||
ных |
режимах |
k < k0l значение |
k0 определяется |
по имеющимся |
данным для материала диска (с учетом принятого допуска на иеупругую деформацию).
Учитывая характер изменения температуры по радиусу, не трудно установить, что минимум подынтегральному выражению в правой части неравенства (5.38) доставляют следующие значе
ния k и функции t± (т): |
|
|
|
|
|||
в области диска |
1 > |
р > |
р* — при стационарном |
режиме работы |
|||
|
|
|
|
k < \ \ |
tx(т) = |
0; |
(5.39) |
в области |
диска |
0 < |
р с |
р;!; — на переходном режиме |
|||
|
|
|
|
k = 1; |
t\ (т) = |
i*. |
(5.40) |
Здесь Р* = |
j [ 1 — (1 — А) -|-] • |
|
|
|
Согласно соотношениям (5.39) и (5.40) скорости неупругой деформации в центральной части диска в предельном цикле ста новятся отличными от нуля в момент наибольшего перепада тем ператур tx (т), т. е. в процессе нагрева, когда напряжения от центробежной нагрузки в этой части диска суммируются с тепло
выми. |
В |
периферийной |
части |
|
|
|||||
диска |
они |
становятся |
отличными |
|
|
|||||
от нуля при стационарном режиме |
|
|
||||||||
работы, т. е. за счет ползучести. |
|
|
||||||||
С увеличением |
длительности |
ра |
|
|
||||||
боты, а также при снижении тем |
|
|
||||||||
пературных |
перепадов |
в |
период |
|
|
|||||
нагрева величина р.,. уменьшается, |
|
|
||||||||
при определенных |
условиях |
она |
|
|
||||||
становится равной нулю. В по |
|
|
||||||||
следнем случае |
минимум |
подын |
|
|
||||||
тегральной функции в соотношении |
|
|
||||||||
(5.38) |
отвечает |
одному |
(стацио |
0 |
0,25 0,50 0,75 р/р0. |
|||||
нарному) |
режиму |
работы |
для |
|||||||
|
|
|||||||||
всех точек диска, что соответствует Рис. |
5.11 |
133
возникновению в предельном цикле совместной деформации по лзучести.
На рис. 5.11 приведены результаты расчета условий односто роннего накопления деформации с использованием уравнения (5.38) и соотношений (5.39), (5.40) при k — 1,0; 0,75; 0,5, что отве чает возрастающей расчетной длительности работы конструкции. Штриховые линии на рисунке определяют условия возникновения знакопеременного течения, найденные на основе предложенной в § 22 для рассматриваемого случая аппроксимации свойств мате риала на различных этапах цикла (при тех же суммарных дли тельностях для ползучести).
Каждая из диаграмм приспособляемости на рис. 5.11 опре деляет такие значения параметров цикла, превышение которых при заданной длительности работы приводит к сравнительно бы строму увеличению накопленных деформаций или их размахов. Если точка, положение которой определяется параметрами рабо чего цикла, окажется вне области приспособляемости, построенной при k = 1, повторные нагружения сразу же (после стабилизации цикла) приведут к циклической (непрекращающейся) пластической деформации (знакопеременной, односторонней или их сочетанию). При расположении рабочей точки внутри данной области иеупругая деформация остается ограниченной (как правило, не превы шает существенно допуска, по которому проведена схематизация реальной диаграммы деформирования) и увеличивается лишь, начиная с такого числа циклов, при котором диаграмма приспо собляемости, построенная для соответствующей длительности (зна чения k), прошла бы через данную точку. В зависимости от того, на каком участке границы области приспособляемости окажется при этом рабочая точка (штриховая линия, сплошная наклонная или вертикальная), реализуется «смешанное» (вязкопластическое) знакопеременное течение, одностороннее накопление деформации с каждым циклом (также смешанного характера), либо (в послед нем случае) вязкое течение, т. е. совместная деформация ползу чести на этапе стационарной работы.
Заметим, что с ростом длительности выдержек (или темпера туры при стационарном режиме) протяженность вертикального участка диаграммы приспособляемости увеличивается (рис. 5.11), следовательно, расширяется диапазон значений параметров цикла, когда несущая способность конструкции ограничивается усло виями работы на стационарном режиме, реализацией своеобраз ного «предельного равновесия» в условиях ползучести.
Г л а в а 6
Некоторые задачи приспособляемости пластин
Поведение пластинок и оболочек за пределами упругости, их несущая способность представляют значительный интерес для многих областей техники. Расчету пластинок и оболочек по пре дельному равновесию посвящена довольно обширная литература. Необходимо отметить, что фундаментальные статическая и кинема тическая теоремы соответствующей теории были впервые сформу лированы и применены к расчету пластинок в СССР еще в 30-е годы А. А. Гвоздевым [5]. В дальнейшем ряд задач о несущей способ
ности |
пластинок |
был рассмотрен А. А. Ильюшиным [23], |
А. Р. |
Ржаницыным |
[45] и другими авторами. Обстоятельные об |
зоры работ, посвященных проблеме упругопластического состо яния пластинок и оболочек, даны в монографии Ходжа [81 ] и в статье Г. С. Шапиро [64].
Существенный интерес представляет также поведение пластин
иоболочек при повторных нагружениях. Первые решения нетри виальных задач приспособляемости пластинок были получены около 10 лет назад [7 ] с помощью приближенных (кинематического
истатического) методов, использующих, в частности, преобразо ванную формулировку теоремы Койтера (см. гл. 1). Применение
симплекс-метода (линейное программирование) к расчету условий прогрессирующего формоизменения круглых пластинок рассма тривалось впервые в статьях, обзор которых был дан в работе [7 ]. Задача приспособляемости формулировалась при этом на основе статической теоремы. Несколько позднее были разработаны соот ветствующие методы с использованием соотношений преобразо ванной кинематической теоремы [13]. В работе [58] к решению задач приспособляемости круглых пластинок был применен прин цип максимума Понтрягина. В работах, в которых использованы методы линейного программирования и принцип максимума, обобщенные переменные вводились в соответствии с процедурой, изложенной в гл. 2. Напомним, что их введение имеет смысл лишь при определении условий прогрессирующего разрушения; анализ условий знакопеременного течения выполняется непосред ственно в напряжениях, и применение каких-либо специальных математических методов здесь не требуется.
135
В этой главе приводятся решения некоторых, представля ющих практический интерес, задач о приспособляемости пластин. В то же время эти решения иллюстрируют применение ряда мето дов. Круглая пластина, подверженная повторным механическим и тепловым воздействиям, при простых условиях опираиия была рассмотрена в [7].
§ 26. Круглая пластинка, защемленная по краю
Рассмотрим условия возникновения прогрессирующего изгиба защемленной по краю пластинки, нагруженной равномерным
давлением |
(0 < |
р < |
р *) |
при |
переменном |
температурном поле *: |
|||
|
|
|
|
t (х, |
Q = |
/, (т) |
+ 1, (т) (с■ |
+ 1 Е2) , |
(6.1) |
где |
0 « |
|
(т) < <*; £ = |
£ |
(— 1 « £ < 1 ) . |
|
|
||
|
Соответствующие термоупругие напряжения в пластинке |
[40] |
|||||||
|
|
|
|
<т£>= |
= |
q ( 1 - 6£ — ЗС*), |
(6.2) |
||
где |
q |
aEt (т) . |
о < |
я < |
<7*- |
|
|
|
|
6 ( 1 - Ю ’ |
|
|
|
||||||
|
|
|
|
|
Воспользуемся результатами анализа предельного равновесия данной пластинки [81 ], принимая приращения пластических прогибов за цикл пропорциональными соответствующим ско ростям:
0 |
\ |
.1 + In* -) |
( |
о |
< |
р |
< |
± |
) ; |
(6 .3 ) |
Aw = Aw0 |
( 1 ----— |
|
|
|
|
|
|
|
|
|
Aw = —Aw0 -г -)п,р— |
( - |
^ |
< |
p |
« |
l ) |
, |
(6 .4 ) |
||
|
0 |
1НIn X |
|
|
|
|
|
|
|
|
D
где * = -£-, b — неизвестный пока радиус окружности, разделя
ющей пластинку на две области, в которых реализуются различные режимы течения.
С помощью соотношений типа
__ 1 |
d2w ф __ |
1 dw |
Xr~ ~ W |
Up* ’ |
W ТЩ> |
определим приращения составляющих кривизны:
Дх, |
Дw0 |
|
х |
|
Дх, = 0 |
( о « р < ^ - ) ; |
|
p { l + \ n x ) |
|||||
|
|
|
|
|||
Дхф = |
— Дх, |
|
Аю0 |
|
|
|
|
Я2 Ра (| + 1па-) |
( т < р < 1 ) ' |
||||
|
|
|
|
(6 .5 )
(6 .6)
(6.7)
* Здесь и в дальнейшем ось £ (г) направлена вверх (положительное направ лен;::* прогиба — вниз).
Приращения пластической де формации за цикл определяются из выражений
Аег = |
—£А ДхЛ= 0; |
Деф= —£Л Дхф |
||
|
|
|
|
(6.8) |
(для упрощения |
штрихи, |
означа- |
||
ющие пластическое |
деформирова |
|||
ние, здесь и далее опускаются). |
||||
Рассмотрим |
вначале програм Рис. 6.1 |
|||
му |
нагружения, |
при |
которой |
|
теплосмены воздействуют |
на пластинку при р = р* = const. |
|||
Условие (1.65), записанное для переменных тепловых напря |
||||
жений с учетом выражений (6.2)—(6.8), дает |
||||
при 0 < р с - |
|
|
|
сгijxAet-/o — о\рт Авф Oj |
(6.9) |
при — < Р < 1 |
|
O ijx A tJ jy ’ O — СГфX А б ф - { - о Гх А е г — 0 . |
(6 . 10) |
Во внутренней части пластины область, где термоупругие напряжения совершают догрузку, охватывает слои, соответству
ющие заштрихованным участк_ам |
эпюры |
|
напряжений |
при |
|||||
Офт =h 0 |
(рис. 6.1, |
Со = |
-у = 2 |
|
) - |
В |
ее |
наружной |
части |
в связи |
с тем, что |
<J $ |
= о^т» |
а Авф = |
—Авг, |
область догрузки |
|||
вообще |
отсутствует. |
|
|
|
|
|
|
|
|
Уравнение (1.67) в данном случае принимает вид |
|
||||||||
|
|
1 |
|
1 |
|
|
|
|
|
|
J р;| Awp dp = |
| М0A ^p dp— |
|
|
|||||
|
о |
|
о |
|
|
|
|
|
|
|
— A J a ^ d t |
J Дефр ф + |
-^ -| d (^ |
} |
|p=i. |
(6.11) |
|||
|
—1 |
|
о |
|
|
|
|
|
|
Пределы первого интеграла в правой части разбиваются в соот ветствии с выражениями (6.3), (6.4). Последний член представляет собой (с точностью до постоянного множителя) диссипацию энергии в кольцевом пластическом шарнире, который образуется в за делке, этот член соответствует выражению (1.70) *. Работа термо упругих напряжений в кольцевом шарнире (в заделке), где имеет место разрыв перемещения, равна нулю, так как угол поворота
* Здесь и в дальнейшем (если не будет специальной оговорки) температурная зависимость предела текучести для упрощения не учитывается.
137
получает приращение в моменты времени, когда эти напряжения отсутствуют.
После выполнения вычислений и некоторого округления усло вие прогрессирующего изгиба защемленной пластины приведем к виду
где р0 = |
11,26 |
— давление, отвечающее предельному равно |
весию |
[81 ]; |
<7° = о5/4 — интенсивность теплосмен, предельная |
по условию знакопеременного течения при р = 0.
На основании кинематической теоремы параметр х, входящий в решение, находится из условия минимума параметра р*. Это
дает уравнение |
|
|
|
Зх2— 2 In * — 5 = — 2 -|jp, |
(6.13) |
с учетом которого из выражения (6.12) получим |
|
|
|
• ^ = 0,535хг. |
(6.14) |
|
Ро |
|
Как следует из (6.13), с увеличением интенсивности теплосмен |
||
^в интервале 0 < |
< 1^ область пластины, в |
которой реали |
зуется распределение прогибов, определяемое выражением (6.3),
постепенно расширяется, при = 1 она распространяется на
всю пластину.
При достаточно большой интенсивности теплосмен программа нагружения р = р*, при которой давление сохраняется постоян ным и равным своему максимальному значению, перестает быть наиболее неблагоприятной для приспособляемости. Такая ситуация
возникает при у > 0,827, когда возникает область (0,827 < |
р < |
, в которой режим течения (6.6) сочетается с условием М |
< |
< 0 (рис. 6.2, [55]) и, следовательно, максимум произведения
Мфр Ахф отвечает нагрузке р = 0 (индексом р отмечен изгибающий момент от действия давления). Для определения условия прогрес сирующего разрушения с учетом этой особенности основное урав нение (1.67) следует переписать в таком виде:
1
J |
М0Дх,,р dp— /г j J |
а%х Ае,рр ф + - ^ | - | - |
|f)==1 = 0 , (6.15) |
где |
<ТфХ— объемлющее |
распределение суммарных напряжений |
от нагрузки и температурного поля. Пределы интеграла в левой части разбиваются в соответствии с изменением аналитических выражений для Ахф, Деф и aj,x.
Учет влияния изменений нагрузки относится к случаю, когда
138
-% >0,5, при этом значение |
м |
|||
\ . |
|
теперь |
(1+МРЯ2 |
|
— |
[определяемое |
~в---- |
||
уже не из уравнения (б. 13), |
|
|||
а |
из другого |
уравнения, |
0 |
|
получаемого |
путем |
мини |
|
|
мизации решения |
(6.15)] |
|
||
сравнительно |
мало |
отли |
|
|
чается от единицы: 0,827 < |
|
|||
< |
-j < 1. Поэтому фактиче- |
Рис. 6.2 |
ски условие приспособляемости изменяется по сравнению с условием (6.11) лишь незначительно (отличие порядка 3% в величине— ),
ипрактического значения в данном случае поправка не имеет. Рассмотрим решение аналогичной задачи, отличающейся тем,
что температурное поле характеризуется линейным (а не пара болическим) законом изменения температуры по толщине пластины
t (т , о = *0М + 0 1 W, о с tx (Т) < /*. |
(6.16) |
В этом случае термоупругие напряжения в защемленной пла стине приводятся к не изменяющимся вдоль радиуса изгибающим моментам [55]:
м№ = М $ = f o - S 1 = 7 ^ (° < <7< <7*). (617)
и при чисто изгибиом механизме разрушения обобщенные пере менные вводятся тривиальным способом.
Приведем решение данной задачи с использованием различ
ных методов. |
в строгой постановке (задача линейного |
||
Статический метод |
|||
программирования). |
Требуется |
определить при |
заданном </* |
|
max |
/)* = ? |
(6.18) |
|
м°г,м1 |
|
|
когда не зависящие от времени (здесь остаточные) изгибающие
моменты, связанные между собой уравнением |
равновесия |
•щ; (РЛ4?) — М% = 0, |
(6.19) |
принадлежат области, ограниченной фиктивной поверхностью взаимодействия. Эта область определяется неравенством
max j^max (| М°г + М^рХ+ |
М$т I, |
|
||
I м% + M !;J X + < у т |, |
I м ° — м% + |
м \% — /И(ррТ ^ |
м'$>— |
|
— М & |
|)] « M o , |
M o = |
o sh \ |
(6.20) |
соответствующим критерию текучести Треска (1.16).
139
Здесь индексами р, q отмечены изгибающие моменты от меха нической нагрузки и тепловые соответственно. Учитывая, что для принятых здесь диапазонов изменения этих величин в соот
ветствии с рис. |
6.2 |
для всех сечений |
пластинки |
|
||
|
|
|
min (М\%— МфрТ) = |
О, |
|
|
|
|
|
X |
|
|
|
а в соответствии |
с |
(6.17) |
|
|
|
|
|
— Mrqx — 0, min М$ х = |
min М ^х = О, |
|
|||
из условия |
(6.20) |
X |
|
X |
|
|
получим |
|
|
|
|||
— М0« М ф < Мо— т а х Ж ф ^ — max Мф^, |
|
|||||
|
|
|
X |
|
X |
|
— М0— min Мгр\ < М? < |
Л10— max М г%;{ |
|
||||
|
|
|
т |
|
X |
|
— Мо~ |
min (Mrfx— Л1фрТ) < |
М?— М° с Af0. |
(6.21) |
|||
Уравнение |
|
X |
представим в виде |
|
|
|
(6.19) |
|
|
||||
|
|
|
р |
|
|
|
|
|
|
рМ° — J M%d$ |
|
(6.22) |
|
|
|
|
О |
|
|
|
(постоянная интегрирования равна нулю, так как при р = 0 момент Mr ограничен). Правую часть уравнения (6.22) представим в виде конечной суммы
р Ж с = Д M%,dh |
(6.23) |
где dj = у (Дру_г + Др;) при / > 1 и / > i;
^ 1 = у APi; ^ = у Ар; при / = i; Др/ = Ру — р/_1-
Накладывая ограничения (6.21) лишь в узловых точках и исклю чая из них с помощью уравнения (6.23)j радиальные изгибающие моменты, после некоторых преобразований систему ограничений можно получить в виде, удобном для применения симплексметода. В качестве примера приведем эту систему для интервала
« р |
< ]/"-[§■ ( У |
^ |
«=«0,628, ] |
/ |
s==s0,827), в котором, |
как видно |
из рис. 6.2, |
Л4фрТ > 0, М |
< |
0. В результате не |
|
сложных |
преобразований |
получим |
|
|
|
|
|
|
1 — -Ц-р?) + 2 — ? > 0 ; |
||
|
|
2 |
X'di |
^ ° ’ |
|
140
|
|
7 r S * / d / |
+ |
p ( 1 - - r | p ? ) > 0; |
|
(6.24) |
||||
|
|
|
/=1 |
i |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
/^i |
|
|
|
|
|
|
|
|
|
2 xldi ~ |
Xi~ |
it W'+ 1» 0. |
|
||||
Здесь введены |
обозначения: |
|
|
|
|
|
|
|||
* = 1 |
+ |
< L ; |
p = |
M |
l i L k l l . |
~a= |
Jk_ |
(6.25) |
||
л |
1 |
^ |
Л4Л’ |
^ |
|
lfifTc/i2 |
’ |
4 |
8as ' |
В расчете, который выполнялся на ЭВМ, число узловых се чений было принято равным двенадцати (при неравномерном шаге, чтобы более точно учесть изменение знаков моментов от внешней нагрузки согласно рис. 6.2); размер матрицы системы ограничений был 13X61. Результаты решения будут рассмотрены в конце параграфа.
Кинематический метод в |
строгой постановке. |
Требуется |
определить |
|
|
min |
д|; = ? |
(6.26) |
А Е ф , Д е "
при ограничениях (1.28), (1.69), (1.57), (2.8), которые принимают
вид
1 1
|
J р Aw dp > |
J (М^ Axr + |
М%^ ДХф) dp; |
(6.27) |
||
|
о |
о |
|
|
|
|
d(Да;) |
= - Я*р Акф; |
|
= |
—Ахг; (Дда)р=1 = |
0; (6.28) |
|
3fT~ |
|
|||||
|
Ак, |
|
Мта |
Э Ф |
( « ) |
|
|
|
ЛИ? |
|
|||
|
|
|
|
|
||
|
Диф = |
|
дф[ « Ю |
Р а > 0. |
(6.29) |
|
|
|
|
дм1 |
|||
|
|
|
|
|
|
|
Здесь Ф(Л4?;714ф, ) = |
0 — фиктивная |
поверхность |
взаимодей |
ствия, определяемая условием (6.20) при сохранении в нем знака равенства. Отметим, что в дальнейшем — при дискретизации рассматриваемой задачи — разрывы перемещений будут заменены более или менее быстрым (в зависимости от принятой дискретной схемы), но непрерывным их изменением. Поэтому в соотношении (6.27) отсутствуют члены, учитывающие работу в разрывах.
141