CC BY
35
УДК 539.3: 534.1
С. П. Иванов
РАСЧЕТ ТОННЕЛЕЙ НА ОСНОВЕ ВАРИАЦИОННОГО МЕТОДА В. З. ВЛАСОВА
Эксплуатация транспортных систем связана с экологической безопасностью. Поэтому проектирование и расчет тоннелей является актуальной задачей. Разрабатывается методика расчета тоннелей с учетом нелинейной диаграммы деформирования материала. Получены основные уравнения и представлен пример расчета.
Введение. В связи с экологической безопасностью и перегруженностью наземного транспорта возникает острая необходимость проектирования подземных транспортных систем. Наиболее экологически безопасным является метрополитен. Тоннели метрополитена короткой длины можно представить как призматические системы, контактирующие с упругой средой. В случае, когда отдельные участки испытывают сжатие в продольном направлении, то эти участки тоннеля необходимо рассматривать как пространственные задачи. Если оболочка тоннеля достаточно протяженная, то ее можно рассматривать как плоскую задачу теории упругости.
Целью работы является разработка методики расчета оболочки тоннелей.
Интерпретация результатов. Расчетная схема предлагается в виде однослойной модели Власова-Леонтьева [1, 2]. Как указывает В.З.Власов, данная модель является более совершенной, чем винклеровская модель, вследствие того, что она способна «распределять» нагрузку. В основу разрабатываемой методики расчета закладывается модель нелинейно-упругого материала оболочечной системы. Экспериментальные данные таких материалов, как бетон, различные сплавы, композиты, показывают, что зависимость между интенсивностями напряжений ст и деформаций е1 можно принять в виде кубического полинома
Е ■в1 -Ех в], (1)
где Е - начальный модуль упругости, Е1 - постоянная, учитывающая степень нелинейности материала (принимаются на основании экспериментальных данных).
Учитываем гипотезы Кирхгофа-Лява. Полагаем, что направляющие тензоров напряжений и деформаций совпадают. Используем известные соотношения между деформациями и перемещениями:
е* = е- 2гх; е = е- ^; ^ = в„- , (2)
„„„ дп ду дп ду д1^ д1^ д1^ хде в — — • в — — • в —--\--• у =-• у =-• у =-
дх д, д, дх дх д, дхд,,
Величины и, V и ^ соответствуют компонентам вектора перемещений точки К срединной поверхности пластинчатой системы в продольном х, поперечном ^ и нормальном г направлениях (рис.1).
Рис. 1. Схема участка оболочки тоннеля в упругой среде
Двумерная задача сводится к одномерной. Перемещения выражаются по В.З.Власову [1] и принимаются в виде разложений:
Чх,я) = £ иг(хМ (,); Чхя) = (х)¥к(я);
г к
м>(х,,) = £ (х)/ (,), (г = 1,2,3,...,т;к,й = 1,2,3,...,п).
(3)
Выражения интенсивности деформаций е^ и объемной деформации в с учетом гипотез Кирхгофа-Лява и сжимаемости материала (о2 = 0, = 0, в82 = 0) можно записать следующим образом:
42
3
е = ■
' 2(1 + v)'
Къ-е,)2 + (е,)2 + (е, -ех)2 +-е2а
1 - 2v
в = ~-(ех +е,X
1 - V
а деформация
е= —
V
1 -V
(ех +е,).
Составим удельную энергию изменения объема и формы [3]
Ф = 1 Кв2 + 2}(1 + v) ■ йег, 2 3 0
где К = Е/[3(1 - 2v)] - модуль объемного сжатия.
Работа, отнесенная к единице площади поверхности оболочки, равна
5/2
А = ^Фйг,
-5/2
(4)
(5)
(6) (7)
(8)
где 5 - толщина элемента оболочки.
Следуя методу В.З. Власова [1], призматическую систему принимаем дискретно-континуальной. Выбор координатных функций м (,), ¥к(^), / (5) осуществляется по виду деформированного состояния системы. Искомые функции и (х), V (х) и (х) являются обобщенными перемещениями и подлежат определению из
й
решения задачи. В дальнейших записях функций переменные, указанные в скобках, опускаем.
Составляем полную энергию для системы
П = Я ^ + qx + qs + qz ^^ (9)
I к d
где дх - интенсивности внешних нагрузок, действующих на систему в продоль-
ном, поперечном и нормальном направлениях.
Из условий совместности деформаций в местах соединений пластин оболочки можно принять при ё = к
Wd = V. (10)
Углы между пластинами учитываются при выборе координатных функций. Определим минимум функционала (9), используя уравнения Эйлера-Лагранжа [4]:
d дЕ дЕ
= 0;
dx dü. dü
, ! ' (11)
d2 dF d dF dF n
--т-- +--Г--= 0,
dx2 dVk dx dVk d¥k
где F - подынтегральная функция (11), штрихи означают обычные производные от функций по переменной х.
Полагаем, что прогибы пластин, составляющих систему, совпадают с осадкой упругой среды. Развернув уравнения (11) и присоединив работу реактивных давлений
0осн. /^Л осн. ^ [- л
j , Qv упругой среды [1J соответственно в продольном и нормальном направлениях:
= £4U; QVT = 2£Р1Ук -£slV, (12)
i k k
получим уравнения равновесия.
Принимая, что деформированное состояние соответствует моменту потери устойчивости системы в упругой среде, в точках поперечного сечения системы все обобщенные поперечные нагрузки Qv выразим через параметр продольной силы Р и кривизну образующей срединной поверхности оболочки [1]
Qv = P£ ~hkV'k , (13)
k
приходим к системе обыкновенных нелинейных дифференциальных уравнений продольно-поперечного изгиба и устойчивости оболочки в упругой среде:
Yi £ aUi-X bjü -£ CjkVk= фj
i i k
P ~ _ 2y v _
- Yi £ ehkVkV + £ Ь ~hk + rhk +(mhk--—ahk Wk- £ shkVk + (14)
V TG 1 + v 1 -v т
+£ chüi+ Qv= фь,
где Qv =1 +
s
т~_7 1 o. — _ 1 o . — _ 1 o . _ e
bJi _ bJi + aji; rhk _ rhk Phk; shk _ Y\shk shk; 7\ _ '
Gji; hk hk ^ / nk' hk /1 hk ^ hk' / i 2 \ -
g g g (1 — v )
E л E H л
a0i = E° VjV.ds; Pl =7EH!4 fkfhds;
J 2(1 + v0)Hm¡ J 12(1 + v0) s
s O _ E0 f(V, — Vo)
s Ub
{Jfkfhds + ^(1 Vo) Jfkfhds}.
(15)
hk TT 2\ I!*/ ^h J./ Uh
Hm (1 — V0) S 6
В уравнениях и в коэффициентах приняты следующие обозначения: Е0, у0 - соответствуют модулю деформации и коэффициенту Пуассона упругой среды; выражения Фj и Фь в (15) учитывают физическую нелинейность материала оболочки [5].
Для решения данной задачи составлена программа на языке Фортран. Интегрирование нелинейных дифференциальных уравнений ведется численным методом Рунге-Кутта. Поиск недостающих краевых условий нелинейной задачи проводится итерационным методом типа Ньютона.
Для определения критической нагрузки применяется шаговый метод. За ведущий параметр принимается продольная нагрузка. Момент расхождения итерационного процесса принимается за критическое состояние оболочки.
В качестве примера рассчитывалась на сжатие с кручением оболочка тоннеля замкнутого поперечного сечения, контактирующая с упругой средой. Геометрические параметры оболочки следующие: 8/а=1/48 - отношение толщины стенки оболочки к поперечному размеру, 1/а=15 - отношение длины к наибольшему поперечному размеру, а/в=1,5 - отношение большей стороны поперечного сечения к меньшей. Отношение деформируемого слоя упругой среды к поперечному размеру системы равно Н/а=2,5, а отношение Е0^=0,001. Степень физической нелинейности принималась Е1/£^=105 . Сжимающие силы прикладываются к торцам оболочки (рис. 2). Края оболочки опираются на диафрагмы, которые считаются абсолютно гибкими из своей плоскости и абсолютно жесткими в ней.
(и/5)-10Ч
УЗ 5"
S0 45"
Q-ra"'^ „ & 1
& «О -yir т / / \ i i
/¿//^ | 1 1 1 1 i » 1 i 1 1
0 i Z 2,6 ¿9 3 3,2
a-io3
Рис. 2. Относительное перемещение угловой точки контура поперечного сечения а-а оболочки в зависимости от действия нагрузки: 1, 2 - по линейной теории соответственно без учета и с учетом упругой среды; 3 - с учетом физической нелинейности и упругой среды
Вывод. На основании проведенных исследований можно сделать заключение, что полученные уравнения позволяют рассчитывать на устойчивость оболочки тоннеля с учетом нелинейного деформирования материала.
Список литературы
1. Власов, В.З. Тонкостенные пространственные системы / В.З. Власов. - M.: Госстройиздат, 1958. -502 с.
2. Власов, В.З. Балки, плиты и оболочки на упругом основании / В.З. Власов, Н.Н. Леонтьев. - М.: ГИФМЛ, 1960. - 491 с.
3. Лукаш, П.А. Основы нелинейной строительной механики / П.А. Лукаш. - М.: Стройиздат, 1978. -204 с.
4. Смирнов, В.И. Курс высшей математики / В.И. Смирнов. - М. - Л.: ГИТЛ, 1957. - 627 с.
5. Иванов, С.П. Расчет нелинейных пластинчатых систем вариационным методом В.З. Власова / С.П. Иванов // Известия вузов. Строительство. - 2002. - №6. - С. 18 - 23.
Поступила в редакцию 10.08.07.
