Механика деформируемого твёрдого тела
УДК 539.376
ЧИСТЫЙ ИЗГИБ БАЛКИ В УСЛОВИЯХ ПОЛЗУЧЕСТИ ИЗ РАЗНОСОПРОТИВЛЯЮЩЕГОСЯ МАТЕРИАЛА
А. М. Локощенко, К. А. Агахи, Л. В. Фомин
Московский государственный университет им. М. В. Ломоносова,
НИИ механики МГУ,
119992, Россия, Москва, Мичуринский пр-т, 1.
E-mails: [email protected], [email protected], [email protected]
Решены задачи о чистом изгибе балок прямоугольного и круглого поперечных сечений при установившейся ползучести с учётом различных характеристик материала при растяжении и сжатии. Определяющие соотношения установившейся ползучести приняты в виде дробно-степенных функций. Приводится сравнение решений задач для балок квадратного и круглого сечений с одинаковыми значениями осевых моментов инерции.
Ключевые слова: балки, изгиб, установившаяся ползучесть, разносопротивляе-мость материала, дробно-степенная модель, квадратное сечение, круглое сечение.
Введение. В настоящее время опубликовано довольно много работ, посвященных изучению чистого изгиба балок в процессе установившейся ползучести. Проведенный обзор некоторых из них [1-13] послужил поводом для дальнейшего исследования этой задачи. Задачи чистого изгиба балок при ползучести для материалов, имеющих одинаковые свойства на растяжение и сжатие, рассмотрены, например, в работах [1, 3-6, 10-13]. В этих работах в качестве определяющих соотношений в основном используются степенные и экспоненциальные функции. В работах [2, 7-9] рассматривается чистый изгиб балок в процессе установившейся и неустановившейся ползучести, в которых учитывается разносопротивляемость материала растяжению и сжатию. При этом для описания ползучести при изгибе авторы используют как степенные [2,7], так и экспоненциальные [9] определяющие соотношения. Во всех перечисленных работах, за исключением работ [12,13], используется феноменологический подход к описанию процесса ползучести. В работах [12,13] предложена методика решения одномерных краевых задач изгиба балки на основе структурной модели среды.
В настоящей статье исследуется чистый изгиб балок с различными формами поперечных сечений в процессе установившейся ползучести. В отличие от рассмотренных работ, для учета разносопротивляемости материала растя-
Александр Михайлович Локощенко (д.ф.-м.н., проф.), зам. директора, НИИ механики МГУ. Камилла Абдул Гусейн кызы Агахи (к.ф.-м.н., доц.), и. о. зав. лаборатории, лаб. прочности и ползучести при высоких температурах НИИ механики МГУ. Леонид Викторович Фомин, ведущий инженер, НИИ механики МГУ.
жению и сжатию принимается сингулярная дробно-степенная модель с различными значениями пределов кратковременной прочности при растяжении и сжатии [14,15].
Проводится сравнение полученных решений для балок круглого и квадратного поперечного сечений с равными значениями их осевых моментов инерции при действии одинакового изгибающего момента.
1. Чистый изгиб балки прямоугольного поперечного сечения. Рассмотрим чистый изгиб балки в условиях высокотемпературной ползучести. Материал балки при соответствующей температуре имеет разные значения пределов кратковременной прочности при растяжении аы > 0 и сжатии <752 < 0 (отношение пределов прочности а = —аь2І(Тьі- Примем, что сечение балки имеет две оси симметрии (Ох и Оу), изгибающий момент М действует в плоскости уОг. Высота сечения балки по оси Оу равна Н, ширина Ь, длина I удовлетворяет неравенствам I » Н, I » Ь.
В этом случае гипотеза плоских сечений имеет вид
Р = Х(У~Уо), (!)
где р = рг — скорость деформации ползучести; Х~ скорость изменения кривизны балки; у — координата, отсчитываемая от срединной линии балки (—0,5Н ^ у ^ 0,5Н)\ уо — координата нейтральной поверхности, на которой отсутствуют напряжения (сг(уо) = <тг(уо) = 0). Смещение нейтральной поверхности от оси симметрии (у = 0) балки при ползучести в рассматриваемой задаче происходит вследствие разносопротивляемости материала растяжению и сжатию.
Определяющие соотношения ползучести в одноосном случае примем в виде, предложенном С. А. Шестериковым и М. А. Юмашевой [14,15], в котором зависимость р от напряжения а имеет дробно-степенной вид:
р = а(—= а У. (2)
V л /(пі,. — п)(ҐТ — пі,- ) /
Таким образом, напряжённо-деформированное состояние изогнутой балки в любой момент времени определяется осевыми напряжениями а = и (у) и осевыми деформациями ползучести р = р(у,і). Мгновенными упруго-пластиче-скими деформациями пренебрегаем.
Уравнения равновесия балки имеют следующий вид:
/уо Г-Н/2 гуо гН/2 ^
{(Т-)(1у + / {(Т+)(1у = 0, / ((Т-)у(1у+ {(Т+)у(1у=—, (3)
-я/2 3уо 7-Я/2 3уо 0
где <т_, (7-1- — напряжения в сжатой и растянутой зонах балки соответственно. Введём следующие безразмерные переменные:
а =—а ь2/&ы, а = а/аы, 1=1А д> = 4М дг = 2ЛГ - _ _ н_<1х - =2У_
ЪН2аЬ1 ’ ЪНаЬ1 ’ Х <Ш 2 АсМ’ У Н'
С учётом введённых безразмерных переменных определяющее соотношение (2) преобразуется к следующему виду:
ф
сії
и
л/(1 — а)(а + а)
(4)
На рис. 1 представлены графики зависимости модуля скорости деформации |ф/(Й| от модуля безразмерного напряжения |<т| (4) при п = 1.
Согласно гипотезе плоских сечений (1) и определяющему соотношению (2) в безразмерных переменных имеем
ф
СІЇ
а2
(1 - а) (а + а)
га/2
= §(5-Л). (5)
Преобразуем последнее выражение: ¡.2
а
Рис. 1. Графики зависимости \с1р/(Щ от \а\ при а = 1,5: 1 — область растяжения; 2 — область сжатия
(1 - а) (а + а)
С& 1,1 - а.
\ 2/га
(у-уо)) =В(у)- (6)
V СІІ
Представим выражение (6) в виде квадратного уравнения относительно безразмерного напряжения а:
(1 + В)а2 + В (а — 1)<т — аВ = 0. (7)
В результате выражение для безразмерных напряжений принимает следующий вид:
-В{а - 1) ± л/В2(а - I)2 + 4(1 + В)аВ
а+'~ -----------------2(ГТб)-----------------• (8)
Система уравнений равновесия (3) в безразмерных переменных приводится к виду
/Уо /•!
[а-)с!у+ / {(т+)(1у = 0; М =
■1 •'УО
ГУО /•!
/ (а-)уйу+ / (а+)уйу- (9)
'-1 Но
Подставив выражения а_|_ и <т_ из (8) в (9) и учитывая (6), получим систему двух уравнений относительно уо и X-
В таблице в столбцах 1 и 3 приведены результаты решения системы (9) при М = 0,5, а = 1,5 для некоторых значений п. На рис. 2 сплошными линиями изображены вычисленные согласно (6), (8), (9) эпюры напряжений в балке при этих же значениях констант. На рис. 3 приведены зависимости смещения нейтральной оси у о от показателя степени п при а = 1,5.
2. Чистый изгиб балки круглого поперечного сечения. В настоящем пункте рассматривается чистый изгиб при ползучести балки круглого поперечного сечения радиуса К, изготовленной из материала с разными значениями преде-
Результаты численного решения систем уравнений для балок прямоугольного (столбцы 1 и 3) и круглого (столбцы 2 и 4) сечений
те Уо X X
1 2 3 4
1 -0,045 -0,038 0,696 0,803
3 -0,061 -0,049 0,148 0,172
15 -0,077 -0,061 0,00001 0,0000131
Т-, г. ^ Рис. 3. Графики зависимости смещения
Рис. 2. Эпюры безразмерного напряжения „ _
_ ^ нейтральной оси уо от показателя те для
а для балок прямоугольного и круглого се- г
„ 1 n ir 1 г балок прямоугольного 1) и круглого 2)
чении при те = 1, 3, 15 и а = 1,5 4 ' „ 4 '
поперечных сечении
лов прочности при растяжении <751 и сжатии (Т52 (отношение пределов прочности по-прежнему равно a =
= —vb2/vbi)- Задача решается с учётом гипотезы плоских сечений. Примем, что сечение балки имеет две оси симметрии (Ох и Оу), изгибающий момент действует в плоскости yOz. Высота сечения балки по оси у равна 2R, ширина b = 2Rcos(p,
Lp — угол между радиусом и осью Ох (см. рис. 4), длина I удовлетворяет неравенству I R.
Как и в п. 1, напряжённо-де-формированное состояние изогнутой балки в любой момент времени определяется осевыми напряжениями <7 = и(у) и осевыми деформациями ползучести р = р(у, t). Мгновенными упруго-пластическими деформациями по-прежнему пренебрегаем.
В этом случае, как и в п. 1, гипотеза плоских сечений имеет вид
Рис. 4. Круглое и квадратное поперечное сечение балок с одинаковыми осевыми моментами инерции
Р = х(У~Уо), y = R sm<p, y0 = Rsm<p0. Введём безразмерные переменные
a = — оь2І(Уь\і д’ = о’/о’ы> t = tA;
(10)
При рассмотрении изгиба балок прямоугольного и круглого сечений все уравнения (5), (7), (8) совпадают, различие наблюдается только в уравнениях равновесия:
где сг+;_ выражаются через В (у) по формуле (8), а В (у) с учётом (6), (10) в безразмерных координатах принимает вид
Таким образом, в системе (11) величины и сро являются неизвестными.
3. Сравнение решений задач об изгибе балок круглого и квадратного поперечного сечения. Проведём сравнение решений задач об изгибе балок круглого и квадратного поперечного сечения (как частный случай прямоугольного поперечного сечения при Ъ = Н) при одном и том же изгибающем моменте с учётом равенства осевых моментов инерции. Из равенства осевых моментов инерции следует, что
Отсюда получим соотношение между радиусом К круглого поперечного сечения и стороной Н квадратного поперечного сечения:
Соотношение размеров круглого и квадратного поперечных сечений балок показано на рис. 4. Сравнивая выражения для безразмерных изгибающих моментов в случае квадратного М и круглого М поперечных сечений, получаем
В п. 1 для задачи изгиба при ползучести балки прямоугольного поперечного сечения был произведен расчёт для М = 0,5. Аналогично выполним расчёт для балки круглого сечения, нагруженной моментом М = 0,336.
В таблице в столбцах с номерами 2 и 4 приведены результаты решения
чений п. На рис. 2 штриховыми линиями представлены эпюры распределения напряжений ¿г для балки круглого поперечного сечения.
Из рис. 2 видно, что при п = 1 распределение напряжений близко к линейному. С ростом показателя п наблюдается перераспределение напряжений к характерному для идеально-пластического материала. Наряду с изменением эпюр напряжений в зависимости от показателя степени п наблюдается
/
/
7Г/2
7Г/2
(11)
Я4 тгЯ4
12 4
Я = (Зтг)1/4Я и 1,752Д.
М и 0,672М.
системы уравнений (11) при М = 0,5, М = 0,336, а = 1,5 для некоторых зна-
смещение нейтральной оси (сг(уо) = 0) от оси симметрии балки за счет раз-носопротивляемости материала.
На рис. 3 приведены графики зависимости смещения нейтральной оси уо от показателя степени п для балок круглого и квадратного поперечных сечений.
Сравнение эпюр распределения безразмерных напряжений и при различных значениях п для балок круглого и прямоугольного поперечных сечений (при одинаковых осевых моментах инерции) под действием одного и того же изгибающего момента показало, что величина этих напряжений по абсолютной величине больше для балки круглого поперечного сечения. Это можно объяснить тем, что вследствие гипотезы плоских сечений выдерживать основную часть нагрузки должны волокна, наиболее удаленные от плоскости симметрии, а у балки круглого сечения в этой области ширина меньше.
Из рис. 3 видно, что для любого п смещение нейтральной оси для балки круглого поперечного сечения меньше.
Легко показать, что в частном случае материала балок с одинаковыми свойствами (a = 1) в аналогичной задаче чистого изгиба смещения нейтральной оси не происходит (уо = 0).
Также необходимо отметить, что используемое дробно-степенное определяющее соотношение наиболее адекватно позволяет описывать процесс ползучести, так как оно ограничивает уровень допустимых напряжений пределами кратковременной прочности при растяжении и сжатии.
Заключение. Решены задачи о чистом изгибе балки при ползучести с учётом различия свойства материала при растяжении и сжатии. Определяющее уравнение принято в виде сингулярной дробно-степенной зависимости скорости ползучести от напряжения, мгновенные деформации не принимаются во внимание. Решения получены для прямоугольного и круглого поперечных сечений балок. На основании эпюр распределения напряжений по сечению балок проведено сравнение полученных результатов.
Представленные эпюры напряжений подтверждают гипотезу о смещении нейтральной оси балки при чистом изгибе в случае разносопротивляющегося материала.
Показано влияние показателя п в дробно-степенной зависимости на характер распределения напряжений.
Работа выполнена при поддержке РФФИ (проекты № 11-08-01015, 11-08-00007).
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. Качанов Л. М. Теория ползучести. М.: Физматгиз, 1960. 456 с. [Kachanov L. М. Theory
of Creep. Moscow: Fizmatgiz, 1960. 456 pp.]
2. Работное Ю.Н. Сопротивление материалов. М.: Физматгиз, 1962. 456 с.
[Rabotnov Yu. N. Strength of Materials. Moscow: Fizmatgiz, 1962. 456 pp.]
3. Работное Ю. H. Ползучесть элементов конструкций. М.: Наука, 1966. 752 с.
[Rabotnov Yu. N. Creep of Structural Elements. Moscow: Nauka, 1966. 752 pp.]
4. Boyle J. Т., Spence J. Stress Analysis for Creep. London: Buttenworth, 1983. 296 pp.; русск.
пер.: Бойл Дж., Спенс Дж. Анализ напряжений в конструкциях при ползучести. М.:
Мир, 1986. 360 с.
5. Радченко В. П., Еремин Ю. А. Реологическое деформирование материалов и элементов конструкций. М.: Машиностроение-1, 2004. 264 с. [Radchenko V.P., Eremin Yu. А.
Rheological Deformation and Failure of Materials and Structural Elements. Moscow: Mashinostroenie-1, 2004. 264 pp.]
6. Сорокин О. В., Самарин Ю.П., Одинг И. А. К расчёту ползучести балок при изгибе// Докл. АН СССР, 1964. Т. 157, №6. С. 1325-1328; англ. пер.: Sorokin О. V., Samarin Yu. Р., О ding LA. Calculation of the Creep of Beams During Bending// Sov. Phys., Dokl.. Vol. 9. Pp. 728-733.
7. Никитенко А. Ф., Соснин О. В. Изгиб балки из материала с разными характеристиками ползучести при растяжении и сжатии // Проблемы прочности, 1971. №6. С. 67-70; англ. пер.: Nikitenko A. F., Sosnin О. V. Beam flexure for materials possessing different creep characteristics in tension and compression // Strength of Materials, 1971. Vol. 3, no. 6. Pp. 693-697.
8. Горев Б. В. К расчёту на неустановившуюся ползучесть изгибаемого бруса из материала с разными характеристиками на растяжение и сжатие / В сб.: Динамика сплошной среды: Сб. научн. тр. Вып. 14. Новосибирск, 1973. С. 44-51. [Corev В. V. То unsteady creep calculation of a flexible rod made of material with different characteristics for tension and compression / In: Dynamics of Continuous Media: Collected scientific papers. 14. Novosibirsk, 1973. Pp. 44-51].
9. Леллеп Я. Установившаяся ползучесть балок в случае материалов с разными характеристиками на растяжение и сжатие // Тр. Тартус. ун-та (Tartu Ulikooli toimetised), 1975. Т. 355, №15. С. 245-252. [Lellep Ja. Steady-state creep of beams in case of materials with different characteristics at tension and compression// Tr. Tartus. Un-ta., 1975. Vol. 355, no. 15. Pp. 245-252].
10. Еремин Ю. А., Кайдалова Л. В., Радченко В. П. Исследование ползучести балок на основе аналогии структуры уравнения состояния материалов и элементов конструкций // Машиноведение, 1983. №2. С. 67-74. [Eremin Yu. A., Kaydalova L. V., Radchenko V. P. The study of creep of beams based on the analogy of the structure of the equation of state of materials and structural elements // Mashinovedenie, 1983. no. 2. Pp. 67-74].
11. Локощенко А. М., Печенина H. E., Шестериков С. А. Долговечность цилиндрического бруса при чистом изгибе // Изв. вуз. Машиностр., 1988. № 9. С. 9-13. [Lokoshenko А. М., Pechenina N. Е., Shesterikov S. A. The durability of a cylindrical rod under pure bending // Izv. vuz. Mashinostr., 1988. no. 9. Pp. 9-13].
12. Басов М. В. Применение структурной модели стержневого типа для решения задачи изгиба и разрушения балки в условиях ползучести/ В сб.: Актуальные проблемы современной науки. Естественные науки: Тр. 3-й Междунар. конф. молодых учёных и студентов. Часть 1. Математика. Механика (30 сентября - 2 октября 2002 г.). Самара: СамГТУ, 2002. С. 47-48. [Basov М. V. Application of a structural model of rod type to solving the problem of bending and fracture of a beam under creep conditions / In: Actual Problems of Modem Science. Natural Sciences. Part 1. Samara: SamGTU, 2002. Pp. 47-48].
13. Андреева E.A., Агафонов А. А. Решение краевой задачи изгиба балки из пластически разупрочняющегося материала на основании структурной модели / В сб.: Актуальные проблемы современной науки. Естественные науки: Тр. Междунар. форума молодых учёных (9-й Междунар. конф.). Часть 1-3. Математика. Математическое моделирование. Механика (20-23 ноября 2008 г.). Самара: СамГТУ, 2008. С. 168-178. [Andreeva Е. A., Agafonov A. A. The solution of the boundary value problem of rod bending made of plastic softening material on the basis of the structural model / In: Actual Problems of Modern Science. Natural Sciences. Part 1-3. Samara: SamGTU, 2008. Pp. 168-178].
14. Шестериков С. А., Юмашева М. А. Конкретизация уравнения состояния в теории ползучести// Изв. АН СССР. Мех. те. тела, 1984. №1. С. 86-91. [Shesterikov S.A., Yumasheva М. A. Specification of Equation of State in Creep Theory// Izv. AN SSSR. Mekh. Tverd. Tela, 1984. no. 1. Pp. 86-91].
15. Шестериков С.А., Юмашева М. А. Вариант уравнения состояния при ползучести и его приложения/ В сб.: Вопросы долговременной прочности энергетического оборудования/ Труды ЦКТИ. Вып. 246. Ленинград: НПО ЦКТИ, 1988. С. 74-77. [Shesterikov S.A., Yumasheva М. A. Variant of state equations for creep and its
applications / In: Voprosy dolgovremennoiy prochnosti energeticheskogo oborudovaniya/ Trudy TsKTI. 246. Leningrad: NPO TsKTI, 1988. Pp. 74-77].
Поступила в редакцию 07/1/2012; в окончательном варианте — 22/1/2012.
MSC: 74К10; 74С99
PURE BENDING OF BEAMS MADE OF MULTIMODULAR BEHAVIOR MATERIAL UNDER CREEP CONDITIONS
A. M. Lokoschenko, K. A. Agakhi, L. V. Fomin
М. V. Lomonosov Moscow State University,
Research Institute of Mechanics,
1, Michurinsky pr., Moscow, 119192, Russia.
E-mails: [email protected], [email protected], lef [email protected]
The problem of pure bending of beams with rectangular and circular cross-section is solved at steady-state creep, taking into account the different characteristics of the material in tension and compression. Constitutive relations for steady-state creep are taken in the form, of fractional functions. A comparison of solutions of problems for beams with square and circular cross-section with the same values of the axial moments of inertia is given.
Key words: beams, bending, steady-state creep, multimodulus behavior of material, fractional model, square cross-section, circular cross-section.
Original article submitted 07/1/2012; revision submitted 22/1/2012.
Alexander M. Lokoschenko (Dr. Sci. (Phys. & Math.)), Deputy Director, Research Institute of Mechanics. Kamilla A. Agakhi (Ph. D. (Phys. & Math.)), Head of Laboratory, Lab. of Strength and Creep at High Temperatures, Research Institute of Mechanics. Leonid V. Fomin, Leading Engineer, Research Institute of Mechanics.