Измерение магнитострикции ферромагнетика с помощью тензодатчика

Удмуртский Государственный Университет
Кафедра Физики Твердого Тела
Лабораторная работа № 4
Измерение магнитострикции ферромагнетика с помощью тензодатчика
г. Ижевск
Бреган Андрей гр.18-31 1998 год.
Цель работы: определение продольной магнитострикции никеля в зависимости от амплитуды напряженности магнитного поля.
Теория.
§ 1 Введение
Данная работа посвящена изучению поведедения ферромагнетиков в магнитном поле.
Хотя магнитное взаимодействие является малой поправкой к электрическим обменным силам, обусловливающим самопроизвольную намагниченность, тем не менее, они играют решающую роль во всем сложном комплексе явлений технического намагничивания. Поэтому выяснение физической природы магнитного взаимодействия в ферромагнетиках имеет не только теоретическое значение, но необходимо и для ясного понимания механизма тех физических процессов, которые обусловливают всю практическую ценность явления ферромагнетизма.
Напомним, что ферромагнетиками называются вещества, в которых магнитные моменты ориентированы вдоль выделенного направления.
Монокристаллы ферромагнетиков анизотропны в магнитном отношении. В качестве примера магнитокристаллической анизотропии на рис.1 приведены кривые намагничивания I(Н) монокристалла кобальта, снятые вдоль гексагональной оси (ось с) и перпендикулярно к ней (ось а). Как видно из рисунка, если магнитное поле H || c, то достаточно приложить поле в несколько сот эрстед для того, чтобы намагнитить кристалл до насыщения. При Н ( с насыщение достигается только при Н ( 104 Э.
Наиболее резко магнитная анизотропия, проявляется в кристаллах гексагональной симметрии (Со, Tb, Dy). Из анализа кривых I(Н), снятых по различным кристаллографическим направлениям, следует, что в ферромагнитных монокристаллах существуют направления, называемые «осями легкого намагничивания» (ОЛН), и направления, называемые «осями трудного намагничивания» (ОТН). Известно, что минимум свободной энергии магнитокристаллической анизотропии достигается, когда намагниченность ориентирована вдоль ОЛН. Для поворота Is из этих направлений требуется затрата определенной работы, которая приводит к росту энергии магнитной или магнитокристаллической анизотропии. Энергией магнитокристаллической анизотропии называют ту часть энергии кристалла, которая зависит от ориентации вектора намагниченности относительно кристаллографических осей. [pic] Рис.1. Кривые намагничивания I(Н) монокристалла кобальта, снятые вдоль гексагональной оси (ось с) и перпендикулярно к ней (ось а).
В случае кобальта эта энергия минимальна, если намагниченность направлена вдоль оси с (при комнатной температуре). При вращении намагниченности Is от оси с энергия анизотропии увеличивается с увеличением угла ( между осью c и направлением Is, достигает максимума при (=90°, т. е. при Is ( с, и затем уменьшается до первоначального значения при ( =180°.
§ 2. Спонтанная магнитострикция и ее вклад в магнитную анизотропию
При возможных изменениях ориентации самопроизвольной намагниченности в кристалле изменяются равновесные расстояния между узлами
решетки. Поэтому возникают самопроизвольные магнитострикционные
деформации. т.е. Опр. При перемагничивании ферромагнетика имеет место магнитное взаимодействие элекектронов, которое влияет на межатомное расстояние, вызывая деформацию кристаллической решетки, что сопровождается изменением линейных размеров тела и появлением соответствующей магнитоупругой энергии. Это явление называется магнитострикцией. В частном случае кубического кристалла в отсутствие внешних напряжений свободная энергия магнитного и упругого взаимодействия (с точностью до шестых степеней в направляющих косинусах вектора Is и вторых степеней тензора магнитострикционных напряжений), равна сумме энергии магнитокристаллической анизотропии fa, упругой энергии fупр и магнитоупругой энергии fму:
fa((i ,ei j)= fa((i ,ei j)+ fупр.((i ,ei j)+ fму. ((i ,ei j) (1)
1) Можно феноменологическим путем получить выражение плотности fa энергии магнитной анизотропии, раскладывая эту энергию в ряд по степеням направляющих косинусов вектора намагниченности (i относительно осей симметрии кристалла. Сначала найдем выражение fa для кобальта, имеющего гексагональную решетку с ОЛН - с, для которого (i =( = cos (Is,с) = cos (. Для гексагональной решетки, обладающей центром симметрии, операция замены ( на - ( должна оставлять энергию инвариантной относительно такого преобразования симметрии. Следовательно, в разложении останутся только члены с четными степенями а, т. е.
fa=K1((2 + K2((4 + ......
(2)
где K1((2 и K2((4 и т. д. - параметры магнитной анизотропии; fa чаще записывают в следующем виде:
fa = K1 sin2(+ K2 sin4(+..., (3)
где K1 и K2 называют 1-й и 2-й константами магнитной анизотропии. Энергия анизотропии кристаллов гексагональной системы в общем случае должна зависеть от азимута (. Но эта зависимость является очень слабой, и ею обычно пренебрегают. Для кубических кристаллов, таких как Fe, Ni, энергия анизотропии выражается в функции направляющих косинусов ((1, (2, (3) намагниченности Is относительно трех ребер куба:
((1=cos(Is, [100]); (2=cos(Is, [010]); (3=соs(Is, [001]). (4)
Энергия анизотропии должна быть такой функцией (1 , (2 , (3, которая оставалась бы инвариантной при преобразованиях симметрии кубического кристалла.
В кубическом кристалле плоскости типа [100] являются плоскостями симметрии. Зеркальное отражение вектора Is в такой плоскости должно оставлять функцию fa((1, (2, (3) инвариантной. Отражение, например, в плоскости (100) заменяет (1 на - (1,оставляя (2 и (3 неизменными. Аналогично зеркальное отражение в плоскостях (010) и (001) изменяет знаки соответственно у (2 и (3. Следовательно, функция fa((1, (2,(3) должна быть инвариантной относительно преобразований
(i ( - (i ( i = 1,2,3)
(5)
Кубический кристалл имеет также плоскости симметрии типа {110}. Отражение в этих плоскостях соответствует преобразованиям
(i ( - (j ( i( j = 1,2,3)
(6)
Первым членом разложения энергии анизотропии кубического кристалла по степеням (1 , (2 , (3, удовлетворяющим требованиям симметрии (5,6), является (21 + (22 + (23 , но этот член разложения всегда равен единице и, следовательно, не описывает эффекта анизотропии.
Следующий член (четвертого порядка относительно (i), (41 + (42 + (43 может быть приведен к виду
(41 + (42 + (43 = 1- 2((21(22+(22(23+(21(23) (7)
так как ((21 + (22 + (23)2 = 1. Далее, член шестого порядка приводится к виду
(61 + (62 + (63 = 1- 3((21(22+(22(23+(21(23)+3(21(22(23 (8)
так как ((21 + (22 + (23)3 = 1.
Энергия анизотропии на единицу объема кубического кристалла с точностью до членов шестого порядка относительно (i представляется в виде линейной комбинации
fa=K1((21(22+(22(23+(21(23)+K2(21(22(23 (9)
Часто членом K2(21(22(23, который обычно меньше первого члена в (9), пренебрегают. Тогда:
fa=K1((21(22+(22(23+(21(23)
(10)
Знаки констант анизотропии K1 и K2 и их относительная величина определяют то кристаллографическое направление, которое в данном кристалле будет «легким». Если К1>0, то первый член в (9) минимален при направлении намагниченности вдоль осей [100], [010], [001], которые в этом
случае являются осями легкого намагничивания.
Если К1j) ;
4) S ((2i-1/3)((2j-1/3) = S (2i(2j – 1/3 (i , j=1,2,3, i>j) ;
f(0)му.= – [B21/(C11 – C12)][ 2/3 – 2 S (2i(2j] - B22/C44 S (2i(2j ,
(i , j=1,2,3, i>j) (22)
f(0)упр.= Ѕ C11 [B21/(C11 – C12)2] [ 2/3 – 2 S (2i(2j] + Ѕ C44 B22S (2i(2j+
+C12[B21/(C11 – C12)2] S (2i(2j – 1/3 ,
(i , j=1,2,3,i>j) (23)
Подставляя (21) и (23) в (1) и учитывая (10), (18) и (19), получим следующее выражение для плотности анизотропной части магнитной энергии кристалла при отсутствии упругих внешних напряжений:
f =(K1+?K1)S (2i(2j (i , j=1,2,3, i>j) (24)
где добавка ?K1 к первой константе анизотропии, обусловленная спонтанной магнитострикционной деформацией равна
?K1= [2B21/(C11 – C12)] + [B22/C44] – [C11{B21/(C11 – C12)2}]+
+ [Ѕ C44 (B22/C244)]+[ C12 {B21/(C11 – C12)2] =
= [B21/(C11 – C12)] – Ѕ[B22/C44]=9/4 (2100(C11-C12) – 9/2(2111C44
(25)
Как видно из (24), вид зависимости плотности энергии от направляющих косинусов не изменился, но константа анизотропии благодаря спонтанной деформации решетки увеличилась.
§2. Физическая природа естественной магнитной анизотропии.
В первых работах Акулова магнитное взаимодействие в ферромагнитных кристаллах с микроскопической точки зрения трактовалось чисто классическим путем. Квантовомеханическая трактовка была дана в работах Блоха и Джентиля. Классическую теорию температурной зависимости констант магнитной анизотропии развили Акулов и Зинер, исходя из представления о том, что около каждого узла решетки можно выделить области ближнего магнитного
порядка с не зависящими от температуры локальными константами анизотропии. Локальные мгновенные намагниченности этих областей из-за теплового движения распределены хаотически и образуют среднюю намагниченность всего кристалла. Отсюда удается определить связь между температурным ходом констант анизотропии и намагниченности в виде
Kn(T)/Kn(0) = [Is (T)/Is (0)]n(2n+1) ,
(26)
где n – порядок константы. Таким образом, мы приходим к универсальной зависимости K1? I3s и K2? I10s. Pезультат (26) получается в приближении теории молекулярного поля . Микроскопические трактовки этой проблемы даны в работах Ван - флека и Канамори. В основе всех расчетов по микроскопической теории магнитной анизотропии лежит учет магнитного взаимодействия между спиновыми
и орбитальными магнитными моментами электронов, принимающих участие в ферромагнетизме. В общем случае оператор магнитной энергии складывается из трех членов.
Hмагн.=U1+U2+U3
(27)
где U1 — оператор, соответствующий движению электронов относительно ионов решетки,— спин-орбитальная энергия; U2 — оператор магнитной энергии, возникающей вследствие относительного движения самих электронов, —орбитальная энергия; U3 — оператор энергии магнитного взаимодействия спиновых магнитных моментов электронов — спиновая энергия (в первом приближении имеет вид дипольного взаимодействия). Эффект орбитального взаимодействия U1 и U2 проявляющийся в случае изолированных атомов в образовании тонкой структуры спектральных линий приводит к появлению «внутренних магнитных полей» порядка 105 э. С другой стороны, «эквивалентное магнитное поле» анизотропии ферромагнетиков, определяемое величиной поля, при котором достигается насыщение в монокристалле вдоль труднейших направлений намагничивания, оказывается порядка 102 э и лишь в редких случаях (Со, пирротин) достигает 103—104 э. Объяснение этого несоответствия заключается в том, что в отличие от атомов, где орбитальные моменты отличны от нуля (за исключением s-состоянии), в ферромагнитных кристаллах (например, в d-металлах и сплавах), как показывают измерения гиромагнитного эффекта, средний орбитальный магнитный момент по кристаллу почти всегда практически равен нулю. Поэтому в первом приближении эффект спин-орбитальных энергий U1 и U2 также равен нулю. Отличный от нуля эффект получается лишь во втором и более высоких приближениях. Что же касается спиновой части магнитного взаимодействия U3, которая хотя и дает отличный от нуля эффект в первом приближении, но тем не менее не обеспечивает наблюдаемый на опыте порядок величины эффективных «полей» благодаря своей малости . Несмотря на отсутствие законченной квантовой трактовки магнитного
взаимодействия в ферромагнетиках, в этой области имеются известные
успехи. Так, например, удалось объяснить правильный порядок величины констант магнитной анизотропии. В частности, без всяких дополнительных соображений из теории следует, что в кубических кристаллах (Fe, Ni) константы анизотропии должны быть меньше по абсолютной величине, чем в случае гексагональных кристаллов (Со, пирротин). Это вытекает из свойств симметрии кубических кристаллов, в которых первое приближение для дипольной энергии U3 и второе приближение для орбитальных энергий U1 и U2 не приводит к зависимости свободной энергии кристалла от ориентации его намагниченности относительно кристаллографических осей. Для получения этой зависимости надо рассматривать следующие приближения, в то время как в гексагональных решетках анизотропия получается и в первом приближении для U3, и во втором для U1 и U2.
Остановимся несколько подробнее на микромеханизме явления естественной кристаллографической магнитной анизотропии. Поскольку
в создании самопроизвольной намагниченности ферро- и антиферромагнетиков основную роль играют электронные спины, то микроскопическая энергия, ответственная за магнитную анизотропию, должна зависеть от состояния этих спинов в кристалле, а также отражать симметрию распределения спиновой и зарядовой (орбитальной) плотности в кристалле. Наиболее простым является механизм магнитного дипольного взаимодействия спинов.
К сожалению, однако, учет лишь дипольного межэлектронного взаимодействия не может, как правило, объяснить наблюдаемую на опыте величину энергии магнитной анизотропии. Другой из упомянутых выше механизмов заключается в связи между
спином и орбитальным движением электронов [например, описываемой
членами U1 и U2 гамильтониана (27)]. Киттель дает следующее наглядное объяснение физического механизма магнитной анизотропии из-за спин-орбитальной связи.
В основу своего объяснения они кладут общепризнанное положение, что, само появление этой анизотропии обусловлено совместным действием спин- орбитальной связи, частичного замораживания орбитальных моментов неоднородными кристаллическими полями и орбитальным обменным взаимодействием соседних атомов. Таким образом, самопроизвольная намагниченность кристалла «чувствует» ионную решетку через орбитальное движение магнитных электронов. Спины, участвующие в намагниченности, взаимодействуют с орбитальным движением с помощью спин-орбитальной связи, а орбитальное движение связано с решеткой
полем лигандов. Микроскопическая энергия, возникающая благодаря этому механизму, может быть в свою очередь двух типов: 1) спин-орбитальная связь, которая зависит от спиновых состояний
двух или более ионов-носителей магнитного момента (парная модель
магнитной анизотропии); 2) связь, зависящая от спинового состояния только отдельных ионов
(одноионная модель магнитной анизотропии). Последний механизм оказывается наиболее близким к реальной ситуации, которая имеет место в неметаллических антиферро - и ферримагнетиках, в которых магнитноактивные ионы находятся в окружении магнитно-нейтральных анионов. Под действием поля лигандов, симметрия которого определяется типом кристалла, происходит расщепление уровней магнитного иона. В результате основному состоянию в зависимости от структуры кристаллической решетки будут соответствовать различные типы уровней, что приводит к магнитной анизотропии кристалла с магнитным порядком.
§ 3. Магнитострикция при техническом намагничивании
Известно, что в процессе технического намагничивания происходит смещение границ доменов и вращение вектора Is. Рассмотрим, как эти процессы влияют на изменение длины кристалла с положительной константой магнитной анизотропии K1. Пусть внешнее магнитное поле параллельно оси [110] и в исходном состоянии объемы Vi доменов, намагниченных вдоль шести направлений легкого намагничивания, равновелики: V0100= V0I00 =V00I0 =V0001 =V00I 0 =V000I =1/6 V, где V – объем кристалла. а) Смещение 180° доменных границ. При этом домены, намагниченные вдоль направлений [100] и_{010], поглощаются доменами,намагниченными в направлениях [100] и [010]. Изменения длины при смещении 180° доменных границ не происходит. После того как смещение этих границ заканчивается, объемы
доменов равны
V0I00= V00I0 =0 При этом средняя намагниченность кристалла
I= 2/3 (2) - Ѕ Is
б) Смещение 90° границ. При этом домены, намагниченные
вдоль направлений [100] и [010], поглощают домены, намагниченные в направлениях [001] и [001]. В конце этого процесса объемы доменов равны V100=V010=V/2,V001=V00I=0 и средняя намагниченность кристалла
I = (2) - Ѕ Is
Отноносительное изменение длины, вызванное смещением 90° границ, будет при этом равно
?l/l=1/3[((l/l)[001] – ((l/l)[100]]=1/3[((l/l)[001] – ((l/l)[010]] = (1/3)(3/4) (100 = ј (100
(28) гле
((l/l)[001] ((1 = (2 = 0, (3 = 1 ),
((l/l)[100] ((2 = (3 = 0, (1 = 1 ),
((l/l)[010] ((1 = (3 = 0, (2 = 1 ),
Отметим, что формула (28) определяет лишь конечное изменение длины, соответствующее намагниченности I= (2) - Ѕ Is, при которой заканчивается процесс смещения 90° границ при условии, что вклад в намагниченность от процессов вращения еще пренебрежимо мал.
в) Вращение. Если процессы смещения 180 и 90 границ заканчиваются в слабых магнитных полях, при которых вклад в намагниченность от процессов вращения пренебрежимо мал, то можно считать, что при дальнейшем росте магнитного поля вращение векторов намагниченности доменов к оси [110] происходит в плоскости (001) .Пусть ( - угол между направлением намагниченности доменов и полем. Тогда намагниченность вдоль поля Н равна I=Is cos (. Относительное удлинение вдоль оси [110] будет
(( l/l)[110] = - 1/2(100 + 3/4(100 ((21+(22) + 3/2(111(1(2 ,
где (1 = cos((/4-(), (2 = sin((/4-() , откуда
(( l/l)[110] = ј (100 + ѕ (111 (2cos2( -1) ( ј (100 + ѕ (111 (2(I/Is)2 -1).
Метод измерений и описание установки
На исследуемые в работе ферромагнитные образцы, наклеены проволочные тензодатчики, изготовленные из тонкой константановой проволоки. Константан имеет очень малый температурный коэффициент омического сопротивления и пригоден для выполнения прецизионных измерений. В результате явление магнитострикции при намагничивании образца изменяются и линейные размеры тензодатчика, что приводит к изменению его омического сопротивления. По величине изменения омического сопротивления можно судить о величине магнитострикции. Коэффициент пропорциональности изменения омического сопротивления тензодатчика называют коэффициентом тензочуствительности.
На рис.2 приведена электрическая схема установки для измерения магнитострикции. Схема собрана в виде стенда с клеммами для подключения внешних приборов, а именно: источника постоянного напряжения 10В. и микровольтметра постоянного тока. Через клеммы также подключается тензодатчик исследуемого образца и источник постоянного магнитного поля (соленоид). Измерительный проволочный тензодатчик ИТД является одним из плеч измерительного моста, во второе плечо которого включен компенсационный тензодатчик КТД, наклеенный на стальную полоску, находящуюся вне магнитного поля. Два других плеча измерительного моста образованы прецизионными резисторами R1 и R2 . Потенциометр R3 служит для грубой компенсации моста, а реоход R4 для точной компенсации. К диагонали измерительного моста через подвижный контакт реохода и клеммы x3,x4 “мкВ” Подключается внешний микровольтметр. Резистор R5 задает ток питания датчиков. На стенде размещен также переменный резистор R6 для регулировки тока в соленоиде и, соответственно, амплитуды магнитного поля в нем. Конструктивно-технологические особенности выполнения измерений магнитострикции.
Измерение магнитострикции является прецизионным процессом, в силу малости измеряемой величены. Даже у никеля, имеющего одну из самых больших величин магнитострикции, она составляет всего 30x10-6. Поэтому при питании тензодатчиков током порядка 5мА, разбаланс моста составляет в зависимости от величены ( от 1 до 10 мкВ, что соответствует изменению сопротивления тензодатчика 10-4-10-3 Ом. Для сравнения: величина термоЭДС контакта медь- конктантан – 30 мкВ/град С, а переходные сопротивления контактов могут достигать величин десятых долей Ома, т.е. ,по крайней мере, два мешающих фактора существенно превосходят измеряемую величину. Основные способы отстройки от мешающих факторов использованные в стенде:
-размещение контактов индуцирующих термоЭДС в соседних плечах моста для компенсации ее;
-использование скользящего контакта реохода (контакта с большим переходным сопротивлением) в высокоомной потнциалометрической цепи;
-обеспечение стабильного температурного режима за счет изготовления соленоида с большим запасом по мощности (во избежание его нагрева) и удаление его от измерительной схемы.
Перечень приборов и оснастки.
Для выполнения работы необходимы: - стенд контроля, - амперметр постоянного тока, - микровольтметр постоянного тока или гальванометр, - источник постоянного напряжения.
Порядок выполнения работы.
В качестве образцов работе используется ленты из никеля (образец №1) и лента из электротехнической стали марки 3414 толщиной 0,35 мм, вырезанной вдоль направления проката, совпадающего с направлением магнитной и кристаллографической текстуры (образец №2) и поперек направления проката, т.е. поперек направления магнитной текстуры (образец №3). Порядок выполнения измерений следующий:
1. Подготовить к измерениям рабочее место.
1. Подключить к клеммам “мкВ” микровольтметр постоянного напряжения к клеммам “U пит” источник постоянного напряжения, установив на нем напряжение 10 В. К клеммам “соленоид” подключить последовательно соедененные амперметр постоянного тока и соленоид. К клеммам “ИТД” подключить выводы тензодатчика, наклеенного на образце №1.
2. Установить переключатель пределов измерения микровольтметра в положение максимальной величине напряжения. Переключатель пределов измерения амперметра установить в положение “2А”, а ручку резистора регулировки тока соленоида в крайнее левое положение, соответствуещее минимальному току.
1. Выполнить измерения магнитострикции.
1. Включить в сеть микровольтметр и источник питания и дать ему прогреться в течении 5-10 мин.
2. С помощью резисторов “компенсация грубо” и реохода “компенсация точно” скомпенсировать измерительный мост – добиться минимальных показателей микровольтметра. Постепенно увеличивая чувствительность микровольтметра, довести ее до предела “100 кВ”.
3. Изменяя ток в соленоиде от минимального значения до 2А записать показания вольтметра, соответствуещее 5-6 точкам значений тока в соленоиде.
4. Выполнить измерения в том же порядке для образцов 2 и 3.
Примечания: 1. Перед заменой образца тензодатчика выключить источник напряжения, а ручки управления приборов установить в положение соответстстаующее п.1.2.
2. В случае, если не удается установить стрелку микровольтметра на нулевое положение в процессе компенсации моста, допускается производить отсчет показаний относительно полученного (ненулевого) значения шкалы.
3. Замену образов производить при выключенном источнике питания микровольтметра.
3. Обработка результатов измерений.
1. Рассчитать величину магнитострикции, соответствующую показаниям микровольтметра.
2. По номиналам резисторов, указанным в схеме на рис.2 рассчитать ток тензодатчика.
3. Зная разбаланс моста в мкВ и ток через тензодатчик рассчитать разбаланс моста в Омах и определить соответствующую этому разбалансу магнитострикцию, используя соотношение:
(=?l/l=((?R/R),
где ?R – разбаланс моста в Омах.
R = 200 Oм – сопротивление тензодатчика,
( = 2 - коэффициент тензочувствительности.
4. Построить зависимость ( =f(H) для всех образцов (постоянная соленоида К = 100 Э/А = H/J, где Н – амплитуда магнитного поля,
J – ток соленоида).
5. Объяснить различия в зависимостях ( =f(H) для образцов из электротехнической стали, вырезанных вдоль и поперек проката.
[pic]
Рис. 2. Электрическая схема установки для измерения магнитострикции.
Литература.
1. Вонсовский С.В. Магнетизм: «Наука», 1971, 2. Ивановский В.И. Физика магнитных явлений, Издательство Московского университета, 1981 3. Г.Рейнбот. Магнитные материалы и их применение. Перевод с немецкого Л.
Энергия, 1974, с. 375-380.
Контрольные вопросы.
1. Что такое магнитострикция ? Какова ее величина ? Чем явление магнитострикции отличается от пъезоэффекта ? 2. Вследствие чего возникают магнитострикционные деформации ? 3. Какова природа магнитной анизотропии ? 4. Какие оси в кубическом кристалле будут ОЛН, а какие ОТН ? Почему ? 5. Какие процессы влияют на изменение длины кристалла с положительной константой магнитной анизотропии K1 ? 6. Зная плотность свободной энергии можно сосчитать всю термодинамику системы. Чему равна плотность свободной энергии деформированного кристалла, намагниченного до насыщения? 7. Как зависят константы магнитной анизотропии от вектора намагниченности в приближении молекулярного поля? 8. Вектор спонтанной намагниченности Is вращается в плоскости (100) монокристалла кубической системы, ( - угол между Is и крисаллографичеким направлением [001]. Вывести выражение для магнитострикции в функции угла ( в направлении [111].