|
发表于 2006-3-18 09:43:47
|
显示全部楼层
双频感应淬火的计算与应用
摘 要:提出一种新的通用性更好的计算方法。该方法应用于实际生产的感应加热装置上,不仅可设计特殊需要的温度分布,而且可对温度场进行分析和模拟。试验结果表明双频加热淬火比单频加热淬火获得更好的表面硬化质量。尽管通常双频加热淬火比传统的单频加热淬火费用更昂贵,但在某些情况下证实还是节约的,在特殊情况下是唯一的工艺方法。 4 f+ Q0 e* v/ }8 Q9 }4 D5 B7 k
关键词:双频感应淬火;温度场;计算
" r$ O, t" R1 J- |- l' Q
8 U$ B9 ]6 `+ w2 W/ V) d1 b引 言
4 R0 e9 Z- b/ X% W3 e& r0 K
. S3 P3 p3 c6 r+ P) ^3 z! q( N, r双频淬火装置如图1,其包含两只并排的感应器分别供给不相同频率的电流。上感应器通入较低的频率,工件移动使加热深度较深(低频率);而下感应器通入较高的频率,工件移动使加热深度较浅,此深度相当于所需的硬化层厚度。
1 z, R# [; \& `- }6 Q
8 A- p$ s; Y9 V+ G0 [3 o1 工件; 2,3 频率分别为f1和f2的感应器; 4,5 磁芯(导磁体)
G" ^! a* x9 s2 `4 w+ d图1 双频淬火装置简图. r) k, f6 f' }
选择两只感应器的功率使上感应器加热工件到居里温度,下感应器在相同时间内沿着顺磁性的工件加热到淬火温度。要设计双频淬火装置的各种参数,有必要作系统的电磁计算,但在专门的文献中缺乏适当的计算方法介绍[1]。 6 K# f5 N. U! M
) I" a/ l! k$ [1 轧辊双频感应淬火的数学分析 ; N" X! V, c9 w, U
- B" B* E$ h* s" q
1.1 简化的假设和计算模型
* X6 R3 i. u9 O
% T% X; G8 b5 g求解双频感应淬火过程的电磁问题是一个复杂的数学问题,困难在于系统零件的有限尺寸、工件电磁参数的非线性变化以及感应器移动速度对特征现象的影响。充分考虑所有的这些因素,则解法复杂,以致于无法实用。于是在对解答的精确性影响不大的前提下,引入三类简化的假设: $ H0 C4 ?- n) s/ P9 L3 o t
: ~2 F r) [# @& y(1)有关工件-感应器系统的几何形状假设: 1 |. I0 i, S9 O' k; `7 [: a g7 k' X
$ H* _ V8 c. @0 @+ s1 `/ G
a.以无限长的轧辊代替有限长度的加热工件;
* ?: U% \2 @: n+ Ib.以有限高度和微小厚度的环代替有限尺寸感应器;
0 d# L( Y# B2 G# N) h, Ac.一个无限长的圆柱体代替两个有限尺寸感应器的磁芯。
+ l# G4 R; t7 C5 I: x" o' b2 X& F, o( H* {0 r7 R; z3 C
(2)有关工件电磁参数的假设:
% ?8 D* H; m$ N
! M1 ?2 `, s' _2 h* na.淬火加热过程可以分成两个阶段。第一阶段,铁磁性工件仅仅由频率f1的感应器在一定厚度内加热到居里温度(这厚度近似等于穿透深度即趋肤深度);第二阶段,频率为f2的感应器加热的工件可分两层,外层是顺磁性材料(μr=1),而内层是仍保持高导磁率的铁磁性材料,参数μ和σ在每一层的体积内作为一常数。 , f" w$ f. P: ?" V# \) s
b.两层中导磁率可由众所周知的方法[2]为基础计算得到,为简化起见内层的导磁率均作为无穷大; - Q0 A+ y5 b) G
c.把频率f1感应器加热得的表面铁磁性层厚度与频率f2感应器加热得的表面顺磁性层厚度看作相等。
* \: n& r. A% U& }' L8 E/ E H/ p6 [2 {; ]7 U* k
(3)为了避免积分方程式的问题而假设: ! }/ \1 \9 W2 o7 b1 [4 y$ _! t* d( } G
! d/ V% T4 Z8 Qa.不计感应器的移动速度; 7 _7 U6 E# k( v9 o& |# R1 h
b.把工件沿着坐标线Z=0分开,这样把每半只轧辊均看成无限长,照图2的两种情况分别计算,然后将得到的介于适当的区内结果进行合并。
8 p! w' P# E9 @# P0 V+ T8 `2 @: \" R8 |; u) ^
图2为考虑所有这些假设的计算模型。 0 U. Q/ M L9 X4 i( @
. i1 c/ r/ }, f" K+ `; G图2 工件“分开”后的计算模型
1 K" b- E4 d4 M2 k! m# A5 \1.2 矢热微分方程式的通解
6 D! M3 `6 i4 k9 b R0 R6 n& m3 x% H; B- N, u
从麦克斯韦尔(Maxwell)方程式计算推出矢势的微分方程式,在圆柱体坐标中有下列形式[3、4]:
8 O! Z; n2 U+ Q( L . o4 x0 ]' |9 U! Q1 a
式中:
5 ~3 Y& B/ F* ^' t8 {+ f5 jm2=jωμσ; / F- w* _0 f" D
ω=2πf电流脉动数(即角频率);μ——介质的导磁率;σ——介质的电导率 8 A6 i! e1 f; _( a% A. t
6 z' ^; m; o( U( Y2 O# r8 i( W电磁场矢量可用A以下列形式表示:
; l9 r( m' J. T: r+ x * {4 a* Y0 P7 d* e7 p* ?
式中: ' t8 J% V( E. d
E——电动势; - l. i/ R+ ] g/ \9 g( C8 K
B——磁感应强度;
! O E$ |7 d3 b4 H+ @: kJ——涡流密度 $ ^/ E/ G6 H1 T& f0 z& v; k
: q$ J# L j: ~方程式(1)的特解可由分离变量法求出,为 ) Z6 r, q. ?7 q5 W' p/ [& _( y
( H* v' s; ~5 z9 SA(r,z)=[C(k)I1(pr)+D(k)K1(pr)]]×[F(k)coskz+G(k)sinkz] (4)
1 H7 s" V( e/ v0 g0 K# j! D% @5 w. n9 S" ?4 R
式中:C(k),D(k),F(k),G(k)——积分常数 # M. x* B: d2 C( y4 Y% d
I1,K1,I0,K0——修正的贝塞尔(Bessel)函数 0 e; ?( g" Q" R1 h
k——分离变量常数
. f' t( z) }8 M9 ~2 Z$ O& e: M: A
4 w$ L) e! t) L9 [* v& pp=(k2+m2)1/2 (5) # U. g; X5 u8 t9 j) ~5 F
: _3 G. a: n+ X% ?/ R
方程式(1)的通解是(4)式的所有特解形式的总和,于是: 9 o; s% F+ P( j9 }
' c3 }% I" u. H( a) U9 `( X在图2的特殊计算区间中,(6)式的解有多种形式。事实上从应用的观点出发,最本质的是工件的外层(区域Ⅱ),仅仅这个区域解的结果如下: 0 h9 `* ]* L: a- u
, J" p2 W4 _0 [) _+ X
式中:p1=k2+jω1μ21σ21 8 \8 g- j" _) h7 ?- a
8 j' F4 v( L' T1 j; k3 R3 zp2=k2+jω2μ22σ22 , M" O0 m y# {/ X% }
; H% w( g$ ~6 ^# t$ ^# h2 {
符号Ⅱ表示计算的Ⅱ区域,而记号1和2分别表示频率f1[图2(a)的情况]和f2[图2(b)的情况]感应器的有关量。 2 b/ y5 J5 n! I1 [ t l
) y, p! y# d3 o ?/ J4 D" w表达式(7)中出现的积分常数可用矢势并符合图2的边界条件求出,有 ; ~+ a' p/ w% R5 |" V8 Z; f
' ]3 J/ A8 l; F" t( h7 b6 G; O* Q
式中:μi——在i区域中的导磁率 ; x( B' M7 r) H
+ q9 r. q6 r3 [3 a# Z2 [0 l& U4 u式中:I1,I2——频率f1和f2的感应电流; * H2 m, W! s- V3 q* P
n1,n2——频率f1和f2的感应器匝数
2 f- |) o7 \! i: K( bH,h——感应器1、2分别对应的工件长度(见图2)
" M+ }+ A$ ?( Y. ]. A
: S% r. C; r9 L4 {基于边界条件(9)算出积分常数,于是求得工件的矢势。 3 ^. i! M' q8 O0 z; @
3 i: h4 I, ]5 Q式中:N1=K1(kR3)I0(kR4)+I1(kR3)K0(kR4)
3 C* O5 U) i$ IN2=K0(kR2)I0(kR4)-I0(kR2)K0(kR4)
4 g4 @8 l8 k i3 BN3=K1(kR2)I0(kR4)+I1(kR2)K0(kR4)
" C* e+ u, r. `0 H2 J" }$ O8 YN4=I1(p1R2)K0(p1R1)+K1(p1R2)I0(p1R1)
! g: g/ d" y9 GN6=I0(p1R1)K0(p1R2)-K0(p1R1)I0(p1R2)
, F/ |) M: m2 q4 e8 F e/ dM4=I1(p2R2)K0(p2R1)+K1(p2R2)I0(p2R1)
. }# h9 `* G& k; x, gM6=I0(p2R2)K0(p2R1)-K0(p2R2)I0(p2R1) $ Y. f; T# m: r H' z
7 ~2 |/ f; a3 I" h4 Y) d
x1=p1/k;x2=p2/k (12) ! d' d8 a4 f4 U3 T; I' h4 d8 A" h
( p6 L4 G0 w. X6 O
μr1,μr2——在图2(a)和(b)情况中的相对导磁率
8 A. Z& d9 k M" T6 w' q) Q' _9 u H6 m6 K+ t" p& S! W8 H
1.3 工件中磁感应强度、涡流密度和有效功率密度 6 L: s/ t/ l2 q4 H: [& w* m
- L7 n7 w2 z" S# ]8 h% e区域Ⅱ中与(3)式相一致的磁感应强度有二个分量BrII和BzII,按式(11)可得下式:
$ R! ~/ q. v, S+ j8 Y 3 b8 k0 ^( \! {) T& T9 ^6 {! I
由于矢势(11)式的微分无实际意义,这也就没有给出磁感应强度的表达公式,但轧辊表面(r=R2)的磁感应强度分量(Bz)通过(13)式借助于计算机算出,结果示于图3中。 & X: N; h' {- z
9 h7 @$ |5 a \+ ?$ `工件上的涡流密度按(3)式可得:
$ b5 F$ C: I+ L/ i, Q& V8 I* g* ~ ( p& X+ P3 d" x& F' _
轧辊表面(r=R2)的涡流密度按上式计算,结果示于图4中。
. u! V$ [) s; l9 ^& t, H( W% ?/ c: S- W7 F. h/ M
工件表面的有效功率密度可按下式计算: + m6 J$ G; g3 w
' X+ p% f% P+ V- g' S' L7 y
由(14)得:
& v( G9 g# p" O' t, \ |$ l 1 X! g i8 X3 p2 X* x/ `, `8 P
按式(16)计算的有效功率密度示于图5中。 9 ^. v# _$ I, ~3 L0 W# c
1 y$ a" {1 m7 L8 j# l# |0 U9 q/ h2 轧辊双频淬火试验
: R I: b: ~" r% Q0 X! T
' d$ n& x5 q( b( W根据上述计算,上感应器选用50 Hz,下感应器保持1000 Hz频率。轧辊双频淬火试验的同时,对轧辊的温度分布进行了测量,结果见图6. ) {& @2 w5 }3 E; ^
; E8 i3 L& u7 k5 j9 J1,2,3,4,5,6离表面不同的距离的温度分布
4 i' H7 U2 J. U! o图6a 双频感应淬火时的温度分布
. _$ b8 F9 x( @8 N o! w' g& v
4 |9 \% I N: V" r: Z2 V( _(离轧辊表面的距离) 9 b8 T: B' @2 [1 M3 P
A淬火加热始瞬间温度;B通过50 Hz感应器后的温度; . F [4 S- Y/ L& Y! X; ^0 z' o
C通过1000 Hz感应器后的温度;D喷雾的瞬间温度 " |, H! c0 z2 u% P0 h
图6b 轧辊的截面Ⅰ的温度分布
8 J' Y- ?: ]! P" W2 r由图6a可知,淬火过程的温度的分布有三个最高点,第一个最高点(t=3 min)相应于预热后的温度,第二个最高点(t=21 min)相应于通过50 Hz感应器后的温度,第三个最高点(t=25 min)相应于通过1000 Hz感应器后的温度。图还表明,适当的淬火温度开始于工件温度250 ℃后的第16分钟,精确来说大约加热仅10分钟,第27分钟后喷雾,处理即结束。淬火中轧辊的温升平均速度大约是1.2 ℃/sec,辊身截面的温度分布见图6b.与单频淬火得的温度分布相比,可以看出:轧辊表面在深度30 mm内的温度差,后者约210 ℃,而前者仅120 ℃,从而保证得到更厚的硬化层和较小的硬度梯度(图7). ; b% H! a1 {% G/ {( p
, t) C) h0 p! ]9 @" P(离轧辊表面的距离)
1 T2 Y' z' v: Z. _& l1一次预热的单频感应淬火;
5 V) i" z0 _7 }1 _9 o$ f; n& Y2一次预热的双频感应淬火 3 S8 s* X4 T( x% V T2 _% W8 g
图7 轧辊的截面硬度分布
/ e2 h% x) Z8 B9 ]- a4 E! L1 T4 E0 T' Q3 结 论
4 S* u- K& z& X
5 s, [" [& @1 o/ f& c: O+ g1 x- r1.双频感应淬火的上述计算方法具有较强的通用性和实用性。
' }2 d; o: \9 M$ X' C' S |2.双频淬火较单频淬火能使工件获得更好的性能(表面硬度,截面硬度梯度、硬化深度),并提高工效和降低能耗。 |
|