Математическое моделирование эффекта «течь перед разрушением»

Бесплатный доступ

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

Сквозной дефект, разрушение, математическое моделирование

Короткий адрес: https://sciup.org/148198983

IDR: 148198983

Текст научной статьи Математическое моделирование эффекта «течь перед разрушением»

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

Анализ вариантов развития поверхностных дефектов на стенках трубопроводов. Определение условий разрушения конструкции трубопровода возможно на основе учета таких свойств материала, как вязкость разрушения, действующих напряжений и размеров трещиноподобных дефектов. Будем полагать, что на внутренней поверхности возникает несквозное повреждение типа трещины полуэллиптической формы глубиной l и полудлиной a (рис. 1).

кI = ^fn y

При этом для определения параметра Y, зависящего от геометрии тела, формы трещины и вида напряженного состояния, можно воспользоваться формулами, приведенными в литературе [13]. На основе решения задачи о поле напряжений вокруг эллиптического дефекта, полученного Грином и Снеддоном, а также Ирвином, в работах [4-6] приведено выражение для коэффициента интенсивности напряжений кI =

^[n i ф

24 l sin ф + -ycos ф I

a2         J

Здесь Ф – эллиптический интеграл второго рода:

П (     2   /2       Л12

a — l ф =   1---—sin ф dф

0 1 a 2 J

где l и a определяются из рисунка 1.

Величину Ф можно представить в виде ряда

Рис. 1. Схематизация несквозной трещины в стенке трубы

n 1 a 2 - I 2    3 ( a 2 - I 2 )

I

2 [ 4 a 2     64 ( a 2 J

Даже при значениях l / a , близких к нулю, третий член в квадратных скобках дает вклад, относительная величина которого оценивается 5%. Поэтому в большинстве случаев им пренебрегают. При этом будем иметь

Зависимость коэффициента интенсивности напряжений от приложенных напряжений и характерного размера дефекта имеет вид

и

3п n 12 ф=--+--у

88 a 2

^П1  (1

KT =            sin ф + —

I                          2

31l2

п + n^v 88 a 2

, 1 74 cos2 ф I

. (6)

графиками, приведенными на рис. 3. С учетом (9) и (10) получим

Соотношения (2) и (6) с небольшими изменениями можно применять к полуэллиптиче-ским поверхностным повреждениям – трещинам (рис. 2).

Рис. 2. Часть стенки трубопровода с осевым повреждением

K],    .= 1,12ai n i

1(ф=П2 )           \ Q .

;

Ku    = 1,12 ^1— l

I ( ф = 0 )     , Q a

.

Заметим, что в случае полукругового дефекта при l / a =1 и без учета пластических деформаций, т.е. при Q=Ф= π /2 будем иметь

к I = —^Jni n

.

Это справедливо при любых значениях угла φ . Отношение значений коэффициентов интенсивности напряжений при φ = π /2 и φ =0, определяемых равенствами (11) и (12), равно

Величина K I , как это нетрудно установить, меняет свои значения вдоль фронта трещины. На конце малой полуоси ( φ = π /2) интенсивность напряжений максимальна, на конце большой полуоси ( φ =0) она принимает наименьшее значение. Эти значения равны

K

1 ( Ф=^ 2 ) ) = £

K Ч ф П ) =

KI (ф=0) =

^ni

Ф . ;

е4П

Ф

Поверхностная трещина оказывает более сильное влияние на прочность, чем внутренняя. Большее влияние поверхностных трещин на прочность элемента конструкции объясняется отсутствием связей с наружной стороны [7]. Учет этого обстоятельства производится с помощью корректировочного коэффициента [4], равного примерно 1,12 и учитывающего наличие свободной тыльной поверхности. Наличие в вершине трещины пластической деформации приводит к тому, что трещина ведет себя так, как если бы она была немного длиннее, чем на самом деле [7, 8]. В этом случае получим окончательное выражение для K I в виде

K 1 ( ф = 0 )      v l ,                                (14)

для поверхностныхи внутренних повреждений

1,12^Л1 Г . 2     l2    2 )1/4

K =-------- sin ф + cos ф\

I Q I          a 2        ) , (9)

где

Q = I Ф 2 - 0,212 ^ 2 2

*             /ст т                 (10)

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

Состояние перед разрушением может развиваться по схеме, изображенной на рис. 4. Предположим, что поверхностный дефект полу-эллиптической формы с размерами 2 a и l расположен на внутренней поверхности трубопровода и ориентирован нормально по отношению к направлению кольцевых напряжений (рис. 2). Из-за агрессивности среды и циклических нагрузок эта трещина будет расти. На основании (15) следует, что поверхностный дефект полуэллипти-ческой формы растет в направлении, параллельном малой оси эллипса, где φ = π /2, до тех пор, пока дефект не примет полукруглую конфигурацию. Далее трещина будет расти как непрерывно распространяющаяся полуокружность до тех пор, пока не выйдет на внешнюю поверхность

трубопровода. Трещина будет расти, распространяясь на оставшиеся неразрушенные объемы (заштрихованная зона) по всей толщине (рис . 4).

Рис. 4. Схема роста полуэллиптического поверхностного повреждения в полукруглое

После достижения трещиной внешней поверхности трубопровода происходит утечка перекачиваемого продукта, что делает реальным обнаружение этой трещины. Однако существует и вероятность нестабильного развития трещины. При наличии сквозной трещины длиной 2 a ее края под действием внутреннего давления имеют тенденцию к повороту, т.е. выгибаются наружу, как это показано на рис. 5. Ниже это обстоятельство учитывается путем введения коэффициента увеличения интенсивности напряжений Фолиаса [9].

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

Условия реализации критерия «течь перед разрушением». В работах Ирвина и Сро-ули [10, 11] был предложен упрощенный критерий «течь перед разрушением». Он основан на предположении о том, что нестабильность, предшествующая разрушению, в трещине размера 2 t возникает при напряжении, равном пределу текучести σ т . Предполагается также, что при таком высоком напряжении имеет место плоское напряженное состояние. Следовательно, остановка может произойти при условии, что интенсивность напряжений при плоском напряженном состоянии удовлетворяет условию

K c ^ ® T V n ( t + r p )

Здесь r p – радиус зоны пластической деформации в вершине трещины. Значение r p принималось равным

l а l=t и σ=σт. Тогда или

K c2 i®. T

Рис. 5. Схема выгибания зоны трубопровода с трещиной

2 n

,

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

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

K Ict = 1,12 M K

.

В этом выражении используются следующие обозначения: KIct – коэффициент интенсивности напряжений для материала с трещиной, распространяющейся в направлении толщины (рис. 6); MK – коэффициент увеличения интен- сивности напряжений Кобаяши [12], учитывающий приближение фронта трещины к свободной поверхности. Для определения величины MK можно воспользоваться графиком, приведенным на рис. 7.

провода (рис. 5), коэффициент интенсивности напряжений следует задать соотношением

pR

K I = M F ^9 Пла = —

1 + 1,16— n a Rt J

Здесь M F – коэффициент увеличения интенсивности напряжений, который был получен теоретически Фолиасом [19]. Смысл этого коэффициента рассмотрен выше. Его величина определяется согласно (20) по формуле

Рис. 6. Схема реализации критерия «течь перед разрушением»

M, = . 1 + 1,61 —

F            Rt

.

Окружное напряжение σ φ в тонкостенной конструкции равно pR / t . Добавляя к этому напряжению величину p , мы учитываем давление внутри трещины. Используя разложение функции Ф в ряд (?4.4), можно получить для давления, необходимого для процесса нестабильного роста поверхностного дефекта, следующее выражение:

Согласно уравнению (20) давление p 2 , вызывающее нестабильное развитие сквозной трещины, задается следующим соотношением :

p 2 =

K Ic

R

Mv — Паа

F

P l =

П

3 +- K I а

Ict

1,12 M K

R

t 7

nl

Рис. 7. Коррекция Кобаяши, учитывающая близость свободной поверхности к фронту трещины

Здесь K I c – вязкость разрушения материала с трещиной, распространяющейся в продольном направлении (см. рис. 6).

Остановка трещины произойдет в случае, если давление, необходимое для распространения сквозной трещины длиной 2 a больше, чем давление, необходимое для нестабильного роста поверхностного повреждения, глубина которого l . Следовательно, критерий утечки до разрушения следует из неравенства p 2 > p 1 , которое с помощью (19) и (20) можно представить в виде

K Ic

R

MT~ Паа

F t

8 П ( 3 + a2 I K ’ct 8 a

>--^-^—

1,12 mk | 1 + R I v a

K

Критериальное выражение (23) можно упростить, если пренебречь единицей по сравнению с отношением R / t . Далее, величину коэффициента Фолиаса M F можно принять равной 1. Это справедливо в случае, когда в выражении (21) величина

Поверхностный дефект перерастет в сквозную трещину с размером 2 a . Для сквозной трещины, имеющей длину 2 a в стенке трубо-

1,61 — <<  1 Rt

В этом случае на величину давления p 2 коррекция M F не влияет. С учетом сказанного (23) можно представить в виде

K

Ic

K

Ict

1 п С l2 ) --3 + -Т I

9 M a 2 K

l

a

Для поверхностных дефектов, характеризуемых малым отношением l / t , M K ≈1. В этом случае правую часть выражения (25) вычислить несложно; графическое изображение данного неравенства представлено на рис. 8 (верхняя кривая). Для пологих дефектов при l / t ≈1 значение M K ≈2. Поэтому нижняя линия на рис. 8 начинается при значении отношения K Ic / K Ict , составляющем примерно 0,5 от соответствующего значения для верхней линии. Для полукруглых поверхностных дефектов M K близко к единице независимо от l / t . Таким образом, в правой части диаграммы обе линии совпадают, что дает возможность найти нижнюю границу для l / t≈ 1. Условие утечки до разрушения более легко выполняется в материалах, обладающих значительной анизотропией. На практике анизотропия вязкости разрушения редко бывает больше, чем 2. Это обстоятельство затрудняет подбор условий «течь перед разрушением» для трещин, глубина которых меньше величины, определяемой отношением l / t≈ 1.

Рис. 8. Область выполнения критерия «течь перед разрушением»

Ситуация существенным образом меняется в случае, когда сквозная трещина находится в плоском напряженном состоянии, а в районе поверхностного дефекта – повреждения имеет место плоская деформация. Последнее весьма вероятно, так как из-за изогнутых берегов трещины пластическое ограничение велико. С учетом сказанного по оси ординат на рис. 8 следует откладывать отношение Kc/KIct, где Kc – вязкость разрушения материала для данной толщины элемента конструкции, т.е. t. Отношение Kc/KIct остается достаточно большими вплоть до сравнительно значительных толщин t. В этих случа- ях добиться реализации утечки до разрушения возможно.

В работе [13] приведены результаты исследования влияния варьирования химического состава сплава и примесей в нем в пределах стандартного состава на вязкость разрушения крупногабаритных прессованных полуфабрикатов. Эти данные свидетельствуют о том, что варьирование состава металла дает различные значения отношения K Ic / K Ict в пределах 1,061,58. Это подтверждает возможность соблюдения условий реализации критерия «течь перед разрушением». В работе [6] подтверждается высказанное выше предположение о возможности замены отношения K Ic/ K Ict в условии (25) отношением K c/ K Ict. Значения этих отношений достаточно велики вплоть до сравнительно больших толщин t . В этом можно убедиться, если воспользоваться представлением K c в виде эмпирического соотношения Ирвина, приведенного в [5]:

Ic

1 + 1,4

t I ° T )

Правомерность использования отношения K c / K Ict в условии (25) подтверждается и в работе [6]. «Оказывается, – утверждает автор, – что в любом случае для устойчивого образования таким путем сквозной трещины обычно требуется наличие прочностной анизотропии, когда K c / K Ict >1. Если толщина стенки достаточно мала, то данное условие «течи перед разрушением» может с хорошей степенью точности быть удовлетворено даже для изотропного материала».

Список литературы Математическое моделирование эффекта «течь перед разрушением»

  • Панасюк, В.В. Экспериментальное определение вязкости разрушения конструкционных материалов: Рекомендации по стандартизации определения вязкости разрушения/В.В. Панасюк, А.Е. Андрейкив, С.Е. Ковчик//Физ.-хим. механика материалов. -1976. -Т. 12, № 2. -С. 10-17.
  • Прикладные вопросы вязкости разрушения/Пер. с англ. под ред. Б.А. Дроздовского. -М.: Мир, 1968. -552 с.
  • Расчеты и испытания на прочность в машиностроении. Методы механических испытаний металлов. Определение характеристик вязкости разрушения (трещиностойкости) при статическом разрушении: Методические указания РД 50-260-81. -М.: Госстандарт СССР, 1982. -55 с.
  • Броек, Д. Основы механики разрушения: Пер. с англ./Д. Броек. -М.: Высш. школа, 1980. -368 с.
  • Херцберг, Р.В. Деформация и механика разрушения конструкционных материалов: Пер. с англ./Р.В. Херцберг. -М.: Металлургия, 1989. -575 с.
  • Хеллан, К. Введение в механику разрушения: Пер. с англ./К. Хеллан. -М. МИР, 1988. -364 с.
  • Васильченко, Г.С. Практическое применение механики разрушения для оценки прочности конструкций/Г.С. Васильченко, П.Ф. Кошелев. -М.: Наука, 1974. -148 с.
  • Браун, У. Испытания высокопрочных металлических материалов на вязкость разрушения при плоской деформации: Пер. с англ./У. Браун, Дж. Сроули. -М.: Мир, 1972. -246 с.
  • Folias, E.S. A finite crack in a pressured cylindrical shell/E.S. Folias//Int. J. Fracture Mech. -1965. -V. 1. -P. 104-113.
  • Irwin, G.R. Progress in the development of crack toughness fracture tests/G.R. Irwin, J.E. Srawley//Materialprüfung. -1962. -V. 4. -P. 1-11.
  • Irwin, G.R. Fracture of pressure vessels/G.R. Irwin//Material for missile and spacecraft. -Mc Graw-Hill, 1963. -P. 204-229.
  • Kobayashi, A.S. Approximate stress intensity factor for an embedded elliptical crack near two parallel free surfaces/A.S. Kobayashi, M. Zii, L.R. Hall//Int. J. Fracture Mech. -1965. -V. 1. -P. 81-95.
  • Микляев, П.Г. Кинетика разрушения//П.Г. Микляев, Г.С. Нешпор, В.Г. Кудряшов. -М.: Металлургия, 1979. -279 с.
Еще
Статья научная