Научная статья на тему 'МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ ЦИЛИНДРИЧЕСКОЙ ОБОЛОЧКИ МЕМБРАННОГО ПОКРЫТИЯ С ПОДКРЕПЛЯЮЩИМ ЭЛЕМЕНТОМ'

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ ЦИЛИНДРИЧЕСКОЙ ОБОЛОЧКИ МЕМБРАННОГО ПОКРЫТИЯ С ПОДКРЕПЛЯЮЩИМ ЭЛЕМЕНТОМ Текст научной статьи по специальности «Строительство и архитектура»

CC BY
49
7
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
МАТЕМАТИЧЕСКАЯ МОДЕЛЬ / ДВУМЕРНАЯ АППРОКСИМАЦИЯ / ДВУМЕРНАЯ ИНТЕРПОЛЯЦИЯ / ГЕОМЕТРИЧЕСКИЙ ИНТЕРПОЛЯНТ / ЦИЛИНДРИЧЕСКАЯ ОБОЛОЧКА / НАПРЯЖѐННО-ДЕФОРМИРОВАННОЕ СОСТОЯНИЕ / ПОДКРЕПЛЯЮЩИЙ ЭЛЕМЕНТ

Аннотация научной статьи по строительству и архитектуре, автор научной работы — Конопацкий Евгений Викторович, Шпиньков Владимир Александрович, Бездитный Андрей Александрович

В работе на примере моделирования напряженно-деформированного состояния цилиндрических оболочек мембранных покрытий с подкрепляющим элементом рассмотрен подход к многомерной аппроксимации, прототипом которого послужил метод наименьших квадратов. Отличительной особенностью предложенного подхода является отсутствие необходимости составления и решения системы линейных алгебраических уравнений для определения полиномиальных коэффициентов аппроксимирующей функции. Вместо этого для минимизации суммы квадратичных отклонений между исходными данными и расчётными используются быстродействующие численные алгоритмы поиска экстремальных значений, реализованные в программном пакете MS Excel в виде функции «Поиск решения». В результате моделирования получаются уравнения, содержащие натуральные значения факторов. А в случае использования нелинейных координатных сеток или нормирования можно легко перейти к натуральным значениям факторов путём замены переменных. Из полученных результатов видно, что предложенный подход к аппроксимации многомерных экспериментальных данных является в достаточной степени гибким и эффективным инструментом, но вместе с тем обладает недостатками, присущими классическому методу наименьших квадратов в части возникновения незапланированных осцилляций между узловыми точками аппроксимации. Поэтому в работе приведен пример использования геометрической теории многомерной интерполяции для решения тех же задач моделирования, но с использованием геометрических интерполянтов. Как видно из результатов, в конкретном случае модели, полученные на основе геометрической теории многомерной интерполяции, наиболее точно отображают характер протекания процесса и потому являются более предпочтительными по отношению к моделям, полученным с помощью двумерной аппроксимации. Вместе с тем аппроксимационные модели получены в виде явных функций, а интерполяционные - в параметрическом виде.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

MATHEMATICAL MODELING OF CYLINDRICAL SHELL STRESS-DEFORMED STATE OF MEMBRANE COATING WITH A REINFORCING ELEMENT

This work considers an approach to multidimensional approximation with regard to modeling the stress-strain state of cylindrical shells of membrane coatings with a reinforcing element. The proposed approach is characterized by the absence of the need to compose and solve a system of linear algebraic equations, in order to determine the polynomial coefficients of the approximating function. Instead, in order to minimize the sum of squared deviations between the initial data and those calculated, high-speed numerical algorithms to define extreme values are used. They were obtained in the MS Excel software package in the form of the “Find solution” function. The proposed approach to the approximation of multidimensional experimental data is a flexible and effective tool. Nevertheless, it has certain disadvantages inherent in the classical least square method in terms of the occurrence of unplanned oscillations between the nodal points of the approximation. Therefore, this paper provides an example of the use of geometric theory of multidimensional interpolation, in order to solve the same modeling problems, but using geometric interpolants. As can be seen from the results, in a specific case, the models obtained on the basis of the geometric theory of multidimensional interpolation most accurately reflect the nature of the process. In this regard, they are more preferable in relation to models obtained using the two-dimensional approximation. Approximation models are obtained in the form of explicit functions, and interpolation models are obtained in parametric form.

Текст научной работы на тему «МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ ЦИЛИНДРИЧЕСКОЙ ОБОЛОЧКИ МЕМБРАННОГО ПОКРЫТИЯ С ПОДКРЕПЛЯЮЩИМ ЭЛЕМЕНТОМ»

Инженерная геометрия и компьютерная графика. Цифровая поддержка жизненного цикла изделий

УДК 004.925.8: 004.942 DOI: 10.14529/build220406

МАТЕМАТИЧЕСКОЕ МОДЕЛИРОВАНИЕ НАПРЯЖЕННО-ДЕФОРМИРОВАННОГО СОСТОЯНИЯ ЦИЛИНДРИЧЕСКОЙ ОБОЛОЧКИ МЕМБРАННОГО ПОКРЫТИЯ С ПОДКРЕПЛЯЮЩИМ ЭЛЕМЕНТОМ

Е.В. Конопацкий1, В.А. Шпиньков2, А.А. Бездитный3

1 Нижегородский государственный архитектурно-строительный университет, г. Нижний Новгород, Россия

2 Донбасская национальная академия строительства и архитектуры, г. Макеевка, Донецкая Народная Республика

3 Севастопольский филиал РЭУ им. Г.В. Плеханова, г. Севастополь, Россия

В работе на примере моделирования напряженно-деформированного состояния цилиндрических оболочек мембранных покрытий с подкрепляющим элементом рассмотрен подход к многомерной аппроксимации, прототипом которого послужил метод наименьших квадратов. Отличительной особенностью предложенного подхода является отсутствие необходимости составления и решения системы линейных алгебраических уравнений для определения полиномиальных коэффициентов аппроксимирующей функции. Вместо этого для минимизации суммы квадратичных отклонений между исходными данными и расчётными используются быстродействующие численные алгоритмы поиска экстремальных значений, реализованные в программном пакете MS Excel в виде функции «Поиск решения». В результате моделирования получаются уравнения, содержащие натуральные значения факторов. А в случае использования нелинейных координатных сеток или нормирования можно легко перейти к натуральным значениям факторов путём замены переменных. Из полученных результатов видно, что предложенный подход к аппроксимации многомерных экспериментальных данных является в достаточной степени гибким и эффективным инструментом, но вместе с тем обладает недостатками, присущими классическому методу наименьших квадратов в части возникновения незапланированных осцилляций между узловыми точками аппроксимации. Поэтому в работе приведен пример использования геометрической теории многомерной интерполяции для решения тех же задач моделирования, но с использованием геометрических интерполян-тов. Как видно из результатов, в конкретном случае модели, полученные на основе геометрической теории многомерной интерполяции, наиболее точно отображают характер протекания процесса и потому являются более предпочтительными по отношению к моделям, полученным с помощью двумерной аппроксимации. Вместе с тем аппроксимаци-онные модели получены в виде явных функций, а интерполяционные - в параметрическом виде.

Ключевые слова: математическая модель, двумерная аппроксимация, двумерная интерполяция, геометрический интерполянт, цилиндрическая оболочка, напряженно-деформированное состояние, подкрепляющий элемент

Введение

На сегодняшний момент развитие строительства большепролетных мембранных покрытий остается актуальным в связи со значительным числом возводимых объектов спортивного, социально-культурного и общественного назначения [1-4]. Ввиду этого вопрос надежности проектируемых конструкций, а в особенности мембранно-

го покрытия, включающегося в совместную работу при поперечном изгибе с подкрепляющими элементами постели, остаётся актуальным. Проектируемые сегодня конструкции отличаются значительным разнообразием геометрии поверхности: плоские, положительной, отрицательной, нулевой гауссовой кривизны, первоначально плоские или с заданной стрелой провисания и т. п. Одной из осо-

бенностей работы таких конструкций является не в полной мере исследованный вопрос в части совместного взаимодействия тонкостенной мембранной оболочки и подкрепляющего элемента конструкции под действием поперечной нагрузки. Теоретические исследования и внедрение конструктивных решений для большепролетных мембранных конструкций изложены в работах [5-8].

С учётом исследований [9, 10] закономерности изменения напряжений, возникающие в тонколистовой мембране и подкрепляющем элементе жёсткости под действием поперечной нагрузки, определяются коэффициентами к, ¿2, ^ , которые характеризуют степень влияния подкрепляющего элемента на его совместную работу с мембранной оболочкой. В результате были получены экспериментальные данные, представленные в [11]. При этом возникает задача анализа, оптимизации и обобщения полученных экспериментальных данных для последующей разработки методики расчёта напряжённо-деформированного состояния тонкостенных оболочек с различной кривизной поверхности.

Постановка задачи. Разработать математические модели для коэффициентов к^ = / (х2),

характеризующих напряженно-деформированное состояние мембранных покрытий с различной кривизной поверхности, где х1 - функция жестко-стных характеристик, а Х2 - кривизна поверхности тонкостенной оболочки.

Метод

Базовыми инструментами математического моделирования для обработки и анализа экспериментальных данных являются интерполяция и аппроксимация. При этом наибольшее распространение в инженерной практике получил аппроксимаци-онный метод наименьших квадратов [12-14], который лежит в основе регрессионного анализа и заключается в минимизации суммы квадратичных отклонений между исходными данными и расчётными:

К( ^ - ^ )2 ^ mn-

где - исходные значения экспериментальных данных; - расчётные значения, полученные в результате аппроксимации.

В соответствии с методом наименьших квадратов минимизация целевой функции в виде сум-

ставлением и решением системы алгебраических уравнений. Однако, как следует из того же математического анализа, в случае с функцией многих переменных исследование на экстремумы (в нашем случае - минимум) необходимо проводить не только внутри, но и с учётом границ исследуемой области, что в значительной степени усложняет как математическое описание, так и программную реализацию такого решения. Исходя из этого в работе предлагается использовать быстродействующие численные алгоритмы поиска экстремальных значений, реализованные в программном пакете MS Excel в виде функции «Поиск решения», по аналогии с [15].

Традиционно критерием адекватности полученных моделей, свидетельствующих о степени их соответствия исходным экспериментальным данным, служит коэффициент детерминации R2, который определяется следующей зависимостью:

КZ - Z, )2

R2 = 1 - -

_\ 2

К(z, -z)

где г - среднеарифметическое значение исходных экспериментальных данных.

Коэффициент детерминации К2 для модели принимает значения от 0 до 1. Считается, что чем ближе значение К2 к 1, тем сильнее зависимость, а при К2 = 1 достигается функциональная зависимость между аппроксимирующей функцией и исходными данными.

Тогда получается, что максимальное значение коэффициента детерминации К2 на интервале от 0 до 1 достигается при минимальных значениях суммы квадратов регрессионных остатков

^ - ) . Исходя из вышеизложенного, пред-

1=1

лагается в качестве целевой функции использовать непосредственно коэффициент детерминации К2, но находить при этом его максимум с помощью функции «Поиск решения» на интервале значений от 0 до 1.

Для первых двух функций ¿1 = /1 (Х2) и

¿2 = /2 (х\, Х2) используются исходные данные

непосредственно из таблицы. Чтобы увеличить точность итоговой модели по отношению к исходным данным, применительно к факторам

Х1 и Х2 ,

мы квадратов регрессионных остатков zi - zi) определяющих функцию к3 = /3 (хьх2) , было

; I

использовано нормирование, т. е. каждое отдельное значение было поделено на максимальное. В результате получим следующие значения факторов для построения модели ¿3 = /3 (Х1, Х2)

происходит методами математического анализа. В общем случае задача минимизации решается методом математического анализа функции многих переменных, который подразумевает определение частных производных с последующим со-

(см. таблицу).

,=1

,=1

n

Нормированные значения факторов для построения модели k3 — f3 (x1,x2)

Х\ Х2

Натуральные значения Нормированные Натуральные значения Нормированные

фактора значения фактора фактора значения фактора

8,89 1,00 0 0

1,11 0,13 0,008 0,259

0,33 0,04 0,011 0,345

- - 0,016 0,515

- - 0,032 1

Обратное вычисление натуральных значений факторов осуществляется с помощью следующих формул:

Х = Зя , х2 , (1)

1 8,89 2 0,032

где хш и %2я - нормированные значения факторов Х1 и Х2.

В результате проведения вычислительных экспериментов были подобраны следующие аппроксимирующие функции, обеспечивающие высокие значения коэффициента детерминации:

k^ — aX\ + ^^Х^Х^ + С^Хх + "2X2 + ^2, б , г. 6 , 5 ,j5, 4 ,/-4

*2 -I- П. "V. "V- Ч- //- "V- С> -

(2)

7 _ б .7 б . 5 .7 5 i 4 _1_ /" 4

k3 = a3X2HXlH + °3 X2H + С3Х2 HX1H + "3 X2H + е3 X2HX1H + f3X2

+Яз X2HX1H + h3X23H + i3X2HX1H + j3X2H + k3 X1HX2H + l3X2H + да3'

Для реализации метода наименьших квадратов были использованы возможности табличного процессора MS Excel. При этом с помощью команды «Поиск решения» были найдены значения полиномиальных коэффициентов, обеспечивающие для каждой модели максимальные значения коэффициента детерминации на интервале от 0 до 1.

В результате получены следующие математические модели в виде полиномиальных зависимостей:

k1 — 2, 55x2X2 - 0,046x2 + 3,057x1x2 + 0,80 1x1 + 2, 925x2 + 0,789;

k2 — 0,058x1 -0,009x1x2 -0,633x1 + 0,02x2 +1,502; k3 — 258,033x6Hx1H + 546,743x6H -3,064x^Hx1H -0,754x^H --260,239x4Hx1H -543,407x^H -258,591x3Hx1H -551,504x3H + +340,309x2Hx1H + 686,824x2H - 74,744x1Hx2H -107,228x2H + 6,672.

Для всех трёх моделей коэффициенты детерминации варьируются в пределах R2 — 0,997...0,999.

" 2

Поскольку общая дисперсия z, - Z) является величиной постоянной, аналогичный результат

г—1

можно получить, если минимизировать непосредственно сумму квадратов регрессионных остатков

" 2

z - zt) . Однако в этом случае решение нужно искать на всём интервале значений, который для каж-

г—1

дой отдельной функции будет различным. Поэтому предпочтительней использовать коэффициент детерминации R2, значения которого всегда находятся в интервале от 0 до 1.

Результаты и обсуждения

Полученные в результате моделирования коэффициенты детерминации говорят об очень высокой достоверности математических моделей. Однако чтобы убедиться в этом наглядно, выполним графическую визуализацию полученных моделей (рис. 1).

Рис. 1. Графическая визуализация математических моделей, полученных с помощью двумерной аппроксимации

Визуализация первой модели подтверждает достоверность модели по отношению к исходным данным. А вторая и третья модели, как видно из рис. 1, несмотря на высокие значения коэффициента детерминации, склонны к незапланированным осцилляциям, которые никак нельзя обосновать, исходя из условий протекания исследуемого процесса. В частности, коэффициенты ¿2 и ¿3 не могут быть отрицательными. Исходя из этого, для второй и третьей модели необходим поиск других решений.

Эту же задачу можно решить с использованием геометрической теории многомерной интерполяции [16-19], которая предусматривает использование геометрических интерполянтов для моделирования многофакторных процессов и явлений. Под геометрическим интерполянтом понимается геометрический

объект, проходящий через наперёд заданные точки [20], координаты которых соответствуют исходным экспериментальным данным. Для построения геометрических интерполянтов необходимо разработать геометрическую схему и описать её аналитически с помощью дуг кривых, проходящих через наперёд заданные точки. Аналитическое описание геометрических интерполянтов выполняется с помощью математического аппарата «Точечное исчисление» [21-23].

Геометрическая схема моделирования напряженно-деформированного состояния цилиндрических оболочек мембранных покрытий с подкрепляющим элементом (рис. 2) представляет собой геометрический алгоритм построения поверхности отклика, проходящей через 15 наперёд заданных точек. Он включает определение опорных линий

Рис. 2. Геометрическая схема моделирования напряженно-деформированного состояния цилиндрических оболочек мембранных покрытий с помощью двумерной интерполяции

A1B1C1D1E1.

A2 B2C2 D2 E2.

A3B3C3D3E3 и обра-

зующей линии M-1M2Mз геометрического интер-полянта.

На основании геометрической схемы (см. рис. 2) разработан вычислительный алгоритм аналитического определения искомой поверхности отклика, который включает 3 однотипных точечных уравнения для определения каждой из направляющих линий:

л(-4 13-3 13-2 2 -3^ -3 64_2 2 16 _

М ; = Л,\ и--и и +— и и - ии + В, 16м и--и и +--ии

] ] I 3 3 I ] I 3 3

+C(-12м 3м + 40м 2м 2 - 12мм 3

)+d{

16 _3 64 _2 2 3 .

—м м--м м + 16мм +

3 3

т—г | _3 13 _2 2 13—3 4

+E, -м м +--м м--мм + м

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

3 3

где и - текущий параметр точечного уравнения, который изменяется от 0 до 1; и = 1- и - дополнение параметра и до 1; 1 <у< 3.

Далее следует определение образующей линии поверхности отклика: М = V (у - 0,1 5у) М1 + 8, 842УМ2 + V (у - 6, 692У ) М3,

где V - текущий параметр точечного уравнения, который изменяется от 0 до 1; V = 1 - V - дополнение параметра V до 1.

Обратим внимание, что разработанная геометрическая схема и полученный на её основе вычислительный алгоритм являются одинаковыми для всех трёх моделей к, к2 и к3, но с разными исходными данными. Это даёт возможность посредством покоординатного расчёта получить 3 отдельные модели, которые аналитически определяются с помощью 3 различных систем параметрических уравнений. Чтобы не перегружать статью, приведём систему параметрических уравнений только для первой модели:

X = 2,927v2 -11,492v + 8,894

x2 = -0,006м4 + 0,096м3 -0,113м2 + 0,056м

k1 =(16,278v2 -16,492v-1,067)м4 + (-15,863v2 -9,31v + 28,373)м3 + +(1,267v2 + 33,066v -37,213)м2 + ^3,582v2 -20,109v +17,547)м --0,231v2 - 2,799v + 4,06.

Выполним визуализацию геометрических моделей для сравнения интерполяционных моделей с аппроксимационными (рис. 3).

Рис. 3. Визуализация геометрических моделей, полученных с помощью двумерной интерполяции

Из рис. 3 видно, что в данном конкретном случае модели, полученные на основе геометрической теории многомерной интерполяции, наиболее точно отображают характер протекания процесса и потому являются более предпочтительными по отношению к моделям, полученным с помощью двумерной аппроксимации. Вместе с тем аппроксимационные модели получены в виде функций, заданных в явном виде, в то время как интерполяционные модели, полученные на основе точечных уравнений и вычислительных алгоритмов на их основе, представлены в параметрическом виде, что затрудняет их дальнейшее использование в инженерной практике.

Выводы

Описанный в работе подход по реализации метода наименьших квадратов с помощью поиска решений в табличном процессоре MS Excel является в достаточной степени гибким и эффективным инструментом. Его преимуществами являются простота реализации, возможность использования натуральных значений факторов и возможность обобщения на многомерное пространство, что расширяет существующий инструментарий моделирования многофакторных процессов и явлений. Но также он обладает недостатками, присущими классическому методу наименьших квадратов в части возникновения незапланированных осцил-ляций между узловыми точками аппроксимации. Поэтому в работе для конкретных моделей реализована двумерная интерполяция на основе геометрической теории многомерной интерполяции, которая показала более качественные результаты не только по отношению к исходным данным, но и по отношению к характеру протекания процесса, что позволило избежать незапланированных осцилляций, которые были ярко выражены на аппроксимационных моделях напряженно-деформированного состояния цилиндрических оболочек мембранных покрытий с подкрепляющим элементом.

Литература

1. Мущанов, В.Ф. Мембранные системы большепролетных покрытий с конструктивным предварительным напряжением / В.Ф. Мущанов,

A.Н. Оржеховский, Д.С. Коровкина // Вестник Донбасской национальной академии строительства и архитектуры. - 2019. - № 4(138). - С. 38-43.

2. Ларионов, И.В. Основные достоинства применения мембранных пространственных конструкций покрытия / И.В. Ларионов // Поколение будущего: взгляд молодых ученых - 2021: сборник научных статей 10-й Международной молодежной научной конференции, Курск, 11-12 ноября 2021 г. - Курск: Юго-Западный государственный университет, 2021. - С. 210-212.

3. Веселов, В.В. Анализ обрушения каркаса при демонтаже покрытия спортивно-концертного комплекса «Петербургский» / В.В. Веселов // Обследование зданий и сооружений: проблемы и пути их решения: материалы XI науч.-практ. конф., Санкт-Петербург, 18-19 ноября 2021 г. -СПб.: СПбПУ, 2021. - С. 5-13.

4. Исследование живучести мембранного покрытия большепролетного сооружения / Г.М. Кравченко, Е.В. Труфанова, Ю.М. Боженко-ва, Д.А. Суслопаров // Строительство и архитектура. - 2019. - Т. 7, № 2. - С. 10-14. DOI: 10.29039/article_5d4c0a3c3ff755.70793988.

5. Еремеев, П.Г. Пространственные тонколистовые металлические конструкции покрытий / П.Г. Еремеев. - М.: Изд-во Ассоц. строит. вузов, 2006. - 560 с.

6. Расчет и проектирование пространственных металлических конструкций / Е.В. Горохов,

B.Ф. Мущанов, Я.В. Назим, И.В. Роменский. - Макеевка: ДонНАСА, 2012. - 560 с.

7. Конструкции стационарных покрытий над трибунами стадионов / Е.В. Горохов, В.Ф. Мущанов, Р.И. Кинаш и др. - Макеевка, 2008. - 404 c.

8. Моисеенко, М.О. Алгоритм расчета пластин с начальным прогибом с учетом податливости контура / М.О. Моисеенко, О.Н. Попов // Научный альманах. - 2017. - № 12-2(38). - С. 43-50.

9. Мущанов, В. Ф. Учет совместной работы тонколистовой мембраны с подкрепляющими элементами стабилизирующей системы / В.Ф. Мущанов, В.А. Шпиньков // Металлические конструкции. - 2016. - Т. 22, № 2. - С. 79-89.

10. Мущанов, В.Ф. Уточнённая оценка совместной работы тонколистовой мембранной обшивки с подкрепляющим элементом / В.Ф. Муща-нов, В.А. Шпиньков // Металлические конструкции. - 2018. - Т. 24, № 3. - С. 133-141.

11. Особенности совместной работы подкрепляющего элемента с пологой цилиндрической безмоментной оболочкой при действии поперечной нагрузки / В.Ф. Мущанов, Е.В. Конопацкий, В.А. Шпиньков, А.А. Крысько // Металлические конструкции. - 2021. - Т. 27, № 2. - С. 97-118.

12. Голованчиков, А.Б. Аппроксимация экспериментальных данных методом наименьших квадратов и методом наименьших относительных квадратов / А.Б. Голованчиков, К.Д. Минь, Н.В. Шибитова // Энерго- и ресурсосбережение: промышленность и транспорт. - 2019. -№ 1(26). - С. 42-44.

13. Герасимов, С.А. Метод наименьших квадратов и метод площадей: что и когда лучше? / С.А. Герасимов // Учебная физика. - 2021. - № 3. -С. 20-25.

14. Конопацкий, Е.В. Геометрический смысл метода наименьших квадратов / Е.В. Конопацкий // Вестник компьютерных и информационных технологий. -2019. - № 9(183). - С. 11-18. DOI: 10.14489/vkit.2019.09.pp. 011-018.

15. Konopatskiy, E. V. Study of high-strength steel fiber concrete strength characteristics under the influence of elevated temperatures using mathematical modeling methods / E. V. Konopatskiy, S.N. Mashtaler, A.A. Bezditnyi // IOP Conf. Series: Materials Science and Engineering. - 687 (2019). - P. 022040. DOI: 10.1088/175 7-899X/68 7/2/022040.

16. Конопацкий, Е.В. Геометрическая теория многомерной интерполяции / Е.В. Конопацкий // Автоматизация и моделирование в проектировании и управлении, 2020. - № 1(7). - С. 9-16. DOI: 10.30987/2658-6436-2020-1-9-16.

17. Konopatskiy, E.V. Geometric modeling of multifactor processes and phenomena by the multidimensional parabolic interpolation method / E.V. Konopatskiy, A.A. Bezditnyi // Journal of Physics: Conference Series: XIII International Scientific and Technical Conference "Applied Mechanics and Systems Dynamics", Omsk, 05-07 November 2019. -Omsk: Institute of Physics Publishing, 2020. -P. 012063. DOI: 10.1088/1742-6596/1441/1/012063.

18. Конопацкий, Е.В. Принципы построения компьютерных моделей многофакторных процессов и явлений методом многомерной интерполяции / Е.В. Конопацкий // Программная инженерия: методы и технологии разработки информационно-вычислительных систем (ПИИВС-2018): сб. науч. тр. II Междунар. науч.-практ. конф., Донецк, 14-15 ноября 2018 года. - Донецк: ДонНТУ,

2018. - С. 309-318.

19. Конопацкий, Е.В. Подход к построению геометрических моделей многофакторных процессов и явлений многомерной интерполяции / Е.В. Конопацкий // Программная инженерия. -

2019. - Т. 10, № 2. - С. 77-86. DOI: 10.1758 7/prin. 10.77-86.

20. Геометрическое моделирование адаптивных алгебраических кривых, проходящих через наперёд заданные точки / Е.В. Конопацкий, И.В. Селезнев, О.А. Чернышева [и др.] // Вестник компьютерных и информационных технологий. -2021. - Т. 18. - № 9(207). - С. 26-34. DOI: 10.14489/vkit. 2021.09.pp. 026-034.

21. Балюба, И.Г. Конструктивная геометрия многообразий в точечном исчислении: дис. ... д-ра техн. наук: 05.01.01 / И.Г. Балюба. - Макеевка: МИСИ, 1995. - 227 с.

22. Балюба, И.Г. Точечное исчисление / И.Г. Балюба, Е.В. Конопацкий, А.И. Бумага. - Макеевка: Донбасская национальная академия строительства и архитектуры, 2020. - 244 с.

23. Балюба, И.Г. Точечное исчисление. Историческая справка и основополагающие определе-

ния / И.Г. Балюба, Е.В. Конопацкий // Физико-техническая информатика (СРТ2020): материалы 8-й Междунар. конф., Пущино, Московская обл., 9-13 ноября 2020 года. - Нижний Новгород: Автономная некоммерческая организация в области информационных технологий «Научно-исследовательский центр физико -технической информатики», 2020. - С. 321-327. DOI: 10.3098 7/conferencearticle_5fd755c0adb 1 d9. 27038265.

Конопацкий Евгений Викторович, доктор технических наук, профессор кафедры инженерной геометрии, компьютерной графики и автоматизированного проектирования, Нижегородский государственный архитектурно-строительный университет (Нижний Новгород), e.v.konopatskiy@mail.ru.

Шпиньков Владимир Александрович, ассистент кафедры теоретической и прикладной механики, Донбасская национальная академия строительства и архитектуры (Макеевка), licvova@gmail.com.

Бездитный Андрей Александрович, кандидат технических наук, доцент кафедры экономики, финансов и учета, Севастопольский филиал РЭУ им. Г.В. Плеханова (Севастополь), bezdytniy@gmail.com.

Поступила в редакцию 21 июня 2022 г.

DOI: 10.14529/build220406

MATHEMATICAL MODELING OF CYLINDRICAL SHELL STRESS-DEFORMED STATE OF MEMBRANE COATING WITH A REINFORCING ELEMENT

E.V. Konopatskiy, e.v.konopatskiy@mail.ru

V.A. Shpinkov2, licvova@gmail.com

A.A. Bezditnyi3, bezdytniy@gmail.com

1 Nizhny Novgorod State University of Architecture and Civil Engineering, Nizhny Novgorod, Russia

2 Donbas National Academy of Civil Engineering and Architecture, Makeyevka, Donetsk People's Republic

3 Sevastopol branch of the Plekhanov Russian University of Economics, Sevastopol, Russia

This work considers an approach to multidimensional approximation with regard to modeling the stress-strain state of cylindrical shells of membrane coatings with a reinforcing element. The proposed approach is characterized by the absence of the need to compose and solve a system of linear algebraic equations, in order to determine the polynomial coefficients of the approximating function. Instead, in order to minimize the sum of squared deviations between the initial data and those calculated, high-speed numerical algorithms to define extreme values are used. They were obtained in the MS Excel software package in the form of the "Find solution" function. The proposed approach to the approximation of multidimensional experimental data is a flexible and effective tool. Nevertheless, it has certain disadvantages inherent in the classical least square method in terms of the occurrence of unplanned oscillations between the nodal points of the approximation. Therefore, this paper provides an example of the use of geometric theory of multidimensional interpolation, in order to solve the same modeling problems, but using geometric interpolants. As can be seen from the results, in a specific case, the models obtained on the basis of the geometric theory of multidimensional interpolation most accurately reflect the nature of the process. In this regard, they are more preferable in relation to models obtained using the two-dimensional approximation. Approximation models are obtained in the form of explicit functions, and interpolation models are obtained in parametric form.

Keywords: mathematical model, two-dimensional approximation, two-dimensional interpolation, geometric interpolant, cylindrical shell, stress-strain state, a reinforcing element

References

1. Mushchanov V.F., Orzhekhovskiy A.N., Korovkina D.S. [Long-Span Membrane Coating Systems with Constructive Prestress]. Proceeding of the Donbas National Academy of Civil Engineering and Architecture, 2019, no. 4(138), pp. 38-43. (in Russ.)

2. Larionov I.V. [Main Advantages of Using Membrane Spatial Roof Structures]. Pokolenie budushchego: vzglyad molodykh uchenykh - 2021: sbornik nauchnykh statey 10-y Mezhdunarodnoy molodezhnoy nauchnoy kon-ferentsii [Generation of the Future: the View of Young Scientists - 2021: Collection of Scientific Articles of the 10th International Youth Scientific Conference] Kursk, 2021, pp. 210-212.

3. Veselov V.V. [Analysis of the Collapse of the Frame During the Dismantling of the Covering of the Sports and Concert Complex "Peterburgsky"] Obsledovanie zdaniy i sooruzheniy: problemy i puti ikh resheniya: Materia-ly XI nauchno-prakticheskoy konferentsii [Inspection of Buildings and Structures: Problems and Ways of Their Solution: Proceedings of the XI Scientific-Practical Conference] St. Petersburg, 2021, pp. 5-13. (in Russ.)

4. Kravchenko G.M., Trufanova E.V., Bozhenkova Yu.M., Susloparov D.A. [Study of a Durability of the Membrane Layer of a Long-Span Superstructure]. Construction and Architecture, 2019, vol. 7, no. 2, pp. 10-14. (in Russ.). DOI: 10.29039/article_5d4c0a3c3ff755.70793988.

5. Eremeev P.G. Prostranstvennye tonkolistovye metallicheskie konstruktsii pokrytiy [Spatial Thin-Sheet Metal Structures of Coatings]. Moscow, Association of Construction Universities Publ., 2006. 560 p.

6. Gorokhov E.V., Mushchanov V.F., Nazim Ya.V., Romenskiy I.V. Raschet i proektirovanie prostranstven-nykh metallicheskikh konstruktsiy [Calculation and Design of Spatial Metal Structures]. Makeevka, DonNACEA Publ., 2012. 560 p.

7. Gorokhov E.V., Mushchanov V.F., Kinash R.I., Shimanovskiy A.V., Lebedich I.N. Konstruktsii statsio-narnykh pokrytiy nad tribunami stadionov [Constructions of stationary coverings over stadium stands]. Makeevka, 2008. 404 p.

8. Moiseenko M.O., Popov O.N. [Algorithm for Calculating Plates with Initial Deflection Taking into Account the Compliance of the Contour] Science Almanac, 2017, no. 12-2(38), pp. 43-50. (in Russ.). DOI: 10.17117/na.2017.12.02.043

9. Mushchanov V.F., Shpinkov V.A. [Accounting of Joint Work of Thin-Sheet Membranes with Fastening of the Stabilizer System] Metal constructions, 2016, vol. 22, no. 2, pp. 79-89. (in Russ.)

10. Mushchanov V.F., Shpinkov V.A. [Updated Assessment of Joint Work of Thin-Sheet Membrane Coating with a Reinforcing Element] Metal constructions, 2018, vol. 24, no. 3, pp. 133-141. (in Russ.)

11. Mushchanov V.F., Konopatskiy E.V., Shpinkov V.A., Krys'ko A.A. [Features of the Joint Operation of the Reinforcing Element with a Flat Cylindrical Momentary Shell Under the Action of a Transverse Load] Metal constructions, 2021, vol. 27, no. 2, pp. 97-118. (in Russ.)

12. Golovanchikov A.B., Minh C.D., Shibitova N.V. [The Approximation of Experimental Data Using the Least Squares Method and the Least Relative Squares Method] Energo- i resursosberezhenie: promyshlennost' i transport [Energy and Resource Saving: Industry and Transport], 2019, no. 1(26), pp. 42-44. (in Russ.)

13. Gerasimov S.A. [Least-square method and area-method: what's better?]. Educational Physics, 2021, no. 3, pp. 20-25. (in Russ.)

14. Konopatskiy E.V. [Geometric meaning of least squares method]. Vestnik komp'yuternykh i informatsion-nykh tekhnologiy [Herald of Computer and Information Technologies], 2019, no. 9(183), pp. 11-18. (in Russ.). DOI: 10.14489/vkit.2019.09.pp.011-018.

15. Konopatskiy, E.V., Mashtaler S.N., Bezditnyi A.A. [Study of high-strength steel fiber concrete strength characteristics under the influence of elevated temperatures using mathematical modeling methods] IOP Conf. Series: Materials Science and Engineering, 2019, no. 687, 022040. DOI: 10.1088/1757-899X/687/2/022040.

16. Konopatskiy E.V. [Geometric Theory of Multidimensional Interpolation]. Avtomatizaciya i modelirovanie v proektirovanii i upravlenii [Automation and Modeling in Design and Management], 2020, no. 1(7), pp. 9-16. (in Russ.). DOI: 10.30987/2658-6436-2020-1-9-16.

17. Konopatskiy E.V., Bezditnyi A.A. [Geometric modeling of multifactor processes and phenomena by the multidimensional parabolic interpolation method] Journal of Physics: Conference Series: XIII International Scientific and Technical Conference "AppliedMechanics and Systems Dynamics", Omsk, Institute of Physics Publishing, 2020, 012063. DOI: 10.1088/1742-6596/1441/1/012063.

18. Konopatskiy E.V. [Principles the Construction of Computer Models of Multifactor Processes and Phenomena's by Multidimensional Interpolation] Programmnaya inzheneriya: metody i tekhnologii razrabotki informat-sionno-vychislitel'nykh sistem (PIIVS-2018): Sbornik nauchnykh trudov II Mezhdunarodnoy nauchno-prakticheskoy konferentsii [Software Engineering: Methods and Technologies of Information and Computer Systems Development: Collected Papers], Donetsk, DonNTU, 2018, pp. 309-318. (in Russ.)

19. Konopatskiy E.V. [Approach to the Construction of Geometric Models of Multifactor Processes and Phenomena by the Method of Multidimensional Interpolation] Programmnaya inzheneriya [Software Engineering], 2019, vol. 10, no. 2, pp. 77-86. (in Russ.). DOI 10.17587/prin. 10.77-86.

20. Konopatskiy E.V, Seleznev I.V., Chernysheva O.A., Lagunova M.V., Bezditnyy A.A. [Geometric Modeling of Adaptive Algebraic Curves Passing Through Predetermined Points] Vestnik komp'yuternykh i informat-sionnykh tekhnologiy [Herald of Computer and Information Technologies], 2021, vol. 18, no. 9(207), pp. 26-34. (in Russ.). DOI: 10.14489/vkit.2021.09.pp.026-034.

21. Balyuba I.G. Konstruktivnaya geometriya mnogoobraziy v tochechnom ischislenii: dis. dokt. tekhn. nauk [Constructive geometry of manifolds in point calculus. Doct. sci. diss.] Makeevka, 1995. 227 p.

22. Balyuba I.G., Konopatskiy E.V., Bumaga A.I. Tochechnoye ischisleniye [Point Calculus]. Makeyevka, Donbasskaya natsional'naya akademiya stroitel'stva i arkhitektury Publ., 2020. 244 p.

23. Balyuba I.G., Konopatskiy E.V. [Point Calculus. Historical Background and Basic Definitions]. Fiziko-tekhnicheskaya informatika (CPT2020): Materialy 8-y Mezhdunarodnoy konferentsii, Pushchino, Moskovskaya obl., 09-13 noyabrya 2020 goda. Nizhniy Novgorod: Avtonomnaya nekommercheskaya organizatsiya v oblasti informatsionnykh tekhnologiy "Nauchno-issledovatel'skiy tsentr fiziko-tekhnicheskoy informatiki" [Physics and Technology Informatics (CPT2020): Proceedings of the 8th International Conference], 2020, pp. 321-327. (in Russ.). DOI: 10.30987/conferencearticle_5fd755c0adb1d9.27038265.

Received 15 June 2022

ОБРАЗЕЦ ЦИТИРОВАНИЯ

Конопацкий, Е.В. Математическое моделирование напряженно-деформированного состояния цилиндрической оболочки мембранного покрытия с подкрепляющим элементом / Е.В. Конопацкий, В.А. Шпиньков, А.А. Бездитный // Вестник ЮУрГУ. Серия «Строительство и архитектура». - 2022. - Т. 22, № 4. - С. 57-65. DOI: 10.14529/build220406

FOR CITATION

Konopatskiy E.V., Shpinkov V.A., Bezditnyi A.A. Mathematical modeling of cylindrical shell stress-deformed state of membrane coating with a reinforcing element. Bulletin of the South Ural State University. Ser. Construction Engineering and Architecture. 2022, vol. 22, no. 4, pp. 5765. (in Russ.). DOI: 10.14529/build220406

i Надоели баннеры? Вы всегда можете отключить рекламу.