В течение достаточно длительного периода времени авторами проводились исследования поведения различных оболочечных конструкций, испытывающих воздействие силовых и тепловых нагрузок, а также высокотемпературной водородной коррозии [2, 6, 7]. При этом коррозионное воздействие водорода на стали в настоящее время изучено достаточно хорошо, и на эту тему имеется множество публикаций [1, 4, 5]. В этих работах показано, что влияние водорода существенным образом сказывается на изменении вида диаграммы деформирования углеродистых сталей, при этом снижаются прочностные и пластические свойства сталей, и происходит их охрупчивание. Особенно сильно ухудшаются длительные прочностные характеристики сталей.
Целью настоящей работы является обобщение существующих подходов к учету влияния высокотемпературной водородной коррозии и практическая реализация этих подходов в составе разработанной авторами методики решения комплексной задачи по оценке прочности и долговечности стальных оболочек вращения с учетом необратимых деформаций и повреждаемости материалов [2, 6, 7].
Обобщенная модель химического взаимодействия
Предполагается, что влияние водорода, контактирующего с материалом конструкций, проявляется по истечении определенного времени, называемого инкубационным периодом tинк. После завершения инкубационного периода начинается интенсивное обезуглероживание материала конструкции, сопровождающееся ухудшением его механических свойств. При этом механические свойства материалов могут существенно изменяться. Этот период называется периодом активных химических превращений tкр. По завершении указанного периода дальнейших изменений механических свойств сталей под воздействием водорода не происходит.
Влияние водорода на свойства материала в данной модели представляется в виде дифференциального уравнения для параметра химического взаимодействия водорода с материалом конструкции µ, изменяющегося от m0 до mКР, принадлежащего интервалу (0; 1) и характеризующего степень поражения материала вследствие водородной коррозии [5]:
(1)
с условиями:
t = 0, m = m0(p, T, wС);
t = tинк, m = mП(p, T, wС); (2)
t = tкр, m = mКР(p, T, wС).
Здесь k(p, T, wС) – коэффициент, учитывающий влияние давления p, температуры T и степени поврежденности wС на кинетику химических превращений, приводящих к деградации физико-механических свойств материала.
Решение уравнения (1) имеет вид [5] (кривая 1 на рис. 1):
, (3)
g0 = (1 – m0) / m0; (4)
; (5)
gП = (1 – mП) / mП, (6)
где m, u, B, z – константы, экспериментально определяемые для различных марок стали. Значения m0, mП и mКР определяются по методике, изложенной в работах [3, 5], и связаны с анализом кинетики выделения метана при взаимодействии водорода с углеродом, содержащемся в материале стенки.
Рис. 1.
При этом учете влияния уровня напряжений на длительность инкубационного периода и скорость высокотемпературной коррозии вводится зависимость коэффициента m от s в виде:
. (7)
Здесь: m0, G, g – константы материала; sЭ – эквивалентное напряжение, si – интенсивность нормальных напряжений, sср – среднее напряжение.
Распределение давления по толщине оболочки вычисляется [5]:
для плоской стенки:
, (8)
для оболочки:
, (9)
где: z – плоская координата; r – сферическая координата; h – толщина оболочки; rВ – внутренний радиус оболочки; rН – наружный радиус оболочки; PВ – внутреннее давление; PН – наружное давление.
Упрощенная модель химического взаимодействия
Так как t < tинк изменений механических свойств материалов еще не происходит, а при t > tкр уже не происходит [5], то логистическую кривую Ферхюльста можно аппроксимировать трехзвенной линией (на рис. 1 линия 2).
В этом случае влияние водорода на свойства материала будет описываться в виде кусочно-линейного закона для параметра химического взаимодействия водорода µ, изменяющегося от 0 до 1 и характеризующего степень поражения материала вследствие водородной коррозии:
(10)
Длительность инкубационного периода, как показано в работе [4], определяется с помощью следующего соотношения:
tинк= m p–u exp(B/T), (11)
где p – парциальное давление водорода, МПа; T – температура; m, u, B – некоторые константы, экспериментально определяемые для различных марок стали.
Длительность периода активных химических превращений определяем следующим образом [3, 5]:
tкр = λ tинк, (12)
где: λ – параметр, являющийся мерой скорости химического взаимодействия водорода с конструкцией, показывающий, во сколько раз время до завершения периода изменений механических свойств материала отличается от продолжительности инкубационного периода, т.е. λ = tкр/tинк. Для стали 20 этот параметр принимает значения в пределах от 5 до 10 [3, 5].
В том случае, если известно, что процесс обезуглероживания конструкции значительно короче времени до начала разрушения, можно использовать еще более простые модели воздействия водородосодержащей среды. В частности, замену логистической кривой, описывающей плавное изменение механических характеристик, на ступенчатое их изменение от начального до конченого состояния, по прошествии либо инкубационного периода, либо периода активных химических превращений.
Рис. 2.
На рис. 2 схематически показан примерный характер изменения временного сопротивления стали 20 вследствие воздействия водорода. Здесь sв и sв* – соответственно значения временного сопротивления материала в исходном и обезуглероженном состоянии. Другие механические характеристики материала (E, sт, d, y) изменяются в процессе обезуглероживания аналогичным образом [3, 5].
При использовании варианта 1 считается, что механические характеристики ступенчато изменяются по истечении периода времени, соответствующего tкр, а в варианте 2 это изменение происходит сразу же после завершения инкубационного периода, то есть при tкр = tинк. При этом вариант 1 дает завышенную оценку времени обезуглероживания, а вариант 2 – заниженную. Кроме того, в варианте 2 не требуется нахождения параметра λ .
Учет влияния напряжений на скорость коррозии и распределение давлений осуществляется так же, как и в обобщенной модели.
Модель оценки кинетики движения фронта обезуглероживания
Предполагается, что воздействие водорода на материал начинается с поверхности оболочки, непосредственно контактирующей с водородосодержащей средой, и проявляется в обезуглероживании материала. По мере проникновения водорода в материал этот процесс распространяются в глубь конструкции с образованием области обезуглероженного материала. Поверхность, разграничивающую области материалов в исходном и обезуглероженном состоянии, будем называть фронтом обезуглероживания.
Кинетика перемещения фронта обезуглероживания будет определяться следующими выражениями [5]:
для пластины:
, (13)
для оболочки:
, 
, (14)
где z – глубина обезуглероживания, отсчитываемая от поверхности контакта оболочки с водородом; tфронта – время, для которого определяется глубина обезуглероживания.
Эта модель построена на основе упрощенной модели химического взаимодействия, целью ее разработки была возможность использования при аналитических расчетах, однако она может использоваться только для некоторых несложных моделей пластин и оболочек и неспособна учитывать нестационарный характер изменения тепловых, силовых нагрузок и давления водорода, а также влияние уровня напряжений на скорость движения фронта обезуглероживания.
В действительности моделей воздействия водородосодержащей среды значительно больше, наиболее распространенные из них подробно рассмотрены в работах [3, 5]. Здесь же приведены только те модели, которые могут использоваться в составе разработанного авторами программного комплекса по расчету напряженно-деформированного состояния оболочечных конструкций [2, 6, 7].
Технически учет влияния водорода на прочностные свойства конструкции реализован следующим образом. Для каждого из материалов конструкции введены механические характеристики двух его состояний – исходного и обезуглероженного.
Процесс нагружения оболочки будем разбивать на ряд малых по времени этапов Dt, в течение которых эти параметры можно принять постоянными, причем длительность этапов также может изменяться в зависимости от скорости изменения давления и температуры. Оболочку разбиваем на малые элементы, напряженно-деформированное состояние, а также степени поврежденности и обезуглероживания которых можно принять однородными. На каждом этапе нагружения и для каждого элемента оболочки будет вычисляться приращение Dm(Dt) по формулам (3) – (6), если используется обобщенная модель химического взаимодействия и (10) если упрощенная. Если процесс обезуглероживания отсутствует (инкубационный период не завершен), то для расчета параметров напряженно- деформированного состояния берутся исходные механические характеристики материала. Если в материале уже происходит процесс обезуглероживания, но он еще полностью не завершен, то значения механических характеристик материала на этой стадии расчета вычисляются путем линейной интерполяции в зависимости от текущего значения m(t) между исходным и обезуглероженным состояниями. Если процесс обезуглероживания завершен, то для расчета берутся характеристики обезуглероженного материала.
В случае использования упрощенной модели химического взаимодействия со ступенчатым переходом к обезуглероженному состоянию для каждого элемента оболочки по формуле (12) вычисляется значение tкр, при наступлении которого и осуществляется переход к обезуглероженному состоянию.
С использованием этого подхода авторами разработана методика, позволяющая исследовать историю изменения осесимметричного упругопластического напряженно-деформированного состояния однослойных и многослойных оболочек вращения с учетом повреждаемости материалов при ползучести и высокотемпературной водородной коррозии, а также оценивать их несущую способность и долговечность. Использование этой методики позволит решать целый ряд важных технических задач, таких как расчет на прочность, оценка несущей способности, долговечности и остаточного ресурса, а также задачи проектирования оптимальных конструкций, применяемых в энергетическом и химическом машиностроении.
Рецензенты:
Богданов Евгений Павлович, профессор, доктор технических наук, профессор кафедры «Информационные системы и технологии» ФГОУ ВПО Волгоградская государственная сельскохозяйственная академия, г. Волгоград.
Овчинников Игорь Георгиевич, доктор технических наук, профессор, ФГОУ ВПО «Саратовский государственный технический университет имени Гагарина Ю. А.», г. Саратов.
Библиографическая ссылка
Белов А.В., Неумоина Н.Г., Поливанов А.А., Попов А.Г. К ВОПРОСУ ОБ ИСПОЛЬЗОВАНИИ РАЗЛИЧНЫХ МОДЕЛЕЙ УЧЕТА ВОЗДЕЙСТВИЯ ВОДОРОДОСОДЕРЖАЩЕЙ СРЕДЫ ПРИ РАСЧЕТЕ НЕСУЩЕЙ СПОСОБНОСТИ СТАЛЬНЫХ ОБОЛОЧЕЧНЫХ КОНСТРУКЦИЙ // Современные проблемы науки и образования. – 2012. – № 6. ;URL: https://science-education.ru/ru/article/view?id=7854 (дата обращения: 19.04.2024).