基于Python仿真载流圆线圈的磁场分布

1、理论基础

R是圆环的半径：
x是距离圆环中心的距离：
$$dB=\frac{{\mu }_0}{4\pi }\frac{Icos\alpha dl}{r^2}$$
$$B=\frac{({\mu }_{0}I{R}^2)}{2(x^2+R^2{)}^{\frac{3}{2}}}$$

仿真载流圆线圈的磁场分布：
设场点P的位置为 ，电流元位置为 ，电流元矢量 ,由此可得场点P相对于电流的位置矢量为：
利用行列式计算Idl×r ，可写为
设载流线圈通过坐标原点垂直于x-y平面放置，电流元Idl沿z轴正向，设x= 0，则y=Rsin(θ) , z=Rcos(θ)，$$z_0=0,dx=0,dy=Rcos(\theta )d_{\theta },dz=Rsin(\theta )d_{\theta }$$
即：
$$dl=[d_x{d}_y{d}_z]=[0Rcos(\theta )d_{\theta }Rsin(\theta )d_{\theta }]$$
$$r=[x_0-x,y_0-y,z_0-z]$$
通过Python求叉乘

import numpy as np
import matplotlib.pyplot as plt
from pylab import meshgrid,  arange, streamplot, show
import mpl_toolkits.mplot3d
from mpl_toolkits.mplot3d import Axes3D

import sympy as sym

dx, dy, dz, x0, x, y0, y, z0, z, R,theta = sym.symbols('dx, dy, dz, x0, x, y0, y, z0, z, R,theta')  # 声明符号变量，否则会认为没有定义

z0=0
dx=0
x=0
y=sym.simplify(R*sym.sin(theta))
z=sym.simplify(R*sym.cos(theta))

dy=sym.diff(y,theta)
dz=sym.diff(z,theta)
dl=[dx, dy,dz]
r=[x0-x,y0-y, z0-z]
dlcr=np.cross(dl,r) 
print(dlcr)          

x=sym.Symbol('x')
f1=sym.sin(x)/x
f2=(1+x)**(1/x)
f3=(1+1/x)**x
#三个参数是 函数，变量，趋向值
lim1=sym.limit(f1,x,0)
lim2=sym.limit(f2,x,0)
lim3=sym.limit(f3,x,sym.oo)
print(lim1,lim2,lim3)

求矢量的叉乘积可得到：

r的大小为：

由毕奥－萨伐尔定律
$$dB⃗=\frac{{\mu }_0}{4\pi }\frac{(Idl⃗\times r⃗)}{r^3}$$
有：

2、算法仿真

在篇文章里有详述

3、实验结果

通过编程可得到其磁感线的分布：