В данной работе изложена разработанная авторами методика решения комплексной задачи по оценке прочности и долговечности стальных оболочек вращения с учетом необратимых деформаций, повреждаемости материалов вследствие ползучести и высокотемпературной водородной коррозии [1, 2, 5, 7]. Эта задача решалась в термовязкоупругопластической постановке с применением соотношений теории неизотермических процессов упругопластического деформирования элементов твердого тела по траекториям малой кривизны, линеаризованным методом дополнительных деформаций. Учет влияния повреждаемости материалов вследствие ползучести на процесс ее деформирования осуществляется путем введения параметра повреждаемости в соотношение термовязкопластичности, для определения которого используется хорошо зарекомендовавшее себя кинетическое уравнение повреждаемости в форме, предложенной Ю.Н. Работновым.
Коррозионное воздействие водорода на стали в настоящее время изучено достаточно хорошо, на эту тему имеется множество публикаций, в которых показано, что влияние водорода существенным образом сказывается на изменении вида диаграммы деформирования углеродистых сталей, при этом снижаются прочностные и пластические свойства сталей и происходит их охрупчивание, ухудшаются длительные прочностные характеристики сталей. Вместе с тем возможность комплексного учета повреждаемости материалов вследствие температурной ползучести и водородной коррозии почти нигде не рассмотрена. Таким образом, основной и наиболее сложной задачей исследования является определение условий перехода механических характеристик сталей из исходного состояния в обезуглероженное. Для этого используется несколько подходов, подробно изложенных авторами в публикации [4]. Рассмотрим некоторые из них.
Обобщенная модель химического взаимодействия
Предполагается, что влияние водорода, контактирующего с материалом конструкций, проявляется по истечении определенного времени, называемого инкубационным периодом tинк. После завершения инкубационного периода начинается интенсивное обезуглероживание материала конструкции, сопровождающееся ухудшением его механических свойств. При этом механические свойства материалов могут существенно изменяться. Этот период называется периодом активных химических превращений tкр. По завершении указанного периода дальнейших изменений механических свойств сталей под воздействием водорода не происходит.
Влияние водорода на свойства материала в данной модели представляется в виде дифференциального уравнения для параметра химического взаимодействия водорода с материалом конструкции µ, изменяющегося от m0 до mКР, принадлежащего интервалу (0; 1) и характеризующего степень поражения материала вследствие водородной коррозии [4]:
(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) имеет вид логистической кривой Ферхюльста [4] (кривая 1 на рис. 1):
, (3)
g0 = (1 – m0) / m0; (4)
; (5)
gП = (1 – mП) / mП, (6)
где m, u, B, z – константы, экспериментально определяемые для различных марок стали. Значения m0, mП и mКР определяются по методике, изложенной в работе [4], и связаны с анализом кинетики выделения метана при взаимодействии водорода с углеродом, содержащемся в материале стенки.
Рис. 1
При этом для учета влияния уровня напряжений на длительность инкубационного периода и скорость высокотемпературной коррозии вводится зависимость коэффициента m от s в виде:
. (7)
Здесь: m0, G, g – константы материала; sЭ – эквивалентное напряжение, si – интенсивность нормальных напряжений, sср – среднее напряжение.
Распределение давления по толщине оболочки вычисляется [4]:
для плоской стенки:
, (8)
для оболочки:
, (9)
где: z – плоская координата; r - сферическая координата; h – толщина оболочки; rВ – внутренний радиус оболочки; rН – наружный радиус оболочки; PВ – внутреннее давление; PН – наружное давление.
Упрощенная модель химического взаимодействия
Установлено, что при t < tинк изменений механических свойств материалов еще не происходит, а при t > tкр уже не происходит [4]. В этом случае логистическую кривую Ферхюльста можно аппроксимировать трехзвенной линией (на рис. 1 линия 2), и влияние водорода на свойства материала будет описываться в виде кусочно-линейного закона для параметра химического взаимодействия водорода µ, изменяющегося от 0 до 1 и характеризующего степень поражения материала вследствие водородной коррозии:
(10)
Длительность инкубационного периода, как показано в работе [4], определяется с помощью следующего соотношения:
tинк= m p–u exp(B/T), (11)
где p – парциальное давление водорода, МПа; T – температура; m, u, B – некоторые константы, экспериментально определяемые для различных марок стали.
Длительность периода активных химических превращений определяем следующим образом [4]:
tкр = λ tинк, (12)
где: λ – параметр, являющийся мерой скорости химического взаимодействия водорода с конструкцией, показывающий, во сколько раз время до завершения периода изменений механических свойств материала отличается от продолжительности инкубационного периода, т.е. λ = tкр/tинк. Для стали 20 этот параметр принимает значения в пределах от 5 до 10 [4].
В том случае, если известно, что процесс обезуглероживания конструкции значительно короче времени до начала разрушения, можно использовать еще более простые модели воздействия водородосодержащей среды, в частности замену логистической кривой, описывающей плавное изменение механических характеристик, на ступенчатое их изменение от начального до конечного состояния.
Учет влияния напряжений на скорость коррозии и распределение давлений осуществляется так же, как и в обобщенной модели.
На основе упрощенной модели химического взаимодействия построена также модель оценки кинетики движения фронта обезуглероживания, целью разработки которой была возможность использования при аналитических расчетах. Однако она может использоваться только для некоторых несложных моделей пластин и неспособна учитывать нестационарный характер изменения тепловых, силовых нагрузок и давления водорода, а также влияние уровня напряжений на скорость движения фронта обезуглероживания.
В действительности моделей воздействия водородосодержащей среды значительно больше, наиболее распространенные из них подробно рассмотрены в работе [4]. Здесь же приведены только те модели, которые могут использоваться в составе разработанного авторами программного комплекса по расчету напряженно-деформированного состояния оболочечных конструкций [1, 2, 5, 7].
Технически учет влияния водорода на прочностные свойства конструкции реализован следующим образом. Для каждого из материалов конструкции введены механические характеристики двух его состояний – исходного и обезуглероженного. Если процесс обезуглероживания отсутствует (инкубационный период не завершен), то для расчета параметров напряженно-деформированного состояния берутся исходные механические характеристики материала. Если в материале уже происходит процесс обезуглероживания, но он еще полностью не завершен, то значения механических характеристик материала на этой стадии расчета вычисляются путем линейной интерполяции в зависимости от текущего значения m(t) между исходным и обезуглероженным состояниями. Если процесс обезуглероживания завершен, то для расчета берутся характеристики обезуглероженного материала.
В частности, с применением данной методики были проведены расчеты составной двухслойной конической оболочки, осевое сечение которой приведено на рисунке 2. Исследуемая конструкция состоит из трех участков. По форме меридиана координатной поверхности все три участка оболочка являются коническими. Левый торец оболочки жестко заделан, правый свободен. Толщина первого участка 40 мм (конусность 30°), третьего – 30 мм (конусность 10°). Второй участок является переходным (конусность 1,575°). Давление водорода осуществляется с внутренней стороны и постоянно во времени.
Внутренний слой, обращенный к водородосодержащей среде, изготовлен из материала 12Х18Н9Т, а внешний — из материала сталь 20. Давление водорода принималось равным 2,5 МПа. При этом в материале оболочки возникают как упругие деформации и деформации ползучести, так и мгновенные пластические. Расчеты приведены для температуры 500°С. Константы материала в соотношениях были взяты из [4], которые для стали 20 при данных условиях нагружения принимают следующие значения: k = 1,49 × 10–5 (МПа)*u, ч; u = 1,73; B = 13500; l = 1–10. Согласно данным, приведенным в работе [4], механически сплав 12Х18Н9Т под воздействием водорода практически не изменяется.
Рис. 2
Влияние уровня напряжений на скорость процесса обезуглероживания в данном расчете не учитывалось, поскольку для этого требуются дополнительные исследования.
В работе исследовалось влияние защитного слоя 12Х18Н9Т фиксированной толщины на длительную прочность оболочки, а также распределение меридиональных напряжений, деформаций зон пластичности и повреждаемости в моменты времени до разрушения, перед разрушением и непосредственно после него. Некоторые результаты приведены на рисунках 3–6.
Рис. 3
На рисунке 3 приведены кривые, показывающие распределение меридиональных (а) и окружных (б) напряжений на внутренней поверхности. Кривая 1 соответствует началу процесса нагружения (упругое решение), кривая 2 – моменту времени t = 1 ч (упругопластическое решение), кривая 3 – моменту времени t = 380 ч (сразу после начала разрушения). На этих графиках видно перераспределение напряжений по мере развития пластических деформаций и накопления повреждений. При этом наиболее нагруженными точками оболочки являются точки внутренней поверхности.
На рисунке 4 показан фрагмент осевого сечения оболочки для начала нагружения (а) и непосредственно перед разрушением (б), в котором наиболее интенсивно происходит накопление повреждений. В начале процесса нагружения на внутреннем слое появляется зона пластичности, размеры которой постепенно увеличиваются, и к моменту разрушения она составляет порядка 80% осевого сечения оболочки. Кроме того, во втором участке начиная с момента времени 260 ч также появляется небольшая зона пластичности. В третьем участке оболочки к моменту начала разрушения зоны пластичности не были выявлены. Начало процесса разрушения происходит в точке внутренней поверхности с меридиональной координатой S = 0,457 м (на рис. 4 место начала разрушения обозначено стрелкой).
Рис. 4
Рис. 5
Рис. 6
На рисунке 5 показано распределение меридиональных (а) и окружных (б) деформаций во внутреннем слое оболочки, а на рисунке 6 – степень поврежденности материала внутреннего слоя. Видно, что наибольшие значения полных деформаций развиваются в точках внутренней поверхности, близких к координате S = 0,457 м, где также накапливаются самые большие повреждения и впоследствии возникает разрушение.
Таким образом, исследуя историю изменения напряженно-деформированного состояния с учетом перемещения фронта обезуглероживания и повреждаемости материалов при ползучести и используя соответствующие критерии мгновенной и длительной прочности, можно оценить ее несущую способность и долговечность с учетом воздействия всех вышеназванных факторов.
Рецензенты:
Богомолов А.Н., д.т.н., профессор, проректор по научной работе ФГБОУ ВПО «Волгоградский государственный архитектурно-строительный университет», г. Волгоград;
Богданов Е.П., д.т.н., профессор кафедры «Информационные системы и технологии» ФГБОУ ВПО «Волгоградский государственный аграрный университет», г. Волгоград.
Библиографическая ссылка
Белов А.В., Поливанов А.А., Неумоина Н.Г. ИССЛЕДОВАНИЕ ВОЗДЕЙСТВИЯ ВОДОРОДА НА НЕСУЩУЮ СПОСОБНОСТЬ МНОГОСЛОЙНЫХ СОСТАВНЫХ ОБОЛОЧЕК В УСЛОВИЯХ ТЕМПЕРАТУРНОЙ ПОЛЗУЧЕСТИ // Современные проблемы науки и образования. – 2015. – № 2-1. ;URL: https://science-education.ru/ru/article/view?id=20790 (дата обращения: 16.04.2024).