CC BY
29
Раздел VIII. Естественные науки
УДК 534.4
А.И. Жорник, В.А. Жорник, П.А. Савочка
ВЛИЯНИЕ ЭНЕРГИИ РАЗРУШЕНИЯ НА РАЗВИТИЕ КОЛЬЦЕВОЙ ТРЕЩИНЫ В ЦИЛИНДРЕ ПРИ ЕГО ОХЛАЖДЕНИИ
Рассчитываются температурные поля и термоупругие и напряжения в сплошном цилиндре с внешней кольцевой трещиной, вызванные охлаждением цилиндра и источником теплоты, который моделирует мгновенное выделение энергии на берегах трещины за счет ее роста при охлаждении цилиндра. Показано, что величины температурных полей и напряжений, вызванных источником теплоты, ничтожно малы по сравнению с термоупругими напряжениями, связанными с охлаждением цилиндра и вызывающими мгновенный рост кольцевой трещины.
Цилиндр; неорганическое стекло; кольцевая трещина; мгновенный источник теплоты; охлаждение; термоупругие напряжения.
A.I. Zhornik, V.A. Zhornik, P.A. Savochka FRACTURE ENERGY INFLUENCE ON RING-SHAPED CRACK DEVELOPMENT IN A CYLINDER UNDER COOLING
Temperature fields and thermoelastic stresses in a solid cylinder with an outer ring-shaped crack caused by cylinder cooling and heat source which models an instant energy liberation on crack banks because of crack's growth at cylinder cooling are computed. It is shown that the values of these temperature fields and stresses are negligibly small as compared with thermoelastic stresses connected with cylinder cooling and causing instant growth of the ringshaped crack.
Cylinder; inorganic glass; ring-shaped crack; instantaneous heat source; cooling; thermoelastic stress.
Введение. В процессе нагрева или охлаждения различных элементов конструкций в них возникают нестационарные градиенты температур, которые обусловливают температурные напряжения. Эти напряжения, усиливаясь вблизи трещиноподобных дефектов, могут вызывать разрушение изделия.
В данной работе рассматривается задача об определении температурных напряжений в сплошном цилиндре с внешней поперечной соосной кольцевой трещиной, расположенной в центре цилиндра, термонапряженное, состояние которого связано как с процессом охлаждения нагретого цилиндра внешней средой, так и с процессом его разрушения, вызванного распространением в нем кольцевой трещины.
Постановка задачи. В работе авторов [1] исследовалось напряженное состояние вблизи кончика внешней поперечной соосной кольцевой трещины радиуса rb, расположенной в середине сплошного цилиндра радиуса rc, длины I, нагретого до постоянной температуры T0 и охлаждаемого путем теплообмена в воде температуры в. Причем теплообмен торцов цилиндра со средой происходил симметрично с коэффициентом теплообмена ао, а цилиндрической поверхности с коэффици-
ентом теплообмена аь Окончательное решение задачи свелось к нахождению зависимости коэффициента интенсивности напряжений (КИН) К, управляющего ростом кольцевой трещины, от времени t для различных размеров трещины (радиусов фронтов трещины гъ).
Рассматривался конкретный случай охлаждения, когда условия охлаждения цилиндра были таковы, что интенсивность теплообмена между торцами и средой была весьма велика. Кроме того, длина цилиндра была значительно больше радиуса (Игс = 25), а исследуемый процесс кратковременный. Поэтому влияние охлаждения торцов на развитие трещины было незначительным. В связи с этим коэффициент О) был взят неограниченным (граничное условие первого рода). Теплообмен между цилиндрической поверхностью и средой задан был критерием Био
Бг = °- г = 90.
На рис. 1 приведены графикы зависимости безразмерного КИН К* от безразмерного времени Ев при различных размерах кольцевой трещины а для этого случая охлаждения.
Рис. 1. Зависимость КИН К от времени Ев для различных размеров трещины а
На этом рисунке
К. Щ (1- у) ^ (|)
' аЕ (Т„- в) ,
где КI - КИН; Е - модуль Юнга; V - коэффициент Пуассона; ат - коэффициент термического расширения материала цилиндра; ¥в = — - критерий Фурье;
Гс
а =------- - температуропроводность; Лт - теплопроводность; ру - плотность;
Русу
су - удельная теплоемкость материала цилиндра; а = гъ/гс - относительный размер фронта кольцевой трещины, К1С - трещиностойкость (постоянная материала), при достижении которой КИНом КI трещина начинает распространяться [2]. Нужно
ТУ'*
заметить, что согласно (1) безразмерная трещиностойкость К1С не является постоянной материала, а зависит еще и от разности температур Т0 - в Это означает,
что разрушение цилиндра зависит не только от материала, но и от разности начальной температуры нагрева цилиндра Т0 и температуры охлаждающей среды в.
Из рисунка видно, что трещина размера а = 0,95 достигает трещиностойко-
Т/' * т/' *
сти К/с = 0,31 (максимальное значение К7) в момент времени Ев= 0,015 (точка
а на рисунке), и далее растет скачком (вертикальные стрелки) до значения а =
0,85. Все более крупные трещины (с меньшим радиусом гъ) на рисунке, обозначенные штриховыми линиями, в этот момент времени имеют меньший К1 . Следовательно, в области разрушения 0,85 < а < 0,95 мгновенно выделяется некоторое количество теплоты [3].
Цель работы. Целью работы является расчет температурного поля и соответствующих ему температурных напряжений, вызванных процессом разрушения, а также оценка влияния их на температурное поле и температурные напряжения, вызванные охлаждением цилиндра средой, и, следовательно, на процесс распространения трещины.
Теоретическое обоснование. Рассмотрим задачу о напряженном состоянии конечного сплошного цилиндра, описанного выше, за исключением того, что внутри цилиндра в плоскости кольцевой трещины в пределах га < г < гъ, 2 = 1/2 в момент времени Ев' >0 возникает мгновенный источник теплоты q(г, 2, () - количество теплоты, выделяемой мгновенно единицей площади
д(г, г, і) = д\%(г - га)-х(г - гь )]^| г -1 '\8(Ро - Ро'),
(2)
где х(г) - функция Хевисайда, 5(х) - 5-функция Дирака.
Решение для температурного поля в этом случае получено с помощью интегральных преобразований Ханкеля по г и Фурье по г и имеет вид:
Т* (г, £ /2, Ро) =
Т (г, £ /2, Ро )-в
= кк Ак (Ро У
(
. V
V ’о Віс
Л
^к А-к • А-к
- ООЭ —— + 81П
2
2
. (3)
Радиальные, окружные, осевые и касательные напряжения в плоскости расположения трещины (2 = 1/2), которые вызываются этим температурным полем, имеют вид:
О-* {г 1 Рз)-агг(г,£/2,Ро)(1 -у)_ уу Ап— (ро)Уп
ггV , 2, 1 атЕТ-в) к к 2 2 г2
уп + Хк ^2
г
хІ7уі
г
Уп-V го J
+
ХЛ_ т
2 У С
Уп£
г
\
Уп~
V го J
+
+
У1 (Уп )
Г°хкІі\^ хкГ-° | + Вії с \хкГо
(і - 2у)ї с I Х—Г- 1 + Хк-£І11 ХкГ-
+
Ві
—
1 + 2(1 -уУ
2 2 Хк го
i:\xk
п=1 к=1
X
2
Г
г
2
о
Х
с
к
е
£
+
2^(V-1)-х; '-р у01 X,
( Ві'2
\ х.г
V к с
2 + хк
2(1 -у)к I
*7)-“- 1іI х^
') г X, V '
хк хк . хк -соя — + мп —
Віп
(4)
* ( ' я 1 СТ I Г —, ГО I =
М 2 )
Сw(г,'/2,Го)(1-у) _ ^ Апк(Го)у
атЕ(Т0 -в)
2
Уп + хк^2
2 2 Л
Уп +
ХкГс
2
п' )
Уп '
Л
' И
Уп“
V гс)
г
-
г
Уп-
V Гс )
+
+ -
Л1 (Уп ) _ х‘ГіІ' (^ )+ВіІ 0 (хкгг 1 (1 - ^ 0(хк- ] + ± - X
Ві 10 (хк7- ]" Г, , 2(1 -V)'2 ] хІК \ І2 (хк 7)
Ві'
( 2у- 1)-
\
'2
Ві -2 + х,2
V Гс )
1о\ хк-
хкВі +
(Ві -2
V хк гс
2 + Хк
2(1 -V)
/,|хк| №(хк -
хк хк хк
к - СОЭ — + 81П —
Ві
(5)
С
г— Гп1 С(г,' / 2, Го)(1 -у) = -УУ Апк (Го)Уп \ У Л
Vг,2, ] аТЕ(Т0 -в) , 9 г2 1УЛ0
п=1 к=1 2 . 2 ^_с_
Уп + Хк ^ 2
' И
Уп-
V - )
+
Л1 (Уп )
2(у-2)10( хк- )-Хк-Іі1 хк-
+
Ві
12 [ хк -
1 + 2(1 -V)
2 2 хк гс
І2 хк
-
х —
X
к
і
г
с
X
о
I
2
2
X
X
г
г
с
X
I
і
г
с
2
2
*
Г
е
Ві (2у- 1)-
Ві
т+хк
І
х,г
—
2
-
с
с
2
£
і
с
Ві -
х, г
V к с
2 + хк
2(1 -V)
І І |хк
СОБ —— + Б1П —
V Ві0 2 2 )
(6)
* і 'і—' ) с (г,' / 2,^о)(і-V) А, (^о)у Г г ( г ,
С_| г-,¥о I = / лч ^ = -ЕЕ п п 1- ^х^іI Уп 7 1 +
2
а^(^ - в)
у + х .-Ст
у п к £ 2
і
£
+
гс Л1 (Уп )
хкІ1 хк гг)+Ві 7101 хкг
+
1+
2'2 (1 -V)
22
і' кт
+•
Ві — ( 2у- 1)-
2
Л
2
Ві т + хк
V гс )
х
х^Ві +
V хк гс
2 + хк
2(1 -V)
где Апк(Ео) имеет вид:
Апк (Гп )= 4Ві 2е
\Ро (У
х,.
х
Ві
Б1П — + СОБ -2
V Ві0
б1п х*. - СОБ х^ + 1
+
(7)
+
дх(і'° - Го !) _з_<
'(Т0 -в)РгСК Л1 (Уп )
У;7 +х42^ IГо'
Г
х
Л
к
хк
хк
СОБ — + Б1П
V Ві0 2 2 )
-Л
Уп
г
-
с V с ) с
Уп
V гс )
2 Л
А
1 хк
1 + —
V Ві 0 )
х2 1 -х
V Ві 0 )
БІП 2 хк 2
2
2х
+--------Б1П х
Ві,
0
х
(8)
где хк - корни уравнения
г
х
к
£
£
£
х
х
х
п
г
х
к
і
г
с
0
£
£
х
х
к
к
2
2 . 2 г
Уп +хк~2
і
X
г
гь
г
а
Ь
X
X
X
X =
X2 - ВІІ 2 хБіп
(9)
уп - корни уравнения
Уи 1 Ы = ви 0 (у); (10)
а а0
в1 = — г , в10 = — I - критерии Био; аь ао - коэффициенты теплообмена ме-
А А
жду цилиндрической поверхностью, торцами и средой. Предполагается, что теплообмен между торцами цилиндра и средой симметричный.
Численный расчет. Численный расчет температурного поля и осевых температурных напряжений проводился по (3) и (6) на ЭВМ в программе МаШсаЛ При этом бралось 100 членов ряда, поскольку при малых временах, которые имели место в задаче, необходимо учитывать их большое число. Расчет проводился только для осевых температурных напряжений, что связано с симметрией задачи относительно плоскости кольцевой трещины и с гипотезой о том, что рост трещины вызывают максимальные нормальные к ее берегам напряжения. Этими напряжениями для кольцевой поперечной трещины являются осевые напряжения. Расчет проводился для цилиндра из неорганического стекла (А = 0,5852 Вт/м-К, ру = 3,03-103 кг/м3, су = 8,33-102 Дж/кг-К, Е = 6-104 МПа, V = 0,23, а = 87-10-7 1/К, КС = 0,71 МПа-м1/2) радиусом гс = 2-10-3 м и длиной I = 5-10-2 м. Коэффициент теплообмена согласно [4] а! = 26-103 Вт/м2К.
На рис. 2 для случая Q =
Я
Кт0 -0)ру
= 10-
приведено распределение
температуры по поперечному сечению цилиндра (р = г/гс) в плоскости расположения трещины (х/1 = 0,5) в различные моменты времени Ев (сплошные линии) и для Q = 0 (штриховые линии).
Рис. 2. Распределение температурного поля по сечению цилиндра в плоскости расположения трещины в различные моменты времени Ев
Из этого рисунка видно, что мгновенный источник Q = 10" в момент времени Ев' = 0,015 резко поднимает температуру на берегах мгновенно проросшей трещины (0,85 < а < 0,95, т/1 = 0,5), создавая на границах кольцевого источника тепла большие градиенты температур, в то время как отсутствие такого источника (штриховые линии) не вызывает таких градиентов.
2
На рис. 3 изображено распределение осевых температурных напряжений по сечению цилиндра р = г/гс в плоскости расположения кольцевой трещины в различные моменты времени Ев, рассчитанные по (6) для Q = 0 (штриховые линии) и Q =10-2 (сплошные линии).
‘2 ----------------------------------------------------------------------------------------------------------------------------------------------- о
0 0.2 0.4 0.6 0.8 1 о
Рис. 3. Распределение осевых напряжений по сечению цилиндра в плоскости расположения трещины в различные моменты времени Ев
Из этого рисунка видно, что в области действия источника интенсивности Q = 10-2 возникают термоупругиие напряжения, причем вблизи границ источника напряжения разных знаков. Внутри источника за счет нагрева - напряжения сжатия (отрицательные), а вне границ в более холодной области - напряжения растяжения (положительные). Для отсутствующего источника Q = 0 таких напряжений нет.
ТУ*
Анализ зависимостей К от времени Ев для различных размеров трещины а при интенсивности Q = 10-2, аналогичный анализу зависимостей, приведенных на рис. 1 для Q = 0, в данной работе не проводится. По-видимому, эти зависимости будут гораздо более сложными, чем изображенные на рис. 1, и, следовательно, трещина будет расти гораздо более сложным образом. Тем не менее проведем оценку интенсивности источника Q для случая прорастания трещины в неорганическом стекле, который отражает процесс роста трещины, приведенный на рис. 1. Согласно данным [3] максимальная энергия разрушения для неорганического стекла равна д = 7,5 Дж/м2. Считая, что вся энергия разрушения превращается в теплоту,
„ д
найдем интенсивность источника у = ^----------------г------ с учетом физико-
~ в)рУ°У
механических постоянных цилиндра и его размеров, приведенных выше. Получаем, 5 1
что Q = 5,9-10 -------. Величину Т0 - в найдем из (1) и указанных выше физико-
70 ~в
ТУ*
механических постоянных с учетом, что К1С = 0,31. В этом случае Т0 - в = 86 К. Тогда Q = 6,9-10-7.
Расчет распределения температур и осевых температурных напряжений в плоскости расположения трещины для Q = 6,9-10-7 по (3) и (6) показывает, что эти распределения совпадают со значением Q = 0.
Заключение. Из проведенного выше анализа следует, что влияние энергии разрушения на температурное поле и на вызванные этим полем температурные напряжения, возникающие в результате охлаждения цилиндра средой, ничтожно
мало. Однако действием источников тепла, обусловленных другими причинами (прохождение электрического тока по проводникам, ядерные и химические реакции и др.), пренебрегать в общем случае нельзя. Для этого необходимо использовать результаты, полученные в данной работе.
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Жорник В.А., Савочка П.А. Хрупкое разрушение цилиндрических тел при нестационар-
ном тепловом воздействии // Известия ЮФУ. Технические науки. - Таганрог: Изд-во ТТИ ЮФУ, 2009. - № 4. - С. 205-211.
2. Irwin G.R. Analysis of stresses and strains near the end of a crack traversing a plate //
J. Appl.Mech., 1957. - Vol. 24, № 3. - P. 361-364.
3. ShandE. // J. Amer. Cer. Soc. - 1961. - Vol. 44. - P. 21.
4. Кутателадзе С.С. Основы теории теплообмена. - Новосибирск: Наука, Сиб. отделение,
1970. - 659 с.
Статью рекомендовал к опубликованию д.ф-м.н., профессор Г.В. Куповых. Жорник Александр Иванович
ГОУ ВПО "Таганрогский государственный педагогический институт".
E-mail: [email protected] г. Таганрог, ул. Инициативная, 48.
Тел.: 88634601807.
Жорник Виктория Александровна
E-mail: [email protected].
Савочка Петр Анатольевич
Технологический институт федерального государственного автономного образовательного учреждения высшего профессионального образования "Южный федеральный университет" в г. Таганроге.
E-mail: [email protected].
347928, г. Таганрог, пер. Некрасовский, 44.
Тел.: +78634371603.
Zhornik Aleksandr Ivanovich
Taganrog State Pedagogical Institute.
E-mail: [email protected]
48, Initsiativnaya Street, Russia, Taganrog, Russia.
Phone: +78634601807.
Zhornik Viktoriya Aleksandrovna
E-mail: [email protected].
Savochka Petr Anatolievich
Taganrog Institute of Technology - Federal State-Owned Autonomy Educational Establishment of Higher Vocational Education “Southern Federal University”.
E-mail: [email protected].
44, Nekrasovskiy, Taganrog, 347928, Russia.
Phone: +78634371603.
