CC BY
9
2017 Математика и механика № 46
УДК 539.3
DOI 10.17223/19988621/46/8
Д.Ю. Гаврилов
НАПРЯЖЕННОЕ СОСТОЯНИЕ БРУСА ПРЯМОУГОЛЬНОГО СЕЧЕНИЯ, ВЫПОЛНЕННОГО ИЗ МАТЕРИАЛА ПОРИСТОЙ СТРУКТУРЫ ПРИ ЧИСТОМ ИЗГИБЕ
Статья посвящена задаче рационального распределения пористости по высоте сечения нагруженного элемента конструкции при чистом изгибе. Решение такой задачи позволит обеспечить необходимую несущую способность при снижении материалоемкости конструкции. Целью исследования является разработка методики решения задачи чистого изгиба бруса прямоугольного поперечного сечения, выполненного из материала пористой структуры с учетом технических ограничений производства, и получение численного решения задачи.
Ключевые слова: пористость, чистый изгиб, брус.
Рассмотрим чистый изгиб бруса, выполненного из материала (стали) пористой структуры с переменной по высоте сечения пористостью р. Ограничимся случаем упругого деформирования, полагая, что максимальные напряжения в брусе ст^ не превышают предела текучести материала стТ . Известно, что упругая характеристика материала (модуль Юнга Е) является функцией пористости материала. Экспериментальные данные, описанные в [1, 2], представлены для пористого железа, изготовленного из порошка марки ПЖМ. Полученные двусторонним прессованием в пресс-формах брикеты спекали при температуре 1470 К в среде водорода в течение 2 ч. После механической обработки образцы подвергали отжигу для снятия наклепа в поверхностном слое и восстановления окислов, образовавшихся в материале от окончания спекания до завершения шлифования. Отжиг производили при температуре 1170 К в среде водорода в течение 1 ч.
На основе представленных выше экспериментальных данных зависимости Е(р) и стТ (р) могут быть представлены полиномами:
Е = а1 + а2р + а3р2 ; (1)
стТ = с1 + с2р + с3 р2 . (2)
Сглаживая заданную функцию методом наименьших квадратов (МНК), коэффициенты а^, с^ найдем с помощью минимизации отклонения сглаживающей функции от заданных точек в некотором среднеинтегральном смысле.
Конечным результатом МНК будут являться значения соответствия стТ , Е пористости в границах от 0 до 0.4. Максимальное значение пористости выбрано из соображения технологических возможностей по изготовлению конструкций.
На рис. 1 представлены графики функций соответственно (1) и (2) при значениях коэффициентов, приведенных в таблице.
Значения коэффициентов аппроксимирующего полинома, МПа
a1 209285,7143 Cl 196,0368975
a2 -535000 C2 -645,964465
a3 321428,5714 C3 627,5555095
Рис. 1а. Графики зависимости модуля Юнга (а) и нормальных напряжений (b) от пористости
Fig. 1a. Diagrams of Young's modulus (a) and normal stresses (b) as a function of porosity
Положим теперь, что пористость р меняется по высоте сечения у:
р = р( у), 0 < у < к/2,
где у - координата высоты полусечения бруса, к - высота полусечения бруса.
Тогда Е = Е(у), т.е. модуль Юнга есть функция координаты высоты (у) поперечного сечения бруса.
Задача изгиба, таким образом, сводится к задаче изгиба бруса, выполненного из неоднородного материла с переменной по высоте сечения упругой характеристикой Е = Е(у) и переменным пределом текучести стг = стг (у).
Обозначим изгибающий момент на брусе через М. При чистом изгибе бруса соблюдается гипотеза плоских сечений, т.е. сечения после деформаций остаются плоскими и нормальными к оси бруса, а геометрические размеры неизменными.
Соответственно закону Гука при растяжении-сжатии для изогнутого бруса определим закон изменения напряжений по высоте сечения:
ст = Ее , (3)
где е = етах - максимальная деформация изогнутого бруса.
Величина изгибающего момента в сечении определится по формуле
/ 2
M = 21 abydy,
(4)
где Ь - ширина поперечного сечения бруса
Цель исследования - подобрать такой закон распределения модуля Юнга Е и, следовательно, пористости р по высоте сечения бруса, чтобы получить максимально возможный изгибающий момент при ограничениях на напряжение ст<стт и на пористость ртт < р < ртах, задаваемую техническими возможностями производства.
Данная постановка является сложной задачей линейного программирования, для решения которой используем приближенный подход. Суть приближенного решения заключается в следующем.
Разбиваем полусечение бруса по высоте сечения на п элементов: у}-,
] = 1, 2,..., п, где у;- - расстояние от горизонтальной оси сечения до соответст-
Н
вующего слоя с шагом, 5 = — (рис. 2). В пределах каждого слоя получаем все
2п
характеристики постоянными, меняющимися от слоя к слою.
En, pn
EP Pj
E3, Рз
E2, Р2
El, Pi
У
E
П
Рис. 2. Схема разбиения сечения бруса по высоте Fig. 2. Partitioning scheme of a beam cross section along the height
Соответственно получим
У = 5, У; = У}-1 + 5, ] = 2, п.
Максимальные напряжения в слоях будут равны
ст; = Е]8; ~ат; , ; = 1 2,- •., п. (5)
Модуль Юнга Е = Е(у;-) зависит от расстояния от оси сечения до слоя у;-, при этом имеем ограничение по напряжению:
СТ<СТГ ( У; ). (6)
Поскольку имеем Е(р), аТ (р) в виде зависимостей (1), (2), то из условия (6) получаем
8 ; {а1 + а2 Р; + а3 Р; 2 } < ^ + С2 Р ^ + ^ Р ^ 2 ,
где значения коэффициентов а и с приведены в таблице, р; - пористость в ;-м слое.
Соответственно получаем условие для определения р;■ на высоте у.
0 ао 8 со а 8 С
р2 + р, 2 _ 2 + 1 ; _ 1 < 0. (7)
аз 8; С3 аз 8; С3
Обращаем неравенство (7) в равенство и решаем его как обыкновенное квадратное уравнение, находя положительный корень. Найденное значение пористости следует согласовать с условием ртт < р; < ртах . В случае выхода искомого
значения р ; за указанные пределы принимаем граничное значение пористости.
Найдя пористость р ; , находим Е; .
Величина изгибающего момента определится по формуле
М « 2£ а;.Ь5 У; = 2£ Е;8}Ь5 у} . (8)
;=1 ;=1
В качестве примера примем А = 20 см, п=20, 8тах = 0.000952381, Ь = 10 см. Получим значение пористости для каждого элементарного слоя и представим в виде эпюры на рис. 3. Эпюра распределения нормальных напряжений изображена на рис. 4.
Суммарный изгибающий момент, воспринимаемый брусом с распределенной по высоте сечения пористостью, будет равен 3.01 кНм. Стоит отметить, что увеличение числа слоев разбиения бруса будет приводить к уточнению значения вычисленного изгибающего момента.
Для оценки роли пористости при изгибе бруса сравним значение воспринимаемого изгибающего момента в случае изменяющейся по высоте сечения пористости с моментом в случае постоянной (средней) пористости при тех же нагрузках.
Расчетные формулы для момента будут следующие:
Для средней пористости:
Еср = а1 + а2 рср + а3 рср2 , а; = Еср8; ,
М(1) « 2Еср £ 8;Ь 5 У; , (9)
] =1
где Еср - значение модуля Юнга, соответствующее средней пористостирср = 0.33.
Для переменной пористости
Еср = E ( p, ).
M(2) « 2^ Е, ejb 5 y .
(10)
j=1
В случае если равномерно распределить среднюю пористость рср = 0.33 по
всему поперечному сечению, получим момент равный 1.844 МНм. Эпюра нормальных напряжений в этом случае отражена на рис. 5.
р, %
Страц, МПа
стср, МПа
> ^,8 0.075=^1 40 40 ДН95 J<1S7 —^99 3 tss ¿55 -/43 )
131jT 43^= 55 J 53 f 61J <
Рис. 3. Эпюра значений рационально подобранной пористости
Fig. 3. Plot of values of rationally chosen porosity
Рис. 4. Эпюра значений нормальных напряжений при рационально подобранной пористости Fig. 4. Plot of values of normal stresses at rationally chosen porosity
Рис. 5. Эпюра значений нормальных напряжений при средней по сечению пористости
Fig. 5. Plot of the values of normal stresses at average porosity in cross section
Анализируя результат, видим, что момент при изменяемой по высоте сечения пористости выше, чем момент при средней пористости на 63 %.
Таким образом, аналитически подтверждено, что при варьировании пористости по высоте сечения можно повысить несущую способность бруса при его изгибе.
ЛИТЕРАТУРА
1. Кашталян Ю.А. Характеристики упругих материалов при высоких температурах. Киев: Наукова думка, 1970. 112 с.
2. Белов С.В. Пористые металлы в машиностроении. М.: Машиностроение, 1981. 247 с.
Статья получена 11.07.2016 г.
Gavrilov D.Yu. (2017) THE STRESS STATE OF A RECTANGULAR CROSS SECTION BEAM MADE OF A MATERIAL WITH POROUS STRUCTURE AT PURE BENDING. Tomsk State University Journal of Mathematics and Mechanics. 46. pp. 64-69
DOI 10.17223/19988621/46/8
The article is devoted to the problem of rational distribution of porosity along the height of the cross section of a loaded structural element at pure bending. We consider pure bending of the beam which is made of a material (steel) with a porous structure of variable porosity along the height of cross section. The research is limited to the case of an elastic deformation under the assumption that the maximum stress in the beam does not exceed the yield stress of the material. It is known that Young's modulus is a function of material porosity.
The porosity is assumed to be variable along the height of the cross section. Therefore, Young's modulus is also a function of the height coordinate of the beam cross section.
The problem of bending is reduced to the problem of bending of a beam which is made of a heterogeneous material with a variable elastic characteristic along the height of cross section and variable yield stress.
According to Hooke's law, at the tension-compression state of bent beam, the law of stress variation along the height of cross section is determined, which makes it possible to calculate the value of the bending moment.
The aim of the study is to select the distribution law of both Young's modulus and, hence, porosity along the height of beam cross section in order to obtain the greatest possible bending moment under restrictions on the stress and porosity.
A numerical solution of the problem corroborates that the bending moment is 63% higher at rationally chosen porosity than that at average porosity. The results of calculations indicate an optimization of the porosity distribution over the cross section.
Keywords: porosity, pure bending, beam.
GAVRILOV Danila Yur'evich (Saratov State Technical University named after Yuri Gagarin, Saratov, Russian Federation)
E-mail: [email protected]
REFERENCES
1. Kashtalyan Yu.A. (1981) Kharakteristiki uprugikh materialov pri vysokikh temperaturakh [Characteristics of elastic materials at high temperatures]. Kiev: Naukova dumka.
2. Belov S.V. (1981) Poristye metally v mashinostroenii [Porous metals in mechanical engineering]. Moscow: Mashinostroenie.
