Механика композитных материалов 2 1980
..pdfуже отмечалось, связан с высокой анизотропией свойств этого материала. Однако этот же материал несет и самые большие потери в удельной мас совой энергоемкости. Плотность балласта в рассматриваемой конструк ции маховика при оптимальных размерах обода и балласта незначи тельно влияет на энергоемкость: с уменьшением плотности балласта объ емная энергоемкость маховика несколько снижается, а массовая — возрастает. Энергия накапливается в основном ободом маховика; вклад балласта в энергоемкость маховика существен лишь у органопластика.
Для окончательного решения вопроса о целесообразности и возмож ности практического применения балласта необходимо оценить абсолют ную величину оптимального давления р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посл= |
2еУЕ |
|
R kf — hR i3 |
( l + ^ ± ! i ) |
+ |
|||||||||||||||
3(1 —v) |
|
----- |
|
|
|
|||||||||||||||
|
|
|
|
|
|
/Ощ-i3—KR\6 \ |
|
|
2r3 |
/ |
|
|||||||||
|
2еуЕ [ |
R i 3 \ V» 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 ’ |
Напряжения в самой оп равке являются результатом
совокупного |
действия на |
нее |
||
всех |
образованных |
слоев |
от |
|
вержденного материала: |
|
|||
|
. . |
2еуЕ |
|
|
<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, |
, |
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'