CC BY
17
УДК548.4:536.48
О ГЕОМЕТРИЧЕСКОМ ПАРАМЕТРЕ В УРАВНЕНИИ ИНТЕНСИВНОСТИ ГЕНЕРАЦИИ ДИСЛОКАЦИИ
М.И. Слободской, JI.E. Попов (Томск) Abstract
From the results of computer simulation of planar dislocation loop emission by Frank-Read source of further evolution of the dislocation loop in a random field of point-like obstacles is estimated of the geometric parameters in the equation rate of generation a dislocations and they is compared with corresponding of analytically calculations
Введение
Один из возможных подходов к изучению пластической деформации кристаллических тел как целостного явления в достаточно общих условиях предложен в цикле работ Л.Е. Попова с соавторами (библиография, например, в [1-5]) Отталкиваясь от атомно-дислокационной концепции упрочнения и отдыха авторы сформулировали модели пластичности кристаллических тел, основным ядром которых является система дифференциальных уравнений эволюционного типа (уравнения баланса деформационных дефектов). Очень важно, что каждый параметр таких моделей имеет вполне определенный физический или прозрачный геометрический смысл и может быть вычислен, либо может быть указан интервал его возможных значений. Это убедительно проиллюстрировано в ряде публикаций (ссылки в [1-5]). Однако с геометрическим параметром в уравнении баланса дислокаций, которому посвящена настоящая работа, (необходимые пояснения - ниже) дело обстоит сложнее
Остановимся на самом элементарном случае: рассмотрим эволюцию только одной изолированной дислокационной петли и временно перейдем на чисто геометрический язык (длина дислокационной петли - периметр Р соответствующей замкнутой плоской кривой; площадь, заметенная дислокационной петлей, - площадь S, ограниченная кривой). Тогда плотность дислокаций р = Р, кристаллографический сдвиг a=bS, а изменение плотности дислокаций (периметра) за некоторое время At представимо в виде
Ap = AP*r!(L)AL, (1)
где L - некоторая линейная характеристика замкнутой плоской кривой, b -модуль вектора Бюргерса, Г, - геометрический параметр, определяемый формой кривой, которая, в свою очередь, может изменяться вместе с L. Г, (L) и L таковы, что Гj (L) • L = Р (для окружности, если за L принять ее диаметр, Г] ~п). Выбор L, вообще говоря, произволен и определяется соображениями удобства, но для сравнения результатов по различным классам линий вопрос о его выборе надо решить однозначно. Тем более, что иногда одним и тем же термином называют разные характеристики плоских кривых и множеств на плоскости (например, диаметр). В дальнейшем мы придерживаемся терминологии, устоявшейся в теории метрических пространств и берущей свое начало в интегральной геометрии [6] - [8]: А(ф) - ширина
множества л в направлении ф (плоской замкнутой кривой - его границы) - это длина ортогональной проекции Н на прямую, параллельную направлению ф; Б = тах{Д(ф)},
ф
К-тт{А(ф)}, соответственно, диаметр и широта множества Е. В случае выпуклых
ф
замкнутых кривых [6]
Р = |д(ф)с1ф (2)
()
и поэтому среднее значение ширины <Д(Ф)) = Р/ п. Очевидно, что смысл последнего равенства в том, что средняя ширина множества (А(ф))-Ь является диаметром окружности, равновеликой по периметру границе данного множества. Именно эта величина и выбрана в качестве линейной характеристики замкнутой кривой (дислокационной петли). Для невыпуклых множеств равенство (2) надо рассматривать как приближенное Далее,
Да = Ь ■ ДБ ~ Ь ■ Г2 (Ь) • ЬАЬ, (3)
г де Ь - модуль вектора Бюргерса, Г2 (I.) - параметр формы кривой, в общем случае отличный от Г| (1.). хотя подходящим выбором линейной характеристики Ь для определенного класса линий эти геометрические параметры можно сделать одинаковыми. (Для окружности, в случае выбора за I. ее диаметра, Г2 (Ь) = л: / 2 ; если же за I. принять радиус окружности, то Г[ ~Г2 = 2- л ) При ДЬ—» 0 комбинируя (1) и
(3), находим
с!р с1Р Г1(Ь) ёР П(Ь)
.----=------'-----::> Ф = Ь— (Ф - - - - ) (4)
с1а ЬсШ Г2(Ь)ЬЬ ёв’ 1 Г2(Ь)' v ’
Из (4) следует, что в общем случае параметр Ф может изменяться в довольно широких пределах: от отрицательных значений, так как периметр замкнутой плоской кривой может убывать с ростом площади, ограниченной этой кривой, - до бесконечности, если принимать во внимание, что дислокация при определенных условиях имеет дробную размерность (линия фрактального типа) [17] Равенство (4) справедливо для любых дифференцируемых (кусочно-дифференцируемых) плоских кривых. Если ограничиться только выпуклыми кривыми, то оно принимает вид
Р (1Р
ф = —-- (5) 7Т • с! о
Как (4), так и (5) практически не пригодны без каких-либо дополнительных предположений для вычисления параметра Ф в силу того, что в достаточно общих условиях неизвестна зависимость Р(8) или зависимости Р и 8 от какого-либо вспомогательного параметра в процессе эволюции кривой. Поэтому, для получения конкретных результатов нужно либо:
а) ввести дополнительные предположения;
б) использовать результаты имитационного моделирования процессов генерации и эволюции дислокационной петли,
в) возможно, что-то можно получить из принципа “длины и площади” [11, с. 354].
Приближение выпуклых самоподобных кривых
Предположим, что с изменением средней ширины замкнутой плоской кривой 1. она остается все время подобной самой себе - ее тип и параметры формы не изменяются (выполняется преобразование подобия, а параметр(ы) формы - инвариант этого преобразования). Тогда из (4) имеем
Ф
271-8
(6)
Таким образом, для вычисления параметра Ф имеем две формулы (4) - общие или, что тоже, (5) и (6)- справедливые только в сформулированных выше предположениях. Значение Ф, подсчитанное по (4) будем обозначать Фа;г, по (6) - Фш, Конечно же, при выполнении предположения самоподобия Фш,= ф<|,[ - Ф.
Таблица 1
Значения параметра Ф для некоторых классов кривых
N Классы линий Ф
1 Окружности 2
2 Треугольники Рельё [6] п-угольники (п=2к+1) я /(я --у/з ) * 2,2288 я/[я-П tg(я/2n)] => 2
3 Прав многоугольники (2п/я) tg(тc /п) => 2
4 Квадраты 8/я * 2,5465
5 Прямоугольники 2(1+А,)2/(яХ)
6 Ромбы 8/(тс вша)
7 Параллелограммы 2(1 +Х)2/(яА. віпа)
8 Эллипсы (8Е2(є))/(л2(1-е2У15)
9 Эпициклоиды [7] т - 1 - кардиоиды 32(1-1/(ш+2)) /я2 64/(3я2) «2.1615
10 Гипоциклоиды [7] т = 3 - Штейнера [7] т = 4 - астроиды [7] 32(1 + 1/ (т-2))/ я2 64/я2 «6,4846 48/я2* 4,8634
В табл. 1 представлены результаты вычислений параметра Ф по формулам (4),
(6). Здесь приняты следующие обозначения: п - число сторон многоугольника (для многоугольников Рельё - нечетное), А, - отношение большей стороны к меньшей для параллелограммов и прямоугольников, а - острый угол параллелограмма или ромба, £
- эксцентриситет эллипса, Е(х) - полный эллиптический интеграл второго рода [18], т -модуль эпи(гипо) циклоид - отношение радиусов неподвижного и катящегося кругов
[7] (в данном случае предполагается целым числом). Фигуры 1-8 выпуклые, из них 1-2 относятся к так называемым множествам постоянной ширины [6,8]. Для кривых 9-10 не выполняется предположение выпуклости, что в данном случае несущественно.
а б
Рис 1 Зависимости геометрического параметра Ф от инвариантных относительно преобразования подобия параметров некоторых плоских фигур: а) - нумерация
кривых соответствует их нумерации в табл. 1; б) - зависимости для семейства параллелограммов (квадраты, ромбы, прямоугольники). Разные кривые отвечают различным значениям а, нижняя линия с символами соответствует семейству эллипсов
с отношением полуосей X
Данные, представленные в табл..1 и на рис.1, в комментариях не нуждаются. Отметим лишь следующее. Для эхш(гипо)циклоид с модулем ш известные формулы дифференциального исчисления для периметров и площадей плоских кривых дают
Рт= 8Я( 1+1/т) , 8т= лЯ2(1±1/ т)(1±2/т).
Переходя к пределу при фиксированном Я (радиусе неподвижного круга),
2
имеем Ыт Рт -8Я , 1лт 8Пг тсК В этом предельном случае эпи(гипо)циклоиды
1П->сю III—>00
«стремятся» к окружности радиуса Я. Но длина этой «окружности» 8Н (вместо 2лЯ), а площадь, ограниченная ею, как и для обычной окружности - тЛ2. Видимо здесь имеет место промежуточный вариант между классическими и «береговыми» [17] линиями
Конечно же, табл.1 можно продолжить. Но даже, если ограничиться только алгебраическими кривыми, то, как известно, кривых второго порядка 4, третьего -1875, четвертого - 84551... [7]. Поэтому крайне необходимо иметь общие, не зависящие от класса кривых, оценки геометрического параметра Ф. Так как Ф может принимать отрицательные значения и стремиться к бесконечности, то общих двухсторонних неравенств для Ф (без каких-либо ограничений на кривые), по-видимому, предложить нельзя. Поэтому не случайно, что в классической монографии по неравенствам [9] и в методах, развиваемых в [10], нет результатов, непосредственно применимых к подобным оценкам. Но если ввести какие-либо ограничения на форму кривых, скажем, в виде определенной комбинации их широты К и диаметра Б , то из классических неравенств интегральной геометрии [6,8]
Р2 >4л8 , л/38>К2 ,Р>тгК,28>КО,
Р>2(02-К2)112 + 2Кагсяп(К/Б),
Р<2ф2 -К2)1/2 + 20агсяп(К/В),
28<К(Ог -К2)172 +В2 агсяп(К/0) , 48<2КР-тгК2
можно получить достаточно точные двухсторонние оценки параметра Ф. Ограничения на форму кривых, конечно же, должны соответствовать конкретному содержанию рассматриваемой физической задачи (анизотропия, специфика барьеров, различия в краевых и винтовых дислокациях...). К сожалению, полная система неравенств остается пока (1983г.) неизвестной. “Нахождение для каждого из неравенств множеств, для которых выполняется равенство, - очень трудная математическая задача’'[6].
Оценки геометрического параметра Ф по результатам имитационного моделирования эволюции дислокационной петли в поле случайно расположенных точечных препятствий
При имитации на ЭВМ генерации и эволюции дислокационной петли [12, 14 -16, 19] в плоскости кристаллографического скольжения, содержащей точечные [13] случайно распределенные препятствия, в определенные промежутки времени, кроме всего прочего, регистрируются периметр дислокационной петли и её площадь (площадь, по которой прошел кристаллографический сдвиг). Эти данные позволяют проследить за геометрическим параметром Ф для эволюционирующей по своим законам (а не в предположении самоподобия) дислокационной петли. При этом имеется возможность вычислять как Ф^, так и Фш1 и тем самым судить о том, как сильно влияет невыполнение предположения самоподобия на приближенный подсчет параметра Ф.
В работе [14] рассматривалось термоактивируемое расширение дислокационной петли в плоскости кристаллографического скольжения со случайно распределенными однотипными препятствиями малой прочности под действием приложенного напряжения. Отмечалось, что в стабильных конфигурациях дислокационные петли по форме близки к окружностям, то есть можно считать, что предположение самоподобия выполняется. Вычисленные значения геометрического параметра Ф (как по формуле
(4) так и по (6)) изменяются незначительно: от 2,0 до 2,81, 2,0 - 2,85, 2,0 - 2,64 в разных ЭВМ-экспериментах. При этом средняя ширина дислокационных петель изменялась в первом эксперименте от 28,4 до 101,8 , во втором - 71,7 - 121,2 и от 101,4 до 146,2 средних расстояний между ближайшими стопорами - в третьем. Зависимости Ф от Ь не являются монотонными.
Значение Ф , вычисленное по результатам работы [15] по формуле (6) для конфигурации, когда дислокационная петля отшнуровывается от источника Франка -Рида, составило 2,21, что сопоставимо со значением Ф для кардиоиды (см. табл. 1).
Иная ситуация в случае эволюции дислокационной петли от источника Франка -Рида, когда в плоскости залегания источника присутствует спектр препятствий с существенно различными прочностями и среди них есть атермические препятствия [16]. Ряд последовательных конфигураций дислокационной петли в этом случае представлен на рис.2.
Рис. 2 . Эволюция дислокационной петли от источника длиной 6 /л под действием внешнего напряжения т* Номера в углах рисунков соответствуют номерам дислокационных конфигурации из табл.2. В квадрате сетки каждого рисунка в среднем по 225 стопоров
Препятствия в плоскости кристаллографического скольжения, за исключением первого блока, в котором представлена эволюция дислокационного сегмента до конфигурации потери устойчивости - конфигурации, после которой все процессы происходят в над барьерном режиме (сегмент теряет устойчивость), не изображены, их
в каждом из представленных фрагментов - более 200 тысяч.. Замкнутые маленькие кривые (часть из них на рисунке выглядит жирными точками) - «острова» незавершенного кристаллографического сдвига [19], оставшиеся за фронтом скользящей дислокации и представляющие собой труднопреодолимые участки плоскости кристаллографического скольжения, которые дислокация обходит по механизму Орована. Конфигурация 11 - дислокационная конфигурация в момент замыкания сегмента-источника в дислокационную петлю (момент срабатывания источника). Вставка в нижней части рисунка («хвост» от источника) - конфигурация, оставшаяся после отшнуровьшания дислокационной петли. Его «геометрия» существенно влияет на генерацию источником последующих петель, вплоть до запирания источника. Количественная информация о дислокационных конфигурациях, частично представленных на рис 2, сведена в табл. 2
Таблица 2
Некоторые характеристики дислокационных петель, представленных на рис. 2
№ Число узлов Периметры Площади Пройдено стопоров
РЬ* РО* ЭЬ* 80* всего перереза нием по мех. Орована
1 3 4 5 6 7 9 11 17
2 70 90,9 4,19 365,3! 0,32 315 301 14
3 586 652,8 158,53 7570,11 140,49 6893 6452 441
4 906 992,6 428,03 16070,7; 365,19 14963 13902 1061
5 1412 1566,8 699,79 25153,5! 570,00 23609 21924 1685
6 1487 1680,5 791,48 27387,6; 649,39 25732 23869 1863
10 1520 1715,1 1126,02 33326,8! 1061,46 32324 28885 3439
11* 1473 1655,3 1209,47 33772,6І 1262,84 32630 29110 3520
12 1563 1768,9 1908,32 44685,4! 2905,74 41782 37238 4544
20 2045 2279,4 2048,84 52024,2! 2961,56 50205 43681 6524
21 2542 2830,0 2208,18 59337,4| 3028,89 57201 50263 6938
В табл.2 в колонке «число узлов» приведено число стопоров, с которыми в данный момент времени контактирует дислокационная петля, где РЬ* и РО* -периметры петли и всех островов незавершенного кристаллографического сдвига в данной конфигурации в единицах среднего расстояния между ближайшими стопорами в плоскости кристаллографического скольжения 1„ ; 8Ь* и 80* - площади, охватываемые дислокационной петлей и островами незавершенного кристаллографического сдвига в единицах (1„)2
Уже простой взгляд на рис.2 убеждает, что предположение самоподобия эволюционирующей дислокационной петли существенно нарушено и формулу (6) для вычисления геометрического параметра использовать, вообще говоря, нельзя. Это же следует и из данных, представленных на рис.З, на котором, в целях сравнения результатов, все же приведены значения Фц^.
Рис.З. Зависимости геометрических характеристик дислокационных конфигураций от средней ширины дислокационной петли (пояснения в тексте)
Зависимости (а) непосредственно регистрируются в машинных экспериментах, (б) - результаты численного дифференцирования кривых (а). В первом и во втором (разница в длинах источников - 4 18 и 12 15 соответственно) блоках - результаты по «до критическим» конфигурациям (блок 1 на рис.2). Зависимости, представленные в третьем и четвертом блоках рисунка, построены по конфигурациям дислокационных петель от соответствующих источников. Их в каждой ячейке по две - «короткие» кривые (незаполненные символы) и «длинные» (темные кружочки). Первые из них построены без учета островов незавершенного кристаллографического сдвига, при построении вторых в периметр дислокационной петли включены их периметры, а из площадей, ограниченных дислокационной петлей, вычтены площади островов (в этих
областях не прошел кристаллографический сдвиг). Как и следовало ожидать, значительна разница между и Фш[ Встречаются конфигурации с отрицательными значениями геометрического параметра; они как раз приходятся на те моменты, когда в результате эволюции дислокационной петли у нее сокращается периметр при некотором возрастании площади (петля «округляется»).
Заключение
Конечно же, усреднение по достаточно представительной группе дислокационных петель значительно сгладит различие теоретических и ЭВМ-экспериментальных оценок рассматриваемого в статье геометрического параметра Кроме того, при количественном описании интенсивности генерации дислокаций в процессе пластической деформации необходимо учитывать и длительность деформирующего воздействия; в простейшем варианте хотя бы два предельных случая: 1) когда время деформирующего воздействия меньше времени элементарного акта кристаллографического сдвига - формирования одной дислокационной петли (импульсное нагружение), 2) когда время деформирующего воздействия превышает время образования дислокационной петли и даже время формирования зоны сдвига, то есть время формирования всей серии дислокационных петель, испускаемых источником до его остановки (статический случай). После этого реально подняться на следующий структурный уровень в описании пластической деформации - зону сдвига. Но э го предмет последующих исследований и публикаций
Библиографический список
1. Попов Л.Е., Конева Н А., Терешко И В. Деформационное упрочнение упорядоченных сплавов. - М.: Металлургия, 1972 - 256 с.
2. Попов Л.Е., Кобытев В С , Ганзя Л.В Теория деформационного упрочнения сплавов. - Томск: Изд-во Том. ун-та, 1981. - 176 с.
3. Попов Л.Е., Кобытев В С , Ковалевская Т А. Пластическая деформация сплавов -М.: Металлургия, 1984,- 182 с.
4. Попов Л.Е., Пудан Л.Я., Колупаева С.Н. и др. Математическое моделирование пластической деформации. - Томск: Изд-во Том. ун-та, 1990 - 185 с.
5. Колупаева С.Н., Старенченко В.А , Попов Л.Е Неустойчивости пластической деформации кристаллов - Томск: Изд-во Том ун-та, 1994 - 300 с
6. Сантало Л. Интегральная геометрия и геометрические вероятности. - М.: Наука, 1983 - 359 с.
7. Савелов А.А. Плоские кривые. Систематика, свойства, применения - М.: Физматгиз, 1960,- 294 с
8 Амбарцумян Р.В., Мекке И., Штоян Д. Введение в стохастическую геометрию. - М : Наука, 1989.-401 с
9. Харди Г Г., Литтльвуд Д.Е., Полна Г. Неравенства. - М.: ИЛ, 1948 - 456 с.
10.Маршалл А., Олкин И. Неравенства: теория можаризации и ее применения. - М : Мир, 1983 - 575 с.
11 Математическая энциклопедия: Гл. редактор И.М.Виноградов. Т.2 Д - Коо - М.: Советская энциклопедия, 1979 - 1104 с. (ст. «длины и площади принцип», с 364)
12.Слободской М И., Матющенко А.В., Голосова Т.Н. Алгоритмизация имитации образования дислокационной петли источником и процесса ее эволюции в плоскости кристаллографического скольжения со случайно распределенными дискретными препятствиями // Математическое моделирование систем и процессов
- 1995 №3. - С. 88-96.
13.Зайцев С И Моделирование движения дислокаций через точечные препятствия // Дефекты в кристаллах и их моделирование на ЭВМ - Л. Наука, 1980. - С 178-191
14 Слободской М И., Кобытев В.С., Попов Л.Е Зависимость средней площади, заметаемой дислокационной петлей после одной термической активации// Изв. вузов. Физика. - 1985. №3.- С. 119-120
15 Слободской М.И., Голосова Т.Н., Попов Л.Е. Источник дислокаций в поле дискретных стопоров// Изв. вузов. Физика. - 1990. №12,- С. 20-24.
16Голосова Т.Н., Слободской МИ., Попов Л.Е. Моделирование источника дислокаций в поле активируемых и неактивируемых дискретных препятствий // Изв. вузов. Физика. - 1992. №10. - С. 20-24.
ПОлемской А.И, Флат А.Я. Использование концепции фрактала в физике конденсированной среды // УФН. - 1993. - Т 163 №12,- С. 1-50
18 Справочник по специальным функциям. Под редакцией М.Абрамовича и И.Стиган -М : Наука, 1979,- 832 с.
19 Слободской МП, Матющенко А.В. Имитационное моделирование генерации дислокационной петли в поле случайно расположенных дискретных препятствий // Математическое моделирование систем и процессов - 1996. №4. - С. 88-95.
