Solution of the boundary-value problem of torsion for solid and hollow cylindrical specimens made of the Steel 45 and AMG-6M alloy under short-term steady-state creep conditions

Full Text

Введение. Повышение требований к прочности и надежности элементов конструкций, работающих в условиях повышенных температур, определяет необходимость исследования процесса ползучести при различных видах напряженного состояния для широкого класса материалов. Так, для цилиндрических изделий типичными испытаниями на ползучесть и длительную прочность являются испытания сплошных стержней при растяжении и кручении, а для трубчатых образцов - при комбинированном действии осевой силы, крутящего момента и внутреннего давления. При ползучести необходимо учитывать физическое разрушение материала вследствие ползучести - длительную прочность. Прогнозирование длительной прочности требует существенных экспериментальных затрат и связано с аппроксимацией экспериментальных данных при всех видах напряженного состояния. Многие задачи прогнозирования длительной прочности конструкций при сложном напряженном состоянии решаются на основе концепции эквивалентных напряжений [1-6]. Данный подход подробно изложен в монографии [7], однако его существенным недостатком является то, что за рамками исследования находится кинетика напряженно-деформированного состояния в процессе ползучести конструкций. Второй подход связан с решением краевых задач ползучести и длительной прочности на основе феноменологических уравнений состояния материала с использованием критериев разрушения. Проведение экспериментальных исследований на ползучесть при сложном напряженном состоянии связано со значительными техническими трудностями, и основная часть испытаний проводится, как правило, при одноосном растяжении цилиндрических образцов. Поэтому желательно иметь реологические модели, построение которых осуществляется лишь по одноосным кривым ползучести и которые бы адекватно описывали процесс ползучести при сложном напряженном состоянии. Примером такой реологической модели может служить энергетический вариант теории ползучести и длительной прочности, описанный в работе [8]. Теоретических методов решения краевых задач в теории ползучести при сложном напряженном состоянии, позволяющих в квадратурах получить зависимости для компонент тензоров напряжений и деформаций, на сегодняшний день относительно немного. В частности, в работе [9] с использованием модели ползучести [8] была решена краевая задача о реологическом деформировании и разрушении толстостенной трубы в условиях внутреннего давления и осевого растяжения и выполнена проверка адекватности метода расчета экспериментальным данным по длительной прочности толстостенных образцов из сталей 20, 12ХМФ и ЭИ694. Здесь же проведено сравнение и с результатами расчетов по концепции эквивалентных напряженных состояний. Для решения задачи о кручении цилиндрических конструкций в условиях ползучести в работе [10] предложен приближенный метод характеристической точки. В его основе лежит существование особой характеристической точки, в окрестности которой реализуется практически неизменное во времени напряженное состояние. При таком подходе все характеристики деформированного состояния определяются в характеристической точке, и необходимость решать краевую задачу отпадает. Этот подход успешно применялся 508 Решение краевой задачи о кручении сплошных и полых цилиндрических образцов . . . для решения разнообразных задач ползучести конструкций [11, 12]. В работе [13] разработан численный метод для расчета кинетики напряженно-деформированного состояния цилиндрических образцов при комбинированном растяжении и кручении. Однако проверка адекватности данного метода, выполненная в [13], носит частный характер, поскольку она была проведена лишь для одного материала (алюминиевого сплава Д16Т), при этом предполагалось, что материал проявляет анизотропию свойств ползучести при чистом растяжении и чистом сдвиге. В связи с вышеизложенным целью настоящей работы является дальнейшая проверка реологической модели [8] и методики расчета [13] с использованием экспериментальных данных по кратковременной установившейся ползучести стали 45 при температурах 725 и 740 ℃ [11] и сплава АМГ-6М при температуре 450 ℃ [14] в условиях одноосного растяжения и при кручении сплошных и полых стержней из этих материалов. 1. Методика решения краевой задачи о ползучести цилиндрического образца при кручении. Изложим методику решения краевой задачи о реологическом деформировании цилиндрического образца с радиусами R1 и R2 в условиях ползучести при кручении, описанную в работе [13]. Отметим, что данная методика применима как для сплошных, так и для трубчатых цилиндрических образцов (при этом для сплошных образцов полагаем R1 = 0). Рассмотрим цилиндрический образец, к торцу которого приложен крутящий момент M . Уравнение равновесия элемента стержня имеет вид R2 2π τ (r, t)r2 dr = M, (1) R1 где r ∈ [R1 , R2 ] - расстояние до оси стержня, τ - касательное напряжение. Уравнение совместности деформаций в соответствии с гипотезой прямых радиусов имеет вид γ(r, t) = rθ(t); γ(r, t) = γ e (r, t) + γ p (r, t), (2) где γ - полный угол сдвига, γ e (r, t) = τ (r, t)/G - упругая компонента деформации сдвига (G = E/ (2(1 + µ)) - модуль сдвига), γ p - сдвиговая деформация ползучести, θ - относительный угол закручивания. Соотношения (2) дают τ (r, t) + γ p (r, t) = rθ(t), G откуда получаем τ (r, t) = G rθ(t) - γ p (r, t) . (3) Подставляя соотношение (3) в (1) и проводя необходимые операции интегрирования, получаем формулу для нахождения угла закручивания: θ(t) = M 2π + GJ J R2 γ p (r, t)r2 dr, (4) R1 509 Ц в е т к о в В. В. где J - момент инерции относительно оси стержня: R2 J = 2π r3 dr = R1 π(R24 - R14 ) . 2 Вычислив по формуле (4) величину угла закручивания, можно определить распределение касательного напряжения по формуле (3). Таким образом, алгоритм решения краевой задачи о кручении следующий: сначала определяется начальное упругое напряженное состояние при t = 0: γ p (r, 0) = 0; τ (r, 0) = M r, J а затем определяются искомые величины в условиях ползучести (при t > 0) по схеме: (4) (3) t→ - γ p (r, t) --→ θ(t) --→ τ (r, t), (5) где γ p (r, t) рассчитывается по выбранной теории ползучести, а цифры над стрелками означают формулы, по которым вычисляется конкретная величина. 2. Реологическая модель. В соответствии со схемой (5) для решения краевой задачи необходимо уметь рассчитывать компоненту деформации ползучести γ p . Для этого можно использовать любую теорию ползучести, которая бы адекватно описывала экспериментальные данные. Для рассматриваемых материалов при испытаниях на кратковременную ползучесть присутствует лишь вторая стадия ползучести. Поэтому в качестве реологической модели можно использовать теорию установившейся ползучести, которая в одноосном случае имеет вид (6) p˙ = cσ m . Здесь p - осевая деформация ползучести; σ - осевое напряжение, c и m - параметры материала при данной температуре. Согласно [8], теория (6) на случай сложного напряженного состояния обобщается следующим образом: 1 3 p˙ij = cS m-1 σij - δij σll , 2 3 (i, j, l = 1, 2, 3), (7) где δij - символ Кронекера, S - интенсивность напряжений, по повторяющемуся индексу l осуществляется суммирование. Тогда для частного случая - чистого кручения - реологическая модель будет иметь вид √ √ γ˙ p = (c 3)(τ 3)m . Параметры c и m определяются по экспериментальным кривым ползучести в условиях одноосного растяжения при σ = const. Для каждого k-того испытания при напряжении σk = const определяется скорость деформации 510 Решение краевой задачи о кручении сплошных и полых цилиндрических образцов . . . ползучести p˙k (например, графически). Далее осуществляется аппроксимация полученных данных (σk , p˙k ) степенной зависимостью (6). 3. Оценивание параметров степенной зависимости на основе модифицированного метода наименьших квадратов. Изложим методику идентификации параметров c и m степенной аппроксимации (8) yˆ(x) = cxm в соответствии с «классическим» подходом, рекомендованным в [15]. Исходными данными для построения модели (8) служат табличные данные (xk , yk ) при k = 1, 2, . . . , N , где N - объем экспериментальных данных. Согласно методу наименьших квадратов (МНК), оценки параметров c и m должны определяться из условия минимизации функционала - суммы квадратов отклонений теоретических значений yˆ(xk ) от табличных значений yk : N (yk - yˆ(xk ))2 → min . Qост (c, m) = k=1 В этом случае система нормальных уравнений формируется на основе вычисления частных производных ∂Qост /∂c и ∂Qост /∂m. Нелинейность функциональной зависимости (8) приводит к нелинейности системы нормальных уравнений относительно параметров c и m, что существенно усложняет процедуру их вычисления. В связи с этим в большинстве расчетных методик применяется метод линеаризации. Логарифмируя обе части соотношения (8), получим ln yˆ(x) = ln c + m ln x. Обозначив uk = ln yk , vk = ln xk , u ˆ = ln yˆ, v = ln x, A = ln c, B = m, получим линейную по параметрам A и B модель, описывающую табличные данные (vk , uk ): (9) u ˆ(v) = A + Bv. Минимизируя для модели (9) «невязку» нейную систему нормальных уравнений вида N k=1 (uk -u ˆ(vk ))2 , получаем ли- ˆ N vk = N uk , N Aˆ + B k=1 k=1 N ˆ N v 2 = N vk uk , Aˆ k=1 vk + B k=1 k k=1 ˆ m ˆ из решения которой с учетом введенных обозначений cˆ = exp(A), ˆ = B находим оценки параметров исходной степенной зависимости (8): cˆ = exp m ˆ = N N k=1 N N N 2 k=1 ln xk - k=1 ln xk · k=1 2 N N 2 N k=1 ln xk - k=1 ln xk N N N k=1 ln xk ln yk - k=1 ln xk · k=1 ln yk 2 N N 2 N k=1 ln xk - k=1 ln xk ln yk · ( ( ) ) . ln xk ln yk ; (10) 511 Ц в е т к о в В. В. Однако оценки, полученные по формулам (10), оказываются достаточно грубыми, поскольку при формировании системы нормальных уравнений минимизировалась не величина N N 2 (yk - cxm k ) , (yk - yˆ(xk ))2 = k=1 k=1 а «невязка» N N 2 (ln yk - ln c - m ln xk )2 . (uk - u ˆ(vk )) = k=1 k=1 Несмотря на простоту этого вопроса, аппроксимацию зависимости (8) можно улучшить и тем самым существенно уменьшить величину погрешности отклонения расчетных данных от экспериментальных (и как следствие этого - повысить точность расчетной методики решения краевой задачи установившейся ползучести), воспользовавшись идеями работ [16, 17]. Представим табличные экспериментальные данные в виде (11) yk = cxm k + ek , где через ek обозначены остатки, которые носят случайный характер. Прологарифмируем обе части выражения (11): ln yk = ln (cxm + ek ) = ln cxm k 1+ ek cxm k = A + B ln xk - ln 1 - ek , yk где A = ln c, B = m. При малых значениях ek величина |ek /yk | < 1 и имеет место разложение в ряд Маклорена ln 1 - ek yk =- ek 1 ek - yk 2 yk 2 - 1 ek 3 yk 3 - ..., используя которое, из предыдущего выражения получим ln yk = A + B ln xk + где eˆk = ek + eˆk , yk e2k e3 + k2 + · · · = ek + O(e2k ) ≈ ek . 2yk 3yk Окончательно имеем yk ln yk = Ayk + Byk ln xk + eˆk . (12) Оценки коэффициентов модели (12) находим из условия минимизации функционала N N e2k k=1 512 N eˆ2k ≈ k=1 (yk ln yk - Ayk - Byk ln xk )2 → min . = k=1 Решение краевой задачи о кручении сплошных и полых цилиндрических образцов . . . Отсюда можно получить линейную систему нормальных уравнений Aˆ Aˆ N 2 ˆ k=1 yk + B N 2 k=1 yk ln xk N N 2 2 k=1 yk ln yk , k=1 yk ln xk = N N 2 2 2 ˆ +B k=1 yk ln xk k=1 yk ln xk = ln yk . С учетом введенных обозначений получаем расчетные формулы для оценок параметров модели (8): N N N 2 2 2 2 k=1 yk ln xk - k=1 yk ln xk · k=1 yk 2 N N N 2 2 2 2 k=1 yk · k=1 yk ln xk - k=1 yk ln xk N 2· N 2 ln x ln y - N 2 ln x · N 2 y y y k k k k=1 k k=1 k k=1 k k=1 yk ln yk 2 N N N 2 2 2 2 k=1 yk · k=1 yk ln xk - k=1 yk ln xk cˆ = exp m ˆ = N k=1 yk2 ln yk · ) ( ( ) . ln xk ln yk , (13) Как будет показано ниже, по сравнению с МНК-оценками (10), полученными для линеаризованной модели, использование оценок (13) позволяет существенно уменьшить погрешность отклонения расчетных и экспериментальных данных для исходной степенной аппроксимации (8): ∆= р 2 N э k=1 yk - yk N э 2 k=1 yk 1/2 · 100%, (14) где N - количество реализаций; ykэ и ykр - экспериментальные и расчетные данные соответственно. 4. Результаты расчетов. Методика (5) решения краевой задачи реализована численно «шагами» по времени с малым шагом дискретизации по временной и пространственной координатам. Приращения компонент тензора деформаций ползучести вычислялись методом Эйлера. При реализации предложенного метода все производные аппроксимировались соответствующими разностными соотношениями, а интегралы вычислялись по квадратурным формулам численного интегрирования. Далее приведем полученные результаты моделирования для стали 45 и сплава АМГ-6М. 4.1. Сталь 45. В работе [11] для стали 45 приведены результаты испытаний на кратковременную ползучесть при температурах 725 и 740 ℃ в условиях одноосного растяжения и при кручении сплошных (R1 = 0) стержней. Значения упругих констант материала при данном уровне температур следующие [11]: E = 1.7 · 105 МПа, µ = 0.3. На рис. 1 точками показаны экспериментальные значения осевой деформации εz при одноосном растяжении образцов из стали 45 при температуре 725 ℃ [11, рис. 5]. На основании этих данных для каждого напряжения были найдены экспериментальные значения скорости деформации ползучести p˙э , представленные в табл. 1. Согласно методике, изложенной в п. 3, по формулам (13) были определены значения параметров реологической модели (7): c = 7.046 · 10-9 (МПа)-m , m = 4.2, в то время как оценки, полученные по формулам (10), следующие: c = 1.263 · 10-9 (МПа)-m , m = 4.648. Для сравнения в табл. 1 также приведены результаты теоретических значений скорости установившейся ползучести p˙р с параметрами c и m, полученными по формулам (10) и (13) соответственно, а также величины оценки погрешности ∆ для 513 Ц в е т к о в В. В. скоростей деформаций в соответствии с формулой (14). Анализ этой информации позволяет сделать вывод, что использование оценок (13) дает более точные оценки параметров модели установившейся ползучести. Данные расчета деформации ползучести с использованием оценок (13) приведены на рис. 1 сплошными линиями, штриховыми линиями показаны результаты расчетов [11]. Для экспериментов на кручение в работе [11] приведены значения интенˆ Однако для прогносивности напряжений Sˆ в характеристической точке R. зирования характеристик напряженно-деформированного состояния по схеме (5) необходимо знать значения крутящего момента. Согласно [10], положение характеристической точки применительно к кручению стержней с до- 1 Axial strain, εz 0.15 3 2 0.1 4 0.05 5 0 0.5 1 1.5 Time, t (hours) 2 2.5 Рис. 1. Значения осевой деформации при растяжении образцов из стали 45 при температуре 725 ℃: точки - экспериментальные данные ([11, рис. 5]); сплошные линии - расчетные данные по модели (6), (13); штриховые линии - расчетные данные [11]; 1 - σ = 53.96 МПа; 2 - σ = 49.05 МПа; 3 - σ = 43.16 МПа; 4 - σ = 39.24 МПа; 5 - σ = 25.51 МПа [Figure 1. Values of axial strain under tension for specimens made of the Steel 45 at a temperature of 725 ℃: dots - experimental data (from [11, Fig. 5]); solid lines - calculated data by model (6), (13); dashed lines - calculated data by [11]; 1 - σ = 53.96 MPa; 2 - σ = 49.05 MPa; 3 - σ = 43.16 MPa; 4 - σ = 39.24 MPa; 5 - σ = 25.51 MPa] Таблица 1 Зависимость скорости установившейся ползучести p˙ от приложенного напряжения σ = const при растяжении стали 45 [Dependencies for steady-state creep rate p˙ on applied load σ = const under tension for the Steel 45] no. σ (MPa) p˙э (hr)-1 p˙р (hr)-1 [11] p˙р (hr)-1 (10) p˙р (hr)-1 (13) 1 2 3 4 5 53.96 49.05 43.16 39.24 25.51 1.35709 · 10-1 8.19814 · 10-2 5.33530 · 10-2 3.82695 · 10-2 4.04542 · 10-3 1.54378 · 10-1 9.32407 · 10-2 4.78672 · 10-2 2.92890 · 10-2 3.42310 · 10-3 1.42002 · 10-1 9.11787 · 10-2 5.03315 · 10-2 3.23177 · 10-2 4.36347 · 10-3 1.33035 · 10-1 8.91421 · 10-2 5.21031 · 10-2 3.49124 · 10-2 5.71554 · 10-3 14.1 7.6 5.0 5.45313 · 10-1 2.92874 · 10-1 1.56127 · 10-1 8.31005 · 10-2 1.28745 · 10-2 3.51548 · 10-1 2.35042 · 10-1 1.51227 · 10-1 9.28795 · 10-2 1.42293 · 10-2 5.01939 · 10-1 2.64972 · 10-1 1.31610 · 10-1 6.07196 · 10-2 3.09341 · 10-3 14.0 28.6 8.1 T = 725 ℃ ∆, % T = 740 ℃ 1 2 3 4 5 ∆, % 514 58.86 53.96 49.05 44.15 29.43 -1 5.13033 · 10 2.24865 · 10-1 1.46536 · 10-1 5.41082 · 10-2 1.80545 · 10-2 Решение краевой задачи о кручении сплошных и полых цилиндрических образцов . . . статочной точностью можно определить как координату пересечения эпюр, соответствующих упругому и установившемуся распределениям касательных напряжений: 4 4 4 ˆ = 3R2 1 - (R1 /R2 ) = 3 R2 - R1 R 3 3 3 4 1 - (R1 /R2 ) 4 R2 - R1 (15) ˆ = 3R2 /4). Найдя положение харак(для сплошных цилиндров R1 = 0 и R теристической точки, из упругого распределения касательных напряжений определяем крутящий момент: M= J τˆ π(R24 - R14 ) τˆ 2π(R23 - R13 ) = = τˆ, ˆ ˆ 2 3 R R (16) √ ˆ 3 - значение касательного напряжения в характеристической где τˆ = S/ точке. Радиусы стержней и значения крутящих моментов, вычисленные по формуле (16), приведены в табл. 2. Расчетная схема (5) позволяет определить величину относительного угла закручивания θ = θ(t). Для сравнения результатов расчета с экспериментальными данными и данными расчета [11] при кручении для найденного θ определялась интенсивность деформаций в характеристической точке: ˆ R εˆi (t) = √ θ(t) 3 (17) На рис. 2 приведены значения интенсивности деформаций в характеристической точке - экспериментальные (точки), расчетные данные, полученные в работе [11] (штриховые линии) и найденные по схеме (5) с использованием модели (7) и оценками параметров (13), и по формуле (17) (обозначены сплошными линиями). Здесь и далее под «экспериментальным» значением интенсивности деформаций при кручении будем понимать расчетноэкспериментальную величину, определяемую из соотношения (17). Также в [11] для одного эксперимента при температуре 725 ℃ приведена зависимость угла закручивания для образца радиуса R2 = 9.988 мм при Таблица 2 Значения параметров при экспериментальных исследованиях на кручение сплошных стержней из стали 45 [Parameter values of experimental study of torsion of the solid cylinders made of the Steel 45] no. R2 (mm) ˆ (mm) R 1 2 3 4 9.995 10 9.988 9.975 7.496 7.5 7.491 7.481 Sˆ (MPa) τˆ (MPa) M (N · mm) 28.32 24.92 22.66 14.73 59222.3 52193.9 47278.4 30611.1 28.32 25.49 59240.1 53316.1 T = 725 ℃ 49.05 43.16 39.24 25.51 T = 740 ℃ 1 2 9.996 9.996 7.497 7.497 49.05 44.15 515 Ц в е т к о в В. В. 0.1 Strain intensity, εˆi 0.08 Рис. 2. Значения интенсивности деформаций в характеристической точке при кручении образцов из стали 45 при температуре 725 ℃: точки - экспериментальные данные ([11, рис. 5]); сплошные линии - расчетные данные по методике (5); штриховые линии - расчетные данные [11]; 1 - Sˆ = 49.05 МПа; 2 - Sˆ = 43.16 МПа; 3 - Sˆ = 39.24 МПа; 4 - Sˆ = 25.51 МПа 2 1 0.06 0.04 3 4 0.02 0 0.5 1 1.5 Time, t (hours) 2 Torsion angle, θ (rad/mm) 0.01 Рис. 3. Изменение угла закручивания для стержня из стали 45 при температуре 725 ℃ (M = 50521.5 Н · мм): точки - экспериментальные данные ([11, рис. 4]); сплошная линия - расчетные данные по методике (5); штриховая линия - расчетные данные [11] 0.008 0.006 [Figure 3. Torsion angle for the specimen made of the Steel 45 at a temperature of 725 ℃ with M = 50521.5 N · mm: dots - experimental data (from [11, Fig. 4]); solid line - calculated data by method (5); dashed line - calculated data by [11]] 0.004 0.002 0 0.2 0.4 0.6 Time, t (hours) 0.8 0.16 1 2 3 Axial strain, εz 0.12 0.08 4 0.04 0 516 [Figure 2. Values of strain intensity at the charachteristic point under torsion for specimens made of the Steel 45 at a temperature of 725 ℃: dots - experimental data (from [11, Fig. 5]); solid lines - calculated data by method (5); dashed lines - calculated data by [11]; 1 - Sˆ = 49.05 MPa; 2 - Sˆ = 43.16 MPa; 3 - Sˆ = 39.24 MPa; 4 - Sˆ = 25.51 MPa] 5 0.3 0.6 0.9 Time, t (hours) 1.2 Рис. 4. Значения осевой деформации при растяжении образцов из стали 45 при температуре 740 ℃: точки - экспериментальные данные ([11, рис. 6]); сплошные линии - расчетные данные по модели (6), (13); штриховые линии - расчетные данные [11]; 1 - σ = 58.86 МПа; 2 - σ = 53.96 МПа; 3 - σ = 49.05 МПа; 4 - σ = 44.15 МПа; 5 - σ = 29.43 МПа [Figure 4. Values of axial strain under tension for specimens made of the Steel 45 at a temperature of 740 ℃: dots - experimental data (from [11, Fig. 6]); solid lines - calculated data by model (6), (13); dashed lines - calculated data by [11]; 1 - σ = 58.86 MPa; 2 - σ = 53.96 MPa; 3 - σ = 49.05 MPa; 4-σ = 44.15 MPa; 5 - σ = 29.43 MPa] Решение краевой задачи о кручении сплошных и полых цилиндрических образцов . . . Рис. 6. Значения осевой деформации при растяжении образцов из сплава АМГ-6М при температуре 450 ℃: точки - экспериментальные данные ([14, рис. 3]); сплошные линии - расчетные данные по модели (6), (13); штриховые линии - расчетные данные [14]; 1 - σ = 19.6 МПа; 2 - σ = 14.7 МПа; 3 - σ = 12.25 МПа; 4 - σ = 9.82 МПа; 5 - σ = 7.5 МПа [Figure 6. Values of axial strain under tension for specimens made of the AMG6M alloy at a temperature of 450 ℃: dots - experimental data (from [14, Fig. 3]); solid lines - calculated data by model (6), (13); dashed lines - calculated data by [14]; 1 - σ = 19.6 MPa; 2 - σ = 14.7 MPa; 3 - σ = 12.25 MPa; 4 - σ = 9.82 MPa; 5 - σ = 7.5 MPa] 0.12 1 Strain intensity, εˆi 0.09 0.06 2 0.03 0 0.25 0.5 0.75 Time, t (hours) 1 0.012 1 2 0.009 Axial strain, εz Рис. 5. Значения интенсивности деформаций в характеристической точке при кручении образцов из стали 45 при температуре 740 ℃: точки - экспериментальные данные ([11, рис. 6]); сплошные линии - расчетные данные по методике (5); штриховые линии - расчетные данные [11]; 1 - Sˆ = 49.05 МПа; 2 - Sˆ = 44.15 МПа [Figure 5. Values of strain intensity at the charachteristic point under torsion for specimens made of the Steel 45 at a temperature of 740 ℃: dots - experimental data (from [11, Fig. 6]); solid lines - calculated data by method (5); dashed lines - calculated data by [11]; 1 - Sˆ = 49.05 MPa; 2 - Sˆ = 44.15 MPa] 3 0.006 4 0.003 5 0 0.5 1 Time, t (hours) 1.5 Strain intensity, εˆi Рис. 7. Значения интенсивности деформа0.012 ций в характеристической точке при кручении образцов из сплава АМГ-6М при темпе1 ратуре 450 ℃: темные точки - экспериментальные данные для сплошных образцов 0.009 ([14, рис. 3]); светлые точки - экспериментальные данные для трубчатых образцов ([14, рис. 3]); сплошные линии - расчетные данные по методике (5); штриховые линии - 0.006 2 расчетные данные [14]; 1 - Sˆ = 19.6 МПа; 2 - Sˆ = 12.25 МПа; 3 - Sˆ = 9.82 МПа; 4 - 3 Sˆ = 7.5 МПа 0.003 [Figure 7. Values of strain intensity at 4 charachteristic point under torsion for specimens made of the AMG-6M alloy at a temperature of 450 ℃: dark dots - experimental data for solid specimens (from 0 0.3 0.6 0.9 1.2 1.5 [14, Fig. 3]); white dots - experimental Time, t (hours) data hollow specimens (from [14, Fig. 3]); solid lines - calculated data by method (5); dashed lines - calculated data by [14]; 1 - Sˆ = 19.6 MPa; 2 - Sˆ = 12.25 MPa; 3 - Sˆ = 9.82 MPa; 4 - Sˆ = 7.5 MPa] 517 Ц в е т к о в В. В. кручении моментом M = 50521.5 Н · мм. Результаты показаны на рис. 3. Здесь погрешность отклонения расчетных данных для угла закручивания θ от экспериментальных данных вычислялась по формуле ∆θ = р 2 N э k=1 θk - θk N э 2 k=1 θk 1/2 (18) · 100%, где N - количество экспериментальных значений; θkэ и θkр - экспериментальные и расчетные значения угла закручивания соответственно. Величина погрешности (18) для расчетных данных по методу настоящей работы составила 3.6%, а для расчетных данных [11] - 8.4%. Аналогичные исследования для стали 45 проводились при температуре 740 ℃. По одноосным экспериментальным кривым ползучести (см. рис. 4) с использованием методики, изложенной в п. 3, были найдены оценки параметров реологической модели (7) по формуле (13): c = 5.085 · 10-14 (МПа)-m , m = 7.342. Экспериментальные данные и результаты расчетов на растяжение и кручение сплошных образцов приведены на рис. 4 и 5 соответственно (обозначения те же). Значения параметров при экспериментальных исследованиях приведены в табл. 2. 4.2. Сплав АМГ-6М. В работе [14] приведены результаты испытаний на кратковременную ползучесть сплава АМГ-6М при температуре 450 ℃ в условиях одноосного растяжения и при кручении сплошных и толстостенных цилиндрических образцов. Упругие константы для данного материала следующие: E = 0.7 · 105 МПа, µ = 0.33. С использованием экспериментальных одноосных кривых ползучести [14, рис. 3] для каждого уровня напряжений определялась скорость деформации ползучести; эти данные приведены в табл. 3. Далее согласно методике, изложенной в п. 3, были найдены значения параметров реологической модели (7) с использованием оценок (13): c = 1.039 · 10-7 (МПа)-m , m = 4.262. Результаты расчетов при растяжении представлены на рис. 6. Здесь же приведены экспериментальные и расчетные данные из [14]. Величины погрешностей расчетных данных для скоростей деформаций ползучести по данным работы [14] и по оценкам (10) и (13) для метода наименьших квадратов по отношению к экспериментальным данным приведены в последней строке табл. 3. Таблица 3 Зависимость скорости установившейся ползучести p˙ от приложенного напряжения σ = const при растяжении сплава АМГ-6М при температуре 450 ℃ [Dependencies for steady-state creep rate p˙ on applied load σ = const under tension for the AMG-6M alloy at a temperature of 450 ℃] no. σ (MPa) p˙э (hr-1 ) p˙р (hr-1 ) [14] p˙р (hr-1 ) (10) p˙р (hr-1 ) (13) 1 2 3 4 5 19.6 14.7 12.25 9.82 7.5 3.34605 · 10-2 9.97811 · 10-3 3.97138 · 10-3 2.00677 · 10-3 7.17912 · 10-4 3.67606 · 10-2 8.53492 · 10-3 3.99199 · 10-3 1.98536 · 10-3 7.71509 · 10-4 3.10171 · 10-2 9.83745 · 10-3 4.75132 · 10-3 1.96568 · 10-3 6.70320 · 10-4 3.34675 · 10-2 9.81964 · 10-3 4.51440 · 10-3 1.75919 · 10-3 5.57711 · 10-4 10.2 7.3 1.8 ∆, % 518 Решение краевой задачи о кручении сплошных и полых цилиндрических образцов . . . Таблица 4 Значения параметров при экспериментальных исследованиях на кручение сплошных и трубчатых цилиндрических образцов из сплава АМГ-6М при температуре 450 ℃ [Parameter values of experimental study of torsion of solid and hollow cylinder specimens made of the AMG-6M alloy at a temperature of 450 ℃] no. R1 (mm) R2 (mm) ˆ (mm) R Sˆ (MPa) τˆ (MPa) M (N · mm) 1 2 3 4 0 8 0 8 10 11 10 11 7.5 9.657 7.5 9.657 19.6 12.25 9.82 7.5 11.32 7.07 5.67 4.33 23700.3 12131.6 11874.3 7427.5 Значения параметров при испытаниях на кручение сплошных и толстостенных цилиндрических образцов приведены в табл. 4. Результаты расчетов интенсивности деформаций в характеристической точке показаны на рис. 7: темными точками отмечены испытания на кручение сплошных образцов, светлыми точками - испытания трубчатых образцов. Выводы. Как следует из информации, представленной на рис. 2, 3, 5, 7, наблюдается соответствие данных расчета на кручение цилиндрических образцов по реологической модели при сложном напряженном состоянии (7) и методике решения краевой задачи (5) с экспериментальными данными, а также с результатами расчетов по методу характеристической точки [11, 14]. В отличие от отмеченных работ [11, 14], в которых осуществлялась аппроксимация всех экспериментальных значений интенсивности деформаций лишь в одной характеристической точке, предложенная методика позволяет по известным одноосным кривым ползучести осуществлять расчет напряженно-деформированного состояния для любых пространственно-временных сечений цилиндрического образца. В качестве примера на рис. 8 и 9 приведены расчетные распределения касательных напряжений в различные временные сечения для сплошного и полого цилиндров соответственно. Из информации, представленной на рис. 8 и 9, следует, что с течением времени происходит перераспределение касательных напряжений по радиусу 8 Shear stress, τ (MPa) Рис. 8. Расчетные эпюры касательных напряжений τ = τ (r, t) при кручении сплошного цилиндра из сплава АМГ-6М при температуре 450 ℃ при Sˆ = 9.82 МПа в различные моменты времени: 1 - t = 0; 2 - t = 0.04 ч; 3 - t = 0.1 ч; 4 - t = 1.5 ч [Figure 8. Calculated diagram for shear stress τ = τ (r, t) under torsion of a solid specimen made of the AMG-6M alloy at a temperature of 450 ℃ with Sˆ = 9.82 MPa at different time points: 1 - t = 0; 2 - t = 0.04 hr; 3 - t = 0.1 hr; 4 - t = 1.5 hr] 6 4 3 2 0 2 1 4 2 4 6 Radius, r (mm) 8 10 519 Ц в е т к о в В. В. Рис. 9. Расчетные эпюры касательных напряжений τ = τ (r, t) при кручении полого цилиндра из сплава АМГ-6М при температуре 450 ℃ при Sˆ = 12.25 МПа в различные моменты времени: 1 - t = 0; 2 - t = 0.01 ч; 3 - t = 0.02 ч; 4 - t = 1.5 ч [Figure 9. Calculated diagram for shear stress τ = τ (r, t) under torsion of an hollow specimen made of the AMG-6M alloy at a temperature of 450 ℃ with Sˆ = 12.25 MPa at different time points: 1 - t = 0; 2 - t = 0.01 hr; 3 - t = 0.02 hr; 4 - t = 1.5 hr] Shear stress, τ (MPa) 8 7 2 1 3 6 4 8 9 10 11 Radius, r (mm) вследствие деформации ползучести. Следует также отметить, что расчетные эпюры касательных напряжений действительно пересекаются в окрестности характеристической точки, теоретическое положение которой было получено расчетным путем по формуле (15). Таким образом, новая методика идентификации параметров даже для простейшей модели установившейся ползучести позволяет существенно снизить погрешность аппроксимации первичных кривых стационарной ползучести (см. рис. 1, 4, 6 и табл. 1, 3). Это, в свою очередь, позволило снизить погрешность решения краевой задачи для угла закручивания стержня (см. рис. 3). Однако в целом расчет напряженно-деформированного состояния посредством решения краевой задачи и методом характеристической точки дают близкие результаты и хорошо согласуются с экспериментальными данными. Конкурирующие интересы. Заявляю, что в отношении авторства и публикации этой статьи конфликта интересов не имею. Авторская ответственность. Я несу полную ответственность за предоставление окончательной версии рукописи в печать. Окончательная версия рукописи мною одобрена.

About the authors

Vitaly V Tsvetkov

Samara State Technical University

Email: vi.v.tsvetkoff@mail.ru
Postgraduate Student; Dept. of Applied Mathematics & Computer Science 244, Molodogvardeyskaya st., Samara, 443100, Russian Federation

References

Supplementary files