Full Text

1. Постановка задачи Трубопровод внутренним радиусом R, толщиной стенки h свободно провисает на участке длиной L (рис. 1, а). Начало продольной координаты x принимается в середине этого участка. Участки -∞ < x < - L/2 и L/2 < x < ∞ находятся в грунте, свойства которого принимаются одинаковыми. Реакция грунта моделируется деформацией системы упругих пружин. При прогибе w(x) трубопровода со стороны системы пружин возникает распределенная сила, равная qsw, а при продольном перемещении u(x) трубопровода - распределенная сила tsu. Здесь через qs и ts обозначены жесткости системы пружин в поперечном и продольном направлениях. Реакция грунта на изгиб трубопровода более подробно рассматривается в [1-4]. В них, а также в [5-10] дается анализ прочности и устойчивости. Учитывается влияние давления p транспортируемой жидкости на изгиб трубопровода. Теоретическое и экспериментальное изучение этого влияния содержится в [11, 12]. Скорость V движения жидкости, а также изменение p на длине L не учитывается. Такое допущение может быть приемлемым, если ρf V2 << p, где ρf - плотность жидкости. Это условие выполняется, например, при данных p = 70 бар = 7·106 кг/м·с2, ρf = 900 кг/м3, V = 10 м/с. В газовых трубопроводах это условие выполняется практически всегда. Принимается допущение о повороте поперечного сечения так, что нормаль к осевой линии остается нормалью при изгибе и сохраняет свою круговую форму. Рассматривается статический линейный изгиб трубопровода. Нелинейный изгиб изучается, например, в работах [13-15]. Первая из них посвящена анализу сильного изгиба пролета трубопровода при подъеме его сосредоточенной силой, а вторая - описанию изгиба сверхпроводящего кабеля, представляющего трубу с концентрическими слоями из разных материалов для прокачки охлаждающего жидкого азота. В [15] рассмотрен нелинейный изгиб двухслойного трубопровода с определяемой из решения задачи длиной пролета (модель укладки трубопровода на морское дно). Рис. 1. Статический изгиб трубопровода Fig. 1. The static bending of pipeline Приведенная модель трубчатой балки широко используется для описания поведения трубопровода [2, 4, 8, 10]. Рассматриваемая расчетная схема справедлива не только в случае трубопровода и грунта, но и в случае конструкций трубопровод-емкость, где более податливая сплошная среда уменьшает концентрацию напряжений в месте соединения трубопровода со стенкой емкости (рис. 1, б). Различные подходы в обратных задачах упругой деформации стержней, балок, пластин рассматриваются, например, в [16, 17]. К самым последним публикациям по теме относятся работы [18-27], в которых изучаются свободные, вынужденные и параметрические колебания. Поведение трубы при ударном внутреннем давлении рассматривается в [27]. При этом влияние грунта не учитывается. В [24] дается обзор исследований. Линейное уравнение статического изгиба имеет вид [2, 11] (1) где E, ρ - модуль упругости и плотность материала; q - распределенные внешние силы на трубопровод, которые состоят из собственных весов жидкости и материала трубопровода единичной длины; g - гравитационное ускорение. Для x > L/2 в правой части (1) добавляется реакция грунта, равная -qsw. В линейном приближении в уравнении изгиба (1) отсутствует продольное перемещение, однако более строгий учет внутреннего перепада давления приводит к появлению такой зависимости. Величина P в (1) состоит из части, связанной с взаимодействием разности внутреннего и внешнего давлений p - ps и изменения угла поворота осевой линии (кривизны), и части, обусловленной кольцевой и осевой силами Nθ, Nx в стенке трубы под действием p - ps [11]. Здесь ps - давление на стенки трубы со стороны грунта. Первая часть π[R2 p - (R + h)2 ps] ≈ πR2 (p - ps ). В дальнейшем будем пренебрегать величиной по сравнению с радиусом R. Как видно из (1), изгибная деформация трубы зависит от условия ее растяжения-сжатия. Но в линейной постановке продольное растяжение-сжатие отделяется от задачи изгиба и не зависит от прогиба. Рис. 2. Продольная и кольцевая силы в стенке трубы Fig. 2. Longitudinal and annular forces in the pipe wall Для получения второй части P рассмотрим продольную силу 2πRNx, где Nx - усилие на продольную полоску единичной ширины стенки трубы толщиной h, как было определено выше (рис. 2). Из закона Гука , . (2) Окружная сила Nθ на единичную ширину кольца . Исключив деформацию в (2), находим . (3) Таким образом, с учетом (3) в уравнении (1) получаем (4) Так как в пределах 0 ≤ x ≤ L/2 на трубу не действуют силы в осевом направлении, то из уравнения равновесия 2πRdNx/dx = 0 следует постоянство Nx и в линейной постановке постоянство εx = du/dx. Отсюда εx = C и u = Cx + G (0 ≤ x ≤ L/2). Для участка x > L/2 деформацию εx определяем из равенства сил 2πRhE(dεx /dx) и tsu, или , , решение которого u = A exp(αx) + B exp(-αx) должно удовлетворять условиям u = 0 (x→∞), а также равенства продольных перемещений и сил левой и правой частей трубопровода в точке x = L/2. Определив A, B, C, G из указанных условий, выражения (4) приводим к виду (5) Из (5) видно, что учет радиального расширения трубы под действием внутреннего давления приводит к уменьшению продольной силы P по всей длине (когда это расширение не учитывается, в (5) ν не входит, поэтому в этом случае результат можно получить, положив ν = 0). Учет продольного перемещения на участке x > L/2 увеличивает здесь значение P, но уменьшает на участке 0 ≤ x ≤ L/2. В случае неглубокого залегания трубопровода в грунте и высокого внутреннего давления (ps << p) выражение P для обоих участков является одинаковым P = πR2p(1 - 2ν ). Такое допущение используется в [20]. С учетом (5) уравнению (1) придаем вид (6) 2. Случай высокого внутреннего давления в трубопроводе и неглубокого залегания в грунте Отбрасывая в (6) члены, содержащие ps, и вводя обозначения (7) второе уравнение (6) запишем в виде . (8) Для участка 0 ≤ ξ ≤ 1/2 в (8) отсутствует член 4β4w. Граничные условия состоят в равенстве нулю угла поворота и перерезывающей силы в средней точке трубы: , (9) и в равенствах прогибов, углов поворота, изгибающих моментов, перерезывающих сил в месте перехода от провисающей части (0 ≤ ξ ≤ 1/2) трубы к участку в грунте (ξ > 1/2): , , , . (10) Кроме того, ставятся условия , . (11) Отношение веса трубопровода q к жесткости пружин qs в поперечном направлении представляет собой равномерную по x осадку трубы в грунте. Решения уравнения (8), а также уравнения без члена 4β4w, удовлетворяющие условиям (9) и (11), имеют вид , (12) , . Таким образом, при удалении от провисающей части трубопровода прогиб стремится к равномерной осадке. Удовлетворяя (12) условиям (10), получаем , , , , , , . Зависимость относительного прогиба W = (24d/q)w в точке ξ = 0 от относительной жесткости грунта для различных значений относительного внутреннего давления μ показана на рис. 3. Рис. 3. Зависимость относительного прогиба W в середине (ξ = 0) пролета трубопровода от относительной жесткости грунта β Fig. 3. The dependence of W relative bending in the middle (ξ = 0) of the pipeline span on the relative stiffness of β soil С увеличением относительной жесткости β грунта происходит уменьшение относительного прогиба, причем тем быстрее, чем меньше относительное внутреннее давление μ. С увеличением μ возрастает относительный прогиб W. Таким образом, внутреннее давление вызывает дополнительный изгиб трубопровода. С уменьшением жесткости грунта происходит стремительный рост прогиба. То значение относительной жесткости грунта, когда прогиб возрастает неограниченно, назовем критическим βcr. Оно определяется из уравнения ∆(β, μ) = 0 или . (13) В отсутствие внутреннего давления (μ = 0) из (13) следует βcr = 0. Это означает, что неограниченное возрастание прогиба в линейной постановке может быть только при исчезающе малой жесткости грунта. При значении μ << 1 имеем βcr ≈ μ. Отметим, что эта оценка остается справедливой даже при нарушении условия μ << 1. Так, при μ = 0,3, 0,5 из (13) имеем βcr ≈ 0,303, 0,516. Другой предельный случай β→∞ соответствует большой жесткости грунта по сравнению с жесткостью на изгиб трубы. Тогда из (13) следует cos(μ/2) = 0 или μcr = π. С учетом (7) имеем значение критического внутреннего давления , , где PE - эйлерово значение продольной сжимающей силы в случае защемленных концов трубы (w = dw/dξ = 0 при ξ = -1/2, 1/2). Для большого, но конечного значения β примем μ = π - ε, где ε << π. Тогда из (13) находим ε ≈ 2π/(2 + β). Поэтому в этом случае μ ≈ πβ/(2 + β) и . Если, например, β = 8, то pcr снижается на 20 % по сравнению с предельным случаем абсолютной жесткости грунта. В случае отсутствия транспортируемой жидкости или малого внутреннего давления (μ = 0) выражение для относительного прогиба имеет вид (14) Зависимость W(β) по первому выражению (14) при ξ = 0 также показана на рис. 3. 3. Определение реакции грунта на изгиб трубопровода Определение параметра qs или относительной жесткости грунта β может быть осуществлено формально по замеру прогиба w трубопровода или деформации крайних волокон, например, в центре пролета (ξ = 0). Для этого могут быть использованы выражения (12), (14). Однако такой способ практически нереализуем. Здесь предлагается способ догружения сосредоточенной силой Q в средней точке ξ = 0 и замеров соответствующего дополнительного прогиба или деформации в этой же точке. Такая обратная задача, очевидно, сопряжена с повышенными требованиями к экспериментальным данным и с вопросом о ее корректности. Эти вопросы здесь не обсуждаются. Снова рассмотрим случай p >> ps, когда справедливо уравнение (8). Тогда относительно прогиба w, вызванного силой Q, имеем уравнения , (15) , где δ(ξ - 0) - дельта-функция. Граничные условия и условия сопряжения те же (9)-(11), но в них функция w(ξ) означает дополнительный прогиб. Решения уравнений (15), удовлетворяющие условиям (9) и (11), имеют вид , (16) , где . После удовлетворения условиям (10), находим , , , . В связи с тем, что определение параметра β по замеренным значениям прогиба или деформации трубы по выражениям (16) весьма затруднительно, рассмотрим более простой случай нулевого перепада давления (μ = 0). Тогда решения (16) приобретают вид , (17) . Зная экспериментальное значение прогиба [w] в точке ξ = 0, можно определить параметр , (18) , , (19) откуда qs = 4dβ4. Пусть надземная часть стальной трубы L=25 м, а ее внутренний диаметр 2R и толщина стенок h равны соответственно 0,309 и 0,008 м. E = 2·1011 Па, при этом d ≈ 51,27 Па. Определим прогиб трубы под действием груза Q = 500 H при жесткости грунта qs ≈ 16611,89 Па (β = 3). В соответствии с первой формулой в (17) получим, что прогиб w ≈ 0,13 м. В обратной задаче определяем жесткость грунта qs, если известно, что под действием дополнительного груза Q = 500 H прогиб трубопровода возрастает на [w] = 0,13 м. По формулe (19) W ≈ 0,25 м, а по формуле (18) β ≈ 3,04, поэтому qs ≈ 17589,61 Па (или qs /E = 8,79·10-8). Анализ показывает, что малые ошибки замера значения [w] приводят в обратной задаче к значительным ошибкам в определении значения жесткости грунта qs. Например, при [w] = 0,125 м значение qs ≈ 19712,52 Па (β ≈ 3,13), а при [w] = 0,135 м значение qs ≈ 15795,79 Па (β ≈ 2,96). На рис. 4 приведены значения β от W также при ненулевых µ (этот график построен при задании β). Рис. 4. Зависимость относительной жесткости грунта β от относительного прогиба трубопровода W для различного относительного давления транспортируемой жидкости µ Fig. 4. The dependence of the relative stiffness of soil β on the relative bending of W pipeline for different relative pressures of the transported fluid μ В силу принятого допущения о нормали к изогнутой оси поперечного сечения трубы продольная деформация нижнего волокна наружной поверхности ε = - (R + h)(d2w/dx2). С учетом выражения (17) при ξ = 0 получаем . E Рис. 5. Зависимость относительной жесткости грунта β от относительной продольной деформации трубопровода E Fig. 5. The dependence of the relative stiffness of soil β on the relative longitudinal deformation of pipeline E Зависимость относительной жесткости грунта β от безразмерной величины E = = (16L2d/Q (R + h)) [ε] показана на рис. 5. Заключение В идеализированной системе труба-транспортируемая жидкость-грунт определяется прогиб трубопровода в зависимости от геометрических и механических характеристик трубы, от плотности и давления жидкости, упругой реакции грунта. Учитывается влияние на изгиб взаимодействия внутреннего перепада давления и кривизны осевой линии, а также возникающего продольного усилия в результате осесимметричной деформации трубы. Рассмотрен случай большого внутреннего давления в трубопроводе и неглубокого его залегания в грунте, когда задача существенно упрощается. С увеличением внутреннего перепада давления происходит рост прогиба пролета трубопровода. Определяются критическое значение жесткости грунта, ниже которого в линейной задаче происходит неограниченное возрастание прогиба, а также значение внутреннего перепада давления, выше которого происходит такое возрастание. В обратной задаче по экспериментально определенным значениям прогиба или деформации трубы определяется значение жесткости грунта.

About the authors

M A Ilgamov

Mavlyutov Institute of Mechanics RAS

A A Yulmukhametov

Mavlyutov Institute of Mechanics RAS

References

Statistics

Views

Abstract - 277

PDF (Russian) - 95