CC BY
12
УДК 539.3
МЕТОД ПОСЛЕДОВАТЕЛЬНЫХ ПРИБЛИЖЕНИЙ В ЗАДАЧЕ РАЦИОНАЛЬНОГО РАСПРЕДЕЛЕНИЯ ПОРИСТОСТИ ПРИ ЧИСТОМ ИЗГИБЕ БРУСА ПРЯМОУГОЛЬНОГО СЕЧЕНИЯ
Шляхов СМ., Гаврилов Д.Ю.
Саратовский государственный технический университет имени Гагарина Ю.А., Саратов, Россия
Статья посвящена задаче рационального распределения пористости по сечению бруса при чистом изгибе. Решение такой задачи позволит обеспечить необходимую несущую способность при снижении материалоемкости конструкции. Целью исследования является подобрать рациональный закон распределения пористости по прямоугольному сечению бруса при технических ограничениях производства.
Ключевые слова: пористость, чистый изгиб, прямоугольник, рациональность, брус, нормальные напряжения, модуль Юнга.
DOI: 10.22281/2413 -9920-2017-03-01-46-51
Рассмотрим чистый изгиб бруса, выполненного из материала (стали) пористой структуры с переменной по сечению пористостью. Интерес к рассмотрению задачи с постановкой именно такого вопроса обуславливается предположением о теоретической возможности варьирования механических свойств нагружаемой конструкции исходя из возникающих в сечениях напряжений. Известно, что механические характеристики материала (модуль Юнга Е и предел текучести аГ) являются функциями пористости материала [1, 2]. Практические результаты показывают, что при повышении уровня пористости р значение модуля Юнга Е снижается, а при снижении уровня пористости - повышается. Аналогичная зависимость просматривается и при анализе пары - предел текучести аГ - пористость р.
Ограничимся случаем упругого деформирования, полагая, что максимальные напряжения в брусе сттах не превышают предела текучести материала (7Г. На основе экспериментальных данных для пористой стали, приведенных в табл. 1, зависимости Е(р) и стт (p) могут быть представлены полиномом.
Таблица 1
Экспериментальные зависимости значений нормальных напряжений и модуля Юнга _от пористости_
Пористость ат, МПа Пористость Е, МПа
0 200 0 210000
0,12 116 0,1 160000
0,21 95 0,2 110000
0,31 59 0,3 80000
0,37 43 0,4 50000
0,43 32 0,5 20000
Сглаживая заданную функцию методом наименьших квадратов (МНК) получим полиномы:
Е = а1 + а2 р + а3 р2, (1)
аГ = Ь + Ь2 р + Ь3 р2. (2)
Конечным результатом МНК будут являться значения соответствия аГ, Е пористости в границах от 0 до 0,4.
На рис.1, а,б приведены графики функций (1) и (2) соответственно, при значениях коэффициентов (МПа), приведенных в табл.2 в интервале 0 < р < 0,4 .
Таблица 2
Значения коэффициентов полиномов, полученных методом наименьших квадратов
Коэффициент Значение коэффициента Коэффициент Значение коэффициента
а1 209285,7143 196,0368975
а 2 -535000 ь2 -645,9644653
а 3 321428,5714 Ъ3 627,5555095
а)
Е МП а
250000
20МИ0
150000
ЮТОМ)
50000
О
б)
гг. МП а
V.
* ♦
— Л ■ А ' ♦
О 0,03 0,06 0,09 0,12 0,14 0,16 0,16 0,21 0,23 0,26 0,29 0,31 0,32 0,34 0,35 0,37 0,39 0,4 0,42 0,43
Рис. 1. Графики зависимости значений модуля Юнга и нормальных напряжений от
пористости: а - Е; б - оТ
Положим теперь, что пористость р переменна по высоте бруса.
Р = Р (у),0 < у <
Тогда Е = Е(у) т.е. модуль Юнга есть функция координаты у (рис.2).
0,03 0,06 0,06 0,1 0,13 0,16 0,16 О,г 0,23 0,26 0,29 0,3 0,32 0,35 0,37 0,4 0,42 0,45 0,47
ъ/
/2
Задача изгиба, таким образом, сводится к изгибу бруса, выполненного из неоднородного материла с переменной по сечению упругой характеристикой Е = Е (у) и переменным
пределом текучести аГ = аГ (у).
Обозначим изгибающий момент на брусе через М, а высоту сечения - h. При чистом изгибе бруса работает гипотеза плоских сечений, т.е. сечения после деформаций остаются плоскими и нормальными к оси бруса.
Соответственно закону Гука при изгибе определим закон изменения напряжений по высоте сечения:
а = Ее . (3)
Величина максимальных деформаций е постоянна для данного сечения, а значение Е = Е() зависит от расстояния от оси сечения yj. При этом имеем ограничение
а <°г(у). (4)
Величина изгибающего момента в сечении определится по формуле
М « 2 а]Ь]5у] = 2±Е]е]Ъ]5у] (5)
1=1 1=1
Целью исследования является подобрать такой закон распределения модуля Юнга Е и, следовательно, пористости р по сечению бруса, чтобы получить максимально возможный изгибающий момент при ограничениях на напряжение а < аГ и на пористость
ртт < р < ртах , (6)
задаваемых техническими возможностями производства.
Данная постановка является сложной задачей линейного программирования, для решения которой используем приближенный подход. Суть приближенного решения заключается в следующем.
Разбиваем все сечение бруса по высоте сечения на п элементов с наружными высотами у 1 а = 1, 2, ..., п) с шагом
^=А,
2п
где h - высота бруса (рис.2).
£тах
ЕП,Р„
« .Р:
Е-:рг
Е
Ж
Я
СМ
Рис.2. Схема разбиения поперечного сечения на элементы
Соответственно получим
У1 =^; у1 = Уу-1 + 8, 1 = 2, 3, ..., п. Максимальные напряжения в слоях будут равны
= Е1£1 - О, ] = 1 2 п (7)
Для решения задачи применим метод последовательных приближений.
Зададимся вначале нулевой пористостью и максимальной деформацией £тах . По формуле (7) найдем значение нормальных напряжений для каждого элемента поперечного сечения. По диаграмме (рис.1, б) находим предельное значение пористости, соответствующее полученному значению нормальных напряжений, полагая о1 = оТ. Найденное значение пористости следует согласовать с условием (6). Примем в нашем случае ртп = 0, ртах = 0,4, т.е. в случае выхода искомого значения р за пределы наложенных ограничений следует принять граничное его значение.
Определяем величину изгибающего момента по формуле
М=2 . (8)
1=1
По полученному значению пористости на каждом слое сечения по формуле (1) найдем Е = Е^ и повторяем расчет.
Будем циклично продолжать описанные выше действия до тех пор, пока не выполнится условие
М+1 - Мг
М+1
•100% < 1%.
В качестве примера примем И = 20 см, п = 20, £тах = 0,000952381 , Ь = 10 см. Получим
значение пористости для каждого элементарного слоя и определяем момент. Результаты всех приближений сведем в табл.3.
Таблица 3
Сводные значения вычисленных значений изгибающих моментов на каждом приближении
№ приближения Значение изгибающего момента, кНм М , - М Значение г+1 г -100,% М,+1 '
5,7205
1 4,6802 18,19
2 4.1285 11,79
3 3,8071 7,79
4 3,5668 6,31
5 3,4283 3,88
6 3,3410 2,54
7 3,2339 3,20
8 3,2201 0,43
При оценке полученных результатов, мы видим, что достигли необходимого результата на 8-м приближении.
Окончательно функции рационального распределения р(у) и о(у) отражены на эпюрах (рис. 3).
Таким образом, описанный метод решения является достаточно точным и объективным методом решения задачи. Варьируя значение £тах, получим функциональную зависимость Мизг (£тах, р), из которой найдем расчетное значение £тах, соответствующее требуемому изгибающему моменту.
Научно-технический вестник Брянского государственного университета, 2017, №1
Р, % а, МП а
о
Рис. 3. Эпюры распределения пористости и соответствующих значений нормальных напряжений по поперечному сечению бруса
Список литературы
1. Кашталян, Ю.А. Характеристики упругих материалов при высоких температурах / Ю.А. Кашталян. - Киев: Наукова думка, 1970. - 112 с.
2. Белов, С.В. Пористые металлы в машиностроении / С.В. Белов. - М: Машиностроение, 1981. - 247 с.
Об авторах
Шляхов Станислав Михайлович - доктор физико-математических наук, профессор кафедры «Теория сооружений и строительных конструкций» Саратовского государственного технического университета имени Гагарина Ю.А.
Гаврилов Данила Юрьевич - аспирант кафедры «Теория сооружений и строительных конструкций» Саратовского государственного технического университета имени Гагарина Ю.А., [email protected].
METHOD OF CONSECUTIVE APPROXIMATIONS IN THE PROBLEM OF RATIONAL DISTRIBUTION OF POROSITY AT THE PURE BEND OF THE BAR OF
RECTANGULAR SECTION
Shlyakhov S.M., Gavrilov D.Yu.
Saratov State Technical University named after Yuri Gagarin, Saratov, Russian Federation
This article is devoted to a task of rational distribution of porosity on bar section in case of a net bend. The solution of such task will allow to provide the necessary bearing capability in case of decrease in a material capacity of a design. A research purpose is to pick up the rational distribution law of porosity for the rectangular section of a bar in case of technical restrictions of production.
Keywords: porosity, net bend, rectangle, rationality, bar, normal tension, Jung's module. DOI: 10.22281/2413 -9920-2017-03-01-46-51
References
1. Kashtalyan Yu.A. Kharakteristiki uprugikh materialov pri vysokikh temperaturakh [The characteristics of elastic materials at high temperatures]. Kiev, Naukova dumka, 1970. 112 p.
2. Belov S.V. Poristye metally v mashinostroenii [Porous metals in mechanical engineering]. Moscow, Mashinostroenie, 1981. 247 p.
Authors' information
Stanislav M. Shlyakhov - doctor of Physical and Mathematical Sciences, Professor at Saratov State Technical University named after Yuri Gagarin.
Danila Yu. Gavrilov - graduate student at Saratov State Technical University named after Yuri Gagarin, [email protected].
Дата публикации (Date of publication): 25.03.2017
