УДК 378.14.015.62:378.147.227:539.3
Ю.С. Холодняк, А.В. Периг, И.А. Матвеев Yu.S. Holodnyak, A.V. Perig, I.A. Matveev
Донбасская государственная машиностроительная академия, г. Краматорск, Украина Donbass State Engineering Academy, Kramatorsk, Ukraine
СОВЕРШЕНСТВОВАНИЕ МЕТОДИКИ ПРОЧНОСТНОГО РАСЧЕТА ДВУТАВРОВОЙ БАЛКИ И ЕЕ ИЗЛОЖЕНИЯ В КУРСАХ ТЕХНИЧЕСКИХ ДИСЦИПЛИН
I-BAR STRENGTH COMPUTATION METHODOLOGY PERFECTION FOR TEACHING OF ENGINEERING DISCIPLINES
Исследован характер изменения эквивалентных напряжений по высоте сечения двутавровой балки, установлены места их максимума во взаимосвязи с параметрами двутавра и внутренними силовыми факторами, действующими в сечении. Проанализированы условия, при которых наибольшие эквивалентные напряжения в сечении двутавровой балки не превышают их допустимой величины. Предложена упрощенная графическая методика комплексной оценки прочности двутавровой балки по нормальным, касательным и эквивалентным напряжениям, которая рекомендуется студентам и преподавателям технических вузов, а также производственным специалистам в области прочностных расчетов.
Ключевые слова: двутавровая балка, эквивалентные напряжения, прочность, методика расчета, факторное пространство.
The equivalent stress distribution along the height of an I-bar cross-section has been studied. The zones of maximum equivalent stress in conjunction with the I-bar geometric parameters and power load at the section have been found. The conditions for maximum equivalent stresses have been specified to not exceed the allowable values. An improved method for the comprehensive assessment of I-bar strength against the normal, tangential and equivalent stresses has been proposed. This computational approach may be recommended to lecturers and students of technical universities as well as engineering specialists in the field of strength computations.
Keywords: I-bar, equivalent stresses, strength, design procedure, quotient space.
В нынешнюю эпоху стремительного развития техники особое значение приобретает подготовка вузами высокообразованных специалистов, владеющих солидным комплексом теоретических знаний и практических навыков, позволяющих создавать все более сложные образцы новой техники. Немаловажную роль в этом призвано сыграть совершенствование методик изложе-
ния базовых дисциплин инженерной подготовки, наполнение их новыми знаниями, углубляющими существующие теоретические представления.
Одними из основных в курсе сопротивления материалов, закладывающем фундамент инженерных расчетов, и ряде других смежных дисциплин являются вопросы прочности массово использующихся в машиностроении, строительстве и ряде других отраслей промышленности конструктивных элементов стержневой формы, работающих в условиях плоского изгиба и именуемых балками.
Указанные вопросы давно уже относятся к классическим, казалось бы, полностью исчерпанным в теоретическом плане и не нуждающимся в дополнительном изучении. Однако, как показывает практика прочностных расчетов, неясные вопросы здесь все же остаются.
Считается общепризнанным [1-4 и др.], что для обеспечения надежной эксплуатации балки должны соблюдаться условия ее прочности по нормальным (1), касательным (2) и эквивалентным (3) напряжениям:
°тах ^ м, (1)
ттах ^[т], (2)
^ [о]. (3)
Методики нахождения максимальных нормальных отах и максимальных касательных ттах напряжений, а также определения соответствующих допускаемых напряжений [о] и [т] достаточно полно и подробно изложены
в учебной литературе и в дополнительном анализе не нуждаются, чего нельзя сказать о методике нахождения максимальных эквивалентных напряжений в балках - аЭЭаХ, прежде всего потому, что места их действия в большинстве случаев неизвестны. Так, авторы упомянутых и многих других источников, приводя примеры полной проверки балок на прочность, рассматривают двутавровые балки и отмечают, без убедительных доказательств, что наибольшие эквивалентные напряжения в их поперечных сечениях действуют на стыках ребра с полками. По другим формам поперечных сечений подобных сведений практически нет.
В ходе исследований по данном вопросу [5, 6], ранее проведенных с участием авторов применительно к балкам прямоугольного и круглого поперечных сечений, было установлено, что наибольшие эквивалентные напряжения в зависимости от соотношения величин изгибающего момента М и поперечной силы Q, действующих в сечении, имеют место либо в крайних верхних и нижних слоях волокон балки, либо в волокнах нейтрального слоя, т.е. в местах действия отах и ттах. Их величины при соблюдении условий прочности (1) и (2) никогда не превышают [о], т.е. условие прочности по эквива-
лентным напряжениям (3) в балках названных сечений выполняется автоматически.
Настоящая работа является логическим продолжением упомянутых исследований. Цель работы - изучить характер изменения эквивалентных напряжений по высоте сечения двутавровой балки, увязать положение и величины их максимумов с внутренними силовыми факторами, действующими в сечении, и на этой основе внести соответствующие коррективы в методику силового расчета указанных балок.
На рис. 1 показана форма сечения двутавровой балки (со стандартными обозначениями) и эпюры нормальных а и касательных т напряжений по высоте этого сечения [7].
Рис. 1. Распределение нормальных и касательных напряжений по высоте сечения
двутавровой балки
Как известно, нормальные напряжения изменяются по линейному зако-
й
а,
шах
Ь
ну:
(4)
а касательные - по квадратичному:
(5)
или
где М и Q - соответственно изгибающий момент и поперечная сила, действующие в сечении; ]х - момент инерции сечения относительно оси х (перпендикулярной силовой плоскости); Sx(з,) - статический момент относительно оси х части сечения, лежащей выше уровня у.
Эквивалентные напряжения в точках сечения балки определяются по третьей или четвертой теориям прочности:
=Л
а2 + ат2,
(6)
где а = 4 - по третьей теории прочности и а = 3 - по четвертой.
Принимая во внимание скачкообразное уменьшение т при переходе от ребра к полке (см. рис. 1), нетрудно сделать вывод о том, что наибольшие эквивалентные напряжения возникают в пределах ребра двутавра, но не в его полках.
Подставив в формулу (6) выражения (4) и (5), получим
M
\
2
J У
V J X J
+ а
( QsXy ^
dJr
(7)
Выразим теперь Sxy через статический момент половины сечения дру-
тавра Sx:
SXy ) = Sx -1 • dy2.
x x 2
(8)
Подставив (8) в (7), получим закон изменения оэкв по высоте ребра двутавра:
( M Л
Ту
V J x J
2
+ а
Ж ( Sx- 2dy 2
(9)
С применением открытой системы компьютерной алгебры GNU CAS Maxima исследуем функцию (9) на экстремум. Возьмем первую производную от оэкв по y и приравняем ее к нулю:
2 • м2у 2 -а<22 (Sx-0,5 •dy2 )у
dll
= 0.
M
2y2 aQ2 (Sx -0,5 • dy2 )
(10)
2r2
d21
Решив уравнение (10), найдем значения у, при которых функция (9) может иметь экстремумы:
Jextr = 0, ±Oq(11)
Как видно из (11), число действительных корней уравнения (10) может быть различным в зависимости от знака подкоренного выражения. При а CQ2 > M2 Cd, или (M/Q) < (а CSX/d)0,5, подкоренное выражение в (11) по-
ложительно и число корней - три (один - на нулевом уровне и два симметричных относительно оси x). При а CQ2 CSX < M2 Cd, или (M/Q) > (а CS/d)0,5, подкоренное выражение в (11) или отрицательно, или равно нулю; в этом случае имеем один действительный корень: yextr = 0.
Таким образом, у каждого из номеров двутавра существует свое отношение M/Q, при достижении и превышении которого зависимость (9) имеет только одну возможную точку экстремума. Назовем это отношение критическим:
(M Л laSr
Q
, , d
У кр
Так, при а = 3 указанное отношение для двутавра № 10 равно 123,8 мм, для двутавра № 30 - 351,7 мм, для двутавра № 60 - 610,5 мм.
Чтобы выяснить, действительно ли является нулевой уровень местом экстремума и какой именно экстремум здесь имеет место (максимум или минимум), вычислим значение второй производной в этой точке от оэкв по у:
d2a 22 d -aQ2 Sx)
(0 ) = *-------. (12)
dy2 W QJxSxyia
Из (12) явствует, что при M2 Cd > a CQ2 CSX, или при (M/Q) > (M/Q)Kp, на нулевом уровне - минимум, так как выражение (12) положительно.
При M 2 СЙ < а CQ2 CSX, или при (M/Q) < (M/Q)Kp, выражение (12) отрицательно и в этом месте - максимум.
При M2 Ш = а CQ2 CSX, или при (M/Q) = (M/Q)Kp, выражение (12) равно нулю, поэтому для установления вида экстремума в данном случае прибегнем к помощи производных более высоких порядков.
Дальнейшие расчеты показывают, что третья производная от оэкв по у, как и вторая, при критическом отношении M/Q на нулевом уровне равна нулю, а четвертая - положительна:
d4%(0) = 3■'Ja'Qd >0, da SXJX
поэтому и в данном случае здесь - минимум [8].
Таким образом, при (M/Q) > (M/Q)Kp функция (9) имеет на нулевом уровне минимум, а при (M/Q) < (M/Q)Kp - максимум.
Выясним теперь вид экстремумов в точках с координатами
Jextr = ±0^2 • ad ( Sx -M 2 d ) •
(13)
Взяв вторую производную от оэкв по у и подставив в нее значения выражения (13), получим
Из (14) следует, что при а Щ2 CSx > M2 Cd, или при (M/Q) < (M/Q)Kp, это выражение положительно. Следовательно, на указанных уровнях функция (9) имеет минимумы.
Таким образом, при (M/Q) < (M/Q)Kp функция (9) имеет максимум на нулевом уровне и два минимума на уровнях с координатами, определяемыми по формулам (13). Необходимо отметить, что эти координаты зависят как от параметров профиля, так и от отношения M/Q в его сечении, в чем легко убедиться, преобразовав (13) к следующему виду:
При (M/Q) = 0 (т.е. при M = 0 и Q Ф 0) точки минимума оэкв максимально удалены от оси х:
Для двутавра № 10 это удаление составляет 101,1 мм, для двутавра № 30 -287,2 мм, для двутавра № 60 - 498,5 мм, что во всех случаях значительно больше половины высоты профиля ( h/2).
При увеличении отношения M/Q места минимума оэкв приближаются к оси х и при (M/Q)^(M/Q)Kp сливаются с местом максимума (у = 0):
Учитывая изложенное, можно констатировать, что в случае, когда (M/Q) < (M/Q)Kp, опасными по эквивалентным напряжениям могут быть как
dy2 eXtt d/x^M(2(2SX"-M2d]
(14)
нулевой уровень (у = 0), так и уровень точки K (yk = (h/2) - t), расположенной на границе ребра и полки (см. рис. 1). Поэтому важно определить такое (граничное) значение M/Q, при котором величины оэкв на обоих этих уровнях одинаковы. Тогда при (M/Q) < (M/Q)TV опасным будет нулевой уровень, а при (M/Q) > (M/Q)rp - уровень точки K.
Значение (M/Q)rp найдем из условия
Подставив (16) и (17) в (15) и решив полученное уравнение относительно M/Q, получим
При а = 3 величина (Ы/0)1.р составляет: для двутавра № 10 - 118,2 мм, для двутавра № 30 - 330,2 мм, для двутавра № 60 - 559,5 мм, т.е. везде
В случаях, когда (M/Q) > (M/Q)Kp, т.е. когда функция (9) имеет минимум на нулевом уровне, опасным по эквивалентным напряжениям является уровень точки K. Обратим внимание, что в данных случаях, как и в предыдущем, отношение M/Q превышает граничное.
Указанное обстоятельство позволяет сформулировать для всех трех возможных случаев единых принцип нахождения мест действия наибольших эквивалентных напряжений в сечении двутавровой балки: если отношение M/Q в сечении не превышает граничного значения, то а™х действует на нулевом уровне, а если превышает - на уровне точки K; при равенстве M/Q и (M/Q)IV наибольшие эквивалентные напряжения действуют одновременно на обоих указанных уровнях.
Сформулированный принцип иллюстрируется рядом кривых, изображенных на рис. 2 и являющихся графиками функции (9), построенными для двутавра № 30 при следующих исходных данных: а = 3; Q = 105 Н;
M = 30,703П06 Нйм (кривая 1), 33,020П06 Нйм (кривая 2), 34,115П06 Нйм (кривая 3), 35Д70П06 Нйм (кривая 4), 37,0000Ю6 Нйм (кривая 5).
°экв (°) = °KK
(15)
При этом оба значения оэкв определим по формуле (9):
(16)
(17)
(M/ Q)rp < (M/Q)Kp.
Теперь выясним, при каких условиях аэкв не превышают [а], т.е. когда выполняется условие прочности по эквивалентным напряжениям (3).
Начнем с нулевого уровня. В этом случае
а ЭкГ О?1 = Лт„„. (18)
106
105
104
103
102
101
100
99
/
" иэкв’ МПа /
Л
/4
^2
у, мм
0 20 40 60 80 100 120 ук
Рис. 2. Влияние отношения М/О на характер изменения эквивалентных напряжений по высоте двутавра № 30: (М/О)1 = 307,03 мм; (М/02 = 330,20 мм;
(М/03 = 341,15 мм; (М/04 = 351,70 мм; (М/05 =
= 370,00 мм
Если в сечении балки выполняется условие прочности по касательным напряжениям [2], то с учетом (18) имеем
а эах
Тэах Ф]. (19)
у/а
Но [т], как известно, связано с [а]:
|т] = М (20)
>/а
Тогда из (19) и (20) следует, что -ЭХ ^ [-], т.е. что здесь автоматически выполняется и условие прочности по эквивалентным напряжениям (3).
Перейдем теперь к уровню точки К. Запишем выражение для оэкв на этом уровне:
/ \2
( МУк Л
V 3 х /
+ а
= [а]
Ук Лх [а]
+ й ж
* (-1 • луі Лх [а] V х 2 ^
(21)
Расчеты показывают, что значение корня в (21) для любого номера двутавра может быть больше единицы, даже если величины М и Q не превышают своих максимально возможных (допустимых) значений:
[М ] = Wx [а], й]х [а]
[й] =
Sr^¡а
(22)
(23)
Следовательно, условие прочности по эквивалентным напряжениям (3) выполняется в данном случае не при любых значениях М и Q, а только при таких, которые не выходят за определенные границы.
Для установления таких границ приравниваем (21) к величине максимального эквивалентного напряжения, еще удовлетворяющего условию (3):
_тах г_-. азкЯ = №
получим
[Ч
Ук
Лх [а]
+й
—( 8х -1 • йуК Лх [а] V х 2 ^
= [а].
(24)
Отсюда - равенство единице квадратного корня в выражении (24), а значит, и его подкоренного выражения:
(
Ук
2
л м
Введем обозначения:
* (\-1 • іуК
Лх [а] V х 2 ^
Ук = 1
= 1.
Л[а] Р’
7а(х - 0,5 • йуК) й]х [а]
(25)
(26)
(27)
Подставим (26) и (27) в (25), получим:
= 1.
Это и есть уравнение линии, ограничивающей геометрическое место точек с координатами М и й, где выполняется условие прочности (3), т.е. дуги эллипса с полуосями Р и у:
При построении указанной линии необходимо ограничить величины М и Q их допустимыми значениями (22) и (23):
Введение ограничений (28) и (29) заменяет часть дуги эллипса прямолинейными отрезками, что позволяет исключить из ограничиваемого факторного пространства все те точки, значения М и Q в которых хотя и удовлетворяют условию прочности (3), но нарушают условия (1) и (2). Таким образом, факторное пространство становится универсальным, пригодным для проверки сразу трех названных условий прочности.
Проиллюстрируем изложенное примером полной проверки прочности двутавровой балки № 30 Ух = 7080П04 мм4, уК = 139,8 мм, й = 6,5 мм, а = 3, Sx = 268П03 мм3, ^х = 472П03 мм3, [о] = 150 Н/мм2), представленной на рис. 3 [7].
Для построения факторного пространства выполним необходимые расчеты, согласно которым Р = 75,966 кНй; у = 194,910 кН; [М] = 70,800 кНй; [2] = 148,715 кНй.
Факторное пространство для приведенной балки (заштрихованная площадь) представлено на рис. 4, где координаты пронумерованных точек являются значениями поперечных сил и изгибающих моментов в соответствующих сечениях балки (рис. 3): 21 = 81,5 кН; М1 = 0 кНй; 22 = 18,5 кН; М2 = 63 кНй; 23 = 73,5 кН; М3 = 63 кНй; 24 = 73,5 кН; М4 = 4,1 кНй; 25 = 73,5 кН; М5 = 44,1 кНй; 26 = 73,5 кН; М6 = 0 кНй; 27 = 0 кН; М7 = 66,4 кНй.
Ук й]х [а]
М <[М], й <[й].
(28)
(29)
Рис. 3. Схема балки с эпюрами поперечных сил й и изгибающих моментов М
Поскольку ни одна из обозначенных точек не выходит за границы факторного пространства (см. рис. 4), то в каждом из соответствующих сечений балки должны одновременно выполняться условия прочности по нормальным (1), касательным (2) и эквивалентным (3) напряжениям. Убедимся в этом.
Рис. 4. Факторное пространство для приведенной балки
Вертикальные координаты всех указанных точек не превышают [М], значит, во всех соответствующих сечениях выполняется условие прочности (1).
Горизонтальные координаты этих точек не превышают [Q], поэтому в соответствующих сечениях балки выполняется и условие прочности (2).
В сечениях 1, 4 и 6 отношения значений M и Q меньше граничного (330,21 мм), поэтому, как отмечалось раньше, условие прочности (3) в них выполняется автоматически, так как выполнено условие (2).
В остальных сечениях - 2, 3, 5 и 7 - отношения M/Q превышают граничное, но поскольку соответствующие точки не выходят за границу факторного пространства, очерченную дугой эллипса, то и в этих сечениях выполняется условие прочности (3). Нетрудно таким же образом убедиться в выполнимости всех условий прочности в любых других сечениях приведенной балки.
Факторное пространство, изображенное на рис. 4, пригодно для полной проверки прочности балок данного профиля (№ 30) при любом варианте их нагружения.
Выполнив подобные построения для всех других номеров двутавра, выпускаемых промышленностью, можно создать солидную справочную базу для прочностных расчетов этих балок по упрощенной графической методике, изложенной выше. Такая методика может быть полезна как студентам и преподавателям технических вузов, так и производственным специалистам в области прочностных расчетов.
Таким образом, ходе исследования установлена взаимосвязь мест действия наибольших эквивалентных напряжений с внутренними силовыми факторами в сечении двутавровой балки. Показано, что такие напряжения могут возникать на уровне центра тяжести сечения, если отношение величин изгибающего момента и поперечной силы в этом сечении меньше определенной (граничной) величины, зависящей от параметров двутавра и принятой теории прочности, а также на стыке ребра с полкой, если указанное отношение превышает граничное; при равенстве этого отношения и его граничной величины наибольшие эквивалентные напряжения возникают сразу на обоих названных уровнях.
Доказано, что в случае максимума эквивалентных напряжений на уровне центра тяжести сечения условие прочности по этим напряжениям выполняется автоматически, если в сечении выполнимо условие прочности по касательным напряжениям.
Установлено, что выполнимость условия прочности по эквивалентным напряжениям при их максимуме на стыке ребра с полкой зависит от комбинации величин изгибающего момента и поперечной силы в сечении балки. Это условие выполняется, если указанная комбинация в виде соответствующей точки факторного пространства с координатами «изгибающий момент -поперечная сила» не выходит за его границу, очерченную дугой эллипса с полуосями, определяемыми параметрами двутавра, допускаемым нормальным напряжением и принятой теорией прочности.
По результатам исследования предложена упрощенная (графическая) методика комплексной оценки прочности двутавровой балки, основанная на использовании построенного факторного пространства с соответствующими ограничениями, вытекающими из условий прочности по нормальным и касательным напряжениям. Она может быть рекомендована к применению как в учебных, так и в практических целях.
Работы по данной тематике целесообразно продолжить в направлении интегрирования предложенной методики в соответствующую компьютерную технологию прочностных расчетов.
Список литературы
1. Сопротивление материалов: учебн. для техн. вузов / Г. С. Писаренко [и др.]; под ред. Г.С. Писаренко. - Киев: Вища шк., 1986. - 775 с.
2. Дарков А.В., Шниро Г.С. Сопротивление материалов: учебник для техн. вузов. - М.: Высшая школа, 1989. - 624 с.
3. Цурпал И.А. Краткий курс сопротивления материалов: учеб. пособие. -Киев: Вища шк., 1989. - 311 с.
4. Матвиенко Ю.Г., Сапунов В.Т. Сопротивление материалов в задачах и решениях: учеб. пособие. - М.: Эдиториал УРСС, 2001. - 88 с.
5. Холодняк Ю.С., Подлесный С.В., Роменский ЕЮ. Прочностные расчеты в машиностроении и необходимость учета в них эквивалентных напряжений // Вісник Донбаської державної машинобудівної академії. - 2006. -№ 1Е (6). - С. 74-80.
6. Особенности прочностного расчета балок и их учет в изложении технических дисциплин / Ю.С. Холодняк [и др.] // Качество образования -управление, сертификация, признание: сб. науч. работ междунар. науч. метод. конференции / под общ. ред. С.В. Ковалевского; ДГМА. - Краматорск, 2011. - С. 463-469.
7. Збірник розрахунково-графічних завдань з курсу «Опір матеріалів» (для студентів всіх механічних спеціальностей денної форми навчання) / Л.В. Кутовий, В. А. Овчаренко, Ю.С. Холодняк [та ін.]; ДДМА. - Краматорск,
2007. - 220 с.
8. Смирнов В.И. Курс высшей математики. - СПб.: БХВ-Петербург,
2008. - Т. 1. - 624 с.
Получено 1.11.2012
Холодняк Юрий Сергеевич - кандидат технических наук, доцент, Донбасская государственная машиностроительная академия (84313, Украина, Донецкая область, г. Краматорск, ул. Шкадинова, 72, е-mail: texmex@dgma. donetsk.ua).
Периг Александр Викторович - кандидат технических наук, старший преподаватель, Донбасская государственная машиностроительная академия (84313, Украина, Донецкая область, г. Краматорск, ул. Шкадинова, 72, е-mail: [email protected], [email protected]).
Матвеев Иван Анатольевич - студент, Донбасская государственная машиностроительная академия (84313, Украина, Донецкая область, г. Краматорск, ул. Шкадинова, 72, е-mail: [email protected]).
Holodnyak Yuri Sergeevich - Candidate of Technical Sciences, Associate Professor, Donbass State Engineering Academy (84313, Ukraine, Kramatorsk, Shkadinova st., 72, е-mail: [email protected]).
Perig Aleksandr Viktorovich - Candidate of Technical Sciences, Senior Lecturer, Donbass State Engineering Academy (84313, Ukraine, Kramatorsk, Shkadinova st., 72, е-mail: [email protected], olexander.perig@ gmail.com).
Matveev Ivan Anatolyevich - Student, Donbass State Engineering Academy (84313, Ukraine, Kramatorsk, Shkadinova st., 72, е-mail: matveev.ivan. [email protected]).