CC BY
9
№ 2 (95)
A, UNI
те)
UNIVERSUM:
технические науки
февраль, 2022 г.
ИНЖЕНЕРНАЯ ГЕОМЕТРИЯ И КОМПЬЮТЕРНАЯ ГРАФИКА
АВТОМАТИЧЕСКАЯ ЛИНЕАРИЗАЦИЯ ВЫПУКЛЫХ ГИПЕРПОВЕРХНОСТЕЙ И НЕСУЩАЯ СПОСОБНОСТЬ ОБОЛОЧЕК
Махмудов Максуд Шералиевич
докторант,
Бухарский инженерно-технологический институт, Республика Узбекистан, г. Бухара E-mail: shm. maxmudov@mail. ru
AUTOMATIC LINEARIZATION OF CONVEX HYPERSURFACES AND THE BEARING CAPACITY OF SHELLS
Makhmudov Maksud Sheralievich
Doctoral candidate, Bukhara Institute of Engineering and Technology, Republic of Uzbekistan, Bukhara
АННОТАЦИЯ
В данной работе решается оптимизационная задача: отыскивается минимум числа граней вписанного или описанного полиэдра при заданной точности решения задачи предельного равновесия. Основу работы составляет автоматическое построение выпуклых полиэдров вокруг поверхности, а также автоматическое формирование матрицы задачи линейного программирования.
ABSTRACT
In this paper, an optimization problem is solved: the minimum number of faces of an inscribed or described polyhedron is found for a given accuracy of solving the limit equilibrium problem. The basis of the work is the automatic construction of convex polyhedra around the surface, as well as the automatic formation of the matrix of the linear programming problem.
Ключевые слова: гиперсеть, метод конечных разностей, интерполяция, краевые условия, принадлежность, точка обвода, узлы, полный дифференциал, непрерывность, эрмитовый сплайн, пространство.
Keywords: Hyper net, finite difference method, interpolation, boundary conditions, affiliation, bypass points, nodes, complete, differential, continuity, Hermitian spline, spaces.
В настоящей статье рассмотрены теоретические результаты, применены к анализу прочности тонкостенных оболочек, используемых в строительстве, машиностроении и других отраслях техники в качестве несущих поверхностей. Среди оболочек покрытий преобладают железобетонные оболочки постоянной и переменной толщины. Современные нормы проектирования строительных конструкций из железобетона допускают неупругое поведение материала и развитие пластических деформаций. Последние возникают при одновременном в различных направлениях (в пологих оболочках обычно четыре - места независимых внутренних усилий).
Условие перехода материала оболочки в состояние имеет вид уравнения, связывающего все силовые факторы с некоторой константой материала. Такое условие называют условием пластичности (текучести), в пространстве внутренних усилий оболочки оно
имеет вид выпуклой замкнутой гиперповерхности второго порядка [6]. Эта поверхность используется в расчетах несущей способности оболочек. Эффективность того или иного метода расчета, в свою очередь, связана со способом линеаризации гиперповерхности текучести.
Механическое поведение и несущая способность оболочек. Теперь рассмотрим условия текучести. Статическая и кинематическая формулировка задач предельного равновесия оболочек показывает, что условие текучести играет важную роль в определении верхних и нижних оценок предельной нагрузки.
В общем случае в каждой точке оболочки действуют восемь внутренних усилий: нормальные силы Ых, Ыу, сдвигающие силы Ыху, поперечные силы Qx,Qy, изгибающие моменты Мх,Му, а также крутящие моменты Мху. Обычно влияние поперечных сил на переход материала в пластическое
Библиографическое описание: Махмудов М.Ш. АВТОМАТИЧЕСКАЯ ЛИНЕАРИЗАЦИЯ ВЫПУКЛЫХ ГИПЕРПОВЕРХНОСТЕЙ И НЕСУЩАЯ СПОСОБНОСТЬ ОБОЛОЧЕК // Universum: технические науки : электрон. научн. журн. 2022. 2(95). URL: https://7universum. com/ru/tech/archive/item/13145
№ 2 (95)
A UNI
/ш. те;
UNIVERSUM:
технические науки
февраль, 2022 г.
состояние невелико, поэтому часть из восьми перечисленных усилий имеет вид:
Р, (Ых, Му, Мху , Мх, Му, Мху) < К,, (1)
где К1 - константа, зависящая от предела текучести 60 материала при простом растяжении (сжатии).
Уравнение (1) в пространстве внутренних усилий Ых,Ыу,Ыху,..., Мху описывает замкнутую выпуклую гиперповерхность, окружающую начало координат. Существуют различные виды конкретных условий текучести (1), переменные Ых,Ыу,..., Мху входят в них, как правило, от второй степени. Вид зависит от свойств материала, от формы оболочки и других факторов. Существенной особенностью условий текучести (1) является то, что функция Р1 связывает все шесть силовых факторов [7].
Современная теория предельного равновесия оболочек позволяет исследователю выбрать условие текучести в готовом виде. Другая возможность состоит в том, чтобы получить это условие в пространстве внутренних усилий Ых,Ыу,..., Мху, исходя из условия текучести, сформулированного в выражениях <5Х ,<5у,<5ху. Переход от пространства напряжений к пространству усилий представляет собой достаточно трудную задачу, часто невыполнимую из-за невозможности проинтегрировать условие в напряжениях по толщине оболочки.
Рассмотрим вначале частный случай - безмо-ментную постановку задачи. При Мх = Му =
Мху = 0, получаем
dN
dN
дх + ду
'ху _
0;
dNy ~ду d2Z
+
dN
ху
дх
- 2 N
дх2™х 2 дхдуПхУ
= 0;
d2Z
(2)
где Ъ (х; у) - уравнение срединной поверхности в декартовых координатах;
q - интенсивность распределенной нагрузки, нормальной к срединной поверхности.
На область, занятую оболочкой, нанесем регулярную сетку с квадратной ячейкой, для узла с номером I; j получим из первого и второго уравнений (2), воспользовавшись односторонними разностями:
(Nx)i-1J - (NX)U + (Nxy). ._i - (Nxy).. = 0 (Ny)iJ-1 - (Ny) + (Nxy)._1. - (Nxy).. = 0
(3)
В третьем уравнении (1) коэффициенты при усилиях Ых,Ыу, Ыху приведены к срединной поверхности, постоянные либо изменяющиеся от узла к узлу. Важно отметить, что система, образованная третьим уравнением (2) и уравнениями (3), является алгебраической системой линейных уравнений.
Обозначив кривизны Кх = ^ ; Ку = ;
ду
d2z
Кху = —— и переходя к безразмерным усилиям, получаем:
пх = ЫхЫ-1ш, пу = ЫуЫ-1; пху = ^уЫ-1; Ыо = (ВокК
где К — тощина оболочки, получим систему линейных уравнений для узла сетки номером:
(Кх nx)itj - 2(Кхупху). . + (Ку пу). . + qt>j = 0 (nx)i-ij - (пх)и + (nXy)i j_1 - (nxy\j = 0 (4) (ny)ij-i - (ny) + (nxy) + (nxy) = 0 .
Уравнения вида (4) могут быть записаны для узла сеточной области.
В более общем случае, то есть в моментной теореме оболочек, для представления вторых производных моментов в третьем уравнении (4) воспользуемся вторыми разностями. Для точки с номером I; j получим:
д2 т,
дх2
д2 Шх дх2
дх2
= Ах 2((mx)i-ij - 2(mx)ij + (mx)i+u)1 = Ау-2((ту). - 2(ту). + (ту). )
(5)
= (4АхАу)-1(тХу).- .-1 + (тху).. - КД - (тху) ^
Здесь Ах, Ау - шаг сетки в каждом из направлений, величины тх,ту,тху — безразмерные изгибающие моменты, введенные:
тх = МХМ0 1; ту = МуМ0 1;
Lxy = 1V1XV1V10 ; ^0
mxv = МХУМ-1; Mo = 4-1&h2.
все алгебраические уравнения, полученные путем сеточной дискретизации, линейны относительно усилий пх,пу, ...,тху.
Условия закрепления кроев позволяют сразу указать значения внутренних усулий в точках сетки, совпадающих с контуром, для свободного опирания краев при:
Подстановка соотношений (5) в (2) позволяет третье уравнение заменить соответствующим алгебраическим выражением. Здесь важно заметить, что
X = ± a,Nx = 0,МХ = 0, 0 и т.п y = ± a,Ny = 0,МУ = 0, Nx Ф0 и т.п.
2
У
№ 2 (95)
A UNI
/ш. те;
UNIVERSUM:
технические науки
февраль, 2022 г.
Пусть для определенности срединная поверхность оболочки есть эллиптический параболоид:
Z=
2а2
( х 2 + у2),
(6)
где 2а - размер стороны оболочки в плане; $ - стрела подъема.
Вычисляя кривизну срединной поверхности, найдем:
„ d2z с
-2
ах2
ay2 J
(7)
^ g2z
= 0 .
Дополнительно обозначим y = f • а-1 -пологость, e = к • f-1 - относительную толщину, р = q в-1 - безразмерную интенсивность равномерной нагрузки. Примем для конкретного примера y = 0,2; e = 0,055.
В линейные алгебраические уравнения, полученные описанным способом, входит величина р — интенсивность распределенности нагрузки, а также внутренние имя nx,ny, ...,rnxy. В соответствии с теоремой о нижней границе предельной нагрузки получаем статическую формулировку задачи: необходимо, варьируя поля внутренних усилий пх, пу,...., rnxy, найти minp.
Так как все названные величины входят в уравнения равновесия в первой степени, для отыскания можно было бы применить эффективный метод линейного программирования, гарантирующий получение решения законченного числа шагов. Ниже представлен метод линеаризации выпуклых гиперповерхностей применительно к задачам несущей способности жесткопластических оболочек.
Конкретный вид условия пластичности зависит от свойств материала, от формы срединной поверхности и других особенностей задачи, среди свойств пластического материала, влияющих на вид и структуру, условия пластичности, особое значение имеет его «знак чувствительности», то есть зависимость
величины предела текучести от знака действующего усилия (растяжение и сжатие). В частном случае материал может одинаково сопротивляться сжатию и растяжению (сталь, алюминий, другие пластичные материалы), общему же случаю отвечает соотношение:
р = б- (6+)-1,
где б-, 6+- пределы текучести материала при сжатии и растяжении.
Для оболочек из равносопротивляющегося материала наибольшее распространение получило условие пластичности Мизеса:
N2 - Nxwy + щ + зм2у +
/х2
"X
+12Ь-2(М2 - Мху + М2 + 3М2у) < N2, (8)
являющееся гиперэллипсоидом в пространстве внутренних усилий.
Для того чтобы построить эффективную процедуру расчета несущей способности жесткопластиче-ских оболочек на основе линейно-программных методов, необходимо взамен криволинейных условий текучести (8) получить их линейные приближения в виде тел, ограниченных участками гиперплоскостей (полиэдрами). Проблема такой заме ног (проблема линеаризации) содержит следующие последовательности.
1. Методика линеаризации должна быть пригодна как для регулярных, так и для песочно-регу-лярных поверхностей.
2. С целью контроля погрешности при аппроксимации гиперповерхностей полиэдрами необходимо всякий раз строить как вписанные, так и описанные полиэдры.
3. Размер задачи, объем вычислений и расход компьютера заметно зависят от числа граней полиэдра и резко возрастают с увеличением их числа. При неполном количестве гиперграней полиэдра объем задачи невелик, но и точность результатов также невелика. Поэтому возникает вопрос о выборе оптимального числа граней, позволяющего получить заданную точность вычислений при наименьших затратах машиной времени.
Список литературы:
1. Ахмедов Ю.Х. Автоматическая аппроксимация односвязных гиперповерхностей полиэдрами применительно к расчетам несущей способности оболочек покрытий: Автореферат. - Киев, 1984. - 208 с.
2. Волков В.Я. Многомерная исчислительная геометрия : монография / В.Я. Волков, В.Ю. Юрков. - Омск : Изд-во ОмГПУ, 2008. - 244 с.
3. Гумен Н.С., Гумен B.C. Геометрическое моделирование некоторых многопараметрических систем химической технологии. - Киев : Вища школа, 1977. - 108 с.
4. Махмудов М.Ш. Формирование линеаризованной дискретной модели условия пластичности / Polish Science Journal. - Warsaw : Sp. z o. o, 2021. - Issue 1 (34). Part 2. - 112 p.
5. Махмудов М.Ш. Элементы гиперсетей и их взаимопринадлежность // Polish Science Journal. - Warsaw, 2020. -Issue 9 (30). - P. 70-73.
6. Ольшак В., Савчук А. Неупругие поведение оболочек. - М. : Мир, 1966. - 144 с.
7. Рассказов А.О., Дехтярь А.С. Предельное равновесие оболочек. - Киев : Вища школа, 1978. - 152 с.
• 7universum.com
UNIVERSUM:
, ТЕХНИЧЕСКИЕ НАУКИ_Февраль. 2022 г.
8. Филиппов П.В. Начертательная геометрия многомерного пространства и ее приложение. - Л. : Изд. ЛГУ, 1979. - 280 с.
9. Makhmudov M.Sh, Akhmedov Y.H. Use Of E4 Space // Describing A Graph-Analytical Representation Of Multi-Factor Events And Processes. - Vol. 09, Issue 09. - P. 194-197.
10. Makhmudov M.Sh. Construction Of A Hypernet With Use Finite Difference Method In E4 Space // Journal NX 6.11 A Multidisciplinary Peer Reviewed Journal. - 2020. - P. 238-239.
№ 2 (95)
