Моделирование напряженно-деформированного состояния вантового пролета конструкций, служащих опорой линии газопровода
Л.Н. Панасюк, В. С. Тюрина, К.А. Пастарнак Донской Государственный Технический Университет
Аннотация: В статье представлены описание расчетной модели и некоторые результаты моделирования напряженно-деформированного состояния вантовых конструкций с учетом динамического воздействия. Отражены некоторые аспекты моделирования общей потери устойчивости в геометрически нелинейной постановке.
Висячий однопролетный вантовый переход через реку Энжерук запроектирован с основным пролетом Ь=164м, начальная стрела провеса Г=15м [1-3]. С одной стороны несущий трос подвешен к пилону (выполняется по схеме «качающейся опоры»), закрепленному пятью оттяжками. С другой стороны - непосредственно растяжками прикреплен к неподвижной массивной опоре пятью оттяжками. Оттяжки несущего каната располагаются на расстоянии 27 и 41м от оси пилона и несущего каната [4, 5]. Полная схема одноцепного перехода представляется в виде 27+164+41м с 4-я несущими канатами диаметром 38мм. На пролетной части висячего перехода прокладывается стальная прямошовная труба газопровода диаметром 325х9мм. Пролетная часть висячего перехода крепится к несущим канатам с помощью подвесок переменной длины из круглой стали диаметром 24мм с шагом 5м. Для восприятия боковых ветровых и сейсмических нагрузок пролетная часть расчаливается 6-ю парами оттяжек с диаметром каната 19мм под углом к горизонтальной поверхности и к оси трубопровода. Также парой поперечных диагональных оттяжек диаметром 19мм крепится пилон первой опоры от ветра и сейсмических нагрузок. Под газовой трубой проложен монорельс (двутавр № 18), по которому может перемещаться монтажная тележка для обслуживания газопровода.
Расчетная схема - пространственная стержневая система [6, 7]. Работа пилонов, трубопровода, монорельса и основного каната моделировалась пространственными геометрически нелинейными (теория третьего порядка) стержневыми элементами. Несущие и ветровые оттяжки пилонов, а также вертикальные элементы подвеса трубы к несущему тросу моделировались стержневыми элементами типа «вант» с учетом односторонней работы (только на растяжение). При моделировании опор пренебрегаем осадкой фундаментов, устраиваемых в жестком скальном основании.
Высокая гибкость конструкции, а также зависимость общей жесткости от достигнутых деформаций обуславливает необходимость учета этапности возведения. Для реализации этого эффекта построение полной модели системы должно проводиться в несколько этапов.
На первом этапе рассматривается расчет троса на собственный вес. Приблизительно задается форма провеса и выполняется итерационный нелинейный расчет по уточнению провеса от собственного веса. В конце этапа определена геометрия несущего троса, деформация от собственного веса и усилия в тросе. На втором этапе корректируется геометрия троса и к его элементам добавляются усилия предварительного напряжения, определенные как усилия в тросе при его расчете на собственный вес. На третьем этапе строится полная расчетная схема сооружения. Этот этап можно считать окончательным в построении топологии модели (рис. 1-2).
Параметры расчетной модели - 1146 конечных элементов, 1034 узлов, порядок системы разрешающих уравнений 6204, нагружений 17, комбинаций нагружений 23.
Рис.1. - Пространственные виды конечно-элементной модели перехода Т.к. выполняется расчет в нелинейной постановке, то нельзя использовать алгоритм определения расчетных сочетаний усилий [8-10]. Для получения диапазона изменения усилий в элементах конструкции необходимо выполнить серию отдельных расчетов на разные комбинации нагружений.
Особый интерес представляет расчет на общую потерю устойчивости. Для гибких вантовых систем «стандартные», предусмотренные в программных комплексах методы определения критической нагрузки общей потери устойчивости часто не приводят к искомому результату. Более того, получаемые в стандартном расчете результаты достаточно абсурдны. Действительно, полученные значения критического параметра нагрузки -небольшие по модулю отрицательные числа. Т.е. по стандартному расчету потеря устойчивости должна происходить при противоположных направлениях нагрузки. Однако, это вполне закономерный результат стандартного расчета. Дело в том, что в алгоритме расчета на общую нелинейную потерю устойчивости на параметр нагрузки умножаются все нагрузки в комбинации, в том числе и собственный вес. Согласно алгоритму нелинейного расчета определяется минимальный параметр нагрузки, при котором система становится неустойчивой. Минимальным по модулю является наиболее близкий к нулю отрицательный параметр нагрузки. При этом вертикальные тяжи, связывающие трубу и несущий трос, начинают работать на сжатие. Т.к. в модели предусмотрена односторонняя работа тяжей только на растяжение, то при появлении сжимающих напряжений подвески выключаются из работы. Трос начинает работать изолированно от трубы. Его собственный вес троса становится отрицательным, основные напряжения становятся сжимающими, что и определяет общую потерю устойчивости при отрицательном значении параметра нагрузки (рис. 2).
Рис.2. - Примеры форм потери устойчивости при «стандартном» нелинейном
расчете
Для устранения этого эффекта достаточно относить параметр ко всем нагрузкам, за исключением собственного веса конструкции. Например, при расчете устойчивости на ветровую нагрузку нужно постепенно увеличивать только ветровую нагрузку. В результате получим, во сколько раз следует увеличить значение ветровой нагрузки, чтобы система потеряла устойчивость.
Критерием потери устойчивости является сингулярность функционала полной потенциальной энергии Лагранжа. При переходе через точку бифуркации вторая вариация функционала Лагранжа равна нулю, что приводит к вырожденной касательной матрицы жесткости. Практически получить точное значение параметра невозможно, поэтому в закритической
точке касательная матрица содержит отрицательные диагональные элементы, что и будет являться численным критерием перехода через точку потери устойчивости.
На рис. 3-7 показаны некоторые результаты моделирования работы вантового перехода.
Комвиндшии £
•JX. min- О.Омм так— О.Омм dV: min- -95.3мм так- 14.1мм dZ: miri- -437.7мм тах- 75.5мм
Рис.3. - Деформации конструкции - комбинация 2 (вес несущих конструкций и воздействие от подвесной конструкции для обслуживания перехода)
Рис.4. - Деформации конструкции - комбинация 21 (собственный вес и сейсмическое воздействие)
• Максимальный прогиб 405мм возникает при равномерном нагреве. Он обусловлен естественным удлинением несущего троса.
• Прогиб 400мм возникает при действии монтажной тележки.
• Прогиб 376мм возникает при сейсмическом воздействии по первому опасному направлению.
• Максимальные прогибы без учета температуры 400мм -воздействие ремонтной тележки в центре пролета и 376мм при действии сейсмической нагрузки.
• Относительно исходного монтажного состояния (прогиб под действием собственного веса 280мм) максимальные дополнительные прогибы от ветровой нагрузки 125.
• Максимальные горизонтальные смещения из плоскости конструкции 1530мм и 1455мм - при колебаниях конструкции от ветрового воздействия. При сейсмических воздействиях максимальные горизонтальные смещения 315мм.
• Максимальные смещения верха пилонов: пилон 1 - 18.3мм (сейсмика) и пилон 2 - 20мм (сейсмика). Отклонение по горизонтали центральной части 63мм - ветровое воздействие.
Рис.5. - Огибающие эпюры изгибающего момента в несущем тросе
Рис.6. - Огибающие эпюры крутящего момента в несущем тросе
Мах N=21.2515 кН (Elem N 370), Min N=0 кН (Elem М 339) _ Г
Mill................ I М III
Рис.7. - Огибающие усилия в элементах вертикальной подвески
References