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

книги / Несущая способность конструкций при повторных нагружениях

..pdf
Скачиваний:
4
Добавлен:
12.11.2023
Размер:
12.63 Mб
Скачать

Предположим,

что

на

первом

этапе с ростом температуры

общего

нагрева диска

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 и функции (т):

 

 

 

 

в области диска

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

Соседние файлы в папке книги