книги из ГПНТБ / Жаров Г.Г. Судовые высокотемпературные газотурбинные установки
.pdfГлава VIII
ТЕМПЕРАТУРНЫЕ НАПРЯЖЕНИЯ В О Х Л А Ж Д А Е М Ы Х У З Л А Х ТУРБИНЫ
§60. Особенности оценки прочности охлаждаемых узлов турбины
При создании новых конструкций охлаждаемых узлов высокотемпературных ГТУ вопросы термоупругости приоб ретают существенное значение. В первую очередь это следует от нести к проточной части газовых турбин.
Помимо температурных напряжений, возникающих вследствие неравномерности подвода тепла по профилю и высоте охлаждаемых узлов (лопаток, дисков) во времени, создается дополнительный источник напряжений, вызванный разностью температур между га зом и хладагентом. Для оболочковых лопаток разность температур является менее существенной, так как тонкая оболочка почти мгно венно прогревается, а для лопаток с внутренними охлаждающими каналами этим обстоятельством пренебрегать нельзя. При измене нии режимов работы влияние обоих факторов усиливается.
Неравномерность нагрева влечет за собой изменение физико-ме ханических свойств материала и возникновение дополнительных разностей температур, что сопровождается неодинаковым тепловым расширением отдельных частей конструкций. Помимо этого, одни частицы препятствуют свободному расширению других, что приво дит к появлению в теле температурных напряжений. Величину и характер температурных напряжений необходимо знать, так как при существующих уровнях и разностях температур в охлаждае мых узлах турбин они играют существенную роль. В сочетании с ме ханическими температурные напряжения могут вызывать появление трещин и разрушение конструкций.
Некоторые материалы при быстром появлении напряжений ста новятся хрупкими и не выдерживают теплового удара. При повто ряющихся изменениях температуры возникают термоусталостные явления в элементах турбин. При длительном воздействии темпе ратурных напряжений может появиться существенная пластическая деформация. Наличие пластических деформаций приводит к пере распределению напряжений по сечению лопатки. Эти напряжения можно определить методом переменных параметров упругости. При определенных условиях значительные деформации могут привести
к тому, что лопатки будут касаться корпуса, а знакопеременные пластические деформации повлекут за собой термоусталостные явле ния. Все это способствует интенсивному или полному разрушению конструкции. У лопаток оболочкового типа это приводит к выпу чиванию пластины, а у сплошных лопаток — к появлению трещин. Помимо возникновения температурных напряжений вследствие изме нения температуры, связанного с подводом тепла, дополнительные температурные напряжения могут создаваться в результате про цесса деформирования тела. В этом случае наблюдается образование и движение тепловых потоков внутри тела, возникновение в нем связанных упругих и тепловых волн, термоупругое рассеивание энер гии и др.
Таким |
образом, увеличение температуры |
газа перед |
турбиной |
||||
с целью |
повышения |
эффективности |
газотурбинной |
установки в |
|||
конечном |
итоге ведет |
к |
снижению |
прочностных |
характеристик |
||
материалов и появлению |
дополнительных |
напряжений, |
которые |
||||
могут стать причиной разрушения высоконапряженных |
элементов |
||||||
газовой турбины. |
|
|
|
|
|
|
Применение различных систем охлаждения узлов газовой турбины позволяет в некоторой степени сократить температурные напряжения путем снижения температуры узлов и выравнивания температурных полей благодаря рациональному выбору геометрии каналов охлаж дения, типа хладагента и других элементов. Однако все это создает дополнительные трудности при оценке температурных напряжений путем расчетов, которые связаны с определением пространствен ного характера изменения напряжений во времени. Кроме того, ре лаксация, возникающая вследствие текучести и ползучести мате риалов, в еще большей степени усложняет задачу, так как в этом слу чае создается перераспределение напряжений. К тому же измене ние температуры материала по лопатке вызывает изменение модуля упругости и коэффициента линейного расширенияНа основе ги потезы плоских сечений можно подсчитать напряжения в неравно мерно нагретой лопатке с учетом зависимости £ и а от температуры:
|
+ |
Е { |
|
|
(272) |
где М ц , |
Мг |
— суммарные |
изгибающие моменты |
от действия ЦБС |
|
|
|
и газовых |
сил; |
|
|
Е, |
a, |
t—модуль |
упругости, коэффициент |
линейного расши |
|
|
|
рения и температура в точке сечения с координа |
|||
|
|
тами і |
И Т]. |
|
Первая часть формулы определяет напряжения, вызванные дей ствиями внешних сил, вторая — напряжения, вызванные неравно мерностью распределения температур. За оси координат принимаются
главные центральные осп термоупругой жесткости. Начало коорди нат называется термоупругпм центром сечения. Положение термо упругого центра сечения и осей определяют аналогично положению центра тяжести сечения. Отличие состоит в том, что в интеграл по площади сечения входит модуль упругости, переменный по сечению. Формула дает хорошее совпадение с действительным распределением напряжений неравномерно нагретых лопаток, но требует сложных расчетов и применения ЭЦВМ.
Поскольку напряженное состояние охлаждаемых узлов турбины зависит от многих факторов, и прежде всего от температурных, це лесообразно рассмотреть решение задачи термоупругости.
§ 61. Основные уравнения |
|
термоупругости |
|
Термоупругость занимается вопросами |
равно |
весия тела как термодинамической системы, взаимодействие |
кото |
рой с окружающей средой заключается лишь в механической |
работе |
внешних сил и теплообмена. Температура тела может изменяться как в результате подвода тепла от внешнего источника, так и вслед ствие самого процесса деформирования. Связь между деформацией и температурой устанавливается с помощью законов термодинамики.
В общем случае постановка задачи термоупругости |
заключается |
|||||||||||||||
в определении |
шестнадцати |
функций |
координат |
и |
времени |
[30]: |
||||||||||
— |
шести |
компонентов |
тензора |
напряжения |
at/-; |
|
|
|
||||||||
— |
шести |
компонентов |
тензора |
деформации |
ги; |
|
|
|
||||||||
— |
трех компонентов вектора перемещения Uf, |
|
|
|
||||||||||||
— |
температуры |
Т. |
|
|
|
|
|
|
|
|
|
|
|
|
||
Эти |
функции |
должны |
удовлетворять: |
|
|
|
|
|
|
|||||||
— |
трем |
уравнениям |
|
движения |
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
о"./, / = |
Р"ї; |
|
|
|
|
|
|
(273) |
|
— |
уравнению теплопроводности |
|
|
|
|
|
|
|
||||||||
|
|
|
Тц |
1-Т'- |
|
Г п |
( 3 |
\ г Э |
Д Ит |
e'kk |
= 0, |
|
|
(274) |
||
где а = — |
|
коэффициент |
температуропроводности; |
|
|
|||||||||||
— |
шести |
соотношениям |
между |
напряжениями и |
деформациями |
|||||||||||
|
|
аи |
= 2цв4 / + |
[Ukk |
- |
(3% + |
2ц.) ат (Т - |
Т0)] |
&и; |
|
(275) |
|||||
— |
шести |
соотношениям между деформациями |
и |
перемещениями |
||||||||||||
|
|
|
|
|
г |
а = 4 - |
|
( |
" |
" |
• |
+ |
( |
2 7 |
6 ) |
при начальных и граничных условиях, заданных определенным обра зом.
В приведенных уравнениях для упрощения записи применяются индексные обозначения и правило суммирования по повторяю щимся индексам, принятые в тензорном анализе:
|
— оси х, у , г в декартовой системе координат обозначаются че |
|||||
рез |
ха |
(а — 1, 2, |
3) или xt {і = |
1, 2, 3); |
||
|
— |
вектор перемещения щ в упругом теле означает вектор с ком |
||||
понентами иг\ |
и 2 ; и3\ |
|
|
|||
|
— напряженное и деформированное состояние упругого тела |
|||||
определяются |
соответственно тензорами второго ранга оС1- и e[f (if — |
|||||
= |
1,2, 3). Символы Оц и &ц- означают |
величины с девятью компонен |
||||
тами. |
Например, |
Оц означает |
любой |
компонент матрицы |
ОцСГі3 а1 3
^21^22^23
— повторяющийся индекс означает суммирование от 1 до 3
аН = ^11 "!~ СТ22 ~Г °"зз'>
— символ Кронекера |
oL! |
обладает |
следующими свойствами: |
|
a t ! = |
1 (при |
/ = |
/); |
|
а ( / |
= |
0 (при |
i=hj); |
— дифференцирование по определенной координате обозначается запятой на уровне индексов с одновременным индексным обозна чением соответствующей координаты:
|
// |
_iffL- |
Р |
- |
ЭЧц . |
|
|
и |
Ч - дх; ' |
е'/.« |
- |
дхкдхі ' |
|
— |
коэффициенты |
Ляме |
|
|
|
|
|
X — |
|
v£ |
|
_ |
Е |
|
Л - |
(1 +v)(l -f- 2v)' |
^ ~ |
2(1 +v ) • |
||
От |
общей системы |
уравнений |
термоупругости можно перейти |
к двум уравнениям, описывающим в векторной форме при опреде
ленных начальных |
и граничных |
условиях изменение в пространстве |
|||||||
и во времени поля |
деформаций |
и температурного |
поля: |
|
|||||
r iA2 « + |
(A, -f- ц.) grad div и — (ЗА, + |
2\л) сст grad Т = ри" = |
|
||||||
= |
V*T — ~T' |
— Т°{ЗХ |
+ 2 ^ |
" т divа" = |
0. |
(277) |
|||
Вектор перемещения |
« можно |
разложить на потенциальную и |
|||||||
соленоидальную |
части |
|
|
|
|
|
|
||
|
|
|
« = grad(t> + r o t I , |
|
(278) |
||||
где Ф — скалярный |
потенциал; |
|
|
|
|
|
|||
А — векторный |
потенциал. |
|
|
|
|
|
После подстановки выражения (277) в (278) получим систему уравнений, решение которой можно представить в виде следующих выражений:
У2 Ф — Дг Ф" — ( З Х ; - ' ; У °т ( Т _ 7 ) = о; |
(279) |
|
с - |
Л. -г ^[L |
|
У2Л
V2 Т 1- Г — Г ° ( 3 ^
2 |
A, -J- 2fl . |
c i = |
— , |
— 0; |
|
(280) |
|
; 2 |
| 1 ) У2 Ф' = 0; |
(281) |
|
2 |
= |
ц |
|
с2 |
— , |
|
где |
СІ |
— скорость |
|
распространения |
упругой безвихре |
||||
|
|
вой волны |
(волны расширения); |
|
|
||||
|
с 2 |
— скорость |
распространения |
упругой волны |
иска |
||||
|
|
жения и поворотов, |
вызывающей |
изменение фор |
|||||
|
|
мы без |
изменения |
объема; |
|
|
|||
Т0 = const—температура |
тела в |
напряженном |
состоянии, при |
||||||
|
|
котором |
Ф = |
0. |
|
|
|
|
|
Исключив Т из выражения (279) и (281), получим уравнение для |
|||||||||
функции |
Ф |
|
|
|
|
|
|
|
|
|
~ 4 - І) Г - і ' |
Щ* ~ |
|
|
™ = о . |
№ |
где с£ — удельная объемная теплоемкость при постоянном тензоре деформации.
Если пренебречь взаимодействием поля деформаций и темпера турного поля, можно получить представление общего решения ди намической задачи термоупругости, в котором скалярный потенциал
Ф и векторный потенциал А определяются из уравнений
• 5 Л = 0,
где
,2 |
. „ я |
і |
а 2 |
с
В поставленной задаче поле деформаций определяется в зависи мости от температурного поляТакая задача называется связанной. Решение связанной задачи термоупругости в общем случае пред ставляет значительные математические трудности. Для приближен ного решения этой задачи используют вариационный принцип.
В зависимости от условий теплообмена связанную задачу можно упростить.
1. При неравномерном нагреве, обусловленном внешним источни ком тепла, можно считать, что температурное поле не зависит от вызываемых им деформаций. Тогда в уравнении теплопроводности (274) не учитывается член механической связи
Г„ (& + 2ц) а т ,,-
лt-kli
изадача термоупругости становится несвязанной. В этом случае поле деформаций определяется независимо от температурного поля. Это
позволяет решить дифференциальное уравнение теплопроводности,
не определяя поле |
деформаций. |
2. Инерционный |
член рщ в уравнении (273) означает, что опре |
деление температурных напряжений в случае изменяющихся с те чением времени температурных полей является по существу не ста тической, а динамической задачей. К динамическим можно отнести также задачу определения температурных напряжений при тепловых ударах, которые имеют место в газовых турбинах при пуске. Как показали исследования, проведенные в работе [48], динамические напряжения имеют существенную величину только при скачкообраз ном изменении температуры поверхности или окружающей среды. Но такое изменение температуры является математической идеали зацией, которая физически не может осуществиться. Все реальные
тепловые |
удары имеют |
конечную |
скорость изменения |
температуры, |
|||
т. е. температуры среды |
Тг достигают конечного значения за отлич |
||||||
ный от нуля интервал времени AT. С увеличением безразмерного |
|||||||
времени |
максимум |
напряжений |
быстро |
уменьшается. |
При т = 3 |
||
(что соответствует |
1СГ12 |
с для углеродистой |
стали) этот максимум |
||||
составляет лишь 14% его значения при т 0 |
= |
0. Таким образом, даже |
при столь малой продолжительности нагрева максимум динамических напряжений снижается до 86%. Эти результаты показывают, что при встречающихся в действительности соотношениях повышение напряжений вследствие динамических эффектов не имеет практи ческого значения. Поэтому температурные напряжения в упругом теле в определенный момент времени можно определить исходя из состояния температурного поля в рассматриваемый момент без учета сил инерции, соответствующих движению частиц тела при перемен ном тепловом расширении, что равносильно отбрасыванию в уравне нии (273) члена рщ.
Постановка задачи термоупругости с учетом рассмотренных до пущений называется квазистатической.
При решении отдельных задач термоупругости в квазистати ческой постановке обычно принимают в качестве основных неизвест
ных компоненты |
вектора перемещения и.ц или |
компоненты |
вектора |
напряжения ап-. |
В соответствии с этим различают постановку |
задачи |
|
в перемещениях и в напряжениях. |
|
|
|
Постановка задачи в перемещениях. Для |
постановки квазиста |
тической задачи термоупругости в перемещениях используется пер вое уравнение системы (277). Отбрасывая инерционный член ри'ї
и считая, что объемные силы в теле отсутствуют, получаем основное уравнение
liYhi - f (X -\- |л) grad div и — (ЗА, -1- 2ц) а т grad Т = 0. |
(284) |
Функция Т предполагается известной из решения задачи теплопро водности.
Общее решение уравнения (284) можно представить в виде суммы общего решения однородного уравнения и„ и частного решения неоднородного уравнения и^.
Общее |
решение |
однородного |
уравнения иа |
было найдено |
|||
П. Ф. Папковичем |
|
|
|
|
|
|
|
|
«0 = |
4(1 —v)B |
— pad (В? — В0), |
(285) |
|||
где В — гармонический |
вектор, |
удовлетворяющий |
уравнению |
||||
|
|
|
V2B = 0; |
|
|||
В0—гармонический |
|
скаляр, |
удовлетворяющий |
уравнению |
|||
|
|
|
V-Яо = |
0; |
|
||
/•— |
радиус-вектор. |
|
|
|
|
|
|
Частное решение |
ит |
имеет |
вид |
[32] |
|
||
|
|
|
йт |
= |
grad Ф, |
(286) |
|
где Ф удовлетворяет уравнению |
Пуассона |
|
|||||
|
|
™=Т=І«т{Т-Т0). |
|
(287) |
|||
Функция Ф называется термоупрушм потенциалом |
перемещения. |
||||||
Порядок |
решения сводится к |
следующему. |
|
По известному температурному полю находят решение уравнения (287). Вычисляют температурные напряжения, отвечающие потен циалу Ф, которые не удовлетворяют заданным условиям на поверх
ности. На |
полученное |
решение накладывается решение уравне |
ний теории |
упругости, |
соответствующее Т — 0, чтобы выполнить |
заданные условия на поверхности. Иногда квазистатическая задача термоупругости сводится к задаче изотермической теории термо упругости, если величину— (ЗА,— 2ц.) аг grad Т рассматривать как вектор плотности объемной силы, а за поверхностные силы принять равномерное нормальное к поверхности растяжение— (ЗА. -|- 2ц.)ссгх
Х(Т-Т0).
Постановка задачи в напряжениях. При постановке задачи термо упругости в напряжениях [30] решение сводится к нахождению шести функций Оц, удовлетворяющих трем уравнениям равновесия
а,;, І = 0; |
(288) |
— шести уравнениям совместности деформаций |
в напряжениях |
(1 -h v) а,-/, па + а'ш, а + £ ( , ' ~ V ) <*тТ: « А / + ЕатТИ = 0; (289)
— трем граничным условиям |
|
аип, - 0, |
(290) |
где п,- — направляющие косинусы, определяющие положение внеш
|
ней нормали к поверхности тела. |
(N — связная), то |
|||
Если |
рассматриваемая |
область многосвязна |
|||
ее посредством (N— 1) разрезов можно превратить в |
односвязную, |
||||
для которой справедливы |
соотношения |
(288) — |
(290). |
|
|
Зная |
напряжения, по |
зависимостям |
(288) — (290) |
с помощью |
|
соотношения |
|
|
|
|
|
|
e„ = |
- п Й - 8,/ W |
T - |
Т0) 6tJ |
(291) |
можно определить деформации, а затем и перемещения.
§ 62. Методы расчета температурных напряжений
Классический метод расчета температурных напряжений по существу состоит в интегрировании дифференциаль ных уравнений с дополнительными температурными членами, что эквивалентно решению вопроса о напряженном состоянии при дей ствии фиктивной поверхностной и объемной нагрузки. Преимущество этого метода состоит в том, что когда удается получить аналити ческое решение, появляется возможность широкого исследования напряженного состояния тела. Основным недостатком этого метода следует считать то, что при сложной системе нагрузок, имеющих место в узлах современных газовых турбин, решение задач с помощью классических методов значительно затруднено. Поэтому в настоящее время существуют решения лишь сравнительно несложных задач
отемпературных напряжениях.
Метод, основанный на обобщении теоремы о взаимности работ.
Этот метод заключается в определении напряженного состояния
в упругом теле под действием неравномерного температурного поля |
||
и сводится к задаче изотермической теории упругости о напряжен |
||
ном состоянии упругого тела под действием |
сосредоточенной |
силы. |
Рассмотрим два напряженных состояния |
упругого тела, |
из ко |
торых первое характеризуется |
напряжениями |
а(1-, |
деформациями |
£ц и перемещениями щ, возникающими под действием внешних сил F, |
|||
/ и температурного поля Т, а второе — напряжениями |
а],-, деформа |
||
циями г}/ и перемещениями и\, |
возникающими |
под действием внеш |
|
них сил F'lf'i и температурного |
поля 7". |
|
|
Работа сил первого состояния на перемещениях второго состоя ния с помощью формулы Остроградского—Гаусса при использова нии уравнений равновесия и граничных условий в напряжениях может быть представлена так:
L b 2 |
= J" Ftu'i dV + |
j fiiii dQ = |
j otj&tj dV. |
(292) |
|
v |
a |
v |
|
Аналогично можно получить выражение для работы сил второго состояния на перемещениях первого
Цл= \ a'ueudV. |
(293) |
v |
|
Из уравнений (292) и (293) с помощью соотношений (275) и (291)
находим |
|
|
L 1 2 - L 3 l = а т Jv |
\(Ґ - T'o)Gkk - (Т - Т0) о!.,] dV. |
(294) |
Зависимость (294) обобщает теорему о взаимности работ для слу чая статической и квазистатической задач термоупругости. Ее можно использовать для нахождения перемещений, возникающих в опре
деленной точке тела при неравномерном нагреве. |
|
|
|
||||
Пусть |
F[ = О, ft = 0, 7" = Т 0 , а система |
внешних сил |
F't и |
f\ |
|||
сводится |
к сосредоточенной единичной силе, |
приложенной |
в точке |
||||
х\ и направленной |
параллельно осп хп |
которая вызывает |
в точке |
х |
|||
напряжения atj (х\; |
х{) и деформации |
е ( / (xjj; х%). В этом |
случае |
из |
уравнения (294) можно получить зависимость для определения пе
ремещений в точке |
от действия температурного поля |
Т: |
|||||||||
|
щ (4) |
= |
(ЗЯ, + |
2ц) ат\(Тv - |
Т0) e?kl{ ( 4 |
хъ) dV |
= |
||||
|
|
|
= ат |
\ (Т - |
Т0) okk |
(4, |
* 6 ) dV. |
|
(295) |
||
|
|
|
|
|
v |
|
|
|
|
|
|
Если модуль упругости Е и коэффициент Пуассона |
v (или А. и |
||||||||||
ц.) зависят |
от температуры, т. е. являются заданными функциями |
||||||||||
координат |
Х|, то |
зависимость |
(295) |
принимает |
вид |
|
|||||
|
щ (4) |
= |
J "(ЗА + 2ц) сст |
(Т - |
70 ) йи (4, |
Ч) dV |
= |
||||
|
|
|
v |
|
|
|
|
|
|
|
|
|
|
|
= |
v\ а т |
(Т - |
Т0) |
akk |
(4, |
л-6) dV. |
|
(296) |
Определив щ (xl) по формуле (296), можно найти компоненты тензоров напряжения и деформации по зависимостям (275) и (276).
Рассмотренный метод позволяет исследовать температурное на пряженное состояние тела при различных начальных и граничных условиях и учесть зависимость £ и v от температуры путем сведе ния задачи термоупругости к задаче изотермической теории упру гости. С помощью этого метода можно легко учитывать наличие внеш ней нагрузки. Метод удобен для решения, так как является особым методом интегрирования совокупности системы дифференциальных уравнений второго порядка в частных производных. Недостатком метода является то, что он требует наличия готовых решений изо термических задач для тел, подверженных действию сосредоточен ных сил.
Следует отметить, что рассмотренный метод позволяет во многих случаях получить экспериментальное решение вопроса о напряжен ном состоянии тела [44].
Вариационные методы. Вариационные методы во многих случаях являются эффективными для приближенных расчетов температурных
напряжений. |
W0 на единицу объема тела в |
|
|
|
Энергия деформации |
напряженном |
|||
состоянии может быть записана в следующем виде |
[10]: |
|
||
2№0 = aifit; = °х.А-л- + |
Ъуу%у -І- *zPa + охуъхц + Охггхг |
+ |
оугги2. |
(297) |
Полная энергия деформации тела |
|
|
|
|
|
W=\wvodV. |
|
|
(298) |
Зная энергию деформации и силы, действующие на тело, можно вывести ряд теорем и принципов: принцип наименьшей работы; тео рему о минимуме потенциальной энергии; теорему Кастильяно и др.
Рассмотрим два принципа, предполагая, что тело находится под действием поверхностных сил и объемных сил при известном темпе ратурном поле.
Обобщение вариационного |
уравнения |
Лагранжа |
для случая за |
дачи термоупругости. Пусть телу сообщены виртуальные переме щения би,-, удовлетворяющие всем кинематическим граничным усло виямНа основании уравнений равновесия можно записать
|
|
|
|
\{olhj |
+ Ft)buidV |
= Q. |
|
(299) |
|
|
|
|
|
v |
|
|
|
|
|
|
По формуле Остроградского—Гаусса |
(294) выражение приводится |
|||||||
к виду |
|
|
|
|
|
|
|
||
|
|
\ OjjiijdUi dQ + |
f Ffiut dV—\ |
ац8еиdV |
= 0, |
(300) |
|||
|
|
Q |
|
|
V |
|
V |
|
|
где |
fi |
— замкнутая |
поверхность. |
|
|
|
|||
|
Введем специальные термодинамические функции: |
|
|||||||
|
— |
плотность |
свободной |
энергии |
|
|
|
||
|
|
|
|
F — V — TS; |
|
|
(301) |
||
|
— |
плотность |
термодинамического |
потенциала |
Гиббса |
|
|||
|
|
|
|
|
G = |
F — оИЕИ. |
|
|
(302) |
Из |
выражений |
(301) |
и (302) |
имеем |
|
|
|
- & = - « « • |
<зм> |
17* |
259 |