CC BY
58
УДК 539.3
Моделирование процесса деформирования и оценка несущей способности системы грунт — тонкостенная конструкция*
Р.А. Каюмов1, Ф.Р. Шакирзянов1, С.С. Гаврюшин2
1 КГАСУ, 420043, Казань, Российская Федерация, ул. Зеленая, д. 1.
2 МГТУ им. Н.Э. Баумана, 105005, Москва, Российская Федерация, 2-я Бауманская ул., д. 5, стр. 1.
Modeling of the deformation process and evaluation of the bearing capacity of a thin-walled structure in the ground
R.A. Kayumov1, F.R. Shakirzyanov1, S.S. Gavryushin2
1 Kazan State University of Architecture and Engineering, Zelenaya st., 1, 420043, Kazan, Russian Federation.
2 Bauman Moscow State Technical University, building 1, 2-nd Baumanskaya str., 5, 105005, Moscow, Russian Federation.
Г(Т|1 e-mail: [email protected], [email protected], [email protected]
Оценка состояния, прочности и долговечности тонкостенных полиэтиленовых трубопроводов в грунте в настоящее время является актуальной проблемой. В статье рассмотрена задача численного анализа предельной нагрузки, вызывающей осадку грунта, взаимодействующего с размещенной в нем полиэтиленовой гофрированной трубой. Для описания процесса деформирования использована упруго-вязко-пластическая модель системы грунт — тонкостенная конструкция, учитывающая дилатансию и изменения механических характеристик грунта во времени, а также влияние температуры на механические характеристики полиэтилена. Приведены результаты численных экспериментов, выполненные с использованием разработанного авторами комплекса конечно-элементных программ. В результате исследования выявлены закономерности влияния геометрических и механических характеристик на процесс осадки и предельную нагрузку для системы грунт—полиэтиленовый трубопровод.
Ключевые слова: механика грунтов, ползучесть, пластичность, предельная нагрузка, дилатансия, метод конечных элементов.
О Estimating the strength and durability of thin-walled plastic pipes in the ground is of particular interest at present. This paper deals with the numerical analysis of the ultimate load causing displacements of the ground in which a polyethylene corrugated pipe is located. An elastic-viscous-plastic model was developed to describe deformations of a thin-walled structure in the ground. The model takes into account the dilatancy and changes in mechanical characteristics of the ground over time and the effect of temperature on the mechanical properties of polyethylene. Numerical experiments were performed using a finite element code developed by the authors. The results of study show the influence of the geometrical and mechanical characteristics of the plastic pipe on the displacements of the ground. The ultimate load causing this phenomenon is determined.
Keywords: soil mechanics, creep, plasticity, ultimate load, dilatancy, finite element method.
* Работа выполнена при финансовой поддержке РФФИ на 2013-2014 г. (проект № 13-01-90708 мол_рф_нр).
Сильная изношенность коммуникаций ЖКХ, тонкостенных трубопроводов в грунте требует оценки их состояния, прочности и прогнозирования долговечности и, в случае необходимости, замены их на трубы из современных, прочных и долговечных материалов. Однако в материалах типа полиэтилен от перепада температур сильно меняются прочностные и жесткостные характеристики, что может обусловить значительное изменение напряженно-деформированного состояния (НДС) и несущей способности системы грунт — тонкостенная конструкция. Поэтому необходима оценка влияния уменьшения жесткости и прочности полиэтилена при нагревании на несущую способность и осадку полиэтиленовой трубы в грунте.
С другой стороны, система грунт — тонкостенная конструкция сама изменяется с течением времени от деформаций пластичности, ползучести, дилатан-сии, что может привести к значительному изменению НДС и несущей способности [1-5]. Поэтому их также необходимо учитывать. В данной работе эта проблема рассматривается для полиэтиленовой гофрированной трубы низкого давления, находящейся в грунте.
Цель работы — определить несущую способность системы грунт — тонкостенная конструкция.
Запишем вариационное уравнение принципа возможных перемещений
¡¡¡оу Ьг^йУ = ¡¡{ О, ЬщйУ + Ц Ры5щйБ, (1)
V V X
где Q¡, Р„1 — объемные и поверхностные нагрузки.
Связь между деформациями и перемещениями будем считать линейной, используя соотношения Коши:
Ц = 2 (и,) + и» ).
(2)
Прочность грунта оценим по критерию Друк-кера — Прагера [6]
Р = аЬ+4ь -к = 0,
(3)
+ т2 + т2 + т2
Ог
= [(ах - О у )2 + (О у -О? )2 + (О? -Ох )2 ■
+ 6(тХу +т2уг +т1х)]1/2.
(4)
^ = 8у + И* аа. (5)
& 1 + АхЕ0 ' 1 + Л2гс к Здесь ц, п, а, р — реологические коэффициенты; А1, А2 — параметры упрочнения; №к1 — тензор, связывающий компоненты тензора напряжений и девиатор скорости деформаций ползучести.
Пластические деформации грунта описываются ассоциированным законом течения [7]
йХ дБ
йб йб доц
(6)
Предельную нагрузку и осадку системы оценим по теории упруго-вязко-пластического деформирования с учетом температуры, так как в полиэтиленовых трубопроводах при колебаниях температуры воды могут возникать значительные деформации. Температурные деформации определим из закона линейного температурного расширения
4 =РЛТ, (7)
где р — коэффициент температурного расширения; ДТ — перепад температуры.
Тогда полная деформация будет определяться выражением
еу =Е!) +£р + г\] + Е!), (8)
в котором скорости пластических деформаций 6Р и скорости деформаций ползучести 6с определяются из соотношений (5), (6), а упругие деформации ё из нелинейных соотношений:
11 ^ 1
г': =1---\ъ0вц +—аа;
1 1 3К 20) 0 1 20 1
аа = (а: +02 +Стэ)/3.
(9)
Здесь К — объемный модуль; О — модуль сдвига, которые зависят от уровня деформированного состояния. Зависимость (9) принимается в виде дробно-линейного соотношения, предложенного С.П. Тимошенко:
ао = 3К (ео) £(,; т = С(у; )у;,
(10)
где коэффициенты а и к зависят от сцепления С и угла внутреннего трения ф;
/2 = 6[(ах - ау )2 + (ау - аг)2 + (аг х)2 ] +
где
К =
Е0о,
; О = -
Оп т.
3(1- 2v)as + Е0 е0' т$ + О0 у ¡'
Прочность полиэтилена оценим по критерию Мизеса [6]
Ползучесть грунта учтем по теории течения с упрочнением, где тензор скорости деформаций ползучести имеет вид [7]
ст5, х5 — экспериментальные механические характеристики; Е0 — начальный модуль сжатия; О0 — начальный модуль сдвига.
При моделировании процессов деформирования грунта необходимо учитывать и дилатансию, т. е. изменение объема грунта от сдвига:
в^ = • (11)
1 + й2 у;
Здесь й1, й2 — экспериментальные механические характеристики; Y¡ — интенсивность сдвига. Тогда объемную деформацию грунта можно представить в виде суммы:
ее = е0 (12)
где 6° — объемная деформация, обусловленная все-
сторонним давлением; Ео — объемная дилатансионная часть деформации, вызванная сдвигом.
Прочность полиэтилена значительно ниже, чем у стали, а коэффициент линейного температурного расширения, наоборот, в 15-20 раз больше [8-9]. Поэтому учтем изменение прочностных и жесткостных характеристик в зависимости от температуры через их линейную аппроксимацию:
wnp04H
(Т) = ДТСТдт + ^прочнС^О )>
Одг = [аПрочн(7о)- СГпрочн №)]/ AT; -^прочн(Т) - АТЕДГ + £Прочн(То); -Едг = [-Епрочн(То ) - £прочн№)] / AT.
(13)
Сама задача определения НДС системы грунт — тонкостенная конструкция решается методом последовательных догружений, где на всех шагах нагружения решается вытекающее из (1) уравнение равновесия в приращениях [9]. Согласно (7), приращения напряжений
{Да} = [D]{Ase}.
(14)
Здесь [D] — матрица упругих характеристик; {Две} — вектор приращений упругих деформаций,
{Две} = {Дв} - {Дв^} - {Двс} - {Asd} - {Дв7}, (15)
где {Дв^}, {Двс} и {Дв^} определяются из соотношений (6), (5), (11), а {Двт} — из соотношения (7). Тогда разрешающее уравнение расчета по упруго-вязко-пластической модели принимает вид [10, 11]
Ш [Ae]T [Dep №}dV = Ш {AQ}T [Su]dV +
V V
+ JJ [AP }T [bu]dS + JJJ {Aed }T [D][bu]dV +
S V
+ JJJ {Aec }T [D][bu]dV + JJJ {AeT }T [D] [bu]dV. (16)
VV
Здесь [Dep] — упруго-пластическая матрица.
Из уравнения (16) определяются приращения перемещений {Дм}, а затем полные перемещения, деформации и напряжения:
{uk+1} = {uk} + {Дм}; {sk+1} = {sk} + {Дв}; {ak+1} = {ak} + {Да},
где k — номер шага по времени. За критерий достижения предельного состояния принимается условие резкого увеличения приращений.
Рассмотрим задачу оценки прочности трубопровода из полиэтилена низкого давления (ПЭНД), закрепленного на неподвижных бетонных опорах и находящегося в грунте (рис. 1). Диаметр трубы D = 0,12 м, относительный диаметр D/t = 17 (t — толщина трубы); глубина залегания h = 1,5 м; перепад температур AT = 30 °С, внутреннее давление в трубопроводе qr = 10 кг/см2. Геометрические, механические и силовые характеристики системы взяты из работы [11].
Рис. 1. Расчетная схема для трубопровода, закрепленного на неподвижных бетонных опорах
Поверочный расчет на прочность полиэтиленового трубопровода проведен с учетом всех действующих на него сил. Для этого использовался разработанный авторами комплекс программ МКЭ на основе 8-узлового изопараметрического трехмерного конечного элемента [10, 11]. Конечно-элементная модель участка трубы приведена на рис. 2. НДС сечения грунта, полученное с учетом всех факторов представлено на рис. 3. Изо-полями показаны уровни напряженности (3). Напряжения вдоль оси г для выбранного участка системы представлены на рис. 4.
J X
1
1
— 0.9
— 0.8
— 0.7
— 0.6
— 0.5
hz 0.4
Рис. 2. Конечно-элементная модель участка трубы (а) и ее поперечное сечение с изополями равного уровня интенсивности напряжения (б)
Рис. 3. НДС сечения грунта
#б [б51] 2014 известия высших учебных заведений. машиностроение
23
Рис. 4. Изополя напряжения для участка системы грунт — тонкостенная конструкция
Влияние толщины трубы на осадку системы грунт — тонкостенная конструкция представлено на рис. 5, а, температуры на предельную нагрузку — на рис. 5, б.
Выводы
1. Предложенный подход моделирования процесса деформирования и потери несущей способности системы грунт — тонкостенная конструкция позволяет решать задачу расчета трубы в грунте в более точной постановке, учитывать влияние различных факторов и изменение механических характеристик системы с течением времени.
t, см
Рис. 5. Зависимость осадки системы грунт — тонкостенная конструкция от толщины трубы (а) и предельной нагрузки от температуры (б)
2. Расчетным путем установлено, что увеличение температуры трубы уменьшает несущую способность системы.
Литература
[1] Бережной Д.Б., Паймушин Б.н. Математическое моделирование этапов строительства сложных сооруже-
ний по трансформирующимся расчетным схемам. Наукоемкие технологии, 200S, т. б, № 8-9, с. 59-б4.
[2] Бережной Д.в., Оагдатуллин М.К., султанов Л.У Моделирование деформирования обделки тоннеля
метрополитена, расположенной в грунте сложной физической природы. Вестник Казанского технологического университета, 2013, т. 1б, № 9, с. 2S0-2SS.
[3] Баженов в.Г., Ломунов в.К., Осетров с.Л., Павленкова Е.в. Экспериментально-расчетный метод иссле-
дования больших упруго-пластических деформаций цилиндрических оболочек при растяжении до разрыва и построение диаграмм деформирования при неоднородном напряженно-деформированном состоянии. Прикладная механика и техническая физика, 2013, т. S4, № 1(317), с. 11б-124.
[4] султанов Л.У, Бережной Д.в. Математическое моделирование несущей способности грунтовых
насыпей. Вестник КГТУ им. А.Н. Туполева, 2013, № 1, с. 117-124.
[5] Гаврюшин с.с., Барышникова О.О, Борискин О.Ф. Численные методы в динамике и прочности
машин. Москва, Изд-во МГГУ им. н.Э. Баумана, 2012. 492 с.
[6] Abrate S. Criteria for yielding or failure of cellular materials. Journal of Sandwich Structures and Materials,
2008, vol. 10, pp. S-S1.
[7] Ивлев Д.Д. Механика пластических сред. т. 1: Теория идеальной пластичности. Москва, Физматлит,
2001. 448 с.
[8] Лавров И.Г. Уточнение математической модели напряженного состояния полиэтиленовых труб
для расчета при различных температурах. Фундаментальные исследования, 2007, № 1, с. 44-4S.
[9] "Öiomas H.R., Missoum H. "Öiree-dimensional coupled heat, moisture, and air transfer in a deformable unsatu-
rated soil. International Journal for Numerical Methods in Engineering, 1999, vol. 44, no. 7, pp. 919-943.
[10] Каюмов P.A., шакирзянов Ф.Р. Моделирование поведения и оценка несущей способности системы
тонкостенная конструкция — грунт с учетом ползучести и деградации грунта. Учен. зап. Казан. ун-та. Сер. Физ.-матем. науки, 2011, т. 1S3, № 4, с. б7-75.
[11] Каюмов P.A., шакирзянов Ф.Р., шевченко с.Ю. Оценка несущей способности системы конструк-
ция — грунт. Известия КГАСУ, 2012, № 4(22), с. 49б-501.
[12] ГОст 1б338-85. Полиэтилен низкого давления. Москва, Изд-во стандартов, 2008. 3S с.
References
[1] Berezhnoi D.V., Paimushin V.N. Matematicheskoe modelirovanie etapov stroitel'stva slozhnykh sooruzhenii po
transformiruiushchimsia raschetnym skhemam [Mathematical modeling of the stages of the construction of complex structures by transforming design schemes]. Naukoemkie tekhnologii [Science Intensive Technologies]. 2005, vol. 6, no. 8-9, pp. 59-64.
[2] Berezhnoi D.V., Sagdatullin M.K., Sultanov L.U. Modelirovanie deformirovaniia obdelki tonnelia
metropolitena, raspolozhennoi v grunte slozhnoi fizicheskoi prirody [Simulation of deformation of tunnel lining subway, located in the soil complex physical nature]. VestnikKazanskogo tekhnologicheskogo universiteta [Herald Kazan Technological University]. 2013, vol. 16, no. 9, pp. 250-255.
[3] Bazhenov V.G., Lomunov V.K., Osetrov S.L., Pavlenkova E.V. Eksperimental'no-raschetnyi metod issledovaniia
bol'shikh uprugo-plasticheskikh deformatsii tsilindricheskikh obolochek pri rastiazhenii do razryva i postroenie diagramm deformirovaniia pri neodnorodnom napriazhenno-deformirovannom sostoianii [Experimental and computational method of studying large elastoplastic deformations of cylindrical shells in tension to rupture and constructing strain diagrams for an inhomogeneous stress-strain state]. Prikladnaia mekhanika i tekhnicheskaia fizika [Journal of Applied Mechanics and Technical Physics]. 2013, vol. 54, no. 1(317), pp. 100-107.
[4] Sultanov L.U., Berezhnoi D.V. Matematicheskoe modelirovanie nesushchei sposobnosti gruntovykh nasypei
[Mathematical simulation of the dirt fills load-carrying capacity]. Vestnik KGTU im. A.N. Tupoleva [Vestnik Kazan State Technical University A.N. Tupolev]. 2013, no. 1, pp. 117-124.
[5] Gavriushin S.S., Baryshnikova O.O, Boriskin O.F. Chislennye metody v dinamike i prochnosti mashin
[Numerical methods in dynamics and strength of machines]. Moscow, Bauman Press, 2012. 492 p.
[6] Abrate S. Criteria for yielding or failure of cellular materials. Journal of Sandwich Structures and Materials,
2008, vol. 10, pp. 5-51.
[7] Ivlev D.D. Mekhanika plasticheskikh sred. Teoriia ideal'noi plastichnosti [Mechanics of plastic media. The theory
of perfect plasticity]. Moscow, Fizmatlit publ., vol. 1, 2001. 448 p.
[8] Lavrov I.G. Utochnenie matematicheskoi modeli napriazhennogo sostoianiia polietilenovykh trub
dlia rascheta pri razlichnykh temperaturakh [Refinement of the mathematical model of stress state of polyethylene pipes for the calculation at different temperatures]. Fundamental'nye issledovaniia [Fundamental research]. 2007, no. 1, pp. 44-45.
[9] Thomas H.R., Missoum H. Three-dimensional coupled heat, moisture, and air transfer in a deformable
unsaturated soil. International Journal for Numerical Methods in Engineering, 1999, vol. 44, no. 7, pp. 919-943.
[10] Kaiumov R.A., Shakirzianov F.R. Modelirovanie povedeniia i otsenka nesushchei sposobnosti sistemy
tonkostennaia konstruktsiia — grunt s uchetom polzuchesti i degradatsii grunta [Behavior modeling and evaluation of the carrying capacity of thin — walled structures, taking into account soil creep and soil degradation]. Uchenye zapiski Kazanskogo Universiteta. Ser. Fiziko-matematicheskie nauki [Proceedings of Kazan University. Physics and Mathematics Series]. 2011, vol. 153, no. 4, pp. 67-75.
[11] Kaiumov R.A., Shakirzianov F.R., Shevchenko S.Iu. Otsenka nesushchei sposobnosti sistemy konstruktsiia —
grunt [Evaluation of bearing capacity of the soil — structure]. Izvestiia Kazanskogo Gosudarstvennogo Arkhitekturno-Stroitel'nogo Universiteta [News of the Kazan State University of Architecture and Engineering]. 2012, no. 4(22), pp. 496-501.
[12] GOST16338-85. Polietilen nizkogo davleniia. [Standard 16338-85. Low-pressure polyethylene]. Moscow,
Standards publ., 2008. 35 p.
Информация об авторах
КАюмов Рашит Абдулхакович (Казань) — доктор физико-математических наук, профессор, зав. кафедрой «Сопротивление материалов и основы теории упругости». КГАСУ (420043, Казань, Российская Федерация, ул. Зеленая, д. 1, e-mail: [email protected]).
ШАКиРзянов Фарид Рашитович (Казань) — кандидат физико-математических наук, ассистент кафедры «Сопротивление материалов и основы теории упругости». КГАСУ (420043, Казань, Российская Федерация, ул. Зеленая, д. 1, e-mail: [email protected]).
ГАвРюШин Сергей Сергеевич (Москва) — доктор технических наук, профессор, зав. кафедрой «Компьютерные системы автоматизации производства». МГТУ им. Н.Э. Баумана (105005, Москва, Российская Федерация, 2-я Бауманская ул., д. 5, стр. 1, e-mail: [email protected]).
Статья поступила в редакцию 01.04.2014
Information about the authors
KAYuMOV Rashit Abdulkhakovich (Kazan) — Dr. Sc. (Phys. Math.), Professor, Head of «Strength of Materials and Fundamentals of the Theory of Elasticity» Department. Kazan State University of Architecture and Engineering (KSUAE, Zelenaya str., 1, 420043, Kazan, Russian Federation, e-mail: [email protected]). SHAKIRZYANOV Farid Rashitovich (Kazan) — Cand. Sc. (Phys. Math.), Assistant of «Strength of Materials and Fundamentals of the Theory of Elasticity» Department. Kazan State University of Architecture and Engineering (KSUAE, Zelenaya str., 1, 420043, Kazan, Russian Federation, e-mail: [email protected]). GAVRYuSHIN Sergey Sergeevich (Moscow) — Dr. Sc. (Eng.), Professor, Head of «Computer Systems of Automated Production» Department. Bauman Moscow State Technical University (BMSTU, building 1, 2-nd Baumanskaya str., 5, 105005, Moscow, Russian Federation, e-mail: [email protected]).
