Обобщённая модель неупругого деформирования и разрушения толстостенной сферической оболочки при ползучести Текст научной статьи по специальности «Физика»

Механика деформируемого твердого тела

УДК 539.376

В.П. Радченко, М.Н. Саушкин, С.Н. Кубышкина

ОБОБЩЕННАЯ МОДЕЛЬ НЕУПРУГОГО ДЕФОРМИРОВАНИЯ И РАЗРУШЕНИЯ ТОЛСТОСТЕННОЙ СФЕРИЧЕСКОЙ ОБОЛОЧКИ ПРИ ПОЛЗУЧЕСТИ

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

Постановка задачи

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

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

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

Исходной информацией для построения обобщенной реологической модели служат кривые ползучести конструктивного элемента в координатах "обобщенное перемещение — обобщенная нагрузка", полученные при стационарных внешних воздействиях, а также обобщенная диаграмма упругопластического деформирования. При этом возможны два способа получения этой информации:

1) проведение натурного эксперимента над конструктивным элементом;

2) решение соответствующей краевой задачи при стационарных внешних воздействиях. Описание неупругого реологического деформирования ряда элементов конструкций в пределах первых двух стадий ползучести при помощи обобщенной модели на основе указанных двух способов получения исходной информации было рассмотрено в [1—5]. В отмеченных работах показана эффективность такого подхода для решения некоторого класса задач. Однако вопросы

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

Естественное обобщение результатов работ [1-5] на случай третьей стадии и разрушения с учетом пластической деформации было предпринято в [6] для описания реологического неупругого деформирования и разрушения растягиваемого резьбового соединения и толстостенной трубы под действием внутреннего давления [7].

Целью настоящей работы является развитие идей [6, 7] на случай толстостенной сферы под действием внутреннего давления, при этом получение необходимой информации для построения соответствующей обобщенной модели сферы основано на решении соответствующей краевой задачи при постоянных значениях внутреннего давления.

Реологическая модель и критерий разрушения материалов при сложном напряженном состоянии

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

В качестве основной реологической модели в настоящей работе используется модель, предложенная для одноосного напряженного состояния в [8-9] и обобщенная на случай сложного напряженного состояния в [10].

Основной вариант определяющих соотношений имеет вид

е- = е, + 9, + Р, ;

/ = 0 Г УУ /9 =

Г УУ

Рц = и„ + ;

при

о, ,,,----------а

2

£ ап

* 1 = с| о-

а

-о,,---а

2 1 2

, (?) = Х ик ($);

а

(1)

(2)

(з)

'я[а( - оПр )-1 • В - / ] [.] • В > 0; (4)

з 1 0

0, [...] • В £ 0, при 2 К --о 2 > апр;

(5)

(6)

(7)

^({)=Х(*);

к

V, (I)=(1+т / а) - тт (/1*1 ()+/ (о+/«));

// = ГУУ

О \П2 1 -5 1 а

к\ * а

а

,0, [.]к > °;

[...]к > °;

О;

= О° (1 +

(Ъ = г(е2 °,}д,} + «(5° КЛ.

(8)

(9)

(1°)

Здесь е,, е-, д,, р- — полная, упругая, пластическая деформации и деформация ползучести

соответственно; и.., V-,,

V V

— вязкоупругая, вязкопластическая и вязкая составляющие де-

1

*

формации ползучести; а у, а0 — соответственно компоненты истинного и номинального тензоров напряжений; E, m — упругие константы материала; E2, S, S0 — соответственно интенсивности тензоров пластической деформации, истинных и номинальных напряжений;

1, а, п — константы модели, описывающие диаграмму мгновенного упругопластического

деформирования; апр — предел пропорциональности; 1к, ak, Ьк, с, n2, тх, а*— константы

модели, при помощи которых описывается первая и вторая стадия ползучести материала и ее обратимая после разгрузки часть; mk, т”к — коэффициенты Пуассона для обратимой и необратимой компонент деформации ползучести; pq, pk — соответственно активные пластические и вязкопластические деформации, которые можно было бы наблюдать при отсутствии пуассо-новского сужения материала; r(E2) и a(S 0) задаются степенными аппроксимациями вида

y(E2) = g • Em2; a(S0) = a, • Sm3, (11)

где g, m2, a,, m3 — константы модели, контролирующие процессы разупрочнения материала

при пластической деформации и деформации ползучести соответственно; о — скалярный па-

раметр поврежденности. В формулах (4)-(7) использованы следующие обозначения:

B = ^3а n -2а° -апр^sign^| а w -2а0j; а0 = ап + а22 + азз. (12)

Расчет пластической qу и вязкопластической vу деформаций осуществляется в главных

осях, поэтому суммирования по индексу v в формулах (3), (8), (12) не выполняется. Очевидно, что при записи (1)—(12) использовалась гипотеза соосности тензоров напряжений и деформаций. Модель (1)-(12) описывает процесс неупругого деформирования с изотропным разупрочнением. Методика определения всех параметров модели (1)—(11) представлена в [8-9]. Она осуществляется на основании одноосных кривых ползучести, полученных при а0 = const и доведенных до разрушения, и диаграммы упругопластического деформирования.

Для прогнозирования времени разрушения материала t* используется критерий разрушения энергетического типа вида [8-10]

n(t.)=fSA+tfSA=1, (13)

W J A J A (S0)

где A и A (S0) соответственно критические величины работ разрушения истинного напряжения на пластической деформации и на деформации ползучести. При этом материал находится в не разрушенном состоянии, если W(t) < 1, и разрушается при выполнении (13). В общем случае

величина A (S0) имеет степенную аппроксимацию вида

Ac (S0 ) = aA • S0mA,

где aA, mA — константы материала, которые могут быть определены по результатам одноосных испытаний по методике, изложенной в [8-9].

Таким образом, построение модели (1)—(12) для описания неупругой деформации при сложном напряженном состоянии и критерия разрушения (13) не требует дополнительных экспериментальных затрат, так как все параметры модели определяются по результатам одноосных испытаний. В частности, в работах [10-11] приведены константы используемой модели для некоторых материалов.

Следует отметить, что согласно (3)-(4) пластические деформации описываются такими же по структуре уравнениями, как и вязкопластическая компонента v деформации ползучести (8),

т.е. также развивается по времени. При этом величина 1 много больше max{l} и скорость

к

деформации пластичности на порядок выше скорости деформации ползучести. При таком предположении за то время, когда пластическая деформация при заданном тензоре напряжений достигнет, согласно (3)-(4), асимптотического значения, накопленная деформация ползучести за это же время будет пренебрежимо малой по сравнению с пластической деформацией, т.е. p У (t) много меньше q у (t). Такой подход к описанию пластических деформаций соответствует

так называемым эндохронным теориям пластичности [12-13]. Обоснованность применения такого подхода приведена в [10]. Соотношения (3)-(4) задают вариант теории пластичности без

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

Решение краевой задачи о неупругом реологическом деформировании и

разрушении толстостенной сферической оболочки

Изложим теперь методику решения краевой задачи для толстостенной сферической оболочки аналогично [11], [14], где решена задача о неупругом деформировании толстостенного цилиндра.

Предположим, что толстостенная сферическая оболочка с внутренним радиусом Я1 и внешним Я2 под действием внутреннего давления р(Г) находится в условиях неупругого деформирования. Задача рассматривается в сферических координатах г, в, у, В силу симметрии задачи деформации сдвига и касательные напряжения равны нулю, и она решается в главных осях, причем ев = еу, ив = иу.

Полные деформации ei ( = г, в, у) описываются соотношениями (1)—(12).

Уравнение совместности деформаций в силу симметрии задачи записывается в виде

г + ев = ег. (14)

дг

Из уравнений (1) и (14) имеем соотношение

де

г^ + ев - ег = я (г, 0, (15)

дг

где функция я(г,г) задается следующим образом:

0Рв + °<7в

я(г,0 = рг + дг - рв - д- - г\-^- + | • (16)

Из закона Г ука (2), записанного в главных осях

= И0 - 2 ив

Е

(і - т)л° - тв

Е

Єу =е-

выражается величина

Г 5е- + - Л )=.

ог Е

Дифференцируя второе уравнение (17) и подставляя его в (18), получаем

г (1 - - гтИ + (1 + /4ив - и0) = Е • я (г, г),

дг дг

из которого с учетом уравнения равновесия

находим дифференциальное уравнение относительно сг° :

дг г дг 1 -т г

Решение уравнения (3.7) с граничными условиями

при фиксированном значении времени ґ записывается в виде

(17)

с учетом которой уравнение (15) принимает вид

\ив - и0 )= Я(г, г). (18)

г + 2л° = 2и- (19)

ог

0Л° + 4 Ост0 = 2Е Я(гґ) (2П)

Л° в =- р(0, Л° в = ° (21)

1г =Н.1 1г =л,

ев =

Л°(г • ґ) =

2Е (х,ґ) , 2Е 1 г 3 , , ,

—?----\ ---------ах —г- х я(х,ґ)ах- р(ґ) -

з(1 - т) х з(1 - т)г3 { ^ ’

[ К3 \ 2Е % (х, ґ) , 2Е 1 г 2 , ч ,

1л 3 -?—\ Г——- ах —т—г Гх я(х, ґ) ах- р(ґ)

[) 3(1 - ті х 3(1 - М)г3 [ ^ ^

(г3 -Пз\пз

\кг-

С учетом (22) из (19) находится выражение для ав :

авв(гґ) = зГ3Е ) Г0х + з(Е ) Л Гх2Я(X ґ) °х- р(ґ) 3(1 - т) х з(1 - т)г [

(г( + І [ \ 2Е І3я(х, ґ) е 1 і3 з ( ґ) а (ґ) ^

^ зт-т) іг^зт-т гз Iх я(х ґ) *- р(,)

(зз)

Соотношения (22) и (2з) задают распределение напряжений а°, а- по радиусу в зависимости от времени. Нетрудно видеть, что при я(г, ґ)=° они совпадают с упругим решением для толстостенной сферы.

Численная реализация расчета кинетики напряженно-деформируемого состояния во времени и разрушения толстостенной сферической оболочки под действием внутреннего давления р в условиях ползучести осуществлялась по хорошо известному методу "шагами по времени", аналогично [11], [14]. Временной интервал разбивался на малые отрезки времени [ґг-, ґі+1] с шагом А ґі, внутри которого внутреннее и внешнее напряженные состояния считались постоянными и соответствующими моменту ґ = ґі • Приращения компонент деформаций в конце отрезка при ґ = ґі+1 находились по формулам (1)-(11) по методу Эйлера. При этом интегралы в формулах (22) и (2з) вычислялись по соответствующим квадратурным формулам численного интегрирования при каждом значении времени ґ = ґі , а для производных по радиусу использовались их конечно-разностные аппроксимации. Время до разрушения в соответствии с (1з) определялось следующим образом: расчет продолжался до того момента времени ґ = ґ», при котором в какой-либо точке сферы выполнялось бы условие ^(ґ») = 1.

В случае упругопластического деформирования расчет происходит также "шагами по времени" (как при ползучести), но здесь время можно рассматривать как некоторое "внутреннее время", играющее роль параметра нагружения. При этом сначала рассчитывается упругое решение при ґ=°, а затем строится упругопластическое решение на основании (1)-(4), (9)—(1з) по аналогии с деформацией ползучести.

Следует отметить, что во всех формулах (14)-(2з) используются номинальные напряжения, тогда как реологические и упругие деформации рассчитываются в соответствии (1)—(1з) с использованием истинных напряжений.

Численное решение краевой задачи

В качестве примера выполнено модельное численное исследование неупругого деформирования и разрушения толстостенной сферической оболочки из сплава ЭИ 698 с внутренним радиусом [1 = °.2° м и внешним — [ 3 = °.28 м при действии внутреннего давления р и температуре Т=7°° ° С. Параметры модели взяты из работы [1°] и представлены в табл. 1 и 2.

Используя разработанную выше методику, была решена серия краевых задач о ползучести и разрушении сферической оболочки при постоянных давлениях, а также задача о упругопластическом разрушении сферической оболочки, при этом для упругопластической задачи в определяющих соотношениях (1)—(1 °) было положено ру = °. В качестве примера на рис.1

сплошными линиями приведены расчетные зависимости для перемещения от времени на внутреннем радиусе сферы, вызванные деформациями ползучести и пластическими деформациями для ряда постоянных значений давления. На рис.2 приведена расчетная диаграмма "упругопластического деформирования" сферы в координатах "перемещение внутреннего радиуса — давление". Информация, представленная на рис.1 и 2, являлась исходной для построения обобщенной модели сферической оболочки.

/

Значения параметров модели (1)—(13) для описания

деформации пластичности сплава ЭИ 698 при Т=700 ° С_____________________

Ипр . МПа Е -105, МПа а, МПа щ П1 Х1, МПа 1 А,, МДж/м 3 Ш2

500.3 1.52 1.18-10-7 1.995 3.77 -10-3 284.3 0

Таблица 2

Значения параметров модели (1)—(13) для описания деформации ползучести сплава ЭИ 698 при Т=700 ° С

и*, МПа £ 4, ч-1 ак х10-4 Ък х10-4 С х10-5 П2 шх а, МПа-1-Ш ш3 а, МПа1~ША ША

490.5 1 0.2 2.96 4.44 2.51 2.9 10.96 9.56-103 -2.03 12.2 0

Р и с. 1. Кривые стационарной ползучести толсто- Р и с. 2. Диаграмма упругопластического де-стенной сферы: сплошные линии — численный экспе- формирования толстостенной сферы в коорди римент; штриховые линии — расчет по модели (24)-(25). натах "радиальное перемещение е — внутрен Цифры: 1 — р0=210 МПа; 2 — р0=220 МПа; 3 — р0=250 МПа; 4 — р0=270 МПа; 5 — р0=290 МПа.

нее давление р 0": сплошная линия — численный

расчет по модели

эксперимент; штриховая линия (24)-(25).

Р и с. 3. Кривые ползучести толстостенной сферы при нестационарных режимах нагружения: сплошные линии — численный эксперимент; штриховые линии — расчет по модели (5.1)—(5.2). Цифры: 1 — р0=190 МПа; 2 — р0=250 МПа; 3 — р0=210 МПа; 4 — р0=230 МПа; 5 — р0=200 МПа.

1 1; 1 1 ... |!

/ и V 1 4 г/

1

31 1,*

Р и с. 4. Кривые ползучести толстостенной сферы при нестационарных режимах нагружения: сплошные линии — численный эксперимент; штриховые линии — расчет по модели (5.1)—(5.2). Цифры: 1

— р0=190 МПа; 2 — р0=280 МПа; 3 — р0=220 МПа; 4

— р0=0 МПа; 5 — р0=210 МПа.

Построение обобщенной реологической модели для сферической оболочки

Если в качестве внешней "обобщенной силы" взять внутреннее давление, а в качестве "обобщенного перемещения" — перемещение на внутреннем радиусе, то, как это следует из рис. 1 и 2, наблюдается полная аналогия в эволюции сферы и одноосного образца при растяжении в процессе ползучести и упругопластического деформирования. В связи с этим обобщенная модель реологического деформирования и разрушения толстостенной сферы как целого по аналогии с одноосным напряженным состоянием [8-9] принимается в виде

Ч .

e = e + e + p,

e=

G

er =

\ q(t) < q^,

[4a(q(t) - q^)n1 -ep(t)] [...] > °

|0, [■■■] < 0, q(t) > qпр ,

p = u + v + w ; i(t) = £ uk (t), Uk (t) = l

(24)

(г. > n2

ak - uk (t)

v q* 0

v(t) = Z vk(t),

к

vk(t) =

Ik

n2

q

q*

- vk (t) [■■■] < 0,

[•■■]> 0;

W(t) = c

q*

Ч = ч0 (1 + ^

= удё р + ачр.

Здесь е — полное обобщенное перемещение; ё, ёр — упругая и пластическая компоненты е; р — компонента обобщенного перемещения, вызванная ползучестью; и, V, ^ — вязкоупругая, вязкопластическая и вязкая составляющие р; О — интегральная (либо локальная) упругая податливость в заданной точке; ч0 и ч соответственно номинальное и истинное значения обобщенной нагрузки, чпр — предел пропорциональности на диаграмме упругопластического

деформирования в координатах e — q; a, n1, l, 1k

ak, bk,

q* —константы модели.

Детальный анализ данных показал, что в общем случае g = y(ep) и a=a(q0) и для них можно использовать степенные аппроксимации вида

g = g (eP)”2, a = aiq”3- (25)

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

В качестве критерия разрушения (локального или интегрального) конструктивного элемента используется соотношение вида

■ q(t) de p (t) f q(t) dp(t)

A*p

A* (q0)

= 1.

(26)

Здесь A*p = const, A*c (q0) = aAqmA ■

В рассматриваемом случае роль обобщенной силы q играет давление р, роль обобщенного перемещения e — радиальное перемещение на внутренней поверхности сферы.

Как отмечалось выше в качестве исходной информации для построения уравнений (24) и критерия (25) служат кривая упругопластического деформирования сферы в координатах e —qo и серия кривых ползучести при q0 = const вплоть до разрушения, которые в настоящей работе получены численным решением соответствующей краевой задачи и представлены на

k

n, c, ті

рис.1 и 2 сплошными линиями. Используя эту информацию по методике, изложенной в [8], были определены все параметры модели (24)-(26), значения которых представлены в таблицах 3 и 4.

Таблица 3

Значения параметров модели (24)-(26) для описания деформации пластичности

толстостенной сф ерической оболочки из сплава ЭИ 698 п ри T=700° C

q^,МПа E ■ 1Q 5, МПа а, МПа щ n1 g , МПа 1 A*, МДж/м 3 m2

212.1 6.13 6.86 -1Q-9 2.178 3 -1Q-3 15.39 -Q.6Q9

Таблица 4

Значения параметров модели (24)-(26) для описания деформации ползучести толстостенной сферической оболочки из сплава ЭИ 698 при Т=700 ° С

q*, М Па к 4 Т1 ак х 1Q-5 bk х 10-5 С х10-8 П2 m1 a1, МПа~l~m m3 a1 , МПа-1-mA mA

1 7.75 10-4 4.29 8.58

2QQ 2 4.99 10-3 2.35 4.7 3.5 3.343 6.188 75.487 -Q.923 Q.9 Q

3 1.68 10-1 2.Q7 4.15

Для сравнения на рис.2 штриховыми линиями показаны рассчитанные по предложенной модели (24)-(26) диаграммы упругопластического деформирования в координатах р0 — е (р0— истинное давление). На рис.1 штриховыми линиями представлены результаты расчета кривых ползучести по модели (24)-(26) для стационарных режимов нагружения при р0 = const. Адекватность модели (24)-(26) с данными решения краевой задачи была выполнена для нестационарных режимов нагружения. На рис.3-5 сплошные линии соответствуют данным численного решения краевой задачи о ползучести толстостенной сферической оболочки (численный эксперимент), а штриховые — данным расчета по модели (24)-(26). Как следует из представленных рисунков, наблюдается соответствие данных численного эксперимента и расчета по модели (24)-(26).

Приведенный пример свидетельствует о том, что в ряде случаев связь обобщенных перемещений с обобщенными нагрузками может быть сформулирована соотношениями типа (24)-(26). Достоинством соотношений типа (24)-(26) является то, что существенно понижается размерность исходной краевой задачи и, как следствие этого, на порядок уменьшается количество вычислений. Краевую задачу при этом необходимо решать лишь для нескольких стационарных режимов нагружения при p0(t) = const. При переменных режимах нагружения нет необходимости в решении краевой задачи, а достаточно воспользоваться " одномерными" соотношениями (24)-(25).Это немаловажный фактор с точки зрения оценки вычислительной устойчивости расчетных алгоритмов, так как при решении полных краевых задач с учетом времени приходится многократно пересчитывать одни и те же итерационные процедуры (реальные элементы конструкций, например, в энергетической промышленности в условиях ползучести эксплуатируются до 104 ^105 часов).

Р и с. 5. Кривые ползучести толстостенной сферы при нестационарных режимах нагружения: сплошные линии — численный эксперимент; штриховые линии — расчет по модели (24)-(25). Цифры: 1 — /|=190 МПа, 2 — р=270 МПа, 3 — р=210 МПа, 4 — р=230 МПа, 5 — р=220 МПа, 6 — р=280 МПа.

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

ее ползучести.

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

1. Самарин Ю.П., Еремин Ю.А. Метод исследования ползучести конструкций // Проблемы прочности. 1985. №4. С.40-45.

2. Самарин Ю.П., Клебанов Я.М. Обобщенные модели в теории ползучести конструкций. Самара: Поволж. отд-ние инженер. акад. РФ. Самарский гос. техн. ун-т, 1994. 197с.

3. Еремин Ю.А., Кайдалова Л.В., Радченко В.П. Исследование ползучести балок на основе аналогии структуры уравнений состояния материала и элементов конструкций // Машиноведение. 1983. №2. С.67-74.

4. Самарин Ю.П., Радченко В.П. О решении краевых задач механики сплошных сред методами теории управления // Механика и прикладная математика. Тула: Приокское кн. изд-во, 1988. С.3-8.

5. Радченко В.П. Компактное представление напряженно-деформированного состояния в точке при решении краевых задач реологии / Прочностные и динамические характеристики машин и конструкций. Пермь: Перм. ГТУ, 1994. С.46-55.

6. Радченко В.П., Кубышкина С.Н. Об одном подходе к описанию реологического деформирования и разрушения резьбовых соединений // Изв. вузов. Машиностроение. 1998. №4-6. С.31-35.

7. Радченко В.П., Кубышкина С.Н. Об одном подходе к описанию реологического деформирования и разрушения толстостенной трубы // Математическое моделирование м краевые задачи: Тр. седьмой межвуз. конф. Ч.1. Самара: СамГТУ, 1997. С .111-116.

8. Радченко В.П. Энергетический вариант одноосной теории ползучести и длительной прочности // Журн. прикл. мех. и техн. физики. 1991. №4. С.172-179.

9. Радченко В.П., Кичаев Е.К. Феноменологическая реологическая модель и критерий разрушения металлов при одноосном напряженном состоянии // Проблемы прочности. 1991. №11. С. 13-19.

10. Радченко В.П. Математическая модель неупругого деформирования и разрушения металлов при ползучести энергетического типа // Вестн. СамГТУ. Сер. Физ.-мат. науки. Вып.4. 1996. С.43-63.

11. Радченко В.П., Кубышкина С.Н. Математическая модель реологического деформирования и разрушения толстостенной трубы // Вестн. СамГТУ. Сер. Физ.-мат. науки. Вып.6.1998.С.23-34.

12. Кадашевич Ю.И., Мосолов А.Б. Эндохронные теории пластичности: основные положения, перспективы развития // Изв. АН СССР. МТТ. 1989. №1. С.161-168.

13. МосоловА.Б. Эндохронная теория пластичности. М.: Институт проблем механики АН СССР, 1988.С.44.

14. Миронова С.Н. Решение задачи упругопластического деформирования и разрушения толстостенной трубы на основании эндохронной теории ползучести //Вестн. СамГТУ. Сер. Физ.-мат. науки. Вып.4. 1996.С.85-92.