ISSN 1992-6502 (Print) 2016. Т. 20, № 2 (72). С. 17-22
Ъыьмт QjrAQnQj
ISSN 2225-2789 (Online) http://journal.ugatu.ac.ru
УДК 669.295:539.4
Прогнозирование усталостной прочности стержня с выточкой из наноструктурного титанового сплава
при симметричном циклическом изгибе В. С. Жернаков1, Т. Н. Мардимасова2, В. Р. Акбашев3
1 zhvs@mail.rb.ru, 2 toma_mard@mail.ru, 3 vadja92@mail.ru
ФГБОУ ВО «Уфимский государственный авиационный технический университет» (УГАТУ)
Поступила в редакцию 20.06.2016
Аннотация. Проведено исследование усталостной прочности элементов конструкций из объемных наноструктурных материалов при многоцикловом нагружении. На основании описанной методики рассчитана усталостная прочность стержней с выточкой, выполненных из наноструктурного титанового сплава Ti-6Al-4V, при циклическом изгибе.
Ключевые слова: усталость, усталостная прочность, наноструктурный сплав, кривая Веллера.
Важным требованием для ответственных деталей в авиации и двигателестроении является высокая усталостная прочность. Современные методы аналитического расчета усталостной прочности конструкций сложной формы характеризуются высокой трудоемкостью. Одним из решений данного вопроса является применение численного расчета, что позволяет уменьшить трудоемкость и повысить производительность проектных работ [1].
В данной работе рассмотрена методика численного расчета усталостной прочности элементов конструкций. Приведены результаты исследования предела выносливости стержня с выточкой.
При расчете многоцикловой усталостной прочности конструкций процесс нагружения происходит при упругом деформировании, т.е. зависимость между напряжениями о, деформациями в, в теле конструкции является линейной.
Методика расчета усталостной прочности элементов конструкций при многоцикловом нагружении представляет собой ряд последовательных этапов.
На первом этапе, с целью определения полей напряжений огу и деформаций ву, выполняется решение статической задачи при максимальных значениях внешних нагрузок.
Исходными данными для расчета являются геометрические характеристики модели элемента конструкции, сеточный аналог, упругие константы материала (Е, ц), а также кривая деформирования материала.
На втором этапе для оценки усталостной прочности конструкции рассчитанное объемное
напряженное состояние сводится к эквивалентному одноосному опр. Такое приведение основано на критериях прочности в условиях статического нагружения [2]. Ниже представлены выражения, используемые для приведения амплитудных и средних напряжений цикла нагружения [3, 4]:
гипотеза пластичности - наибольших касательных напряжений Треска-Сен-Венана:
ъПа =ст1а"ст3а
°m ~1m 3 m
(1)
гипотеза пластичности Губера-Мизеса-Генки -потенциальной энергии формоизменения:
^ (°1a -a2a) +(°2a-°3a) +(°3a-°1a)
(2)
,2
,2
ат=^[а1ш-а2ш ] +[а2ш-а3ш ] +|°3ш -а1ш ]
где 01а, 02а, оза, 01т, 01т, озт,- главные напряжения тензоров амплитудных и средних значений напряжений.
Отметим, что при использовании соотношений (1) и (2) приведенное среднее напряжений цикла принимает только положительные значения. В случае необходимости учета знака напряжений применяется следующее выражение:
°ар (°1 а "°2 а)2 +(°2 а "°3 а)2 +(°3 а "°1а)2 ^(ча)
m (°1m "°2m)2 +(°2m "°3m)2+(°3m "°1m)2 xsign(°1m)
2
где sign (x
(x )=
1, x > 0 -1, x < 0
- функция знака
числа.
Для сведение асимметричного цикла напряжений к эквивалентному симметричному циклу используются соотношения для амплитуды эквивалентного симметричного цик-[5] на основании зависи-
ла напряжений □
мости Гудмана:
7пр a
а
пр "т
аПЧ
,атР^
(4)
0,
апр <0 m
где опч -предел прочности материала.
Таким образом, результатом приведения ха-
нр нр
симметричному являются следующие характе-
ристики: а t
экв. п ,ат =0 .
Для учета конструктивно-технологических факторов (концентрация напряжений, масштабный фактор, методы упрочнения деталей) применяется коэффициент пересчета предела выносливости К/. Коэффициент К/ - это отношение предела выносливости объекта О-1Д к пределу выносливости стандартного образца О-1 при одинаковой асимметрии цикла.
Таким образом, предел выносливости объекта исследований определяется следующим выражением:
а-1Д = Kfа-1 •
(5)
Выражение для амплитуды цикла напряжений при симметричном цикле с учетом характеристик имеет вид [6, 7]:
ааД = aaj Kf •
(6)
При асимметричном нагружении характеристики цикла определяются по уравнениям
[7]:
1° аД =° а/К/ [ ° тД =° т
(7)
В соответствии с ГОСТ 25.504-82 [8] коэффициент снижения предела выносливости определяется по формуле:
Kf = ч
K.
а
Kda KFa
—1
(8)
где Ка - эффективный коэффициент концентрации напряжений, К^ - коэффициент влияния абсолютных размеров поперечного сечения, К - коэффициент влияния шероховатости поверхности, К - коэффициент влияния поверхностного упрочнения.
Результатом вышеуказанных преобразова-
ний являются следующие
ки: а:
шахД"
тщД
экв. аД ' ° тД
„экв. п
характеристи-R= -1.
Следующим этапом расчета является вычисление в каждой точке тела значение числа циклов N при которых происходит усталостное разрушение.
При этом значение N определяется с помощью заданной кривой усталости материала для симметричного цикла изменения нагружения.
Кривая усталости материала описывает зависимость уровня амплитудных напряжений от числа циклов до разрушения. Одной из характеристик цикла является коэффициент асимметрии, который определяется следующим выражением:
R = а . а
min / шах
(9)
Кривая задается набором точек
°ah, Nh
h = 1, q
, при этом а ,< а , ., h = 1, q-1 ah ah +1
Значения Оан, N принимаются исходя из справочных данных, основанных на результатах проведенных экспериментов [7].
Описание кривой усталости на интервалах
между указанными точками И = 1, q проводится по следующим соотношениям [7]: линейная
N =
(Nh - Nh+1) <
+ аah+1Nh - l3ahNh+1
аah+1 - аah полулогарифмическая
lg N -
(lg Nh - lg Nh+1 К + 0ah+1 "8 Nh - аah "g Nh
0ah+1 - °h
и двойная логарифмическая интерполяция
lg N = -
lNh - lgNh+1 )lg°a + 'g°ah+1 "8Nh - '80ah >8Nh+1
lg 0ah+1 - lg0ah
(10)
(11)
,(12)
где Оа - значение амплитуды напряжения, для которого надо найти долговечность N.
В том случае, если Оа < Оа1, полагается, что долговечность N= N6, где N6 - база эксперимента [7].
1
а
a
экв.
a
a
(
При задании кривой усталости необходимо учитывать, какой вид напряженного состояния возникает в теле исследуемого объекта (одноосное напряженное состояние, либо чистый сдвиг) и, исходя из этого, подбирать необходимую кривую, полученную при эксперименте с аналогичным характером напряженного состояния.
Для определения вида напряженного состояния используется параметр в, который определяется выражением:
в =
V °3,
|°1 II °1|
Р1Л °1
(13)
где 01, <52, оз - главные напряжения тензора напряжений в некоторой точке, полученного в результате статического расчета. Введенный параметр принимает значения из отрезка [-1; 1]. В зависимости от значения параметра в напряженное состояние соответствует одному из следующих видов [7]:
• в = -1 - чистый сдвиг (используются кривые усталости, полученные при циклическом кручении);
• в = 0 - одноосное напряженное состояние (используются кривые усталости, полученные при циклическом изгибе или растяжении -сжатии в зависимости от вида распределения напряжений);
• в = 1 - чистое двухосное растяжение (сжатие) (используются кривые усталости, полученные при растяжении (сжатии)).
Основными характеристиками усталостной прочности конструкции являются коэффициент запаса по долговечности и коэффициент запаса по амплитудам напряжений
Согласно ГОСТ 23.207-78 [8], коэффициент запаса по долговечности определяется соотношением базовой долговечности Ы к текущему значению параметра N [7]:
N п5Р
N
(14)
где N = 10' ;
Коэффициент запаса по амплитудам напря-а
жений п определяется отношением амплитуды напряжения <а1 для стандартного образца, при базовом количестве циклов Ы, к текущему
экв.
значению приведенного параметра аа [19]:
-А
(15)
пБР = аа 1 / ааД .
Условием усталостного разрушения элементов конструкции является превышения приведенного амплитудного напряжения значения
Оаэкв при данном числе циклов для стандартного
образца: аэккв ) а 1 .
аД а 1
В данной работе проведен расчет усталостной прочности стержней с выточкой при симметричном циклическом изгибе. Численное решение статической задачи было достигнуто с помощью метода конечных элементов, реализованного с использованием программно-вычислительного комплекса ЛЫБУБ 14.5 [9].
Материалом стержней являлся нанострук-турный титановый сплав Ti-6Al-4V (£=0,83*105 МПа, v=0,32, ог=1250 МПа, опч=1370 МПа). Геометрические размеры стержня с выточкой представлены на рис. 1.
б
Рис. 1. Геометрические размеры стержня (а) и выточки (б)
Расчетная схема представлена на рис. 2 и представляет собой схему изгиба с кручением. В точке А стержень имеет жесткое закрепление. В точке В прикладывалась нагрузка изменяющаяся по симметричному циклу Я= —1.
Рис. 2. Расчетная схема для стержня с выточкой при изгибе
Приведение объемного напряженного состояния к эквивалентному одноосному
Для приведения рассчитанного объемного напряженного состояния к эквивалентному одноосному, использовалась гипотеза пластичности Губера-Мизеса-Генки - потенциальной
а
энергии формоизменения: (Equivalent (VonMises)).
Сведение асимметричного цикла напряжений к эквивалентному симметричному
Так как значение нагрузки изменялось по симметричному циклу, то учет асимметрии цикла не требовался.
Корректировка значений эквивалентных амплитудных напряжений для учета конструктивно-технологических факторов
Учет конструктивно-технологических факторов проводился при помощи коэффициента пересчета предела выносливости Kf , определенным по следующей формуле:
Kf =ч
к„
1
Kda KFa
—1
(16)
Чувствительность материала к масштабному фактору и местным напряжениям вычисляли при помощи выражения (5.20):
K„
2а„
K
da
1+1 88,3
G Га L
(17)
Градиент напряжений. Увеличение О снижает чувствительность материала к местным напряжениям. Градиент напряжений О , подобно теоретическому коэффициенту концентрации напряжений, зависит от формы тела и условий нагружения, и определяется законом изменения напряжений в окрестности очага концентрации:
ё = ^ + 1 = + = 27Д4, (18)
г с ОД 2,7
где й - диаметр стержня в зоне концентратора напряжений, г - радиус скругления на дне выточки, ф рассчитывалось по формуле:
Ф = = ^^=0.32. (19)
4iiz
ИГ \ IX
— 2
Протяженность концентратора Ь:
Ь = л*й= л*2.7 = 8.48 мм, где й - диаметр стержня на дне выточки. Учитывая вышеизложенное, вычислено
(20)
(21)
где ао - теоретический коэффициент концентрации напряжений, V о - коэффициент чувствительности материала к концентрации напряжений и масштабному фактору.
Таким образом, определившись со значениями необходимых коэффициентов, проводился расчет коэффициента пересчета предела вынос-
Kf= ч
K„
V Kda KFa
= 1
3,43+y-l I = 0,29 , (22)
где Kfc - коэффициент влияния шероховатости поверхности. Рабочую поверхность стержня подвергали шлифовке и механической полировке с обеспечением шероховатости Ra0,63 мкм, что соответствует Kfc= 1, K— коэффициент влияния поверхностного упрочнения. Так как стержень не подвергался никаким технологическим операциям по упрочнению, то K= 1.
Расчет числа циклов N, при которых происходит усталостное разрушение конструкции
На рис. 3 приведена кривая усталости для гладких образцов из наноструктурного титанового сплава Ti-6Al-4V.
Рис. 3. Кривая усталости для наноструктурного титанового сплава Ti-6Al-4V
По завершению подготовки исходных данных и решения статической задачи проводился расчет усталостной прочности.
ливости
Kt
Рис. 4. Поля осевых напряжений а? и деформаций е? в рабочей части стержня, возникающих при изгибе
В качестве примера на рис. 4 приведены поля напряжений о2 и деформаций 82 в стержне с выточкой, возникающих при изгибе силой Кщ^ 20,5Н.
Максимальные растягивающие напряжения возникают на внешних слоях в зоне концентратора = 372 МПа, а максимальные сжимающие -на внутренних и равны о2 = 372 МПа. Аналогичный характер наблюдается и по деформациям, при этом на внешних слоях 82 составляет 0,000117, а на внутренних - 0,000118.
Рис. 5. Кривая усталости для стержня с выточкой при симметричном изгибе
На рис. 5 приведены кривые усталости для стержня с выточкой из наноструктурного титанового сплава Ti-6Al-4V,1 полученные численным решением и экспериментальным путем.
ВЫВОДЫ
Анализ результатов численного расчета усталостной прочности стержней с выточкой из наноструктурного титанового сплава Ti-6Al-4V выявил, что предел выносливости о-1д соответствует значению а^д = 222 МПа. Экспериментальные исследования усталостной прочности стержней показали, что а^д достигает уровня
220 МПа. Таким образом, отклонение расчетных данных от экспериментальных составляет всего 0,9%. Однако, при числе циклов Ы= 104 разрушение детали при численном расчете происходит при уровне напряжений цикла в 252 МПа, а при эксперименте - 325 МПа. Отклонение составляет 22,5%. Таким образом, заметим, что при приближении значения амплитуды напряжений к малоцикловой области, разница между результатами численного решения и полученными экспериментальным путем увеличивается. Однако, в пределах многоцикловой области (Ы > 105), отклонение не превышает 8%.
СПИСОК ЛИТЕРАТУРЫ
1. Гречищев Е. С., Ильяшенко А. А. Соединения с натягом: Расчеты, проектирование, изготовление. М.: Машиностроение, 1981. 247 с. [E. S. Grechishev, A. A. Ilyashenko. Interference fit: Calculations, design, manufacture. MOSCOW: Mashinostroenie, 1981. 247 p. ]
2. Допуски и посадки: справочник. В 2 ч./ В. Д. Мягков, М. А. Палей., 6 изд., перераб. и доп. Л.: Машиностроение. Ленингр. отд-ние, 1982. [Tolerance and fit: Manual./ V. D. Myagkov, M. A. Palei., 6 edit. Leningrad: Mashinostroenie, 1982].
3. Зенкевич О. С.. Метод конечных элементов. М.: Мир, 1975. 540 с. [O. S. Zenkevich. Finite element method. Moscow: Mir, 1975. 540 p. ]
4. Жернаков В. С. Сопротивление материалов - механика материалов и конструкций. Уфа: УГАТУ. 2012. 494 с. [ V. S. Zhernakov. Strength of materials - Mechanics of materials and structures. Ufa: UGATU, 2012. 494 p. ]
5. Жернаков В. С. Влияние напряженно-деформированного состояния деталей из объемных наномате-риалов на их усталостную прочность /Семенова И. П., Ермоленко А. Н.//Вестник УГАТУ. 2009. Т. 12, № 2(31). С. 62-68. [ V. S. Zhernakov, I. P. Semenova, A. N. Ermolenko," Effect of stress-strain state of parts of bulk nanomaterials on their fatigue strength" (in Russian), in Vestnik UGATU, vol. 12, no. 2 (31), pp. 62-68, 2009. ]
ОБ АВТОРАХ
Жернаков Владимир Сергеевич, зав. каф. сопротивления материалов, д-р техн. наук, проф., засл. деят. науки РФ. Чл.-корр. АН РБ. Автор 8 монографий, 2 учебников и 8 учебных пособий*, более 250 статей, 90 авторских свидетельств и патентов, им подготовлено 12 канд. и 8 д-ров техн. наук. Область науч. интересов - решение фундаментальных и прикладных задач теории упругости, пластичности, ползучести и мех. разрушения применительно к основным элементам авиац. конструкций. Мардимасова Тамара Николаевна, канд. техн. наук, доц. Автор 5 учебных пособий, свыше 100 публикаций, 4 авторских свидетельств и патентов. Область науч. интересов -решение задач механики процессов упруго-вязкого пластического деформирования конструкций. Акбашев Вадим Ринатович, асп. каф. СМ, маг. техн. и тех-нол. (УГАТУ, 2015).
METADATA
Title: Prediction of fatigue strength rod with notching made of nanostructured titanium alloy under symmetric cyclic bending. Authors: Zhernakov V.S., Mardimasova T.N., Akbashev V.R.1
Email: 1 vadja92@mail.ru. Language: Russian.
Source: Vestnik UGATU (scientific journal of Ufa State Aviation Technical University), vol. 20, no. 2 (72), pp. 17-22, 2016. ISSN 2225-2789 (Online), ISSN 1992-6502 (Print). Abstract: A study of the fatigue strength of structural elements of the bulk nanostructured materials with highcyclic loading. Based on the described method is designed fatigue strength recessed rods made of nanostructured titanium alloy Ti-6Al-4V, under cyclic bending. Key words: fatigue, fatigue strength, nanostructured alloy, Wohler curve.
About authors:
Zhernakov Vladimir Sergeyevich, Head of the Department of Strength of Materials, Doctor of Technical Sciences, Professor, Honored Scientist of the Russian Federation. Corresponding Member of the Academy of Sciences of RB. Author of 8 monographs, 2 textbooks and teaching aids 8, more than 250 arti-cles,90 patents, they trained 12 candidates and 8 doctors of technical sciences. Research interests - the solution of fundamental and applied problems, the theory of elasticity, plasticity, creeping and fracture mechanics as applied to the basic elements of aircraft structures.
Mardimasova Tamara Nikolaevna, candidate of Technical Sciences., associate professor. Author of 5 textbooks and over 100 Publications, 4 patents. Research interests - the solution of problems of mechanics processes viscoelastic plastic deformation structures.
Akbashev Vadim Rinatovich, graduate student Department of Strength of Materials USATY. Master of Engineering and Technology (USATU 2015).