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

книги / Температурные напряжения и малоцикловая усталость

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

материала при заданной температуре. Для каждой выбранной долговечности размах деформации при пределе выносливости вычитают из размаха упругой деформации на рис. 4.15, а и соот­ ветствующую точку наносят на рис. 4.15, б. Эти точки (рис. 4.15, б) не могут образовать прямую линию, поскольку они получены вычитанием постоянной величины из линейного соотношения на рис. 4.15, а, тем более, что линейная зависимость для кривых усталости в двойных логарифмических координатах фактически точно не выполняется. Однако отклонение от линейности неве-

и Г*

0,100

А

0,010 __I __1

 

 

V

 

 

 

 

4/

а

 

Г

 

 

 

 

 

ЛяА

О

 

 

 

 

 

Ё _

1

 

Е

г

 

 

1

 

^----

 

ю’ 10* 10} 10* 105 Ю6циклов

0)

 

Рис. 4.15. Иллюстрация

 

метода перестроения кри­

0,001

вой усталости

в

упругих

деформациях

с

учетом

10' 10* 10* 10* 10* 10 6 циклов

предела выносливости

6)

 

 

 

лико, и полученные точки могут быть легко аппроксимированы прямой по методу наименьших квадратов (на рис. 4.15, б ли­ ния А'В').

3.

Размах

упругой деформации

реконструируется

исходя из

линии

А'В',

полученной на

рис. 4.15, б. Для

того чтобы опре­

делить ДееЬ

к ординатам А 'В '(Де*/

— 2е_1)

прибавляют 2е_ь

Перечерченная

кривая для

примера

показана линией

А’В'С’

на рис!

4.15, в.

Как уже отмечалось,

эта линия не может быть

прямой, однако она очень близка к ней и хорошо аппроксимирует исходную линию АВС на рис. 4.15, а (штриховая на рис. 4.15, в). Таким образом, кривая на рис. 4.15, в хорошо аппроксимирует зависимость для упругой составляющей и в то же время учиты­ вает такое важное свойство материала, как существование пре­ дела выносливости.

Рис. 4.15, в дает возможность записать уравнение кривой циклического деформирования/учитывающее упругое деформи­ рование при напряжениях ниже предела выносливости. Для заданных значений долговечности размах напряжений можно определить по кривой А’В'С' при известном значении модуля

упругости, а размах полной деформации — сложением

упругой

и пластической составляющих деформации (по А'В’С'

и ЭЕ)

(линии ИЕ на. рис. 4.15, в и а идентичны). Кривая

циклического

деформирования для 315,6° С дана на рис. 4.16 вместе с кривыми

для других температур,

определенных этим же

способом.

4.3.2.

Преобразование кривой

циклического

деформирования

в текущую

кривую

деформирования,

необходимую для

расчета

в общем случае. Кривая циклического деформирования связывает

размах напряжения с размахом деформации, но не характеризует

изменение напряжения деформации в любой период нагружения.

 

 

 

 

 

 

Вместе с тем, в общем слу­

 

 

 

 

 

 

чае необходима именно такая

 

 

 

 

 

 

информация и поэтому целе­

 

 

 

 

 

 

сообразно

разработать мето­

 

 

 

 

 

 

дику

перехода

от

 

кривой

 

 

 

 

 

 

циклического

 

деформирова­

 

 

 

 

 

 

ния в размахах к текущей

 

 

 

 

 

 

кривой циклического дефор­

 

 

 

 

 

 

мирования. Как и в разделе

 

 

 

 

 

 

4.2.2,

криваядеформирования

 

 

 

 

 

 

для каждого

участка

нагру­

 

 

 

 

 

 

жения

будет

определяться

 

 

 

 

 

 

для среднего значения темпе­

 

 

 

 

 

 

ратуры этого

участка тела.

 

 

 

 

 

 

 

При малом числе

циклов

 

 

 

 

 

 

нагружения

кривая

дефор­

 

 

 

 

 

 

мирования

в

 

любой

точке

 

 

 

 

 

 

может

изменяться

сложным

0

I

<■

6

9

ю&е,%

образом (для стали 4340 по-

казано

на

рис.

3.3).

Если,

Рис. 4.16. Кривые циклического дефор-

например, приложена цикли-

мироваиия для

нескольких значений тем-

ческая деформация

0,015, то

 

 

пературы

 

 

в первом цикле размах напря­

 

 

 

 

 

 

жений

оказывается

равным

267 кгс/мм2; в

последующих циклах

размах напряжений

посте­

пенно уменьшается пока примерно при 800 циклах не достиг­ нет 183 кгс/мм2. Этот размах напряжения (182 кгс/мм2) сохра­ няется на большей части долговечности, т. е. от 800 до 2800 циклов, когда материал разрушается. Если испытания при цикли­ ческом деформировании прервать после 100 циклов и построить кривую деформирования, то она будет отличаться от кривой, полу­ ченной после 800 циклов. Хотя в детали можно проследить распре­ деление напряжений и деформаций в каждом цикле с учетом изме­ нения кривой деформирования (предполагая, что это изменение известно или может быть определено), однако такое усложнение в большинстве инженерных расчетов, не оправдано. Поэтому целе­ сообразно построить текущую кривую деформирования для пе­ риода стабилизации размахов, который охватывает большую часть долговечности. Кроме того, как отмечалось выше, особен­ ности формы кривой деформирования не имеют существенного

192

значения при определении текущего размаха деформации, по­ скольку даже грубое предположение об инвариантности упругой деформации приводит к хорошим результатам и, таким образом, для большинства приложений пригоден приближенный подход.

В качестве примера рассмотрим рис. 4.17, которому соответ­ ствует кривая циклического деформирования, приведенная на рис.'4.16 для 315,6° С. Предположим, например, что после цикли­ ческого деформирования в течение 1000 циклов при постоянном размахе деформаций 0,015 статическая кривая деформирования

будет ОРВС1 на рис. 4.17.

 

 

Точка

В

определяет дефор­

 

 

мацию 0,0075 (0,015/2)

и на­

 

 

пряжение 91 кгс/мм2 (182/2);

 

 

эти значения

соответствуют

 

 

точке

диаграммы

цикличе­

 

 

ского

деформирования

при

 

 

размахе

деформации

0,015.

 

 

Если

кривая

ОРВС} неизве­

 

 

стна, то ее заменяют кривой

 

 

ОАВС, в которой ОА опре­

 

 

деляется модулем упругости,

 

 

а АВС является горизонталь­

 

 

ной линией, проходящей че­

 

 

рез известную точку В. Об­

 

 

ласти

сжатия

соответствует

Рис. 4.17. Замена фактической кривой

линия ОА'В'С', аналогичная

линии

ОАВС

при

растяже­

деформирования

кривой идеальной пла­

стичности, проходящей через одну точку,

нии.

Замена

эта

возможна

соответствующую

кривой циклического

лишь

в

том

случае,

когда

деформирования

размах

деформации

равен

 

 

0,015; если это другая величина, то предел текучести должен быть принят равным напряжению, соответствующему этому раз­ маху. Например, если размах деформации равен 0,045, то уровень напряжения на кривой деформирования по рис. 4.17 должен

быть принят

116 кгс/мм2

вместо 91 кгс/мм2, поскольку

размах

напряжения,

связанный

с размахом деформации 0,045,

равен

232 кгс/мм2.

 

 

 

Возникает вопрос относительно выбора кривой (деформирова­ ния, когда размах деформации, полученный при расчете по методу постоянства упругой деформации, меньше деформации, соответ­ ствующей пределу усталости. В этом случае горизонтальный участок кривой деформирования на рис. 4.17 может быть про­ веден на уровне напряжения, равного пределу усталости. Напри­ мер, если в положении О (рис. 4.8) размах деформации, полу­ ченный по методу постоянства упругой деформации, равен 0,002 (соответствующая амплитуда напряжения по кривой цикличе­ ского деформирования.равна 19 кгс/мм2) и предел усталости равен 27,1 кгс/мм2, то кривая деформирования имеет горизонтальный

участок на уровне 27,1 кгс/мм2. Далее будет построено решение без тех. ограничений, которые обусловливаются идеальной пла­ стичностью, хотя и принимается, что кривые деформирования

вкаждой точке не имеют участка упрочнения. Это связано не только с тем, что кривые деформирования отличаются в каждой точке, но и с тем, что при итерациях кривые деформирования пере­ страиваются таким образом, чтобы обеспечить размах деформации

вкаждой точке, соответствующий размаху напряжения по неидеа­ лизированной кривой циклического деформирования. Поскольку на кривой АВ рис. 4.17 представляет интерес только точка В, которая совпадает с точкой действительной диаграммы деформи­ рования ОРО., можно ожидать, что идеализация кривой на стадии, не входящей в окончательное решение, не будет иметь существен­ ного значения.

Следует также рассмотреть вопрос относительно построения кривой деформирования с учетом влияния двухосности напря­ женного состояния. Простейший метод состоит в следующем: кривую деформирования в точке выбирают после предваритель­ ной итерации, в процессе которой определяют компоненты дефор­

мации, при этом используют рекомендации раздела 3.5.2. Опре­

делив

разность

деформаций ех — е2,

е2 — е3 и е3 — б!

(где е1(

в2 и

в3 — три

главные

деформации)

и размахи

этих

разно­

стей

[ А (ех — е2), А (е2 — е3) и А (е3— ех) ], по уравнению (3.33)

находят размах эквивалентной деформации, по

методу

раз­

дела 3.5.2— размах напряжения. Половина этого

размаха

нап­

ряжения

принимается

в качестве предела текучести (АВС на

рис.

4.17).

Значение его изменяется от итерации к итерации,

по­

скольку размах деформаций в предыдущей итерации исполь­ зуется для определения кривой деформирования в последующей.

4.3.3. Определение долговечности. Определение долговечности при заданном размахе эквивалентной деформации рассмотрено в разделе 3.5.2. Строго говоря, до установления повторяющихся от цикла к циклу параметров нагружения (раздел 3.6) может возникнуть средняя деформация. Однако такая средняя дефор­ мация в типичных задачах термического или механического нагру­ жения является ничтожно малой и ее можно не учитывать. Следует также заметить, что рассмотренный в разделе 3.5 метод, относя­ щийся к циклической долговечности при сложном напряженном состоянии, основан на неподтвержденной гипотезе о возможности использовать данные, полученные при одноосном напряженном состоянии, и формулу Мизеса для пересчета на двухосное напря­ женное состояние, в то же время фактических данных, подтверж-. дающих справедливость этой гипотезы, получено еще недоста­ точно. Впрочем, как уже отмечалось, этот метод может быть легко заменен другим, если он окажется более обоснованным экспери­ ментально.

4.3.4. Решение на основе метода инвариантности упругой де­ формации. В этом случае упругое решение является предельно

194

простым, поскольку может быть получено для диска, например, по уравнениям (4.7) при допущении о равенстве нулю пластиче­ ских деформаций. Вычисления легко могут быть выполнены для каждого приращения нагрузки, указанного на рис. 4.13, путем подстановки температуры, скорости вращения и нагрузки на ободе в уравнения (4.7). Кольцевые деформации в центре диска и на ободе приводились на рис. 4.14. Результаты подробного решения даны на рис. 4.18, на котором указаны деформации и разности деформаций на ободе.

Рис. 4.18. Деформации на ободе, вычисленные по уравнениям упругости

4.3.5. Решение по методу циклической пластичности. Любому пластическому циклическому расчету предшествует определение размаха эквивалентной деформации по методу инвариантности упругой деформации (см. раздел 4.3.4). Если условия нагруже­ ния (см. рис. 4.13) подставить в уравнения (4.7), в которых все значения пластической деформации равны нулю, то 'в каждой точке диска могут быть определены три компонента упругой де­ формации. На рис. 4.18, а показаны результаты расчета для обода на рис. 4.8). Как отмечалось, определение отдельных деформа­ ций несущественно, значительно более важно знать разности ком­ понентов деформаций. Разности деформаций (ег — е0, ее — ег и е2 — ег) на ободе показаны на рис. 4.18, б. Эквивалентный раз­ мах деформации в каждой точке может быть определен по] размахам разностей деформации, для обода формула для вычисления эквивалентной деформации приведена на рис. 4.18, б. Особо отметим, что в расчете использованы размахи разностей деформа­ ций, а не наибольшие их значения. Например, для вычисления

соответствующего компонента А (е0 — ег) необходимо исполь­ зовать положения 6 и 2 приращения нагрузки, в то же время для получения компонента Д (ег — е,) приращения— 4 и 0. Как показано на рис. 4.18,6, размах эквивалентной деформации на ободе составляет 0,00615. Аналогично определен размах эквива­ лентной полной деформации в других точках, он показан штри­ ховой линией на рис. 4.19. Эти размахи деформаций и данные рис. 4.16 затем использованы для вычисления кривых деформи-

0

0,2

0,0

0 ,6

0,8

1,0

г

 

 

 

Рис.

4.19.

Распределение

эквивалентного

Рис. 4.20. Кривые деформн-

 

размаха деформаций в диске

 

рования

в

первом прибли­

рования

 

каждой точке,

 

схеме,

жении

для

циклического

В

ПО

упругопластического расчета

г

 

на

рис.

. , _ ’

 

круглой плоской плиты с ра-

приведеннои

4 .1 7 .

Кривые де-

ДИ7 ЛЬНЫМ

перепадом темпе-

формирования показаны на рис. 4.20,

ратур (для различных радиу-

на котором дана

кривая деформирова-

сов,

см. рис. 4.8)

ния лишь в первом квадранте. По оси ординат отложено отношение напряжений к модулю упругости, зна­

чение которого можно найти на рис. 4.9. Однако при переменном модуле упругости уравнения (4.7) получаются чрезвычайно гро­ моздкими, поэтому постоянный модуль упругости принимают и для кривой деформирования (ниже он всюду равен 1,97 X X 104 кгс/мм2).

Поскольку для каждой точки тела определяют единственную кривую деформирования, порядок вычислений оказывается до­ вольно простым, однако решение становится практически воз­ можным только при использовании быстродействующих вычисли­ тельных машин, так как при решении методом циклической пла­ стичности требуется выполнить огромный объем вычислений. При этом обработка исходных данных существенно отличается от той, которая связана с применением ручного счета. Например, кривая деформирования не задается графиком, а вычисляется. Размах напряжения определяется по кривой деформирования

196

графически по известному размаху деформации, при использо­ вании быстродействующих вычислительных машин эта задача решается одним из методов последовательных приближений, на­

пример,

методом

Нью­

 

тона—Рафсона.

 

Для

 

этого

в

машину

вво­

 

дится

вспомогательная

 

программа

для

опреде­

 

ления

размаха

напря­

 

жения по конкретному

 

размаху

 

деформации.

 

Аналогично

предусмат­

 

ривается вспомогатель­

 

ная

 

программа

для

 

спрямления

кривой ус­

 

талости

по

методу наи­

 

меньших квадратов (см.

 

раздел

4.3.1

 

и

рис.

 

4.15, в).

Различные де­

 

тали

 

счета,

известные

 

операторам

быстродей­

 

ствующих вычислитель­

 

ных машин,

здесь

под­

 

робно

 

не

 

рассматри­

 

ваются.

 

 

 

стадии,

 

Основные

 

 

входящие в процесс сче­

 

та,

можно

сгруппиро­

 

вать

 

следующим

обра­

 

зом:

 

 

 

 

 

 

 

 

1.

В машину вводят

 

основные данные, кото­

 

рые характеризуют ма­

 

териал

и

условия тер­

 

мического

и

механиче­

 

ского

нагружения.

 

2.

Определяют

кри­

 

вую деформирования в

 

каждой

точке и выпол­

 

няют расчет

по теории

 

пластических прираще­

 

ний

 

при

постепенном

Рис. 4.21. Блок-схема расчета по методу цикли­

нагружении и разгрузке

ческой пластичности

для

такого числа

цик­

 

лов, которое необходимо, чтобы последующие циклы не изменяли ни одну переменную, выбранную как критерий достижения уста­ новившегося состояния. Для описанного здесь метода в качестве критерия используют эквивалентную пластическую деформацию

за цикл, и установившееся состояние считается достигнутым, если разница в деформациях в двух соседних циклах во всех точках не

превышает

0,01%.

,

3.

Как

только

это состояние достигнуто, вычисляют размах

деформации

в

каждой

точке для последнего цикла нагружения.

По размахам деформации, определенным для установившегося состояния, вычисляют новые кривые деформирования в каждой точке и эти кривые сравнивают с кривыми, использованными ранее. Если они совпадают во всех точках с точностью 0,1%, то задача считается решенной. Если этого не происходит, то этап 2 повторяют. Для новых кривых деформирования вычисляют новые размахи деформаций до тех пор, пока кривые деформирования в двух последовательных приближениях не станут удовлетворять

указанному

значению критерия

сходимости.

 

 

рис. 4.21,

Блок-схема

указанных

операций

приведена на

а различные стадии вычислений

подробно описаны

в

табл. 4.1.

 

 

4.1. Программа основного

расчета по методу

 

 

 

 

 

 

циклической пластичности

 

 

 

1. Считывание данных:

 

 

 

 

 

 

 

 

 

а) числа приращений нагрузки за цикл;

 

 

 

 

 

б) распределения температур при каждом приращении;

 

 

в) центробежного слагаемого

в каждом приращении

1(1 — |х2)/8Е] раз2# 2;

г)

исходного предела текучести в каждой точке вдоль

радиуса;

зависимости

д)

постоянных,

определяющих

свойства

при

растяжении,

в

от средней температуры;

котором

выбран

предел

усталости

ъ

 

е)

числа циклов, при

напряжений

2.

Расчет

постоянных

для уравнения,

связывающего

размах

с размахом деформаций:

 

 

 

 

 

 

 

 

 

а) характеристик материала при растяжении в каждой точке по радиусу;

б)

коэффициентов и показателей степени

кривых усталости,

выраженных

в упругих и пластических деформациях:

Дел = (Л')7 и Дер = М (И )\

где Оуу» М и 2 рассчитаны по формулам (3.10), если имеются данные по усталости, или вычислены по формулам (3.25)—(3.28);

в) деформации, соответствующей пределу усталости в каждой точке Аб_ 1 =

= в/Е (А/-1)7;

 

 

 

 

 

что

г) новых постоянных X и ф методом наименьших квадратов, так

 

 

А&е1 — Дев1 =

X (Л ^

для

10 ^

^

106.

 

3.

Начало вычислений по новой кривой деформирования

 

 

 

Е/*р = 8

—В—0.

 

 

 

Эту стадию проходят в начале задачи, а также всякий раз, когда определяют

новые

кривые деформирования

для

установившихся размахов деформаций,

не согласующихся с принятыми пределами

текучести.

 

 

4.

- Начало

вычислений в

новом

цикле

 

 

 

 

 

е^р/цикл = размах

(ег — а Т) =

раЗМах

(ео — аТ) =

 

 

= размах

(е2 — а Т) = размах

(ег -—не) размах

(е2 — ег) =

О.

5. Приращение нагрузки:

а)

установление

распределения температур для приращения 'нагрузки;

б)

установление

скорости вращения для приращения нагрузки, *

6.Решение задачи пластичности в приращениях для конкретной нагрузки

инакопленных пластических деформаций;

а) для любого изменения нагрузки сначала предполагается, что

Ав/-р = Дв©р = 0;

б) решение уравнения (4.2) методом итераций для определения вг, е©, е2, ДеГр, Деор и применение способа, подобного приведенному в разделе 2.4.2;

в) в.записи новых условий нагружения заменяют исходную кривую дефор­ мирования

е* =

ег - Е

дЧ Р;

ее =

ео ~

2

д е0р‘»

е2 = ег- Е

ДЧр

(см. раздел 2.5.5.) и используют е^,

и

е* в уравнении (2.9).

7. Проверка сходимости по Де2Р во всех точках по критерию прекращения итераций при каждом приращении нагрузки

I Аб/р

I

Д б / р

 

Если сходимость не достигнута, то следует п. 6, б, если достигнуто, то п. 8. В обсуждаемой задаче величина критерия равна 0,0001, она может быть увели­

чена для

сокращения времени вычислений.

8. Для значений пластических деформаций в конце каждого приращения

нагрузки

вычисляют дополнительно значения, данные в пп. а—в1. В конце по­

следнего приращения нагрузки в цикле (возврат к начальному условию) опре­ деляют величины по п. г с рассмотрением всех значений Бицикле и вычислением разностей максимальных и минимальных величин:

 

 

 

я) Дв2Р =

— Двгр — Дб0 р,

 

 

е^/цикл =

е^/цикл +

Де,.р;

 

 

б)

^

Аегр == 2

^ ' ъгр +

&&гр\

 

 

 

2

Ае0р = 2

А'вор + Де0Р;

 

 

 

^

Де2Р =

^

Д'в2Р +

Де2р;

 

 

 

^

Дв/р =

2

^ 'в*р

в)

стг и <т0

[уравнение (2.8)]; а/, е/и е/р

[уравнения (2.3), (2.9), (2.14)1;

г)

размахн

вг— а 7\

в©

— а 7\

е2 — а 7\ гг — в©, в© — е2, е2— ег.

9.

Проверка завершения цикла

нагружения:

а) если цикл завершен, то дальше следует п. 10; б) если цикл не завершен, то следует повторить расчет п. 5

10. В

конце цикла

нагружения проверяют достижение установившегося

состояния,

используя в

качестве критерия е,-р/цикл

 

 

К - р / ц и к л — е | р / ц и к л |

 

 

е ,р / ц и к л

Если эта

разность

положительная,

то следует

п. 4.

Если

отрицательная,

то п.11.

 

 

 

 

 

 

 

 

11. Вычисляют новый предел текучести для наибольшего размаха в каждой

точке на радиусе.

 

 

 

 

 

 

и

«

-

2от

М /2 от

А

\ г / Ф

в случае, когда

Наибольший размах

равен —

I - ~ -

• Де_1 1

 

2а.

 

 

2о>р

то в этой точке обеспе­

он больше —2^

• Если наибольший размах меньше —

 

чено упругое состояние.

 

 

 

 

 

 

12. Проверка сходимости по от

 

 

 

 

 

Если эта разность положительная, то расчет повторяют с п.З. Если

отрица­

тельная, то далее следует п.13. В этой задаче критерий

принят равным

0,001.

13.

Запись

выведенного ответа.

относятся к

предыдущему приближе­

Замечание. Величины со штрихами

нию.

 

 

 

 

 

 

 

 

4.3.6.

Анализ

результатов.

Важный

результат расчетов при­

веден на рис. 4.19. Штриховой линией, как уже отмечалось, пока­ зан размах эквивалентной деформации, который получен из упругого расчета деформаций и использован для построения кри­ вой деформирования в каждой точке в исходном приближении. Сплошная кривая соответствует размаху эквивалентной деформа­ ции при расчете по методу циклической пластичности и показы­ вает отличие размахов эквивалентной деформации в этих двух случаях. Это отличие очень мало, и сходимость итераций по существу достигается после проведения шести приближений. Различия между промежуточными итерациями почти незаметны в сравнении с теми, что получены в шестом приближении. Малые изменения эквивалентного размаха деформаций от приближения к приближению означают, что кривые деформирования, опреде­ ленные по деформациям, полученным из упругого расчета, почти совпадают с точными. Однако, как будет показано в этом па­ раграфе, это не относится к составляющим деформаций в пласти­ ческой области; только эквивалентный размах деформаций, кото­ рый определяет кривую деформирования и циклическую долго­ вечность, практически не меняется. Также показано, что инва­ риантность результата имеет место лишь тогда, когда эквивалент­ ный размах деформации определяют рассматриваемым выше спо­ собом; при использовании других методов расчета необходимо проводить дополнительный анализ.

Результаты расчета показаны также на рис. 4.22 и 4.23, на которых приведены, соответственно, компоненты деформации и разности этих компонентов на ободе в зависимости от приращения нагрузки. Результаты упругого анализа, данные на рис. 4.18, приведены на рис. 4.22, а и 4.23, а для сравнения. На рис. 4.22, б—г показано изменение компонентов деформации в нескольких

200