книги / Температурные напряжения и малоцикловая усталость
..pdfющее приближение прибавляет к решению один член из разло жения е* в степенной ряд. При достаточно большом числе прибли жений можно получить решение требуемой точности. Для дости жения инженерной точности может оказаться достаточным дватри приближения в зависимости от выбранного интервала для переменной х.
2.3.4. Сходимость. Критерием точности метода и окончания расчета является получение при последующих приближениях неизменяющегося решения. Графически понятие о сходимости иллюстрируется рис. 2.4. В данном случае рассмотрена простая
а) |
5) |
в) |
Рис. 2.4. Иллюстрация сходимости при ис пользовании метода последовательных при ближений для решения двух однородных линейных уравнений
задача — два совместных линейных алгебраических уравнений. Процесс приближений состоит из произвольного выбора значения для у, по которому для одного из линейных уравнений опреде ляется х, далее это значение х используют во втором уравнении для определения нового значения у, затем возвращаются к пер вому уравнению для получения нового значения х. Процесс по» вторяют поочередно с двумя уравнениями для получения новых значений х и у. Этот процесс при решении рассмотренных здесь уравнений очень прост, однако аналогичный подход применяется и тогда, когда должно быть решено очень большое число совмест ных уравнений.
Результаты расчетов показаны графически на рис. 2.4.
В случае а начальной точкой является А и результаты после довательных приближений указаны штриховой линией. На рис. 2.4 видно, что последовательные приближения сходятся к точке пересечения двух прямых, которая является точным решением. Случай б соответствует результатам последовательных приближе ний для другого направления прямых, описываемых линейными уравнениями. Здесь точное решение достигается по спирали вокруг
точки пересечения с постоянно перемещающимся радиусом. Слу чай в показывает ориентацию прямых, для которых последова тельные приближения не приводят к точному решению. Это слу чай, когда обе линии расположены под углом 45° к осям. После довательные приближения повторяют одни и те же значения, но не определяют точку пересечения.
Случай г показывает результаты при неблагоприятной ориен тации прямых, которая дает расходящееся решение. Как показано штриховой линией и стрелками, последовательные решения ухо дят все дальше от точного решения, изображаемого точкой пере сечения линий, т. е. решение расходится.
Для простых задач можно установить правила, которые опре деляют заранее, будет ли данный метод последовательных при ближений давать сходящееся или расходящееся решение. Однако в большинстве сложных случаев, встречающихся в задачах пол зучести и пластичности, разработать такие правила нелегко. В одних случаях последовательные приближения сходятся к ре шению; в других — решение расходится. Применение этого метода не вызывает опасений, так как процесс никогда не схо дится к неверному решению. Если приближения сходятся, то всегда к правильному решению, а расходимость быстро выяв ляется в процессе последовательных приближений. Обычно ока зывается возможным добиться сходимости решения за счет изме нения некоторых деталей процесса.
2.4.ПРИМЕНЕНИЕ МЕТОДА ПОСЛЕДОВАТЕЛЬНЫХ ПРИБЛИЖЕНИЙ К ЗАДАЧАМ ПЛАСТИЧНОСТИ
2.4.1.Приближения по напряжениям и деформациям одновре менно. Применение метода последовательных приближений к ре шению задач пластичности впервые описано в связи с конечно разностным расчетом дисков газовых турбин [2.8]. Решение незначительно отличается от описанного выше. Для усовершен ствования метода использована комбинация напряжений и дефор маций, полученных из первых приближений (рис. 2.5). Первый шаг — расчет упругих напряжений без учета пластических де формаций. Эквивалентные напряжения определяются по уравне нию (2.3). Если в какой-либо точке напряжение превышает предел пропорциональности материала при температуре в этой точке, то возникают пластические деформации и для определения
напряжений с учетом пластических деформаций необходимо ис пользовать метод последовательных приближений.
Допустим, например, что эквивалентное напряжение в п-й точке тела соответствует точке А на продолжении линейного участка кривой деформирования (см. рис. 2.5). Поскольку точка А расположена выше предела упругости (точка В), должны появиться пластические деформации. Напряжение и деформацию необходимо откорректировать так, чтобы они попали на реальную
92
диаграмму, характеризующую материал. Как и в начальной точке, полная деформация для реальных значений напряжений предполагается равной деформации в точке А. Значения напря* жения и деформации в данном состоянии соответствуют линии постоянной деформации, проходящей через А или С. Эквивалент ная пластическая деформация в п-й точке ерл дана отрезком СП. Значение егрл и е0рл можно
получить, |
используя величину |
|
|
||||||
ерл вместе с огп, о0л и а1п по |
|
|
|||||||
упругому |
расчету |
и применяя |
|
|
|||||
уравнения (2.6). |
|
|
|
|
|
||||
|
Используя уравнения равно |
|
|
||||||
весия |
и |
совместности, |
можно |
|
|
||||
подсчитать новые значения агп, |
|
|
|||||||
а0п и аш. По уравнению совме |
^ |
|
|||||||
стности |
полная |
деформация |
|
||||||
принимается равной упругой де- |
$ |
|
|||||||
формации, соответствующей |
не- |
^ |
|
||||||
известным напряжениям, |
плюс |
^ |
|
||||||
пластическая деформация, опре- |
$ |
|
|||||||
деленная |
при |
первом прибли |
|
|
|||||
жении. Новое значение а1ппре |
|
|
|||||||
вышает |
значение |
напряжения |
|
|
|||||
в точке П, например, его вели |
|
|
|||||||
чина |
может |
соответствовать |
|
|
|||||
точке |
Е. |
Хотя |
напряжения |
|
|
||||
в точке Е и деформация |
СП |
|
|
||||||
удовлетворяют условиям равно |
|
|
|||||||
весия и совместности, они соот |
|
|
|||||||
ветствуют точке Р, |
не лежащей |
|
|
||||||
на |
кривой |
деформирования, |
р нс. 2.5. Кривая деформирования при |
||||||
в результате |
чего |
физические |
одноосном растяжении, иллюстрирую- |
||||||
условия, |
определяемые |
свой- |
щая процесс определения правиль- |
||||||
ствами материала, еще не удов- |
пых значени® пластических дефор- |
||||||||
|
|
|
т-г |
|
** |
^ |
|
мзции |
|
летворяются. Поскольку любое |
|
|
|||||||
значение ерл, меньшее, чем СП, |
|
чем в точке Е, СЮ является |
|||||||
дает обычно значения а1п больше, |
|||||||||
нижним пределом ерл. Величина а1п, |
рассчитанная по деформации |
||||||||
ерл, |
соответствующей СИ, слишком велика, и прирост деформа |
ции ЕС, соответствующий этому,а1п, равен верхнему пределу ерл. Истинная величина ерл лежит между СП и ЕО, поэтому среднее значение между СП и ЕО, показанное, как НК, будет по пред положению хорошим приближением.
Новые значения сгп, ст0л и а1п можно рассчитать с исполь зованием величины НК, ъм, а также напряжений в точке Е а[п и агЛ и а0л, полученных в предыдущем приближении. Предположим, что это новое значение а.п лежит в точке М. Так как напряжение в точке М выше, чем в точке Н, приращение НК- величины ерл
становится новым нижним пределом врп. Далее, так как напря жение в точке М. меньше, чем в точке Е, соответствующее прира щение МЫ есть новый верхний предел ер„, "и процесс может быть повторен опять со средним значением, лежащим между МЫ и НК. Аналогично, если расчет по деформации врп, равной НК, дает а1п В точке Р, то НК будет являться новым верхним пределом и РС} — новым нижним пределом. Если же величина о1п соответствует точке Р, то НК станет новым верхним пределом врю но- Сй остается нижним пределом. Процесс повторяют до тех пор, пока
Е,гЮ3
о о,г о,о о,б о,в г/н
<$1 ,кгс/ммг________________ ______________________________
/Г_т |
--- — |
» - |
■— |
|
|
1 |
1с. |
0,012 |
0,016 |
е |
|
О |
0ОО0 0,000 |
|
|||
|
5) |
|
|
|
|
Рис. 2.6. Распределение температур |
(а), кривая деформирования (б) |
||||
и эквивалентной полной деформации |
(а) для закаливаемого длинного |
цилиндра
значение врп, использованное в расчете, и значение врп, соответ ствующее вычисленному о1п, не окажутся равными.
2.4.2. Приближения по деформациям. Форма кривой деформи рования в пластической области такова, что небольшие изменения напряжений вызывают существенные изменения деформации. Поэтому в ряде случаев желательно не вводить напряжения, а учи тывать только деформации. Остановимся на методе, который не посредственно вытекает из шести ранее перечисленных пунктов решения методом последовательных приближений и связан с при ближениями по деформациям. Для пояснения рассмотрим нагре тый длинный сплошной цилиндр с радиусом Р, хотя применение метода эффективно и во многих других случаях [2.61. Определе ние распределения температур как функции времени является, задачей теплопроводности и здесь не разбирается. Предположим, что температура в любой момент известна и надо определить рас пределение напряжений и деформаций. На рис. 2.6 показано распределение температур и кривая деформирования материала. Для простоты предположим, что эта диаграмма справедлива при
94
всех температурах; применение метода для случая, когда кривая деформирования зависит от температуры, рассмотрено ниже.
Полную деформацию выразим через пластические деформации. Как всегда, используем уравнение равновесия, совместности и кривую деформирования. Задача является осесимметричной, по этому уравнения равновесия и совместности представлены соот ношениями (1.18) и (1.20) при толщине, равной единице (анализ сделан для слоя единичной толщины по центру цилиндра), и © = = 0. Кроме того, используем условия плоской деформации — равенство осевого перемещения по всему сечению, за исключением
торца, е2 = соп$1:. Пусть цилиндр свободен по концам, так что
я
суммарное осевое усилие по сечению | а/йг равно нулю.
|
|
|
|
|
* |
о |
|
|
|
Эти условия приводят к соотношениям |
|
|
|||||||
8° = Т ^ |Г "7=~ I аТг йг + |
2(1-и ) |
“ Т 1 |
+ |
®еЛ гйг + |
|||||
|
о |
|
|
|
|
|
о |
|
|
|
, |
|
1 |
2р |
г |
|
* . . |
|
|
|
|
Г е/-р — ееР |
|
||||||
|
+ |
2 |
(1 - | 1) |
.) |
г |
йГ+ |
|
|
|
|
|
|
|
|
о |
|
|
|
|
6Г = |
— 8, |
|
1+1* |
а Т |
1— 2ц- |
|
|
||
|
1—ц |
егр + |
|
||||||
|
|
’0 ' |
1— 1* |
|
ГР |
|
|||
|
+ |
^ |
|
Г |
^ |
, |
г + г С . |
|
|
|
|
^ |
|
||||||
&г~ |
Я2 |
|
I аТг Лг— \ (егр + е6р)гйг |
(2.19) |
|||||
о |
|
о |
|
|
|
||||
|
|
|
|
о |
|
|
о |
|
|
|
|
|
|
|
|
я |
|
|
|
|
+ |
т ( г = й ' - к \ ( е'р + е°р>гс1г- |
|
Соотношения (2.19) решаются последовательными приближе ниями. Нулевое приближение полных деформаций получаем, предполагая, что егр и е0р равны нулю. Тогда полную эквивалент ную деформацию вычисляем по уравнениям (2.9), эквивалентную пластическую деформацию определяем по рис. 2.6, в и первые приближения для егр и е0Р — по уравнению (2.11). Эти величины егр и е0р подставляем в (2.19) и получаем новое приближение для ег, е0 и е2. Процесс повторяем столько раз, сколько необходимо
для получения желаемой степени сходимости. Напряжения можно вычислить по основным зависимостям для напряжений и дефор маций (2.8). Расчеты для первого приближения показаны в табл.
2.2.
2.2. Расчет деформаций в длинном сплошном цилиндре
п |
гЩ |
аГЛО» |
|
|
евР |
ег Л03 |
е0 .1О» |
||
|
% |
ур. |
(2.19) |
ур. (2.19) |
|||||
1 |
0 |
9,5 |
|
0 |
0 |
9,457 |
9,457 |
||
2 |
0,75 |
9,5 |
- |
0 |
'0 |
9,462 |
9,452 |
||
3 |
0,80 |
9,29 |
|
0 |
0 |
9,085 |
9,439 |
||
4 |
0,85 |
8,81 |
|
0 |
0 |
8,240 |
9,393 |
||
5 |
0,90 |
7,68 |
|
0 |
0 |
6,236 |
9,299 |
||
6 |
0,95 |
7,06 |
|
0 |
0 |
5,274 |
9,109 |
||
7 |
1,0 |
6,83 |
|
0 |
0 |
5,045 |
8,911 |
||
л |
е2-Ю3 |
ег Ю3 |
|
«в»-10' |
V |
I 0:, |
|
||
УР. (2.19) |
УР- |
(2.19) |
ур. (2.11) |
||||||
|
(рис. 2.7, о) |
ур. |
(2.11) |
||||||
1 |
8,912 |
0,363 |
|
0 |
|
0 |
0 |
||
2 |
0 |
0,363 |
|
0 |
|
0 |
0 |
||
3 |
0 |
0,310 |
|
0 |
|
0 |
0 |
||
4 |
0 |
0,668 |
|
0 |
|
0 |
0 |
||
5 |
0 |
1,926 |
|
0,91 |
—0,907 |
0,545 |
|||
6 |
0 |
2,494 |
|
1,40 |
—1,397 |
0,754 |
|||
7 |
0 |
2,578 |
|
1,48 |
—1,477 |
0,738 |
|||
♦ |
Первое приближение. |
|
|
|
|
|
|
Результаты различных приближений представлены на рис. 2.7, на котором видно, что для получения решения необходимо отно сительно немного приближений (сплошная линия — нулевое при ближение, штриховая — первое, штрих-пунктирная — четвертое
ипятое).
2.4.3.Быстрота сходимости. Проблема сходимости имеет боль шое значение в методе последовательных приближений этого типа. Для решения, которое получается из уравнений (2.19), можно показать, что условия сходимости удовлетворяются; однако точное решение можно получить при достаточно большом числе прибли жений [2.41. При наличии быстродействующих ЭЦВМ большой объем вычислений не может служить существенным ограничением. Если необходимо максимально уменьшить объем необходимых расчетов, можно применить по крайней мере два способа: 1) на чинать в первом приближении с распределения напряжений или
деформаций, как можно более близкого к точному решению; 2) вводить улучшения в полученные решения в процессе после довательных приближений, не полагаясь только на первое при ближение.
Предположение о том, что распределение упругих деформаций является правильным, приводит к хорошим результатам в боль шинстве задач с температурными нагрузками, поскольку в них окончательное распределение оказывается мало отличающимся от первоначально принятого [2.61. Если приложена внешняя нагрузка или возникнет локализованная пластическая деформа ция, когда остальной объем тела является упругим, предположе ние о сохранении упругого распределения деформаций не может
Рис. 2.7. Распределение напряжений (а) и деформаций (б) в сплошном цилиндре после различного числа приближений
[2.6 ]
служить наилучшим исходным приближением. В этом случае для получения лучшего исходного приближения можно исполь зовать упрощенный анализ, экспериментальные данные или сочетание предварительного расчета и эксперимента. Полезным может оказаться предварительный расчет, в котором учет пласти ческих деформаций ограничен лишь малой областью. При этом приближения сходятся быстро, хотя точность решения ограничена. Однако это решение может служить хорошим исходным прибли жением для решения, охватывающего большую пластическую область и более совершенного, чем в случае применения упругого решения в качестве исходного приближения.
Второй способ.можно применять, если рассматривать скорость изменения деформации в данной точке от приближения к при ближению с экстраполяцией до скорости, равной нулю. Например, можно делать три или четыре приближения, вычерчивая для каждой точки график приращения деформаций за приближение в зависимости от принятой деформации и экстраполируя его до деформации, обеспечивающей приращение, равное нулю. Поле ченное распределение деформаций можно использовать как новое исходное приближение для выполнения нескольких дополнитель ных приближений. Этим способом-в работе [2.61 решена задача
о вращающемся диске за 12 приближений с такой же точностью, которая получается за 40 приближений обычным методом.
2.4.4. Решение уравнений другими методами. В рассмотренном примере уравнения удалось проинтегрировать непосредственно, так как члены, содержащие пластическую деформацию, рассма тривались как известные величины. В ряде случаев непосред ственное интегрирование невозможно, и тогда после оценки пла стических деформаций из предшествовавшего приближения при ходится использовать метод конечных разностей* коллокации или другие приближенные приемы решения дифференциальных уравнений. При этом могут быть использованы многие приемы решения задач термоупругости, описанные выше, причем пласти ческие деформации, определяемые в каждом приближении, ■ должны входить в соотношения между напряжениями и деформа циями.
2.5. УПРУГОПЛАСТИЧЕСКИЙ РАСЧЕТ С ИСПОЛЬЗОВАНИЕМ ТЕОРИИ ТЕЧЕНИЯТЕОРИЯ ТЕЧЕНИЯ КАК РАЗВИТИЕ ДЕФОРМАЦИОННОЙ ТЕОРИИ
2.5.1. Ограничения деформационной теории. Как отмечалось, деформационная теория использует соотношение между напря жениями и деформациями. Во многих приложениях применение этой теории может быть эффективным, хотя ее исходные положе ния не всегда справедливы. Большинство проблем возникает в тех случаях, когда происходит «разгрузка», что имеет место при уменьшении напряжения в какой-либо точке тела, нагружен ной за пределами упругости; при разгрузке напряжения изме няются в пределах упругости, а пластическая деформация сохра-. няется, так что линейное соответствие между напряжением и де формацией не выполняется. Окончательную'Величину деформа ции можно определить при рассмотрении сложной истории на гружения, а не окончательного состояния. Известные трудности могут появиться также, если в процессе нагружения изменяется отношение главных напряжений. Таким образом, теория течения должна быть использована преимущественно тогда, когда усло вия применимости деформационной теории не выполняются, т. е.
не осуществляется простое нагружение.
Ниже рассмотрено применение этой теории к задачам о тем пературных напряжениях. Во многих случаях термического на гружения возникает нестационарное распределение температур, вызывающее в некоторых областях тела пластические деформации с разгрузкой до достижения окончательного распределения тем ператур. При стационарном распределении температур процессы ползучести постоянно изменяют распределение напряжений, отно шения главных напряжений, иногда также приводят к разгрузке. В- этом случае следует использовать метод приращения.
В следующих разделах изложены основы расчета методом при ращений, представленным как развитие деформационной теории. При этом оказывается возможным не только достичь лучшего понимания сложных процессов деформирования, но и получить конкретные решения, основанные на несколько ограниченной теории.
2.5.2. Поведение материала при одноосном напряжении. Для лучшего понимания основных уравнений теории течения жела тельно рассмотреть поведение материала при нагружении и раз грузке после однократного, по крайней мере, превышения пре дела текучести. Рассмотрим вйачале одноосное напряженное со стояние, распространив затем этот анализ на сложное состояние.
|
На рис. 2.8 линией ОАВСИР показана |
|
|
||||||
кривая деформирования материала. Если |
|
|
|||||||
при нагружении или разгрузке не превы |
|
|
|||||||
шается предел текучести А, то материал |
|
|
|||||||
остается упругим, и все процессы идут по |
|
|
|||||||
линии ОА\ полная разгрузка всегда воз |
|
|
|||||||
вращает материал к его -начальному со |
|
|
|||||||
стоянию О, независимо от истории нагру |
|
|
|||||||
жения или разгрузки. Если, однако, ак |
|
|
|||||||
тивное |
|
нагружение |
осуществляется |
до |
|
|
|||
точки Д |
и затем происходит разгрузка, |
|
|
||||||
то |
она |
осуществляется |
не |
по ДСВЛО, |
|
|
|||
а |
по |
линии ДД"Д'" |
Если |
напряжение |
|
|
|||
снижается до Д " и затем |
вновь увеличи |
|
|
||||||
вается, то нагружение идет по линии |
|
|
|||||||
Д"ДЕ; |
материал приобретает новые свой |
Рис. 2.8. Поведение ма |
|||||||
ства: предел текучести (Д) |
и кривую де |
териала |
при нагружении |
||||||
формирования Д"'Д"ДЕ (фактически вбли |
и |
разгрузке |
|||||||
зи Д происходит незначительное откло |
|
|
|||||||
нение |
этой кривой, и кривая несколько отличается от линии Д/7, |
||||||||
этот эффект обычно |
не учитывается в упругопластических рас |
||||||||
четах). |
Кроме того, |
предполагается, |
что |
материал при сжатии |
и растяжении имеет одинаковую кривую деформирования, и ее из менение вследствие пластических деформаций в этих случаях одинаково.
Таким образом, линия 0"М""Д"" изображает исходную кривую деформирования при сжатии (начало координат сдвинуто в точку О"", так что ОД'" = Д'"0""), а линия Д "'Д ""Е""— кривую сжатия после предварительного растяжения или сжатия до напряжения Д. При одноосном нагружении после однократной пластической деформации ер, деформация при последующем дей ствии напряжения а независимо от его направления — дальней
шее |
нагружение (например, |
аР) |
или разгрузка |
(например, |
аг,') |
— составит: |
|
|
|
|
е = ер + |
- ^ + |
Дер, |
(2.20) |
где Дер — некоторая новая пластическая деформация. Произой дет ли новая пластическая деформация или нет, зависит от конеч
ного напряжения |
и новой |
кривой |
деформирования О В"В"'Р"" |
|||
Если |
напряжение |
увеличилось до |
аР, то |
приращение |
пластиче |
|
ской |
деформаций |
составит |
Дер; |
если |
напряжение |
снизилось |
до Б ", то пластической деформации не возникнет, даже если напряжение таково, что вызывает пластическую деформа цию в исходном состоянии. Когда напряжение достигает значения ар-" то деформация
2.5.3 Поведение материала при сложном напряженном состоя* нии. Этот подход, имеющий определенное значение при анализе задач с разгрузкой при одноосном нагружении, еще более суще ственен при рассмотрении случаев плоского и объемного напря женного состояния. При сложном напряженном состоянии кри терием наличия или отсутствия последующих пластических де формаций является характер изменения эквивалентного напря жения. Если пластические деформации в теле уже произошли, и конечное эквивалентное напряжение равно а0 (см. рис. 2.8), то зависимость между эквивалентным напряжением и эквива лентной деформацией изображается линией Б"'Б''Б Р . В том случае, когда изменение компонентов напряженного состояния привело к уменьшению эквивалентного напряжения, пластиче ская деформация не увеличивается: материал разгружается по линии Ой "Б"' и упругие деформации определяются законом Гука. Если же изменение компонентов увеличивает эквивалентное напряжение до о>, то пластическая деформация увеличивается на Дер. Компоненты пластической деформации, соответствующие этой эквивалентной пластической деформации, определяются соответственно компонентами ах, ау и ог напряженного состояния
в |
точке Р. Пусть Агхр — приращение |
пластической деформации |
|
в |
направлении х. Используя уравнение (2.6) для общего случая |
||
объемного напряженного состояния, |
получим |
|
|
|
Ле*р = 4"Г [° х ~~ 1 ^ |
+ а*)] • |
(2-2*) |
Компонент полной деформации в направлении х при дости жении точки Р может быть получен с учетом: 1) теплового рас ширения; 2) упругих деформаций; 3) всех предыдущих пластиче ских деформаций вследствие начальной истории нагружения до напряженного состояния в О; 4) приращения пластической де формации при нагружении материала от точки Б до точки Р. Таким образом, если все первоначальные пластические деформации в точке Б, получаемые при какой-либо сложной истории нагруже-
100