Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2013. № 3 (32). С. 77-86
УДК 539.376; 539.434
НАПРЯЖЁННО-ДЕФОРМИРОВАННОЕ СОСТОЯНИЕ ЦИЛИНДРИЧЕСКОГО ОБРАЗЦА ИЗ СПЛАВА Д16Т В УСЛОВИЯХ ОСЕВОГО РАСТЯЖЕНИЯ И КРУЧЕНИЯ ПРИ ПОЛЗУЧЕСТИ
В. П. Радченко, В. В. Цветков
Самарский государственный технический университет, Россия, 443100, Самара, ул. Молодогвардейская, 244.
E-mails: [email protected], [email protected]
Предложена методика 'расчёта 'реологического деформирования и разрушения цилиндрических образцов в условиях ползучести для трёх видов напряжённого состояния: растяжение, чистое кручение, совместное действие растягивающей нагрузки и крутящего момента. Методика базируется на теории ползучести и длительной прочности энергетического типа. Выполнены расчёты для всех трёх видов напряжённого состояния для сплошных и полых цилиндрических образцов из алюминиевого сплава Д16Т при температуре 250 С. Проведено сравнение расчётных данных с соответствующими экспериментальными данными при каждом виде напряжённого состояния. Наблюдается соответствие расчётных и экспериментальных значений. Показано, что приложение крутящего момента к растягиваемому осевой нагрузкой образцу интенсифицирует процесс ползучести и приводит к снижению показателей длительной прочности. Наблюдается существенное перераспределение осевых и касательных напряжений по радиусу в зависимости от времени. Приведены оценки погрешностей отклонения расчётных данных от экспериментальных значений.
Ключевые слова: ползучесть, длительная прочность, цилиндрический образец, растяжение, кручение, растяжение с кручением, метод расчёта, экспериментальные данные, сплав Д16Т.
Ужесточение температурно-силовых режимов нагружения, снижение материалоёмкости и использование всех предельных возможностей материала сопротивляться неупругому деформированию и разрушению — центральная задача машиностроения, авиастроения и других отраслей промышленности. Это требует развития методов расчёта элементов конструкций в условиях ползучести с учётом разупрочнения для различных видов напряжённого состояния. Целью данной работы является исследование напряжённо-деформированного состояния сплошного цилиндрического образца в условиях ползучести при совместном действии кручения и растяжения на основе энергетического варианта реологического деформирования и разрушения [1], а также сравнительный анализ данных расчёта по двум различным моделям, представленным в работах [1,2].
Для экспериментальной проверки всех последующих исследований использовались экспериментальные данные по ползучести и длительной прочности для сплава Д16Т при T = 250 С, приведённые в работе [2] для случаев
Владимир Павлович Радченко (д.ф.-м.н., проф.), заведующий кафедрой, каф. прикладной математики и информатики. Виталий Владимирович Цветков, студент, каф. прикладной математики и информатики.
растяжения, кручения и совместного растяжения и кручения (первая стадия ползучести у данного материала отсутствует).
1. Одноосная реологическая модель ползучести и длительной прочности.
В качестве основной модели используется энергетический вариант теории ползучести и длительной прочности [1]. Уравнения состояния этой модели без учёта первой стадии ползучести имеют вид
p = cam, a = ао(1 + w), w = аар.
(1)
Здесь р — деформация ползучести; ш — параметр повреждённости; а и а0 — истинное и номинальное напряжения соответственно; с, т, а — подлежащие определению параметры модели.
В качестве критерия разрушения принимается величина
Q(t) = Г
о
adp
Акр
(2)
t, час
Рис. 1. Экспериментальные (точки) и расчётные кривые ползучести при растяжении для сплава Д16Т при температуре 250 °С: 1 — сто = 68,7 МПа; 2 — сто = 73,6 МПа; 3 — сто = 78,48 МПа; 4 — ст0 = 88,3 МПа; 5 — ст0 = = 98,1 МПа
где Акр — критическая величина работы истинных напряжений на деформации ползучести, по достижении которой в момент времени t = t* происходит разрушение материала и Q(t*) = 1.
В общем случае параметры а и Акр имеют степенные аппроксимации а = = aiama, Акр = AiamA, где аь ша, А1 и тл — параметры модели.
В модели (1) и (2) в отличие от модели работы [2] во всех соотношениях используется не номинальное напряжение ао, а истинное a.
Идентификация параметров c, т, а1, та, А1, тл хорошо известна [1] и осуществлялась по следующей схеме.
1. С использованием прямолинейных участков кривых ползучести (на начальной стадии w ~ 0 и a = ао) методом наименьших квадратов при пяти ао = const определялись параметры степенной аппроксимации c и m первого соотношения (1).
2. Интегрирование первого соотношения (1) при ао = const и начальных условиях p(0) = 0, w(0) = 0 даёт зависимость деформации ползучести от ао:
P(t) = -
1
таао
ln 11 - arncam+1t|
(3)
и величину работы истинного напряжения на деформации ползучести:
Г 1 A(t) = adp = -оа
(1 - arncam+1t)-1/m - 1
(4)
Параметр а для каждой кривой ползучести, представленной на рис. 1, определялся численно из условия прохождения функции (3) через последнюю точку (t*,p*). При этом уравнение (3) разрешалось относительно а при заданных c, m, ст0, t*, р* = p(t*) и выполнялась аппроксимация параметра а степенной зависимостью а = а^т™.
3. Далее для каждой кривой при ст0 = const по формуле (4) рассчитывалась критическая величина работы разрушения Акр = A(t*) и известных а, m, c и строилась степенная зависимость Акр = Aict0™a.
В результате для модели (1), (2) получены следующие параметры: c = = 2,4 ■ 10-15 (МПа)-т, m = 5,405, а1 = 0,91 (МПа)"та-1, ma = -0,43, Ai = 19,72 (МПа)-1, тл = 0.
На рис. 1 сплошными линиями показаны результаты расчёта по модели (1), (2), штриховыми линиями — расчёт по методике работы [2], точками — экспериментальные данные [2] по ползучести сплава Д16Т при чистом растяжении при T = 250 С. Анализ информации, представленной на рис. 1, позволяет сделать вывод о сопоставимости экспериментальных данных и результатов расчёта по моделям, представленным в работах [1,2].
2. Методика расчёта при чистом кручении. Согласно [1], уравнения состояния (1), (2) при введении скалярного параметра повреждённости на случай сложного напряженного состояния для изотропной ползучести обобщаются следующим образом:
Pij = 2cSm-i (CTij - 1 Okk j, aij = ст0(1 + w), w = а^стур,. (5)
Здесь pj —тензор деформации ползучести; ст,, ст0, —соответственно компоненты истинного и номинального тензоров напряжений; S, So — соответственно интенсивности тензоров истинных и номинальных напряжений; c, m — константы модели; а^0) задаётся степенной аппроксимацией а^0) = а^О"".
Критерий разрушения энергетического типа для сложного напряжённого состояния имеет вид
Q(t*)= Г' =1, (6)
J0 Акр
где Акр задаётся аппроксимацией Акр = A1 S0™A, А1 и m-л — параметры модели.
Все параметры модели (5), (6) для случая изотропной ползучести могут быть определены по результатам одноосных испытаний на растяжение (их значения приведены выше).
Конкретизируем соотношения (5), (6) для частного случая — чистого кручения. Для данного вида напряжённого состояния в работе [2] приведены результаты испытаний, выполненных на тонкостенных образцах из рассматриваемого сплава Д16Т. Соответствующие экспериментальные данные представлены на рис. 2 точками. В [2] отмечается, что сплав Д16Т обладает анизотропией свойств ползучести при растяжении и кручении. Поэтому все параметры модели (5), (6) в условиях кручения будут отличаться от параметров, полученных при растяжении.
Частный случай модели (5), (6) при чистом кручении имеет вид
Yp = ^3с*(^3т )m Т = То(1 + w), w = a*TYp,
(7)
где 7Р — угол сдвига в ползучести; т, то — истинное и номинальное касательные напряжения; с*, ш*, а* — параметры.
Выражение для критерия разрушения при кручении имеет вид
г *
600 t, чах:
о
Td7p Акр
= 1,
(8
Рис. 2. Экспериментальные и расчётные кривые ползучести для сплава Д16Т на кручение при температуре 250°С: 1 — то = = 33,94 МПа; 2 — т0 = 36,3 МПа; 3 — т0 = = 39,7 МПа; 4 — то = 45,3 МПа;
где t* — время разрушения при кручении, АКр — критическая величина работы при кручении.
В силу анизотропии свойств материала Д16Т при растяжении и кручении параметры модели для кручения с*,
m*, a* = a1S0 a, AKP = A (S0 = v3t0) определяются по результатам
испытаний в координатах «yp - t» при т = const.
Решение (7) при т0 = const и начальных условиях Yp(0) = 0, w(0) = 0 имеет вид
Y
p —
1
T0m*a*
ln
1 - a*m*c*^^3T0)m*+1t
(9)
а для работы получаем
1
A(t)=j,Td7P=f1 -
a*m*c*
m*+1, \-1/m*
1
Для дальнейших исследований и конкретизации параметров модели (7) используются экспериментальные данные по ползучести при чистом кручении в координатах «yp - t» при Т0 = const для материала Д16Т при T = = 250С, приведённые в [2] и представленные точками на рис. 2.
1. С использованием прямолинейных участков кривых ползучести методом наименьших квадратов при четырёх значениях Т0 определялись параметры с* и m* степенной аппроксимации первого соотношения (7) (здесь можно считать, что т = Т0).
2. Далее при известных значениях m* и с* и условии прохождения функции (9) через точку разрушения (t*,Y*) для каждого значения т0 определялся параметр a* и выполнялась аппроксимация a* зависимостью a* = a*(^3T0)ma.
3. После определения параметров с*, m* и a* для каждого т0 = const определялась работа разрушения АКр = A(t*) и строилась степенная зависимость АКр = А*(\/3т0)тА.
t
В итоге получены следующие параметры: с* = 2,1 ■ 10 15 (МПа) m*, m* = = 5,808, а* = 118,3 (МПа)-^-1, m% = -1,726, A* = 214,6 (МПа)-тА-1, mA = -0,668.
На рис. 2 сплошными линиями показаны результаты расчёта 7p = 7p(t) по модели (7), (8), точками — экспериментальные данные [2], штриховыми линиями — расчётные данные по модели [2].
3. Методика решения краевой задачи о кручении стержня кругового поперечного сечения при ползучести. Рассмотрим стержень кругового поперечного сечения радиуса R, к которому приложен крутящий момент M = M(t).
Уравнение равновесия элемента стержня имеет вид
,-к
Г R
2W то(г, t)r2dr = M(t), (10)
J о
где г € [0, Л] — расстояние от элемента площади сечения до нейтральной оси, то (г, ¿) —номинальное напряжение.
Уравнение совместности деформаций в соответствии с гипотезой прямых радиусов принимает вид
7(г,*) = г0(£), 7 (г, ¿) = 7е (г, ¿) + 7 (г,*), (11)
где 7(г, ¿) — полный угол сдвига, 7е (г, ¿) = т0 (г, ¿)/С — угол сдвига в упругости, С = Е/(2(1+ —модуль сдвига материала, 7(г,¿) —угол сдвига в ползучести, #(£) —относительный угол закручивания.
При расчёте весь процесс деформирования разбивается на отрезки , ¿7+1] длиной Д^, j = 0,1, 2,..., в пределах которых все характеристики напряжённого состояния считаются постоянными. Тогда для j-того отрезка времени можно записать
Y(r, tj+i) = Y (r, tj) + ^G^ + A7j (r),
где
A7j(r) = 7(r, tj+i) - 7(r, j^ AT0j(r) = то(r, tj+i) - To(r j). (12) Из (10) и второго соотношения (12) получаем Г R
2W AToj(r)r2dr = M(tj+i) - M(j). (13)
о
Соотношения (11) дают
Atoj (r) G
+ Ayj (r) = rA%,
где Д07- = #(¿7+0-#(¿7). Умножая последнее равенство на 2пг2^г, интегрируя полученное в пределах от 0 до Л и используя (13), получаем
= (М(¿7+1) - М(¿7)) + Ё4 I Д7Р(г)г2^г. (14)
81
Легко видеть, что
Дтад(г) = С(тАв3 - 7))
(15)
Далее получаем все искомые характеристики по следующей схеме. Приращения неупругих деформаций при фиксированном г на отрезке [¿7 ,¿7+1] определяются по формулам (14), (15). При этом неупругие деформации и параметр повреждённости имеют нулевые начальные условия. Величины 7Р, то, ш, 9 при £ е [¿7,¿7+1] определяются по рекуррентному соотношению
Р (¿7 + 1)= Р (¿7 )+ДР, ,
где Р — любая из этих величин.
Расчёт ползучести стержня осуществляется до тех пор, пока выполняется условие
П(*) = Г
о
< 1.
в, ММ г
0,01
А?(\/3то)тА
Если при каком-либо значении радиуса стержня выполняется условие = = 1, то предполагается, что при £ = происходит разрушение образца.
В качестве примера был рассмотрен случай неупругого деформирования и разрушения цилиндрического образца кругового поперечного сечения с радиусом К = 10 мм из сплава Д16Т при температуре Т = 250 °С. Экспериментальные данные для указанного цилиндра приведены в [2]. Все параметры реологической модели (7), (8) приведены выше, а значения модуля Юнга и коэффициента Пуассона следующие: Е = 56000 МПа, ^ = 0,3.
На рис. 3 сплошными линиями представлены расчётные данные по модели (7), (8), (10)-(15), точками — экспериментальные значения погонного угла закручивания 9 = 9(£) для трёх величин крутящего момента М = {69,3; 77,05; 92,45} кН ■ мм. Здесь же штриховыми линиями показаны расчётные данные по модели из работы [2].
На рис. 4 и 5 представлено распределение номинальных и истинных касательных напряжений соответственно по радиусу стержня в различные моменты времени при М = 92,45 кН ■ мм.
Как следует из информации, представленной на рис. 3, наблюдается удовлетворительное соответствие расчётных и экспериментальных значений деформации ползучести в координатах «относительный угол закручивания 9 -
Нр^2 .....1..... 1./ : ♦ / /
■ * /
0 200 400 600 г, час
Рис. 3. Экспериментальные и расчётные значения погонного угла закручивания стержня из сплава Д16Т при температуре 250С: 1 — М = 69,3 кН • мм; 2-М = 77,05 кН • мм; 3 — М = 92,45 кН • мм
то • 10 , МПа
го • Ю-1, МПа
з2
4 ^^
1/ 5
0
6
р, мм
Рис. 4. Распределение номинального касательного напряжения по радиусу стержня в различные моменты времени при М = 92,45 кН • мм: 1 — Ь = 0 ч; 2 — Ь = 38 ч; 3 — Ь = 76 ч; 4 — Ь = 114 ч; 5 — Ь = 152 ч
2 4 6 8 Р>мм
Рис. 5. Распределение истинного касательного напряжения по радиусу стерж-различные моменты времени при
ня
M
в различные моменты = 92,45 кН • мм: 1 — t = 0 ч;2 — t = 38 ч; 3 — t =76 ч; 4 — t =114 ч; 5 — t =152 ч
время t» при различных значениях крутящего момента. Из рис. 4 и 5 следует, что при ползучести происходит существенное перераспределение номинального и истинного касательных напряжений по радиусу в зависимости от времени.
4. Решение краевой задачи о реологическом деформировании цилиндрического образца при совместном растяжении и кручении. Рассмотрим реологическую модель ползучести при совместном действии растяжения и кручения с учётом отмеченной анизотропии свойств для этих видов напряжённого состояния. Исходя из общей модели при сложном напряжённом состоянии и одномерных моделей частного вида (1) и (7) для случая анизотропной ползучести определяющие соотношения можно принять в следующем виде:
p = cS--1o, y = 3c*S—*-1т, о = oo(1 + ш),
т = то(1+ ш), So = ^о2 + 3Т0, S = So(1 + ш), (16)
ш = aop + a*TYp, а = a1(S0)-a, а* = а\ (S0)—,
где p — осевая деформация ползучести, yp — угол сдвига в ползучести, о и 00 — истинные и номинальные нормальные напряжения, т и T0 — истинные и номинальные касательные напряжения, S и S0 — интенсивности истинных и номинальных напряжений; все остальные параметры определяются моделями (1), (2) и (7), (8).
Соответственно изменится и критерий разрушения. Расчёт осуществляется до тех пор, пока выполняется условие
п(*)= i t^dp- + Г^-г < i. (17)
J0 A1(S0)- J0 A*(S0)— v 7
Очевидно, что модели (1), (2) и (7), (8) являются частными случаями модели (16), (17) при T0 = 0 и 00 = 0 соответственно.
Рассмотрим решение краевой задачи о ползучести полого толстостенного цилиндрического образца с внутренним радиусом R1 и внешним R2 при совместном действии растягивающей нагрузки Q и крутящего момента M.
Запишем уравнения равновесия и совместности деформаций с использованием гипотезы плоских сечений и прямых радиусов:
2п Г" сто(т, 1)тйт = Я(1), + р(т, I) = ео(Ь), (18)
то(т, г)
/Яг &
где т е ,К2] —расстояние от элемента площади до оси цилиндра, 9{Ь) — относительный угол закручивания, е0(г) — осевая деформация образца.
Решение задачи (18), (19) на основе реологической модели (16), (17) осуществляется численно «шагами» по времени аналогично случаю чистого изгиба.
В качестве примера рассмотрим результаты расчёта ползучести полого цилиндрического образца из сплава Д16Т (Т = 250 °С) при К1 = 5 мм, К2 = = 10 мм и различных комбинациях Q и М, а также сравним их с соответствующими экспериментальными данными, представленными в работе [2]. На рис. 6 и рис. 7 сплошными линиями показаны расчётные значения осевой деформации ползучести р = р(г) при различных значениях растягивающей нагрузки Q и крутящего момента М, на рис. 8 приведена кинетика перераспределения номинальных напряжений сто = сто(т) и касательных напряжений то = то(т) в различные моменты времени при Q = 18,48 кН, М = 47,2 кН ■ мм.
'Яг
Я
2п [ 2то(т,г)т2(1т = М(г), 7р(т,г) = тб(г), (19)
Рис. 6. Осевая деформация ползучести полого цилиндрического образца из сплава Д16Т (Т = 250 С) при осевой нагрузке ^ = 16,17 кН и различных значениях крутящего момента: 1 — М = о, 2 — М = = 28,9 кН • мм, 3 — М = 43,3 кН • мм, 4 — М = = 5о,6 кН^мм; точки — экспериментальные данные [2]; сплошные линии — расчёт по модели (16)—(19); штриховые линии — расчёт по модели [2]
Рис. 7. Осевая деформация ползучести полого цилиндрического образца из сплава Д16Т (Т = 250 С) при осевой нагрузке ^ = = 18,48 кН и различных значениях крутящего момента: 1 — М = о, 2 — М = 47,2 кН • • мм; точки — экспериментальные данные [2]; сплошные линии — расчёт по модели (16)—( 19); штриховые линии — расчёт по модели [2]
Рис. 8. Распределение номинальных напряжений а0 = 00 (г) (а) и номинальных касательных напряжений т0 = т0(г) (б) в различные моменты времени при ^ = 18,48 кН, М = 47,2 кН • мм: 1 — I = 0; 2 — I =100 ч; 3 — I = 200 ч; 4 — I = 250 ч
Q, кН M, кН • мм С ч tp, ч tp, ч [2]
16,17 0 760 825 827
16,17 28,9 416 453 430
16,17 43,3 303 255 323
16,17 50,6 223 187 247
18,48 0 365 371 368
18,48 47,2 143 132 151
Анализ этих данных позволяет сделать вывод, что приложение крутящего момента к нагруженному осевой нагрузкой цилиндрическому образцу существенно интенсифицирует процесс ползучести и приводит к перераспределению во времени не только касательных напряжений по радиусу, но и нормальных напряжений. С целью сопоставления результатов расчётов в таблице для каждой комбинации растягивающей осевой силы Q и крутящего момента M приведены экспериментальные t* и расчётные значения tp времени разрушения по модели (16)—(19) и по модели работы [2]. Ошибка отклонения расчётного времени разрушения от экспериментального вычислялась по формуле
N t* _t*
д = N Е ■ 100%, N i=1 i
где N = 6 — число вариантов расчёта, и для модели (16)-(19) составила величину 9,8%, а для модели, предложенной в [2], —8,4%. В целом наблюдается удовлетворительное соответствие расчётных данных по предложенной авторами методике экспериментальным данным и данным расчёта по модели работы [2].
БИБЛИОГРАФИЧЕСКИЙ СПИСОК
1. В. П. Радченко, Ю. А. Ерёмин, Реологическое деформирование и разрушение элементов конструкций. М.: Машиностроение-1, 2004. 264 с. [V. P. Radchenko, Yu. A. Eremin, Rheological Deformation and Fracture of Materials and Structural Elements. Moscow: Mashinostroenie-1, 2004. 264 pp.]
2. О. В. Соснин, Б. В. Горев, А. Ф. Никитенко, Энергетический вариант теории ползучести. Новосибирск: Институт гидродинамики СО РАН СССР, 1986. 95 с. [O. V. Sosnin,
B. V. Gorev, A. F. Nikitenko, Energy variant of creep theory. Novosibirsk: Inst. of Hydrodynamics, USSR Acad. of Sci., 1986. 95 pp.]
Поступила в редакцию 24/VI/2013; в окончательном варианте — 24/VII/2013.
MSC: 74D10
THE STRESS-STRAIN STATE OF CYLINDRICAL SAMPLE FROM ALLOY D16T UNDER AXIAL TENSION AND TORSION CREEP
V. P. Radchenko, V. V. Tsvetkov
Samara State Technical University,
244, Molodogvardeyskaya st., Samara, 443100, Russia.
E-mails: [email protected], [email protected]
The method for calculation of cylindrical sample rheological deformation and fracture in creep conditions for three types of stress state: tension, pure torsion, the combined effect of tensile load and torque is offered. The procedure is based on the theory of creep and creep rupture strength of energy type. Calculations are performed for all three types of stress state for solid and hollow cylinder tests from aluminum alloy D16T at 250 C. Comparison of the calculated data with the corresponding test data for each type of stress state is conducted. It shows the agreement between the calculated and experimental values. The intensification of creep and decreasing of creep rupture strength factors after the application of torque to the specimen under the axial tension are established. The substantial redistribution of axial and shear stresses along the radius depending on the time is observed. The estimates of errors of deviation of calculated data from the experimental values are given.
Key words: creep, creep rupture strength, cylindrical sample, tension, torsion, tension with torsion, the method of calculation, test data,, alloy D16T.
Original article submitted 24/VI/2013; revision submitted 24/VII/2013.
Vladimir P. Radchenko (Dr. Phys. & Math. Sci.), Head of Dept., Dept. of Applied Mathematics & Computer Science. Vitaliy V. Tsvetkov, Student, Dept. of Applied Mathematics & Computer Science.