Реальный соленоид: точный расчет поля и визуализация результатов

4 часа назад Прокомментировать Просмотры: 7

Несколько соленоидов с током

Поставим перед собой следующую задачу: пусть задан соленоид высотой 2h, имеющий внешний радиус R2 и внутренний радиус R1 (где R1 < R2). В центральную полость помещен цилиндрический стальной сердечник аналогичной высоты 2h и радиуса R1. Поверх сердечника плотно, виток к витку и строго перпендикулярно оси, намотан круглый медный провод диаметром d, причем эта величина пренебрежимо мала по сравнению с габаритами самого соленоида. По обмотке пропускается постоянный ток I (также стоит проанализировать случай переменного тока I(t)), циркулирующий по часовой стрелке. Требуется вычислить индукцию магнитного поля B(x,y,z, t) во всех ключевых зонах: внутри системы, за ее пределами и непосредственно в толще сердечника.

Решение:Разместим начало координат O правой декартовой ортогональной системы XOYZ точно в геометрическом центре цилиндра.Выполним чертеж:

Перейдем к цилиндрической системе координат:Пусть орт i

ориентирован в сторону проекции рассматриваемой точки на плоскость XOY,

орт j направлен по часовой стрелке ортогонально i в пределах плоскости XOY, а орт k

совпадает с положительным направлением оси OZ. Принятые орты формируют правую тройку векторов.

Закон Био-Савара-Лапласа для единичного витка — Закон Био‑Савара‑Лапласа для единичного витка

Внутри сердечника (при r0 < R1 и ‑h < z0 < h) магнитная проницаемость равна μ стали, внутри самой обмотки соленоида (при R1 < r0 < R2 и ‑h < z0 < h) принимается μ меди, в то время как за пределами этих элементов (при r0 > R2 или выходе за границы по z0: z0 < ‑h либо z0 > h) средой считается воздух с проницаемостью μ воздуха.

Индукция магнитного поля демонстрирует обратную кубическую зависимость от толщины провода, приближенно пропорционально 1/d^3.

Иными словами, уменьшение сечения проводника влечет за собой мощное усиление магнитного поля.

При переходе к переменному току вида I = I(t), результирующая индукция B(t) будет синхронно следовать изменениям управляющего сигнала I(t).

Ниже представлен скрипт на языке Python, разработанный при поддержке нейросети для графической демонстрации пространственного распределения поля:

import numpy as np import matplotlib.pyplot as plt from mpl_toolkits.mplot3d import Axes3D from matplotlib import cm Физические константы mu0 = 4 np.pi 1e-7 # Тл·м/А Параметры соленоида R1 = 0.01 # внутренний радиус (м) R2 = 0.02 # внешний радиус (м) h = 0.04 # половина высоты (м) d = 0.005 # диаметр провода (м) I = 0.001 # сила тока (А) Функция для численного интегрирования def calculate_B(r0, z0): Константа для расчетов const = mu0 * I / (np.pi * d ** 2) # Пределы интегрирования r_min = R1 - d / 2 r_max = R2 + d / 2 # Интеграл I1 I1 = 0 for r in np.linspace(r_min, r_max, 100): for z in np.linspace(-h, h, 100): denom = ((r - r0) ** 2 + (z - z0) ** 2) ** 1.5 I1 += r * (z0 - z) / denom # Интеграл I2 I2 = 0 for r in np.linspace(r_min, r_max, 100): for z in np.linspace(-h, h, 100): denom = ((r - r0) ** 2 + (z - z0) ** 2) ** 1.5 I2 += r * (r - r0) / denom # Компоненты поля Br = const * I1 Bz = const * I2 return Br, Bz Создаем сетку для визуализации r = np.linspace(0, 2 R2, 50) z = np.linspace(-2 h, 2 * h, 50) R_grid, Z_grid = np.meshgrid(r, z) Рассчитываем поле Br_grid = np.zeros_like(R_grid) Bz_grid = np.zeros_like(Z_grid) for i in range(R_grid.shape[0]): for j in range(R_grid.shape[1]): Br_grid[i, j], Bz_grid[i, j] = calculate_B(R_grid[i, j], Z_grid[i, j]) Визуализация fig = plt.figure(figsize=(14, 8)) Радиальная компонента plt.subplot(121) plt.pcolormesh(R_grid, Z_grid, Br_grid, cmap=cm.coolwarm) plt.colorbar(label="Br (Тл)") plt.title('Радиальная компонента магнитного поля') plt.xlabel('r (м)') plt.ylabel('z (м)') plt.axis('equal') Осевая компонента plt.subplot(122) plt.pcolormesh(R_grid, Z_grid, Bz_grid, cmap=cm.coolwarm) plt.colorbar(label="Bz (Тл)") plt.title('Осевая компонента магнитного поля') plt.xlabel('r (м)') plt.ylabel('z (м)') plt.axis('equal')

plt.tight_layout() plt.show()

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

Таким образом, в рамках данного исследования были выведены аналитические соотношения для расчета магнитного поля многовиткового соленоида при постоянном и переменном питании путем интегрирования отклика одиночного витка по всему объему обмотки, а также подготовлен эффективный инструмент визуализации на Python.

Литература: