Ответ на вопрос
Коротко — причины и модели.
Причины появления вихревых токов
- Закон Фарадея: изменение магнитного потока через замкнутые контуры в пластине индуцирует электрическое поле
\(\nabla\times\mathbf{E}=-\dfrac{\partial\mathbf{B}}{\partial t}\).
- Движение проводника в поле даёт сдвиг магнитного потока (присутствует мотонная ЭДС) через член \(\mathbf{v}\times\mathbf{B}\) в законе Ома для движущегося проводника
\(\mathbf{J}=\sigma(\mathbf{E}+\mathbf{v}\times\mathbf{B})\).
- Направление индуктированных токов подчиняется закону Ленца: создаваемые токи противодействуют изменению потока.
- Неоднородное поле и геометрические замыкания (плоские области, окна, края) формируют замкнутые локальные пути тока → вихревые токи концентрируются в зонах градиента поля или плотных проводящих областях.
- Частотный/динамический эффект: при высокой частоте токи «забираются» в поверхность — скин‑эффект (см. ниже).
Ключевые формулы (для расчёта полей, сил и потерь)
- Ома (движущийся проводник): \(\mathbf{J}=\sigma(\mathbf{E}+\mathbf{v}\times\mathbf{B})\).
- Сила на объём: \(\mathbf{f}=\mathbf{J}\times\mathbf{B}\); общая сила/момент: \(\mathbf{F}=\displaystyle\int_V \mathbf{J}\times\mathbf{B}\,dV\).
- Плотность тепловыделения (Джоуля): \(p=\mathbf{J}\cdot\mathbf{E}=\dfrac{\mathbf{J}^2}{\sigma}\).
- Магнитная диффузия (квазистационарно, \(\varepsilon\) пренебрежимо): \(\nabla^2\mathbf{B}=\mu\sigma\,\dfrac{\partial\mathbf{B}}{\partial t}\).
- Скин‑глубина (гармонический): \(\delta=\sqrt{\dfrac{2}{\mu\sigma\omega}}\).
- Простой «стержень/рамка» (оценка тормозящей силы): для подвижной длины \(L\) в поле \(B\) индуцированная ЭДС \(\mathcal{E}=B L v\), ток \(I=\dfrac{\mathcal{E}}{R}\), сила \(F=I L B=\dfrac{B^2 L^2 v}{R}\), мощность \(P=I^2R=\dfrac{B^2 L^2 v^2}{R}\).
- Для гармонических полей с индуктивностью \(L\): \(I(\omega)=\dfrac{\mathcal{E}(\omega)}{R+j\omega L}\), средняя мощность \(P=\dfrac{|\mathcal{E}|^2 R}{2(R^2+(\omega L)^2)}\).
Как смоделировать влияние (практический план)
1. Выбрать приближение
- Магнито‑квазистационарное (дисперсия смещение невелико): дискредитирует \(\varepsilon\), решают для векторного потенциала \(\mathbf{A}\) или поля \(\mathbf{H}\).
- Если частота низкая и толщина мала — применять тонкоплёночное (thin‑sheet) приближение / поверхностные токи.
2. Решение электромагнитной части
- Уравнения: \(\nabla\times(\nu\nabla\times\mathbf{A})+\sigma\frac{\partial\mathbf{A}}{\partial t}=\mathbf{J}_{\text{source}}-\sigma\mathbf{v}\times(\nabla\times\mathbf{A})\) (формализм A‑вариант), либо H‑формулировка.
- Граничные условия: внешние поля, симметрии, замыкание токов.
- Инструменты: FEM (COMSOL, Ansys Maxwell, Opera), BEM/аналитика для простых геометрий.
- Тонкие элементы/ламинаты: использовать эффективное удельное сопротивление или специальные элементы для тонких пластин.
3. Вычисление сил и моментов
- Вычислить \(\mathbf{J}\) и затем \(\mathbf{F}=\int_V \mathbf{J}\times\mathbf{B}\,dV\). Для роторных систем интегрировать момент \(\tau=\int_V \mathbf{r}\times(\mathbf{J}\times\mathbf{B})\,dV\).
- Для быстрых оценок использовать эквивалент R‑L‑контур и формулы выше; при малых скоростях поведение ≈ вязкое: \(F\approx -c v\).
4. Тепловое сопряжение
- Источник в уравнении теплопереноса: \(Q(\mathbf{x},t)=\dfrac{\mathbf{J}^2}{\sigma}\).
- Тепловое уравнение: \(\rho c_p \dfrac{\partial T}{\partial t}-\nabla\cdot(k\nabla T)=Q\).
- Учесть температурную зависимость \(\sigma(T)\), теплообмен с окружением.
5. Механическая связь
- Подключить вычисленную ЭДС/силы в уравнения движения: \(m\ddot{x}=-\int_V (\mathbf{J}\times\mathbf{B})_x\,dV + F_{\text{внеш}}\) (или вращательный эквивалент).
- Для динамики использовать явную/неявную интеграцию, учесть обратную связь поля от положения.
6. Верификация и упрощения
- Сравнить с аналитикой для простых случаев (скользящий стержень, круглая пластина в однородном поле).
- При необходимости учитывать скин‑эффект и магнитное насыщение (\(\mu(B,T)\)).
Практические инженерные выводы и меры управления эффектом
- Торможение/демпфирование: вихревые токи превращают кин. энергию в тепло — их можно использовать (электромагнитные тормоза, демпферы) или уменьшать.
- Нагрев: рассчитывается через \(Q=\int_V \mathbf{J}^2/\sigma\,dV\); важна теплоотдача и температурная зависимость \(\sigma(T)\).
- Стабилизация: в малых колебаниях вихревые токи дают вязкую демпфировку \(F\propto -v\) — полезно для пассивной стабилизации (маятники, гиростабилизаторы).
- Снижение нежелательных вихревых токов: ламинирование, прорези/слоты, повышение удельного сопротивления (композиты), уменьшение толщины (меньше площади замыкающих контуров), магнитная экранировка, изменение геометрии поля.
- Управление для тормозов: проектировать R/L и геометрию так, чтобы достичь желаемой зависимости силы от скорости и теплового режима.
Короткое правило масштаба
- Силы и потери масштабируются примерно как \(B^2\), пропорциональны объёму активной проводящей ткани и зависят линейно от проводимости и квадратично от скорости в простых случаях: \(P\sim \dfrac{B^2 L^2 v^2}{R}\).
Если нужно, могу привести конкретную простую модель (прямоугольная пластина: геометрия, выражение \(R\), оценка силы и тепла) или подсказать настройки для FEM (формулировка, сетка, тонкий слой).
Еще