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

Механика композитных материалов 2 1980

..pdf
Скачиваний:
6
Добавлен:
15.11.2022
Размер:
11.16 Mб
Скачать

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

Для окончательного решения вопроса о целесообразности и возмож­ ности практического применения балласта необходимо оценить абсолют­ ную величину оптимального давления р0пт и сравнить ее с прочностью однонаправленных композитов на сжатие поперек волокон (Пг~). Рас­ четы показывают, что для исследованных композитов ропт имеет порядок 1000 кгс/см2 при яр = Зн-10, в то время как Пг~ для этих же композитов лежит в диапазоне 660—2500 кгс/см2. Здесь следует отметить, что отно­ сительно небольшие изменения размеров несущего диска в районе т**\ (см. рис. 2) связаны с весьма существенными изменениями р0пт (это видно из рис. 1). Например, при использовании легкого балласта опти­ мальный размер несущего обода несколько увеличивается и при яр = 0,25 Ропт равен примерно 500 кгс/см2 для всех материалов. Если р0пт>Пг~ для балласта или композита, то следует выбирать значение давления, равное допустимому значению радиальных напряжений, и по графику рис. 1 определять соответствующее ему оптимальное значение т\. Энер­ гоемкость обода при выбранном т\ можно определить по графику рис. 2. Энергоемкость балласта вычисляется по формуле (4), в которую под­ ставляется значение, определяемое выбранным давлением (формула (1).

4. Рассмотрим случай, когда в качестве балласта используется диск из изотропного материала с модулем упругости EQ{2\ много меньшим, чем окружной модуль EQР) материала обода, и соединенный не только с ободом, но и с валом. Таким образом, диск выполняет роль не только балласта, но и центра или ступицы маховика; это приводит к существен­ ному упрощению изготовления маховика и его балансировки. В ряде случаев в качестве центра может служить оправка, на которую наматы­ вался обод маховика. Оценим возможности повышения энергоемкости маховика при использовании балласта такого типа.

Поскольку E BW ^ E QW, то при определении напряженного состояния в диске можно с достаточной точностью принять, что перемещения на поверхностях контакта диска с валом и ободом равны нулю. В этом слу­ чае безразмерные напряжения в балластном диске можно получить, ис­

пользуя зависимости из [8] в следующем виде:

 

 

 

 

аг(2)

 

пр

[ (1 +v) (tn^ + m2)

+

а г<2>= pV(1)(D2fe2(3 + Vf.0(l))

8(3+v,-e(1))

 

m\2m2

 

 

 

(5)

 

+ (1 — v)

P2

 

 

 

 

 

 

 

 

 

ае(2)=-

сте( 2)

 

 

2^

_

руР)(й262(3 + Vr0(1))

8 (3 + Vr0(1))

[ (1 + v) (m{2+ tn )

 

 

 

 

 

 

— (1—v)mi2m2/p2— (l + 3v)p2],

 

 

где v — коэффициент Пуассона материала балласта.

 

и обода

Из формул (5) следует, что на границе балластного диска

 

возникают сжимающие напряжения

 

 

 

 

 

Пр(т\2— т2)

(2)

" ( 2 )

 

(6)

 

 

 

 

= VO\,

 

4(3 + vre(1))

|р= 77Ц

p = m ,

mt

10,00 0,420

7,00 0,450

5,00 0,477

3,00 0,518

1,00 0,596

0,50 0,640

0,25 0,668

10,00 0,380

7,00 0,413

5,00 0,444

3,00 0,488

1,00 0,570

0,50 0,610

0,25 0,641

Отношение удельных энергоемкостей*

Стеклопластик

1,02 (0,77+0,25)

0,20 (0,15+0,05)

1,07 (0,82+0,25)

0,24 (0,19 + 0,05)

1,10 (0,86+0,24)

0,29 (0,23 + 0,00)

1,14 (0,92 + 0,22)

0,37 (0,30+0,07)

1,14 (1,00 + 0,14)

0,56 (0,49 + 0,07)

1,12 (1,01+0,11)

0,68 (0,62 + 0,06)

1,08 (1,01+0,07)

0,78 (0,74+0,04)

Боропластик

0,83 (0,68+0,15)

0,19 (0,16 + 0,03)

0,89 (0,74+0,15)

0,24 (0,20+0,04)

0,94 (0,78+0,16)

0,27 (0,23 + 0,04)

1,00 (0,85+0,15)

0,32 (0,27+0,05)

1,07 (096+0,11)

0,56 (0,50+0,06)

1,06 (0,98 + 0,08)

0,68 (0,63+0,05)

1,04 (0,99 + 0,05)

0,79 (0,75 + 0,04)

* См. примечание к табл. 2.

0,388

0,426

0,461

0,511

0,601

0,644

0,683

0,478

0,512

0,543

0,588

0,679

0,724

0,765

Отношение удельных энергоемкостей*

Углепластик

0,80 (0,65 + 0,15)

0,16 (0,13+0,03)

0,87 (0,71+0,16)

0,20 (0,16 + 0,04)

0,94 (0,76+0,18)

0,23 (0,19 + 0,04)

1,03 (0,85+0,18)

0,31 (0,25+0,06)

1,14 (0,99 + 0,15)

0,51 (0,45 + 0,06)

1,14 (1,03 + 0,11)

0,65 (0,59 + 0,06)

1,11 (1,04+0,07)

0,77 (0,72 + 0,05)

Органопластик

1,02 (0,66 + 0,36)

0,10 (0,06 + 0,04)

1,13 (0,74 + 0,39)

0,13 (0,08 + 0,05)

1,24 (0,84 + 0,40)

0,16 (0,11 + 0,05)

1,37 (0,97+0,40)

0,23 (0,16 + 0,07)

1,49 (1,19+0,30)

0,42 (0,33 + 0,09)

1,46 (1,23+0,23)

0,55 (0,46 + 0,09)

1,37 (1,21+0,26)

0,67 (0,60 + 0,07)

Напряжения, возникающие на границе вала и балластного диска, — рас­ тягивающие и численно равны напряжениям (6). Эпюра радиальных напряжений монотонна, а окружные напряжения имеют максимум:

 

= 8(3 + Vre(l)) [ ( ‘ + V )

v) ( l + 3v)], (7)

*

l l (1—v)

 

где p -

!

 

Общий размер маховика m (отношение радиуса вала к наружному радиусу обода) теперь уже не варьируется. Полагая его равным 0,1, по­ лучим для приведенного давления на внутренней поверхности обода вы­ ражение

__ пр(п1г 0,01)

Р = 4(3+ V,G(1)) '

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

Результаты анализа энергоемкости маховика в виде обода со сплош­ ным диском — балластом между ободом и валом — представлены в табл. 3. Видно, что для всех композитов максимальная удельная объем­ ная энергоемкость достигается в случае, когда плотности балласта и материала обода близки. При этом ее значения практически совпадают с величинами, достигнутыми при оптимальном использовании «чистого» балласта. Удельная массовая энергоемкость существенно возрастает с уменьшением плотности балластного диска; при минимальных ее значе­ ниях массовая энергоемкость практически совпадает с максимальной величиной, полученной при использовании чистого балласта. Таким об­ разом, центр целесообразно изготавливать из материала с плотностью не выше материала обода. При этом удается сохранить на достаточно высоком уровне обе характеристики энергоемкости.

Выводы. 1. Использование балласта обеспечивает существенный при­ рост удельной объемной энергоемкости маховика лишь в случае изготов­ ления его обода из органопластика — максимальный прирост ~ 50% ; для стекло-, боро- и углепластика максимальное увеличение энергоем­ кости незначительно (от 7 до 14%).

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

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

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

СП И С О К Л И Т Е Р А Т У Р Ы

1.Гулиа Н. В. Маховичные двигатели. М., 1976. 170 с.

2.Портнов Г Г., Кулаков В. Л. Исследование энергоемкости маховиков из компо­

зитов, изготовленных намоткой. — Механика полимеров, 1978, № 1, с. 73—81.

3.Post R. F., Post S. F. Flywheels. — Sci. Amer., 1973, vol. 17, December, p. 19—23.

4.Reedy E. D., Gerstle F. P. Design of spoked-rim composite flywheels. — Proc.

1977 Flywheel

Technol. Symp. San

Francisco, 1977, p. 99— 110.

5. Селезнев Л. H„ Портнов Г. Г. Хордовая намотка дисков лентами из компози­

тов. — Механика полимеров, 1977, № 6, с. 998—

1001.

6. Knight

С. Е., Kelly /. I. jr., Huddleston R. L., Pollard R. E. Development of the

«Bandwrap»

flywheel. — Proc.

1977 Flywheel

Technol. Symp. San Francisco, 1977,

P 7? Chamis С. C., Kiraly L. I. Rim-spoke composite flywheels: Detailed stress and vibration analysis. — Proc. 1975 Flywheel Technol. Symp. Lawrence Livermore Lab.,

1975,

p. 110— 116.

__

8.

Лехницкий С. Г

Анизотропные пластинки. М., 1957. 463 с.

Институт механики полимеров

Поступило в редакцию 20.09.79

АН Латвийской ССР, Рига

 

МЕХАНИКА КОМПОЗИТНЫХ МАТЕРИАЛОВ, 1980, № 2, с. 300—303

УДК 678.02:678.067.5

Л.В. Клычников, С. П. Давтян, Р. А. Турусов, С. И. Худяев,

Н.С. Ениколопян

ВЛИЯНИЕ УПРУГОЙ ОПРАВКИ НА РАСПРЕДЕЛЕНИЕ ОСТАТОЧНЫХ НАПРЯЖЕНИЙ ПРИ ФРОНТАЛЬНОМ ОТВЕРЖДЕНИИ СФЕРИЧЕСКОГО ОБРАЗЦА

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

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

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

Пусть к некоторому моменту времени фронт отверждения продви­ нулся от R] и R2. С лой I твердый, а зона отверждения находится в тон­ ком слое AR = R3— R2. Будем считать, что твердый слой I способен лишь к упругому деформированию с постоянным модулем Юнга Е и коэффи­ циентом Пуассона v. Упругие постоянные оправки Е\, v\. Полагаем, что в слое II реализуются упругие деформации с теми же Е и v, что и в слое

I,

и деформации усадки

еу: гт1= ет1\ еет = еФ1= е©1;

Б гп

= егп + еу; ее11= ефп = в011+ еу.

 

Очевидно, что для

каждого последующего

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

Уравнения равновесия в напряжениях для случая сферической симметрии при отсутствии объемных сил имеют вид:

Рис. 1. Модель отверж­ дения сферы на упругой оправке: фронт иниции­ руется от внутренней по-

dor 2

(аг- а е ) = 0.

(1)

- ^ - + у

Соотношение Коши:

 

 

du

и

(2)

гг у у ;

е0= еф= — ,

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

ег=-^(сг, —2vae); е0 = еф= -i-[ae(l - v ) - v a r].

(3)

Общее решение (1) — (3) для напряжений записывается следующим об­ разом:

г , Cs .

С2

(4)

°Г=С1+7 Г '

ae= C|_2rr

Здесь Сь Сг — константы, определяемые из следующих граничных усло­ вий для слоев I, II и оправки

r = R 0:

r = R i:

r = R 2:

r = R3:

о II .°« О

ar°= ar1; и°=и1; e0°=eeI;

Or1= aru; и1= ии\ ее1= еви + Еу \

агп = 0.

(5)

(6)

17)

(8)

(5)и (8) означают отсутствие усилий на свободной поверхности оправки

ивнешней (подвижной) поверхности отвержденного материала. Первые равенства в (6) и (7) отражают условия непрерывности нормальных уси­ лий на границе слоев, а вторые — непрерывность радиальных перемеще­

ний. Из (5) — (8) с учетом (3), (4) получим:

с

i = _ 2 e ^ _____R3з - 7 ? 2 3 .

1 =

2 еуЕ

Я3* - Я 23

» 3.

 

1 3(1 —V) /?33 —Я/?!3 ’

2

3(1 —V)

«33- W

' ’

г п ___ 2б„£

Rj3- l R i z

 

2еуЕ

R i'-IR ,1

1

3(1 —v)

tf33- W

 

’ 2

~ 3 ( l - v )

tf33- W

* 3 '

 

г П .

Т(у б1^ -(Л з3- е д

 

 

 

 

3 ( l - v ) _____________

 

 

 

'

(1 -Я о 3/Я|3)(Лз3-М ?13) V

b

 

 

 

^ ^ т - ( « з 3- Л 23)

 

 

 

 

C2» =

3(1 —v)_____________ /,

,vr> 3

 

 

( l - / ? o W ) ( * 3 3- W r

j °

 

 

 

 

 

l+i?o3/(2/?i3)- ( l - v i ) - v i +

 

E i

 

 

(2v-l)-^r-

 

 

l - f l o W

 

 

 

 

 

 

 

l+ t f 03/(2tfi3)

(1-Vi) -

'

i

7 Г

 

 

 

 

VI+ Y (V+1) —

 

l - ^ o W

Напряженное состояние произвольно выбранного слоя с номером k\ определяется совокупным взаимодействием как со всеми последующими твердыми слоями, так и со всеми предыдущими, как единой твердой частью сферы. Напряжения от взаимодействия со всеми предыдущими слоями есть:

0гпРед =

2еуЕ

R h f—kRi3

2еуЕ

3(1 - v )

3-X R r

«в пред= - 3 ( Г 7 о х

 

 

X Д*,+15-М?13 '

2г3 /

Напряжения в этом же слое от действия одного из последующих слоев

с номером k2 {k2>k\)

определяются:

 

 

 

2еуЕ

 

,

2гуЕ

Л»г+1з - а д (

,

Я,3 \

т :

X

0|'

= з (1—v)

 

 

 

00

3 (1 —v)

 

 

 

 

 

^л2+13— 3

/

R I3

\

 

 

 

 

X Rk2+id—hR\J \

l+ X ~2^ )

 

 

 

Rk — радиус отвержденного слоя с номером k.

Суммарные напряжения в рассматриваемом слое от совокупного взаи­ модействия как с предыдущими слоями, так и со всеми последующими

с номерами от k x+ 1 до k2 имеют вид:

 

 

 

 

 

 

 

 

 

 

 

 

о hi =

 

 

 

2еуЕ

 

R k f—hRi3

 

 

 

 

 

 

 

 

 

 

сггпред _|_ агпосл =

 

 

 

 

 

-

(

'

1-

- Т

-

)

>

+

з72Гу£-Гх

 

 

 

 

3(1 - v )

 

 

 

 

 

 

 

 

 

Д*1+13-ЯД !3

 

 

Г6

 

 

 

3 (1 — V)

 

 

 

 

 

«2

R,l+i3- R k 3

 

 

 

 

 

 

 

 

 

 

 

 

(

b

- ^

)

 

 

 

 

 

 

 

 

 

 

 

X

LRk+i~^R\3 ;

R h i^ f^ R h i+ u

 

 

(9)

 

 

 

 

 

Ai+1

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

а0А1= аепред+а0посл=

УЕ

 

R kf — hR i3

( l + ^ ± ! i )

+

3(1 —v)

 

-----

 

 

 

 

 

 

 

 

 

/Ощ-i3—KR\6 \

 

 

2r3

/

 

 

2еуЕ [

R i 3 \ R h + i3— R h3

 

 

 

 

 

 

 

 

+- 3(1 —v)l

/\ 1 ~HAi ■2r3 )

£

^ft+i3—hRi3 ;

 

 

 

 

 

 

( 10)

 

 

 

 

 

 

fti+i

 

 

 

 

 

 

 

 

 

 

 

 

 

Обозначим

Ru = R\+kA;

A =

Я - Я

 

k = k\,ki + \,... ,k 2.

 

При

A-*-0,

 

 

 

 

 

 

 

n

 

 

 

 

 

 

 

 

 

 

 

 

(n —>-oo)

полагая ^i + &iA = r,

R x+ k2k = R

и

переходя

к интегралам

в (9), (10):

 

. ч

2еуЕ

/

R I3

\

R

 

dga

 

 

 

 

 

 

2еуЕ

 

 

f

 

 

 

 

 

 

 

х

° ' (г) =

з ( Т ^ Г

 

 

 

)

J

 

 

 

 

Оо(г)

 

3(1 —v)

 

 

 

2

\

2гЗ

/

;

|3- я « !

 

 

 

 

 

 

 

 

 

 

получаем окончательно выражения

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

, ,

2еу£

/,

 

 

Я,3

\,

/

 

Я3- « , 3

\

 

 

 

 

 

0,(г)

3(1 - V )

\

 

г3

) П\

 

r 3-

w

 

/ ;

 

 

(

.

2еу£

Г / ,

, ,

Я 13

V

/

tf3- W

\

3

1

 

Rl < r ^ R ,

0е(Г)

3(1 - v )

1 \ 1+Л

2г3

/;П\

г3 — XRi3_ /

 

2

J

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

( 12)

которые дают о> и ао в точке г, возникающие в сфере при фронтальном отверждении ее изнутри с учетом упругой оправки. R — граница фронта отверждения.

Если оправка отсутствует, т. е. Ro->R\, тогда Л-И и полученные выше формулы (11), (12) приобретают вид:

 

,

 

 

2ЪуЕ

/

 

Я,3

\

,

/ Я3-/?,3

\

 

 

 

3(1 - v) \

 

г3

/

П \

r3- ^ i 3

/

/ ч

2еуЕ

Г

/

, ,

Я,3

\

/

R3- R x* \

3

1

00^

3(1 —v)

L

\

2г3

)

П (

г3 —Я13 )

~2 J ’

Напряжения в самой оп­ равке являются результатом

совокупного

действия на

нее

всех

образованных

слоев

от­

вержденного материала:

 

 

. .

уЕ

 

 

<rr™ >W =T ^ vT (1- X ) X

(\-Ro4r3) .

R3- X R t3

 

(l - R o 3IRi3)

(1—

 

g6°™(r)=

2еу£\

(1 - Я )Х

 

 

с5(1 - v )

 

 

 

о - е д а д ,

^ 3- w

.

х

( 1 - а д ^ 1 3)4

о - а д

3 ’

Рис. 2. Кривые распределения радиальных (/)

и

окружных

(2)

напряжений в слое полимера

и

оправке:

£i =

2- 104 кгс/мм2; Vi = 0,27; £ 0 =

= 47 мм; R\ = 50 мм; £ = 2 00 кгс/мм2; v = 0,32; £ = 100 (а), 55 мм (б).

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

Заметим, что в оправке появляются сжимающие напряжения (агопр<С0, с70опр<О ), а слой полимера будет сжат по радиусу (аг< 0 ) и растянут в окружных направлениях (ае>0).

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

 

 

 

С П И С О К Л И Т Е Р А Т У Р Ы

 

 

 

1.

Турусов

Р. А.,

Давтян С. П., Шкадинский К. Г.,

Розенберг

Б. А.,

А ндреев­

ская Г

Ф., Ениколопян Н. С. Механические явления в условиях распространения фронта

отверждения. — Докл. АН СССР,

1979, т. 247, № 1, с. 97.

Розенберг

Б. А.,

Ениколо­

2.

Давтян

С. П.,

Арутюнян

X. А., Шкадинский К. Г.,

пян Н. С. Закономерности отверждения эпоксидных олигомеров диаминами в условиях распространения фронта реакции. — Высокомолекулярные соединения. Сер. А, 1977,

т.19, № 12, с. 2726.

3.Сурков Н. Ф., Давтян С. П., Розенберг Б. А., Ениколопян Н. С. Расчет стацио­

нарной скорости фронта при отверждении эпоксидных олигомеров диаминами. — Докл. АН СССР, 1976, т. 228, № 1, с. 141.

Отделение Института химической физики АН СССР,

Поступило в редакцию 27.11.79

Черноголовка

 

МЕХАНИКА КОМПОЗИТНЫХ МАТЕРИАЛОВ, 1980, № 2, с. 304—307

УДК 678.2:678.067.5

А. И. Цаплин, С. В. Бочкарев

НЕСТАЦИОНАРНОЕ ТЕМПЕРАТУРНОЕ ПОЛЕ ВРАЩАЮЩЕГОСЯ СТЕКЛОПЛАСТИКОВОГО ЦИЛИНДРА, ОБЛУЧАЕМОГО ПАРАЛЛЕЛЬНЫМ ПУЧКОМ ЭЛЕКТРОНОВ

В последнее время в технологическом процессе изготовления изделий цилиндрической формы из полимерных композитных материалов пред­ лагается использовать энергию ускоренных электронов [1], которая по­ зволяет интенсифицировать процесс отверждения. Направленный пучок электронов приводит к разогреву материала, максимальная температура при этом может достигать критического значения, при превышении кото­ рого в полимерном материале происходят деструктивные изменения, свя­ занные с интенсивным газовыделением [2]. Образующиеся пересыщенные газовые растворы вызывают ухудшение механических свойств композита. Кроме того, наоднородное температурное поле приводит к появлению температурных напряжений [3]. В результате возникает задача матема­ тического описания тепловых явлений в композите, облучаемом ускорен­ ными электронами.

1. Рассматривается бесконечный полый цилиндр с заданными внут­ ренним (п) и внешним (г2) радиусами, вращающийся с постоянной уг­ ловой скоростью ш, облучаемый параллельным пучком электронов за­ данной ширины b (рис. 1). Мощность теплового источника Q может ре­ гулироваться энергией электронов Е и током пучка I. Перенос тепловой энергии в цилиндре описывается уравнением, общий вид которого можно найти в [4]. Применительно к рассматриваемой задаче нестационарное уравнение переноса энергии имеет вид:

рс

dt

dt \

I

dt \

X dt

1

д

( 1)

(

 

дг

/ г

дг

г2

д ф

 

 

 

где t — температура; т — время; X — коэффициент теплопроводности; р — плотность; с — теплоемкость; г, ср — радиальная и окружная коор­ динаты. Тепловой источник действует в зоне облучения, характеризуе­ мой углом ф:

Q _ j EI

фо^ф^фо + ф;

I 0

9 < ф о ; ф > ф о + Ф-

Рис. 1. Расчетная схема.

Переходя к логарифмической переменной z = lnr и представляя уравнение (1) к виду, который оно имело бы при постоянном значении коэффициента теплопроводности, получим:

dt

, dt

X*

(

дЧ

дЧ

\

Q

 

~ -Г

J 3

-----------Ту-----

\

2~ + ~ Д Т ' /

/------------

рс

> (2 )

дх

Оф

e2zpс

\

dz2

уф2

 

где коэффициент эффективной теплопроводности

dt дХ dt дХ dz dz дф дф

Х* = Х +

дЧ ( дЧ dz2 ' дф2

Уравнение (2) решается при краевых условиях, включающих началь­ ное распределение температуры

t (т = 0) = t0

и граничные условия третьего рода на внутренней

f=z, —

(/п tс)

(3 )

и наружной поверхностях цилиндра

z = z„

( a h2 + ОС л ) ( / п —tс ) ,

(4)

где осл — коэффициент лучистой теплоотдачи, определяющийся по за­ кону Стефана—Больцмана и зависящий от температур поверхности /п, окружающей среды tc и степени черноты е ссл = 5,67 •10-8е(/п + /с+ + 546)[(/п + 273)2\+ (/с+ 273)2]. Коэффициенты конвективной теплоотдачи a/ии а/(2 считаются заданными.

2. Для численного решения задачи область покрывается регулярной сеткой

7=1, 2,

,N + 1; hz = ~

In —

/ С 2,

,

hq>= дд

 

N

г\

 

 

 

А= 1,2,.

где N, М, L — числа разбиений области соответственно в направлениях радиальной и окружной координат, а также число временных интерва­ лов. Разрешающие соотношения для определения дискретных значений температур в узловых точках определяются на основе экономичной факторизованной разностной схемы [5]. Применительно к уравнению (2) наряду с двумя временными слоями к и /г+1 вводится промежуточный слой k-\-l/2. Обозначая температуру на слоях k, k + l/2, /г+1 соответ­ ственно t, t, i, запишем конечно-разностную аппроксимацию уравне­ ния (2):

0,5/1т + ©

t i , i + 1— /г'Д— 1

^»*г,j

/ /г+ l.j

2 / f ,j+

2hv

e2zi(pc)i,j

'

/*z2

.

ti,i—\ 2/j,j +

\ _|___ Q .

(5)

 

V

'

(рФ.З

 

 

+ co

7i,j+l — /г,j—1 _

e2zi(p6’)i,j

\

2/i,j + /j-l,j

' 0,5/it

2ЛФ

/1г2

 

2/i,j + /i,j+I \ i

Q

( 6)

 

Лф2

'

(Pc)i,j

 

 

Полученные соотношения (5), (6) образуют для всех внутренних узло­ вых точек (£= 1 ,2 ,..., Л7; /=1,2, .. . , М —1) системы алгебраических урав­ нений с трехдиагональными матрицами, решение которых осуществля­ ется методом покоординатной прогонки.

Для граничных узловых точек значения температур определяются из условий (3), (4), записанных в конечных разностях:

ГС+ fiVj'fljV+l,

 

/c + /2,jal,j _

/jv+и

a N + \ ,j :

■i,j

1 +fljv+l,j

(алг + огл

1 ■+ al,i

ai'3 0Lh\hzez'