DOI:10.19476/j.ysxb.1004.0609.2002.s1.005
调幅交变磁场下圆柱形金属液内电磁力分布
上海大学上海市钢铁冶金新技术应用开发重点实验室
上海大学上海市钢铁冶金新技术应用开发重点实验室 上海200072
摘 要:
从Maxwell方程出发 , 推导并计算了连续交变磁场和调幅交变磁场作用下 , 圆柱形金属液中电磁力密度在半径方向上的分布及其随时间的变化情况。发现在连续及调幅磁场情况下 , 电磁力频率均以原输入电磁场频率的两倍变化 , 且在调幅磁场的线性上升段和下降段 , 电磁力幅值的变化近似为二次抛物线形 , 在调幅磁场的恒值段则与连续施加时一样。这为合理设计电磁连铸时的调幅磁场打下基础。
关键词:
中图分类号: TG249.9
收稿日期:2001-10-09
基金:国家自然科学基金资助项目 ( 5 97340 80 , 5 99740 17);
Distribution of electromagnetic force in cylinder shape liquid metal under amplitude modulated alternating magnetic field
Abstract:
The distribution of electromagnetic force along with radius and time in cylinder shape liquid metal under continuous and amplitude modulated alternating magnetic field was deduced and calculated based on Maxwell equation. It is found that the frequency of electromagnetic force is the twice of the magnetic field under both continuous and amplitude modulated condition, and its amplitude changes in a square parabola type while amplitude of magnetic field changes in a linear type.
Keyword:
amplitude modulated magnetic field; electromagnetic continuous casting; electromagnetic force;
Received: 2001-10-09
软接触结晶器电磁连铸技术是一项新的钢坯连铸技术设想。 它将高频电磁场透过结构特殊的结晶器, 作用到结晶器内的液态金属上, 使铸坯的表面及内部质量都得到显著的改善
目前关于电磁场改善铸坯表面质量机理的研究, 主要集中在电磁场对铸坯初始凝固行为的影响上。 研究表明, 在该技术中, 电磁场透过结晶器壁作用在液态金属上, 使之与模壁的接触减少, 实现“软接触”的状态。 这种状态的实现, 使初生坯壳和结晶器之间的保护渣道得以拓宽, 从而减小了因结晶器振动所导致的保护渣道内动态压力的变化
显然, 由电磁场产生的电磁力在金属液中的分布, 以及在施加间断磁场下电磁力在时域上的变化情况, 对于精确控制连铸坯的初始凝固行为, 具有极为重要的作用, 而弄清它们的规律, 具有很大的意义。 到现在为止, 已系统研究过的间断磁场类型仅限于方波的间断磁场, 并且均关注的是在实验条件下结晶器内弯月面的波动行为
1连续交变磁场下的电磁力
计算依据的物理模型如图1所示。 计算中采用众多热模拟实验均采用的液态金属Sn, 其电阻率为0.200 μΩ·m, 磁化率为-45×10-11m·ks
1) 螺旋管线圈为无限长直螺旋管。
2) 不考虑在电磁力作用下金属液形成的类“弯月面”的自由面。
3) 不考虑液态金属容器壁对磁场的屏蔽作用, 因为屏蔽作用仅改变磁场强度的大小, 并不影响其分布。
图1 计算物理模型
Fig.1 Physical model for computation
当水冷螺旋管线圈中通以一定频率ω的交变电流时, 就会在液态金属的侧面产生一个大小为B0cosωt的磁场, 其中B0在线圈中没有铁磁性物质的情况下正比于线圈中的最大电流Imax。 同时, 液态金属内部磁感应强度的分布满足方程
由于本问题的对称性及长直螺旋管的特性, 可知B在径向及角向的分量均为0, 即Br=Bθ=0, 令轴向分量Bz (r, t) =Re[Bz (r) ·exp (jωt) ], 代入控制方程中, 分离变量后得到:
此方程为零阶的标准Kelvin方程, 又称为变型的Bessel方程
和
得
综上所述, 在液态金属液内磁感应强度的分布为
式中 B0为与最大电流相对应的金属液侧面的磁感应强度,
由交变的电磁场所产生Lorentz体积力密度F由电流密度和磁感应强度的叉积得到, 其只有径向的分量Fr, 即
F=J×B=[Fr, 0, 0]=[F2ω+Fave, 0, 0] (6)
其中
式中
C0=ber
C
bei1 (kr) ]2+[bei0 (kR) ·ber1 (kr) -
bei0 (kR) ·ber1 (kr) ]2
C
bei0 (kR) ]2+[bei (kr) ·ber0 (kR) -
ber0 (kr) ·bei0 (kR) ]2
ζ=cos-1{[ber0 (kr) ·ber0 (kR) +
bei0 (kr) ·bei0 (kR) ]/C2}
θ=cos-1{[ber0 (kR) ·ber1 (kr) +
bei0 (kR) ·bei1 (kr) ]/C1}-π/4
在上述诸式中, beri (x) 和beii (x) 为i阶的Kelvin函数,
图3所示为B0=0.008 T, r/R=0.9时, 在电磁场频率为100 Hz和1 000 Hz时计算的电磁力密度随时间的变化情况。 由图可知, 在各种频率的电磁场下, 电磁力密度均以磁场频率的两倍频率变化。
图2 各种频率下电磁力在半径方向上的分布
Fig.2 Distributions of electromagnetic force along with radius in several frequencies
图3 连续施加电磁场时电磁力随时间变化
Fig.3 Electromagnetic force changes with time under continuous magnetic field
2 调幅磁场下电磁力
在施加连续的恒幅电磁场的基础上, 人们提出了施加幅值随时间变化的调幅电磁场的构想, 相当于在一较高频率的磁场上 (称为载波) 附加一频率很低周期性调幅波 (称为调制波) 。 仍取图1所示的物理系统, 如果载波磁场可写为B0cosω0t, 调制波写为f (t) 的形式, 则调幅磁场的数学描述为B0cosω0t·f (t) , 其中f (t) 为一周期函数。
一般地, 当f (t) 满足Dirichlet收敛定理时, 周期为2l的函数f (t) 可以写成Frourier级数的形式, 将之代入B0cosω0t·f (t) 变换后可以得到调幅磁场的级数表现形式为
式中
(n=0, 1, 2, …)
(n=1, 2, 3, …)
可以看出, 调幅磁场除了幅值为a0B0/2的原有频率ω0的谐波成分外, 还有频率分别为ω0+nπ/l, ω0-nπ/l的两个系列的谐波, 且与原载波磁场同相位时幅值为anB0/2, 与原载波磁场相位相差+π/2或-π/2时幅值为bnB0/2。 不妨令ω1n=ω0+nπ/l, ω2n=ω0-nπ/l, B00=a0B0/2, B1n=anB0/2, B2n=bnB0/2, 则式 (8) 可以简写成
对于式 (9) 中的每一项, 在液态金属内部, 均满足方程 (1) , 同理, 对每一项按式 (2) ~ (5) 解出相应的磁感应强度的分布, 叠加后即得到调幅磁场在液态金属内产生的磁感应强度的分布情况:
式中
同样, 对应于式 (10) 中的每一项, 都会感应出一相应的电流密度, 叠加后就得调幅磁场作用下在液态金属中感生的电流密度:
最后得到调幅磁场在液态金属内产生的Lorentz力的分布为
F (r, t) =J×B={Fr, 0, 0}
={Jθ (r, t) ·Bz (r, t) , 0, 0} (12)
以上的分析对于任何波形、 任何频率的调幅磁场均适用, 所以具有一定的普遍性。 当采用较规则的调幅波形时, 式 (10) 和 (11) 中的待定系数可以很容易地求得。 依据上面的推导, 分别计算了方波、 三角波、 锯齿波、 梯形波和正弦波5种情况, 它们的相关特性列于表1中。
正弦波调幅磁场与其它波形有显著的区别, 它只包含频率为ω1, ω2的两种谐波, 且其幅值均为B0/2。 作者选择梯形波调幅磁场进行计算, 因为它既有三角波和锯齿波调幅磁场一样的线性上升和下降部分, 还有和矩形波一样的恒值部分。 同前, 取B0=0.008 T, r/R=0.9, 梯形波的特性如表1中所示。 计算过程是先根据式 (10) 计算出0.9倍半径处的磁感应强度大小, 再据式 (11) 算出该点处的电流密度, 即可以根据式 (12) 得到该点的电磁力密度。 图4所示为载波磁场频率为20 Hz和100 Hz的梯形波调幅磁场下电磁力密度在调制波一个周期内的变化情况。 由图可知, 在调幅磁场的恒值部分, 电磁力密度随时间的变化情况与连续磁场时一样。 但在磁场的线性上升和下降部分, 电磁力密度的变化不再是线性的了, 其包络近似为二次抛物线。 同样的计算表明, 在更高频率时的情况与图4中的规律完全一样。 由图还可知道, 在高频载波磁场的作用下, 金属液中所受到的电磁力以两倍于其的频率的正值和负值间变化, 但正值的绝对值要大于负值的。 对电磁连铸中金属液中说, 有意义的是电磁力在一个变化周期内的“平均作用效果”。 设电磁力的变化周期为T, 定义在交变磁场下一个周期的平均作用力
调幅磁场在液态金属中产生电磁力的这种分布特性, 为正确设计电磁连铸中磁场施加方式提供了依据。 比如, 在连铸结晶器振动过程中, 因为保护渣道内的压力导致初始凝固壳的周期性形变, 产生振痕
表1 五种调幅磁场及其数学特性
Table 1 Five kinds of amplitude modulated magnetic field and their mathematical feature
Parameter | Square wave | Triangle wave | Zigzag wave | Trapezoid wave | Sinusoidal wave |
B00 |
|||||
B1n |
0 | 0 | 0 | ||
B2n |
0 | 0 | 0 | ||
ω1n |
ω0+ (4n-2) π | ω0+ (4n-2) π | ω0+2nπ | ω0+2nπ | ω1=ω0+2nπ |
ω2n |
ω0- (4n-2) π | ω0- (4n-2) π | ω0-2nπ | ω0-2nπ | ω2=ω0-2nπ |
图4 不同频率梯形波调幅磁场下电磁力在调制波一个周期内的变化
Fig.4 Electromagnetic force changement in a modulated wave cycle under different frequency trapezoid amplitude modulated magnetic field
图5 各种波形调幅磁场作用下金属液内平均电磁力随时间的变化
Fig.5 Change of average electromagnetic force in molten metal with time under several kinds of amplitude modulated magnetic field
3 结论
从Maxwell方程出发, 推导并计算了连续交变磁场和调幅交变磁场作用下, 圆柱形金属液中电磁力密度在半径方向上的分布及其随时间的变化情况。 发现:
1) 在连续及调幅磁场情况下, 电磁力频率均以原输入电磁场载波频率的两倍变化。
2) 在调幅磁场的情况下, 平均电磁力的变化规律与输入的磁场的幅值频率变化相同, 但在调幅磁场幅值的一次上升段和下降段, 电磁力幅值的变化近似为二次抛物线形, 在调幅磁场在恒值段则与连续施加时一样。
参考文献