基于永磁体离心量对磁悬浮热解石墨片光驱转动现象的研究 下载: 801次
1 引言
1842年英国物理学家Earnshaw[1]提出了磁悬浮的概念,之后关于磁悬浮的理论日趋成熟,并应用于各个领域。抗磁悬浮技术是抗磁材料因在外部施加的磁场下受力而克服其自身重力并稳定悬浮于磁场中的一种技术。热解石墨作为抗磁材料,有着体积磁化率高、密度低、质量轻、负磁化率系数大等特点,可以在常温下悬浮于永磁体材料(NdFeB永磁铁等)的磁场下。目前磁悬浮热解石墨已经有了多种应用[2-4]。近几年日本科学家Kobayashi和Abe [5]通过研究发现在圆柱形嵌套的永磁体上方磁悬浮的热解石墨片可以在激光操控下运动。永磁体布置包括一块实心的圆柱状NdFeB永磁铁和一块与其等高的空心环形状NdFeB永磁铁,且实心磁铁以N、S极相反的方式放置于圆环形永磁铁中。热解石墨片的直径略大于实心永磁铁的外径,这可保证热解石墨能稳定悬浮于系统中而不至于滑落。当激光照射到圆形石墨片边缘处时石墨片开始转动,当激光照射至圆形石墨片左侧时石墨片顺时针旋转,照射至右侧时圆形石墨片逆时针旋转。在此之后Niu等[6]研究了方形石墨片在方形永磁体上方磁悬浮的稳定条件,发现方形石墨片在方形永磁体阵列产生的磁场中不具有旋转不变性,这导致石墨片表面的受力分布不均匀,使得石墨片产生具有方向性的旋转,其扭矩与旋转角度呈三角函数关系。
本文对圆形热解石墨片的旋转问题进行了进一步的分析。由之前的研究结果可知,热解石墨片的受力分布不均是造成其旋转的主要原因。由于圆形热解石墨片本身具有旋转不变性,故本文认为是磁场分布的不对称造成了圆形热解石墨片的旋转。在实际装配嵌套的永磁体时,由于制造误差和磁极间的相互吸引,内外磁铁之间会产生一定的离心量,该离心量会造成磁场分布的不对称。当激光沿着石墨片边缘处照射时,磁场分布的不对称性使得石墨片在磁场中受到的合力存在切向分力,该分力造成了圆形石墨片具有方向性的旋转。
2 理论计算
若想要计算抗磁物质磁悬浮的稳定条件,需要分别对外部磁场的空间分布和抗磁物质在外部磁场中的受力分布进行计算。
计算外部磁场空间分布的主要方法有等效磁荷法和等效电流法。电流模型主要用于计算螺旋线圈附近磁场的空间分布,磁荷模型主要用于计算矩形、环形等简单永磁体附近的磁场分布。本文主要计算永磁体产生的磁场分布,故采用等效磁荷法进行计算[7-9]。由磁场分布公式可得永磁体的空间磁场分布为[10]
式中:B为永磁体产生的磁感应强度;H为永磁体产生的磁场强度;M0为永磁体的剩磁强度,由于永磁体沿磁化方向均匀充磁,故只存在面磁荷σm ,M0为常矢量;μ0为空气的磁导率。本文以永磁体上表面中心点为原点,上表面为xoy平面构造坐标系,则磁场强度H可表示为
式中:S为永磁体面磁荷组成的边界面;μp为永磁体磁导率;r为参考系原点到坐标空间计算点的距离。接着由空气与永磁体的边界条件可得
式中:n为永磁体边界面的外法线单位向量。在具体计算中,μ0,M0,μp为已知量,首先将(3)式代入(2)式求得空间磁场强度分布,接着将(3)式代入(1)式求得空间磁感应强度分布。
接着便可以计算抗磁物质在磁场中稳定悬浮的条件。其中抗磁物质的受力分布可表示为
式中:f为抗磁物质单位体积在磁场中受到的力;V为抗磁物质的体积。根据文献[ 11]:
式中:M为抗磁物质在磁场中的磁化强度;
式中:B为永磁体产生的磁感应强度大小;Fx、Fy、Fz分别为抗磁物质所受的力沿x、y、z轴方向的分力;χx、χy、χz分别为抗磁物质沿着x、y、z轴方向的分量。
当抗磁物质稳定悬浮时,其沿坐标轴z轴方向的分力Fz等于其重力G,则有
式中:ρ为抗磁物质的密度;g为重力加速度。将(10)代入(5)式,解得
由(11)式可得抗磁物质稳定悬浮的位置。Bz分量可以通过联立(1)~(3)式求得,接着沿z轴方向取适当的空间步长,通过计算得到dB/dz。
已知抗磁物质的扭矩公式为[12]
式中:rg是抗磁物质中心点到计算点的位置向量。由(12)式可以看出,扭矩主要由两部分构成,其中第一部分
根据抗磁物质的磁场势能公式[11]:
可以计算得到抗磁物质在空间磁场中的势能分布。根据势能分布可以确定抗磁物质的稳定位置。
基于上述公式,对方形热解石墨片在方形永磁体上方的磁悬浮状态进行计算分析,得到了与文献[ 6]相似的结果,这验证了利用等效磁荷法分析石墨磁悬浮的有效性。根据方形石墨片的受力分布规律,发现方形热解石墨片在不同旋转角度区间内的受力分布不对称是造成其旋转的主要原因,这为下面分析圆形热解石墨片的光驱转动现象提供了基础。
3 圆形热解石墨抗磁悬浮以及光驱转动分析
本节对圆形热解石墨片磁悬浮的光驱旋转机理进行分析。在进行圆形热解石墨片磁悬浮研究时,首先需要对嵌套的圆环形永磁体进行装配。为了便于装配,外部圆环磁铁的内径必定大于内部圆环磁铁的外径,由于异性相吸的作用,这两个磁铁之间将产生一定的离心量。该离心量将导致磁铁磁场分布的不对称性。本研究将分析激光照射圆形热解石墨片边缘时,在这种离心量所带来的磁场不对称性下圆形热解石墨片在磁场中的受力分布特性,从而研究圆形热解石墨片发生方向性旋转的规律。
第2节介绍了等效磁荷法的计算公式,本节将使用该方法进行计算。计算模型如
图 1. 圆形热解石墨片仿真模型示意图
Fig. 1. Diagram of simulation model of circular pyrolytic graphite disk
采用等效磁荷法计算永磁体空间磁场分布,如
图 2. 环形嵌套永磁体磁场分布。(a)(d) Bx; (b)(e) By; (c)(f) Bz
Fig. 2. Magnetic field distribution of circular nested permanent magnets. (a)(d) Bx; (b)(e) By; (c)(f) Bz
由
图 3. 圆形石墨片稳定悬浮处受力分布。(a)(c) Fx; (b)(d) Fy
Fig. 3. Force distribution of circular graphite sheet for stable suspension. (a)(c) Fx; (b)(d) Fy
根据Kobayasi的实验现象,激光分别照射在热解石墨片表面边缘处对称位置时,石墨片的旋转方向相反。根据热解石墨材料的物理性质,其磁化率具有各向异性,磁化率数值随温度的升高而下降。当采用60mW功率的激光照射至热解石墨片表面时,石墨片温度变化可根据热传导方程式计算得出,如
图 4. 激光照射热解石墨表面时的温度变化
Fig. 4. Temperature change of pyrolytic graphite surface irradiated by laser
根据(12)、(13) 式,可以计算激光照射下圆形热解石墨片的扭矩以及势能变化。
图 5. 激光照射下圆形热解石墨片扭矩分布图。(a)激光照射下圆形石墨片磁悬浮示意图;(b) Δx=0mm时扭矩分布。
Fig. 5. Torque distribution of circular pyrolytic graphite disk under irradiation of laser. (a) Schematic illustration of magnetic suspension of circular graphite disk under irradiation of laser; (b) torque distribution when Δx=0mm
当环形嵌套永磁体存在离心量且激光沿着热解石墨片表面边缘照射时,石墨片受力分量Fx、Fy如
图 6. Δx=0.15mm时热解石墨片受力分布示意图
Fig. 6. Force distributions of pyrolytic graphite disk when Δx=0.15mm
根据(13)式计算圆形热解石墨片的势能与照射点的关系,如
图 7. Δx=0.15mm时激光照射下热解石墨片扭矩及磁势能分布。(a)势能分布;(b)扭矩分布
Fig. 7. Torque and potential energy distributions of pyrolytic graphite disk under irradiation of laser when Δx=0.15mm. (a) Potential energy distribution;(b) torque distribution
4 结论
首先介绍了等效磁荷法,并通过计算得到永磁体空间磁场分布,并基于此计算了热解石墨抗磁悬浮的稳定条件。接着使用该方法计算分析了方形热解石墨片在方形永磁体上方的磁悬浮稳定条件,分析得到方形热解石墨片在不同旋转角度区域内的所受合力的切向分量方向不同是造成其具有方向性旋转的主要原因,并基于此思路针对圆形热解石墨片的光驱转动现象进行进一步分析。分析发现永磁体装配时产生的离心量会造成磁场分布的不对称,进而导致圆形热解石墨片受力分布不均匀。当激光照射在石墨片表面时,温度的变化引起照射区域石墨磁化率的改变,使得圆形热解石墨片在该处的受力分布发生变化。当环形嵌套的永磁体严格同轴时,激光照射下的圆形热解石墨片所受合力沿半径方向,不存在切向分量,扭矩恒为零。当永磁体存在偏心量时,磁场分布不对称,圆形热解石墨片产生的合力存在切向分量,使得扭矩在(0°,180°)和(180°,360°)两个区间内的数值相反,进而引起了圆形热解石墨片的方向性转动。此研究解释了圆形热解石墨片光驱转动的现象,并为之后进行磁悬浮转盘激光器的研究提供了理论基础。
[1] Earnshaw S. On the nature of the molecular forces which regulate the constitution of the luminiferous ether[J]. Transactions of the Cambridge Philosophical Society, 1842, 7: 97- 112.
[2] Pelrine R. Diamagnetic levitation[J]. American Scientist, 2004, 92(5): 428-435.
[4] Simon M D, Heflinger L O, Geim A K. Diamagnetically stabilized magnet levitation[J]. American Journal of Physics, 2001, 69(6): 702-713.
[6] Niu C, Lin F, Wang Z M, et al. Graphene levitation and orientation control using a magnetic field[J]. Journal of Applied Physics, 2018, 123(4): 044302.
[7] 李世鹏, 张卫平, 陈文元, 等. 抗磁悬浮的原理、特点和应用[J]. 磁性材料及器件, 2011, 42(1): 1-6, 44.
[8] 李景天, 宋一得, 郑勤红, 等. 用等效磁荷法计算永磁体磁场[J]. 云南师范大学学报(自然科学版), 1999, 19(2): 3-5.
Li J T, Song Y D, Zheng Q H, et al. Computation of the magnetic field of permanent magnet with equivlent magnetic charge method[J]. Journal of Yunnan Normal University (Natural Sciences Edition), 1999, 19(2): 3-5.
[9] 苏宇锋, 叶志通, 张坤. 抗磁悬浮石墨转子理论分析与实验[J]. 中国机械工程, 2017, 28(9): 1039-1043.
Su Y F, Ye Z T, Zhang K. Theoretical analyses and experiments of micromachined graphite rotor based on diamagnetic levitation[J]. China Mechanical Engineering, 2017, 28(9): 1039-1043.
[10] Campbell P, Chari M. D’Angelo J. Three-dimensional finite element solution of permanent magnet machines[J]. IEEE Transactions on Magnetics, 1981, 17(6): 2997-2999.
[11] Nguyen J, Contera S, Llorente García I. Magneto-electrical orientation of lipid-coated graphitic micro-particles in solution[J]. RSC Advances, 2016, 6(52): 46643-46653.
[12] Arora R K. Forces and torques on current circuits placed in nonuniform magnetic fields[J]. IEEE Transactions on Magnetics, 1989, 25(3): 2789-2790.
Article Outline
童鑫, 唐锋, 李建郎. 基于永磁体离心量对磁悬浮热解石墨片光驱转动现象的研究[J]. 激光与光电子学进展, 2021, 58(1): 0126001. Tong Xin, Tang Feng, Li Jianlang. Optical Drive Rotation of Magnetically Levitating Pyrolytic Graphite Disk Based on Eccentricity of Permanent Magnets[J]. Laser & Optoelectronics Progress, 2021, 58(1): 0126001.