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

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

Агахи К.А. fKamilla@imec.msu.ru ), Кузнецов В.Н., Ковальков В.К., Шестериков С.А.

Работа выполнена при поддержке гранта Российского фонда фундаментальных исследований № 02-01-00075.

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

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

В работе [1] были предложены определяющие соотношения деформационного типа для ортотропной упругопластической среды, чувствительной к виду напряжённого состояния, обладающей свойством дилатансии, т. е. появления объёмной деформации в результате сдвиговой деформации, и нелинейным изменением объёма.

Для ортотропной среды введём условие пластичности, восходящее к условию Мизеса, которое записывается в виде £ > 1, где смешанный инвариант £ определён следующим образом:

£ = Я + М <»•

здесь - тензор деформаций, Ауы, Я°к1 - материальные тензоры, в - константа

материала. Первое и третье слагаемые представляют собой смешанные инварианты тензора деформаций и материальных тензоров, второе - скаляр, где в = ац8ц, (в -

объёмная деформация), Я = 3Яук15у5к1.

Простейшей структурой тензоров Л и Л°к1 в главных осях ортотропии, позволяющей описать экспериментальные данные, является следующая. Отличными от нуля принимаем компоненты: Л1111,..., Л3333 ,Л212,..., Л3131; Л10111,..., Л°333 , остальные компоненты равны нулю.

Рассмотрим случай, когда объёмная деформация среды нелинейно зависит от напряжений [2] и введём соответствующий "предел объёмной упругости" в"; тогда получаем ещё одно условие появления нелинейных деформаций: в >в*. Введём в области £ > 0 обобщённую "интенсивность деформаций" :

е = £ (2).

Тогда "условие пластичности" примет вид системы неравенств:

/ е > 1

{в>г (3)-

Ясно, что неравенства могут рассматриваться как независимые. Определяющие соотношения примем в виде

= [ст - ата

здесь оц - тензор напряжений, сь.к1 - тензор упругих модулей, а^к1, а°°к1 - материальные тензоры, К - модуль объёмного сжатия а (е) и а0 (в) - материальные функции, учитывающие нелинейность.

Отметим некоторые особенности соотношений (4). Слагаемое соответствует

упругому деформированию; второе слагаемое а( е )аь-к1£ь- содержит функцию нелинейности а (е), аналогично теории малых упругопластических деформаций; слагаемое Као (в) учитывает нелинейное изменение объёма. Член а°ыба (е) связан с

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

Структура соотношений (4) позволяет построить их потенциал, который является естественным обобщением потенциала теории малых упругопластических деформаций [3] и имеет вид:

1 о е в

1 3 с с с с

^ = 2 ^ы^ы— 4 л) еа (е Уе-л]ао (вУв.

4 Л1 у ' Л

1 в

Существование потенциала обеспечивается условиями:

а1111 = Л(Л1Ш — Лв), а1111 = 2ЛЛ1111, а1212 = -сР, а3232 = ЛЛ4 ,

здесь введён инвариант с = 5 Сук1 — 15у5к1 ^, имеющий смысл удвоенного среднего

модуля сдвига (в изотропной среде с = 2^).

Для определения материальных функций а(е), ао(в) можно использовать опыт

на кручение цилиндрического образца, по методике, изложенной ниже.

Предварительно рассмотрим задачу об упругопластическом кручении однородного, трансверсально - изотропного толстостенного кругового цилиндра

(задача А). На торцах цилиндра действует система сил, приводящихся к крутящему моменту М и продольной силе Р, а на внутренней поверхности действует давление О, массовые силы и начальные напряжения отсутствуют. Внешняя боковая поверхность свободна от напряжений. Выберем цилиндрическую систему координат (^в^), в которой ось Ъ совпадает с осью цилиндра, являющейся его осью анизотропии, и плоскость z=0 - его плоскость симметрии.

Цилиндр предполагается «длинным» и, следовательно, деформации и напряжения зависят только от одной координаты г, а перемещения имеют вид:

u = ^г), V = тгс+^(г), w = (5),

где ^г), ^(г), п(г) - неизвестные функции радиуса, а Т и Я - некоторые константы.

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

= С18] -а

а4 = с4484 - а (е)а4484 - Као (в) (6),

=

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

1 - гг, 2 - ее, 3 - zz, 4 - ze, 5 - ге, 6 - гё;

Граничные условия на боковой поверхности имеют вид:

°гг = = °вв=0 при 1 = °"гг = °гг =°вв= 0 при 1 = Г2 (7),

где г1 и г2 - внутренний и внешний радиусы цилиндра соответственно. На торцах выполняются интегральные граничные условия

г2 г2 2п^О22Ыт = Р, 2ёг = М (8).

1 1

Выражая с помощью формул Коши деформации через перемещения, запишем уравнения равновесия в перемещениях, подставляя вместо напряжений их выражения через деформации из определяющих соотношений (6), с учетом формул (5) получим три уравнения для трёх неизвестных функций ^г), ^(г), п(г). Перенося в них нелинейные члены вправо, имеем:

и 41) +1«'(1)-^ = F

?(1) + 1*(1)-У- = С (9),

П(1 = н

где Б, О, Н имеют следующий вид: : 1 - (г);

I 1 и(г) | 1

^ = о(е) и"(г) + -и'(г)—^ +ф'(в)—(а„и'(г) + а12и(г) + а13Я + с11а10) + Ка'0

О = 0(е)а11а ^(г) + ±£'(г+ о'(е)^±

С11 С12

Н = о

-44

а11 а12

С11 С12

44

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

Рассмотрим теперь постановку и решение задачи В об определении функций о,юо (в частности, о восстановлении точной зависимости между касательным

напряжением т и деформацией сдвига у) из опыта на кручение толстостенного или сплошного цилиндрического образца. Эта нестандартная задача ставится следующим образом: в задаче об упругопластическом кручении с растяжением и давлением цилиндра, изложенной выше, в число неизвестных включаются функции о(е), а>0 (0) ; в

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

Таким образом имеем следующую постановку задачи:

1) уравнения равновесия в перемещениях (9);

2) граничные условия на боковой поверхности (7) и на торцах (8);

3) деформация сдвига на внешней боковой поверхности.

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

= С18] -о

о4 = с4484 - о (е)а4484

^6 = (С11 - С12)86 -0(е)(а11 -а12)86 Метод состоит в следующем:

1. Обозначим касательное напряжение в нормальном сечении а2в = а4 через

т

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

о 2М

п

т

2. Находим из соотношений (10) функцию ш0 нулевого приближения по формуле:

а>0 = ш4 = -

С4£4 Т0

(11),

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

*г=2£4 (12)-

3. Решается задача А - стандартная краевая задача о кручении, рассмотренная выше, с известной функцией пластичности ш =ш° и обычными граничными условиями на свободной поверхности, (т.е. при этом исключается условие, задающее на поверхности деформацию сдвига).

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

С44(г24 - Г4) а4

поверхности, которая имеет вид: М0 = тгу

___М4

4Т Г2

Г2 г Зш0(е) ёг А

5. Сравнение экспериментальной зависимости "крутящий момент М - сдвиг у",

т.е. Мэ=Мэ(у) и вычисленной аналогичной кривой М°=М°(у) - показывает погрешность, внесённую предположением о равномерности распределения напряжений по сечению. Разность функций Мэ и М0, умноженная на

нормирующий коэффициент к =---1-, даёт поправку для

П

вычисления зависимости т1=т1(у) в первом приближении, т.е. строится зависимость «касательное напряжение - сдвиг» первого приближения по формуле

т1=т0+к(Мэ-М0).

6. Указанным в пункте 2 способом, с помощью формул (10), (11), по кривой т =т (у) находится функция ш первого приближения. После этого строится первое приближение решения задачи А и т.д.

Таким образом, для (п+1) - го приближения получим тп+1=тп+к (Мэ-Мп).

Сравниваем полученную зависимость тп1(у) с зависимостью предыдущего приближения тп, и если они совпадают с заданной точностью, то итерационный процесс завершается. Функция пластичности шп , для которой вычисленная и экспериментальная зависимости «крутящий момент - сдвиг», в известном смысле, совпадают, может быть принята в качестве искомой истинной материальной функции. Другими словами, зависимость т-у в (п+1)-м приближении представляет собой истинную зависимость, которая, в принципе, могла бы быть получена непосредственно из опыта на кручение весьма тонкостенного образца.

Заметим, что в случае, когда экспериментальная кривая Мэ=Мэ(у) является монотонной и выпуклой, как показывает опыт практических вычислений, удаётся получить на 4-ой - 5-ой итерации совпадение расчетной кривой и экспериментальной с точностью, которую обеспечивает эксперимент.

а4454

Предложенный метод определения истинной зависимости т-у легко обобщается для основной модели (7) с двумя функциями нелинейности о и о0. При этом несколько изменяются шаги 2 и 3 описанной процедуры.

В 2. ищутся функции о и о0, материальные тензоры а,у, а^ нулевого приближения, а затем первого и п-го по методике, описанной в пп. 5 и 6. В 3. решается такая же задача А, но уже для основной модели (7). В заключение приведём решение задачи А методом упругих решений, когда в каждом приближении получается точное аналитическое решение.

Систему нелинейных уравнений (19) будем решать методом последовательных приближений. Для этого введем итерационную схему:

<+1 (г )+ -и'к+1 (г )-^ = Р

г г

¿+1 (г ) +-¿+1 (г)-£^ = Ок (13),

г г

П (г) + = Нк

г

где Рк = Р (ик, ек) Ок = О (¿, ек) Нк = Н (п ек).

Эти обыкновенные дифференциальные уравнения Эйлера, как известно [4,5], имеют точное решение при любой правой части, а именно:

ик+1 (г) = Л\г + 4 - + /к (г) г

¿+1 (г ) = В\г + В\1 срк(г ) (14).

г

П+1 (г) = С + Ск2 1п (г) + gk (г) где /к (г ), <рк(г ), gk (г ) вычисляются по формулам:

1 г: 2; 1 г: 2; г : г;

/к(г)=2{1гРкрек)p, <(г)=2{1гОкрек)p, gk(г) = \р1п I- Нкре)р

г1 V г ? г1 V г ? г \Р ?

Константы Аiк, В*, Сiк (i=1,2), Тк, Як определяются из граничных условий (7), (8), которые с учетом определяющих соотношений (6), формул для перемещений и полученных выражений для и,д,ц примут вид:

Ск+1С13 + ЛГ

г г

-вк+1 Л = А (г), -Вк+1-1 = 1\ (г2), СГ1 -1 = А (г), СГ1 -1 = А (2), 2Лк+\ъ + як+Съ = А,

г1 Г2 Г\ Г2

Тк = 1 м

(Г 4 - Г 4 ) г2

-а441 г 0(ек )г

где обозначено:

Л С,,

44

А (г ) = ш(вк) а як+Ак (а + а)+(а - а) 1- сп/к '(г)- си +кшо (вк),

Ьк2

С11 - С12 > г г ) С44 > г )

г - г

г2 { { {

'¡г ш(е к) 2а 3 А! + 0^13

г > > >

Л(г) + ^ | + «33^

+ а0ш

/

ёг

;)

На первой итерации принимаем ш(е)=0, шо (в) = 0, т.е. имеем упругую задачу. Можно показать, что в этом случае на каждой итерации дк = цк = 0.

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

1

= Ак - А-+И(г), 4 = Ак - 4-+

1 + /к

1

4 = „к 4 = 0 4 = 0 4 = ткг

г г 2

где константы Ак, А^, Тк, Як находятся из системы линейных уравнений:

с,Як +

С13 Кк +(с„ + С|2 )Ак +

г

^^ Ак = ь\-

2с^ Ац + С33 ^^ —

Тк = 1

М

п сл

( - г,4) 4

-а

'2

^ гш (ек )г

Обобщенная интенсивность деформаций е (2) в этих уравнениях определяется выражением

{ {{ А 12 { Ак г ;2"1 /- х2

А - А + / + Ак - ^т + /

Л

> >

-Ч 2

г г

?

+ Л (як)2-Лв(^2А1к + / + + Як^ +

+2Л? <2Ак + / + /1 + ЛЯк +1Л4г2 (Тк )21.

Таким образом, поставленная задача А допускает в каждом приближении аналитическое решение (14), которое и используется в пунктах 4, 6 при решении задачи В.

Частный случай изотропного материала получим, полагая:

С1]к1 (Л, м) =: С1111 = С2222 = С3333 = Л + 2м, С1122 = ... = С3311 = Л, С1212 = ... = С3131 = М

а1111 = а2222 = а3333 = 2М (1 - в), а1122 = ... = а3311 = кмв, а1111 = ... = а3333 = 2М

1 2 0 0

Лш = ... = ^3333 = Л212 = ... = ^3131 ГТ 2 , Л111 = ... = ^3333 ,

е, ^/3ek

здесь Л, м - параметры Ламе, е!1 - предел текучести по деформациям.

4

2

2

Тогда определяющие соотношения (2) принимают вид:

ст1 = 2/[1 -о(е)] е1 + К

оА = 2//1 -о(е )] е.

О-\о0

-а0о

записываются

Правые части уравнений равновесия (9) виде:

Р = о

V I

О = о(е)< ¿(г) + -¿(г\ + о'(е)< ¿(г)-^± , Н = о(е)((г) + П(г)) + о'(е)П(г). V г г ? > г ?

Уравнения для определения констант АД В*, С^ (i=1,2), Тк,Як примут вид:

—+1 + 2 — + /Лк+1 + 2/+' -1 = А

-вк+1\ = А (г), - в2к+4 = А (г2), с^ = А (г), С2кЧ = А (г), 2—Г+—+/як+1 = А

Тк _

м

п

(( - г!4 ) 2

'2

2/ г Зо(ек )г

где обозначено:

А

А

А(г) = гг-гт|г о(ек)|-3/Лк -3/\П(г) +

/к -(г4 Я

Л + а°о(ек) + Коо (вк)- — {/'(г) + +

ёг

е2 =

Обобщенная интенсивность деформаций е определяется выражением 1

е

V - V

лк-ч-+/I

)2 : Лк {;2 ^

Л1 ~ +

V г г ? ?

?

+7 (як )2 - 3 3 2Л+/¿+/+я ) +

+2—01 2Лк + /'к+ ^ | + —0Як +

г'(тк )2)•

Полученные результаты проиллюстрируем численным экспериментом, который был проведён по следующей схеме:

1. Считаем, что материал характеризуется одной материальной функцией

о = о(е) (15),

которая считается известной (заданной).

1

2. Решаем задачу А, т.е. стандартную задачу о кручении цилиндра, свойства материала которого описываются функцией ш = ш (е) для ряда последовательных возрастающих значений момента М, что даёт возможность построить зависимость

у = у( М ) (16),

аналогичную той, что получается экспериментально.

3. Считаем зависимость (16) заданной, а зависимость (15) искомой; для нахождения (15) ставится и решается задача В; результат решения - вычисленная зависимость ш = ш (е) сравнивается с исходной зависимостью (16).

4. Расчёты были проведены для модели изотропного материала с линейным упрочнением (17) [3].

\Оу,у<ух

цг) = -

ш

(е ) =

Оух + О' (-у, ),Г>К 0у,у<у,

(О - О') {у-у,) ,У>У, (17).

Оу

На рис. 1, 2 представлены результаты расчётов, которые дают представление о величине расхождения кривых, изображающих зависимость М ~ у, точную, построенную по решению задачи В (нижние кривые), и приближённую, построенную в предположении равномерного распределения напряжений по сечению (верхние кривые). Рассмотрены случаи, когда отношение модуля упрочнения С к модулю сдвига О составляет 0,5; 0,1; 0,01 для сплошного цилиндра (рис. 1) и для отношения радиусов, равного 1/3 (рис. 2).

Отметим, что для слабого упрочнения влияние толщины стенки можно считать малым.

Рис. 1

Сравнение точной и приближённой зависимостей М - у, найденных в численном эксперименте. Внутренний радиус г1 = 0, внешний радиус г2 = 1,7 см. кривые 1, 2, 3 - в предположении равномерного распределения

касательных напряжений по толщине; 1' ,2' ,3' - точное решение, соответственно для О' = О /2; О /10, О /100.

Сравнение точной и приближённой зависимостей М - у, найденных в

численном эксперименте. Внутренний радиус г1 = 1,2, внешний радиус г2 = 1,7

см. Кривые 1, 2, 3 - в предположении равномерного распределения касательных напряжений по толщине;

1' ,2' ,3' - точное решение, соответственно для О' = О /2; О /10, О /100.

Список литературы

1. Кузнецов В.Н., Агахи К.А. Построение материальных функций и численный метод решения краевых задач пластичности с учётом влияния гидростатического давления. Изв. АН Аз.ССР, сер. физ.-тех. и матем. наук, 1976, № 5, с. 130 - 135.

2. Агахи К.А., Кузнецов В.Н., Шестериков С.А. Об учёте нелинейности объёмной деформации в деформационной теории пластичности. Электронный журнал "Исследовано в России", http://zhurnal.ape.ru/articles/2003/014.pdf, с. 142 - 152.

3. Ильюшин А.А. Пластичность (основы общей математической теории). Наука, М., 1963, 316 с.

4. Степанов В.В. Курс дифференциальных уравнений, М., 1964, стр. 211.

5. Камке Э. Справочник по обыкновенным дифференциальным уравнениям, М., 1968, стр. 100, 137, 449.