Научная статья на тему 'Численный метод определения параметров модели ползучести разупрочняющегося материала'

Численный метод определения параметров модели ползучести разупрочняющегося материала Текст научной статьи по специальности «Математика»

CC BY
184
41
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
ПОЛЗУЧЕСТЬ РАЗУПРОЧНЯЮЩЕГОСЯ МАТЕРИАЛА / РАЗНОСТНЫЕ УРАВНЕНИЯ / ОБОБЩЕННАЯ РЕГРЕССИОННАЯ МОДЕЛЬ / НЕЛИНЕЙНАЯ РЕГРЕССИЯ / СРЕДНЕКВАДРАТИЧЕСКОЕ ОЦЕНИВАНИЕ / STRAIN SOFTENING CREEP MODEL / DIFFERENCE EQUATIONS / GENERALIZED REGRESSION MODEL / NONLINEAR REGRESSION / ROOT MEAN SQUARE EVALUATION

Аннотация научной статьи по математике, автор научной работы — Зотеев Владимир Евгеньевич, Макаров Роман Юрьевич

Тенденции к уменьшению массы машин при улучшении их качества, а также стремление к наиболее полному использованию механических свойств материалов требуют постоянного совершенствования и развития известных методов расчета и анализа напряженно-деформированного состояния материалов в условиях ползучести. В статье рассматривается новый численный метод оценки параметров математической модели ползучести разупрочняющегося материала на основе экспериментальных диаграмм, построенных по результатам испытаний при различных напряжениях. В основе метода лежит обобщенная регрессионная модель, построенная на основе разностных уравнений, описывающих диаграммы ползучести. Полученные соотношения между коэффициентами разностного уравнения и параметрами деформации ползучести позволяют свести задачу параметрической идентификации к итерационной процедуре среднеквадратичного оценивания коэффициентов, линейной на каждом шаге итерации обобщенной регрессионной модели. Проведена апробация разработанного численного метода на пяти экспериментальных кривых ползучести алюминиевого сплава, подтверждающая достоверность полученных соотношений и эффективность численного метода.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

A numerical method for the determination of parameters of the strain softening creep model

The trends of decreasing of weight of machines and increasing of their quality, and intention to the fullest use of mechanic properties of materials demand the development of numerical methods for analysis of the stress-strain state of materials undo the terms of creep. The article discusses the development of new numerical method for determining the parameters of the strain softening creep model. The base of new method is generalized regression model, which was built on the basis of difference equations for describing the creep. The relations between coefficients of difference equation and parameters of the strain softening creep model allow reduce the problem of parametric identification to an iterative procedure for RMS of coefficients of regression model, which is linear. The approbation of numerical method with five creep curves of aluminum alloy is accomplished. The approbation confirms scientific credibility of built relations and efficiency of new numerical method.

Текст научной работы на тему «Численный метод определения параметров модели ползучести разупрочняющегося материала»

Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки. 2016. T. 20, № 2. С. 328-341

ISSN: 2310-7081 (online), 1991-8615 (print) doi: http://dx.doi.org/10.14498/vsgtu1488

Математическое моделирование, численные методы и комплексы программ

УДК 519.651.3:539.376

ЧИСЛЕННЫЙ МЕТОД ОПРЕДЕЛЕНИЯ ПАРАМЕТРОВ МОДЕЛИ ПОЛЗУЧЕСТИ РАЗУПРОЧНЯЮЩЕГОСЯ МАТЕРИАЛА

В. Е. Зотеев, Р. Ю. Макаров

Самарский государственный технический университет, Россия, 443100, Самара, ул. Молодогвардейская, 244.

Аннотация

Тенденции к уменьшению массы машин при улучшении их качества, а также стремление к наиболее полному использованию механических свойств материалов требуют постоянного совершенствования и развития известных методов расчета и анализа напряженно-деформированного состояния материалов в условиях ползучести. В статье рассматривается новый численный метод оценки параметров математической модели ползучести разупрочняющегося материала на основе экспериментальных диаграмм, построенных по результатам испытаний при различных напряжениях. В основе метода лежит обобщенная регрессионная модель, построенная на основе разностных уравнений, описывающих диаграммы ползучести. Полученные соотношения между коэффициентами разностного уравнения и параметрами деформации ползучести позволяют свести задачу параметрической идентификации к итерационной процедуре среднеквадратичного оценивания коэффициентов, линейной на каждом шаге итерации обобщенной регрессионной модели. Проведена апробация разработанного численного метода на пяти экспериментальных кривых ползучести алюминиевого сплава, подтверждающая достоверность полученных соотношений и эффективность численного метода.

Ключевые слова: ползучесть разупрочняющегося материала, разностные уравнения, обобщенная регрессионная модель, нелинейная регрессия, среднеквадратическое оценивание.

© 2016 Самарский государственный технический университет. Образец для цитирования

Зотеев В. Е., Макаров Р. Ю. Численный метод определения параметров модели ползучести разупрочняющегося материала // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2016. Т. 20, № 2. С. 328-341. аок 10.14498/vsgtu1488. Сведения об авторах

Владимир Евгеньевич Зотеев (д.т.н., доц.; zoteev-ve@mail.ru; автор, ведущий переписку), профессор, каф. прикладной математики и информатики.

Роман Юрьевич Макаров (makaroman1@yandex.ru), аспирант, каф. прикладной математики и информатики.

Введение. Исследование поведения материалов в условиях ползучести является важной научно-технической задачей. Тенденции к уменьшению массы машин при улучшении их качества, а также стремление к наиболее полному использованию механических свойств материалов требуют совершенных методов расчета, в которых наиболее полно отражены механические свойства материалов [1]. Однако нелинейность определяющих уравнений ползучести и длительной прочности существенно затрудняет применение аналитических и повышает роль численных методов решения.

При высоких температурах характерными особенностями деформационно-прочностного поведения материалов являются отсутствие начальной стадии ползучести (стадия упрочнения) и сравнительно небольшая продолжительность третьей стадии (стадия разупрочнения) [2]. Вместе с тем именно третья стадия ползучести как стадия, предшествующая разрушению, является определяющей при расчете срока службы проектируемого элемента конструкции по катастрофическому критерию отказа [3]. Большинство существующих методик определения параметров моделей ползучести разупроч-няющихся материалов обладают одним существенным недостатком. Они не обеспечивают достаточную адекватность построенной модели данным эксперимента, вследствие чего обладают низкой помехозащищенностью и точностью оценок параметров модели [3-5]. Для устранения этого недостатка в данной работе предлагается новый численный метод нелинейного оценивания, в основе которого лежат разностные уравнения, описывающие диаграммы ползучести материала при различных постоянных номинальных напряжениях.

Определяющие уравнения. При математическом моделировании процессов деформирования и разрушения материалов используем энергетический подход, предложенный в [5]. Этот подход базируется на принципе суперпозиции упругой, пластической деформаций и деформации ползучести, а также методе разделения деформации ползучести. Основной вариант определяющих соотношений представлен в [5]. Из этих соотношений с учетом отсутствия вязкоупругой и вязкопластической составляющих в деформации ползучести p вытекают уравнения

p = cam, a = а0(1 + ш), ш = aap),

где ао и a — соответственно номинальное и истинное напряжения (ао ^ 0); ш — параметр поврежденности; c и m — константы модели, при помощи которых описывается вторая стадия ползучести; a — параметр модели, контролирующий процесс разупрочнения материала на третьей стадии ползучести. Из этой системы дифференциальных уравнений для нулевых начальных условий можно построить временные зависимости деформации ползучести

p(t,aoj) =--1—ln(l - amca0"+1t), (1)

a0j ma j

которые могут быть использованы для описания третьей стадии ползучести при различных постоянных значениях номинального напряжения aoj. Очевидно, что для третьей стадии ползучести важнейшим участком моделирования является промежуток времени непосредственно перед моментом разрушения материала, когда кривая ползучести практически совпадает с вертикальной асимптотой. Поэтому при выборе критерия «близости» модели

к эксперименту необходимо в первую очередь ориентироваться на величину отклонения прогнозируемого времени разрушения от данных эксперимента. Это означает, что при оценке параметров математической модели (1) следует минимизировать норму разности по временной координате:

INI2 = II* - ¿II2 ^ min,

где t — Nj) -мерный вектор данных эксперимента tk,j, к = 0,1, 2,...,

Nj — 1, вычисленный по M диаграммам ползучести, в каждой из которых взято по Nj точек, j = 1,2,..., M; i — N?) -мерный вектор результатов

вычислений tk,j, к = 0,1, 2,..., Nj — 1, по модели (1) при значениях деформации ползучести pk,j в соответствующих точкам эксперимента (tk,j ,Pk,j), к = = 0,1, 2,... ,Nj — 1, j = 1, 2,..., M. При этом зависимость tk,j от pk,j в явном виде можно получить из выражения (1):

4 ,j=—i1—exp(—maaoj pk,j)],

cm^m a (2)

к = 0,1,2,...,Nj — 1, j = 1,2,...,M.

Таким образом, задача параметрической идентификации модели для описания третьей стадии ползучести сводится к минимизации нелинейной функции от трех переменных c, m и a:

М N-1 f 1 1 2

J = £ {tk,j--m+i" t1 — exP(—ma^ojPk,j m ^ min.

j=1 k=0 cm%- a

Данная задача относится к классу задач нелинейной регрессии, методы исследования которой описаны в [6] и сопряжены со сложностью поиска решения нормальной системы нелинейных уравнений. В данной работе предлагается новый численный метод решения задачи нелинейной регрессии, в основе которого лежат разностные уравнения, описывающие результаты эксперимента. Такой подход позволяет свести решение нелинейной задачи к итерационной процедуре, на каждом шаге которой оценки коэффициентов разностного уравнения вычисляются на основе линейной обобщенной регрессионной модели [7-10]. Численный метод оценки параметров рассматриваемой модели ползучести на основе разностных уравнений по совокупности M экспериментально построенных при aoj = const кривых включает следующие шаги:

- предварительная обработка выборки результатов наблюдений для каждой j-той кривой с целью равномерной дискретизации данных эксперимента по переменной pk,j = Apjк, где Apj — шаг (период) дискретизации;

- построение разностных уравнений, связывающих несколько последовательных значений дискретной модели (2), коэффициенты которых известным образом связаны с параметрами c, m и a модели (1);

- формирование обобщенной регрессионной модели, описывающей результаты эксперимента и линейной относительно коэффициентов разностного уравнения;

- итерационная процедура среднеквадратичного оценивания коэффициентов обобщенной регрессионной модели;

- вычисление оценок параметров деформации ползучести на основе среднеквадратичных оценок коэффициентов обобщенной регрессионной модели;

- статистическая оценка погрешности результатов вычисления и построение доверительных интервалов для каждой из M кривых ползучести.

Построение разностных уравнений и обобщенной регрессионной модели, описывающих результаты эксперимента. На первом шаге численного метода проводится предварительная обработка результатов эксперимента для каждой из M кривых ползучести, построенных при постоянных значениях номинального напряжения ooj = const. Предварительная обработка заключается в сглаживании результатов наблюдений с помощью скользящих средних и формировании равномерной (по переменной pk,j с шагом Apj) выборки экспериментальных данных заданного объема Nj. Как показали численно-аналитические исследования процедуры предварительной обработки, систематическая погрешность от сглаживания на основе скользящих средних в точках (Pk,j ,tk,j) эксперимента на два порядка меньше, чем величина ek,j естественного разброса результатов наблюдений.

В основе построения разностных уравнений, связывающих несколько последовательных значений дискретной модели

ik,j =-- exp(-maCToj Apjk)] ,

сша'т^а (3)

k = 0,1, 2,...,Nj — 1, j = 1, 2,..., M,

лежит методика, описанная в [7]. Из формулы (3) при k = 2, 3,... ,Nj, j = = 1, 2,... ,M имеем следующее:

cma m+1atk,j = 1 — exp(- maa0j Apj k),

cmam+1atk-i,j = 1 — exp[—maa0j Apj (k — 1)] =

= 1 — exp(—maa0jApjk) exp(maa0jApj),

cmam+1 atk-2,j = 1 — exp [—maa0jApj(k — 2)] =

= 1 — exp (—maa0jApjk) exp(2maa0jApj).

Отсюда, переходя к разностям, можно получить соотношение

ln tk-1,j —tk-2,j = a0j Apj\i, k = 2,3,... ,Nj, j = 1,2,...,M, tk,j tk-1,j

где X1 = ma.

При k = 0 и k = 1 из формулы (3), соответственно, получаем

to,j = 0, j = 1, 2,..., M

¿1] =-т+г- [1 - ехР(—таоо] Ар])] , ] = l, 2,...,М.

ето0]+ а

Из последнего соотношения имеем

ехр^оо] Ар^ )отаат+1$1,у = ехр( Л^о/ Ар]) — 1.

Отсюда после логарифмирования получаем

1п [ехр(Л1ст0^- Ар2) — 1] — 1п ¿^ = Л1 оо] Ар] + Л2 + Л3 1п оо], з = 1,2,... ,М,

где Л2 = 1п(ста) и Л3 = т + 1.

Таким образом, имеем систему разностных уравнений, связывающих последовательность Ьк-2], tк-1,у и дискретных значений модели (3), ] = = 1, 2,..., М:

Ь о,2 = 0;

1п[ехр(Л1^о2 Ар]) — 1] — 1п ¿1,2 = Л1ао2 Ар] + Л2 + Л3 1п оо];

1п — 1к-2,] = оо] Ар] Л1, к = 2,3,...,^- — 1. Ь к] £ к-1]

(4)

Коэффициенты в системе разностных уравнений (4) связаны с параметрами деформации ползучести соотношениями

Л1 = та, Л2 = 1п(ста), Л3 = т + 1,

которые позволяют по оценкам коэффициентов Л] вычислять оценки параметров деформации ползучести. Используем соотношения

Ь-к] = ¿к,] + £к,], к = 0,1, 2,..., N — 1, 3 = 1, 2,...,М,

где Ьк] — результаты эксперимента, т.е. значения переменной Ь при аргументе Ар]к для каждой з'-той кривой; £к] —естественный разброс результатов наблюдения относительно модели (3). После линеаризации по переменным 61], £к-2], бк-1] и £к], в первом приближении получаем

1п ¿1 ] = 1п (Ь1 ] — е1 ]) « 1п Ь1 ] — 612,

, Ьк-1] — £ к-2] _ , Ьк-1] — Ьк-2] — (бк-1] — 6к-2]) 1п —---- — 1п ----- ~

Ь к] — £к-1] Ьк] — Ьк-1] — (бк] — 6к-1] ) (5)

Ьк-1] — Ьк-2] , 1 _ ^ '

~ 1п--\--£к—2] —

Ьк] Ьк-1] Ьк-1] Ьк-2]

Ьк] — Ьк-2] 6 . + 1 £ .

(Ьк-1] — Ьк-2]) (Ьк] — Ьк-1]) , Ьк] — Ьк-1] ,

Подставляя формулы (5) в систему уравнений (4), имеем математическую модель, которая в форме разностных уравнений описывает результаты эксперимента (рк] ,Ьк]), к = 0,1, 2,... — 1, 3 = 1, 2,..., М.

Полагая ео,, =0 с учетом нулевых начальных условий

ро,, = 0, ¿о,, = 0,

получаем

( ¿о,, =0;

1п 1п

ехр(Л1ао,Ар,) - 1

¿к— 1,, — ¿к—2,,

¿к,,' ¿к—1,,

+

= Л1^о, Ар, + Л2 + Лз 1п сто, - —^;

= сто, Ар, Л1 -¿к, ¿к—2,,

1

¿1,

1,, — ¿к—2,,') (¿к,,' — ¿к—1,, )

¿к—1,, — ¿к—2,, -ек—1 а

ек—2,7 +

(6)

¿к,, ¿к— 1,

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

где к = 2, 3,..., Ж,- - 1, j = 1, 2,..., М.

Система разностных уравнений (6) лежит в основе формирования обобщенной регрессионной модели. Введем следующие векторы и матрицы:

- Л = [ Л1 Л2 Лз ]Т — трехмерный вектор коэффициентов обобщенной регрессионной модели;

- 6л = [ Ь1 62 ••• 6м ]Т —блочный вектор размера Ж?' х , состоящий из М векторов 6,- размера (Ж, х 1) вида

Ь, =

1п

¿о,,

ехр(Л1СТр,- Ар,-) - 1 1п - ¿о,,

1п

j = 1, 2,.

М;

(7)

1п ¿М, —2,,' - ¿М,—3,,' ¿М — 1,, - ¿М—2,,

^ = [ ^2 ■ ■ ■ ^м ]Т — блочная матрица размера Ж, х 3^),

состоящая из М матриц размера (Ж, х 3) вида

^ =

еТ еТ е1 е2

0 0 0 сто, Ар, 1 1п сто, сто, Ар, 0 0

_ сто, Ар, 0 т

j = 1, 2,

М;

Т

м

— блочный вектор естественного разброса дан-

- ' = [

ных эксперимента относительно модели (остатков) размера х1

состоящий из М векторов е, размера (Ж, х 1):

= [ £о,,' е1,,' е2,,' еМ—1,, ]Т , j = 1 2,

М;

1

0

е

,

- П = [ ПТ ПТ ''' пМ ]Т — блочный вектор «невязки» обобщенной ре-

грессионной модели размера ^^М1 Ж, х 1), состоящий из М векторов П, размера (Ж, х 1):

П, = [ П1,, П2,, Пз,, ■

Пм,, ]Т , j = 1, 2,..., М,

где

Пг,, =

е1,,

П1,, = ео,, = 0, П2,, = ",

¿г— 1,, ¿г—3,,

¿г—2,, ¿г—3,,

"ег-3,, +

(¿г—2,, - ¿г—3,,) (^г—1,, - ¿г—2,,) 1 ^

ег- 2,, -

ег—1,,, г = 3, 4,..., Ж,;

Р

Р1 В В. . . В

В Р2 В. . . В

В В Р3 . . . В

В В В. . . Рм

¿г—1,, ¿г—2,,

блочно-диагональная квадратная мат-

рица линейного преобразования вектора остатков размера Ж, х

х £м=1 Ж,^, состоящая из М матриц Р, размера (Ж, х Ж,), ] = 1, 2,...,

М, и нулевых матриц В. Квадратные матрицы Р, —нижние треугольные трехдиагональные матрицы, элементы которых описываются следующим образом:

- диагональные элементы:

,

р=

г,г

Г

¿г— 1,, ¿г—2,, поддиагональные элементы:

г = 1;

г = 2;

г = 3, 4,..

(9)

,

,

рг,г 1 =

0,

¿г—1,, ¿г—3,,

(^г—2,, - ¿г—3,,) (^г—1,, - ¿г—2,,) '

подподдиагональные элементы: ,1

г = 2;

г = 3,4,

(10)

,

рг,г 2 =

¿г—2,, ¿г—3,,

г = 3,4,

(11)

Остальные элементы матриц Р, равны нулю: рк г = 0 при г > к и г<к - 2, к, г = 1, 2,..., Ж,. ,

1

1

1

Таким образом, с учетом введенных выше векторов и матриц получаем обобщенную регрессионную модель в виде

Ьх = FA + п; П = Ре.

(12)

Итерационная процедура среднеквадратического оценивания коэффициентов обобщенной регрессионной модели. В соответствие с выбранным критерием «близости» математической модели (3) экспериментальным данным ('Pk,j ,tk,j) минимизируется остаточная сумма квадратов

Q = INI2 = ||t - t||2 ^ min,

которую с учетом системы (12) при невырожденной матрице Р можно представить в виде

Q = ||е||2 = ||Р-1ЬХ - Р-1F\H2

или

Я = (Ъх — ГЛ) 1 П-1(ЪХ — ГЛ) = Ъ{ П-1ЪХ — 2Ъ{ О-1ГЛ + Л* Г 1 П-1РЛ,

где О = РРт — матрица размера (^ММ=1 N х ^ММ=1 ^^ . Дифференцируя остаточную сумму квадратов по переменным Л1 и приравнивая нулю результаты дифференцирования, получаем систему нормальных уравнений, которая в векторной форме имеет вид

(F - Wa)tQ-1FA = (F - Wx)1 П-1Ьх

1

(13)

где блочная матрица Шх = [ Ш1 Ш2 ■ ■ ■ Шм ]Т размера {^ММ=1N х 3^, состоящая из М матриц размера (Ж] х 3) вида

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Wj =

aoj Apj

0 0

0 0

1 - exp(-A10oj Apj)

0 0 0

00

j = 1, 2,..., M.

(14)

Так как элементы

мм = a°j A'j

j, 1 - exp(-A10ojApj)

матрицы Wx и элементы

bu = ln

exp(A1aojApj) - 1 11,j

0

0

вектора где 1? = 2 + N, ] = 1, 2,..., М, ¿1 = 2, зависят от величины А1, система уравнений (13), вообще говоря, нелинейная. Однако она легко разрешается относительно вектора коэффициентов А:

А = [(Р - - ^л)тО-1Ьл. (15)

Формула (15) лежит в основе итерационной процедуры среднеквадратичного оценивания элементов вектора А. Рекуррентное соотношение, описывающее итерационный процесс на основе (15), имеет вид

д(^+1) = [(р - ^)ТО-1Р]-1 (Р - г = 0,1,2,3,..., (16)

А Л Л

Т1л(») т.(»)

где матрица и вектор оЛ на г-том шаге вычисляются на основе г-того приближения Л«.

В соответствии с моделью (6) оценку начального приближения Л10) для

первого элемента вектора А = [ А1 А 2 А3 ]Т можно найти из условия минимизации

м N-1

?=1 1=2 " Отсюда получаем следующую оценку:

[1П - Ар, А 10)

I- ¿к*?' ¿к — 1 ?

2

—> ш1п.

Ем^Е -1 ,0? Ар? 1п ^-4к"2'

Л(0) = ^ к=2 0? ^ ¿м(17)

Л1 ЕМ=1(^? - 2),2?Ар2 . (7)

Алгоритм итерационной процедуры среднеквадратичного оценивания коэффициентов обобщенной регрессионной модели может быть описан следующим образом.

1. По формулам (8) и (9)—(11) вычисляются элементы матрицы регрессо-ров Р и матрицы линейного преобразования вектора остатков Р соответственно.

2. По формуле (17) вычисляется начальное приближение Л10) первой компоненты вектора коэффициентов А системы разностных уравнений (6).

3. По формулам (7) и (14) при значении А1 = Л1г) при г = 0 вычисляются элементы вектора Ьл и матрицы производных ^л.

4. Вычисляются правая и левая части системы нормальных уравнений (13).

5. По формуле (16) находится приближение л(г+1) вектора коэффициентов обобщенной регрессионной модели (12).

6. Сравниваются по модулю величины Л1г+1) и Л1г). Если |Л1г+1) - Л1г)| < 5, где 5 > 0 — заданное значение предельной абсолютной погрешности, то процесс вычислений заканчивается и за вектор оценок коэффициентов обобщенной регрессионной модели принимается вектор Л(г+1). Если |Л(г+1) - Л(г)| ^ 5, то, увеличив г на единицу, следует вернуться к шагу 3 данного алгоритма вычислений.

По найденным оценкам коэффициентов системы разностных уравнений (6) оценки параметров модели (1) можно вычислить по формулам

„ ехрЛ2 „ с „ Л1

с = —с-, т = Лз — 1, а = с-. (18)

Л1 Лз — 1

Апробация численного метода оценки параметров модели ползучести разу-прочняющегося материала на основе разностных уравнений. Описанный выше алгоритм численного метода оценки параметров модели (1) был апробирован по данным работы [11] при обработке пяти экспериментальных кривых ползучести алюминиевого сплава Д16Т (пруток диаметром 40 мм) при температуре 250 С при испытаниях на растяжение и сжатие (при постоянных номинальных напряжениях -о = 98.1, 88.3, 78.5, 73.6 и 68.6 МПа.

На этапе предварительной обработки экспериментальных данных были сформированы равномерные по координате р выборки равного объема N3 = 20, ] = 1,2,3,4, 5 со следующими периодами дискретизации: Др1 = 0.00522, Др2 = 0.00572, Дрз = 0.00647, Др4 = 0.00699 и Др5 = 0.00633. При этом установлено, что предварительная обработка результатов эксперимента на основе сглаживания с помощью скользящего среднего практически не меняет исходные диаграммы ползучести.

В процессе реализации описанного выше алгоритма численного метода оценки параметров модели (1) получены следующие промежуточные результаты.

Начальное приближение Л^ первой компоненты вектора коэффициентов системы разностных уравнений, вычисленное по формуле (17), следующее:

5 ^19 _ д„ 1 Ък-1Я—1<к-2,1

Т,1=^к=2 аоз ДРз 1п

Л(о) = ^ ^к=2 о] з у ^ = 0.6471.

1 18 Е5=1 -о] ДР?

В соответствии с алгоритмом была реализована итерационная процедура среднеквадратичного оценивания коэффициентов обобщенной регрессионной модели, в результате после 5 итераций были получены коэффициенты обобщенной регрессионной модели:

1 = 0.861, 2 = —29.708, з = 5.4449.

Отсюда, используя соотношения (18), находим оценки параметров модели (1):

с =ехр(Л^ = 1.46 ■ 10-13, т = Лз — 1 = 4.449, а = = 0.192.

1 з — 1

Таким образом, построенная по экспериментальным кривым ползучести математическая модель ползучести разупрочняющегося материала на примере алюминиевого сплава Д16Т при температуре 250 С имеет вид

= 1.2543 ■ 10-13а^49 — ехР(—0,861-0зРк,з)) . (19)

Исходя из этого в качестве зависимости деформации ползучести от времени будем использовать функцию вида

p(t, aoj) = -

1

0.861ctq.

1п(1 - 2.254 ■ 10-1VQ/49t).

Остаточные суммы квадратов для каждой ^'-той кривой ползучести в относительных единицах составили: 24.2, 8.8, 7.8, 11.3 и 4.9% соответственно. По совокупности всех кривых остаточная сумма квадратов в относительных единицах равна 8.3%. Коэффициент детерминации Л2, который косвенно может служить оценкой адекватности построенной модели экспериментальным данным по совокупности всех кривых ползучести, составляет Я2 = 0.96, что свидетельствует о высокой степени адекватности построенной модели (19) данным эксперимента [11].

На основе статистического анализа результатов математического моделирования найдены оценки погрешности для параметров модели и построены доверительные интервалы для кривых ползучести. На рисунке при различных номинальных напряжениях ио] изображены:

- экспериментальные кривые (точки);

- кривые (сплошные линии), построенные на основе разработанной модели (19);

- границы доверительных интервалов (штриховые линии), вычисленные с использованием неравенства Чебышева при доверительной вероятности в = 0.97.

0.15

^ 0.1 а

•а

СО

а

а) а;

б 0.05

50 100 150 Time, t (hours)

100 200 300 Time, t (hours)

200 400 600 Time, t (hours)

Расчетные доверительные интервалы на основе построенной модели (19) при различных значениях номинального напряжения: ст01 = 98.1 МПа (а), ст02 = 83.3 МПа (b), <7оз = 78.5 МПа (c) [Calculated confidence intervals on the basis of the constructed model (19) at different values of the nominal stress: a01 = 98.1 MPa (a), ст02 = 83.3 MPa (b),

стоз = 78.5 MPa (c)]

Выводы. Таким образом, разработан новый численный метод оценивания параметров модели ползучести, описывающей разупрочнение на третьей стадии, по совокупности кривых ползучести, построенных при обработке результатов эксперимента для различных значений номинального напряжения.

В основе метода лежат разностные уравнения, описывающие результаты эксперимента, сформированная на их базе обобщенная регрессионная модель, коэффициенты которой известным образом связаны с параметрами модели ползучести, и итерационная процедура среднеквадратического оценивания коэффициентов обобщенной регрессионной модели. Такой подход к решению задачи нелинейного оценивания позволяет минимизировать величину отклонения модели от данных эксперимента и тем самым обеспечить высокую точность оценок параметров модели. Результаты апробации нового численного метода при обработке результатов эксперимента в форме кривых ползучести алюминиевого сплава Д16Т при температуре 250 С подтвердили достоверность полученных соотношений и выводов, а также высокую эффективность численного метода в задачах оценивания параметров рассматриваемой модели ползучести. ORCIDs

Владимир Евгеньевич Зотеев: http://orcid.org/0000-0001-7114-4894 Роман Юрьевич Макаров: http://orcid.org/0000-0001-8620-1222

БИБЛИОГРАФИЧЕСКИЙ СПИСОК

1. Малинин Н. Н. Прикладная теория пластичности и ползучести. М.: Машиностроение, 1975. 387 с.

2. Соснин О. В., Любашевская И. В., Новоселя И. В. Сравнительные оценки высокотемпературной ползучести и разрушения конструкционных материалов // ПМТФ, 2008. Т. 49, №2. С. 123-130.

3. Bellenger E., Bussy P. Phenomenological modeling and numerical simulation of different modes of creep damage evolution // International Journal of Solids and Structures, 2001. vol.38, no. 4. pp. 577-604. doi: 10.1016/S0020-7683(00)00042-1.

4. Расчеты и испытания на прочность. Расчетные методы определения несущей способности и долговечности элементов машин и конструкций. М.: ВНИИМАШ, 1982. 90 с.

5. Радченко В. П., Еремин Ю. А. Реологическое деформирование и разрушения материалов и элементов конструкций. М.: Машиностроение-1, 2004. 264 с.

6. Демиденко Е. З. Линейная и нелинейная регрессия. М.: Финансы и статистика, 1981. 302 с.

7. Зотеев В. Е. Параметрическая идентификация диссипативных механических систем на основе 'разностных уравнений/ ред. В. П. Радченко. М.: Машиностроение, 2009. 344 с.

8. Зотеев В. Е., Заусаева М. А. Определение параметров двумерных динамических процессов на основе разностных схем // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2010. №1(20). С. 154-161. doi: 10.14498/vsgtu781.

9. Зотеев В. Е. Математические основы построения разностных уравнений для задач параметрической идентификации // Вестн. Сам. гос. техн. ун-та. Сер. Физ.-мат. науки, 2008. №2(17). С. 192-202. doi: 10.14498/vsgtu623.

10. Зотеев В. Е. Параметрическая идентификация линейной динамической системы на основе стохастических разностных уравнений // Матем. моделирование, 2008. Т. 20, № 9. С. 120-128.

11. Соснин О. В., Горев Б. В., Никитенко А. Ф. Энергетический вариант теории ползучести. Новосибирск: ИГиЛ СО АН СССР, 1986. 95 с.

Поступила в редакцию 21/III/2016; в окончательном варианте — 13/V/2016; принята в печать — 27/V/2016.

3oTeeB B. E., MaKapoB P. K>.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Vestn. Samar. Gos. Techn. Un-ta. Ser. Fiz.-mat. nauki

[J. Samara State Tech. Univ., Ser. Phys. & Math. Sci.], 2016, vol. 20, no. 2, pp. 328-341

ISSN: 2310-7081 (online), 1991-8615 (print) doi: http://dx.doi.org/10.14498/vsgtu1488

MSC: 74S30

A NUMERICAL METHOD FOR THE DETERMINATION OF PARAMETERS OF THE STRAIN SOFTENING CREEP MODEL

V. E. Zoteev, R. Yu. Makarov

Samara State Technical University,

244, Molodogvardeyskaya st., Samara, 443100, Russian Federation.

Abstract

The trends of decreasing of weight of machines and increasing of their quality, and intention to the fullest use of mechanic properties of materials demand the development of numerical methods for analysis of the stress-strain state of materials undo the terms of creep. The article discusses the development of new numerical method for determining the parameters of the strain softening creep model. The base of new method is generalized regression model, which was built on the basis of difference equations for describing the creep. The relations between coefficients of difference equation and parameters of the strain softening creep model allow reduce the problem of parametric identification to an iterative procedure for RMS of coefficients of regression model, which is linear. The approbation of numerical method with five creep curves of aluminum alloy is accomplished. The approbation confirms scientific credibility of built relations and efficiency of new numerical method.

Keywords: strain softening creep model, difference equations, generalized regression model, nonlinear regression, root mean square evaluation.

ORCIDs

Vladimir E. Zoteev: http://orcid.org/0000-0001-7114-4894 Roman Yu. Makarov: http://orcid.org/0000-0001-8620-1222

REFERENCES

1. Malinin N. N. Prikladnaia teoriia plastichnosti i polzuchesti [Applied theory of plasticity and creep]. Moscow, Mashinostroenie, 1975, 387 pp. (In Russian)

2. Sosnin O. V., Lyubashevskaya I. V., Novoselya I. V. Comparative estimation of high-temperature creep and rupture of structural materials, J. Appl. Mech. Tech. Phys., 2008, vol.49, no. 2, pp. 262-266. doi: 10.1007/s10808-008-0036-0.

© 2016 Samara State Technical University. Please cite this article in press as:

Zoteev V. E., Makarov R. Yu. A numerical method for the determination of parameters of the strain softening creep model, Vestn. Samar. Gos. Tekhn. Univ., Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. & Math. Sci.], 2016, vol. 20, no. 2, pp. 328-341. doi: 10.14498/vsgtu1488. (In Russian) Authors Details:

Vladimir E. Zoteev (Dr. Techn. Sci.; zoteev-ve@mail.ru; Corresponding Author), Professor, Dept. of Applied Mathematics & Computer Science.

Roman Yu. Makarov (makaroman1@yandex.ru), Postgraduate Student, Dept. of Applied Mathematics & Computer Science.

3. Bellenger E., Bussy P. Phenomenological modeling and numerical simulation of different modes of creep damage evolution, International Journal of Solids and Structures, 2001, vol.38, no. 4, pp. 577-604. doi: 10.1016/S0020-7683(00)00042-1.

4. Raschety i ispytaniia na prochnost'. Raschetnye metody opredeleniia nesushchei sposobnosti i dolgovechnosti elementov mashin i konstruktsii [Design calculations and strength testing. Calculation methods for the determination of bearing capacity and durability of machine elements and constructions]. Moscow, All-Russian scientific research institute for standardization and certification in the branch of mechanical engineering, 1982, 90 pp. (In Russian)

5. Radchenko V. P., Eremin Yu. A. Reologicheskoe deformirovanie i razrusheniia materialov i elementov konstruktsii [Rheological deformation and fracture of materials and structural elements]. Moscow, Mashinostroenie-1, 2004, 264 pp. (In Russian)

6. Demidenko E. Z. Lineinaia i nelineinaia regressiia [Linear and Nonlinear Regression]. Moscow, Finansy i statistika, 1981, 302 pp. (In Russian)

7. Zoteev V. E. Parametricheskaia identifikatsiia dissipativnykh mekhanicheskikh sistem na osnove raznostnykh uravnenii [Parametric Identification of Dissipative Mechanical Systems Based on Difference Equation], ed. V. P. Radchenko. Moscow, Mashinostroenie, 2009, 344 pp. (In Russian)

8. Zausaeva M. A., Zoteev V. E. Definition of Parameters of 2D Dynamical Processes on the Basis of Difference Schemes, Vestn. Samar. Gos. Tekhn. Univ. Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. & Math. Sci.], 2010, no. 1(20), pp. 154-161 (In Russian). doi: 10.14498/vsgtu781.

9. Zoteev V. E. Mathematical base for difference equations formulation in parametrical identification problems, Vestn. Samar. Gos. Tekhn. Univ. Ser. Fiz.-Mat. Nauki [J. Samara State Tech. Univ., Ser. Phys. & Math. Sci.], 2008, no. 2(17), pp. 192-202 (In Russian). doi: 10.14498/vsgtu623.

10. Zoteev V. E. Parametrical identification of linear dynamical system on the basis of stochastic difference equations, Matem. Mod., 2008, vol.20, no. 9, pp. 120-128 (In Russian).

11. Sosnin O. V., Gorev B. V., Nikitenko A. F. Energeticheskii variant teorii polzuchesti [Energetic variant of the creep theory]. Novosibirsk, Institut Gidrodinamiki, 1986, 95 pp. (In Russian)

Received 21/III/2016;

received in revised form 13/V/2016;

accepted 27/V/2016.

i Надоели баннеры? Вы всегда можете отключить рекламу.