Задача 1. Постройте силовые линии магнитного поля, созданного: 1) одной обмоткой с током; 2) двумя обмотками, по которым текут токи в одном направлении или в противоположных направлениях.

Рассмотрите виток с током радиусом R, лежащий в плоскости XOY. Разобьем его на элементы и введем векторы сонаправленные с током. Расположение этого элемента витка и точки A, в которой определяется индукция магнитного поля, задаются векторами:

По закону Био-Савара-Лапласа для вектора индукции и его проекций можно записать:

При этом используются формулы:

Используется программа ПР-1. Для построения силовой линии задают начальную точку A₀, определяют в ней направление индукции магнитного поля, строят небольшой отрезок в направлении вектора индукции, переходя в новую точку A₁. Затем повторяют все снова. Результаты вычислений - на рис. 1, 2, 3.

Программа ПР-1.

Рис. 1. Магнитное поле обмотки

Рис. 2. Магнитное поле двух обмоток. Токи сонаправлены

Рис. 3. Токи в обмотках направлены противоположно

Задача 2. Постройте силовые линии магнитного поля, созданного тремя параллельными проводниками с током.

Рис. 4. Расчет магнитного поля двух токов

Допустим, два проводника с токами I1 и I2 пересекают плоскость XOY в точках с координатами (x1,y1) и (x2,y2). Индукция магнитного поля в точке A(x,y) может быть найдена по формулам:

Используемая программа ПР-2 случайным образом выбирает точку A₀, определяет проекции вектора индукции и в его направлении строит отрезок единичной длины, находя точку A₁. Затем все повторяется снова. Выполнив 10000 шагов, программа снова случайно выбирает точку A₀ и начинает строить вторую силовую линию и т.д.

Программа ПР-2.

Рис. 5. Магнитное поле трех проводников с током

Задача 3. Создайте компьютерную модель ферромагнетика Изинга. Промоделируйте процесс образования областей спонтанной намагниченности (доменов).

Рассмотрите двумерную модель магнетика. Представьте себе сетку, в каждом из узлов которой находится атом, спин (собственный магнитный момент) которого может быть направлен от нас, к нам или лежать в плоскости рисунка. Энергия взаимодействия соседних атомов друг с другом -Jsi,jsi+1,j, атомы, расположенные далеко друг от друга, не взаимодействуют. Внешнее магнитное поле способствует повороту спина атома и изменяет энергию системы на H или -H. Энергия всей системы равна:

Процесс моделирования состоит в следующем: 1) задается начальная ориентация всех спинов случайным образом; 2) определяется энергия системы; 3) случайно выбирается один атом, и его спин изменяется на противоположный, определяется энергия системы в новой конфигурации; 4) если в результате переориентации спина энергия уменьшилась, то такое изменение конфигурации принимается, в противном случае - отклоняется; 5) результат выводится на экран; 6) возврат к операции 2. Для решения этой задачи используется программа ПР - 6.5, при этом моделируется ситуация, в которой в левом верхнем углу индукция магнитного поля направлена от нас, а в правом нижнем углу - к нам. Используется не квадратная, а шестиугольная сетка, результат вычислений - на рис. 6.5.

Программа ПР-3.

Рис. 6. Модель ферромагнетика Изинга

Задача 4. Ферромагнитный стержень вставили в обмотку из N витков, которая подключена к источнику тока. Увеличивают силу тока i в обмотке, измеряют индукцию B магнитного поля при данной напряженности H=iN. Затем уменьшают силу тока, пропускают ток в противоположном направлении, увеличивая и затем уменьшая его до нуля. Создайте компьютерную модель намагничивания и перемагничивания ферромагнетика, получите кривую намагниченности.

Рассмотрим двумерную модель ферромагнетика, состоящую из прямоугольной сетки N x M (кристаллической решетки), в узлах которой находятся атомы с магнитным моментом p. Изначально магнитное поле отсутствует, ферромагнетик полностью размагничен, магнитные моменты атомов ориентированы произвольно, образуя с осью Ox случайные углы αij. Если ось Ox сонаправлена с силовой линией магнитного поля, то среднее значение индукции:

где H -- напряженность магнитного поля, α_ij -- угол между вектором магнитного момента p и осью Ox. Индукция магнитного поля B_ij в данном узле сетки обусловлена внешним полем обмотки μ₀H и магнитным полем, создаваемым соседними атомами. В проекции на ось Ox:

Со стороны магнитного поля на атом, находящийся в узле (i, j), действует вращающий момент M_VR. При повороте атома возникает момент "упругих сил", препятствующий повороту. В результате действия этих моментов, а также сил внутреннего трения, атом поворачивается на Δα. Запишем уравнения:

Здесь x -- небольшая случайная величина в интервале [-0,1; 0,1], учитывающая хаотические колебания атомов в узлах кристаллической решетки. Давая случайное приращение x углу α, мы "встряхиваем" нашу модель, предоставляем возможность найти более "удачные" значения α. Это важно еще и потому, что в анализируемой модели число атомов существенно меньше, чем в реальном ферромагнитном стержне.

Суть вычислительного эксперимента состоит в плавном увеличении напряженности магнитного поля H и пересчете α_ij у всех атомов с нахождением соответствующих средних значений индукции магнитного поля B. В результате на плоскости B от H получится кривая намагниченности с характерным горизонтальным участком, соответствующим магнитному насыщению, когда магнитные моменты всех атомы ориентированы по полю. Для получения петли гистерезиса необходимо, чтобы напряженность поля H изменялась по гармоническому закону.

Используемая программа ПР-4 содержит цикл по времени t, в котором пересчитываются H, α_ij, B_SR и строится зависимость B(H).

Рис. 7. Результаты моделирования магнитного гистерезиса.

Программа ПР-4.

Задача 5. Методами компьютерного моделирования изучите процесс намагничивания ферромагнетика, получите кривую намагниченности и петлю гистерезиса. Как изменятся результаты, если уменьшить амплитуду колебания тока через обмотку.

Необходимо запустить программу PR-4 и подождать пока на экране строится петля гистерезиса (рис. 7.1). Шаг по времени должен быть достаточно мал. При уменьшении амплитуды колебаний напряженности поля получается петля, изображенная на рис. 7.2.

Задача 6. Методом компьютерного моделирования исследуйте зависимость формы петли гистерезиса от величины магнитного момента p атомов данного ферромагнетика.

Если уменьшить магнитный момент атомов p в несколько раз (от 3,4 до 1,4), то уменьшится индукция магнитного поля, соответствующая насыщению, петля гистерезиса сожмется по вертикали (рис. 7.3).

Тексты программ находятся в zip-архиве, файл gl-61.pas.

Майер, Р. В. Задачи, алгоритмы, программы / Р. В. Майер [Электронный ресурс]. - Глазов: ГГПИ, 2012 // Web-site http://maier-rv.glazov.net .