Three-Dimensional Analysis of T-connections of Cylindrical Shells Considering Stages of Construction
- Authors: Kosytsyn S.B.1, Akulich V.Y.1, Osetinskii L.N.1
-
Affiliations:
- Russian University of Transport
- Issue: Vol 21, No 3 (2025)
- Pages: 181-190
- Section: Analytical and numerical methods of analysis of structures
- URL: https://journals.rudn.ru/structural-mechanics/article/view/46166
- DOI: https://doi.org/10.22363/1815-5235-2025-21-3-181-190
- EDN: https://elibrary.ru/SSMRFB
- ID: 46166
Cite item
Full Text
Abstract
A numerical analysis of the stress-strain state of intersecting cylindrical shells has been performed taking into account various types of nonlinearities (physical, contact, and structural). Intersecting cylindrical shells are considered as part of a large-scale three-dimensional “shell - soil” system. The study identified the most stressed areas of intersecting shells that require special attention during the modeling process. It should be noted that the particular complexity of these models lies in their high computational dimensionality since they include both the structure under consideration and its surrounding soil, which imposes additional requirements on the software packages used to solve such problems. Using modern methods of numerical analysis can significantly improve the quality of modeling and increase the accuracy of the results obtained. In particular, by considering the nonlinear properties of materials, it is possible to more accurately assess the actual behavior of shells under different external influences. In addition, the analysis shows that accounting for construction stages has a significant impact on the distribution of stresses in the intersecting shells. This underscores the need to consider the sequence of construction works, which in turn can contribute to increasing the overall reliability of the structure. Future prospects involve testing and validating the developed numerical analysis techniques on real “shell - soil” systems. The application of the obtained results to the design of load-bearing structures in real underground construction projects is also envisaged.
Full Text
1. Введение Проектирование и возведение уникальных зданий и сооружений в современной строительной практике требует комплексного подхода, включающего точные расчеты и всесторонние исследования поведения несущих конструкций под воздействием различных факторов. Это обусловлено как возрастающей сложностью архитектурных и инженерных решений, так и необходимостью обеспечения высокой надежности, безопасности и долговечности объектов. Важную роль в процессе проектирования играют численные методы, позволяющие моделировать поведение конструкций при различных типах нагрузок: статических, динамических, температурных, сейсмических и др. Важным аспектом численного моделирования является учет стадийности строительства, который позволяет более точно воспроизвести реальный процесс возведения конструкции. Одними из наиболее значимых и ответственных объектов строительства являются уникальные подземные сооружения, такие как станции метрополитена, подземные переходы и транспортные узлы, глубокие тоннели, подземные депо и галереи, а также инженерные сооружения для обеспечения жизнедеятельности метро, включая системы вентиляции, водоотведения и энергоснабжения [1-3]. В процессе строительства подземных сооружений и метрополитенов, помимо основных объектов, таких как перегонные тоннели и станции, важную роль играют вспомогательные сооружения, например тоннельные сбойки. Эти элементы обеспечивают функциональность основных сооружений, но при этом требуют особого внимания при проектировании и расчетах, так как создаваемые ими пересечения создают сложное напряженно-деформированное состояние, требующее детального анализа. Таким образом, необходимо развивать методики численного анализа подобных систем «оболочка - основание», которые могут быть успешно применены для повышения точности расчетов и обеспечения надежности и безопасности строительных конструкций. Универсальный программный комплекс ANSYS Mechanical [4-8], использованный в работе, предоставляет широкие возможности для моделирования статического и динамического напряженно-деформированного состояния сложных пространственных конструкций, включая здания, сооружения и комбинированные системы. Программный комплекс поддерживает учет различных типов нелинейностей, таких как физическая (нелинейные свойства материалов), контактная (взаимодействие между телами) и конструктивная (стадийность возведения и нагружения конструкций). Это позволяет максимально точно воспроизвести реальные условия работы конструкций и оценить их поведение под приложенной нагрузкой. Кроме того, программный комплекс поддерживает использование языка программирования APDL (ANSYS Parametric Design Language), что дает возможность расширять его функциональные возможности. 2. Метод расчета Ранее авторами выполнен численный анализ напряженно-деформированного состояния ортогонально пересекающихся цилиндрических оболочек, взаимодействующих с основанием [9]. В данной работе приведен вид тройникового соединения, где оболочки расположены относительно друг друга в плане под углом 15 градусов. Выбор угла обусловлен приближением к реальным размерам конструкций метрополитена (камеры съездов). Оболочки также разбиты на несколько частей для их поэтапного возведения. Далее приведены основные размеры моделируемых пересекающихся оболочек и окружающего основания, а также параметры моделей материалов для них. Основная цилиндрическая оболочка диаметром D1 = 5,50 м и толщиной t1 = 0,25 м; примыкающая оболочка диаметром D2 = 3,85 м и толщиной t2 = 0,20 м. Оболочка [11-13] задана линейно-упругой моделью материала: модуль упругости Esh = 30 000 МПа, коэффициент Пуассона μsh = 0,2, плотность ρsh = 2300 кг/м3. Основание задано моделью материала Мора - Кулона (O. Mohr, C.A. Coulomb) [14-15] со следующими параметрами: модуль деформации Esoil = 30 МПа, коэффициент поперечной деформации μsoil = 0,3, плотность ρsoil = 2000 кг/м3, сцепление Csoil = 10 кПа, угол внутреннего трения φsoil = 25°. По аналогии с предыдущей работой [10] узел пересекающихся цилиндрических оболочек размещен на глубине 30 м от дневной поверхности. Размеры грунтового массива выбраны из условия затухания напряженно-деформированного состояния грунта и приняты по 5 диаметров большой оболочки слева и справа от нее. На рис. 1 показан общий вид расчетной модели (рис. 1, а) и вид пересекающихся цилиндрических оболочек под углом 15 градусов в плане (рис. 1, б), которые состоят из четырех частей основной цилиндрической оболочки и четырех частей примыкающей оболочки. а б Рис. 1. Расчетная модель: а - общий вид; б - вид тройникового соединения И с т о ч н и к: выполнено В.Ю. Акуличем в программном комплексе ANSYS Mechanical Figure 1. Calculation model: а - general view; б - view of the T-connection S o u r c e: made by V.Yu. Akulich in ANSYS Mechanical software package На рис. 2 показано сравнение места стыка рассматриваемых тройниковых соединений. Сверху - (рис. 2, а) для ортогонально пересекающихся оболочек и снизу (рис. 2, б) для пересекающихся под углом 15 градусов в плане. Следует отметить, что чем меньше угол пересечения оболочек, тем большие размеры приобретает место их стыка, что может привести к сильному увеличению общих габаритов модели и росту необходимых вычислительных ресурсов. а б Рис. 2. Место стыка рассматриваемых тройниковых соединений: а - ортогонально пересекающихся оболочек; б - оболочек, пересекающихся под углом 15 градусов в плане И с т о ч н и к: выполнено В.Ю. Акуличем в программном комплексе ANSYS Mechanical Figure 2. The joint of the considered T-connections: a - orthogonally intersecting shells; б - shells intersecting at an angle of 15 degrees in plan S o u r c e: made by V.Yu. Akulich in ANSYS Mechanical software package На нижней и боковых поверхностях грунтового массива, а также на краях тройникового соединения заданы граничные условия, обеспечивающие геометрическую неизменяемость системы. Эти условия включают фиксацию перемещений в соответствующих направлениях, что исключает возможность смещения конструкций и грунтового основания за пределы расчетной области и позволяет воспроизвести их реальную работу. Нагрузка на систему приложена только от собственного веса конструкций и грунтового массива [11-13]. Моделирование контактного взаимодействия реализовано с помощью контактных пар, расположенных на внешней поверхности тройникового соединения и границе грунтового основания. Зона контакта определяется в ходе решения задачи, поскольку ее точное положение заранее неизвестно и зависит от совокупности факторов, включая приложенные нагрузки, характеристики материалов, граничные условия и другие параметры [14-18]. В процессе расчета взаимодействующие поверхности могут как вступать в контакт, так и разъединяться, что требует учета нелинейного характера их взаимодействия. В целях анализа влияния стадийности строительства на напряженно-деформированное состояние конструкции составлены расчетные случаи с различным количеством этапов возведения тройникового соединения: 1, 2, 4 и 8 стадий. Вначале идет активация основной оболочки, затем примыкающей оболочки со стороны основной. В каждом случае дополнительно предусмотрена нулевая стадия, которая служит для определения бытового состояния грунтового основания до начала строительных работ. Также выполнен расчет без учета стадийности строительства для сравнения результатов. На рис. 3 показан вид одного из промежуточных этапов возведения пересекающихся цилиндрических оболочек, когда основная оболочка D1 = 5,50 м и толщиной t1 = 0,25 м уже установлена (1-4-я части), а примыкающая D2 = 3,85 м и толщиной t2 = 0,20 м установлена частично (5-я часть). На местах установки следующих частей (6-8-я части) примыкающей оболочки на данном этапе активированы только элементы основания грунта. Рис. 3. Вид одного из промежуточных этапов возведения пересекающихся цилиндрических оболочек И с т о ч н и к: выполнено В.Ю. Акуличем в программном комплексе ANSYS Mechanical Figure 3. The view of the intermediate stage of construction of intersecting cylindrical shells S o u r c e: made by V.Yu. Akulich in ANSYS Mechanical software package 3. Результаты и обсуждение По результатам выполненных расчетных случаев проведен сравнительный анализ максимальных эквивалентных напряжений по Мизесу [19-21] с результатами, полученными ранее для ортогонально пересекающихся цилиндрических оболочек. Кривые изменения напряжений в оболочке в зависимости от количества стадий в расчетном случае показаны на рис. 4. Для сравнения выбран расчетный случай с 8 этапами возведения оболочек, пересекающихся под углом 15 градусов в плане (зеленая штриховая линия на графике). Также показано максимальное эквивалентное напряжение по Мизесу в тройниковом соединении без учета стадий. Рис. 4. Максимальные эквивалентные напряжения по Мизесу в пересекающихся оболочках И с т о ч н и к: выполнено В.Ю. Акуличем в программном комплексе ANSYS Mechanical Figure 4. Maximum equivalent stresses (von Mises) in intersecting shells S o u r c e: made by V.Yu. Akulich in ANSYS Mechanical software package Напряжения по Мизесу σe определяются по формуле где - главные напряжения. Максимальные эквивалентные напряжения по Мизесу составляют 71,2 МПа и 158,6 МПа при 8 стадиях расчета в вариантах с ортогональным соединением и с соединением под углом 15 градусов в плане соответственно. Следует отметить, что максимальные напряжения возникают локально (на стыке) и в остальном теле оболочек напряжения значительно ниже. Распределение максимальных эквивалентных напряжений по Мизесу в теле ортогонально пересекающихся цилиндрических оболочек в случае с 8-ю стадиями возведения показано на рис. 5. Максимальные значения напряжений возникают в месте сопряжения оболочек. Распределение максимальных эквивалентных напряжений по Мизесу в теле оболочек в случае с соединением под углом 15 градусов в плане с 8-ю стадиями возведения показано на рис. 6. Максимальные значения напряжений возникают в месте сопряжения оболочек. В таблице представлены максимальные эквивалентные напряжения по Мизесу в теле тройникового соединения для рассмотренных расчетных случаев. а б Рис. 5. Распределение максимальных эквивалентных напряжений по Мизесу в теле ортогонально пересекающихся цилиндрических оболочек: a - место сопряжения оболочек; б - увеличенный вид места сопряжения оболочек И с т о ч н и к: выполнено В.Ю. Акуличем в программном комплексе ANSYS Mechanical Figure 5. Distributions of maximum equivalent stresses (von Mises) in the body of cylindrical shells: a - is the interface of the shells; б - enlarged view of the connection of shells S o u r c e: made by V.Yu. Akulich in ANSYS Mechanical software package Рис. 6. Распределение максимальных эквивалентных напряжений по Мизесу в теле оболочек в случае с соединением под углом 15 градусов в плане. Справа показано место сопряжения оболочек И с т о ч н и к: выполнено В.Ю. Акуличем в программном комплексе ANSYS Mechanical Figure 6. Distributions of maximum equivalent stresses (von Mises) in the body of the shells in the case of a connection at an angle of 15 degrees in plan. The location of the shell connection is shown on the right S o u r c e: made by V.Yu. Akulich in ANSYS Mechanical software package Максимальные эквивалентные напряжения по Мизесу в теле тройникового соединения № Количество стадий возведения Максимальные эквивалентные напряжения по Мизесу, МПа Ортогонально пересекающиеся цилиндрические оболочки 1 Без учета стадий 105,6 2 1 78,2 3 2 72,5 4 4 77,4 5 8 75,2 Цилиндрические оболочки с соединением под углом 15 градусов в плане 6 8 158,6 И с т о ч н и к: выполнено В.Ю. Акуличем в программном комплексе ANSYS Mechanical Maximum equivalent stresses (von Mises) in the T-connection № Number of construction stages Maximum equivalent stresses (von Mises), MPa Orthogonally intersecting cylindrical shells 1 Without stages 105.6 2 1 78.2 3 2 72.5 4 4 77.4 5 8 75.2 Cylindrical shells with a connection at an angle of 15 degrees in plan 6 8 158.6 S o u r c e: made by V.Yu. Akulich in ANSYS Mechanical software package 4. Заключение 1. Полученные результаты подтвердили важность учета стадийности возведения при определении напряженно-деформированного состояния тройниковых соединений, взаимодействующих с основанием. 2. Максимальные эквивалентные напряжения по Мизесу составляют 71,2 МПа и 158,6 МПа при 8 стадиях расчета в вариантах с ортогональным соединением и с соединением под углом 15 градусов соответственно. Следует отметить, что максимальные напряжения возникают локально (на стыке) и в остальном теле оболочек напряжения значительно ниже. 3. Перспективы дальнейших исследований связаны с применением нелинейных материалов оболочки и различными вариантами контактного взаимодействия оболочки и основания.About the authors
Sergey B. Kosytsyn
Russian University of Transport
Email: kositsyn-s@yandex.ru
ORCID iD: 0000-0002-3241-0683
SPIN-code: 9390-7610
Advisor of Russian Academy of Architecture and Construction Sciences, Doctor of Technical Sciences, Professor, Head of Department of Theoretical Mechanics
15, Obraztsova St, GSP-4, Moscow, 127994, Russian FederationVladimir Yu. Akulich
Russian University of Transport
Author for correspondence.
Email: vladimir.akulich@gmail.com
ORCID iD: 0000-0002-9467-5791
SPIN-code: 8428-4636
Candidate of Technical Sciences, Associate professor of Department of Theoretical Mechanics
15, Obraztsova St, GSP-4, Moscow, 127994, Russian FederationLeonid N. Osetinskii
Russian University of Transport
Email: leonid.osetinsckij@yandex.ru
ORCID iD: 0009-0003-8268-7482
SPIN-code: 8229-2270
student, technician at the “Heat and Mass Transfer in Construction” Scientific Research Center
15, Obraztsova St, GSP-4, Moscow, 127994, Russian FederationReferences
- Zolotov A.B., Akimov P.A., Sidorov V.N., Mozgaleva M.L. Numerical and analytical methods for calculating building structures. Moscow: ASB Publ.; 2009. (In Russ.) ISBN 978-5-93093-675-9
- Akimov P.A., Mozgaleva M.L. B-spline wavelet discrete-continual finite element method for the local solution to the two-dimensional problem of the theory of elasticity. Monthly Journal on Construction and Architecture. 2022;17(1):32–41. (In Russ.) https://doi.org/10.22227/1997-0935.2022.1.32-41 EDN IYPOKW
- Perelmuter A.V., Slivker V.I. Calculation models of structures and possibilities of their analysis. Moscow: ASV Publ.; 2020. (In Russ.) ISBN 9785930938067
- Gallager R. The finite element method. Fundamentals. Moscow: Mir Publ.; 1984. (In Russ.) Available from: https://dwg.ru/dnl/3236 (accessed: 02.03.2025).
- Trushin S.I. Finite element method. Theory and problems. Moscow: ACB Publ.; 2008. (In Russ.) ISBN 978-5-93093-539-4
- Basov K.A. ANSYS: user reference. Moscow: DMK Press, 2005. (In Russ.) ISBN 978-5-97060-593-6
- Theory Reference for the Mechanical APDL and Mechanical Applications. ANSYS, Inc. 2009.
- Bate K., Vilson E. Numerical methods of analysis and the finite element method. Moscow: Strojizdat Publ.; 1982. (In Russ.) Available from: https://dwg.ru/dnl/4071 (accessed: 02.03.2024).
- Mangushev R.A., Dyakonov I.P., Polunin V.M., Bashmakov I.B., Paskacheva D.A. Mathematical modeling of the operation of plate elements when working together with a soil base in conditions of flat deformation. Housing Construction. 2024;(11):37–46. (In Russ.) https://doi.org/10.31659/0044-4472-2024-11-37-46 EDN BYRBSC
- Kosytsyn S.B., Akulich V.Yu. Numerical stress analysis of orthogonally intersecting cylindrical shells interacting with soil considering stages of construction. Structural Mechanics of Engineering Constructions and Buildings. 2024; 20(4):303–310. https://doi.org/10.22363/1815-5235-2024-20-4-303-310 EDN TVXXYV
- Klochkov Yu.V., Dzhabrailov A.Sh., Ishchanov T.R., Marchenko S.S., Andreev A.S., Klochkov M.Yu. Finite element analysis of an elliptical cylinder in geometrically nonlinear formulation using vector form of interpolation procedure. PNRPU Mechanics Bulletin. 2022;(1):58–71. (In Russ.) https://doi.org/10.15593/perm.mech/2022.1.06 EDN MYVJBF
- Dzhabrailov A.Sh., Nikolaev A.P., Klochkov Yu.V., Gureeva N.A., Ishchanov T.R. Nonlinear deformation of axisymmetrically loaded rotation shell based on fem with different variants of definitional equations. Izvestiya of Saratov University. Mathematics. Mechanics. Informatics. 2022;22(1):48–61. (In Russ.) https://doi.org/10.18500/1816-9791-2022-22-1-48-61 EDN JHCOIF
- Treshchev A.A. Thin-walled cylindrical shells as an object of experimental determination of mechanical properties of structural materials (are classical representations always acceptable when testing tubular specimens). Structural Mechanics and Structures. 2020;(1):7–18. (In Russ.) EDN CTAVJZ
- Lalin V.V., Le T.K.Ch. Calculation of building structures for several dynamic actions with static consideration of higher vibration modes. Structural Mechanics of Engineering Constructions and Buildings. 2020;16(3):171–178. (In Russ.) https://doi.org/10.22363/1815-5235-2020-16-3-171-178 EDN VSEGWP
- Mozgaleva M.L., Akimov P.A. Localization of solution of the problem for Poisson’s equation with the use of B-spline discrete-continual finite element method. International Journal for Computational Civil and Structural Engineering. 2021;17(3):157–172. (In Russ.) https://doi.org/10.22337/2587-9618-2021-17-3-157-172 EDN MFCJOI
- Yankovsky A.P. Refined model of viscoelastoplastic deformation of reinforced cylindrical shells. PNRPU Mechanics Bulletin. 2020;(1):138–149. (In Russ.) https://doi.org/10.15593/perm.mech/2020.1.11 EDN GEUKCU
- Zenkevich O.K. The finite element method in engineering. Moscow: Mir Publ.; 1975. (In Russ.) Available from: https://djvu.online/file/DtUw9BqXrtZCc (accessed: 02.03.2024).
- Zveryaev E.M., Pyhtyn A.V., Hoa V.D. Spatial problem for rectangular elastic plate. Structural Mechanics and Analysis of Constructions. 2021;(4):2–11. (In Russ.) https://doi.org/10.37538/0039-2383.2021.4.2.11 EDN SUBZUM
- Bakulin V.N. Model for analysis of stress-strain state of three-layer cylindrical shells with rectangular cutouts. Mechanics of Solids. 2022;(1):122–132. (In Russ.) https://doi.org/10.31857/S0572329922010032 EDN UFTFXV
- Li J., Shi Z., Liu L. A scaled boundary finite element method for static and dynamic analyses of cylindrical shells. Engineering Analysis with Boundary Elements. 2019;98:217–231. https://doi.org/10.1016/j.enganabound.2018.10.024
- Zang Q., Liu J., Ye W., Yang F., Pang R., Lin G. High-performance bending and buckling analyses of cylindrical shells resting on elastic foundation using isogeometric scaled boundary finite element method. European Journal of Mechanics — A/Solids. 2023;100:105013. https://doi.org/10.1016/j.euromechsol.2023.105013
Supplementary files










