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

Моделирование связанных процессов в реагирующих средах. Монография

.pdf
Скачиваний:
0
Добавлен:
07.09.2026
Размер:
2 Мб
Скачать
☆
2.3.Механическаямодификацияпорошковойсреды

Выполнение пунктов а — в в модели Нестеренко осущест­вляется в рамках сферически-симметричной модели (рис. 3, в). Центральная часть составной модельной сферы с радиусом c представляет собой недеформируемую область частицы по­рошка на стадии процесса схлопывания поры.
Геометрические параметры модели единичной ячейки Не­стеренко имеют вид
3
'

ddd
''
b α ac αα α
33

000 0 0
;; 2; 1,
222
a
n

'
a
0


(2.10)
где d — диаметр частицы.
Для описания зависимостей предела текучести материала от температуры использовалось соотношение
1
σ TT T T

Tmm
где σ
— характерное значение предела текучести материала
Т1
σ
T
при низких температурах; T
1
0
— температура плавления.
m
TT
,
m
Зависимость вязкости материала от температуры принима­лась в виде

— вязкость расплава.
где η
m
ηη B

exp ,
m
11

TT

m
В краевой задаче (2.3) — (2.8) для модельной составной сферы вместо соотношений (2.7) использовались начальные условия:
a'(0) = a'
aTrT
0
n

00
'
;
00; ,0 .
Выражения (2.8) были заменены выражениями для геометри­ческих параметров модели единичной ячейки Нестеренко, полу­ченными из условий а — в для случая плотной упаковки (2.10).
С применением зависимостей σ
и η от температуры были
T
получены решения краевой задачи, представляющие опреде­ляющие соотношения пористой среды [330]:
81
2.Физико‐математическаямодель
'
b
n

  


Pt σηB
  
ρ cT σηB
ss T m
21 6exp

ρ aa a a
snn n n

'

a

n

' ' '2 '2

(2)1 1;
 

  
21 12exp ;

  
1
  
''33 33'3

baaα rraa
 
nn n

'*
ac αα α

n
Trr
  
Tm
1
TTTrr
mm
  
'

a
n
 
Tr r
Tr TT r
mm
'1; ' ;
00 0 0
d
3
2; .
2
b


00

11


'4

1
'
2
n

a
n

'4
b
n

11
 

13
 

 

13
2
(2.11)
Важно отметить, что в модели Нестеренко в качестве ха­рактерного размера для несвязных порошков вводится размер частиц, а не поры, как в модели Кэрролла — Холта, что позво­ляет адекватно моделировать реальную порошковую смесь, в которой при одних и тех же размерах частиц и пористости ре­альная величина пор может быть различной. Модель единич­ной ячейки Нестеренко учитывает зависимость теплофизиче­ских свойств от температуры и позволяет принимать во вни­мание неоднородность пластического деформирования частиц (в частности, локализацию интенсивных пластических дефор­маций по границам порошковых частиц). В модели Кэрролла — Холта неучет зависимости теплофизических характеристик от температуры приводит
к невозможности адекватно рас­сматривать локализацию зон повышенного тепловыделения, возникновения струйных течений в поверхностных зонах час­тиц. Анализ этих факторов необходим для адекватной оценки реакционной способности реагирующей порошковой смеси в локальных микрообъемах среды, для учета влияния теплового эффекта на изменения механических свойств материала и ус­ловия инициирования механохимических превращений.
Экспериментально
обнаружено, что при ударном нагружении пористого образца возможны существенно непохожие меха­низмы деформирования [291]. В одном случае деформация час-
82
2.3.Механическаямодификацияпорошковойсреды
s
тиц имеет характер местного смятия контактирующих поверхно­стей, в другом — спокойного общего течения, в третьем — рез­кого расплескивания, в основном локализующегося в поверхно­стных слоях частицы, прилегающих к поре. В реальных порош­ковых образцах перечисленные процессы могут реализовываться одновременно в различных зонах. Это определяет существенную разницу динамического нагружения и статического. В
работе
[64] показано, что в отличие от квазистатического при ударном
нагружении порошков в так называемом динамическом режиме поверхность частиц очищается, и окислы концентрируются в от­дельных карманах — исходных порах. При этом в динамическом режиме наблюдается разрушение окисных пленок уже в перед­ней части фронта ударного импульса, и дальнейшая деформация происходит при
частично очищенных поверхностях.
С использованием системы (2.11) в модифицированной
модели может быть найдена удельная энергия W
, диссипи-
д
руемая в окрестности поры на стадии ее схлопывания [64]:
'
b
n
3
c
WTTrdr
v

д

ρ
'
s
a
n
2
0
.
Когда можно пренебречь процессами вязкого течения,
вводится среднее для процесса значения предела текучести
и выражение для Wд запишется в виде
σ
Tcp
2
σ
WHxHxHxx
xx
01
Tcp
 
Д
3
1

ΠΠΠ
010
  
ρ
s
;.
01
1
Π
1

,ln1,
(2.12)
0
Если же основной вклад в уплотнение вносят процессы вязкого затекания пор, то после введения понятия средней вязкости η' получено следующее выражение [64]:
'
b
n
4' ' 2'
W
ηαdαηηP

д
31
ρααα ρt ρ d ρ


'
a
n
0
ss s
. (2.13)
83
2.Физико‐математическаямодель
I
Разница между всей запасенной удельной энергией сжатия W и удельной энергией, диссипируемой в окрестности сфери­ческой поры на стадии ее схлопывания W микрокинетическую энергию: W
= W – Wд. Она расходуется
k
, представляет собой
д
на совершение работы очищения поверхностных контактных слоев частиц и, возможно, уплотнения порошкового материала за счет того, что в поверхностных слоях частиц возникает жидкоподобное поведение материала [64]. С увеличением им­пульса приложенной нагрузки доля диссипированной энергии, потраченной на пластическую деформацию и вязкое затекание пор, в общей энергии ударного перехода уменьшается,
а зна­чит, растет доля энергии ударного импульса, которая может быть потрачена на разрушение поверхностных слоев частиц.
Работу разрушения поверхностных слоев частиц можно
оценить по формуле Ирвина [292]:
212
A πνEK
1
 , (2.14)
*
C
где А
— работа увеличения площади трещины на единицу;
*
— предельное значение коэффициента интенсивности
K
IС
напряжений — трещиностойкость материала.
Таким образом, для увеличения площади трещины от зна­чения S меньшую чем 2А(S
до значения S1 необходимо затратить энергию, не
0
— S0). Если при разрушении поверхности
1
отдельной частицы площадь вновь образованной поверхности увеличивается от 0 до площади поверхности частицы S
, то
ч
энергия, необходимая для «очистки» всей поверхности одной частицы, будет равна
WKS
oi I ч
21
2
πν
E
2
.
Полная энергия, требуемая для разрушения поверхностных слоев со всех частиц слоя в этом случае примет вид
M
21
i
πν
1
WKS
oIiчi
2
i
2
E
i
,
где M — число частиц в слое реакционной ячейки.
84
2.3.Механическаямодификацияпорошковойсреды
5
В работе [51] предлагается шкала давлений, характери­зующая специфическую реакцию порошкового тела на им­пульсное воздействие. Шкала представляет собой последова­тельность формул, позволяющих вычислять для каждого ма­териала ряд пороговых давлений ударного импульса, соответ­ствующих переходу в качественно различные варианты мак­ротечения, режима деформирования и нагрева частиц:
1. Давление, ниже которого невозможен
ударный переход
(по Кэрроллу — Холту):

— предел текучести;
где σ
T
ρ
— начальная плотность порошка;
0
— плотность монолитного материала.
ρ
s
PP σ

sl T
min
2
3
ρ
s
ln ,

ρρ
s

0
2. Нижняя граница реализации «струйного» режима уп-
лотнения порошка (локализации деформации в поверхностных слоях частиц) [52]:
ρ
Hv
0
0,36 ,
P
s
2
1
ρρρ

s
0
где Hv = 2,58σ
— предел прочности.
σ
в
Пороговое давление P
— твердость материала частиц;
в
и плотность ρ за фронтом ударного
s2
импульса связаны ударной адиабатой по Кэрроллу — Холту:

2
P σ
sT
2
3
ρ
ln ,
s

ρρ
s

причем эти два уравнения нужно решать совместно. С учетом того что в рассматриваемой порошковой среде возможны хи­мические превращения, шкала последовательности критиче­ских уровней давлений для компонентов порошковой смеси определяет условия повышения реакционной способности реагирующих компонентов механической активацией.
8
2.Физико‐математическаямодель
6
2.4.Критерийинициирования«турбулентного»режима динамическогоуплотнения
 
В процессе динамического нагружения в локальных мик­рообъемах среды (на мезоуровне в некоторых подобластях среды конечного размера, являющихся внутренними для пред­ставительного объема исследуемого материала) происходит зарождение дислокационной структуры, полос локализованно­го сдвига, вихревых структур и других дефектов, а также раз­рушение межкристаллических и поверхностных слоев, опре­деляющих механическую модификацию порошкового риала и активацию реагирующих компонентов.
Е. Эйкинс и Н. Н. Тадани [293] рассмотрели возмож-
Д. ность реализации сверхбыстрых превращений в ударнонаг­руженной порошковой системе Ni-Al в зависимости от кон­фигурации исходных порошков — формы, размера, распре­деления. Экспериментально обосновано, что одно из важ­ных условий реализации сверхбыстрых превращений — по­явление «турбулентного» течения, способствующего интен­сивному перемешиванию. Этот постулат хорошо согласует­ся с заключениями, сделанными
Я. Хори, который механиз­мом перемешивания, необходимым для реакции, называл образование микроструктуры и разрушение вещества в вих­рях (приводится по работе [293]). Межчастичный сдвиг, возникающий при сплющивании и поперечном течении час­тиц, подобен механизму, предложенному в роликовой мо­дели А.
Н. Дремина [134].
«Турбулентный» режим уплотнения, характеризуемый
развитием вихревых структур и поперечным течением час­тиц, назван Ю.
И. Мещеряковым «катастрофическим» [202].
Автор доказал, что при определенных условиях может про­изойти смена однородного режима деформирования на ре-
мате-
8
2.4.Критерийинициирования«турбулентного»режимадинамическогоуплотнения
7
жим «катастрофический», означающая смену режима энерго­обмена между мезо- и макроуровнями деформирования [202]. При переходе к «катастрофическому» режиму деформирова­ния динамические характеристики среды претерпевают рез­кое пороговое изменение, что приводит к формированию субструктуры материала (фрагментации). Предложенный
И. Мещеряковым критерий перехода от однородного ди-
Ю. намического деформирования к «катастрофическому» [202] получен на основе статистического подхода. Рассматривается распространение одномерной плоской волны, для которой свойственно распределение частиц по скоростям. Показано, что за смену режима структурообразования отвечают стати­стические динамические характеристики среды: дисперсия скорости частиц, средняя скорость частиц в ударной волне, скорость изменения
дисперсии и средней скорости частиц на
мезоуровне.
Уравнения баланса импульса и массы сводятся к уравне-
нию второго порядка с источником в правой части:
2
m
1
,
tx
(2.15)
где С
u
2
u
2
— продольная скорость звука;
0
pl
— массовая скорость, соответствующая пластической
2
2
t
Cx
0
2
pl
u
2
u
2
x
11

2
деформации;
λ, μ — коэффициенты Ламе. Левая часть уравнения — волновой оператор. В случае нестационарного пластического фронта u
pl
может быть представлена в виде суммы равновесной и неравновес­ной частей, каждая из которых определяется равновесной f неравновесной f
частями функции распределения мезочастиц
1
0
и
по скоростям v:
pl
p
,
1
1
pl
нo
fffdvvfudvvfu
.,
101
(2.16)
8
2.Физико‐математическаямодель
Равновесная часть функции распределения представляется максвелловской функцией распределения (по аналогии с ме­ханикой жидкости и газа):
23
nf
0
DDD
21

3
2
uv
ii
,exp
(2.17)
2
D
i
где n — средняя плотность мезочастиц;
u
— i-я компонента массовой скорости;
i
2
— дисперсия i-й компоненты массовой скорости.
D
i
Неравновесная часть функции распределения f
вычисля-
1
ется с использованием кинетического уравнения релаксацион­ного типа и в динамически деформируемом твердом теле име­ет вид
f
p
f
где τ
— время релаксации функции распределения частиц по
р
1
2
v
2
D
0
,
(2.18)
u
скоростям, соответствующее изменению состояния системы частиц, при котором они успевают обменяться импульсами в результате взаимодействия друг с другом [202].
С учетом функции (2.18) выражение для неравновесной части средней массовой скорости примет вид
2
D
pl
u
н
Аналогично, представляя напряжения σ
1
2
. (2.19)
u
pl
в виде суммы равновесной и неравновесной частей, выражая шаровую ком­поненту неравновесной части через дисперсию скорости мезо­частиц D
2
, получим выражение (2.15) в виде
22
D
tu
 
x
.
(2.20)
2
C
0
x
u
u
2
2
t
2
2
2
2
C
o
2
x
D
 
Путем преобразований уравнение (2.20) сводится к виду
88
2.4.Критерийинициирования«турбулентного»режимадинамическогоуплотнения
DD u u DD u DD
 
 
2

124
Cu
0
 
 

222.
uu
 
    
uu x t x uu xxuu
 
2
  
DD DD uDD
     
2
xuu tuu tuu
  
22 2
22 2
(2.21)
Когда
 

, (2.22)
1
DD DD
 
uu uu

уравнение (2.21) сводится к следующему выражению:
2
2
t
u
D
D
u
 
t
.
(2.23)
x
u
Уравнение (2.22) имеет только тривиальные решения. Та­ким образом, при выполнении условия (2.22) волновое движе­ние среды характеризуется резким затуханием массовой ско­рости в главном направлении — в направлении распростране­ния импульса давления. В соответствии с законом сохранения энергии с момента выполнения условия (2.22) диссипативные процессы динамического уплотнения порошкового материала реализуются развитием вихревых структур
и т. п.
Размер мезообъемов (подобластей конечного размера внут­ри рассматриваемого представительного объема реагирующего материала), на которых выполняется условие (2.22), как раз и обусловливает характерный локальный размер возможных яче­ек турбулентности. С этим размером связан размер субструктур продуктов физико-химических превращений, а значит, анализ возможных размеров ячеек турбулентности представляет от­дельный результат вычислительных экспериментов
.
Все вышесказанное позволяет заключить, что критерий (2.22) определяет одно из необходимых условий формирова-
ния субструктуры наноскопических размеров на мезоуровне реагирующего порошкового материала и саму возможность ударного синтеза в процессе динамического уплотнения.
89
2.Физико‐математическаямодель
x
2.5.Структурнаямакромеханикаразрушения.
Дискретныеособенностидинамикиразрушения
твердыхтел
Структурная механика разрушения, предложенная Г. Ней­бером и В. мер инженерного подхода к сложнейшей и высокоразработан­ной проблеме прикладной науки.
Критерий разрушения Нейбера — Новожилова может быть построен на основе следующих базовых принципов:
1. Все твердые тела состоят из пространственных струк-
турных элементов конечного размера.
2. Элементарный акт разрушения есть разрушение одного
структурного элемента.
3. Параметры критерия разрушения разом, чтобы в «канонических» случаях получились результа­ты классической теории разрушения.
Эти принципы применимы для формулировки условия
разрушения.
Рассмотрим, например, случай плоской деформации
образца в условиях статического нагружения. Допустим,
Рис. 4. Элемент структуры
на пространственной шкале
В. Новожиловым [294—297], — характерный при-
выбираются таким об-
что разрушение развивается
d
вдоль направления Ох, яв­ляющегося линией симметрии. Будем считать, что на про­странственной шкале задана структура
, а линейный размер
одной структурной составляю­щей равен d (рис. 4).
90
Соседние файлы в предмете [НЕСОРТИРОВАННОЕ]