автореферат диссертации по энергетике, 05.14.03, диссертация на тему:Расчетное моделирование процессов накопления повреждений корпуса коллектора парогенератора ПГВ-1000 в вероятностной постановке

кандидата технических наук
Алексеев, Петр Викторович
город
Москва
год
2013
специальность ВАК РФ
05.14.03
цена
450 рублей
Диссертация по энергетике на тему «Расчетное моделирование процессов накопления повреждений корпуса коллектора парогенератора ПГВ-1000 в вероятностной постановке»

Автореферат диссертации по теме "Расчетное моделирование процессов накопления повреждений корпуса коллектора парогенератора ПГВ-1000 в вероятностной постановке"

Национальный исследовательский центр «Курчатовский институт»

Институт реакторных материалов и технологий (КЦЯТ)

На правах рукописи

АЛЕКСЕЕВ Петр Викторович

РАСЧЕТНОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССОВ НАКОПЛЕНИЯ ПОВРЕЖДЕНИЙ КОРПУСА КОЛЛЕКТОРА ПАРОГЕНЕРАТОРА ПГВ-1000 В ВЕРОЯТНОСТНОЙ ПОСТАНОВКЕ

Специальность 05.14.03 - ядерные энергетические установки, включая проектирование, эксплуатацию и вывод из эксплуатации

Автореферат диссертации на соискание ученой степени кандидата технических наук

Москва — 2013

9 ЯН,В 2014

005544970

Работа выполнена в Институте реакторных материалов и технологий (КЦЯТ) Национального исследовательского центра "Курчатовский институт"

Научный руководитель: доктор технических наук, профессор

Тутнов Александр Александрович

зам. директора ИРМТ КЦЯТ НИЦ "Курчатовский институт"

Официальные оппоненты: доктор технических наук, профессор

Гетман Александр Федорович

ВНИИАЭС, главный научный сотрудник

кандидат технических наук Ткачев Владимир Валентинович

отдел канальных реакторов ИАС КЦЯТ НИЦ "Курчатовский институт", начальник лаборатории

Ведущая организация: ОАО "ОКБМ Африкантов"

г. Нижний Новгород

Защита состоится_в_часов_минут на заседании

диссертационного совета_в НИЦ "Курчатовский институт"

по адресу: 123182, г. Москва, пл. Академика Курчатова, д. 1

С диссертацией можно ознакомиться в научно-технической библиотеке НИЦ "Курчатовский институт"

Автореферат разослан _

Ученый секретарь диссертационного совета

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ Актуальность темы исследования

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

Кроме того, действующие руководящие документы, используемые при расчетах на прочность оборудования и трубопроводов ЯЭУ, устанавливают предельно допустимую вероятность возникновения аварии, некомпенсируемой

штатной системой безопасности, равную 10 (1/год) на установку. Указанное требование нормативных документов, наряду с необходимостью снятия излишнего консерватизма расчетных обоснований прочности конструкций, в том числе и коллекторов парогенераторов (ПГ) типа ВВЭР-1000, определяет актуальность темы диссертации, направленной на создание расчетного инструмента для доказательства маловероятного (гипотетического) характера отказов, недопустимых с точки зрения обеспечения безопасности ЯЭУ.

Степень разработанности темы исследования

Проблема моделирования повреждаемости элементов конструкции ЯЭУ изучается давно. За это время были разработаны методики и расчетные коды, позволяющие моделировать процессы повреждения элементов конструкций с учетом различных механизмов разрушения, а также оценить их надежность и долговечность. Основные существующие методики подробно рассмотрены в работе. Однако стоит заметить, что большинство из них являются "незавершенными" с точки зрения применимости на практике для обоснования надежности элементов конструкции реакторных установок. Наибольшую практическую ценность представляют собой вероятностные методики моделирования, созданные на основе метода пересечения статистических величин. К таким методикам в частности относится комплекс алгоритмов, реализованный в расчетных кодах МАВР-1.1, МАВР-2.1. Однако, на сегодняшний день, указанные методики устарели в связи с большим прогрессом за последние 20 лет в совершенствовании механики разрушения и изменением нормативной базы. Подводя итог, можно сказать, что, несмотря на множество существующих, на сегодняшний день, методов моделирования повреждаемости и оценки надежности конструкций, завершенного инструмента для обоснования безопасной эксплуатации реакторных установок, отвечающего современным достижениями механики разрушения и требованиям актуальных руководящих документов, не существует.

Цели и задачи работы

Целью проводимых в диссертации исследований является решение следующих задач:

1. Разработка комплекса методик для моделирования повреждаемости корпусов коллекторов ПГ.

Для достижения этой цели потребовалось:

• разработать методику моделирования зарождения новых дефектов в процессе эксплуатации;

• разработать методики вычисления критических размеров трещин в соответствии с требованиями актуальных руководящих документов;

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

2. Разработка и внедрение современного программного средства для оценки вероятности возникновения течи корпуса коллектора ПГ. Для достижения этой цели потребовалось:

• разработать программный комплекс (ПК) на основании созданной методики моделирования повреждаемости;

• провести анализ влияния неопределенности отдельных блоков программы на результаты расчета вероятности возникновения течи корпуса коллектора ПГ;

• провести верификационные исследования программного средства;

• провести апробацию ПК, путем применения его для практического обоснования надежности элементов РУ.

Научная новизна работы

Заключается в:

• разработке методики моделирования процессов зарождения микроповреждений в стенке корпуса коллектора на основе трех возможных механизмов повреждаемости: разрушения вследствие накопления неупругих

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

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

• реализации методики определения критических размеров дефектов при вязком разрушении через локальную и общую пластическую нестабильность;

• разработке программного блока, реализующего метод «теней» при моделировании циклической нагрузки, а также блоки и методики оценки значений коэффициента интенсивности напряжений за пределами диапазона применимости инженерных формул, приведенных в отраслевых нормативных документах;

• создании ПК МАВР-4.1 для расчета вероятности возникновения течи коллекторов ПГ, на основании новых методик и блоков расчета, а также блоков, заимствованных из ПК МАВР-1.1 и МАВР-2.2;

• создании версии МАВР-4.1, реализующей механизм параллельных вычислений на многопроцессорном кластере;

• проведении исследований влияния неопределенности исходных данных на результаты расчета вероятности возникновения течи коллектора, с получением соответствующих зависимостей;

• оценке точности расчетов с помощью программного средства, с учетом погрешностей вычислительных алгоритмов и неопределенности в задаваемых исходных данных;

• проведении верификационных исследований, показывающих соответствие между эксплуатационными данными о выходе из строя коллекторов и глушении

теплообменных трубок ПГ ВВЭР-1000 и расчетными прогнозами, полученными с использованием МАВР-4.1;

• проведении обоснования гипотетического характера события, заключающегося в крупномасштабном разрушении коллектора ПГ ПГВ-100МКП проекта ВВЭР-ТОИ, возникновении некомпенсируемых системой безопасности течей, в соответствии с современными представлениями механики разрушения, требованиями руководящих документов и учетом погрешности вычислений программного средства.

Практическая и теоретическая значимость работы

Результаты данной работы имеют практическую значимость при прогнозировании и обоснованиях надежности и долговечности ПГ. Созданный ПК МАВР-4.1 позволяет оценить вероятность возникновения течи коллектора ПГ на разных этапах его эксплуатации. Эти результаты дают возможность обосновать надежность коллектора ПГ в соответствии с требованием нормативных руководящих документов обеспечить вероятность безотказной работы ЯЭУ не более 10~7 на реактор в год. Оценка неопределенности программного средства, полученная в рамках настоящей работы, позволяет проводить данные обоснования с учетом погрешности результатов расчетов, обусловленной неопределенностью исходных данных и погрешностью компьютерных вычислений.

Результаты исследований по влиянию неопределенности исходных данных на вероятность возникновения течи коллектора ПГ могут быть использованы при проектировании или продлении срока службы ЯЭУ в качестве основания для выбора приоритетных направлений в работах по улучшению прочностных или эксплуатационных характеристик ПГ.

Методология и методы исследований

Все исследования, проведенные в рамках данной работы, проводились с использованием ПК МАВР-4.1. В основе МАВР-4.1 лежат законы и

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

Положения, выносимые на защиту

• Методика моделирования повреждаемости коллектора в вероятностной постановке и ее реализация в программном комплексе МАВР-4.1;

• результаты исследования влияния неопределенности исходных данных на результаты расчетов вероятности возникновения течи;

• результаты верификационных исследований;

• результаты применения МАВР-4.1 для обоснования прочности коллектора ПГ в проекте ВВЭР-ТОИ.

Степень достоверности и апробация результатов работы

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

Результаты исследований и разработки отдельных методик докладывались на различных конференциях и семинарах, а именно: на 13-й и 14-й конференции молодых специалистов по ядерным энергетическим установкам, проведенных ОКБ "Гидропресс" в 2011 и 2012 г.г. в г. Подольске; на 7-м межотраслевом семинаре "Прочность и надежность оборудования" в 2011 г. в г. Звенигороде; на 7-й Российской конференции "Методы и программное обеспечение расчетов на прочность" в 2012 г. в г. Геленджике.

Кроме того, по теме работы было опубликовано 2 статьи в журнале "Атомная энергия", а также выпущен ряд научно-технических отчетов в НИЦ "Курчатовский институт".

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

Работа содержит описание вероятностной методики моделирования накопления повреждений, разработанной для оценки вероятности возникновения течей корпуса коллектора ПГ, а также, описание программного кода МАВР-4.1, созданного на основе данной методики, результаты исследования влияния неопределенности исходных данных на вероятность возникновения течи, верификационные исследования и краткое описание результатов практического применения МАВР-4.1.

Коллектор ПГ типа ПГВ-1000 представляет собой цилиндрический сосуд давления, закрепленный с двух концов в корпусе ПГ. На каждый ПГ приходится по два коллектора: "холодный" и "горячий". Конструктивные особенности коллектора связаны с наличием перфорации под теплообменные трубки, а также "клина" - зоны сопряжения двух цилиндрических поверхностей: корпуса коллектора и корпуса ПГ (рисунок 1).

(6)

"а" - сопряжение корпусов коллектора и парогенератора; "б" - коллектор и зона его перфорации; "в" - схема коллектора и "клина". Рисунок 1. Схема коллектора и "клина".

Проведенные анализы повреждений коллекторов ПГ, выявленных в 1978-1991 годы показали, что развитию трещин способствует, главным образом, статическое и циклическое нагружение конструкции, а также коррозионные процессы, возникающие в зазоре между теплообменной трубкой и корпусом коллектора при недовальцовке трубок. Большинство повреждений было расположено в зоне "клина", что объясняется повышенной концентрацией напряжений в этой зоне. Разработанная методика учитывает результаты данных исследований и эксплуатационно-конструкционные особенности коллектора ПГ.

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

См и

(С) = | (х, фхкр (х, С)<кс1С (1)

с, о

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

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

на ряд равных временных интервалов, и последующие стадии моделирования выполняются последовательно на каждом отрезке времени. Вторым этапом производится моделирование зарождения новых дефектов в процессе эксплуатации. Возможны следующие механизмы накопления повреждений:

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

• термически активированное коррозионное растрескивание,

• коррозионно-усталостное повреждение.

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

(2)

Функция повреждения а>>7, вычисляется путем нелинейного суммирования функций повреждения от воздействия каждого из трех механизмов зарождения трещины: сое, сок и сос - повреждения за счет неупругих деформаций, термически активированнные коррозионные повреждения и коррозионно-усталостные повреждения соответственно:

«V = (<' + + К3 У4 (3)

Коэффициенты т1...т4 подбираются при математической обработке экспериментальных данных, в случае их равенства единице формула (3) отражает линейный закон суммирования повреждений.

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

9

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

со.

= I ' (4)

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

) ^Т (5)

0гк(Сс„Со2, ■ •)

где 1долговечность перемычки (участка перфорированной зоны коллектора,

заключенного между соседними отверстиями), зависящая от химического состава среды и прочих коэффициентов, характеризующих условия нагружения.

Моделирование циклической повреждаемости выполняется в соответствии с Нормами расчета на прочность оборудования и трубопроводов атомных энергетических установок. Выражение для функции циклической повреждаемости записывается в соответствии с принципом линейного суммирования повреждений:

к N.

=2 ж:« ' = 1 01

где Л^ - число циклов /-го типа с одинаковой амплитудой напряжений, асимметрией цикла и температурой; к - общее число различных типов циклов;

Л^ - число циклов г-го типа, необходимое для разрушения материала.

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

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

(7)

где Сл- характеристика материала, - коэффициент интенсивности

напряжений.

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

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

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

2 Вр

КМ-Кщ п

уК'М-К^

аь< 1 (8)

где: Вр - протяженность фронта расчетного дефекта; К,{ь) — распределение {К,) по фронту Ь расчетного (полуэлиптического) дефекта; К)С(Ь) — распределение вязкости разрушения по фронту трещины Ь; Кт^п— константа материала;

где / - параметр механики разрушения (./-интеграл), определяемый численным методом, при решении упруго-пластической задачи, или с использованием нормативных инженерных формул.

Стоит заметить, что одной из особенностей настоящей методики является учет влияния эффектов коротких трещин и двуосного нагружения на вязкость разрушения, в соответствии с методиками расчета на сопротивление хрупкому разрушению корпусов реакторов АЭС с ВВЭР-1000 при эксплуатации, а также, учет влияния перфорации и геометрии поверхности на вязкость разрушения. Было показано, что вязкость разрушения не является табличной характеристикой материала. Существует ряд факторов, существенно влияющих на данную величину. Такими факторами являются эффекты коротких трещин и двуосного нагружения, а также наличие перфорации в стенке корпуса коллектора. Учет

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

Методика предусматривает моделирование по трем механизмам разрушения: хрупкому, хрупко-вязкому и вязкому. Выбор одного из трех механизмов определяет критерий критического подрастания дефектов.

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

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

Л = а(Л)уМ4>, (9)

А0=2пас2/Е' (10)

где А - площадь раскрытия трещины, А^ - упругая составляющая площади раскрытия трещины для бесконечной пластины с центральной трещиной длиной 2с под действием номинального напряжения о, коэффициент а(Я) - поправка на кривизну оболочки, Я - параметр оболочки, X5) - поправка на размер пластической зоны, Е равно модулю упругости Е при плоско-напряженном

состоянии (тонкостенные сосуды) и Е =Е(} — у2) для плоской деформации, V -коэффициент Пуассона.

Описанная методика была реализована в виде ПК МАВР-4.1. Расчеты вероятности возникновения течи можно условно разделить на три этапа:

• расчет поля остаточных напряжений, возникающих при запрессовке теплообменных труб, с учетом технологии запрессовки;

• расчет нестационарных полей температуры и напряжений в корпусе коллектора при проектных режимах эксплуатации, включая гидроиспытания, аварийные ситуации, проектное землетрясение и максимальное расчетное землетрясение;

• расчет вероятности разрушения коллектора в процессе эксплуатации путем раздельного вычисления вероятностей возникновения течей разного размера и вероятности крупномасштабного разрушения (разрыва).

Первые два этапа расчетов выполняются с помощью известных конечноэлементных кодов общего назначения (например, АКБУБ (в ОКБ "Гидропресс"), 1^011-1 (в НИЦ "Курчатовский институт") и др.). На основании этих вычислений формируются файлы, содержащие поля температуры и напряжений (интенсивности или приведенных, а также компонент тензора напряжений в координатах, привязанных к плоскости трещины), которые являются частью исходных данных для МАВР-4.1.

Следует заметить, что вероятностные вычисления связаны с большим количеством исходных и промежуточных данных, что рождает потребность в существенных вычислительных ресурсах. Для сокращения времени машинных вычислений была создана параллельная версия программы МАВР-4.1. Данная версия адаптирована для работы на многопроцессорных вычислительных кластерах, позволяющих рассчитывать объемные массивы одновременно большим числом процессоров. В настоящее время, параллельная версия МАВР-4.1 установлена на многофункциональном вычислительном комплексе НИЦ "Курчатовский институт". Благодаря существенному снижению времени вычислений, впервые появилась возможность проведения полноценных 14

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

В рамках работы были проведены исследования, направленные на получение сведений о влиянии неопределенности исходных данных на вероятность возникновения течи корпуса коллектора ПГ, а также, оценена суммарная неопределенность программного средства. Исследования преследуют две цели: получение представлений о влиянии наиболее значимых исходных характеристик материалов, конструкционных и эксплуатационных факторов на результат расчета вероятности возникновения течей, и оценка неопределенности результатов расчетов МАВР-4.1, обусловленной погрешностью исходных данных. Работа проводилась на базе данных проекта АЭС-2006 (ВВЭР-ТОИ), для течи с условным диаметром ДУ-43. В качестве примера, в диссертационной работе приведены исследования для 4-х наиболее значительных параметров исходных данных, характеризующих работу отдельных блоков программы: начальная дефектность (стандартное отклонение начальной дефектности); параметр, характеризующий докритический подрост трещин (коэффициент ускорения в уравнении роста трещины в водной среде (уравнении Париса)); критический коэффициент интенсивности напряжений, определяющий, главным образом, критическую фазу развития дефектов (матожидание и стандартное отклонение критического КИН); поле напряжений, влияющее на результаты моделирования на каждом этапе. Также было определено влияние нарушений водно-химического режима (ВХР) на вероятность возникновения течи коллектора. В рамках исследования, в рассматриваемый параметр исходных данных вносилось отклонение его от номинального значения, производился расчет вероятности возникновения течи и сопоставление полученного результата с исходным. В качестве характеристики степени влияния исходных данных на результат был применен коэффициент влияния, показывающий, на сколько процентов изменится результат, в случае отклонения рассматриваемого параметра на один процент. Коэффициент влияния вычисляется из отношения:

где Аисх - процент изменения исходных данных относительно "номинального" варианта; Арез - процент изменения результата расчета относительно "номинального" варианта, вычисляемый по формуле:

где АР - величина, на которую изменился результат; Рисх - расчетная вероятность для исходного варианта.

Данный коэффициент был рассчитан для каждого параметра исходных данных для трех возможных ориентаций трещин в корпусе коллектора на протяжении проектного срока службы РУ 60 лет. Краткие сведения о результатах приведены в таблице 1. Минимальные и максимальные значения коэффициентов влияния позволяют не только оценить степень влияния той или иной исходной характеристики на вероятность возникновения течи, но и сделать выводы о слабом изменении коэффициентов влияния на протяжении всего периода эксплуатации.

^■рез = р "100%

|ДР|

(12)

реэ

Таблица 1. Диапазон изменения коэффициентов влияния.

Рассматриваемый параметр Направление трещины Минимальный коэффициент влияния Максимальный коэффициент влияния

Стандартное отклонение начальной дефектности прямая 4,23 4,29

кольцевая 3 4,2

косая 2,5 4,2

Коэффициент в уравнении Париса прямая 1,7 7,1

кольцевая 0,2 3,6

косая 0,1 16

Стандартное отклонение критического КИН прямая 25 70

кольцевая 0,01 0,64

косая 0,01 0,17

Матожидание критического КИН прямая 15,4 17,1

кольцевая 7,5 11,6

косая 0,01 1,3

Напряжения прямая 12 67

кольцевая 0,1 15,9

косая 1,6 12,7

Исследование влияния ВХР на вероятность возникновения течи показали, что при эксплуатации с концентрацией хлора в теплоносителе, равной 300 мкг/кг, на протяжении всего проектного срока службы, вероятность возникновения течи может возрасти на 2,5 порядка. Однако такие значения наблюдаются лишь на начальных этапах эксплуатации старых моделей коллекторов, где вальцовка теплообменных трубок проводилась взрывом и верхняя часть коллектора была заневолена в корпусе ПГ. Подобное нарушение ВХР при эксплуатации

современных ПГ значительно менее существенно влияет на надежность коллектора и достигает отклонения в 2-3 раза на поздних этапах эксплуатации.

Также была вычислена общая неопределенность, вносимая всеми рассмотренными исходными данными, и представлена в виде зависимостей позволяющих оценить изменение неопределенности в определении вероятности возникновения течи в зависимости от срока службы РУ, для трех ориентаций дефектов (рисунок 2).

неопределенность

м-.........=.........г?....... 5.....- - с\.......:........31

2,5 .....................................—"

1Д ЩЙДИ .............:........-.........—...... И ----------------------------------------------

0.5 —----- -4-----------

о_______:_:_|___________________;__________.________________г__________........._________£лет,

О 10 20 30 40 50 60

"/"- для вертикальных трещин; "2"- для кольцевых трещин; "3" - для косых

трещин.

Рисунок 2. Зависимость от времени неопределенности конечного расчетного значения вероятности возникновения течи для трех ориентаций трещины в стенке коллектора

Максимальная неопределенность расчетов по программному средству МАВР-4.1, за счет неопределенности исходных данных составляет 3,87, а погрешность, обусловленная особенностями машинных вычислений, составляет 0,13 (в долях). Суммарная максимальная погрешность МАВР-4.1 составляет 3,88 (в долях).

Диссертационная работа также содержит верификационные исследования МАВР-4.1. Верификация такого рода программ связана с трудностями,

\

ц

2 --1

/

т г*

\3\ г.т >т

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

Проводилось сравнение расчетной наработки на отказ ПГ, полученной с использованием МАВР-4.1, и эксплуатационных данных о повреждаемости коллекторов первого контура. В качестве эксплуатационных были взяты данные о выходе из строя по причине возникновения течи теплоносителя первого контура ПГ ПГВ-1000 в период 1978 - 1991 г.г. Расчетный промежуток процента отказов на разных стадиях эксплуатации, а также эксплуатационные данные о выходе из строя представлены в таблице 2. В последней графе таблицы 2 приведена доля отказов «холодных» коллекторов от всего числа «холодных» коллекторов, с трубчаткой, завальцованной взрывом и находившихся в эксплуатации в рассматриваемый период.

Таблица 2. Сопоставление расчетных и эксплуатационных данных о выходе коллекторов ПГ из строя.

Время Расчетное Минимальная Макси- Эксплуата- Доля

работы, матожидание расчетная мальная ционная доля отказов

лет отказа на доля отказов, расчетная отказов, холодных

интервале (в%) доля (в % от 24 коллекто-

времени, отказов, «холодных» ров (в % от

(в%) (В%) коллекторов) их общего кол-ва)

0-2 12,52 9,29 61,1 41,69 15,63

2-4 9,04 7,45 44,12 33,33 12,50

4-6 6,32 5,19 30,84 8,33 3,12

4-8 5,18 4,16 25,28 16,67 6,25

8-оо 38,32 21,56 100 - 62,5

Из таблицы 2 видно, что эксплуатационные сведения о выходе из строя коллекторов по причине течеобразования попали в расчетный промежуток. Доля отказов от общего числа «холодных» коллекторов, находившихся в эксплуатации в рассматриваемый период, удовлетворительно коррелирует с расчетным матожиданием отказов (второй столбец). Это подтверждает достоверность расчетов по разработанной методике.

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

Также было проведено сопоставление расчетных и экспериментальных данных о выходе из строя теплообменных трубок ПГ в различные годы. Теплообменная трубка также является сосудом давления, и такая модель также может быть рассчитана при помощи МАВР-4.1, как частный случай, где диаметр

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

Сплошная линия — расчетный диапазон глушения трубок. Пунктирная линия -расширение расчетного диапазона за счет погрешности вычислений. Точками

показана эксплуатационная доля глушения в различные годы. Рисунок 3. Сопоставление расчетных и экспериментальных данных о глушении теплообменных труб в различные годы.

Видно, что все значения отказов в различные годы укладываются в расчетный диапазон отказов с учетом погрешности программного средства. Это также подтверждает достоверность результатов моделирования.

В главе 5 диссертации приводятся результаты вероятностного анализа разрушения коллектора ПГ в рамках проектов ВВЭР-ТОИ с исходными данными проекта АЭС-2006, проведенного с использованием ПК МАВР-4.1 Работа проводилась в соответствии с требованием Главного конструктора ВВЭР обосновать, что вероятность реализации аварии, связанной с течью коллектора ПГ из первого ко гп ура во второй с условным диаметром более 43 мм (Ду-43), не превышает величину 10~6 на реактор в год (либо указать технические и/или

глушения, %

1

годы

1974 1976 1978 1980 1982 1984 1986 1988 1990 2000 2002 2004

организационные мероприятия, выполнение которых обеспечит выполнение указанного выше требования).

Оценка вероятности возникновения течи проводилась для двух крайних вариантов распределений начальной дефектности по размерам (варианты № 1 и № 7), постулированных Главным конструктором на основании нормативных требований по входному и эксплуатационному неразрушающему контролю дефектов материала.

Расчеты вероятности разрушения при эксплуатации перфорированной зоны коллектора первого контура ПГ ПГВ-1000МКП показали следующие результаты (таблица 3):

Таблица 3. Результаты расчетов вероятности возникновения течи.

Вариант начальной дефектности Течь > ДУ-43 Течь > ДУ-100

Макс, вероятность течи коллектора за 60 лет №1 4-Ю"4 3-Ю"5

№7 6-Ю6 2 Ю"9

Макс, среднегодовая вероятность течи на блок РУ №1 5,4 10"5 4-Ю-6

№7 8-Ю"7 2,7-10"6

Вероятность крупномасштабного разрушения коллектора за 60 лет №1 2Л0'п

№7 2 10~13

Среднегодовая вероятность разрыва коллектора на блок РУ Для всех постулированных уровней дефектности 4-10"13

Таким образом, исследования показали, что вероятность разрушения корпуса коллектора РУ АЭС-2006, включая образование течи более Ду-43, на протяжении 22

периода эксплуатации 60 лет не превышает 10~б на реактор в год, если начальная дефектность корпуса в районе перфорированной зоны не хуже, чем задаваемая вариантом начальной дефектности № 7.

ЗАКЛЮЧЕНИЕ

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

• На основании разработанной методики создан программный комплекс МАВР-4.1, служащий для расчета вероятности возникновения течей корпуса коллектора на разных этапах эксплуатации реакторной установки.

• Проведено исследование и получены сведения о влиянии неопределенности основных исходных данных на вероятность возникновения течи теплоносителя первого контура.

• Оценена общая погрешность программного средства.

• Проведены верификационные исследования программного комплекса МАВР-4.1 и доказана достоверность результатов моделирования повреждаемости.

• Программный комплекс МАВР-4.1 использован при практических расчетах: проанализирована степень надежности корпуса коллектора парогенератора ПГВ-1000МКП и обоснованы условия неразрушающего контроля дефектов материала, при которых для проектов АЭС-2006 и ВВЭР-ТОИ выполняются требования руководящих документов по безопасности эксплуатации (в части непревышения предельно допустимого уровня вероятности возникновения течей и разрывов). Все обоснования проводились с учетом погрешности программного средства и неопределенности исходных данных.

СПИСОК РАБОТ, ОПУБЛИКОВАННЫХ АВТОРОМ ПО ТЕМЕ ДИССЕРТАЦИИ

Публикации в изданиях, рекомендованных ВАК России:

1. Алексеев, П.В. Расчет вероятности возникновения течи теплоносителя в коллекторе парогенератора ВВЭР по программе МАВР-4.1. /П.В.Алексеев, О.Д. Лоскутов, А.А.Тутнов // Атомная энергия. — 2013. - т.114, вып.1. — с. 17-21.

2. Алексеев, П.В. Вероятностный анализ прочности главного разъема реакторной установки ВВЭР-1000. /П.В.Алексеев, А-др.С. Киселев, Ал-й.С. Киселев и др.// Атомная энергия.-2013.- т.115, вып.З. - с. 147-151.

Материалы конференций:

1. Алексеев, П.В. Методические подходы к оценке влияния перфорации на трещиностойкость корпуса коллектора парогенератора. /П.В.Алексеев, О.Д. Лоскутов, А.А.Тутнов // Материалы 13-й конференции молодых специалистов по ядерным энергетическим установкам. - Подольск, 2011.

2. Алексеев, П.В. Методические подходы к оценке трещиностойкости корпуса коллектора парогенератора. /П.В.Алексеев, О.Д. Лоскутов, А.А.Тутнов// Материалы 7-го межотраслевого семинара "Прочность и надежность оборудования". — Звенигород, 2011 г.

3. Алексеев, П.В. Анализ влияния неопределенности исходных данных на вероятность возникновения течи теплоносителя парогенератора АЭС-2006. /П.В.Алексеев// Материалы 14-й конференции молодых специалистов по ядерным энергетическим установкам. — Подольск, 2012 г.

4. Лоскутов, О.Д. Вероятностный анализ коэффициентов запаса, используемых при расчете на прочность. /О.Д.Лоскутов, П.В.Алексеев, А.А.Тутнов// Материалы 7-й Российской конференции "Методы и программное обеспечение расчетов на прочность". - Геленджик, 2012 г.

5. Алексеев, П.В. Расчет вероятностей возникновения течей теплоносителя в коллекторе парогенератора ВВЭР (ПК МАВР-4.1). /П.В.Алексеев, О.Д. Лоскутов, А.А.Тутнов// Материалы 7-й Российской конференции "Методы и программное обеспечение расчетов на прочность". — Геленджик, 2012 г.

Подписано в печать 25.11.2013. Формат 60x90/16 Печать цифровая. Усл. печ. л. 1,5 Тираж 65. Заказ 95

Отпечатано в НИЦ «Курчатовский институт» 123182, Москва, пл. Академика Курчатова, д. 1

Текст работы Алексеев, Петр Викторович, диссертация по теме Ядерные энергетические установки, включая проектирование, эксплуатацию и вывод из эксплуатации

НАЦИОНАЛЬНЫЙ ИССЛЕДОВАТЕЛЬСКИЙ "КУРЧАТОВСКИЙ ИНСТИТУТ"

к

На правах рукописи

04201456007

АЛЕКСЕЕВ ПЕТР ВИКТОРОВИЧ

РАСЧЕТНОЕ МОДЕЛИРОВАНИЕ ПРОЦЕССОВ НАКОПЛЕНИЯ ПОВРЕЖДЕНИЙ КОРПУСА КОЛЛЕКТОРА ПАРОГЕНЕРАТОРА ПГВ-1000 В ВЕРОЯТНОСТНОЙ

ПОСТАНОВКЕ

Специальность 05.14.03 - Ядерные энергетические установки, включая проектирование, эксплуатацию и вывод из эксплуатации

Диссертация

на соискание ученой степени

кандидата технических наук

Научный руководитель:

д.т.н., проф., А. А. Тутнов

Консультант:

к.т.н., О. Д. Лоскутов

Москва, 2013

СОДЕРЖАНИЕ

1. Обзор методик моделирования повреждаемости конструкционных материалов ...4

1.1 Основные факторы, влияющие на процесс разрушения коллектора парогенератора ВВЭР-1000 и ВВЭР-1200 в процессе эксплуатации...................4

1.1.1. Общие сведения о причинах повреждения конструкции.................4

1.1.2. Механизмы зарождения и развития трещин в коллекторах парогенераторов ВВЭР-1000 и ВВЭР-1200..................................................5

1.2. Модели разрушения конструкционных материалов.......................................8

1.2.1. Общие подходы к моделированию процессов разрушения.............8

1.2.2. Полуэмпирические модели.................................................................11

1.2.3.Структурные модели...........................................................................32

2. Методика расчета вероятности разрушения коллектора парогенератора и описание программного средства МАВР-4.1..................................................................57

2.1. Общие положения............................................................................................57

2.2.Моделирование исходной дефектности материала с учетом пропуска дефектов при НК......................................................................................................59

2.3. Моделирование накопления повреждений....................................................64

2.4. Моделирование докритического роста дефектов..........................................79

2.5. Вычисление вероятности течи........................................................................89

2.6. Определение функции распределения критических размеров дефектов для разных механизмов разрушения.........................................................................108

2.6.1 хрупкое разрушение...........................................................................108

2.6.2. Хрупко-вязкое разрушение..............................................................114

2.6.3. Вязкое разрушение (пластический шарнир)...................................117

2.7. Расчет поправочного коэффициента для учета влияния перфорации на вязкость разрушения.............................................................................................120

2.8. Расчет эквивалентного диаметра течи..........................................................123

2.9. Структура программного средства МАВР-4.1............................................125

3. Анализ влияния неопределенности исходных данных на вероятность возникновения течи корпуса коллектора парогенератора..........................................130

3.1. Цели исследований........................................................................................130

3.2. Описание исходных данных.........................................................................130

3.3. Влияние начальной дефектности материала на вероятность возникновения течи........................................................................................................................136

3.4. Влияние параметров докритического подроста трещин на вероятность возникновения течи..............................................................................................140

3.5. Влияние критических свойств материала на вероятность возникновения течи.........................................................................................................................144

3.6. Влияние напряжений на вероятность возникновения течи.......................150

3.7. Определение неопределенности исходных данных....................................153

4. Верификация программного средства......................................................................155

4.1 .демонстрация правильной работы блоков программы................................156

4.2 сопоставление расчетных и экспериментальных данных о выходе из строя коллекторов ПГ по причине течеобразования....................................................158

4.3. Сопоставление расчетных и экспериментальных данных о глушении теплообменных трубок........................................................................................160

4.4. Неопределенность исходных данных...........................................................161

5. Применение МАВР-4.1 для практического обоснования надежности..................162

Заключение......................................................................................................................167

Литература.......................................................................................................................168

1. ОБЗОР МЕТОДИК МОДЕЛИРОВАНИЯ ПОВРЕЖДАЕМОСТИ

КОНСТРУКЦИОННЫХ МАТЕРИАЛОВ

1.1 основные факторы, влияющие на процесс разрушения коллектора

парогенератора ввэр-1000 и ввэр-1200 в процессе эксплуатации

1.1.1. Общие сведения о причинах повреждения конструкции

Если исключить из рассмотрения выходы из строя конструкций вследствие резких нерасчетных перегрузок, природных катаклизмов, неподдающихся контролю, грубых ошибок при проектировании или эксплуатации или неблагоприятного сочетания перечисленных факторов, то остальные случаи выхода конструкций и машин из строя можно отнести преимущественно к одной из двух больших групп: из-за предельных состояний, связанных с чрезмерным износом трущихся деталей и поверхностей, находящихся в контакте с рабочей или окружающей средой, или же из-за предельных состояний, наступивших в результате постепенного накопления в материале рассеянных повреждений приводящих к зарождению и развитию макроскопических трещин. Коллекторы парогенераторов и прочие сосуды давления можно отнести к классу элементов конструкций, работающих при высоких уровнях общей нагруженности, и степень механической изнашиваемости которых сравнительно мала. Поэтому первая группа причин повреждаемости в данной работе рассматриваться не будет, и все причины достижения предельных состояний будут отнесены ко второй группе. Зачастую накопление повреждений происходит еще на стадиях изготовления, монтажа или подготовки к эксплуатации. Кроме того, существуют работы, целью которых являлось изучение основных причин выхода их строя "холодных" коллекторов парогенераторов ПГВ-1000 [1]. Результаты данных исследований также указывают на деградацию прочностных свойств материала перфорированной зоны коллектора как на одну из основных причин возникновения течи теплоносителя и, как следствие, выхода из строя парогенератора. Деградация является следствием

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

Выделяют несколько основных факторов, вызывающих деградацию прочности оборудования и трубопроводов РЭУ:

• пластическая деформация;

• радиационное повреждение;

• циклическая усталость;

• коррозия

Из данного списка можно исключить радиационное повреждение, как фактор, не имеющий отношения к парогенератору. Несмотря на то, что элементы парогенератора РУ непосредственно контактируют с теплоносителем первого контура, наведенная активность является недостаточной, чтобы вызывать значительное изменение свойств материала. При расчете на прочность элементов парогенератора воздействием облучения также пренебрегают согласно действующим руководящим документам. Остальные перечисленные причины старения металла существенно влияют на надежность и долговечность парогенератора.

1.1.2. Механизмы зарождения и развития трещин в коллекторах парогенераторов ВВЭР-1000 и ВВЭР-1200

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

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

- математический разрез в теле из неповрежденного материала. Трещину считают заданной, а ее размер достаточно большим по сравнению с максимальным размером структуры материала - размером зерна, кристаллита, волокна и т. п. Такие трещины называют макроскопическими (в отличие от микроскопических трещин, размер которых имеет порядок характерного размера структуры материала или менее). Задача состоит в том, чтобы найти закономерности роста трещины при различных свойствах материала и различных процессах нагружения, а также установить условия, при которых этот рост устойчив, т. е. малые приращения нагрузок или малые изменения размеров трещин не приводят к ее интенсивному росту. В действительности процесс разрушения можно разделить на две стадии. Первая стадия - накопление рассеянных повреждений. Зарождение макротрещины и накопление микроповреждений как правило занимает значительную часть общего ресурса (по различным данным от 50 до 90 %). Основными причинами его, актуальными для парогенераторов, является термически активированное коррозионное растрескивание, накопление неупругих деформаций до предельного уровня, определяемого деформационной способностью материала при длительном нагружении и коррозионно-усталостное повреждение. Накопление повреждений зачастую начинается еще до начала эксплуатации конструкции. В частности, значительная часть микроповреждений в корпусе коллектора парогенератора ПГВ-1000 может образоваться на этапе вальцовки теплообменных трубок. Изначально процесс вальцовки производился с помощью взрыва, что, по мнению многих специалистов, явилось доминирующей причиной выхода из строя серии парогенераторов 111 Б-1000 в период 1978-1991 г. по причине повреждения корпуса коллектора. Впоследствии вальцовка взрывом была заменена на гидровальцовку, что положительно сказалось на замедлении процесса образовании рассеянных повреждений и, как следствие, надежности парогенератора в целом.

Вторая стадия - рост и развитие макротрещин. Если в детали или элементе не было начальных технологических трещин, то зарождение первой макроскопической

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

Подтверждением правомерности использования такой схемы повреждения служат многочисленные исследования. Так например, в 2006 г ОКБ «ГИДРОПРЕСС» было проведено исследование темплетов парогенератора №1, 5-го блока НВ АЭС. Фрактографическими исследованиями было показано, что образование дефекта и его развитие на глубину от 2 до 3 мм, происходило, преимущественно, по механизму коррозионного растрескивания под напряжением, вместе с тем имели место и признаки коррозионной усталости. Между начальной стадией развития дефекта и на глубину ~ 2 - 3 мм и стадией его последующего развития имеется ступенька, что свидетельствует о смене механизма разрушения. Последующее развитие дефект происходило преимущественно по механизму коррозионной усталости с наличием участков коррозионного растрескивания. Между первой и второй стадиями развития дефекта имел место значительный инкубационный период. Разворот трещины по мере ее продвижения

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

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

1.2. Модели разрушения конструкционных материалов 1.2.1. Общие подходы к моделированию процессов разрушения

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

взаимодействии с окружающей средой. Современное состояние механики материалов и конструкций (теории упругости и пластичности, строительной механики, механики разрушения и др.), а также прикладных методов расчета машин и конструкций позволяет с большой степенью достоверности предсказывать поведение механических систем, если известны свойства материалов и заданы внешние воздействия. В теории надежности механических систем свойства материалов и воздействий приняты случайными, поэтому поведение объекта также носит случайный характер. Нормативные требования и технические условия эксплуатации накладывают определенные ограничения на эти параметры. Ограничения могут быть сформулированы в виде условия нахождения некоторого случайного вектора, зависящего от времени характеризующего качества объекта, в заданной области, Отказам и предельным состояниям соответствуют выходы этого случайного вектора из области допустимых состояний. Таким образом, основная задача теории надежности — оценка вероятности безотказной работы на заданном отрезке времени - сведена к задаче о выбросах случайных процессов. Соединение методов механики материалов и конструкций с теорией случайных процессов составляет основу теории надежности механических систем [4,5].

Различают два подхода к оценке надежности объекта: полуэмпирический и структурный. В некоторых источниках он носит различные названия. Так, например, в [3] он именуется как системная теория надежности, формально-математический, формально-статистический или феноменологический. Оценка надежности здесь осуществляется на основании обобщения результатов наблюдений и экспериментов и не ставится целью объяснение или полное описание существа явлений на физическом уровне. В частности, модель макроскопической трещины -пример подхода, который не учитывает элементы структуры реальных материалов. Другими примерами служат способы описания процессов накопления повреждений при циклических и длительных нагрузках, основанные на введении мер повреждений. Эти меры не допускают прямой интерпретации на уровне структуры

материала. Более того, мера повреждений вообще не имеет четкого физического истолкования, кроме, может быть, двух ее предельных значений, отвечающих начальному (неповрежденному) состоянию и состоянию полного исчерпания ресурса. Наиболее широко используются полуэмпирический подход при оценке надежности на основании статистики отказов. В настоящее время это особо распространено в связи с высоким уровнем совершенства вычислительной техники, позволяющей быстро обрабатывать, преобразовывать и хранить большое число элементов, требующихся для точного описания отказов объекта. Для конструкций, где получение развернутой статистики затруднено по причине малосерийности и/или крайне редкими фактами отказа наиболее актуален структурный подход к оценке надежности. Модели, созданные на его основе также именуются физико-статистическими [3]. Структур