CN101360983B - 温度推定方法及装置 - Google Patents

温度推定方法及装置 Download PDF

Info

Publication number
CN101360983B
CN101360983B CN2006800512893A CN200680051289A CN101360983B CN 101360983 B CN101360983 B CN 101360983B CN 2006800512893 A CN2006800512893 A CN 2006800512893A CN 200680051289 A CN200680051289 A CN 200680051289A CN 101360983 B CN101360983 B CN 101360983B
Authority
CN
China
Prior art keywords
temperature
value
inferring
mentioned
variation
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN2006800512893A
Other languages
English (en)
Other versions
CN101360983A (zh
Inventor
大野忠久
本桥勇人
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Azbil Corp
Original Assignee
Azbil Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Azbil Corp filed Critical Azbil Corp
Publication of CN101360983A publication Critical patent/CN101360983A/zh
Application granted granted Critical
Publication of CN101360983B publication Critical patent/CN101360983B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01KMEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K7/00Measuring temperature based on the use of electric or magnetic elements directly sensitive to heat ; Power supply therefor, e.g. using thermoelectric elements
    • G01K7/42Circuits effecting compensation of thermal inertia; Circuits for predicting the stationary value of a temperature
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01KMEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K1/00Details of thermometers not specially adapted for particular types of thermometer
    • G01K1/02Means for indicating or recording specially adapted for thermometers
    • G01K1/022Means for indicating or recording specially adapted for thermometers for recording
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01KMEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K17/00Measuring quantity of heat
    • G01K17/02Calorimeters using transport of an indicating substances, e.g. evaporation calorimeters

Landscapes

  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Chemical & Material Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Combustion & Propulsion (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)
  • Control Of Temperature (AREA)
  • Drying Of Semiconductors (AREA)

Abstract

本发明提供一种温度推定方法及装置。温度推定装置(50)具有:存储部,存储关于加热器(3)和工件(1)的热传导模型的参数;温度测量部,测量加热器(3)的温度;工件温度变化量推定部,根据加热器(3)的温度变化量、加热器温度变化量的临时推定值、和工件温度变化量的临时推定值,来推定工件温度变化量的真值,其中,加热器(3)的温度变化量是根据加热器(3)的温度而求得的,加热器温度变化量的临时推定值是基于热传导模型来推定加热器(3)的温度变化量而得到的值,工件温度变化量的临时推定值是基于热传导模型来推定工件(1)的温度变化量而得到的值;以及工件温度推定部,根据所推定的工件温度变化量的真值来推定工件温度的真值。

Description

温度推定方法及装置
技术领域
本发明涉及例如半导体制造装置等的加热冷却处理装置,特别涉及以下温度推定方法及装置,其基于与推定对象之间具有热阻的加热器等温度测定可能点的温度来推定推定对象的温度。
背景技术
在如半导体制造装置那样运行中不能实际测量进行加热冷却处理的被加热冷却物(半导体制造装置时为晶片)的温度的情况下,往往进行如下操作,即、在装置的试运行和调整时预先在被加热冷却物上安装温度传感器,制作关于加热冷却处理部和被加热冷却物的模型,在正式运行时,基于模型边推定被加热冷却物的温度边进行处理。不过,在该方法中,模型和实际的被加热冷却物之间产生偏差时,有产生温度推定误差的问题。
为了解决这个问题,人们考虑在半导体制造装置中应用下述技术,即、使用数值分析法实时进行模型识别(例如,参照日本特许第3331758号公报)。将日本特许第3331758号公报所公开的技术应用在半导体制造装置中时,在不能直接测量想要知道温度的部位(以下称为推定对象)的温度、并且正式运行时安装温度传感器的温度可测量点和推定对象之间的热阻及推定对象的热容量都不清楚的情况下,可以利用递推最小2乘法及梯度法、最大似然估计法来计算求得温度可测量点与推定对象之间的热阻值、推定对象的热容量值,直到根据以热阻值和热容量值为系数的模型而计算出的温度可测量点的温度推定值、和温度可测量点的温度实测值一致,使用所计算出的热阻值来推定推定对象的温度(以下称为已有技术1)。根据这样的已有技术1,即使识别时的模型和实际的推定对象之间存在偏差的情况下,也能精度优良地推定推定对象的温度。
此外,作为使用了其他推定方法的技术,在日本特开2000-88666号公报中公开了推定埋设在地下的管道内电缆的导体温度的技术(以下,称为已有技术2)。在该已有技术2中,测量欲求导体温度的目标管道内部的温度、土壤温度以及所有管道的通电电流值,并基于土壤热阻值求来自周围的管道的热影响,并且根据各管道的通电电流求得导体的热流值,从而求得目标管道附近的土壤温度变化,并根据各深度下土壤原本具有的基底温度和土壤温度变化来求得计算上的土壤温度,根据该土壤温度和目标管道内的导体热流值,按照包含导体温度的解析模型来计算管道内的温度,并比较管道内温度的计算值和管道内温度的实测值,当比较结果为一致时,将解析模型中所使用的导体温度确定为准确的推定结果,当比较结果为不一致时,重新定土壤热阻、基底温度以及计算所使用的热常数,之后再次进行计算。
但是,一边使用数值解析法来识别热阻值一边进行温度推定的已有技术1,存在下述问题,即,在热阻值和热容量值的推定初始值不准确的情况下,在开始准确的温度推定之前需要花费时间,并且在求得系数值之前,推定温度因系数误差而不准确。另外,已有技术1,在热阻和热容量等热常数在加热冷却处理时发生变化的情况下,存在开始准确的温度推定之前需要花费时间的问题。再有,一般来讲,在递推最小2乘法等推定收敛速度快的方法中,存在计算量大的问题。
同样,在已有技术2中,也存在下述问题,即、在土壤的热阻值等热常数的推定初始值不准确的情况下,重新定热常数需要花费时间,在开始准确的温度推定之前需要花费时间。已有技术2为计算管道内电缆的导体温度和土壤的热阻值的例子,温度推定的精度和推定收敛速度也不需要那么高等级,但是,在被加热冷却物的质量因温度的推定精度而发生大的变化的加热冷处理中的温度推定中,提高该精度以及推定收敛速度具有非常重要的意义,即,当不能准确推定被加热冷却物的温度的情况下,就不能进行适当的加热冷却处理。
发明内容
本发明就是为了解决上述问题而做出的,第1目的是提供一种温度推定方法以及装置,即使在不能直接测量推定对象的温度、并且在加热冷却处理时推定对象和温度可测量点之间的热阻发生变化的情况下,也能精度优良且迅速地通过简单的运算处理来推定推定对象的温度。
此外,本发明的第2目的在于提供一种温度推定方法及装置,即使在不能直接推定推定对象的温度、并且在加热冷却处理时推定对象和温度可测量点之间的热阻以及推定对象的热容量发生变化的情况下,也能迅速且精度优良地通过简单的运算处理来推定推定对象的温度。
本发明的温度推定方法,具有:测量步骤,测量与推定对象之间具有热阻的温度可测量点的温度;温度变化量推定步骤,根据上述温度可测量点的温度变化量、温度可测量点的温度变化量临时推定值、和推定对象的温度变化量临时推定值,来推定上述推定对象的温度变化量的真值,其中,上述温度可测量点的温度变化量是根据上述温度可测量点的温度而求得的,温度可测量点的温度变化量临时推定值是基于关于上述温度可测量点和上述推定对象的热传导模型来推定上述温度可测量点的温度变化量而得到的值,推定对象的温度变化量临时推定值是基于上述热传导模型来推定上述推定对象的温度变化量而得到的值;以及温度推定步骤,根据在该温度变化量推定步骤中所推定的温度变化量的真值来推定上述推定对象的温度的真值。
另外,本发明的温度推定装置具有:存储部,存储关于温度可测量点和推定对象的热传导模型的参数,该温度可测量点与上述推定对象之间具有热阻;温度测量部,测量上述温度可测量点的温度;温度变化量推定部,根据上述温度可测量点的温度变化量、温度可测量点的温度变化量临时推定值、和推定对象的温度变化量临时推定值,来推定上述推定对象的温度变化量的真值,其中,上述温度可测量点的温度变化量是根据上述温度可测量点的温度而求得的,温度可测量点的温度变化量临时推定值是基于上述热传导模型来推定上述温度可测量点的温度变化量而得到的值,推定对象的温度变化量临时推定值是基于上述热传导模型来推定上述推定对象的温度变化量而得到的值;以及温度推定部,根据由该温度变化量推定部所推定的温度变化量的真值来推定上述推定对象的温度的真值。
此外,本发明的温度推定方法,具有:测量步骤,测量与推定对象之间具有热阻的温度可测量点的温度;温度变化量推定步骤,根据上述温度可测量点的温度变化量的实测值、上述温度可测量点的温度变化量临时推定值、和上述推定对象的温度变化量第1推定值,来求取上述推定对象的温度变化量第2临时推定值,上述推定对象的温度变化量第2临时推定值是视上述推定对象的热容量值未发生变化时的值,其中,上述温度可测量点的温度变化量的实测值是根据上述温度可测量点的温度而求得的,上述温度可测量点的温度变化量临时推定值是基于关于上述温度可测量点和上述推定对象的热传导模型进行推定而得到的值;上述推定对象的温度变化量第1推定值是基于上述热传导模型进行推定而得到的值;热阻值变化系数推定步骤,根据上述推定对象的温度变化量第2临时推定值和上述热传导模型来推定上述热阻值的变化值;以及温度推定步骤,根据不同采样时刻的上述温度可测量点的温度、和上述热阻值的变化值来推定上述推定对象的温度的真值。
另外,本发明的温度推定装置,具有:存储部,存储关于温度可测量点和推定对象的热传导模型的参数,该温度可测量点与上述推定对象之间具有热阻;温度测量部,测量上述温度可测量点的温度;温度变化量推定部,根据上述温度可测量点的温度变化量的实测值、上述温度可测量点的温度变化量临时推定值、和上述推定对象的温度变化量第1临时推定值,来求取上述推定对象的温度变化量第2临时推定值,上述推定对象的温度变化量第2临时推定值是视上述推定对象的热容量值未发生变化时的值,其中,上述温度可测量点的温度变化量的实测值是根据上述温度可测量点的温度而求得的;上述温度可测量点的温度变化量临时推定值是基于上述热传导模型进行推定而得到的值,上述推定对象的温度变化量第1临时推定值是基于上述热传导模型进行推定而得到的值;热阻值变化系数推定部,根据上述推定对象的温度变化量第2临时推定值和上述热传导模型来推定上述热阻值的变化值;以及温度推定部,根据不同采样时刻的上述温度可测量点的温度、和上述热阻值的变化值来推定上述推定对象的温度的真值。
如以上所说明的那样,根据本发明,通过执行下述步骤,即使在不能直接测量推定对象的温度、并且推定对象和温度可测量点之间的热阻随每个推定对象而产生偏差,或者在加热冷处理中热阻从识别热传导模型时的值发生了变化的情况下,也能精度优良且迅速地推定推定对象的温度,上述步骤是指,测量步骤,测量温度可测量点的温度;温度变化量推定步骤,根据温度可测量点的温度变化量、温度可测量点的温度变化量临时推定值、和推定对象的温度变化量临时推定值,来推定推定对象的温度变化量的真值,其中,温度可测量点的温度变化量是根据温度可测量点的温度而求得的,温度可测量点的温度变化量临时推定值是基于关于温度可测量点和推定对象的热传导模型来推定温度可测量点的温度变化量而得到的值,推定对象的温度变化量临时推定值是基于热传导模型来推定推定对象的温度变化量而得到的值;以及温度推定步骤,根据在该温度变化量推定步骤中推定的温度变化量的真值来推定推定对象的温度的真值。此外,在本发明中,由于即使热阻在温度推定中持续变化,也不受热阻变化的影响,所以能够迅速且精度优良地推定推定对象的温度。另外,在本发明中,由于不需要数值分析处理和统计处理,所以能够将计算量减少到非常少,即使在处理能力比较弱的加热冷却处理装置的温度控制处理器中,也能以足以进行实时处理的程度的简单运算来进行推定处理。
此外,根据本发明,通过执行下述步骤,即使在不能直接测量推定对象的温度、并且推定对象与温度可测量点之间的热阻值和推定对象的热容量值不明确、热阻值和热容量值随每个推定对象而有偏差,或者在加热冷却处理中热阻值和热容量值从识别热传导模型时的值发生了变化的情况下,也能精度优良且迅速地推定推定对象的温度,上述步骤是,测量步骤,测量温度可测量点的温度;温度变化量推定步骤,根据温度可测量点的温度变化量的实测值、温度可测量点的温度变化量临时推定值、和推定对象的温度变化量第1推定值,来求取推定对象的温度变化量第2临时推定值,该推定对象的温度变化量第2临时推定值是视推定对象的热容量值未发生变化时的值,其中,温度可测量点的温度变化量的实测值是根据温度可测量点的温度而求得的,温度可测量点的温度变化量临时推定值是基于关于温度可测量点和推定对象的热传导模型进行推定而得到的值,推定对象的温度变化量第1推定值是基于热传导模型进行推定而得到的值;热阻值变化系数推定步骤,根据推定对象的温度变化量第2临时推定值和热传导模型来推定热阻值的变化值;以及温度推定步骤,根据不同采样时刻的温度可测量点的温度、和热阻值的变化值来推定推定对象的温度的真值。此外,在本发明中,由于即使在温度推定中推定对象和温度可测量点之间的热阻值和推定对象的热容量值持续变化,只要以与它们的变化相比足够快的周期执行运算,既可不受热阻值的变化和热容量值的变化的影响,所以能够迅速且精度优良地推定推定对象的温度。此外,在本发明中,能够迅速且精度优良地计算推定对象与温度可测量点之间的热阻值的变化。再有,在本发明中,只沿时间轴方向进行少数几个采样的数值计算,不需要进行数值解析处理和统计处理,所以能够将计算量减少到非常少,即使在处理能力比较弱的加热冷却处理装置的温度控制处理器中,也能以足以进行实时处理的程度的简单运算,来进行温度推定处理以及热阻值变化系数推定处理。
附图说明
图1是表示本发明的实施例1的加热冷却处理装置的结构的框图。
图2是表示图1的加热冷却处理装置的控制器的温度推定装置的结构例的框图。
图3是表示图2的温度推定装置的动作的流程图。
图4是表示利用图2的温度推定装置得到的工件温度的推定结果和工件温度的实测值的1个例子的图。
图5是表示利用图2的温度推定装置得到的热阻值变化系数的推定结果和热容量值变化系数的实测值的1个例子的图。
图6是表示利用图2的温度推定装置得到的工件温度的推定结果和工件温度的实测值的其他例子的图。
图7是表示利用图2的温度推定装置得到的热阻值变化系数的推定结果和热阻值变化系数的实测值的其他例子的图。
图8是表示利用图2的温度推定装置得到的工件温度的推定结果和工件温度的实测值的其他例子的图。
图9是表示利用图2的温度推定装置得到的热阻值变化系数的推定结果和热阻值变化系数的实测值的其他例子的图。
图10是表示本发明的实施例2的加热冷却处理装置中的控制器的温度推定装置的结构例的框图。
图11是表示图10的温度推定装置的动作的流程图。
图12是表示利用图10的温度推定装置得到的工件温度的推定结果和工件温度的实测值的1个例子的图。
图13是表示利用图10的温度推定装置得到的工件温度的推定结果和工件温度的实测值的其他例子的图。
图14是表示利用图10的温度推定装置得到的热阻值变化系数的推定结果和热阻值变化系数的实测值的1个例子的图。
图15是表示利用图10的温度推定装置得到的热容量变化系数的推定结果和热容量值变化系数的实测值的1个例子的图。
图16是表示利用图10的温度推定装置得到的工件温度的推定结果和工件温度的实测值的其他例子的图。
图17是表示利用图10的温度推定装置得到的热阻值变化系数的推定结果和热阻值变化系数的实测值的其他例子的图。
图18是表示利用图10的温度推定装置得到的热容量变化系数的推定结果和热容量变化系数的实测值的其他例子的图。
具体实施方式
(实施例1)
图1表示了本发明的实施例1的加热冷却处理装置的结构。在图1中,2是承载工件1的加热板,3是安装在加热板2上的加热器,4是安装在加热器3上的温度传感器,5是控制器,6是环境温度传感器。控制器5具有温度推定装置50和控制装置51。控制器5可以由具有CPU、存储装置及接口的计算机、和控制这些硬件资源的程序来实现。CPU按照存储装置所存储的程序执行以下处理。
以下,对温度推定装置50的温度推定原理进行说明。设加热器3的温度为Th、要推定的工件1的温度为Tw、周围的环境温度为Te、从控制器5的控制装置51向加热器3输出的控制器输出值为MV、采样时间(采样周期)为ts,关于加热器温度Th和工件温度Tw,根据关于温度可测量点(加热器3)和推定对象(工件1)的热传导模型,建立式子,进行离散化并整理后,如下所示。
Tw(n+1)=Tw(n)+ts·(a·(Th(n)-Tw(n))-b·(Tw(n)-Te(n)))...(1)
Th(n+1)=Th(n)+ts·(c·MV(n)-d·(Th(n)-Tw(n)))...(2)
在式(1)、式(2)中,Tw(n)为时刻n时的工件温度、Tw(n+1)为从时刻n经采样时间ts后的时刻n+1时的工件温度、Th(n)为时刻n时的加热器温度、Th(n+1)为时刻n+1时的加热器温度。此外,系数a、b、c、d为热传导模型的参数。这些系数a、b、c、d,是在装置的维护(maintenance)中,以时间序列计量加热器温度Th、工件温度Tw以及环境温度Te,通过使用最小2乘法而识别的值。并且加热器温度Th和环境温度Te在维护时和正式运行时的任意一者中都能测量,但是工件温度Tw仅能在维护时测量,在正式运行中不能测量。在这里,工件初始温度可以测量或者可以利用其他位置的温度近似。
此外,设在正式运行时,作为温度可测量点的加热器3和作为推定对象的工件1之间的热阻值,随工件的不同而存在偏差,是不明确的。设正式运行中相对于识别热传导模型时求得的热阻值发生了变化的热阻的变化值(准确来说,为热阻的倒数的变化值)为热阻值变化系数mr,对基于正式运行时的热传导模型的式子进行离散化,整理后的结果如下所示,该热传导模型考虑了热阻值变化系数mr
Tw(n+1)=Tw(n)+ts·(mr·a·(Th(n)-Tw(n))-b·(Tw(n)-Te(n)))...(3)
Th(n+1)=Th(n)+ts·(c·MV(n)-mr·d·(Th(n)-Tw(n)))...(4)
式(1)、式(2)为基于识别时的热传导模型来推定工件温度Tw、加热器温度Th的式子,而如上所述,式(3)、式(4)是进而考虑热阻值变化系数mr来推定工件温度Tw、加热器温度Th的式子,利用式(3)、式(4)而求得的工件温度Tw(n+1)、加热器温度Th(n+1)相当于真值。
当时刻n之前的加热器温度Th、工件温度Tw以及环境温度为Te已知,或者可以适当地推定时,对于从时刻n至时刻n+1的2个采样之间的工件温度变化量,基于识别时的热传导模型来进行推定时,则由式(1)变为下式。
Tw(n+1)-Tw(n)=ts·(a·(Th(n)-Tw(n))-b·(Tw(n)-Te(n)))...(5)
将式(5)所示的值称为工件温度变化量的临时推定值。
此外,在考虑正式运行时的每个工件的热阻值变化系数mr时,工件温度变化量由式(3)变成下式。
Tw(n+1)-Tw(n)=ts·(mr·a·(Th(n)-Tw(n))-b·(Tw(n)-Te(n)))...(6)
将式(6)所示的值称为工件温度变化量的真值。
两个采样之间的工件温度变化量的临时推定值和真值之间的差,为式(5)的右边与式(6)的右边之差,如下式所示。
ts·(1-mr)·a·(Th(n)-Tw(n))...(7)
同样,当时刻n之前的加热器温度Th、工件温度Tw以及环境温度Te已知,或者可以适当地推定时,对于从时刻n至时刻n+1的2个采样之间的加热器温度变化量,基于识别时的热传导模型进行推定时,由式(2)变为下式。
Th(n+1)-Th(n)=ts·(c·MV(n)-d·(Th(n)-Tw(n)))...(8)
将式(8)所示的值称为加热器温度变化量的临时推定值。
此外,当考虑正式运行时的每个工件的热阻值的变化系数mr时,加热器温度变化量由式(4)变为下式。
Th(n+1)-Th(n)=ts·(c·MV(n)-mr·d·(th(n)-Tw(n)))...(9)
将式(9)所示的值称为加热器温度变化量的真值。
2个采样之间的加热器温度变化量的临时推定值和真值之间的差,为式(8)的右边与式(9)的右边之差,如下式所示。
ts·(1-mr)·d·(Th(n)-Tw(n))...(10)
在这里,认真观察式(7)、式(10),可知有下式成立。
(工件温度变化量的临时推定值-工件温度变化量的真值)/(加热器温度变化量的真值-加热器温度变化量的临时推定值)=a/d...(11)
即,根据实际的测量值求2个采样之间的加热器温度变化量的真值,计算2个采样之间的加热器温度变化量的真值和临时推定值之间的差,在所计算出的差值上乘以a/d,由此能求得2个采样之间的工件温度变化量的临时推定值和真值之间的差。
基于识别时的热传导模型的2个采样之间的工件温度变化量的临时推定值是可以计算的。当设2个采样之间的加热器温度变化量的真值为ΔTh、2个采样之间的加热器温度变化量的临时推定值为带帽ΔTh,model(以下,同样将在字符上标记的“∧”称为带帽)、2个采样之间的工件温度变化量的临时推定值为带帽ΔTw,model时,2个采样之间的工件温度变化量的真值的推定值带帽ΔTw,可以如下式这样求得。
(算式1)
Δ T ^ w = Δ T ^ w , mode l - ( a d ) · ( ΔT h - Δ T ^ h , mode l ) · · · ( 12 )
利用式(12)的关系,在设自时刻n至时刻n+1的2个采样之间的工件温度变化量的真值的推定值为带帽ΔTw(n)时,时刻n+1时的工件温度的推定值带帽Tw(n+1),可以根据时刻n时的工件温度的推定值带帽Tw(n)如下式这样计算出来。
(算式2)
T ^ w ( n + 1 ) = Δ T ^ w ( n ) + T ^ w ( n ) · · · ( 13 )
此外,在式(6)中,使用工件温度的推定值带帽Tw(n)替代时刻n时的工件温度Tw(n)、使用工件温度变化量的推定值带帽ΔTw(n)替代工件温度变化量Tw(n+1)-Tw(n)时,可以利用下式求得时刻n时的热阻值变化系数mr的推定值带帽mr(n)。
(算式3)
m ^ r ( n ) = 1 ( a · ( T h ( n ) - T ^ w ( n ) ) ) · ( Δ T ^ w ( n ) t s + b · ( T ^ w ( n ) - T e ( n ) ) ) · · · ( 14 )
如图2所示,温度推定装置50具有:存储部500、温度测量部501、工件温度变化量推定部502、工件温度推定部503、以及热阻值变化系数推定部504,存储部500存储预先识别的热传导模型的参数、测量数据和计算数据。
使用图3,对温度推定装置50的动作进行说明。首先,温度测量部501在正式运行前的初始状态下,取得由温度传感器4测量的加热器温度Th、由环境温度传感器6测量的环境温度Te、和由仅可在初始状态使用的温度传感器测量的工件温度Tw。存储部500存储温度测量部501所取得的加热器温度Th、环境温度Te、工件温度Tw(图3步骤S1)。
接下来,开始正式运行,温度测量部501取得当前时刻n+1时的加热器温度Th(n+1)和环境温度Te(n+1)。存储部500存储温度测量部501所取得的加热器温度Th(n+1)和环境温度Te(n+1)(步骤S2)。
然后,工件温度变化量推定部502参照存储部500所存储的值,计算当前时刻n+1的加热器温度Th(n+1)和1个采样前的时刻n的加热器温度Th(n)之间的差,即、计算2个采样之间的加热器温度变化量的真值ΔTh,利用式(8)计算2个采样之间的加热器温度变化量的临时推定值带帽ΔTh,model,利用式(5)计算2个采样之间的工件温度变化量的临时推定值带帽ΔTw,model,根据这些计算出的值,利用式(12)来计算从时刻n至时刻n+1的2个采样之间的工件温度变化量的推定值带帽ΔTw(n)(步骤S3)。存储部500存储该工件温度变化量的推定值带帽ΔTw(n)。
此时,在当前时刻n+1为工件投入后的最初采样时刻的情况下,式(5)、式(8)中的加热器温度Th(n)是在初始状态下测量的温度,工件温度Tw(n)是在初始状态下测量的温度或利用其他位置的温度近似的温度。此外,式(8)中的控制器输出值MV(n)是在时刻n由控制装置51所计算出的值,但是在当前时刻n+1为工件投入后的最初采样时刻的情况下,为初始状态下所输出的值。
接下来,工件温度推定部503,根据存储部500所存储的工件温度变化量的推定值带帽ΔTw(n)和工件温度的推定值带帽Tw(n),利用式(13)计算当前时刻n+1时的工件温度的推定值带帽Tw(n+1)(步骤S4)。存储部500存储该工件温度的推定值带帽Tw(n+1)。并且,在当前时刻n+1为工件投入后的最初采样时刻时,作为工件温度的推定值带帽Tw(n)使用在初始状态下测量的、或者利用其他位置的温度近似的工件温度Tw
然后,热阻值变化系数推定部504,根据存储部500所存储的工件温度变化量的推定值带帽ΔTw(n)、工件温度的推定值带帽Tw(n)、加热器温度Th(n)和环境温度Te(n),利用式(14)计算时刻n时的热阻值变化系数mr的推定值带帽mr(n)(步骤S5)。存储部500存储该热阻值变化系数mr的推定值带帽mr(n)。另外,在当前时刻n+1为工件投入后的最初采样时刻时,环境温度Te(n)为初始状态下所测量的值。
以下,同样地在加热冷却处理装置的动作停止之前(在步骤S6中为是),在每个采样时间反复进行步骤S2~S5的处理。在当前时刻n+1不是工件投入后的最初采样时刻的情况下,加热器温度Th(n)和环境温度Te(n)是1个采样前所取得的值,作为式(5)、式(8)的工件温度Tw(n),使用工件温度的推定值带帽Tw(n)。
控制器5的控制装置51,与这样的温度推定装置50的动作并行地取得温度推定装置50所计算出的时刻n时的工件温度的推定值,计算出控制器的输出值MV,使得该工件温度的推定值和目标温度一致,在时刻n+1时向加热器3输出该控制器输出值MV。
如以上那样,在本实施例中,即使在正式运行中工件1和加热器3之间的热阻值从模型识别时的值发生了变化,也能即时推定工件温度,而不用知道其变化量。此外,在本实施例中,可以利用所计算出的工件温度来即时推定热阻值的变化量。此外,在本实施例中,即使不知道热阻值的变化量,也能迅速地进行温度推定,所以,即使在温度推定中热阻值连续变化的情况下也能即时进行温度推定,也可以利用它来计算连续变化的热阻值。此外,在本实施例中,由于在温度的推定运算中未使用数值解析和统计运算,所以计算负担轻,从而也可以安装处理能力弱的处理器。
此外,在本实施例中,通过执行热阻值变化系数推定步骤,即、根据工件温度的真值和热传导模型来推定工件1和加热器3之间的热阻的变化值,能够迅速且精度优良地计算热阻值的变化。在这种情况下,热阻值的变化,即使是非预期的变化也没有关系。再有,在本实施例中,将此前所求得的热阻值用作固定值,加热器3的温度变小到接近测量误差程度的水平之后,也能精度优良地求得工件温度。
图4~图9表示本实施例的效果。图4、图6、图8是表示本实施例的工件温度的推定结果和工件温度的实测值的1个例子的图。图5、图7、图9是表示本实施例的热阻值变化系数的推定结果和热阻值变化系数的实测值的1个例子的图。在图4、图5的例子中,假定了工件1和加热器3之间的热阻值在正式运行时从模型识别时的值发生了变化,并达到一定值的情况,工件温度的推定值带帽Tw和工件温度的实测值Tw充分一致,同样,热阻值变化系数的推定值带帽mr和热阻值变化系数的实测值mr也充分一致,从而可知能够精度优良地推定工件温度和热阻值。此外,在图6、图7的例子中,假定了工件1和加热器3之间的热阻值在温度推定中持续变化的情况,可知,能够精度优良地推定工件温度、热阻值。此外,在图8、图9的例子中,假定了工件1和加热器3之间的热阻值在温度推定中不规则地持续变化的情况,可知,可以精度优良地推定工件温度、热阻值。
在本实施例中,加热器温度变化量的实测值与基于热传导模型的加热器温度变化量的临时推定值之间的差,和基于热传导模型的工件温度变化量的临时推定值与工件温度变化量的真值之间的差,呈以a/d表示的固定关系。该比值a/d,若返回到原热传导模型中的微分方程式来考虑,则为由加热器温度所代表的温度可测量点与由工件温度所代表的推定对象的热容量比。
此外,在本实施例中,以工件为推定对象,并设加热器为温度可测量点,不过,反之即使在工件温度可测量、加热器温度不可测量的情况下,式(11)的关系也完全相同,所以,可以以ΔTw为工件温度变化量的实测值,如下式所示这样计算出加热器温度变化量的真值的推定值带帽ΔTh
(算式4)
Δ T ^ h = Δ T ^ h , mode l - ( d a ) · ( Δ T w - Δ T ^ w , mode l ) · · · ( 15 )
感应加热装置等,由于在工件侧产生热量,所以可以将工件看做加热源,并且,很多情况下不可测量工件温度。在这种情况下,设工件温度为Th、装置可测量点的温度为Tw,使用式(15)来推定工件温度的方法是适当的。
此外,在本实施方式中,温度测量要素仅为加热器温度、工件温度、环境温度,在正式运行中模型从模型识别时发生的变化只是加热器与工件之间的热阻值的情况下,即使增加加热器和工件间以外的模型要素的数量,也仅仅是利用式(5)、式(6)、式(8)、式(9)所消去的项增加,其结果是,能够导出式(11),所以可以应用相同的方法。即使在如上所述工件温度可测量、加热器温度不可测量的情况下,在能够如感应加热装置等那样将工件看做加热源、并且工件温度不可测量的情况下,或者在控制器的输出与加热器之间引入其他模型要素的情况下,增加加热器和工件间以外的模型要素的数量,也同样能导出式(11)。
再有,在本实施方式中,加热器温度,被模型化为直接受控制器输出值的影响,但是,如下述的式(16)、式(17)所示,也可以在控制器的输出和加热器之间引入其他的模型要素。在式(16)、式(17)中,设模型要素p的温度为Tp
Th(n+1)-Th(n)=ts·(c·(Tp(n)-Th(n))-d·(Th(n)-Tw(n)))...(16)
Th(n+1)-Th(n)=ts·(c·(Tp(n)-Th(n))-mr·d·(Th(n)-Tw(n)))..(17)
替代式(8)、式(9),求式(16)的右边与式(17)的右边之间的差,消去加热器温度和工件温度以外的要素,结果可计算出式(10),其结果是也可导出式(11)。
此外,在本实施例中,在计算热阻值变化系数mr的推定值带帽mr(n)时,由于使用将式(6)变形后的式(4),所以使用工件温度的推定值带帽Tw(n)和工件温度变化量的推定值带帽ΔTw(n),不过,热阻值变化系数mr的推定值带帽mr(n),也可以利用工件温度的推定值带帽Tw(n)和将式(9)变形后的下述式(18)来求得。
(算式5)
m ^ r ( n ) = 1 ( d · ( T h ( n ) - T ^ w ( n ) ) ) · ( c · MV ( n ) - Δ T h ( n ) t s ) · · · ( 18 )
同样,在控制器的输出与加热器之间引入其他模型要素的情况下,热阻值变化系数mr的推定值带帽mr(n),也可以利用将式(17)变形后的下述式(19)来求得。
(算式6)
m ^ r ( n ) = 1 ( d · ( T h ( n ) - T ^ w ( n ) ) ) · ( c · ( T p ( n ) - T h ( n ) ) - Δ T h ( n ) t s ) · · · ( 19 )
同样,在可以如感应加热装置等那样将工件看做加热源、并且不能测量工件温度的情况下,热阻值变化系数mr的推定值带帽mr(n),可以利用以带帽Tw置换式(14)或式(18)的Th、以Th置换带帽Tw后的式子来求得。
(实施例2)
接下来,对本发明的实施例2进行说明。在本实施例中,由于加热冷却处理装置的构成和实施例1相同,所以使用图1的符号进行说明。
首先,对本实施例的温度推定装置50的温度推定原理进行说明。当针对加热器温度Th和工件温度Tw,根据关于温度可测量点(加热器3)和推定对象(工件1)的热传导模型来建立式子,进行离散化并整理后,时刻n+1时的工件温度Tw(n+1)变为上述式(1),时刻n+1时的加热器温度Th(n+1)变成上述式(2)。
设在正式运行时,作为推定对象的工件1的热容量值,对于每个工件,其偏差都不明。设正式运行中相对于热传导模型识别时所求得的热容量值发生了变化的热容量的变化值为热容量值变化系数kr,对基于正式运行时的热传导模型的式子进行离散化,整理后的结果如下所示,该热传导模型考虑了热容量值变化系数kr
Tw(n+1)=Tw(n)+(ts/kr)·(a·(Th(n)-Tw(n))-b·(Tw(n)-Te(n)))...(20)
Th(n+1)=Th(n)+ts·(c·MV(n)-d·(Th(n)-Tw(n)))...(21)
观察式(21)可知,工件1的热容量的变化没有直接呈现在关于加热器温度Th的模型式中。关于工件温度Tw,设时刻n之前的温度已知或可以适当地推定时,关于从时刻n到时刻n+1的加热器温度Th变化的式子成为下式。
(算式7)
T h ( n + 1 ) = T h ( n ) + t S · ( c · MV ( n ) - d · ( T h ( n ) - T ^ w ( n ) ) ) · · · ( 22 )
同样,关于工件温度Tw,设时刻n之前的温度已知或可以适当地推定时,关于从时刻n+1到时刻n+2的加热器温度Th变化的式子成为下式。
(算式8)
T h ( n + 2 ) = T h ( n + 1 ) + t S · ( c · MV ( n + 1 ) - d · ( T h ( n + 1 ) - T ^ w ( n + 1 ) ) ) · · · ( 23 )
式(22)、式(23)中不可测量或者未知的参数,为时刻n时的加热器温度的推定值带帽Tw(n)(以下,同样将在字符上附加的“∧”称为带帽)、时刻n+1时的工件温度的推定值带帽Tw(n+1)。分别使用1个时间点前的加热器温度来推定带帽Tw(n)、带帽Tw(n+1)。即,设使用时刻n+1的加热器温度Th(n+1)而推定的、时刻n时的工件温度的推定值为带帽Tw,rev(n)时,带帽Tw,rev(n)可以根据式(22)如下式所示这样进行推定。
(算式9)
T ^ w , rev ( n ) = ( 1 d ) · ( ( 1 t S ) · ( T h ( n + 1 ) - T h ( n ) ) - c · MV ( n ) ) + T h ( n ) · · · ( 24 )
同样,设使用时刻n+2的加热器温度Th(n+2)而推定的、时刻n+1时的工件温度的推定值为带帽Tw,rev(n+1)时,带帽Tw,rev(n+1)可以根据式(23)如下式所示这样进行推定。
(算式10)
T ^ w , rev ( n + 1 ) = ( 1 d ) · ( ( 1 t S ) · ( T h ( n + 2 ) - T h ( n + 1 ) ) - c · MV ( n + 1 ) ) + T h ( n + 1 ) · · · ( 25 )
此外,可以根据式(20)而得到以下的式(26)。
(算式11)
k ^ r ( n ) · ( T ^ w ( n + 1 ) - T ^ w ( n ) ) = t S · ( a · ( T h ( n ) - T ^ w ( n ) ) - b · ( T ^ w ( n ) - T e ( n ) ) ) · · · ( 26 )
这里,在根据式(24)、式(25)计量了时刻n+2的加热器温度Th(n+2)时,式(26)的左边的带帽Tw(n+1)-带帽T(n)部分,可以推定为带帽Tw,rev(n+1)-带帽Tw,rev(n)。式(26)的右边,带帽Tw(n)在时刻n+1时可以推定为带帽Tw,rev(n),所以,结果是在n+2的时间点,能够计算出时刻n时的热容量值变化系数的推定值带帽kr(n)。
由以上可知,即使在正式运行中热容量值相对于热传导模型识别时所求得的工件1的热容量值的变化程度不明的情况下,也能计算工件温度和工件热容量值。此外,只要模型式的1个步骤与工件热容量值的变化相比,为足够短的间隔,则即使热容量值动态地变化,也可以进行同样的计算。
接下来,在正式运行时,在作为温度可测量点的加热器3和作为推定对象的工件1之间的热阻值、以及对于每个工件而言工件1的热容量值的偏差不明的情况下,对工件温度、工件热容量值、加热器3与工件1之间的热阻值进行推定。设正式运行中相对于热传导模型识别时所求得的热阻值发生变化了的热阻值的变化值(准确来说,是热阻值的倒数的变化值)为热阻值变化系数mr,如前所述,设正式运行中相对于热传导模型识别时所求得的热容量值发生了变化的热容量值的变化值为热容量值变化系数kr,对基于正式运行时的热传导模型的式子进行离散化,整理后的结果如下所示,该热传导模型考虑了这些变化系数。
Tw(n+1)=Tw(n)+(ts/kr)·(mr·a·(Th(n)-Tw(n))-b·(Tw(n)-Te(n)))...(27)
Th(n+1)=Th(n)+ts·(c·MV(n)-mr·d·(Th(n)-Tw(n)))...(28)
时刻n到时刻n+1的2个采样之间的工件温度变化量,基于识别时的热传导模型进行推定时,如下式所示。
Tw(n+1)-Tw(n)=ts·(1·a·(Th(n)-Tw(n))-b·(Tw(n)-Te(n)))...(29)
将式(29)所示的工件温度变化量Tw(n+1)-Tw(n)称为工件温度变化量的第1临时推定值。
此外,在考虑正式运行时的每个工件的热阻值及热容量值与模型识别时的偏差时,2个采样之间的工件温度变化量为下式所示。
kr·(Tw(n+1)-Tw(n))=ts·(mr·a·(Th(n)-Tw(n))-b·(Tw(n)-Te(n)))...(30)
将式(30)所示的工件温度变化量kr·(Tw(n+1)-Tw(n))称为工件温度变化量的第2临时推定值。该第2临时推定值,是工件温度变化量的真值(Tw(n+1)-Tw(n))放大kr倍后的值,是视工件1的热容量值未发生变化(即,设热容量值变化系数kr为1)时的值。此外,该第2临时推定值,也可以考虑为加热器3和工件1之间的热阻值发生变化的值。
2个采样之间的工件温度变化量的第1临时推定值和第2临时推定值之间的差,为式(29)的右边与式(30)右边之间的差,如下式所示。
ts·(1-mr)·a·(Th(n)-Tw(n))...(31)
同样,时刻n到时刻n+1的2个采样之间的加热器温度变化量,在基于识别时的热传导模型进行推定时,如下式所示。
Th(n+1)-Th(n)=ts·(c·MV(n)-1·d·(Th(n)-Tw(n)))..(32)
将式(32)所示的加热器温度变化量Th(n+1)-Th(n)称为加热器温度变化量的临时推定值。
此外,考虑热阻值的偏差时,加热器温度变化量如下式所示。
Th(n+1)-Th(n)=ts·(c·MV(n)-mr·d·(Th(n)-Tw(n)))..(33)
将式(33)所示的加热器温度变化量Th(n+1)-Th(n)称为加热器温度变化量的真值。
2个采样之间的加热器温度变化量的真值与临时推定值之间的差,为式(32)的右边与式(33)的右边之间的差,如下式所示。
ts·(1-mr)·d·(Th(n)-Tw(n))...(34)
在这里,认真观察式(31)、式(34),则可知下式成立。
(工件温度变化量的第1临时推定值-工件温度变化量的第2临时推定值)/(加热器温度变化量的真值-加热器温度变化量的临时推定值)=a/d
…(35)
即,根据实际的测量值来求2个采样之间的加热器温度变化量的真值,计算2个采样之间的加热器温度变化量的真值与临时推定值之间的差,通过在所计算出的值上乘以a/d,可以求出2个抽样之间的工件温度变化量的第1临时推定值与第2临时推定值之间的差。
基于识别时的热传导模型的2个采样之间的工件温度变化量的第1临时推定值是可以计算的。在设2个采样之间的加热器温度变化量的真值为ΔTh(n)、2个采样之间的加热器温度变化量的临时推定值为带帽ΔTh,model(n)、2个采样之间的工件温度变化量的第1临时推定值为带帽ΔTw,model(n)、2个采样之间的工件温度变化量的第2临时推定值为带帽ΔTw,tmp(n)时,如下式所示这样求取带帽ΔTw,tmp(n)。
(算式12)
Δ T ^ w , tmp ( n ) = Δ T ^ w , mode l ( n ) - ( a d ) · ( Δ T h ( n ) - Δ T ^ h , mode l ( n ) ) · · · ( 36 )
时刻n+1时的工件温度的推定值带帽Tw,tmp(n+1),可以根据时刻n时的工件温度推定值带帽Tw(n),如下式所示这样计算得出。
(算式13)
T ^ w , tmp ( n + 1 ) = Δ T ^ w , tmp ( n ) + T ^ w ( n ) · · · ( 37 )
另一方面,在式(30)中,替代时刻n时的工件温度Tw(n)而使用工件温度的推定值带帽Tw(n),作为kr·(Tw(n+1)-Tw(n))=带帽ΔTw,tmp(n),针对mr求解时,时刻n时的热阻值变化系数mr的推定值带帽mr(n)可以利用下式求得。
(算式14)
m ^ r ( n ) = 1 ( a · ( T h ( n ) - T ^ w ( n ) ) ) · ( Δ T ^ w , tmp ( n ) t s + b · ( T ^ w ( n ) - T e ( n ) ) ) · · · ( 38 )
式(38),表示只要时刻n的工件温度推定值带帽Tw(n)似乎可靠,则即使不知道热容量值变化系数kr,也能推定出热阻值变化系数mr
接下来,使用式(33),根据时刻n至时刻n+1的2个采样之间的加热器温度变化量的真值ΔTh(n)=Th(n+1)-Th(n)、和利用式(38)求得的时刻n时的热阻值变化系数mr的推定值带帽mr(n),反向求出时刻n时的工件温度的推定值带帽Tw,rev(n)。
(算式15)
T ^ w , rev ( n ) = ( Δ T h ( n ) t S - c · MV ( n ) ) ( m ^ r ( n ) · d ) + T h ( n ) · · · ( 39 )
接下来,在下一个采样时刻的计量步骤中,同样地使用时刻n+1到时刻n+2的2个采样之间的加热器温度变化量的真值ΔTh(n+1)=Th(n+2)-Th(n+1),反向求出时刻n+1时的加热器温度的推定值带帽Tw,rev(n+1)。
(算式16)
T ^ w , rev ( n + 1 ) = ( Δ T h ( n + 1 ) t S - c · MV ( n + 1 ) ) ( m ^ r ( n ) · d ) + T h ( n + 1 ) · · · ( 40 )
当温度推定的计算周期与热阻值的变化速度相比足够快的情况下,由于时刻n+1时的热阻值变化系数mr的推定值带帽mr(n+1)、和时刻n时的推定值带帽mr(n)为基本相同的值,所以,在式(40)中,替代带帽mr(n+1)而使用时刻n时的热阻值变化系数mr的推定值带帽mr(n)。即,在本实施例中,是设热传导模型式在轮回2个采样周期左右的短暂时间内热阻值不变化,而推定工件温度的。
如式(41)所示那样,求取式(40)所示的时刻n+1时的工件温度的推定值带帽Tw,rev(n+1)、与式(39)所示的时刻n时的工件温度的推定值带帽Tw(n)之差,作为时刻n至时刻n+1的2个采样之间的工件温度变化量的第3临时推定值带帽ΔTw,rev(n),。
(算式17)
Δ T ^ w , rev ( n ) = T ^ w , rev ( n + 1 ) - T ^ w ( n ) · · · ( 41 )
利用式(36)所求得的工件温度变化量的第2临时推定值带帽ΔTw, tmp(n),为带帽ΔTw,tmp(n)=kr·(Tw(n+1)-Tw(n)),以带帽ΔTw, rev(n)替代该(Tw(n+1)-Tw(n)),针对kr进行求解时,时刻n时的热容量值变化系数kr的推定值带帽kr(n)可以如下式所示这样求得。
(算式18)
k ^ r ( n ) = Δ T ^ w , tmp ( n ) Δ T ^ w , rev ( n ) · · · ( 42 )
此外,若设工件温度变化量的第3临时推定值带帽ΔTw,rev(n)大致等于时刻n至时刻n+1的2个采样之间的工件温度变化量的真值的推定值带帽ΔTw时,则时刻n+1时的工件温度的真值的推定值带帽ΔTw(n+1)可以利用反向求解式(41)的下式来计算。
(算式19)
T ^ w ( n + 1 ) = Δ T ^ w , rev ( n ) + T ^ w ( n ) · · · ( 43 )
以下,通过反复进行同样的运算,可以继续计算出热阻值、热容量值、以及温度推定值的各个步骤的推定值。
接下来,对基于以上那样的原理的本实施例的温度推定装置50的动作进行说明,如图10所示,本实施例的温度推定装置50具有:存储部600、温度测量部601、工件温度变化量推定部602、工件温度推定部603、热阻值变化系数推定部604、以及热容量值变化系数推定部605。该存储部600预先存储所识别出的热传导模型的参数、测量数据和计算数据。
使用图11,对本实施例的温度推定装置50的动作进行说明。首先,温度测量部601,在正式运行前的初始状态下,取得由温度传感器4测量出的加热器温度Th、由环境温度传感器6测量出的环境温度Te、由仅可在初始状态下使用的温度传感器测量出的工件温度Tw。存储部600存储温度测量部601所取得的加热器温度Th、环境温度Te、工件温度Tw(图11步骤S10)。
接下来,正式运行开始,温度测量部601取得当前时刻n+1时的加热器温度Th(n+1)和环境温度Te(n+1)。存储部600存储温度测量部601所取得的加热器温度Th(n+1)和环境温度Te(n+1)(步骤S11)。
然后,在当前时刻n+1为运行开始后的最初采样时刻的情况下(步骤S12中为否),进入步骤S13,在当前时刻n+1为运行开始后的第2个以后的采样时刻的情况下(步骤S12中为是),进入步骤S15。
在这里,由于是最初采样时刻,所以工件温度变化量推定部602,参照存储部600所存储的值,计算当前时刻n+1的加热器温度Th(n+1)与1个采样前的时刻n的加热器温度Th(n)之间的差,即、计算2个采样之间的加热器温度变化量的真值ΔTh(n),并利用式(32)计算2个采样之间的加热器温度变化量的临时推定值带帽ΔTh,model(n),利用式(29)计算2个采样之间的工件温度变化量的第1临时推定值带帽ΔTw, model(n),根据这些计算值,利用式(36)计算时刻n至时刻n+1的2个采样之间的工件温度变化量的第2临时推定值带帽ΔTw,tmp(n)(步骤S13)。存储部600存储该工件温度变化量的第2临时推定值带帽ΔTw, tmp(n)。
此时,在当前时刻n+1为运行开始后的最初采样时刻的情况下,式(29)、式(32)中的加热器温度Th(n)为初始状态下所测量的温度,工件温度Tw(n)为初始状态下所测量的温度或者利用其他部位的温度近似的温度。此外,式(32)中的控制器输出值MV(n)是在时刻n时由控制装置51所计算出的值,不过,在当前时刻n+1为运行开始后的最初采样时刻的情况下,为初始状态下所输出的值。
接下来,热阻值变化系数推定部604,根据存储部600所存储的工件温度变化量的第2临时推定值带帽ΔTw,tmp(n)、时刻n时的工件温度的推定值带帽Tw(n)、加热器温度Th(n)以及环境温度Te(n),利用式(38)计算时刻n时的热阻值变化系数mr的推定值带帽mr(n)(步骤S14)。存储部600存储该热阻值变化系数mr的推定值带帽mr(n)。并且,在当前时刻n+1为运行开始后的最初采样时刻的情况下,作为工件温度的推定值带帽Tw(n),使用在初始状态下所测量出的或利用其他部位的温度近似的工件温度Tw
接下来,对经过1个采样周期,当前时刻n+1成为运行开始后的第2个采样时刻时的动作进行说明。步骤S11的处理如上所述。
当前时刻n+2为运行开始后的第2个采样时刻时,工件温度推定部603,根据2个采样之间的加热器温度变化量的真值ΔTh(n+1)、热阻值变化系数推定部604所计算出的热阻值变化系数mr的推定值带帽mr(n)、加热器温度Th(n)和控制器输出值MV(n),利用式(40)来计算1个周期之前的时刻n+1时的工件温度的推定值带帽Tw,rev(n+1)(步骤S15)。
在当前时刻n+2为运行开始后的第2个采样时刻的情况下,利用式(41),根据已经求得的时刻n的工件温度的真值的推定值带帽Tw(n)以及时刻n+1的工件温度的推定值带帽Tw,rev(n+1)来计算工件温度变化量的第3临时推定值ΔTw,rev(n)。另外,根据式(43)求工件温度的真值的推定值带帽Tw(n+1)。存储部600存储工件温度的推定值带帽Tw,rev(n+1)和工件温度的真值的推定值带帽Tw(n+1)。
然后,热容量值变化系数推定部605,根据工件温度变化量推定部602所计算出的工件温度变化量的第2临时推定值带帽ΔTw,tmp(n)和工件温度推定部603所计算出的工件温度变化量的第3临时推定值带帽ΔTw,rev(n),利用式(42)来计算2个周期之前的时刻n时的热容量值变化系数kr的推定值带帽kr(n)(步骤S16)。存储部600存储该热容量值变化系数kr的推定值kr(n)。
以下,同样地在加热冷却处理装置的动作停止之前(步骤S17中为是),在每个采样时刻都反复进行步骤S11~S16的处理。在当前时刻不是运行开始后的最初采样时刻的情况下,加热器温度Th(n)和环境温度Te(n)为1个采样之前所取得的值,作为式(29)、式(32)的工件温度Tw(n),式(38)、式(43)的工件温度的真值的推定值带帽Tw(n),使用利用1个时间点之前的加热器温度Th(n+1)而计算出的值。
控制器5的控制装置51,与这样的温度推定装置50的动作并行地取得温度推定装置50所计算出的时刻n时的工件温度的真值的推定值带帽Tw(n),并计算控制器输出值MV,使得该工件温度的推定值和目标温度一致,并在时刻n+2以后,作为MV(n+i)向加热器3输出(i>1)。
如以上那样,在本实施例中,即使不能直接测量工件1的温度、并且工件1和加热器3之间的热阻值及工件1的热容量值不明、热阻值和热容量值随每个工件1的不同而存在偏差,或者在正式运行中热阻值和热容量值从热传导模型识别时的值发生了变化的情况下,也能迅速且精度优良地推定工件1的温度。此外,在本实施例中,能够迅速且精度优良地计算出工件1和加热器3之间的热阻值的变化以及工件1的热容量值的变化。
此外,在本实施例中,即使不知道热阻值的变化量以及热容量值的变化量,也能迅速且精度优良地推定工件温度,所以,在温度推定中热阻值以及热容量值发生了变化的情况下,也能即时推定工件温度,利用其还能计算出连续变化的热阻值以及热容量值。此外,在本实施例中,在时间轴方向上仅进行少数几个采样的数值计算,在温度推定运算中未使用数值解析和统计运算,所以计算负担轻,从而也可以安装处理能力弱的处理器。
另外,在本实施例中,工件温度变化量的第2临时推定值、与根据在多个采样时刻所推定的工件温度而求得的工件温度变化量的第3临时推定值之比,为从热传导模型识别时发生了变化的、工件1的热容量值的变化值,据此,即使在工件1的热容量值不明、热容量值随每个工件的不同而存在偏差,或者在加热冷却处理中热容量值从热传导模型识别时的值发生了变化的情况下,也能迅速且精度优良地计算出工件1的热容量值。再有,在本实施例中,即使在加热冷却处理装置的温度控制处理器中,也能以足以进行实时处理的程度的简单运算来进行热容量值变化系数推定处理。
图12~图18表示本实施例的效果。图12、图13、图16是表示本实施例的工件温度的推定结果和工件温度的实测值的一个例子的图,图14、图17是表示本实施例的热阻值变化系数的推定结果和热阻值变化系数的实测值的一个例子的图,图15、图18是表示本实施例的热容量变化系数的推定结果与热容量值变化系数的实测值的一个例子的图。并且,图12~图18所示的实测值,全部为通过仿真而计算出的值,不过,在这里为了简化记载而全部称为实测值。
在图12~图15的例子中,假定了工件1和加热器3之间的热阻值以及工件1的热容量值在正式运行中从模型识别时的值发生了变化,并变为固定值的情况,由图可知,工件温度的推定值带帽Tw和工件温度的实测值Tw充分一致,热阻值变化系数的推定值带帽mr和热阻值变化系数的实测值mr充分一致,此外热容量值变化系数的推定值带帽kr和热容量值变化系数的实测值kr也充分一致,从而能迅速且精度优良地推定工件温度、热阻值变化系数、以及热容量值变化系数。
此外,在图16~图18的例子中,假定了工件1和加热器3之间的热阻值以及工件1的热容量值在温度推定中持续变化的情况,由图可知,工件温度、热阻值变化系数、以及热容量值变化系数的各推定值与各自的实测值充分一致,能够迅速且精度优良地推定工件温度、热阻值变化系数、以及热容量值变化系数。
在本实施例中,加热器温度变化量的实测值(真值)与基于热传导模型的加热器温度变化量的临时推定值之间的差,和基于热传导模型的工件温度变化量的第1临时推定值与工件温度变化量的第2临时推定值之间的差,呈以a/d表示的固定关系。
此外,在本实施例中,设工件1为推定对象、加热器3为温度可测量点,但是相反,在工件温度可测量、加热器温度不可测量时,在时刻n时的工件温度推定式(式(44))的基础上,使用1个采样之后的工件温度推定式(式(45)),并进一步使用时刻n+1时的热容量值变化系数kr的推定值带帽kr(n+1)与时刻n时的推定值带帽kr(n)大致相等、时刻n+1时的热阻值变化系数mr的推定值带帽mr(n+1)与时刻n时的推定值带帽mr(n)大致相等的关系,可以如式(46)那样计算出时刻n+1时的加热器温度的真值的推定值带帽Th(n+1)。
(算式20)
T w ( n + 1 ) = T w ( n ) + ( t S k ^ r ( n ) )
· ( m ^ r ( n ) · a · ( T h ( n ) - T w ( n ) ) - b · ( T w ( n ) - T e ( n ) ) ) · · ( 44 )
(算式21)
T w ( n + 2 ) = T w ( n + 1 ) + ( t S k ^ r ( n ) )
· ( m ^ r ( n ) · a · ( T h ( n + 1 ) - T w ( n + 1 ) ) - b · ( T w ( n + 1 ) - T e ( n + 1 ) ) ) · · · ( 45 )
(算式22)
T ^ h ( n + 1 ) = T h ( n ) + t S · ( c · MV ( n ) - m r ( n ) · d · ( T h ( n ) - T w ( n ) ) ) · · · ( 46 )
对于感应加热装置等而言,由于在工件侧产生热量,所以可以将工件看做加热源,并且往往不可测量工件温度。在这种情况下,设工件温度为Tw、装置可测量点的温度为Th,使用式(46)来推定工件温度的方法是适当的。
此外,在本实施例中,温度测量要素仅为加热器温度、工件温度和环境温度,不过,在正式运行中模型从模型识别时发生的变化只是加热器3与工件1之间的热阻值、和工件1的热容量值的情况下,即使在热传导模型式上追加其他的可测量点,式子的未知量的数量也不变,所以,其结果是可以导出式(35),从而可以应用相同的方法。
再有,在本实施例中,以加热器温度直接受控制器输出值的影响的方式进行模型化,不过,如下述的式(47)所示,控制器的输出与加热器之间也可以插入其他的模型要素。在式(47)中,设模型要素p的温度为Tp(n)。只要可以测量温度Tp(n),则此前的讨论内容可以原封不动地使用。
Th(n+1)-Th(n)=ts·(c·(Tp(n)-Th(n))-mr·d·(Th(n)-Tw(n)))..(47)
此外,在实施例1和实施例2中,使用环境温度传感器6来测量环境温度,不过,在环境温度可看做一定值的情况下,不测量环境温度,而使用固定值进行计算,也不会影响计算方法。
另外,在实施例1和实施例2中,没有考虑温度测量值的测量误差,但是,如果将使用以下滤波器来消除了测量误差的结果考虑为温度测量值,则不失讨论内容的一般性,该滤波器是与温度变化倾向和热阻值变化速度相比足够快的滤波器。
另外,在实施例1和实施例2中,主要针对加热处理进行了记载,不过,在冷却处理中也可以应用本发明。
(产业上的可利用性)
本发明能够应用到半导体制造装置等加热冷却处理装置中。

Claims (12)

1.一种温度推定方法,其特征在于,具有:
测量步骤,测量与推定对象之间具有热阻的温度可测量点的温度;
温度变化量推定步骤,根据上述温度可测量点的温度变化量、温度可测量点的温度变化量临时推定值、和推定对象的温度变化量临时推定值,来推定上述推定对象的温度变化量的真值,其中,上述温度可测量点的温度变化量是根据上述温度可测量点的温度而求得的,温度可测量点的温度变化量临时推定值是基于上述温度可测量点和上述推定对象的热传导模型来推定上述温度可测量点的温度变化量而得到的值,推定对象的温度变化量临时推定值是基于上述热传导模型来推定上述推定对象的温度变化量而得到的值;以及
温度推定步骤,根据在该温度变化量推定步骤中所推定的温度变化量的真值,来推定上述推定对象的温度的真值。
2.根据权利要求1所述的温度推定方法,其特征在于,
上述温度变化量推定步骤包括下述步骤:利用上述温度可测量点的温度变化量的实测值与上述温度可测量点的温度变化量临时推定值之间的差、和上述推定对象的温度变化量临时推定值与上述推定对象的温度变化量的真值之间的差存在一定关系,来推定上述推定对象的温度变化量的真值。
3.根据权利要求1所述的温度推定方法,其特征在于,
还具有热阻值变化系数推定步骤,在该热阻值变化系数推定步骤中,根据上述推定对象的温度的真值和上述热传导模型,来推定上述温度可测量点与上述推定对象之间的热阻的变化值。
4.一种温度推定装置,其特征在于,具有:
存储部,存储关于温度可测量点和推定对象的热传导模型的参数,该温度可测量点与上述推定对象之间具有热阻;
温度测量部,测量上述温度可测量点的温度;
温度变化量推定部,根据上述温度可测量点的温度变化量、温度可测量点的温度变化量临时推定值、和推定对象的温度变化量临时推定值,来推定上述推定对象的温度变化量的真值,其中,上述温度可测量点的温度变化量是根据上述温度可测量点的温度而求得的,温度可测量点的温度变化量临时推定值是基于上述热传导模型来推定上述温度可测量点的温度变化量而得到的值,推定对象的温度变化量临时推定值是基于上述热传导模型来推定上述推定对象的温度变化量而得到的值;以及
温度推定部,根据由该温度变化量推定部所推定的温度变化量的真值来推定上述推定对象的温度的真值。
5.根据权利要求4所述的温度推定装置,其特征在于,
上述温度变化量推定部,利用上述温度可测量点的温度变化量的实测值与上述温度可测量点的温度变化量临时推定值之间的差、和上述推定对象的温度变化量临时推定值与上述推定对象的温度变化量的真值之间的差存在一定关系,来推定上述推定对象的温度变化量的真值。
6.根据权利要求4所述的温度推定装置,其特征在于,
还包括热阻值变化系数推定部,该热阻值变化系数推定部根据上述推定对象的温度的真值和上述热传导模型,来推定上述温度可测量点与上述推定对象之间的热阻的变化值。
7.一种温度推定方法,其特征在于,具有:
测量步骤,测量与推定对象之间具有热阻的温度可测量点的温度;
温度变化量推定步骤,根据上述温度可测量点的温度变化量的实测值、上述温度可测量点的温度变化量临时推定值、和上述推定对象的温度变化量第1临时推定值,来求取上述推定对象的温度变化量第2临时推定值,上述推定对象的温度变化量第2临时推定值是视上述推定对象的热容量值未发生变化时的值,其中,上述温度可测量点的温度变化量的实测值是根据上述温度可测量点的温度而求得的,上述温度可测量点的温度变化量临时推定值是基于上述温度可测量点和上述推定对象的热传导模型进行推定而得到的值,上述推定对象的温度变化量第1临时推定值是基于上述热传导模型进行推定而得到的值;
热阻值变化系数推定步骤,根据上述推定对象的温度变化量第2临时推定值和上述热传导模型来推定上述热阻值的变化值;以及
温度推定步骤,根据不同采样时刻的上述温度可测量点的温度、和上述热阻值的变化值来推定上述推定对象的温度的真值。
8.根据权利要求7所述的温度推定方法,其特征在于,
上述温度变化量推定步骤还包括下述步骤:利用上述温度可测量点的温度变化量的实测值与上述温度可测量点的温度变化量临时推定值之间的差、和上述推定对象的温度变化量第1临时推定值与上述推定对象的温度变化量第2临时推定值之间的差存在一定关系,来推定上述推定对象的温度变化量第2临时推定值。
9.根据权利要求7所述的温度推定方法,其特征在于,
还具有热容量值变化系数推定步骤,在该热容量值变化系数推定步骤中,根据上述推定对象的温度变化量第3临时推定值和上述温度变化量第2临时推定值来推定上述推定对象的热容量值的变化值,上述推定对象的温度变化量第3临时推定值是根据在多个采样时刻所推定的上述推定对象的温度而求得的。
10.一种温度推定装置,其特征在于,具有:
存储部,存储关于温度可测量点和推定对象的热传导模型的参数,该温度可测量点与上述推定对象之间具有热阻;
温度测量部,测量上述温度可测量点的温度;
温度变化量推定部,根据上述温度可测量点的温度变化量的实测值、上述温度可测量点的温度变化量临时推定值、和上述推定对象的温度变化量第1临时推定值,来求取上述推定对象的温度变化量第2临时推定值,上述推定对象的温度变化量第2临时推定值是视上述推定对象的热容量值未发生变化时的值,其中,上述温度可测量点的温度变化量的实测值是根据上述温度可测量点的温度而求得的,上述温度可测量点的温度变化量临时推定值是基于上述热传导模型进行推定而得到的值,上述推定对象的温度变化量第1临时推定值是基于上述热传导模型进行推定而得到的值;
热阻值变化系数推定部,根据上述推定对象的温度变化量第2临时推定值和上述热传导模型来推定上述热阻值的变化值;以及
温度推定部,根据不同采样时刻的上述温度可测量点的温度、和上述热阻值的变化值来推定上述推定对象的温度的真值。
11.根据权利要求10所述的温度推定装置,其特征在于,
上述温度变化量推定部,利用上述温度可测量点的温度变化量的实测值与上述温度可测量点的温度变化量临时推定值之间的差、和上述推定对象的温度变化量第1临时推定值与上述推定对象的温度变化量第2临时推定值之间的差存在一定关系,来推定上述推定对象的温度变化量第2临时推定值。
12.根据权利要求10所述的温度推定装置,其特征在于,
还具有热容量值变化系数推定部,该热容量值变化系数推定部根据上述推定对象的温度变化量第3临时推定值和上述温度变化量第2临时推定值来推定上述推定对象的热容量值的变化值,上述推定对象的温度变化量第3临时推定值是根据在多个采样时刻所推定的上述推定对象的温度而求得的。
CN2006800512893A 2006-01-19 2006-12-26 温度推定方法及装置 Active CN101360983B (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
JP2006010898A JP4658818B2 (ja) 2006-01-19 2006-01-19 温度推定方法および装置
JP010898/2006 2006-01-19
PCT/JP2006/325936 WO2007083498A1 (ja) 2006-01-19 2006-12-26 温度推定方法および装置

Publications (2)

Publication Number Publication Date
CN101360983A CN101360983A (zh) 2009-02-04
CN101360983B true CN101360983B (zh) 2010-08-04

Family

ID=38287454

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2006800512893A Active CN101360983B (zh) 2006-01-19 2006-12-26 温度推定方法及装置

Country Status (4)

Country Link
JP (1) JP4658818B2 (zh)
KR (1) KR101026432B1 (zh)
CN (1) CN101360983B (zh)
WO (1) WO2007083498A1 (zh)

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4939456B2 (ja) * 2008-02-21 2012-05-23 株式会社山武 温度推定方法および温度推定装置
JP5578028B2 (ja) * 2010-10-29 2014-08-27 セイコーエプソン株式会社 温度測定装置および温度測定方法
CN102735358B (zh) * 2011-03-29 2016-03-30 株式会社爱德克斯 温度推定装置以及温度推定方法
JP5737331B2 (ja) 2013-05-27 2015-06-17 横浜ゴム株式会社 加硫制御方法および加硫制御システム
DE102013213448B4 (de) * 2013-07-09 2021-12-09 Siemens Aktiengesellschaft Elektronikbaugruppe mit Leistungshalbleiter
JP6123741B2 (ja) * 2014-06-20 2017-05-10 トヨタ自動車株式会社 冷却器
JP5952438B1 (ja) * 2015-01-22 2016-07-13 ファナック株式会社 電動機の温度推定装置
JP6163590B1 (ja) * 2016-05-31 2017-07-12 株式会社日阪製作所 シミュレーション方法、シミュレーションプログラム、及びこのプログラムを内蔵した記憶媒体を含むシミュレーション装置
CN106225943B (zh) * 2016-08-31 2018-12-14 南京理工大学连云港研究院 一种基于瞬态传热理论测量高温石英熔融炉温度的装置及方法
WO2018186172A1 (ja) * 2017-04-03 2018-10-11 シャープ株式会社 電子機器
TWI646315B (zh) 2017-11-15 2019-01-01 財團法人工業技術研究院 溫度估算裝置及溫度估算方法
JP7326344B2 (ja) * 2018-05-24 2023-08-15 アプライド マテリアルズ インコーポレイテッド 空間分解ウエハ温度制御のための仮想センサ
JP7071946B2 (ja) * 2019-06-21 2022-05-19 東京エレクトロン株式会社 プラズマ処理装置
US20220412815A1 (en) * 2021-06-25 2022-12-29 Rosemount Inc. Heat flow-based process fluid temperature estimation system with thermal time response improvement

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1322940A (zh) * 2000-05-09 2001-11-21 丰田自动车株式会社 温度推定方法及其设备
CN1647361A (zh) * 2002-03-20 2005-07-27 丰田自动车株式会社 电动机温度推定装置和电动机控制装置

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP3331758B2 (ja) * 1994-07-25 2002-10-07 株式会社デンソー 排出ガス浄化装置の温度制御装置
JPH07283158A (ja) * 1994-04-11 1995-10-27 Tokyo Electron Ltd 熱処理装置およびその温度制御方法
JP2000088666A (ja) * 1998-09-11 2000-03-31 Fujikura Ltd 地中管路内電力ケーブルの導体温度計算方法及び装置
JP3403160B2 (ja) * 2000-09-27 2003-05-06 東京エレクトロン株式会社 熱処理装置、熱処理装置用制御装置
JP3812884B2 (ja) * 2001-06-27 2006-08-23 日本電信電話株式会社 温度検出方法および温度検出器
JP4285759B2 (ja) * 2003-07-28 2009-06-24 株式会社日立国際電気 基板処理装置及び基板処理方法
JP4712343B2 (ja) * 2003-10-30 2011-06-29 東京エレクトロン株式会社 熱処理装置、熱処理方法、プログラム及び記録媒体

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1322940A (zh) * 2000-05-09 2001-11-21 丰田自动车株式会社 温度推定方法及其设备
CN1647361A (zh) * 2002-03-20 2005-07-27 丰田自动车株式会社 电动机温度推定装置和电动机控制装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
JP特开平7-283158A 1995.10.27

Also Published As

Publication number Publication date
CN101360983A (zh) 2009-02-04
KR20080081027A (ko) 2008-09-05
WO2007083498A1 (ja) 2007-07-26
JP2007192661A (ja) 2007-08-02
KR101026432B1 (ko) 2011-04-07
JP4658818B2 (ja) 2011-03-23

Similar Documents

Publication Publication Date Title
CN101360983B (zh) 温度推定方法及装置
EP3278032B1 (en) Method and system for determining characteristic parameters of a hydraulic network
CN106471733A (zh) 用于确定转子温度的方法和装置、计算机程序、计算机程序产品
KR101086128B1 (ko) 제어 기기 및 전력 추정 방법
CN104535257B (zh) 一种硅压阻温度补偿评估方法
CN111649794B (zh) 热式流量计及流量校正方法
Polom et al. Exploiting distinct thermal response properties for power semiconductor module health monitoring
CN113935253A (zh) 一种基于数据赋权的页岩气井经验产量递减模型拟合方法
Groover et al. A continuing study in the determination of temperatures in metal cutting using remote thermocouples
CN109270460A (zh) 锂离子电池热失控的能量计算方法
CN111208425A (zh) 高精度异步电机系统状态模型的构建方法及异步电机状态检测方法
CN105404329A (zh) 用于对成型模具进行调温的方法和装置
CN109270455B (zh) 基于弱敏集合卡尔曼滤波的感应电机状态监测方法
CN107703761B (zh) 一种气动调节阀粘滞特性参数的估计方法
CN113172119A (zh) 模具工作中折弯处恒温控制系统、控制方法、终端、介质
CN108169565A (zh) 一种用于电导率测量的非线性温度补偿方法
JP2008151739A (ja) 温度推定方法および装置
CN110376241B (zh) 一种振荡式量热反应釜传热因子和系统热容测量方法
Špička et al. The mathematical-physical models and the neural network exploitation for time prediction of cooling down low range specimen
CN109031942A (zh) 一种利用传递函数修正量热仪量热过程延迟环节的方法
Gilev et al. Mathematical model of induction heating with heat transfer of cylindrical body for pressing treatment
CN115307768B (zh) 一种ntc检测时间补偿算法
Takenaga et al. Process Optimizer for adjusting film thickness and in-film dopant concentration at the same time
JPH08132506A (ja) 射出成形用金型の熱伝達係数算出方法
Wu et al. A quick surface heat flux estimation method for typical structures of high-speed aircrafts

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant