CC BY
134
УДК 622.692.4
В. И. Ларионов, П. А. Новиков, А. К. Гумеров
АНАЛИЗ НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ ТРУБОПРОВОДА НА УЧАСТКАХ С КАРСТАМИ
Приведен расчет напряженно-деформированного состояния линейной части магистральных трубопроводов, выполненный методом конечных элементов. Рассмотрены перемещения трубопровода в грунте в вертикальной и горизонтальной плоскостях, вдоль оси трубопровода, а также сопротивление грунта в разных плоскостях, независимых друг от друга.
E-mail: [email protected].
Ключевые слова: трубопровод, магистральный, карсты, деформация,
анизотропный отпор, нелинейные пружины.
Подземные магистральные трубопроводы — это протяженные объекты, которые проходят не только через стабильные районы с устойчивыми грунтами, но и районы с активными проявлениями грунтовых изменений, в том числе через болота, реки, сейсмические разломы, многолетнемерзлые грунты, участки с карстовыми проявлениями, склоны и овраги. Пересекая карстовые полости и провалы, нефтепровод деформируется вместе с грунтом под действием собственного веса и веса грунта, а также под воздействием перепада температуры и внутреннего давления [1]. Схемы состояния трубопровода в зоне карстовых полостей, характерных для термокарстов, приведены на рис. 1, а, в зоне карстовых провалов — на рис. 1, б.
Рассмотрим случай, когда трубопровод проходит над карстовой воронкой.
Расчетная схема модели с указанием всех нагрузок и воздействий представлена на рис. 2.
Решение основывается на системе дифференциальных уравнений
где V и т — вертикальное и продольное смещение трубы; (г) и дг (г) — вертикальная и продольная нагрузки на трубу; Е — модуль упругости металла трубы; N (г) — продольная (осевая) сила, Н; 7у (г) — момент инерции поперечного сечения трубы относительно вертикальной оси у, совпадает с Зх(г).
d2w qz (z)
(1)
dz2 nDES^
Рис. 1. Схема состояния трубопровода в зонах термокарстов (а) и карстового провала (б)
У 91 А д> и
сТГ^тхШ ШШШШШ ГТТТТБ в
AT Т ТIТ1Т Т ТТ Т11 |А
<12 , Ч
1
и и L,
ь
Рис. 2. Расчетная схема участка трубопровода в зоне карста
Первое уравнение системы (1) — это уравнение продольно-поперечног< изгиба упругой балки в плоскости уг), а второе — уравнение, описывающее состояние трубопровода при продольном сдвиге по оси г.
Геометрические характеристики труб, наиболее часто используемые в расчетах магистральных трубопроводов, определяются следующим образом:
^(г) = п5т(В - 5т) (2)
— площадь поперечного сечения стенки трубы;
, = ^ - ев - 25т)4] (3)
64
— осевой момент инерции поперечного сечения трубопровода, где В
— наружный диаметр трубопровода; 5т — толщина стенки трубы.
Рассмотрим расчетную схему (см. рис. 2), на которой обозначено: СВ — расчетный участок длиной Ь; АВ — участок с карстом длиной Ь\; С А и В В — участки вне карстовой зоны длиной Ь2.
Определим основные нагрузки, действующие на трубопровод [3]. Сила, действующая на трубопровод вне области карста,
С1 = Ссв + Сгр,
где ссв — собственный вес трубы с продуктом на единицу длины, равный
Ссв(у) = - ^ - (В - 2^т)2] - - 2^т)2 - ПТи^и (П + 6И) (4)
(здесь 7г, 7н, 7и — удельные веса металла трубы, нефти и изоляции; ¿т, 6и — толщина стенки трубы и изоляции; П — внешний диаметр трубы);
Сгр = -7грВ (к + 0,1073В)
— масса грунта над трубой на единицу длины (7Ф — удельный вес грунта над трубой; к - глубина залегания трубы, т.е. высота от дневной поверхности земли до верхней образующей трубы).
Сила реакции грунта (в исходном состоянии) д2 = -д1. На участке АВ (в области карста) реакция грунта д2 = 0. Сила, действующая на трубу в области карста,
!дсв, если над трубой нет грунта;
Ссв + Сгр, если над трубой есть грунт.
Если нижняя образующая трубы достигнет дна карста, то появится дополнительная сила д4, уравновешивающая силу с3:
С4 = -Сз.
Вне области карста (участки С А и ВП) труба также должна двигаться вниз под действием силы с1 и внутренних сил, появляющихся вследствие изгиба трубы, при этом возникает дополнительная сила, препятствующая этому движению:
Сгр(у) = -Су AvП,
где Су — коэффициент постели грунта при движении трубы вниз; Av — смещение грунта вниз (при смещении вниз величина Av отрицательная); П — внешний диаметр трубы.
Сила сгр(у) не может превышать предельного значения спр(у). Продольная сила, действующая на каждый элемент трубы, определяется по формуле
N.. = (цакц(р) - аАгЕ + Ае2г)Г, где ц — коэффициент Пуассона металла трубы;
_ П - 25т
акц(р) = 2£т
— кольцевое напряжение в трубе; а — коэффициент температурного расширения металла трубы; А£ — температурный перепад, т.е. разность между температурой эксплуатации трубы и температурой замыкания трубопровода; Е — модуль упругости металла трубы;
д 1 - 10
Аегг = —— ¿0
— относительное удлинение ¿-го элемента трубы; 10 — начальная длина элемента трубопровода; I — измененная длина элемента трубопровода; Г — площадь поперечного сечения трубы.
Изменение продольных размеров участков трубопровода вызывает продольное смещение трубы относительно грунта, при этом появляются силы, препятствующие этому смещению:
Сгр(2) = -пС,ш Атг П;
— коэффициент постели при продольном перемещении грунта; Аыг — продольное смещение грунта; П — внешний диаметр трубы.
Сила дгр(2) не может превышать предельного значения дпр(г).
В качестве исходных данных принимается прямолинейный участок трубопровода. Внутри участка строится конечно-элементная сетка с размером элементов, равным к.
Затем организуется итерационный процесс (прогонка) (] = = -1, 0, 2, 3,... ,р) в пределах всех узлов конечно-элементной модели (г = 1, 2,3,..., п). Определяется результирующее смещение каждого узла г с помощью формул, полученных при решении системы (1),
vj+1 = AVi j + vj+i) + BVi j + vj+J + Cj;
wj+1 =0,5 (W_i + wj+i + C'm),
(5)
где коэффициенты определяются следующим образом:
А = 4Е1Х + N к2 В = -Е,1Х
т -тт, Т . „Л 2 ; ™ ЙР 7 I ОЛГ-?
6E Jx + 2Nj h2 6EJX + 2Nj h2 ■
j
%
сj = qyih ; сj = jjh2
Ш
6EJx + 2N,j h2' nDEÖT1
здесь г — номер узла сетки; ] — номер этапа итерации; N — суммарная продольная сила, действующая в г-м узле на ]-й итерации и вычисляемая по формуле
N = (цакц - аАЬЕ + А^г )Г; (6) А ^ = У(^ - 1) 2(к-!+-- 1 (7)
— относительное удлинение оси трубопровода на конечном элементе
2уг
■ „ 7
в г-м узле при ]-и итерации; дуг — суммарная вертикальная нагрузка,
• „ 7
прикладываемая в г-м узле на ]-и итерации; — суммарная продольная нагрузка, прикладываемая в г-м узле на ]-И итерации, причем
дУг = 0; 7 = (^акц - аА1Б(8)
— на участках С А и ВО;
7 = -дсв - 7 = - аАгЕ(9)
— на участке АВ.
В качестве примера рассмотрим первую итерацию
У1 = А1г (у— + ) + В1г + + т1 = 0,5 (т°г-1 + + С^) ;
дуг = дсв дгр;
N = (^акц - аМЕ + Де^)*1; д 1 _ - у1-!)2 + (2^ + т!+1 - т_1)2
Де1г = й 1,
а также вторую итерацию
У2 = А2г (У_1 + У+1) + В2 (У!_2 + У1+2) + С2т;
т2 = 0,5 (т1-! + т!+1 + С2шг) .
На каждоИ последующей итерации смещение Ду7+!, входящее в формулу определения реакции грунта, будет вычисляться из соотношения
ДУ7+! = ДУ^ - , (10)
где Ду7 — смещение, полученное на ]-И итерации в г-м узле, — перемещение, вычисленное на ]-И итерации в г-м узле.
После определения смещений по всем узлам находим по формуле следующую величину:
п
ББ7 = ^(|у7| + |т7|) . (11)
г=1
Сравним значение Б Б7 со значением Б Б7-1, полученным на предыдущем шаге. В случае если полученные значения отличаются между собой на величину Д, меньшую заданной точности е, то расчет заканчивается. В противном случае переходим к следующей итерации.
Повторяя последовательные приближения в зависимости от требуемой точности решения, получаем такое состояние конечно-элементной
модели, когда смещения всех узлов практически не изменяются. Это означает, что получено искомое решение системы (1).
Далее представлена блок-схема, описывающая алгоритм работы программного комплекса, рассчитывающего напряженно-деформируемое состояние (НДС) трубопровода на участке карста (рис. 3).
Ввод исходных данных Характеристики трубопровода; Свойства транспортируемого продукта; Характеристики грунтов; Состояние трубопровода; Размеры участка с карстом
Рис. 3. Алгоритм расчета НДС трубопровода на участке карста
Рассмотрим задачу о расчете НДС трубопровода с наружным диаметром В = 1220мм и толщиной стенки трубы 22мм (см. рис.2). Температурный перепад "эксплуатация-укладка" составляет -20 °С. Длина карстового участка (АВ) равна 50 м; сдвиг вдоль трубопровода составляет 5 м; глубина карста — 5 м; глубина залегания трубопровода — 1м. Модуль упругости материала равен 2,06 • 105 МПа; предел текучести материала трубы — 500 МПа; коэффициент температурного расширения — 1,2 • 10-5; рабочее давление в трубе — 5 МПа. При расчетах использовали удельные веса, Н/м3: стали — 77000; нефти — 8700; изоляции — 11000; грунта - 16800. Коэффициент постели грунта составляет 1,32МПа/м; предельная реакция грунта — 6237 Н/м.
Результаты расчетов приведены на рис. 4. Видно, что в районе карста основные показатели НДС трубопровода составляют: кольцевое напряжение 133,6 МПа; продольное осевое напряжение 89,53 МПа; продольная осевая деформация 4,3462• 10-4; на отметке 110м — максимальное смещение трубы 0,21 м, минимальный радиус изгиба 1,4028 х х 103 м, максимальная кривизна трубы 7,1286 •Ю-4 м-1, максимальный изгибающий момент 2180000 Н-м, максимальное напряжение изгиба трубы 89,58 МПа, минимальное и максимальное суммарные продольные напряжения на стенке трубы — 1,54 МПа и 180,69 МПа, максимальная интенсивность напряжений 162,36 МПа, минимальная и
Рис. 4. Графики зависимости напряжения изгиба (а) и интенсивности напряжений (б) трубы от расстояния
максимальная суммарные деформации стенки трубы 7,4733 • 10-6 и 8,7716 • 10-4; максимальная перерезывающая сила 0,33 МН на отметке 85 м; коэффициент запаса по пределу текучести 3,08.
Выводы: 1. Предложен метод определения НДС магистрального трубопровода на участках с карстами.
2. Разработан алгоритм решения уравнений и предложена его программная реализация.
3. Разработанная модель взаимодействия трубопровода с грунтом на участке с карстом позволяет исследовать и оценить влияние опасных геологических процессов на НДС и надежность участков трубопровода.
Работа, результаты которой использованы в настоящей статье, выполнена по заказу ОАО "АК "Транснефть " при финансовой поддержке Министерства образования и науки РФ.
СПИСОК ЛИТЕРАТУРЫ
1. Ш а м м а з о в А. М., Зарипов Р.М., Чичелов В. А., Коробков Г.Е. Расчет магистральных газопроводов в карстовой зоне. - Уфа: Гилем, 1999. -213 с.
2. Б о р о д а в к и н П. П., Синюков А. М. Прочность магистральных нефтепроводов. - М.: Недра, 1984. - 245 с.
3. АйнбиндерА. Б., Камерштейн А. Г. Расчет магистральных трубопроводов на прочность и устойчивость. - М.: Недра, 1982. -340 с.
Статья поступила в редакцию 25.06.2012 Валерий Иванович Ларионов окончил Военно-инженерную ордена Ленина Краснознаменную академию им. В.В.Куйбышева в 1972г. Д-р техн. наук, профессор, заместитель директора по научной работе НОЦ ИЭС МГТУ им. Н.Э. Баумана. Автор более 200 начных работ в области теории рисков и промышленной безопасности потенциально опасных объектов.
V.I. Larionov (b. 1941) graduated from The Order of Lenin and Red Banner Military Engineering Academy n. a. V.V. Kuibyshev in 1972. D. Sc. (Eng.), professor, deputy director for scientific work of the Scientific and Educational Center of Extremal Situation Study of the Bauman Moscow State Technical University. Author of more than 200 publications in the field of theory of risks and industrial safety of potentially dangerous objects.
Павел Александрович Новиков окончил МГТУ им. Н.Э. Баумана в 2008 г. Научный сотрудник НОЦ ИЭС МГТУ им. Н.Э. Баумана. Автор более 10 научных статей в области промышленной безопасности и эксплуатации опасных производственных объектов.
P.A. Novikov graduated from the Bauman Moscow State Technical University in 2008. Researcher of the Scientific and Educational Center of Extremal Situation Study of the Bauman Moscow State Technical University. Author of more than 10 publications in the field of industrial safety and exploitation of dangerous industrial objects. Айдар Кабирович Гумеров окончил Уфимский государственный авиационный технический университет в 2006 г. Канд. техн. наук, старший научный сотрудник ГУП ИПТЭР Республики Башкортостан. Автор более 30 научных работ в области эксплуатации и промышленной безопасности трубопроводного транспорта. A.K. Gumerov graduated from the Ufa State Aviation Technical University in 2006. Ph. D. (Eng.), senior researcher of the Institute of Problems in Energy Resource Transport of the Republic of Bashkortostan. Author of more than 30 publications in the field of exploitation and industrial safety of pipeline transport.
