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

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

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

Знание хотя бы одного термодинамического потенциала позво­ ляет определить все термодинамические параметры. С учетом (303) и условий на поверхности

 

 

 

о</"•/ =

//

 

 

 

зависимость

(300)

принимает

вид

 

 

 

 

б '\FdV

— \ FiUi

dV — \

dQ]

0.

(305)

 

V

 

V

 

Q

 

 

 

При отсутствии

объемных

н поверхностных

сил

 

 

 

 

б =

\FdV

=

Q.

 

(306)

 

 

 

 

v

 

 

 

 

Уравнение (305) является обобщением начала возможных пере­

мещений Лагранжа для случая упругого равновесия.

 

Обобщение

принципа

минимума

потенциальной

энергии

деформа­

ции для случая задачи

термоупругости-

Если тело подвержено та­

кой произвольной вариации, при которой новые компоненты тензора

напряжения

ац

+ 5а,-у удовлетворяют

уравнениям равновесия, то

 

 

 

 

 

 

6ff,/ i /

= 0.

 

 

 

(307)

Для удовлетворения граничных условий необходимо вариации

поверхностных

 

сил

подчинить

условию

 

 

 

 

 

 

 

 

 

8fi = baunh

 

 

 

(308)

При этом должны быть удовлетворены все условия кинемати­

ческих связей. На

основании

зависимости

(304)

можно записать

 

 

 

 

J ( e / /

+

- ^ - ) 6 a i / d V = 0.

 

(309)

 

 

 

 

v

 

 

 

 

 

 

 

 

Выражая компоненты тензора деформации е£/- через компоненты

вектора перемещения u-L

и

применяя

формулу

Остроградского —

Гаусса, из

(309)

имеем

 

 

 

 

 

 

 

 

 

j

u^OijiijdQ

f щдоц, ,-dV +

J 8G dV = 0.

(310)

 

Q

 

 

 

V

 

 

 

V

 

 

Принимая

во

внимание

зависимости (307) и (308), получаем

 

 

 

 

8(\GdV+

 

с1^а\=0.

 

 

(311)

Если поверхностные силы отсутствуют

или выполняется

условие

их неизменяемости

(б/,- = 0 ) ,

то

 

 

 

 

 

 

 

 

 

 

&\GdV

=

0.

 

 

(312)

Эта формула обобщает принцип минимума потенциальной энер­ гии деформации. Сформулировать этот принцип можно следующим образом [10]: при перемещениях, удовлетворяющих одновременно

граничным условиям и условиям равновесия, потенциальная энергия достигает минимума.

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

имела минимальное

значение.

 

 

 

 

 

 

 

Пусть

потенциальная

энергия

имеет

вид

 

 

Р = \Рг

(х, у,

z,

S,

Sx,

Sy,

Sz,

Sxx,

 

Sm,

 

Sxy, Sxi,

Syz)dV,

(313)

v

 

 

 

 

 

 

 

 

 

 

 

 

 

где

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

S

=

S{x,

y,

z);

 

Sx

=

- ^ ;

 

 

 

 

 

 

c,

_

_dS_ .

-

_

dS

 

 

 

 

 

 

v

~

ду

'

2

dz

'

 

 

Если заменить S

через 5 (x, у, z)

+

є„

(x, у,

z), где n

(x, y, z)

функ­

ция, обращающаяся в нуль

на

границе;

є—произвольный

малый

параметр. Тогда потенциальная энергия Р имеет минимум при є = О

и

dp

 

 

 

 

 

 

 

 

 

 

 

 

 

 

- j - =

0.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

as

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Используя

это условие,

 

можно прейти к выражению

 

 

 

 

 

 

 

 

 

\L(S)n(x,

у, z)dV

= 0,

 

 

 

(314)

где

 

 

 

 

 

V

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

{

'

 

dS

 

дх

'

dSx

ду '

dSy

dz

'

dSz

i "

 

 

 

 

 

 

JP_

 

дРг

 

 

дРг

і JP_

dPx

 

 

 

 

 

 

 

~^

dx2

'

dSxx +

 

ду2 '

dSyy

dz2 '

dSzz

+

 

 

 

 

,

9 (

J

P _

 

дР,

 

d2

_dPj_

 

 

_дР\_\

n

i 4 .

 

 

"г"

 

\

дхду

'

dSxy

r

dxdz

' dSxz

"г" dydz

'

dSyz

J '

v

;

Уравнение (315) — уравнение Эйлера вариационной задачи, определяемой выражением (313). Решив уравнение (315), можно найти функцию 5 (х, у, z). Приближенное решение этого уравнения можно получить, использовав метод Рейли—Рица, который заклю­ чается в следующем. Задаются перемещения

Чп

t chnRM,

(316)

где Rin удовлетворяют тем же граничным условиям, что и"«,-.

Подставив иіп в (313) и выразив Р как функцию aUl, получим выражение, удовлетворяющее принципу минимальной потенциаль­ ной энергии:

откуда можно определить и1п.

Более подробно о других

методах

изложено в соответствующих разделах математики.

 

Преимущество

вариационных

методов заключается прежде всего

в том, что они

позволяют найти приближенное решение

многих

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

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

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

турбины часто можно ограничиться

решением

двумерной задачи.

Прежде чем перейти к анализу

методов

решения двумерных

задач напряженного состояния тел,

рассмотрим

в квазистатической

постановке две типовые плоские задачи термоупругости: о плоской деформации тела и о плоском напряженном состоянии тела.

§ 63. Плоская задача термоупругости

Задача о плоской деформации. Плоская дефор­ мация возникает в длинном цилиндрическом или призматическом теле (лопатке, роторе) при плоском температурном поле Т(х,у,%).

Для плоской деформации характерными являются перемещения

и = и(х, у); v — v(x, у); w = 0,

(318)

при которых

О",хг

Приняв в уравнении (291) ги = ггг = О, получим

v(o^ + o ^ ) - a x £ C T - 7 0 ) .

(319)

Тогда решение задачи о плоской деформации можно свести К решению следующей системы уравнений (объемные и поверхностные силы отсутствуют):

 

 

дохх

+і

даХу = 0 ;

 

 

 

 

дх

'

ду

 

 

 

 

 

дахУ

і

даии

 

 

 

 

дх

1

ду

 

 

 

& х

х =

1ЁТ ( ° Х Х ~

^ У

У >

~'~а т ' ^

~~ Т ° ) '

(320)

НУ

— -щ- ІРуу

— v\csxx)

+ а Т і (Т — То);

 

 

 

р

1 Н- vx

,

 

 

 

 

1

l

и ГТ

 

 

 

 

 

 

 

 

 

 

^Ху

£

 

I ^ХУ

 

 

 

дии,

диw

 

__ 1 / ди .

dv \

~дх~''

гУУ-~ду~'

£ху-~2

\~dy~-r

-fa)

є,-,-- —

; е „ „ _ — ;

при граничных условиях на наружном контуре, которые можно задать либо в напряжениях:

<УххПх + охи = !х(х, у);

«V і * + °УУПУ

=

U (х> У)>

(321)

либо в перемещениях:

 

 

 

и = 8i

(*> У);

 

" -

ёг

(*. у),

 

где

 

 

 

Задача о плоском напряженном состоянии. Плоское напряженное состояние возникает в тонкой пластине (оболочковые лопатки), поверхности которой свободны от внешних усилий. При этом

°ы = ахг = в иг — °-

Решение задачи в этом случае сводится к решению следующей системы:

двхх

 

і

дахУ

_ и-

 

 

~дГ

+

 

~ d y - - U '

 

 

доху

_|_ дауу

_ Q .

 

 

дх

 

'

ду

 

 

 

*хх = 4 " ( f f * *

™ У Н ) +

а т ( г

~ ^о) ;

 

єда = 4 " (°w

w

«

) +

а т ( J '

Т'о);

(323)

263

 

 

 

 

*zz

= —

 

(0-.v.v + Vyy)

+

« г

(T

To)\

 

 

 

 

 

 

 

ди

 

_

 

dv

 

 

_

 

1

/

ди .

dv

\

 

 

 

Є х х ~

IF'

 

гш

~

~ду~ ' &

x b

~

~~2~ \ду

 

/

 

при граничных

условиях

на

наружном

контуре

в напряжениях

и в перемещениях (322).

 

 

 

 

 

 

 

 

 

 

 

 

 

Решение задачи в перемещениях. Решение задачи в этом случае

сводится к

решению

системы

(320),

которую можно

представить

в

следующем

 

виде:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

рФи +

(Л +

^ ) . | _ ( * L +

- * L ) _

{ З Х

 

+

2ц.) « т - g -

=

0; (324)

 

+

(Я. +

, ) ^

(-g- +

 

-

(3* +

2ц) а т

- f -

=

0

с

граничными

условиями

(322). Если

в

(324)

подставить

величины

 

 

 

 

 

к =

~

 

> Ц

 

2(1

-l-v,)

 

 

 

 

 

 

 

 

 

1

 

V j

 

 

а, Е,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(ЗХ +

2р.) а т

 

 

-1

 

 

 

 

 

 

 

 

 

 

 

 

 

- V,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

и заменить

£ х ,

 

v b

а т 1

на

Е,

v,

ост,

то

получим систему уравнений

в перемещениях для плоского напряженного состояния, к которой сводится система (323).

Пользуясь классическим методом, решение системы (324) можно искать как сумму общего решения однородного уравнения и частного решения неоднородного уравнения.

Частное решение соответственно зависимости (286) имеет вид

" т — дх ' ° г ~ ду '

где термоупругий потенциал перемещений определяется из урав­ нений:

— для плоской деформации

 

^

= \ ~ ^ ( Т - Т 0 ) ;

(325)

— для плоского напряженного состояния

 

V2 O = ( l + v ) a T ( 7 , - 7 , 0 ) .

(326)

К частному решению необходимо присоединить общее решение однородной системы уравнений, содержащее необходимое число постоянных интегрирования для удовлетворения граничных усло-

вий (322). Такая постановка задачи пригодна как для одиосвязных, так и для многосвязных тел.

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

 

Постановка

задачи

для

односвязных

тел. В этом случае

систему

уравнений (320) можно свести к неоднородному

бигармоническому

уравнению

 

 

 

 

 

 

 

 

 

^ + • ^ « ^ = 0,

 

( 3 2 7 )

а

систему (323) к

 

 

 

 

 

 

 

 

Vdcp + EaTV2T

= 0

 

(328)

с

граничными

условиями

(при отсутствии поверхностных

сил)

 

 

 

 

Ф = |f=

0-

 

(329)

 

Постановка

задачи

для многосвязных

областей.

Для многосвязных

тел к системе (327)—(329) добавляются граничные условия на вну­

тренних контурах:

akx

+ $ky

+

у,

(330)

Ф =

дф

 

дх

і

ft

dy

да

~

к 'дп

 

P l i

dn

-

На контуре

L/; при

/г =

1,

2, . . ., N,

где N — число внутренних отверстий.

3N

постоянных (aft; рА ; yk)

определяются из удовлетворения 3N условий однозначности для

перемещений и углов

поворота

[30J

 

 

 

 

 

 

 

, д n , . . J O

 

, _. f дТ

 

 

(331)

 

 

Е

J ^ V v d S

+ ccJ

*LdS = 0;

 

_1_

 

 

 

 

 

 

 

 

 

 

Е

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(332)

4 - J

{y-W + xl)^dS

 

+ ar$

^

+ xl)TdS

= 0. (333)

Lk

 

 

 

 

Lk

 

 

 

 

Условия (331) соответствуют плоскому напряженному состоянию.

Если в этих условиях заменить Е и а т

на Ег

и сст 1 , получим

условия

однозначности,

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

плоской деформации.

 

Физическое

толкование

неоднозначности

перемещений

и углов

поворота

в многосвязных

телах.

Неоднозначность

перемещений и

углов поворота

в многосвязных

телах

обусловлена

напряжениями,

возникающими в них не от действия внешних сил, а вследствие обра­ зования особого рода деформаций, называемых дислокациями [30].

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

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

замкнутого контура

L z

в точке

р0

возникают (при отсутствии тем­

пературного

поля)

следующие Аи,

До, Дш,:

 

 

 

 

Аи = и + -

и

= -

уАи2

-

- L }

(х-^-

-

У

~ )

V2cpdS;

(334)

Av = v+-v~

 

 

= л-Асо2 --LJ

{y-£f+

 

хж)

V

V S ;

(335)

 

 

 

 

 

 

 

Ч

 

 

 

 

 

 

 

Дсо2

=

(о+ -

соГ =

-g- \ 4i

^(pdS.

 

(336)

Учитывая

(331)

 

и

(336),

получаем

равенства

 

 

 

 

а т і

 

(х ~д§

У'дії) T

d S

= A " — </AoV-

 

 

 

аЛ

Ж ~ х

ш ) T d

S

= Л у

_

х А с о

-

 

<337)

 

а т {

 

|

^

=

- Дсо г

 

 

k=\,2,...,N,

 

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

В случае плоской деформации в уравнениях (334)—(336) и ра­ венствах (337) величины Е и ат следует заменить на Ех и а т 1 . За­ метим также, что с помощью дислокаций решается важный для оптического метода (экспериментального метода определения тем­ пературных напряжений) вопрос — о переходе от напряжений в мо­ делях к напряжениям в действительной конструкции [39].

§64. Методы решения плоской задачи термоупругости

Методы, использующие функции комплексного переменного. В настоящее время методы, использующие функции комплексного переменного, являются наиболее общими в применении

•к плоской задаче теории термоупругости. Если представить функцию напряжения как разность двух функций

Ф = U — V,

(338)

то решение двумерной задачи можно свести к решению следующей

системы уравнений

 

[10,

39]:

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

У*и =

0,

WV =

k(T

— T0),

 

 

(339)

где

/г =

Еат

— для

плоского

напряженного

состояния;

 

k —

_ т ^

— для

плоской

деформации.

 

 

 

 

 

Функция U может быть выражена

через аналитические

функции

в виде

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

Тогда

 

 

 

2U =

2 Ф і

(г) + г Ф і (z) +

^ 2

(2)

+

аР 2

(2).

 

(340)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

cxx

+ oyy

=

Va(U-V)=

 

4 a

' f

r v

)

,

(341)

 

 

 

 

 

 

 

 

 

 

 

 

 

 

dzdz

 

 

 

где

 

 

 

 

 

z

= x

+

іу;

z =

х

і у;

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

д2

,

_д*_

_

.

д2

 

 

 

 

 

 

 

 

 

 

 

 

 

дх2

'

ду»

 

d z d z

 

 

 

 

 

или,

имея

в

виду

(340),

 

 

 

 

 

 

 

 

 

 

 

а

 

 

охх

+ ови

=

2 [фі (2) +

фі (і)] — £ (Т — Т0),

(342)

 

 

 

,

 

 

. „.

 

d2U

d2U

 

0 .

d2U

 

 

 

 

 

 

 

 

 

 

 

 

 

 

І

"хх

 

І

А Л ;

2

Д У

І

 

 

Д Х Д

У

 

Принимая

во

внимание,

что

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

д2

2/

д2

 

,

д2

 

 

 

 

 

 

 

 

 

дх2

 

ду2

 

\dz2

1

QZ*

 

 

 

 

 

 

 

 

 

2i

9 2

 

 

-

/ д2

 

 

д2

 

 

 

 

 

 

 

 

 

 

 

дхду

 

 

у dz2

 

 

dz2

 

 

 

 

уравнение (343)

можно

преобразовать:

 

 

 

 

 

 

 

 

 

 

Оуу ахх

+ 2 ю ^ = 2 [гфГ (г) +

г|н (г)] —

 

где

Уравнения (342) и (344) определяют напряженное состояние тела,

при этом на функцию Ф і (z)

можно наложить

условия

[10]

Ф І (0) =

Re [-L Ф І (0)] =

0.

(345)

267

Обычно функции фі (z) и і)?! (г) представляют в виде полиномов

 

СО

д„2<<;

СО

/ г -

cPl

(г) = £

^ (г) =; £ 6„z»,

(346)

где cn и Ь п — комплексные

числа.

образом,

Отображающая

функция

г =

со (S) выбирается таким

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

Метод конечных разностей

(метод сеток). Выражения (327) и

(328) можно преобразовать в

конечно-разностные уравнения вида

 

 

Р " { [б (7 + ~

) +

8] ер,?1 -

4 [(1 + у) X

 

 

 

 

Л?*

 

 

 

 

 

 

 

 

 

X (cpffi. , + Ф Ш , /) +

( 1 +

- f ) ( ф ^ - і

+ Ф/'+ж)]

+

 

 

 

+ 2 (Ф £1. у _ ! + ф ^ } ,

-f- Ф # і. / - і + Ф Ш , Л И ) +

 

 

 

+ Y ( Ф Й . ; + Ф Ф . /) + Y

< Ф ? Ж + Ф / ' + Ц +

 

 

' к,^

[7 (Я[±1, ,•+#{#./ + Hit}-i

+

tfft+r-

2(1-|- 7) H[f\

= 0,

(347)

где

 

 

 

 

 

 

 

 

 

 

Граничные

условия

 

 

 

 

 

 

 

 

 

 

/4-і

Ф ^ Ч 2 А ( ^ . ) 4 =

Ф а - Ч Л ,

 

(348)

 

 

фз/г

 

 

 

 

 

 

 

 

 

 

 

где

ф Ж

— значение функции

напряжений

в законтурной

точке

 

 

по определенному

направлению

или у);

 

 

Ф ^ 1 — то же в предконтурной точке;

h — расстояние от предконтурной точки до контура в рас­ сматриваемом направлении;

значение на контуре производной по тому же на-

правлению.

Уравнения (347) с граничными условиями (348) решаются ме­ тодом итераций. Обычно две-три итерации обеспечивают необходимую точность.

Вычислив функцию напряжений Ф-+1, можно определить на­ пряжения по следующим соотношениям:

eV^'-'ff+'W", (349)

« ,

=

' - l ' <

"2

 

;

(350)

+

W

' і і,

/—1

' і і,

і 't + i .

У—1

1 ' ' + ' . /+1

(351)

а затем по (319) найти огг для задачи о плоской деформации. Преимущество метода конечных разностей состоит в том, что

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

§65. Температурные напряжения

вохлаждаемых лопатках турбин

Рабочие лопатки без каналов и с каналами охлаждения. Расчеты температурных напряжений в рабочих лопат­ ках можно выполнять двумя способами:

— определяют только продольные составляющие температурных

напряжений, которые складываются

с

напряжениями растяжения

от центробежных сил и напряжениями

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

и поэтому представляют наибольшую

опасность;

— определяют общее напряженное состояние лопатки, вызванное неравномерным температурным полем.

I . Если предположить, что изменения температуры по высоте лопатки не происходит, то для балки с несимметричным поперечным сечением справедливы формулы для определения напряжений при косом срезе [22]:

o'z!=:aTE(Tcp

— kyX — k 2 y ) — a?ET;

(352)

Tcp

= ±\TdS;

(353)

 

б"