首页 / 专利库 / 数学与统计 / 层间位移 / 一种涂层感应加热方法

一种涂层感应加热方法

阅读:438发布:2024-01-13

专利汇可以提供一种涂层感应加热方法专利检索,专利查询,专利分析的服务。并且本 发明 涉及一种涂层 感应加热 方法,在感应加热过程中,管体、线圈与空气都属于不同介质,具有不同的磁热特性,根据不同介质的材料磁热特性,整个作用域空间被划分为三个个作用区域,感应加热过程中涂层 熔化 的 热能 来源于基体的感应 电流 产生的 焦 耳 热 ,因此,对感应 磁场 与 电场 的耦合分析成为感应加热过程理论分析的核心,在研究加 热管 体材料机械特性保护方面,热应 力 场分析是辅助热场分析的重要手段;优点是:加热速度快,加热 精度 高, 温度 控制稳定,极大的提高了工作效率,节省成本。,下面是一种涂层感应加热方法专利的具体信息内容。

1.一种涂层感应加热方法,包括边界作用域划分,电-磁耦合分析,热-场耦合分析,其特征在于:1)在感应加热过程中,管体、线圈与空气都属于不同介质,具有不同的磁热特性,根据不同介质的材料磁热特性,整个作用域空间被划分为三个个作用区域,分别为Ω1:
管体区,为生热涡流区,材料为导电与导热介质,无外加激励电场,在此区域电流电荷被考虑,自由电荷产生热加热基体,Ω2:线圈区,为交变电流激发区,材料为导电材料,在此区域感应电流激励下产生交变磁场,给基体和提供激发磁场,给基体激发感应电流,Ω3:空气区,为非生热涡流区,材料为非导电介质,无外加激励源,在此区域发生辐射、热对流
2)感应加热过程中涂层熔化热能来源于基体的感应电流产生的焦耳热,因此,对感应磁场与电场的耦合分析成为感应加热过程理论分析的核心,电磁场分析中,首先引用经典电磁理论的麦克斯韦方程组:
▽×D=ρ
▽×B=0


其中,D为电位移,C/m2,B为磁感应强度,Gs,E为电场强度,N/C,H为磁场强度,Oe,J为电流密度,A/m2,ρ为电荷密度,C/m2,
焦耳热产生于自由电荷,电位移密度 不是产生焦耳热的原因,只有在频率达到
10MHz时,才会产生焦耳热,再者,因为电流密度J远远大于电位移密度 所以不给计算在内,麦克斯韦方程组转化为:
▽×B=0    (2-6)

▽×H=J    (2-8)
在此过程中的辅助协调公式如下:
B=μH    (2-9)
J=σE(欧姆定律)    (2-10)
其中,μ为磁导率,σ为电导率,
由磁感应强度B满足等式(2-6)的散度为零的,根据亥姆霍兹定理引入磁向量势A,得到等式
B=▽×A    (2-11)
代入公式(2-11)到(2-7)中,可以得到:

因 满足亥姆霍兹定理,由此引入标量电位势 公式(2-12)可以进一步
推导为:

公式(2-11)与公式(2-13)可以同时满足公式(2-6)与公式(2-7),因此代入辅助协调方程公式(2-9)与公式(2-10)到公式(2-8)中,可得

标量电位势 与磁向量势A不是唯一的,对于不同的应用规范会有不同的对应值,洛伦兹规范和库仑规范常用于解决电磁协调分析方程,在是准静态近似条件下,相对于洛伦兹规范利用库伦规范可极大的简化模型方程,在库伦规范协应用后,标量势与矢量势都可以被表达为标量电位势 与磁向量势A的函数,因此,库伦规范表达式▽·A=0被代入公式(2-14),惩罚公式 引入公式(2-14),得到公式

引入旋度计算公式▽×▽×A=▽(▽·A)-▽2A到公式(2-15),可得

代入公式(2-10)与(2-12)可得

其中,Je为感应电流密度,Js为源电流密度,表达为
对于谐波分析,电磁量是单一频率的振动函数,因此A可以写成:
A=A0ejωt    (2-18)
Js=J0ejωt    (2-19)
其中,J0为源电流密度的振幅,A0为磁势振幅,ω为频率,可表达为ω=2πf,由于构建模型的对称性,外加激励电流可以设计为:


把公式(2-18)与公式(2-19)代入公式公式(2-16),可得

对于轴对称的问题,柱坐标系下的公式(2-22)转化为

以上公式适用于不同介质区域的电磁场,包括磁体、非磁体、导电介质和非导电介质,根据感应加热建模区域建模介质不同,由公式(2-13)、(2-15)与(2-17),各个区域描述方程如下:
管体区Ω1范围内描述方程:



线圈区Ω2范围内描述方程:



空气区Ω3范围内描述方程:
J=0    (2-30)

针对于不同区域边界介质的特性,设定各划分区域的边界条件如下:
Ω1法线方向上:
n·(D2-D1)=σe0    (2-32)
n·(B2-B1)=0    (2-33)

Ω1径线方向上:
n×(E2-E1)=0    (2-35)
n×Hout=i0 or n×(H2-H1)=0    (2-36)
Ω2:法线方向上:
n·(J2-J1)=0    (2-37)
Ω3:法线方向上:
n·(B2-B1)=0    (2-38)
径向方向上:
n×(H2-H1)=0    (2-39)
其中,Di为电流密度,Bi为磁流密度,Ji为电流密度,Ei(i=1,2)为感应涡流电场强度,n为单位法线向量,σe0为自由电荷密度,i0为交界表面法相方向上的电流密度,Hout为磁场密度,方向沿着界面指向外侧,特别是在超高频率下公式n×Hout=i0代替公式n×(H2-H1)=0作为边界条件,
在感应加热过程中管体被磁化后,其表面应力张量的变化即可产生电磁力,电磁力的计算引用公式Q=0.24I2Rt与公式(2-17),可得磁场强度控制方程:
▽×H=Js+Je    (2-40)
为了对方程(2-40)进行求解,需要将方程转化为弱解形式,由此引入引入爱因斯坦张量表达式:

其中,εijk为符号变量,i,j,k为下标,取得下标组合(1,-1,0),根据公式(2-41),两个变量的叉积即为
(ν×ω)i=εijkνjωk    (2-42)
其中,v与w为任意代数向量,对应的旋度公式为:

由公式(2-42)和公式(2-43),方程转化成爱因斯坦张量表达式:

其中,Hk为磁场强度分量,Ai为磁动势分量,A为磁路中的磁动势矢量,σe为磁介质电导率,
对公式(2-44)运用加权残值法积分可得:

其中:ψi为测试函数;VB为磁场的空间积分域。公式(2-45)中包含磁场强度H对位置变量x的散度,以及散度与测试函数ψi的乘积,可以进一步通过分部积分运算性质进行化简,得到

方程(2-46)中乘积项Hkψi对位置变量x的偏导数的体积积分实际上是乘积Hkψi的散度计算,根据散度计算定理,一个矢量V散度的体积分等于V在此体积表面上的法向分量的面积分,所以方程(2-46)可以进一步表示为

其中,n表示包围积分空间域,VB表面的法相向量,其指向向外, 表示沿体积VB的表面进行积分。将方程(2-46)代入式(2-47)可得

公式(2-48)即为爱因斯坦标记下三维磁场分布计算方程。为方便方程的求解,需要将公式(2-48)写成矩阵形式,根据爱因斯坦标记法的运算规则,对方程(2-48)进行矩阵标记的转换:

式中磁场空间中边界处磁场强度的切向分量满足HT=H×n,对于变分变量,磁感应强度B与磁动势A满足下列表达式:
δB=▽×δA    (2-50)
则方程(2-50)的弱解形式可以进一步简化为

公式(2-51)即为弱解形式的电磁场偏微分方程,其中自变量为磁动势矢量A,从公式(2-51)的等号右边可以看到,系统的输入为外加电流密度与系统边界条件的矢量和;
3)在研究加热管体材料机械特性保护方面,热应力场分析是辅助热场分析的重要手段,由此,基于热场分析的热-应力耦合分析成为研究感应加热的重要理论依据,感应加热过程涉及管体热能传导的问题,现应用经典热传导的傅里叶定律,得到如下等式:

其中,ρ为材料密度,kg/m3;C为比热,J/(kg·℃);T为温度时间函数,℃;k热导系数,W/(m·℃);Q为感应涡流产生的热量,J,
感应加热的热能由感应涡流产生,结合焦耳定律,可以得到:

感应加热为非线性瞬态热传导过程问题,在笛卡尔坐标系内,建立三维瞬态热传导控制方程:

感应加热过程基体不具有内热源,因此qv=0,上式转化为Fourier方程形式:

对于工件温度的边界条件可以描述为法向导数,根据热力学第一定律表面的对流和热辐射可以表达如下:

其中,n为边界上指向外侧发的法向向量,Tair为工件的环境温度,℃;h为对流换热系数,W/(m2·℃);εemi为计算辐射的发射率,k为热导率,σb为斯特藩-玻尔兹曼常数或斯特藩常量5.67×10-8W·m-2·K-4,
感应加热过程对于不同区域边界上边界表达式各有不同,其具体边界分析如下:
Ω1管体区与Ω3空气区、Ω1管体区与Ω3空气区在传导热流量为Q时的边界Γ1:

Ω1管体区与Ω3空气区、Ω1管体区与Ω3空气区对流换热为h(T-T0)时的边界Γ2:

所有感应加热划分区域高温辐射换热为εσ(T4-T04)的边界Γ3:

感应加热过程开始阶管体受到热应力和电磁力的影响,当管体温度达到居里温度时,部分失磁的管体只受到热应力作用,管体未失磁部分受到磁场和热应力的双重作用,在感应加热成形完成后,失去感应加热激励电流,管体在迅速降温过程中,温度分布的不均匀性使得基体和涂层受到温度梯度作用,在涂层和基体分别产生一定的应力集中,感应加热热应力场分析主要集中在管体的降温阶段,由此,管体不受到外力和磁场力作用,仅在感应电流加热引起的温度梯度作用下,满足虚功原理条件,在体积为V,边界为S的连续介质中,t+Δt时刻的虚功积分表达式为:
∫Vσ·δε·dV=∫Vq·δu·dV+∫Sp·δu·dS    (2-60)
其中,σ为修正的拉格朗日应力,ε为和应变张量,q为体积力,p为面积力矢量,δu为虚位移,
考虑符合Mises屈服准则的热弹塑性材料本构模型如下式:

式中,为等效应力;σ’为应力偏移量分量;T为温度; 为材料等效塑性应变;对于热弹塑性过程,总的应变增量可分解为弹性应变εije、塑性应变εijp和热应变εijth

若令Dijkl为弹性应力应变关系系数,则从弹性的应力应变关系得出应力变化率

说明书全文

一种涂层感应加热方法

技术领域

[0001] 本发明涉及涂层喷涂领域,具体涉及一种涂层感应加热方法。

背景技术

[0002] 传统的管材面漆喷涂过程中,需要对管材进行离线地炉内加热,但是取出之后,难以使其保温在最佳温度附近,喷涂过程中的温度波动严重严重影响了面漆的附着质量,而且加热炉的加热效率严重影响了后续管材的喷涂效率,限制了生产节拍的提高。而作为一种新型的加热方式,感应加热可以实现无接触加热,并有着加热速度迅速,加热精度高等优点,并被广泛应用于工业生产当中。另外,在喷涂过程中,传统的喷头运动速度不稳定,喷头的启停动作、电网电压的波动影响了涂层的涂覆效果,涂层出现了众多的螺纹面,影响了涂层厚度的均匀性。因此通过研究温度稳定、运动精确的喷涂系统,对于提高钻杆防腐涂层的稳定性和可靠性有着重要的经济和技术意义。

发明内容

[0003] 为解决现有技术的不足,本发明提供一种涂层感应加热方法。
[0004] 本发明的技术方案是:一种涂层感应加热方法,包括边界作用域划分,电-磁耦合分析,热-场耦合分析,1)在感应加热过程中,管体、线圈与空气都属于不同介质,具有不同的磁热特性,根据不同介质的材料磁热特性,整个作用域空间被划分为三个个作用区域,分别为Ω1:管体区,为生热涡流区,材料为导电与导热介质,无外加激励电场,在此区域电流电荷被考虑,自由电荷产生热加热基体,Ω2:线圈区,为交变电流激发区,材料为导电材料,在此区域感应电流激励下产生交变磁场,给基体和提供激发磁场,给基体激发感应电流,Ω3:空气区,为非生热涡流区,材料为非导电介质,无外加激励源,在此区域发生辐射、热对流
[0005] 2)感应加热过程中涂层熔化热能来源于基体的感应电流产生的焦耳热,因此,对感应磁场与电场的耦合分析成为感应加热过程理论分析的核心,电磁场分析中,首先引用经典电磁理论的麦克斯韦方程组:
[0006]
[0007]
[0008]
[0009]
[0010] 其中,D为电位移,C/m2,B为磁感应强度,Gs,E为电场强度,N/C,H为磁场强度,Oe,J为电流密度,A/m2,ρ为电荷密度,C/m2,
[0011] 焦耳热产生于自由电荷,电位移密度 不是产生焦耳热的原因,只有在频率达到10MHz时,才会产生焦耳热,再者,因为电流密度J远远大于电位移密度 所以不给计算在内,麦克斯韦方程组转化为:
[0012]
[0013]
[0014]
[0015] 在此过程中的辅助协调公式如下:
[0016] B=μH(2-9)
[0017] J=σE(欧姆定律)(2-10)
[0018] 其中,μ为磁导率,σ为电导率,
[0019] 由磁感应强度B满足等式(2-6)的散度为零的,根据亥姆霍兹定理引入磁向量势A,得到等式
[0020]
[0021] 代入公式(2-11)到(2-7)中,可以得到:
[0022]
[0023] 因 满足亥姆霍兹定理,由此引入标量电位势 公式(2-12)可以进一步推导为:
[0024]
[0025] 公式(2-11)与公式(2-13)可以同时满足公式(2-6)与公式(2-7),因此代入辅助协调方程公式(2-9)与公式(2-10)到公式(2-8)中,可得
[0026]
[0027] 标量电位势 与磁向量势A不是唯一的,对于不同的应用规范会有不同的对应值,洛伦兹规范和库仑规范常用于解决电磁协调分析方程,在是准静态近似条件下,相对于洛伦兹规范利用库伦规范可极大的简化模型方程,在库伦规范协应用后,标量势与矢量势都可以被表达为标量电位势 与磁向量势A的函数,因此,库伦规范表达式 被代入公式(2-14),惩罚公式 引入公式(2-14),得到公式
[0028]
[0029] 引入旋度计算公式 到公式(2-15),可得
[0030]
[0031] 代入公式(2-10)与(2-12)可得
[0032]
[0033] 其中,Je为感应电流密度,Js为源电流密度,表达为
[0034] 对于谐波分析,电磁量是单一频率的振动函数,因此A可以写成:
[0035] A=A0ejωt  (2-18)
[0036] Js=J0ejωt  (2-19)
[0037] 其中,J0为源电流密度的振幅,A0为磁势振幅,ω为频率,可表达为ω=2πf,[0038] 由于构建模型的对称性,外加激励电流可以设计为:
[0039]
[0040]
[0041] 把公式(2-18)与公式(2-19)代入公式公式(2-16),可得
[0042]
[0043] 对于轴对称的问题,柱坐标系下的公式(2-22)转化为
[0044]
[0045] 以上公式适用于不同介质区域的电磁场,包括磁体、非磁体、导电介质和非导电介质,根据感应加热建模区域建模介质不同,由公式(2-13)、(2-15)与(2-17),各个区域描述方程如下:
[0046] 管体区Ω1范围内描述方程:
[0047]
[0048]
[0049]
[0050] 线圈区Ω2范围内描述方程:
[0051]
[0052]
[0053]
[0054] 空气区Ω3范围内描述方程:
[0055] J=0  (2-30)
[0056]
[0057] 针对于不同区域边界介质的特性,设定各划分区域的边界条件如下:
[0058] Ω1法线方向上:
[0059] n·(D2-D1)=σe0  (2-32)
[0060] n·(B2-B1)=0  (2-33)
[0061]
[0062] Ω1径线方向上:
[0063] n×(E2-E1)=0  (2-35)
[0064] n×Hout=i0or n×(H2-H1)=0  (2-36)
[0065] Ω2:法线方向上:
[0066] n·(J2-J1)=0  (2-37)
[0067] Ω3:法线方向上:
[0068] n·(B2-B1)=0  (2-38)
[0069] 径向方向上:
[0070] n×(H2-H1)=0  (2-39)
[0071] 其中,Di为电流密度,Bi为磁流密度,Ji为电流密度,Ei(i=1,2)为感应涡流电场强度,n为单位法线向量,σe0为自由电荷密度,i0为交界表面法相方向上的电流密度,Hout为磁场密度,方向沿着界面指向外侧,特别是在超高频率下公式n×Hout=i0代替公式n×(H2-H1)=0作为边界条件,
[0072] 在感应加热过程中管体被磁化后,其表面应力张量的变化即可产生电磁力,电磁力的计算引用公式Q=0.24I2Rt与公式(2-17),可得磁场强度控制方程:
[0073]
[0074] 为了对方程(2-40)进行求解,需要将方程转化为弱解形式,由此引入引入爱因斯坦张量表达式:
[0075]
[0076] 其中,εijk为符号变量,i,j,k为下标,取得下标组合(1,-1,0),根据公式(2-41),两个变量的叉积即为
[0077] (ν×ω)i=εijkνjωk  (2-42)
[0078] 其中,v与w为任意代数向量,对应的旋度公式为:
[0079]
[0080] 由公式(2-42)和公式(2-43),方程转化成爱因斯坦张量表达式:
[0081]
[0082] 其中,Hk为磁场强度分量,Ai为磁动势分量,A为磁路中的磁动势矢量,σe为磁介质电导率,
[0083] 对公式(2-44)运用加权残值法积分可得:
[0084]
[0085] 其中:ψi为测试函数;VB为磁场的空间积分域。公式(2-45)中包含磁场强度H对位置变量x的散度,以及散度与测试函数ψi的乘积,可以进一步通过分部积分运算性质进行化简,得到
[0086]
[0087] 方程(2-46)中乘积项Hkψi对位置变量x的偏导数的体积积分实际上是乘积Hkψi的散度计算,根据散度计算定理,一个矢量V散度的体积分等于V在此体积表面上的法向分量的面积分,所以方程(2-46)可以进一步表示为
[0088]
[0089] 其中,n表示包围积分空间域,VB表面的法相向量,其指向向外, 表示沿体积VB的表面进行积分。将方程(2-46)代入式(2-47)可得
[0090]
[0091] 公式(2-48)即为爱因斯坦标记下三维磁场分布计算方程。为方便方程的求解,需要将公式(2-48)写成矩阵形式,根据爱因斯坦标记法的运算规则,对方程(2-48)进行矩阵标记的转换:
[0092]
[0093] 式中磁场空间中边界处磁场强度的切向分量满足HT=H×n,对于变分变量,磁感应强度B与磁动势A满足下列表达式:
[0094]
[0095] 则方程(2-50)的弱解形式可以进一步简化为
[0096]
[0097] 公式(2-51)即为弱解形式的电磁场偏微分方程,其中自变量为磁动势矢量A,从公式(2-51)的等号右边可以看到,系统的输入为外加电流密度与系统边界条件的矢量和;
[0098] 3)在研究加热管体材料机械特性保护方面,热应力场分析是辅助热场分析的重要手段,由此,基于热场分析的热-应力耦合分析成为研究感应加热的重要理论依据,感应加热过程涉及管体热能传导的问题,现应用经典热传导的傅里叶定律,得到如下等式:
[0099]
[0100] 其中,ρ为材料密度,kg/m3;C为比热,J/(kg·℃);T为温度时间函数,℃;k热导系数,W/(m·℃);Q为感应涡流产生的热量,J,
[0101] 感应加热的热能由感应涡流产生,结合焦耳定律,可以得到:
[0102]
[0103] 感应加热为非线性瞬态热传导过程问题,在笛卡尔坐标系内,建立三维瞬态热传导控制方程:
[0104]
[0105] 感应加热过程基体不具有内热源,因此qv=0,上式转化为Fourier方程形式:
[0106]
[0107] 对于工件温度的边界条件可以描述为法向导数,根据热力学第一定律表面的对流和热辐射可以表达如下:
[0108]
[0109] 其中,n为边界上指向外侧发的法向向量,Tair为工件的环境温度,℃;h为对流换热系数,W/(m2·℃);εemi为计算辐射的发射率,k为热导率,σb为斯特藩-玻尔兹曼常数或斯特藩常量5.67×10-8W·m-2·K-4,
[0110] 感应加热过程对于不同区域边界上边界表达式各有不同,其具体边界分析如下:
[0111] Ω1管体区与Ω3空气区、Ω1管体区与Ω3空气区在传导热流量为Q时的边界Γ1:
[0112]
[0113] Ω1管体区与Ω3空气区、Ω1管体区与Ω3空气区对流换热为h(T-T0)时的边界Γ2:
[0114]
[0115] 所有感应加热划分区域高温辐射换热为εσ(T4-T04)的边界Γ3:
[0116]
[0117] 感应加热过程开始阶管体受到热应力和电磁力的影响,当管体温度达到居里温度时,部分失磁的管体只受到热应力作用,管体未失磁部分受到磁场和热应力的双重作用,在感应加热成形完成后,失去感应加热激励电流,管体在迅速降温过程中,温度分布的不均匀性使得基体和涂层受到温度梯度作用,在涂层和基体分别产生一定的应力集中,感应加热热应力场分析主要集中在管体的降温阶段,由此,管体不受到外力和磁场力作用,仅在感应电流加热引起的温度梯度作用下,满足虚功原理条件,在体积为V,边界为S的连续介质中,t+Δt时刻的虚功积分表达式为:
[0118] ∫Vσ·δε·dV=∫Vq·δu·dV+∫Sp·δu·dS  (2-60)
[0119] 其中,σ为修正的拉格朗日应力,ε为和应变张量,q为体积力,p为面积力矢量,δu为虚位移,
[0120] 考虑符合Mises屈服准则的热弹塑性材料本构模型如下式:
[0121]
[0122] 式中,为等效应力;σ’为应力偏移量分量;T为温度; 为材料等效塑性应变;对于热弹塑性过程,总的应变增量可分解为弹性应变εije、塑性应变εijp和热应变εijth[0123]
[0124] 若令Dijkl为弹性应力应变关系系数,则从弹性的应力应变关系得出应力变化率[0125]
[0126] 本发明的有益效果是:加热速度快,加热精度高,温度控制稳定,极大的提高了工作效率,节省成本。

具体实施方式

[0127] 一种涂层感应加热方法,包括边界作用域划分,电-磁耦合分析,热-力场耦合分析,1)在感应加热过程中,管体、线圈与空气都属于不同介质,具有不同的磁热特性,根据不同介质的材料磁热特性,整个作用域空间被划分为三个个作用区域,分别为Ω1:管体区,为生热涡流区,材料为导电与导热介质,无外加激励电场,在此区域电流电荷被考虑,自由电荷产生焦耳热加热基体,Ω2:线圈区,为交变电流激发区,材料为导电材料,在此区域感应电流激励下产生交变磁场,给基体和提供激发磁场,给基体激发感应电流,Ω3:空气区,为非生热涡流区,材料为非导电介质,无外加激励源,在此区域发生辐射、热对流;
[0128] 2)感应加热过程中涂层熔化的热能来源于基体的感应电流产生的焦耳热,因此,对感应磁场与电场的耦合分析成为感应加热过程理论分析的核心,电磁场分析中,首先引用经典电磁理论的麦克斯韦方程组:
[0129]
[0130]
[0131]
[0132]
[0133] 其中,D为电位移,C/m2,B为磁感应强度,Gs,E为电场强度,N/C,H为磁场强度,Oe,J为电流密度,A/m2,ρ为电荷密度,C/m2,
[0134] 焦耳热产生于自由电荷,电位移密度 不是产生焦耳热的原因,只有在频率达到10MHz时,才会产生焦耳热,再者,因为电流密度J远远大于电位移密度 所以不给计算在内,麦克斯韦方程组转化为:
[0135]
[0136]
[0137]
[0138] 在此过程中的辅助协调公式如下:
[0139] B=μH  (2-9)
[0140] J=σE(欧姆定律)  (2-10)
[0141] 其中,μ为磁导率,σ为电导率,
[0142] 由磁感应强度B满足等式(2-6)的散度为零的,根据亥姆霍兹定理引入磁向量势A,得到等式
[0143]
[0144] 代入公式(2-11)到(2-7)中,可以得到:
[0145]
[0146] 因 满足亥姆霍兹定理,由此引入标量电位势 公式(2-12)可以进一步推导为:
[0147]
[0148] 公式(2-11)与公式(2-13)可以同时满足公式(2-6)与公式(2-7),因此代入辅助协调方程公式(2-9)与公式(2-10)到公式(2-8)中,可得
[0149]
[0150] 标量电位势 与磁向量势A不是唯一的,对于不同的应用规范会有不同的对应值,洛伦兹规范和库仑规范常用于解决电磁协调分析方程,在是准静态近似条件下,相对于洛伦兹规范利用库伦规范可极大的简化模型方程,在库伦规范协应用后,标量势与矢量势都可以被表达为标量电位势 与磁向量势A的函数,因此,库伦规范表达式 被代入公式(2-14),惩罚公式 引入公式(2-14),得到公式
[0151]
[0152] 引入旋度计算公式 到公式(2-15),可得
[0153]
[0154] 代入公式(2-10)与(2-12)可得
[0155]
[0156] 其中,Je为感应电流密度,Js为源电流密度,表达为
[0157] 对于谐波分析,电磁量是单一频率的振动函数,因此A可以写成:
[0158] A=A0ejωt  (2-18)
[0159] Js=J0ejωt  (2-19)
[0160] 其中,J0为源电流密度的振幅,A0为磁势振幅,ω为角频率,可表达为ω=2πf,[0161] 由于构建模型的对称性,外加激励电流可以设计为:
[0162]
[0163]
[0164] 把公式(2-18)与公式(2-19)代入公式公式(2-16),可得
[0165]
[0166] 对于轴对称的问题,柱坐标系下的公式(2-22)转化为
[0167]
[0168] 以上公式适用于不同介质区域的电磁场,包括磁体、非磁体、导电介质和非导电介质,根据感应加热建模区域建模介质不同,由公式(2-13)、(2-15)与(2-17),各个区域描述方程如下:
[0169] 管体区Ω1范围内描述方程:
[0170]
[0171]
[0172]
[0173] 线圈区Ω2范围内描述方程:
[0174]
[0175]
[0176]
[0177] 空气区Ω3范围内描述方程:
[0178] J=0  (2-30)
[0179]
[0180] 针对于不同区域边界介质的特性,设定各划分区域的边界条件如下:
[0181] Ω1法线方向上:
[0182] n·(D2-D1)=σe0  (2-32)
[0183] n·(B2-B1)=0  (2-33)
[0184]
[0185] Ω1径线方向上:
[0186] n×(E2-E1)=0  (2-35)
[0187] n×Hout=i0or n×(H2-H1)=0  (2-36)
[0188] Ω2:法线方向上:
[0189] n·(J2-J1)=0  (2-37)
[0190] Ω3:法线方向上:
[0191] n·(B2-B1)=0  (2-38)
[0192] 径向方向上:
[0193] n×(H2-H1)=0  (2-39)
[0194] 其中,Di为电流密度,Bi为磁流密度,Ji为电流密度,Ei(i=1,2)为感应涡流电场强度,n为单位法线向量,σe0为自由电荷密度,i0为交界表面法相方向上的电流密度,Hout为磁场密度,方向沿着界面指向外侧,特别是在超高频率下公式n×Hout=i0代替公式n×(H2-H1)=0作为边界条件,
[0195] 在感应加热过程中管体被磁化后,其表面应力张量的变化即可产生电磁力,电磁力的计算引用公式Q=0.24I2Rt与公式(2-17),可得磁场强度控制方程:
[0196]
[0197] 为了对方程(2-40)进行求解,需要将方程转化为弱解形式,由此引入引入爱因斯坦张量表达式:
[0198]
[0199] 其中,εijk为符号变量,i,j,k为下标,取得下标组合(1,-1,0),根据公式(2-41),两个变量的叉积即为
[0200] (ν×ω)i=εijkνjωk  (2-42)
[0201] 其中,v与w为任意代数向量,对应的旋度公式为:
[0202]
[0203] 由公式(2-42)和公式(2-43),方程转化成爱因斯坦张量表达式:
[0204]
[0205] 其中,Hk为磁场强度分量,Ai为磁动势分量,A为磁路中的磁动势矢量,σe为磁介质电导率,
[0206] 对公式(2-44)运用加权残值法积分可得:
[0207]
[0208] 其中:ψi为测试函数;VB为磁场的空间积分域。公式(2-45)中包含磁场强度H对位置变量x的散度,以及散度与测试函数ψi的乘积,可以进一步通过分部积分运算性质进行化简,得到
[0209]
[0210] 方程(2-46)中乘积项Hkψi对位置变量x的偏导数的体积积分实际上是乘积Hkψi的散度计算,根据散度计算定理,一个矢量V散度的体积分等于V在此体积表面上的法向分量的面积分,所以方程(2-46)可以进一步表示为
[0211]
[0212] 其中,n表示包围积分空间域,VB表面的法相向量,其指向向外, 表示沿体积VB的表面进行积分。将方程(2-46)代入式(2-47)可得
[0213]
[0214] 公式(2-48)即为爱因斯坦标记下三维磁场分布计算方程。为方便方程的求解,需要将公式(2-48)写成矩阵形式,根据爱因斯坦标记法的运算规则,对方程(2-48)进行矩阵标记的转换:
[0215]
[0216] 式中磁场空间中边界处磁场强度的切向分量满足HT=H×n,对于变分变量,磁感应强度B与磁动势A满足下列表达式:
[0217]
[0218] 则方程(2-50)的弱解形式可以进一步简化为
[0219]
[0220] 公式(2-51)即为弱解形式的电磁场偏微分方程,其中自变量为磁动势矢量A,从公式(2-51)的等号右边可以看到,系统的输入为外加电流密度与系统边界条件的矢量和;
[0221] 3)在研究加热管体材料机械特性保护方面,热应力场分析是辅助热场分析的重要手段,由此,基于热场分析的热-应力耦合分析成为研究感应加热的重要理论依据,感应加热过程涉及管体热能传导的问题,现应用经典热传导的傅里叶定律,得到如下等式:
[0222]
[0223] 其中,ρ为材料密度,kg/m3;C为比热,J/(kg·℃);T为温度时间函数,℃;k热导系数,W/(m·℃);Q为感应涡流产生的热量,J,
[0224] 感应加热的热能由感应涡流产生,结合焦耳定律,可以得到:
[0225]
[0226] 感应加热为非线性瞬态热传导过程问题,在笛卡尔坐标系内,建立三维瞬态热传导控制方程:
[0227]
[0228] 感应加热过程基体不具有内热源,因此qv=0,上式转化为Fourier方程形式:
[0229]
[0230] 对于工件温度的边界条件可以描述为法向导数,根据热力学第一定律表面的对流和热辐射可以表达如下:
[0231]
[0232] 其中,n为边界上指向外侧发的法向向量,Tair为工件的环境温度,℃;h为对流换热系数,W/(m2·℃);εemi为计算辐射的发射率,k为热导率,σb为斯特藩-玻尔兹曼常数或斯特藩常量5.67×10-8W·m-2·K-4,
[0233] 感应加热过程对于不同区域边界上边界表达式各有不同,其具体边界分析如下:
[0234] Ω1管体区与Ω3空气区、Ω1管体区与Ω3空气区在传导热流量为Q时的边界Γ1:
[0235]
[0236] Ω1管体区与Ω3空气区、Ω1管体区与Ω3空气区对流换热为h(T-T0)时的边界Γ2:
[0237]
[0238] 所有感应加热划分区域高温辐射换热为εσ(T4-T04)的边界Γ3:
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈