Russian Federation
Belgorod, Russian Federation
UDK 62 Инженерное дело. Техника в целом
Building art consists in the unity of the function of the structure and its form, embodied in the material. The bearing capacity of the material depends entirely on the chosen form. This is expressed in the rational correlation of load and support. At the level of determining the topology of a structure, the process of “enveloping” a force field in natural systems by matter is considered as an analogy. This phenomenon correlates with the technical support of the oblique bend, which is typical, in particular, for the run of the roof of the structure. A channel that is rational under conditions of direct bending loses its effectiveness with oblique bending. The best performance is found in the Z-profile with a vertical wall. The orientation of its shelves corresponds to the principle of material saturation of the areas adjacent to the external force field. But even this profile does not fully satisfy the most effective resistance to oblique bending. The introduction of an inclined wall makes it possible to bring the material of the shelves closer to the external force field. A comparison of the functioning of the mentioned profiles is given on numerical examples, united by the designated cross-sectional area. For a Z-profile with an inclined wall, formulas for geometric characteristics are derived. The angle of inclination of the wall is determined from the condition of transforming an oblique bend into a straight bend, that is, the coincidence of the trace of the force plane with the main axis of the beam section. In the framework of the above studies, we can talk about a decrease in stress by about 80%. It is also oblique bending at nonlinear physical law considererd.
cross-sectional shape of the beam, load-bearing capacity of the material of construction, oblique bend, Z-profile with an inclined wall, Prandtl’s diagram
Введение. Строительное искусство состоит в единстве функции конструкции и ее формы, воплощенной в материале. Несущая способность материала всецело зависит от выбранной формы. Это выражается в рациональном соотнесении нагрузки и опоры [1-5]. Когда строительный материал определен, несущая способность конструкции является функцией ее топологии, геометрии и параметров элементов [6-10]. Этот принцип проектирования известный итальянский зодчий П. Нерви назвал «сопротивляемостью по форме».
Во второй половине XIX века инженер К. Кульман и анатом Г.Х. фон Мейер совместно проводили серьезные исследования влияния механики окружающей среды на трабекулярные структуры. Они положили начало биомеханике несущих конструкций. На результаты их сотрудничества обратил внимание профессор В. Ру, впервые использовавший термин «функциональная адаптация» [11].
Тектоника биологических систем выглядит значительно сложнее, так как «конструкции» включаются в организм, выполняющий качественно иные функции. Форма возникает здесь в результате наложения различных подсистем. Тем не менее закон «борьбы элементов» в организме [11] устанавливает конструктивное сочетание – максимум работы при минимуме материала. Функциональное раздражение создает укрепление нагруженного органа путем возрастания поставки вещества. Отсутствие раздражения дает возможность доставить вещество в другие органы, где, напротив, наблюдается повышение раздражения. Так представляется «обволакивание» материей силового поля, создающее возможность органическим системам приспособляться к воздействию внешних факторов.
Долгие по времени и многократные по частоте внешние воздействия умеренной интенсивности вызывают как функциональную, так и морфологическую перестройку затронутых структур. В то же время рациональная экономичность организмов при снижении уровня внешних воздействий ведет к снижению функциональных возможностей и прочностных показателей в первую очередь опорно-двигательного аппарата.
Конструкционная бионика представляет природные системы в виде объемных моделей. В этом заключается не только средство установления законов структурообразования организмов, но и способ решения насущных практических задач.
Идеи конструкционной бионики могут способствовать совершенствованию теории структурного синтеза, исходящей из принципа стационарного действия. «Обволакивание» материей силового поля можно толковать как противостояние внешним воздействиям со стороны поля внутренних сил в материале, отождествляемых с потенциальной энергией деформации. Сосредоточение материала диктуется «вызовом» внешнего силового поля.
Рассмотренное противостояние нагрузки и конструкции можно рассматривать в рамках принципа «сопротивляемости по форме». Остановимся на общем случае изгиба бруса (косом изгибе) [12-14]. Помимо линейно-упругого деформирования, рассмотрим случай физической нелинейности.
Методика исследования. Вариационное исчисление – универсальный аппарат для творчества в механике. Вариационные принципы представляют собой ее теоретическую основу с энергетическим содержанием, устанавливают рациональные подходы к решению задач. В то же время они имеют большое практическое значение. При сформировавшихся функционалах можно вывести дифференциальные уравнения и естественные граничные условия при таких обстоятельствах, когда непосредственный подход затруднителен. Кроме того, прямые методы дают возможность решить задачу на вариационной основе, не прибегая к дифференциальным уравнениям.
Эти свойства присущи и вариационным принципам синтеза конструкций, сформированным на распространении функционального пространства на конфигурацию и модули упругости материала. При этом потенциальная энергия системы, находящейся в устойчивом равновесии, имеет абсолютный минимум по перемещениям, а функционал в точке стационарности – минимакс, а именно минимум по функциям перемещений максимумов по функциям конфигурации и (или) модулей упругости [15-18].
В работе [15] это демонстрируется на примере проектирования консоли с поперечным сечением в форме двутавра, нагруженной моментом в главной плоскости, при дополнительном условии ограничения прогиба.
Дополнительные условия (в форме уравнений связи) накладываются на переменные функционала вариационной задачи, в которой определяется ее условный экстремум. Следствием стационарности функционала являются уравнения, присущие базовому функционалу (Лагранжа или Кастильяно), уравнения связи, специфические уравнения структурообразования. Последние обуславливают надлежащий в традиционном проектировании критерий оптимальности. Важно заметить, что он приобретает здесь объективный характер, в то время, как в других случаях, кроме оговоренных в теореме Васютинского [19], эта существенная установка не лишена субъективности.
В упомянутом примере использовано свойство двойственности реализации задач на условный экстремум. Вместо исследования функционала с энергетическим содержанием рассматривается минимизация площади поперечного сечения.
В общем случае изгиба бруса такая постановка задачи усложняется в связи с двучленностью формулы напряжений и необходимостью итерационной процедуры расчета. Более простое решение вытекает с позиций «обволакивания» материей силового поля.
В качестве целевой функции избираем экстремальные напряжения в сечении бруса, входящие в условия прочности. Для упрощения постановки задачи останавливаемся на чистом изгибе, вызывающем исключительно нормальные напряжения.
Включенные в поле зрения профили рассматриваются в последовательности их появления в рамках технического прогресса. Для оценки их эффективности сохраняется внешнее силовое поле и объем используемого материала, что выражается в одинаковой площади поперечного сечения бруса. В каждом случае вызывает интерес возможность трансформации косого изгиба в прямой изгиб, что упрощает деформационное состояние элемента конструкции.
Основная часть. Под косым изгибом подразумевается общий случай изгиба, при котором плоскость действия момента не содержит главную ось инерции сечения. Он характерен, в частности, для прогона покрытия сооружения (рис.1). Здесь f-f – след силовой плоскости, z, y – главные центральные оси инерции сечения бруса.
Для коробчатого сечения решается изопериметрическая задача (известна площадь сечения при заданном контуре), в которой определяются стороны внутреннего прямоугольника, при которых достигается минимум абсолютной величины напряжений в контрольных точках сечения.
Исследование на экстремум функции упомянутых напряжений сводится к системе двух уравнений, приводимых к одному нелинейному уравнению. Свидетельством минимума функции напряжений служит положительная разность ее значений для обозначенного решения и сравниваемого варианта в оговоренной окрестности.
Рис. 1. Профили в качестве прогонов покрытия
В продолжение исследований рациональной формы сечений бруса при косом изгибе [20] остановим внимание на профилях, показанных на рис. 1. Заметим, что в отдельных случаях они могут состоять из перфорированных элементов с поперечными связями [21].
Первый тип профиля представлен на рис. 2. В численном примере приняты размеры: b=11,2 см; t=2,5 см; h=20 см; t1=3,3 см. Площадь сечения составляет 122 см2, моменты инерции имеют величины: Iz=6517 см4, Iy=2224 см4. Центр тяжести О находится на расстоянии 4,98 см от наружной грани стенки.
Рис. 2. Профиль типа швеллер
Представив, как обычно, косой изгиб бруса моментом М как одновременный изгиб бруса моментами Мz=Мcosβ и Мy=Мsinβ, вычислим модуль экстремального напряжения при β=300:
Силовая плоскость проходит через II и IV квадранты системы координат zOy. И если во II квадранте расположена большая часть верхней полки длиной 9,52 см, то в IV квадранте – лишь небольшая часть нижней полки длиной 1,68 см. Это обстоятельство приводит к развороту нижней полки вправо и переходу к Z-овому профилю (рис. 3).
Рис. 3. Z-овый профиль с вертикальной стенкой
Оставляем те же размеры и площадь сечения. Моменты инерции имеют величины: Iz=6517 см4, Iy=3589 см4, Dzy= - 3553 см4.
Для нахождения положения главной оси инерции, определяемого углом α0 по отношению к оси z, вычисляем
Значению 2,4 соответствует величина α0=33040/. Главные моменты инерции равны: Iu=9051 см4, Iv=1055см4. Модуль экстремального напряжения при β=300 равен 0,0013М, то есть меньше, чем у сечения швеллера в 2,6 раза. Это является следствием того, что силовая плоскость значительно приблизилась к главной плоскости бруса (угол между ними составляет 3040/).
Если швеллер в качестве прогона «обречен» на косой изгиб, то при Z-овом профиле может произойти трансформация этой сложной деформации в прямой изгиб. Условием для этого служит равенство β= α0.
Дальнейшее развитие идеи «обволакивания» материей внешнего силового поля приводит к рассмотрению Z-ового профиля с наклонной стенкой (рис.4). Это своего рода скачок от качественного преобразования формы сечения, о чем свидетельствуют рис. 2 и 3, к количественному изменению его параметров (за счет длины и толщины стенки и ее смещения). В этом проявляется стремление к достижению лучшего противостояния нагрузке внутреннего силового поля.
Рис. 4. Z-овый профиль с наклонной стенкой
Такого рода сечение можно представить в виде суммы трех параллелограммов. Момент инерции относительно оси z равен
Момент инерции относительно оси y вычислим, пренебрегая площадями
и пренебрегая площадями
В итоге
Центробежный момент инерции имеет выражение:
Оставляем те же размеры и площадь сечения, что и у профиля на рис. 3. Трансформация косого изгиба в прямой изгиб может произойти (β= α0=300) при некотором угле γ. Используя формулу (1), записываем уравнение
Введя в него величины моментов инерции по формулам (2) – (6) и β= 300 получаем уравнение с неизвестным γ (из-за громоздкости не приводится). Для решения трансцендентного уравнения использован графический метод, в результате чего получена величина γ=7030/.
Вычисляем моменты инерции по формулам (2) – (6): Iz=6667 см4, Iy=4282 см4, Dzy= - 2334 см4, затем главный момент инерции:
Модуль экстремального напряжения равен 0,00072М, что в 1,8 раза меньше, чем в случае Z-ового профиля с вертикальной стенкой.
Рассмотрим теперь чистый косой изгиб при упруго-пластических деформациях. Ему соответствуют следующие условия: 1) линейная деформация по направлению оси бруса (х) для любого волокна, взятого в плоскости y=const; 2) кривизна изгиба kx для всех волокон в каком-либо поперечном сечении принимается постоянной, хотя точные кривизны отдельных волокон неодинаковы; 3) используется идеальный упруго-пластический материал с диаграммой Прандтля; 4) внутренними усилиями являются изгибающие моменты Mz и My.
Выведем зависимость между моментами Mz и My и кривизнами kx1 (в плоскости xOy) и kx2 (в плоскости xOz). Используем выражение линейной деформации:
(9)
при записи выражений Mz и My:
, (10)
. (11)
где Es – секущий модуль.
Нейтральная ось (нулевая линия) проходит через центр тяжести сеченя (рис. 5) и имеет угловой коэффициент . В линейно-упругой области нормальное напряжение пропорционально расстоянию от точки до нулевой линии. Граница между линейно-упругой и пластической областями, определяется уравнением
. (12)
где ɛy – деформация, соответствующая пределу текучести . Возможны два случая расположения пластических областей.
Рис. 5. Поперечное сечение бруса: а – случай 1, б – случай 2 косого изгиба
Случай 1 (рис. 5, а). Начальная абсцисса пластической деформации определяется по формуле:
. (13)
Выражения моментов имеют вид:
, (14)
. (15)
В итоге имеем:
(16)
(17)
В частном случае при имеем
. (18)
Первая формула характеризует момент при изгибе в плоскости xOy (с учетом того, что ).
Случай 2 (рис. 5, б). Выражения моментов имеют вид:
, (19)
. (20)
В итоге имеем:
, (21)
. (22)
Нелинейный характер зависимостей (16), (17), (21), (22) нацеливает на использование обратного метода, предусматривающего построение номограмм и . Задаются величины кривизн и по формулам (16), (17), (21), (22) вычисляются моменты Mz и My. Построенные номограммы используются при расчетах напряженно-деформированного состояния бруса, силовая плоскость которого содержит момент М.
Существует и другой путь решения задачи. Задаются характеристики материала: предел текучести , модуль продольной упругости Е, размеры сечения b и h, положение силовой плоскости (угол ) и координаты точки на границе упругой и пластической областей . По формулам (16), (17), (21), (22) вычисляются моменты Mz и My, после чего составляется уравнение , в котором неизвестными являются кривизны . Второе уравнение составляется на основе условия пластичности в заданной точке на границе упругой и пластической областей. Решение системы двух уравнений дает величины и . Затем находим величину и строим кривые: и , которые используются в расчетах напряженно-деформированного состояния бруса..
Момент , отвечающий пределу пластического сопротивления при косом изгибе, вычисляется в предположении равенства по всей площади сечения.
Для случая 1 ( ) имеем:
. (23)
В частном случае при
. (24)
Для случая 2 ( ) имеем:
. (25)
В частном случае при
. (26)
Суждения, касающиеся имеющей место идеализации так называемого пластического шарнира, соотносятся со строгостью применения физического закона в виде диаграммы Прандтля.
Выводы
Идея перехода от косого изгиба к прямому, то есть от сложной к простой деформации, сопряжена с теорией структурного синтеза, что в данном случае выражается в поиске оптимальной конфигурации. Вариационная постановка задачи структурного синтеза сводится к итерационной процедуре расчета. Более приемлемой оказалась вытекающая из нее идея «обволакивания» материей силового поля, имеющая также обоснование с позиций биомеханики.
В качестве целевой функции избраны экстремальные напряжения, зависящие от конфигурации сечения. Целенаправленное сопоставление швеллера и Z-ового профиля выявило значительное преимущество последнего. Завершающий шаг от косого изгиба к прямому – переход в сечении от вертикальной к наклонной стенке. Эта задача решается при исходных размерах и площади сечения, полученных при оптимизации Z-ового профиля с вертикальной стенкой. Достигнуто уменьшение экстремального напряжения в 1,8 раза.
Решена физически нелинейная задача косого изгиба бруса из материала с диаграммой Прандтля. Предложен обратный метод, предусматривающий построение номограмм – зависимостей внутренних усилий от кривизн, которые затем используются для расчетов напряженно-деформированного состояния бруса
1. Majid K.I. Optimum design of structures. London: Newnes - Butterworths, 1979. 238 p.
2. Diaz A.R., Kikuchi N. Solutions to shape and topology eigenvalue optimization sing a homogenization method. Int. J. Numer. Methods Eng., 1992. No 35. Pp. 1487--1502.
3. Rozvany G. I. N., Zhou N., Sigmund O. Topology optimization in structural design. Advances in design optimization. London: Adeli, 1994. Pp. 240-299.
4. Yang R.J., Chahande A.I. Automotive applications of topology optimization. Structural Optimization, 1995. No 9. Pp. 245-249.
5. Cardoso E.L., Fonseca J.S.O Complexity control in the topology optimization of continuum structures. J. of the Bras. Soc. Of Mech. Sci&Eng. 2003. Vol. 25. No. 3. Pp. 293-301.
6. Stromberg L.L., Beghini A., Baker W.F., Paulino G.H. Topology optimization for braced frames: Combining continuum and beam. Column elements. Engineering Structures. 2012. No. 37. Pp. 106-124.
7. Kryzhevich G.B., Filatov A.R. Complex approach to topologic and geometric optimization of ship’s constructions [Kompleksnyj podkhod k topologicheskoj optimizacii sudovyh konstrukcij]. Memoirs of Krylov’s State Scientific Center. 2020. No 1 (391). Pp. 95-108. (rus)
8. Tamrazjan A.G., Aleksejcev A.V. Modern optimization methods of constructive solutions for supporting systems of buildings and structures [Sovremennye metody optimizacii konstruktivnyh reshenij dlya nesushchih sistem zdanij i sooruzhenij]. Bulletin of MSBU. 2020. Vol.15. No. 1. Pp. 12-30. (rus)
9. Mishchenko A.V. Optimization of lattice-nonhomogeneous pivotal construchions on the basis of energy criterion [Optimizacziya strukturno neodnorodnyh sterzhnevyh konstrukcij na osnove energeticheskogo kriteriya]. News of Institutions of Higher Educations. Building. 2021. No. 6. Pp. 20-31. (rus)
10. Locmanova V.S., Vladimirova E.A. Supporting systems optimization of frames on the basis of energy criterion [Optimizaciya opornyh ustrojstv ramnyh sistem na osnove energeticheskogo kriteriya].Traditions, Modern Problems and Perspectives of Building Development. Grodno. 2021. Pp. 141-145. (rus)
11. Roux W. Gesammelte Abhandlungen über Entwicklungsmechanik der Organismen. Bd. 1-2. Leipzig, 1895.
12. Frajman F.I. Optimal design of a thin-shelled beam with box section under oblique bend [Optimal’noe proektirovanie tonkostennoj balki korobchatogo secheniya pri kosom izgibe]. Applied Mechanics. 1971. No 11. Pp. 96-104. (rus)
13. Pochtman Y.M. At question on selection of optimal cross-sections form under oblique bend [K voprosu o vybore optimal`noj formy sechenij pri kosom izgibe]. News of Institutions of Higher Educations. Mechanical Engineering. 1972. No. 2. Pp. 9-12. (rus)
14. Sofronov Y.A. About optimal shape of beam cross-section under oblique bend [Ob optimal`noj forme poperechnogo secheniya balki pri kosom izgibe]. Memoirs of Kazan Aviation Institute. Issue 189. 1975. Pp. 36-42. (rus)
15. Yuriev A.G. Variational principles of structural mechanics [Variacionnye principi stroitel’noj mehaniki]. Belgorod: Publishing house of BelGTASM. 2002. 90 p. (rus)
16. Zinkova V.A. Topology optimization of metallic trusses [Optimizaciya topologii metallicheskih ferm]. Bulletin of BSTU named after V.G. Shukhov. 2015. No. 2. Pp. 37-40. (rus)
17. Zinkova V.A., Yuriev A.G., Peshkova E.V. Designing of tube trusses without gusset plate with joint connections. International Journal of Applied Engineering Research. 2015. No. 5. Vol. 10. Pp. 12391-12398.
18. Yuriev A.G., Smoljago N.A., Yakovlev O.A. Displacements in pivotal systems over elasticity limit [Peremeshcheniya v sterzhnevyh sistemah za predelom uprugosti]. Bulletin of BSTU named after V.G. Shukhov. 2022. No. 3. Pp. 25-31. (rus)
19. Wasiutynski Z. On the congruency of the forming according to the minimum potential energy with that according to the equal strength. Bulletin de L’Academie Polonaise des Sciences. Serie des Sciences Techniques. 1960. Vol. 8. No. 6. Pp. 259-268.
20. Yuriev A.G., Tolbatov A.A., Smoljago N.A., Yakovlev O.A. Rational cross-section of a beam under oblique bend [Raczional`noe sechenie brusa pri kosom izgibe]. Bulletin of BSTU named after V.G. Shukhov. 2017. No. 11. Pp. 60-63. (rus)
21. Romanjuk V.V., Vasilenko V.B. Construction deformationity of perforated purlin with Z-profile under oblique bend [Deformativnost` konstrukczii perforirovannogo progona Z-obraznogo profilya v usloviyah kosogo izgiba]. Bulletin of Byelorussia-Russia University. 2015. No 1 (16). Pp. 111-119. (rus)
22. Pisarenko G.S., Yakovlev A.P., Matveev V.V. Handbook with respect to strength of materials [Spravochnik po soprotivleniyu materialov]. Kiev: Scientific Idea, 1988. 736p. (rus)