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

книги из ГПНТБ / Жаров Г.Г. Судовые высокотемпературные газотурбинные установки

.pdf
Скачиваний:
21
Добавлен:
25.10.2023
Размер:
17.85 Mб
Скачать

Глава 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 —

 

 

_

Е

 

Л -

(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 — связная), то

Если

рассматриваемая

область многосвязна

ее посредством (N1) разрезов можно превратить в

односвязную,

для которой справедливы

соотношения

(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—\

ациdV

= 0,

(300)

 

 

Q

 

 

V

 

V

 

 

где

fi

— замкнутая

поверхность.

 

 

 

 

Введем специальные термодинамические функции:

 

 

плотность

свободной

энергии

 

 

 

 

 

 

 

F — V — TS;

 

 

(301)

 

плотность

термодинамического

потенциала

Гиббса

 

 

 

 

 

 

G =

F оИЕИ.

 

 

(302)

Из

выражений

(301)

и (302)

имеем

 

 

 

- & = - « « •

<зм>

17*

259