К теории флуктуаций роста трещин в теплоэнергетическом оборудовании ТЭС и АЭС

УДК 621.18.021

К теории флуктуаций роста трещин в теплоэнергетическом оборудовании ТЭС и АЭС

Семенов В.К., д-р техн. наук, Беляков А.А., асп., Щебнев В.С., канд. техн. наук.

Предлагается стохастическая математическая модель роста трещин в теплоэнергетическом оборудовании, основанная на уравнении Фоккера-Планка и полуэмпирическом уравнении роста средней длины трещины.

Ключевые слова: стохастическая модель, рост трещины, прогнозирование.

To the theory of fluctuations of growth of cracks in heat power equipment heat power plant and the atomic power station

Semenov V.K., doctor of science, Belyakov A.A. post-graduate student, Shhebnev V.S., candidate of science

In work the stochastic mathematical model of growth of cracks in the heat power equipment is offered. The model is based on the Fokker-Planck equation and the semiempirical equation of growth of average length of a crack.

Keywords: stochastic model, crack growth, forecasting.

Оборудование ТЭС и АЭС работает под постоянным воздействием изменяющейся температуры, давления, циклических нагрузок и пр. Такие условия работы приводят к накоплению в материале оборудования значительных повреждений и развитию макроскопических трещин. Нередко зародыши трещин находятся в материале еще до ввода его в эксплуатацию [1].

Растрескивание стенок теплообменных трубок и других элементов оборудования определяется целым комплексом условий (накопление продуктов коррозии, тепловые и динамические режимы работы, механические воздействия и пр.), многие из которых являются неконтролируемыми. Следовательно, на процесс накопления дефектов следует смотреть как на стохастический и исходить из вероятностных представлений. Многие из величин можно считать определенными, детерминированными, тогда как сам процесс роста трещины - случайный и длина трещины -случайная величина, изменяющаяся во времени непрерывным образом. Причем если интервал времени мал, то и изменение длины трещины тоже мало, т.е. за малое время наиболее вероятны переходы, в результате которых длина трещины изменяется незначительно, а большие изменения длины трещины маловероятны. В таких условиях эволюцию системы можно рассматривать как непрерывный стохастический процесс марковского типа [2]. Это означает, что вероятность перехода системы из одного состояния (начального) в другое (конечное) зависит только от начального состояния и не зависит от тех состояний, которые предшествовали начальному. Последнее предположение является весьма общим и, хотя вначале не может быть доказано, получает обоснование в дальнейшем. Так, будет показано, что уравнения, описывающие кинетику роста тре-

щины на детерминированном уровне, являются по существу следствием сделанного допущения.

Статистическое описание роста трещины будем осуществлять заданием функции распределения вероятностей р(/Д причем р(/,/)с(/ представляет собой вероятность того, что данный элемент имеет трещину, длина которой находится в интервале от / до / + б/. Рассматриваемый нами стохастический процесс является непрерывным по обеим переменным t и /. Такие стохастические процессы подчиняются уравнению Фоккера-Планка [2]

^ -I [ A )p(l. I )]+|, dt dl dl

[B(l )p(l, t)].

Кинетические коэффициенты Л(Г) и В(/) соответственно представляют собой среднее и среднеквадратичное изменение длины трещины за единицу времени:

да да

А(/) = | ц<о(/, ц)ёц, В(/) = -1 ц 2ю(/, ц )ёц.

0 0

Здесь д - длина зоны пластической деформации, возникающей перед вершиной трещины; (/,д)бд -вероятность изменения длины трещины в интервале от / до / + д за единицу времени. Как известно [1], рост трещины связан с развитием зоны пластической деформации перед вершиной трещины, причем д// <<1. Этот процесс носит пороговый характер, т.е. величина д должна быть больше или равна некоторому минимальному значению 5. Согласно сделанным предположениям, функция ю(/,д) быстро убывает с увеличением д, так что основную роль играют малые значения д, т.е. 5. Так как 5 > 0, то рассматриваемый нами стохастический процесс является процессом чистого рождения. В этом случае кинетические коэффициенты

1

связаны простым соотношением В(/) = —5А(/). Учи-

тывая, что 5 является единственным геометрическим параметром задачи развития трещины в бесконечной среде, длину l будем измерять в безразмерных единицах 1/5. При такой нормировке коэффициент A(l) в уравнении Фоккера-Планка умножится на 1/5, а коэффициент B(l) - на 1/52. Нормированные таким образом коэффициенты будут

1

связаны соотношением B(l) = — A(l). Чтобы не вводить новых обозначений, далее сохраним за l, A и В старые обозначения.

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

Сначала найдем уравнение для средней длины трещины. При этом анализ проведем в общем виде, считая коэффициенты А и В разными. С этой целью умножим левую и правую части уравнения (1) на l и проинтегрируем по всевозможным значениям длины трещины от l0 до да. Имея в виду, что при l = l0 и l = да функция распределения p(l, t) = 0 , после несложных преобразований получим

d~dr 1 A

7l0

Индекс l0 означает, что усреднение ведется при фиксированном значении начальной длины трещины. Усредняя по начальным длинам, получим уравнение

=< A

dt

Уравнение (3) имеет очевидный физический смысл: скорость роста трещины определяется средней скоростью потока роста трещин. Однако для того, чтобы им воспользоваться, нужно перейти от среднего потока к потоку от среднего. Для этого разложим A(l) в ряд Тейлора вблизи <l>:

A(l) * A(< l >) + dA(< l >) (l -<l>) + dl

1 d2A(< l >).. . . 2

+-------^—-(l- < l >)2.

2 dl2

Усредняя по l, получим

A(l) * A(< l >) +1 d A(<21 >) A.

2 dl2

При малой дисперсии распределения А =< (/-< / >)2 > << < /2 > вторым слагаемым

можно пренебречь:

« А

На детерминированном уровне описания именно это уравнение должно лежать в основе всех теорий роста трещин в аппаратуре. Заметим,

что в случае линейной зависимости А(/) это уравнение становится точным.

Выведем теперь уравнение для дисперсии распределения:

А =<

Из уравнения (3) имеем

d < l >2 dt

= 2 < l > A(l).

Умножая все члены уравнения (1) на / и интегрируя по /, с учетом (7) и (8) получим

= 2 [< /А

Раскладывая </А(/)> в ряд Тейлора и усредняя по /, найдем

< lA(l) >=< l > A < l > +

[lA(l )k

<l >

•A.

2 6/2

С учетом (5) и (9) уравнение для дисперсии распределения запишется в виде

dA = ( dA

dt ~ I dl

• A + 2B(< l >).

Если А(< / >), В(< / >) явно от времени не зависят (или зависят одинаковым образом), то вместо времени t можно ввести новую переменную:

1ГГ = АТб-[|п А2

6 < / > 6 < / >- -I А

Полученное уравнение является линейным и интегрируется в квадратурах:

A =

A

A

<l >

Ao

dl

Здесь учтено, что в принятых переменных 1

В(/) = — А(/). Аналогичным образом можно найти

уравнения для последующих моментов распределения.

Различные полуэмпирические уравнения роста трещин (уравнения Пэриса, Эрдогана-Формана и др.) основаны на теории Гриффитса, в основу которой положен закон сохранения энергии. Энергия, необходимая для развития трещины, равна энергии упругости, заключенной в объеме трещины. Уравнения роста трещин в бесконечной среде при непрерывной и циклической нагрузке можно записать в следующем виде [1]:

d < l > оп

= а-

dt

Kc

m/2

d < l >=p(Mm <l >m/2

dt

Kc

Здесь а и р - эмпирические коэффициенты пропорциональности, учитывающие геометрию элемента и внешние условия; ст - напряжение в материале; Кс =^2уБ - критический коэффициент интенсивности напряжений, при котором распространение трещины происходило бы со скоростью

упругой волны ~ ст ; Е - модуль Юнга; у - удель-

ная энергия пластического деформирования; т -продолжительность цикла. Как показали эксперименты, эмпирический коэффициент т для большинства материалов и сплавов может принимать значения от 2 до 6 [1]. В частности, для углеродистых сталей т = 4.

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

6/

&

-а(1)(1 - /*)" .

Здесь / * = </>/1 - новая переменная; / - предельная длина трещины; п = т/2; а (0 - эмпирический коэффициент, учитывающий указанные выше свойства материала, условия эксплуатации и его старение:

а

Входящие в (17) коэффициенты должны определяться на основе регрессионного анализа результатов натурного эксперимента.

Для п = 1 рост средней длины трещины определяется экспоненциальным законом / *

Для п > 1 рост трещины происходит по закону

1

/ Ч?) = 1--------!-----т.

[1 + (п - 1)!Ц)]п-1

Здесь принято обозначение Ц?) = |а(?)М .

о

С учетом (16) дисперсия распределения определяется выражением

1

А Р й = -?-

5L 2п -1

а й - безраз-

Здесь Ар - размерная дисперсия, мерная величина.

На рис. 1, 2 представлены графики зависимостей дисперсии распределения от длины трещины при различных значениях числа п.

Анализ графиков показывает, что с ростом числа п максимум дисперсии сдвигается влево, а дисперсия в максимуме уменьшается.

В качестве примера на рис. 3 приведена кривая прогноза роста трещины из тепловой канавки ротора турбины, изготовленного ПО ЛМЗ из стали Р2 [3], при циклической нагрузке. Коэффициенты идентификации функции прогноза определены на основе регрессионного анализа экспериментальных данных по методу наименьших квадратов. Предельная длина трещины составляет 150 мм, а время измерено сотнями циклов.

Рис. 1. Зависимость дисперсии распределения й = А5// от длины трещины при п = 1

Рис. 2. Зависимость дисперсии распределения й = А5// от длины трещины при п = 2

Рис. 3. Функция прогноза роста трещины из тепловой канавки ротора турбины:  - данные эксперимента; -

функция прогноза

Итак, на основании уравнений (6) и (13) можно спрогнозировать среднюю длину трещины, ее флуктуацию Та и время достижения данного состояния:

/ (?) *< /(?) >±,ЩГ).

Список литературы

1. Болотин В.В. Прогнозирование ресурса машин и конструкций. - М.: Машиностроение, 1984.

2. Либов Р. Введение в теорию кинетических уравне- 3. Израилев Ю.Л. Основы теории живучести турбин.

ний. - М.: Мир, 1974. Рекомендации и опыт реализации. - М.: Мин-во энергетики и

электрификации СССР, 1991.

Семенов Владимир Константинович,

Ивановский государственный энергетический университет,

доктор технических наук, профессор кафедры атомных электростанций,

телефон (4932) 385-778.

Беляков Андрей Александрович,

Ивановский государственный энергетический университет, аспирант кафедры атомных электростанций.

Щебнев Владимир Сергеевич,

Ивановский государственный энергетический университет,

кандидат технических наук, заведующий кафедрой атомных электростанций,

телефон (4932) 385-778.