【电磁场计算】磁偶极子在外磁场中的受力

磁偶极子受力计算:从推导到解析
本文详细探讨了磁偶极子在外磁场中的受力计算,通过磁势能的梯度得出受力公式,并澄清了一些文献中的误解。推导过程中指出在无外电流条件下,受力简化为磁矩与磁场梯度的点乘。进一步通过矢量展开,明确了受力表达式的矩阵形式,揭示了与磁场雅可比矩阵的关系。最终得出简洁的受力表达式,并强调了与磁场梯度矩阵的区别。

磁偶极子在外磁场中的受力

背景

最近在总结电磁场计算公式的时候发现,关于磁偶极子受力的计算公式,在不同文献中的说明有点出入,在此特地写一篇博客澄清一下。

受力推导

通过磁偶极子在外磁场中的磁势能求梯度得出其受力:
F = − ∇ ( U ) = − ∇ ( − m ⋅ B ) = m × ( ∇ × B ) + ( m ⋅ ∇ ) B \begin{align*} F &= -\nabla(U) \\ &= -\nabla(-m \cdot B) \\ &= m \times (\nabla \times B) + (m \cdot \nabla)B \end{align*} F=(U)=(mB)=m×(×B)+(m)B
在计算的区域内(磁偶极子所在的微小空间)内可认为无外电流,即 ∇ × B = 0 \nabla \times B=0 ×B=0,故有:
F = ( m ⋅ ∇ ) B (1) F=(m \cdot \nabla)B \tag{1} F=(m)B(1)

以上是标准的计算公式,很多教材中也是这样写的。而有的文献中直接描述为 F F F等于磁矩 m m m与磁场梯度 ∇ B \nabla B B相乘—— F = m ⋅ ∇ B F=m \cdot \nabla B F=mB,这是不对的。

矢量展开

式(1)右侧是磁矩 m m m右点乘梯度算子 ∇ \nabla ,再对磁感应强度 B B B求梯度,展开为:
F = ( m x ⋅ ∂ ∂ x + m y ⋅ ∂ ∂ y + m z ⋅ ∂ ∂ z ) B = [ m x ∂ B x ∂ x + m y ∂ B x ∂ y + m z ∂ B x ∂ z m x ∂ B y ∂ x + m y ∂ B y ∂ y + m z ∂ B y ∂ z m x ∂ B z ∂ x + m y ∂ B z ∂ y + m z ∂ B z ∂ z ] = [ ∂ B x ∂ x ∂ B x ∂ y ∂ B x ∂ z ∂ B y ∂ x ∂ B y ∂ y ∂ B y ∂ z ∂ B z ∂ x ∂ B z ∂ y ∂ B z ∂ z ] [ m x m y m z ] = J ( B ) m \begin{align*} F &=(m_x \cdot \frac{ \partial}{ \partial x } + m_y \cdot \frac{ \partial}{ \partial y } +m_z \cdot \frac{ \partial}{ \partial z } )B \\ &= \begin{bmatrix} m_x \frac{ \partial B_x }{ \partial x } + m_y \frac{ \partial B_x }{ \partial y } + m_z \frac{ \partial B_x }{ \partial z } \\ m_x \frac{ \partial B_y }{ \partial x } + m_y \frac{ \partial B_y }{ \partial y } + m_z \frac{ \partial B_y }{ \partial z } \\ m_x \frac{ \partial B_z }{ \partial x } + m_y \frac{ \partial B_z }{ \partial y } + m_z \frac{ \partial B_z }{ \partial z } \end{bmatrix}\\ &= \begin{bmatrix} \frac{ \partial B_x }{ \partial x } & \frac{ \partial B_x }{ \partial y } &\frac{ \partial B_x }{ \partial z } \\ \frac{ \partial B_y }{ \partial x } & \frac{ \partial B_y }{ \partial y } & \frac{ \partial B_y }{ \partial z } \\ \frac{ \partial B_z }{ \partial x } & \frac{ \partial B_z }{ \partial y } & \frac{ \partial B_z }{ \partial z } \end{bmatrix} \begin{bmatrix} m_x \\ m_y \\ m_z \end{bmatrix} \\ &= J(B) m \end{align*} F=(mxx+myy+mzz)B=

在MATLAB中,绘制磁偶极子(magnetic dipole)的磁场分布图通常涉及到使用数学公式来计算磁场强度,并利用MATLAB的绘图工具如`meshgrid`、`surf`或`quiver`来可视化结果。以下是一个简化的步骤: 1. **定义磁场公式**:磁偶极子磁场公式为 \( \vec{B} = \frac{\mu_0}{4\pi} \frac{3(\vec{m} \cdot \hat{r})\hat{r} - \vec{m}}{r^3} \),其中 \(\vec{m}\) 是磁偶极矩,\(\hat{r}\) 是从磁偶极子到某点的位置向量,\(r\) 是距离,\(\mu_0\) 是真空中的磁导率。 2. **创建网格**:使用`meshgrid`函数创建x、y和z网格,以便对每个网格点计算磁场。 3. **计算磁场**:根据磁场公式,遍历网格点,计算每个点的磁场强度。 4. **绘制磁场**: - 使用`surf`函数创建一个表面图,如果想要更直观地显示磁场方向,可以使用`quiver`函数画出箭头表示磁场方向。 - 或者,你可以只画出等值线,使用`contourf`或`isosurface`函数。 ```matlab % 假设m = [x y z]是磁偶极矩,r = sqrt(x.^2 + y.^2 + z.^2)是距离 [x, y, z] = meshgrid(-10:0.1:10); % 创建网格 r = sqrt(x.^2 + y.^2 + z.^2); m = [1 0 0]; % 假设磁偶极矩在x轴上 B = (mu0/4/pi) * (3*(m.'*r)./r.^3 - m./r.^3); % 计算磁场 % 选择可视化方法 % surf(x, y, z, B) 或 quiver(x, y, z, real(B(:, :, 1)), real(B(:, :, 2)), real(B(:, :, 3))) 显示磁场 % contourf(x, y, B) 或 isosurface(x, y, z, B, 'EdgeColor', 'none') 显示等值面 % 清晰显示坐标轴 view(3) xlabel('X') ylabel('Y') zlabel('Z') % 提示:可能需要设置合适的颜色标度和轴范围以更好地展示磁场强度 ```
评论 3
成就一亿技术人!
拼手气红包6.0元
还能输入1000个字符
 
红包 添加红包
表情包 插入表情
 条评论被折叠 查看
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值