книги / Скоростное деформирование конструкционных материалов
..pdfРис. 1.44. Зависимость давлении от напряжения, при ложенного к конденсатору
экспериментальные зависимости деформации от времени для начальных скоростей дефор мации в диапазоне 330—2030 с-1. Отме чается, что кривые, соответствующие на чальным скоростям деформации, превыша ющим ё0 = 2000 с-1, соответствуют состоя нию материала, предшествующему разру шению.
Среди прогрессивных методов технологии обработки металлов давлением использова ние удара занимает важное место. Исследова
ния высокоскоростного деформирования заготовок ведутся во многих направлениях. Широкое распространение получил метод дефор мации листовых заготовок электрогидравлическим методом через передающую среду. Одна из схем деформации заготовки ударной волной через воду представлена на рис. 1.43.
Между двумя электродами в течение короткого интервала времени возникает высокая температура. Жидкость превращается в пар, который быстро расширяется и создает ударную волну. Удар ная волна, взаимодействуя с заготовкой, деформируется в матрицу.
Наряду с электрогидравлическим методом для высокоскоростной деформации трубчатых образцов используется энергия магнитного поля. Трубчатый образец помещают внутрь катушки соленоида. Энергия, накопленная в конденсаторах, разряжается в течение короткого промежутка времени через катушку-соленоид. В резуль тате образуется магнитное поле, развивающее давление на заготовку.
Давление магнитного поля определяют по формуле |
Р = |
х |
||||
X ( й о о У 1 где |
— коэффициент; (х0 — магнитная |
проницаемость |
||||
материала; |
Нм — напряженность магнитного поля. Напряженность |
|||||
магнитного |
поля определяют по формуле Ни = ( “^ |
) '/~“ |
> гДе |
|||
k0 — постоянный |
коэффициент; с — емкость конденсаторной |
бата |
||||
реи; и — напряжение на обкладках |
конденсатора; |
I — длина |
||||
катушки; г — радиус катушки. На рис. |
1.44 приведена зависимость |
давления, развиваемого магнитным полем, от напряжения, при ложенного к конденсатору.
Наиболее распространенным методом регистрации высокоскоро стной деформации оболочек является метод высокочастотной ре гистрации. Необходимую для съемки частоту кадров определяют
по формуле v = 0,7vr$N, где vr — радиальная скорость |
оболочки; |
||
Р — масштаб съемки; N — разрешающая |
способность оптики. Для |
||
обычных условий съемки v « |
10е кадров в секунду. Следует отме |
||
тить, что важной проблемой в |
методах |
высокочастотной |
регистра |
ции, как и в любых других оптических методах, является |
проблема |
||
защиты оптики от осколков разрушенной оболочки. |
|
31
Г ла ва 2 ОПРЕДЕЛЯЮЩИЕ СООТНОШЕНИЯ ДИНАМИЧЕСКОГО ПОВЕДЕНИЯ МЕТАЛЛОВ
2.1. ДИНАМИЧЕСКИЕ ЗАВИСИМОСТИ НАПРЯЖЕНИЕ - ДЕФОРМАЦИЯ
Рассмотрим зависимости между напряжениями и деформациями для случая одноосного напряженного состояния (стержень). В ра боте [21 содержится обзор исследований, посвященных этому во просу. Авторы обзора отдают предпочтение функциональному под ходу, т. е. заданию связи между напряжениями и деформациями в форме функционала
(0 = F[T',
структура которого определяется из экспериментов. Сложные струк турные механизмы, ответственные за неупругое поведение материала, при таком подходе учитываются интегрально через структуру дина мического функционала пластичности. В настоящей работе разви вается подход, основанный на использовании кинетических (эволю ционных) уравнений, описывающих структурные изменения мате риала. Следует отметить, что во многих случаях кинетический и функциональный подходы эквивалентны.
Теория распространения упругопластических волн в материалах, напряжения в которых зависят от деформации, но не от их скорости, развивалась X. А. Рахматулиным 117] и другими отечественными и зарубежными авторами. В литературе эту теорию часто называют
теорией |
Рахматулина — Кармана — Тейлора. |
В |
ее основе |
лежит |
||||||
предположение о |
существовании |
однозначной |
зависимости |
ах = |
||||||
= / (ei)> |
что позволяет |
свести задачу |
о |
распространении |
волны |
|||||
в стержне к решению волнового уравнения |
|
|
|
|
||||||
|
дЧ |
_ |
2 |
дЧ |
- |
V |
1 |
dgi |
( 2. 1) |
|
|
дР |
~ |
С[ |
дха ’ |
||||||
|
! |
V |
р |
de2 |
|
где и — смещение частиц среды; ct — продольная скорость волны. Теория распространения упругопластических волн в материалах, напряжения в которых не зависят от скоростей деформаций, интен сивно развивалась для исследования одномерных динамических волновых задач, когда реализуется двухосное напряженное или деформированное состояние, а также для исследования распростра нения двухмерных волн. Благодаря своей простоте и доступности эта теория широко применяется для расчетов элементов конструкций
при кратковременных импульсных воздействиях.
Однако многочисленные экспериментальные данные по динами ческому деформированию многих конструкционных материалов
32
указывают на необходимость учета зависимости деформаций от вре мени. Особенно явно эта необходимость проявляется при обработке данных экспериментов, в которых изучаются переходные процессы:
нагружение — разгрузка, |
нагружение — ступенчатое дополни |
тельное нагружение и т. |
д. |
В рассматриваемых в предыдущей главе экспериментальных данных Дж. Липкина и Б. М. Малышева исследуется распростране ние дополнительных импульсов в предварительно пластически де формированном (статически или динамически) материале. Во всех указанных экспериментах скорость распространения добавочного импульса равнялась скорости упругой волны, а не локальной ско
рости (Ci = V ^Oj/prfej), как это следует из теории Рахматулина — Кармана — Тейлора. Кроме того, экспериментальные данные по изучению распространения догрузочных волн высокой интенсив ности показывают, что вторичная волновая картина воспроизводит первоначальную, т. е. впереди распространяется упругий пред вестник, а за ним — пластический фронт.
Результаты описанных выше экспериментов не укладываются в рамки теории Рахматулина — Кармана — Тейлора. Указанные расхождения являются следствием того, что в теории не учитывается влияние скорости деформации материала на процесс образования волн.
К одним из первых работ, учитывающих чувствительность мате риалов к скоростям деформации, следует отнести работы А. А. Ильюшина [6, 7]. Обобщением определяющего уравнения, пред ложенного А. А. Ильюшиным, служит модель вязкоупругого тела
Р п (~ д Г ) ai = ("аГ) 8ь (2 ^
где Я„ и Qn — полиномиальные дифференциальные операторы сте пеней п и т соответственно,
|
Рп |
|
|
|
b |
д> |
|
|
|
|
|
|
|
||
|
|
|
|
|
/= О } dt> |
' |
|
Заметим, что максимальное значение степени |
п зависит от /я, |
||||||
с тем, чтобы избежать противоречий со 2-м законом Ньютона. |
|||||||
Уравнения |
движения для |
среды |
|
|
|
||
|
д2и |
_ |
до1 |
д&х |
ди |
|
(2.3) |
|
Р “д7* |
|
дх |
1 dt |
~~ дх |
|
|
|
|
|
|
||||
совместно с равенством (2.2) приводят |
к уравнениям |
||||||
|
-%г 1Рп (о.)1 = |
-jr - |г |
К » (о.)]. |
|
(2.4) |
||
С помощью |
преобразования |
Лапласа |
|
|
|
||
|
/(/>) = |
|
1 /( 0 *-'’' <и |
|
(2.5) |
||
|
|
|
о |
|
|
|
|
2 МоПборода В. Г1. ■■ДР• |
33 |
Рис. 2.1. Эпюры напряжений в стержне моделей Максвелла (я), Фойхта (б) и че тырехпараметрической (в) для разных моментов времени
уравнение (2.2) приводится к виду |
|
||
Л. 0>) о. = |
(/>)•!£--£-. |
(2.6) |
|
где Рп(р) = 2 aiPl+2, |
Q (/?) = |
2 bjpl. |
|
i= 0 |
|
/=0 |
|
Решение уравнения (2.6) в изображениях по Лапласу таково: |
|
||
6 = С, ехр [— / Р |
п (p)!Qm(р) х] + Сг exp [ / Р л (p)/Q„, (р) х), |
(2.7) |
где константы Q и С2 определяются из граничных условий задачи. Оригинал изображения (2.7) сравнительно просто находится для следующих моделей:
1) модель Максвелла: До<*1 |
-7^ |
= &1 -75J-; |
|||||
2) модель Фойхта: a0ai = |
M i + |
&i |
; |
|
|||
3) четырехпараметрическая модель |
|
|
|||||
„ |
^ |
л |
d oi |
д,2аг |
и |
дег . |
, а2е! |
W |
i + |
a i - g f |
+ а -г -ы г = |
bi |
“ а Г + |
а/2 * |
На рис. 2.1, а—в показаны эпюры напряжений в стержне в раз ные моменты времени соответственно для каждой из моделей.
Можно показать, что система уравнений (2.3) совместно с опре деляющим уравнением Максвелла имеет две действительных харак теристики, т. е. является системой уравнений с частными производ ными гиперболического типа. В случае использования моделей Фойхта и четырехпараметрической модели вязкоупругого тела тип уравнений меняется на параболический. Это видно и из решений, представленных на рис. 2.1.
В табл. 2.1 приведены значения параметров, используемых при построении решений.
В. В. Соколовский и Л. Мальверн постулировали существование статической кривой св = / (ег) и вязкие деформации рассматривали
34
Таблица 2.1
Модель |
п0-ю-«. |
a i |
а 2. Ю " |
Ьв. |
V 10\ |
6,- Ю -*. |
|
С-1 |
ГПа |
ГПл |
НДмм-с) |
||
Максвелла |
0,02 |
1,0 |
0 |
0 |
0,02 |
0 |
Фойхта |
1,2 |
0 |
0 |
1,44 |
0,012 |
0 |
Четырехпараметрнче- |
411 |
197 |
2,32 |
0 |
0,0021 |
0,488 |
ская |
|
|
|
|
|
|
как меру превышения действующих напряжений над статическими. Определяющие соотношения для среды записывали в следующем виде
dajdt = Edejdt — [(ах — ае)//с] (В. В. Соколовский)
dojdt = Edejdt — f[a1 — <JS (E^] (Л. Мальверн),
где |
as |
(Ej) — статическая |
кривая зависимости |
напряжение — де |
|||||
формация; |
ое — статический предел |
текучести материала; |
tc — |
||||||
эмпирическая постоянная. |
|
|
|
|
|
||||
|
Определяющее уравнение типа Л. Мальверна часто используется |
||||||||
в следующей форме |
|
|
|
|
|
||||
|
|
|
dl - £4, = |
-g- [exp |
- |
1 ] • |
(2.9) |
||
В |
соотношение |
(2.9) входят |
два |
неизвестных параметра |
сгс |
||||
и |
tc, |
а |
также |
зависимость |
crs (ех) — статическая диаграмма. |
Нетрудно показать, что с помощью интегрирующего множителя
exp |
\t/tc + |
[Еех — os (е^ ]/<тс} |
дифференциальная |
связь между |
||||
напряжениями и деформациями сводится к функционалу вида |
||||||||
|
|
<*i = |
<*stei) - |
<JCIn [exp ( — |
*~*л ) exp ( — £ei~ gs ) X |
|||
|
|
( ' — |
, р = а М ] ) + . |
) |
x |
|||
- |
r |
[ - |
^ |
“ |
] |
|
( г - |
(2.10) |
|
|
|
|
|
|
|
|
|
где |
of |
= |
Ее1/ при |
t = td, ef = |
Ei (/'). |
|
Таким образом видно, что достаточно простое дифференциальное уравнение (2.9) приводит к функционалу со сложной структурой. При этом следует иметь в виду, что для адекватного описания экспе риментальных данных в широком диапазоне изменения скоростей деформаций правая часть уравнения (2.9) будет иметь вид более сложный, а это в свою очередь приведет к существенному усложне нию структуры динамического функционала пластичности.
Часто при решении динамических задач априори ограничиваются без достаточных причин линейной частью функционала в его раз-
2* |
35 |
ложении в ряд по полилинейным функционалам. Однако для задач динамического неупругого деформирования металлов, когда суще ственна нелинейность вязких и пластических свойств, такой подход неоправдан и возникает проблема создания обоснованной, надежной методики определения структуры динамического функционала пла стичности. Для решения волновых задач такая методика, по-види мому, отсутствует.
С точки зрения числовой реализации задач нестационарного движения среды с волнообразованием, на наш взгляд, удобнее использовать определяющие уравнения в форме неголономных диф ференциальных соотношений типа (2.9). Прежде всего это связано с тем, что уравнения сохранения — неразрывности, импульса и энергии, используются обычно в дифференциальной форме.
В том случае, когда в правую часть определяющего уравнения (2.9) входят некоторые параметры накопления, для которых имеются свои дифференциальные уравнения, представить такое определя ющее уравнение в виде динамического функционала пластичности вида (2.10) вообще говоря не удается.
Н. Кристеску обобщил соотношения (2.8) для случая мгновенной пластичности, когда мгновенная реакция материала на динамическое воздействие неупругая. В соответствии с этим предположением определяющее уравнение для такого материала записывают в виде
di = (p(<Ti, Е1) 81 + ф(а1, |
е:). |
(2.11) |
|||||
Запишем замкнутую систему уравнений движения среды с учетом |
|||||||
(2.11) и исследуем характеристики системы |
|
|
|||||
dv |
__ _1_ |
до\ |
dv |
_ |
д&х |
(2. 12) |
|
dt |
~~~р~ |
дх * |
дх |
~~ ~дГ * |
|||
|
|||||||
4 г |
= Ф(®1. *i).-fr + |
(а*> е‘)- |
|
Для того, чтобы выявить тип системы (2.12), дополним ее следу ющими соотношениями
d o ^ ^ - d x + ^ - d t , |
|
+ |
(2.13) |
‘Ь - Т Г ' Ь + Т Г * -
Расширенную систему (2.12) и (2.13) разрешим относительно производных
дах _ |
ф (daxdt/p) — dvdx — |
dx2 |
|
dt |
~ |
(fdt2/p — dxz |
|
дгх |
|
dojdt/p — dvdx — tydt2/p |
|
dt |
~~ |
qdl2/p — dx2 |
’ |
36
dej _ Idxdv—d^dl/p •+-i|xtt*/pl + [q>dP/p — dx*\ de, |
(2.14) |
||
|
|
[*prf/*/P — <i»aJ |
|
|
|
|
|
dv |
_ |
daxdtlp — i|)d/2/P— dvdx |
|
dx |
~ |
ydP/p —dx* |
|
Из этих уравнений следует, что система уравнений (2.12) обладает
тремя семействами |
характеристик |
|
|
|
|
|
т г |
= |
± |
= |
± с, (О,, е,), |
(2.15) |
|
|
|
= |
0, |
|
|
(2.16) |
причем на характеристиках (2.15) имеют место зависимости |
|
|||||
|
= |
±рс, (CTL ех) rfu + |
ф (ab |
ex) d/, |
(2.17) |
|
a на (2.16) — зависимость |
|
|
|
|
||
|
da\ |
= cp (<jlf ej) det + |
ф (alt |
e j dt. |
(2.18) |
Следовательно, возмущение, возникающее на конце стержня, будет
распространяться как волна со скоростью |
± |
ф ^ |
£l^ ■, отличной |
||
от |
звуковой. Очевидно, что в случае |
Соколовского — Мальверна |
|||
Ф |
(o'!, |
= Е, и скорость распространения |
головного возмущения |
||
будет |
равна скорости звука ct = ± У Е/р. |
Далее |
на конкретных |
числовых расчетах будет показано, что использование соотношений типа (2.11) не описывает ряда экспериментальных данных по рас пространению волн в стержнях из конструкционных материалов. Не описываются, например, явления запаздывания текучести и затухания амплитуды упругого предвестника, а также наличие «плато» пластических деформаций вблизи торца ударяющегося стержня.
В отличие от модели Соколовского — Мальверна, в которой по стулируется существование статической диаграммы as = as (ег), в модели материала с наследственными свойствами постулируется
существование динамической диаграммы ad = cprf (еГ), а последу ющее поведение материала определяется «сползанием» по этой диаг
рамме |
в соответствии |
с функциональным равенством |
|
|
|
t |
|
|
G\ = |
фdМ ) - f К (t - s) CTi (s) ds, |
(2.19) |
|
|
6 |
|
где К |
(t — s) — разностное ядро релаксации, определяющее эффект |
«памяти» материала.
Многочисленные экспериментальные данные по скоростному де формированию конструкционных металлических материалов опи сываются в рамках модели, основанной на дислокационных пред ставлениях [13, 14]. Структура уравнений этой модели имеет вполне определенное физическое обоснование. Тем не менее эту модель
37
Рис. 2.2. Зависимость скорости деформации от деформации
нельзя назвать физической в том смысле, что все фигурирующие в ней параметры определяются из независи мых физических теорий или экспери ментов. Эти параметры определяются из решения обратной задачи и сравне ния с данными механических экспери
ментов по методике, описанной ниже. Таким образом, речь идет о совместном физико-механическом подходе к построению опреде ляющих соотношений.
Определяющее соотношение для скорости неупругих деформаций часто используется в форме Гилмана— Джонстона,
*? = -Ьс, (п, + -§-те?) ехр [ - 2 ( D + - f Лв?/<г,) ] , (2.20)
где п0 — начальная п л о т н о с т ь дислокаций; Ъ — вектор Бюргерса; т — коэффициент размножения (дислокаций); Н — постоянная упрочнения; D — эмпирическая постоянная. Нетрудно видеть, что с ростом неупругих деформаций е'1скорость неупругих деформа ций (при некотором фиксированном at = const) падает после дости жения некоторого критического значения е" = . На рис. 2.2
приведены зависимости ёу от еу для разных уровней напряжений. Физический смысл этого свойства тот, что при росте общей плотности дислокаций «запираются» подвижные дислокации и таким образом происходит упрочнение. Ниже на конкретных числовых расчетах будет показано, что с использованием определяющего соотношения (2.20) удается описать как затухание амплитуды упругого пред вестника, так и образование «плато» постоянных пластических де формаций вблизи соударяющихся поверхностей пластин или стержней.
2.2. ОПРЕДЕЛЯЮ Щ ИЕ СООТНОШЕНИЯ П Р И СЛОЖНОМ Н А Г Р У Ж Е Н И И
Обобщение определяющих соотношений для среды в случае одноосного напряженного состояния (стержень) на случай сложного напряженного состояния представляет собой важную научную и практическую проблему. Исследования такого рода необходимы прежде всего для дальнейшего совершенствования методов расчета конструкций и технологических процессов. При функциональном методе описания динамической пластичности важную исходную информацию о свойствах динамического функционала пластичности могут дать постулат изотропии и принцип запаздывания, предло женные А. А. Ильюшиным.
Ниже используется принцип экстремальности реального про цесса, приводящий к условию ортогональности вектора скорости
неупругих деформаций к поверхности нагружения в пространстве напряжений. Для случая сложного нагружения строятся определя ющие соотношения. При этом в число определяющих параметров среды включаются параметры, характеризующие изменения ее вну тренней структуры. Для этих параметров записываются кинетические (эволюционные) уравнения.
Непосредственным обобщением соотношений Соколовского — Мальверна (2.8) являются соотношения Пэжины для расчета дина мических задач при сложном напряженном состоянии. Приведем один из вариантов определяющих уравнений типа уравнений Пэ
жины: |
|
ёц = -Щ- |
для F < 0. |
= |
(2.21) |
где / 2 = SijSij, G, /С, ц, <7j), Si], Ьц, ё,], F — соответственно модули упругого сдвига и объемного сжатия, коэффициент вязкости, компо ненты тензора напряжений, компоненты девиатора тензора напря жений, компоненты тензора скоростей деформаций, компоненты девиатора тензора скоростей деформаций, F — поверхность теку чести. Для функций Ф (F) используются следующие зависимости:
Ф(Г) = (.
Ф(Р) = (Ф-'Г-
Ф (F) = |
1 |
.Ъ |
|
|
k = }/(oi + 2HWp)/3t
k = (сг. -f- Яб;) j Т,
А = (а, + Я е?)//3 .
где ое, k, Н, Wp, е? — соответственно пределы текучести при простом
растяжении и сдвиге, постоянная упрочнения, работа на пласти ческих деформациях, интенсивность пластических деформаций.
Последовательный термодинамический подход к построению полной системы уравнений движения сплошной среды с учетом релаксации напряжений развивается в работах С. К. Годунова [4]. Вводится удельная внутренняя энергия как функция параметров состояния среды: г,] — компонент тензора эффективных деформаций
иs — энтропии. Связь между компонентами тензора напряжений а1'
итензора ги задается формулами Мурнагана, которые в произволь ной системе координат имеют вид:
<Ji'' = p (6 fi-2 e!i) Д - , |
(2.23) |
39
где р = |
Ро (det||6fJ- — 2е,7 ||)|/2; |
р0 — начальное |
значение |
плот |
ности; 8 |
— символ Кронекера; |
0 (ijj, s) — удельная внутренняя |
||
энергия |
среды; s — удельная энтропия. |
вводится |
метри |
|
Наряду с тензором эффективной деформации |
ческий тензор эффективной упругой деформации gij. Предполагается, что мгновенная реакция среды упругая, а уравнение состояния при нимается в виде U = U (klt k3t s), где klt k2t k3 — коэффициенты удлинений в направлении главных осей тензора деформаций. Фор мулы Мурнагана в главной системе координат
Ро |
W . |
= Л®_ |
9 U . |
||
k.,ka |
дку ’ |
2 Мз |
дк« ’ |
||
|
Ро |
dU , |
T _ d U _ |
(2.24) |
|
3 ' кхк« |
дк3 ’ |
ds |
|||
|
где Т — абсолютная температура.
В этом случае для главных компонент метрического тензора эффективной упругой деформации можно написать gi = 1/jfe? (i = = 1,2, 3). Компоненты тензора эффективной деформации опреде ляются формулами ё,7 = (6,7 — 2gij)l2.
Полная система уравнений движения среды с учетом релаксации, предложенная С. К. Годуновым, запишется в следующем виде:
|
|
~д1 |
|
д+х ^ (FW m ) = |
°Л>»> |
|
|||
|- ( р \и + |
ЭД./2]) + |
-Щ- И |
(U + вд,/2)] = |
(2.25) |
|||||
deki |
I |
. |
deki ____ 1 |
/ |
dvh |
. дщ \ , |
dvm |
, |
|
dt |
+ Vm |
dxm |
2 |
\ |
dxi |
dxh ) + |
eAm"faT + |
||
|
+ |
e'"t'S T |
= |
_ |
I f |
(**' ~ "T" 6|“ ) ’ |
|
||
|
> |
+ |
|
+ |
|
^ |
- 0 . |
|
Совместно с формулами Мурнагана (2.23) система (2.25) является полной. Заметим, что в систему (2.25) не включены уравнения не разрывности и уравнение для энтропии 5. Можно показать, что при выполнении условия
|
|
|
= ° |
(2-Щ |
уравнение неразрывности является следствием этой системы. |
||||
При выполнении условий |
|
|
||
Ь д и /д Ь -Ь д и /д к г ^ |
Л |
Ш / д к г - к 3ди/дк3 ^ |
Л |
|
щ - щ |
^ |
и’ |
>и> |
|
|
k3dU/dk3- k JU /дк! ^ л |
|
||
|
|
Ж=щ. |
> и |
|
система (2.25) является |
гиперболической. |
|
40