Добавил:
Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:

книги / Скоростное деформирование конструкционных материалов

..pdf
Скачиваний:
2
Добавлен:
12.11.2023
Размер:
11.58 Mб
Скачать

Рис. 4.29. Сравнение расчетных (сплош­

6, мм

ные линии) данных по глубине фазовых

переходов с экспериментом

 

Зависимость предела текучести о.т давления принимали линейной. Прямоугольниками, размер кото­ рых соответствует точности опреде­ ления бе и 1»0, нанесены глубины зоны постоянного упрочнения бе

взависимости от скорости удара vu. Сравнение экспериментальных

ирасчетных данных показало

хорошее подтверждение предположения о связи упрочнения с фа­ зовым переходом.

В первом приближении значения кинетических коэффициен­ тов /Г2, Л12можно рассматривать как предельные и их уточнение долж­

но

быть связано с уточнением замеров скорости удара

при v0 ^

^

2000 м/с. Расчеты показывают, что более интенсивная

кинетика

перехода (большие значения /Гг или линейная кинетика) приводят или к параллельному смещению кривой б (и0) вверх или к увеличе­ нию ее наклона. Штриховая линия на рис. 4.29 показывает зависи­

мость

ба (vn)

для линейной кинетики с

коэффициентом L12 =

= 6,45

с/м2

и постоянным значением s* =

0,48 ГПа. Уменьшения

наклона кривой б (и0) можно добиться увеличением М, характери­ зующим увеличение предела текучести s* с давлением, и увеличением показателей нелинейности кинетического уравнения /г12. Такое спе­ цифическое влияние упомянутых трудно измеряемых величин суще­ ственно облегчает исследование кинетики.

В заключение остановимся подробнее на следующем обстоятель­ стве. В первых сериях экспериментов для скоростей удара i>„ = = 2700—3000 м/с глубина зоны постоянного упрочнения получа­ лась существенно меньше (на 3—5 мм), чем показанная прямо­ угольниками на. рис. 4.29, т. е. зависимость 6t. (v0) практически вы­ ходила на насыщение по i>0. Теоретические расчеты с различными вариациями /Тг, п12, Д|2, М не давали такого эффекта. Анализ пока­ зал, что это насыщение кривой Ье (v0) связано с двумерными эффек­ тами, а именно определяется боковой разгрузкой второй волны (Т)12), на которой происходит фазовый переход и которая при доста­ точном ослаблении существенно отстает от первой волны (D,). Про­ стые оценки с данными о двухволновой конфигурации возмущения показывают, что волна может потерять плоскую конфигурацию существенно раньше, чем волна £>12, т. е. боковая разгрузка (если не принять учитывающих это обстоятельство мер) может снять про­ цесс фазовых переходов на некоторой глубине, которые не произо­ шли вследствие бокового разгруження и в связи с этим дополнитель­ ного падения давления.

В связи с этим для больших скоростей удара (у0 = 2700—3000 м/с) был увеличен в полтора раза диаметр мишени (d = 130 мм вместо

5*

131

d = 90 мм), ударяющей пластины и заряда ВВ (но при той же их толщине, что оставляло прежней скорость ударника). Это сразу исключило влияние двумерных эффектов на процесс фазовых пере­ ходов / 2 в центре образца и дало заметно более высокие значе­ ния 6е, которые хорошо согласуются с теоретическими результатами.

4.5. УДАР ЦИЛИНДРИЧЕСКОГО БОЙКА О ПРЕГРАДУ

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

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

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

д2ш

ао

с

J

Р Ж

= г

Eez = (Jz - v a o>

е* = ес А

 

Еее = а0- VCL

(4.26)

где w — радиальные смещения; сге, а2 и е0, е2 — окружные и осевые компоненты тензоров напряжений и деформаций. Система (4.26) при г = const имеет решение

- ( sin7 7 T = ^ ) ,' / 1 - v 2 / c ' ] ’

 

V 1-

v2 ,

 

cit

г / Г ^ я ) ) ] 1

Еёpv

 

Clt

V l — V2 Cl

r\f \ — V*’

132

Из полученных решений видно, что для малых значений t напряжен­ ное состояние близко к одноосному деформированному состоянию. С течением времени осевое напряжение осциллирует с частотой

ci!2nrV 1 — v2 и возрастает до напряжения при одноосном напряжен­ ном состоянии.

Будем рассматривать двумерное нестационарное осесимметричное движение сплошной среды в цилиндрической системе координат г, г (от угла, ввиду осевой симметрии, величины не зависят).

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

дгср

д2(р .

1

_д_ /

дц> \

(4.27)

W

~д*

~

~дг V ~ дГ )

 

причем v = grad <р; г — радиальная координата; г — координата вдоль оси цилиндра. Граничные условия запишем в виде:

= vQпри х = 0, (р = 0 при г = 1.

Начальные условия нулевые. Далее будем иметь -дело с безразмер­ ными величинами.

Применим к уравнению (4.27) преобразование Фурье—Бесселя

по г, Фурье по z и Лапласа по t.

В результате получим ряд

 

* = _ v

2v°J° ^ пГ)

 

 

 

sV i (*„) W + s2+<?2) ’

 

 

где Jv — функция Бесселя, a kn — ее корни, причем

 

Ф* (s, q) = j j ф (t, х) e~st cos qx dt dx.

 

(4.28)

о 0

 

 

 

 

 

Так как давление в среде определяется

равенством р = —

то

из соотношения (4.28) получим

 

 

 

 

р* = 2 [ Бф*/dr =

— V -7-т----1-2—

j

(4.29)

j

{till

+

5“ + Я") k7t

 

 

 

 

 

4»о<7

 

(4.30)

о

4

+

.*S + Я ') skn

 

133

Обращая соотношения (4.29) и (4.30), получим

р = -

1>0У \4 r J 0(kn V i2~ x 2) при t > x ,

(4.31)

 

A-J kh

 

 

i

(4.32)

V

- knx j (T2 - Л'2) 2'Л(Лп/ х 2

.Из решений (4.31) и (4.32) следует, что в области течения суще­ ствуют зоны с растягивающими давлениями.

Ниже рассмотрим задачу о соударении цилиндров с учетом проч­ ностных и релаксационных свойств материала. Уравнения движения в цилиндрической системе координат в безразмерной форме имеют вид

диг

_ даг .

darZ

j

gr

 

стф

 

~dt

 

 

дГ

'

dr

 

1

 

г

 

*

dVz

__

д°гг

|

d °z

|

Orz

Q Q \

dt

~

 

dr

+

dz

+

r

 

где

 

vz =

Vz/Ch

t = C,t/X0,

p = p/po,

vr = vrlch

 

r =

r / %

 

2 =

Z/XQ

Or =

Or/PO'h

0(P = Й ф /ро 7.

 

П2 =

dz/poCJ,

Qrz

O rz IP o ^ h

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

V

п •

П|/ = Eijkfikli

где еА, — компоненты тензора

скоростей деформации; стУ/ — кова-

риантная производная по времени, введенная Зарембой и Яуманном, причем

у . d / ч

Gin№mj Ojmtomit

где (&ki — компоненты тензора завихренности. Такой материал часто называют гипоупругим. Для изотропного материала закон Гука для скоростей упругих деформаций запишется в виде

+ atbij.

Запишем полную систему уравнений, описывающую нестационар­ ное осесимметричное движение гипоупругой среды. После некоторых

134

несложных преобразований определяющие уравнения для гипоупругой среды запишем в следующем виде

а7 = ^ + ( 1 - 2 г ) ^ + ( | _ 2 с !) ^ ,

0 У - < 1 - 2 »

где с = cs/ci — отношение поперечной упругой скорости звука к про­ дольной.

Решение системы (4.33) ищут в области В, представляющей собой

объединение двух областей В2 (0 <

г <

R2) 0 <

z <

Нх) и В2 (0 <

с г «

 

R2; Н2 < z <

Их +

# 2) на полуинтервале

] 0; te] и удовле­

творяет

следующим

граничным условиям

 

 

 

 

аг =

аг2 = 0

при

г =

 

0 < z <

Hlt

г = R>,

 

к

Ht -j-

a2 =

ori = 0

при

z = 0,

0 < r <

/?ь

z =

Hi -{- # 2,

0

г < R2.

Начальные данные при t = 0:

 

 

 

 

 

 

 

 

 

vr (0,

M) = ar(0, M) =

o«(0,

/И) =

 

 

 

 

 

 

 

=

cr(p (0, M) = <Jrz (0, M) = 0

 

 

 

для M

£ В,

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

v:(0,M ) = 0, M £B.Z

 

 

 

 

 

 

 

 

vz(0,

M) = i'0> 0 ,

 

 

 

 

 

На линии контакта цилиндров поставим условие прилипания, т. е. смещения при z = Нх + 0, 0 с г < Rt и г = Нх — 0, 0 < г < Rx непрерывны. Следует отметить, что линия г = 0 является для си­

стемы (4.33). особой. На линии

г = 0,

vr = 0, arz = 0.

случаю плоского

соуда­

Аналогично

рения пластин,

бесконечных в попереч­

ном направлении, от плоскости контакта в противоположные стороны (рис. 4.30) распространяются плоские волны сжатия Г и 2'. Совместно с волнами сжатия от границы соударяющихся плоскостей, име­ ющей форму окружности «?, распростра­ няются волны разгрузки. В сечении пло­ скостью, проходящей через ось симметрии образцов, фронт волны разгрузки . имеет форму полукруга с центром на внешней

Рнс. 4.30. Волновая картина в результате соуда­ рения двух цилиндров

135

3

г

1Z)

-1

-J —

)

Х )г

Рис. 4.31.

Линии равного

уровня а2:

Рис. 4.32. Липни равного уровня ог:

1 — ;0.0125;

2 - 0.0118; 3 —^0.0017

/ —.0,0003; 2’— 0,0059; 3\— 0.002

границе

поверхности

контакта

цилиндров. Очевидно, чтс в обла­

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

На рис. 4.30 штриховые горизонтальные линии соответствуют отраженным от свободных поверхностей ударника и мишени волн разгрузки 1 и 2, движущихся навстречу друг другу. Штриховые ли­ нии 3 и 4, 5 на рис. 4.30 соответствуют переднему и задним фронтам боковой волны разгрузки. В области АВС, ограниченной линиями У, 3 и 4 происходит интерференция продольной, отраженной от свобод­ ной поверхности мишени, и боковой волн разгрузки. На рис. 4.31 показаны линии уровня постоянных напряжений а2 в момент вре­ мени t = 0,106, соответствующий началу отражения от свободной поверхности 2 = 0. Из рис. 4.31 следует, что линии уровня а2 = = const являются прямыми линиями вплоть до фронта волны раз­ грузки и симметричны относительно плоскости контакта. Отраже­ ние фронтов сжатия и боковой разгрузки происходит одновременно. В результате интерференции в окрестности угла 2 = 0, z = Rx линии уровня отличаются от линий уровня в окрестности точки пересечения фронта волны сжатия и боковой поверхности пластин. На рис. 4.32 показаны линии уровня аг = const. Отметим, что ни одна линия <jr = const не пересекает боковую поверхность пластин (так как на боковой поверхности выполнено условие а,. = 0), тогда как линии <J2 = const пересекают эту поверхность.

Практический интерес представляет задача об осесимметричном соударении двух цилиндров разного радиуса. Задача решалась ме­

тодом С.

К.

Годунова при следующих

параметрах

г/„

= 0.025

(100 м/с),

У?! =

0,666, Нх = 0,233,

= 1,

Я2 = 0,433.

В

началь­

ные моменты времени от линии пересечения плоскости контакта и боковой поверхности ударника распространяется волна разгрузки, которая примыкает к плоскому фронту сжатия. На рис. 4.33 показана схема волновой картины в'момент времени t = 0,1065, соответствую­ щий началу отражения волны сжатия от плоскости 2 = 0. Здесь: 1 — плоскость контакта; 2 — фронт волны сжатия, распространяю­

щийся вверх; 3 — фронт волны сжатия,

распространяющийся вниз;

4 — фронт волны разгрузки. В области,

ограниченной линиями 2, 3

и4, движение одномерно. На рис. 4.34 показаны линишуровня oz =

=const в тот же момент времени. На рис. 4.34 видно, что линии уровня а2 = const прямолинейны вплоть до фронта тороидальной волны разгрузки, а вблизи линии пересечения плоскости контакта

136

ибоковой поверхности ударника напряжение <т2 резко падает до нуля. Аналогичным образом можно проанализировать и линии уровня сгг = const. Следует заметить, что линия пересечения поверхности контакта и боковой поверхности ударника является особой. Анало­ гом в квазистатике служит задача о штампе. Однако, если в статиче­ ских задачах эта особенность локализована в некоторой малой об­ ласти, то в случае распространения волн эта особенность является источником волн напряжений, существенно, как это видно на рис. 4.33

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

Учет в расчетах, релаксационных свойств среды может привести к существенным изменениям напряжений и деформаций в мишени. В модели С. К. Годунова релаксационные свойства среды учиты­ ваются путем введения нелинейной вязкости. Результаты расчетов по модели С. К. Годунова показывают, что учет релаксационных свойств среды приводит к уменьшению действующих напряжений на 12—15 % [3].

Рассмотрим задачу о продольном ударе упруговязкопластичного цилиндра по абсолютно жесткой преграде. В качестве определяющих соотношений используются соотношения типа соотношений Пэжины (4.10). Интегрирование системы уравнений движения проводи­ лось по модифицированному методу Лакса—Вендроффа вдоль би­ характеристик системы. Впервые подобная модификация для реше­ ния динамических упругих и упругопластических задач с малыми деформациями использовалась Р. Клифтоном. На рис. 4.35 приве­ дены расчетные (сплошные линии) и экспериментальные зависимости

деформаций

е(р от

отношения глубины распространения волны L

к радиусу цилиндра R. В расчетах использовали следующие пара­

метры: v =

0,29;

р =

8 -103 кг/м3; <т0 = 0,156 ГПа,

Я = 1,9 ГПа,

Г = 35 с-1;

р =

5;

k = 0,25. Точность расчетов

обеспечивалась

проверкой интеграла энергии. Относительная погрешность не пре­ вышала 8 %.

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

Рис. 4.33. Соударение цилиндров рал-

Рис. 4.34.

Линин равного уровня о.:

ЦЫХ радиусов ( / — ЛИНИЯ контакта)

/ — 0.01-4;

;>— 0,91; 3 — 0.001

137

Рис. 4.35. Зависимость окружныхостаточных деформаций от глубины

Рис. 4.36. Удар цилиндра по жесткой преграде (штриховыми линиями пока­ зана пластическая зона)

В сторону, противоположную направлению удара, распространяется пластический фронт. Моделирование удара на ЭВМ в рамках квазигидродинамической модели показывает, что пластический фронт, который вызывает торможение упругой части тела, ближе к жесткой границе преграды, чем это показывает профиль ударника (рис. 4.36). На рисунке видно, что пластическое течение происходит и вблизи оси симметрии цилиндра вследствие схлопывания волн боковой раз­ грузки. Вариация отношения длины цилиндра к диаметру в расчетах показали, что отношение окончательной длины цилиндра к началь­ ной при заданной скорости удара является постоянным. Полное время установления состояния покоя при заданной скорости соударения прямо пропорционально начальной длине цилиндра. Расчеты на ЭВМ подтвердили, что начальная скорость пластической деформации линейно зависит от скорости удара. Расчетные скорости деформаций меняются от 2-104 с"1 при скорости удара 0,6 мм/мкс для коротких образцов. Увеличение температуры при высоких скоростях удара приводит к быстрой релаксации напряжения течения. При скорости удара 0,5 мм/мкс увеличение температуры в алюминии в результате пластического деформирования составляет 500—600 К. Следует от­ метить, что для металлов предел текучести для скоростей деформа­ ции свыше некоторого критического значения не зависит от этих скоростей.

Для среды с релаксационными свойствами в декартовой системе координат запишем замкнутую систему уравнений, предложенную С. К-Годуновым [4]:

 

 

 

Qt

(РУ*) "г ~faT(PyAui ■

 

 

 

dekk

,

fekk_

 

dxt

,

9_

dvk

-R n,

(i-

 

dt

 

dXi

-

 

+

ze'ih‘d ^ ’

dtjj .

 

d t j j ___1

/

duj

[

 

j

dv,t

дик

 

!

 

 

2

\

 

'

ddJ \ ,

dt

h

dxk

dxj

dxi

+

b|fc-^7

+ fc^ i i r

138

где Rfj — релаксационные члены, причем

т0 — время релаксации (функция параметров состояния). Уравне­ ние неразрывности, явно не включенное в систему (4.34), является ее следствием.

Запишем уравнение (4.34) в криволинейной еистеме координат. Получим общее выражение для производной по времени компонент тензора второго ранга. Известно, что для любого постоянного век­

тора v = v'e,- =

v'e,- справедливы соотношения

 

 

 

- vV T L ;

- £ = - V

v * ‘.

(4.35)

где

Г),,,, — символы Крмстоффеля. Составим

скалярную

величину

X =

T j/v uv'i в

момент времени

V, а в момент t: К =

7\jV'V, где

Тц — компоненты некоторого тензора второго ранга. Составим раз­ ность

Я ' -

Я =

7

V V

' -

7 >

V

(V=-

I)- § +- | .0 ( ( / ' -1)1)

=

=

|5 -

тк, V.S* -

tu, Y / )

W

(/' -

I) f 0((f' -

/)") =

=

-

( W

" Г™, -

Т.Х Т;,/) п'п' (I

-

I) + 0((/- -

0

) =

 

 

=

L n‘nl v

-

v + 0

-

w -

 

<4-36>

Здесь величины Тц относятся к базису eh а без крышек — к базису е Таким образом, из (4.36) следует

D T u

Dt

в базисе et и

•Ssh II

*u>

'Ъ с>

 

1

 

- T to V jV *

(4.37)

,D T ij

дТи 1

_ т

v kr'!1-

- T imvkГД

(4.38)

Dt

1

1 к1

 

 

в базисе eh

Для тензора деформаций е;,- справедливы соотношения

Тогда для полной производной по времени в базисе et получим

^ = 4 - + VJ®<) - “ *<«V -

139

Учитывая (4.39), уравнения для тензора скоростей деформаций за­ пишутся в виде

 

 

= -г

 

 

 

 

- *».j W " -

 

 

где gij — компоненты метрического тензора.

= 0, i Ф /, gn‘ =

g 93 =

В цилиндрической системе координат g {j

' = 1, g22 = ' г2,

g22g22 =

1.

£ 22б22 = е2.

 

 

 

 

 

 

 

Покажем, что справедливо соотношение

 

 

 

 

 

 

 

оо

4б22

_

4бф

,

о

JV

 

 

 

(4.40)

 

 

S

 

” 5 Г

~

“ З Г

'

^ е<г

Г

 

 

 

 

 

 

Действительно

 

 

 

 

 

 

 

 

 

 

 

 

 

=

" Ж

( ‘Я ' 6*2) +

“ ' 4

(

4

f22)

+

1'2 4

( 4

Е22) +

 

+

 

= 4

^

)

+

2г 22<!22 4

=

 

+

2f„

 

Тогда для

еф получим

 

 

 

 

 

 

 

 

 

 

 

 

 

 

= (I -

2е„,)

 

 

ч

 

 

 

 

^ > 3.

(4.41)

Аналогично получаются соотношения для ег и ez.

Окончательно запишем полную систему уравнений (безразмерных) для нестационарного движения вязкоупругой среды с цилиндриче­ ской симметрией

4 (рv )+4 (PW)+4 (р»я= 4 (a'S )+4-

dt е^ иг ~ д Г ^ - ^ - д Г -

= П _

9с ^

_

о»

вг — (ег + е2 + вф)/3

 

11

 

дг

Je~dT

Г0------ »

 

А р

_1_ 7,

A L

(1 - 2 e ,) 4 & ~

 

 

 

 

 

 

_ 9с f r y ,

вг - ( е г + Ег -4еа,)/3

(4.42)

 

 

 

 

 

- ^ +

 

°

(1 - % ) Т~ - **~ (РЦ * -+ Кф)/3 ,

 

 

д(Р')

f ~Qf (р™г) +

(Р"Уг) — 0.

 

 

5/

 

140