CN104050320A - 一种电熨斗底板温度分布的计算方法 - Google Patents

一种电熨斗底板温度分布的计算方法 Download PDF

Info

Publication number
CN104050320A
CN104050320A CN201410265935.2A CN201410265935A CN104050320A CN 104050320 A CN104050320 A CN 104050320A CN 201410265935 A CN201410265935 A CN 201410265935A CN 104050320 A CN104050320 A CN 104050320A
Authority
CN
China
Prior art keywords
physical model
iron
shaped pipe
temperature
dull
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.)
Granted
Application number
CN201410265935.2A
Other languages
English (en)
Other versions
CN104050320B (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.)
China Jiliang University
Original Assignee
China Jiliang University
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 China Jiliang University filed Critical China Jiliang University
Priority to CN201410265935.2A priority Critical patent/CN104050320B/zh
Publication of CN104050320A publication Critical patent/CN104050320A/zh
Application granted granted Critical
Publication of CN104050320B publication Critical patent/CN104050320B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • General Induction Heating (AREA)
  • Investigating Or Analyzing Materials Using Thermal Means (AREA)

Abstract

本发明公开了一种电熨斗底板温度分布的计算方法。将电熨斗底板视作截面相同的平板,建立平板物理模型A;将电熨斗加热管视作U形管,建立U形管物理模型B;根据电熨斗加热管与电熨斗底板的位置,将U形管物理模型B平行嵌入到平板物理模型A上建立有限元模型C,根据有限元模型C通过有限元法计算得到平板下表面中心线的中点温度作为关键点温度,计算各时刻下的平板温度分布作为电熨斗的底板温度分布。本发明通过测量所得的参数作为基本条件,避免了电熨斗附加热电偶及导热硅油对底板温度场的影响,能够完整地计算得到任意时刻任意位置底板温度的分布,精度高,操作方便。

Description

一种电熨斗底板温度分布的计算方法
技术领域
本发明涉及一种温度分布的计算方法,尤其是涉及一种电熨斗底板温度分布的计算方法。
背景技术
电熨斗已经成为现代生活必备的家庭工具之一,而且现在的电熨斗不仅仅具有加热功能,还包括其它附加功能如:恒温、蒸汽、喷雾等等,这些功能的实现要有一个基本的条件才能实现,就是电熨斗底板的精确测温,所以电熨斗底板温度分布的检测也是科学研究的一个重要方向。
目前对电熨斗的检测方式还停留在用热电偶进行接触测量,本文创新了一种方法,可以对不同结构的底板温度进行计算,通过与福禄克Ti25热像仪进行比较,温度的误差在5℃以内,为今后设计新的电熨斗以及底板温度检测提供一种新的可行方法。
现有的对电熨斗底板温度检测的方法,存在的缺陷和不足之处如下:
由于接触测温要求测量热电偶放置于电熨斗底板上,一般需要附着上导热硅油,用聚酯胶固定,操作不方便,还要影响电熨斗的底板温度场;检测时只测量电熨斗底板四个点的温度,不能很好反应整个底板的温度状态;存在进一步改进空间。
发明内容
针对背景技术中存在的不足,本发明的目的是提供一种电熨斗底板温度分布的计算方法,简单有效,可以通过测量电熨斗的外观尺寸、加热管尺寸、功率等信息而进行温度预测,减少不必要的实验和材料。
为解决上述问题,本发明的电熨斗底板内嵌有电熨斗加热管,包括以下步骤:
1)将电熨斗底板视作截面相同的平板,建立平板物理模型A,平板物理模型A的参数包括平板的厚度d、底面积S、密度ρ、导热系数k和比热c、室温Tf以及平板对空气的对流换热系数hf,平板的厚度d、底面积S、密度ρ、导热系数k和比热c均与电熨斗底板的属性相同;
2)将电熨斗加热管视作U形管,建立U形管物理模型B,U形管物理模型B的参数包括U形管中心线与U形管一侧直管部分的中轴线之间的距离r1、U形管两侧直管部分的半径r2和U形管两侧直管部分的长度h;
3)根据电熨斗加热管与电熨斗底板的位置,将U形管物理模型B平行嵌入到平板物理模型A上建立有限元模型C,有限元模型C中,U形管物理模型B和平板物理模型A的对称中心线重合,U形管物理模型B嵌入平板物理模型A的深度为U形管两侧直管部分的半径r2
4)根据有限元模型C通过有限元法计算得到平板下表面中心线的中点温度作为关键点温度,计算各时刻下的平板温度分布作为电熨斗的底板温度分布。
所述的步骤4)中计算各时刻下的平板温度分布具体采用以下方法:
设定温度常数T,当该关键点温度大于温度常数T时,将近似热流密度qw代入所述的有限元模型C中计算当前时刻下的平板温度;当该关键点温度小于温度常数T时,将近似热流密度qw等于零代入所述的有限元模型C中计算当前时刻下的平板温度。
所述的温度常数T为电熨斗内的双金属片控温元件的断开温度。
所述的近似热流密度qw根据电熨斗的功率Q和U形管物理模型B通过以下公式计算得到:
q w = Q 2 π r 2 ( π r 1 + 2 h ) · R w
其中,Rw=0.7为电熨斗加热管的加热丝热量传递到电熨斗加热管的外部套管上热量的百分比。
本发明的有益效果是:
本发明是在电熨斗底板温度计算过程中,通过测量所得的参数作为基本条件,本发明方法避免了电熨斗附加热电偶及导热硅油对底板温度场的影响,能够完整地计算得到底板温度的分布,精度高,操作方便。
附图说明
图1为平板物理模型A的示意图。
图2为U形管物理模型B的截面示意图。
图3为有限元模型C的示意图。
具体实施方式
下面结合附图及具体实施例对本发明作进一步详细说明。
本发明中电熨斗底板内嵌有电熨斗加热管,本发明方法包括以下步骤,
1)如图1所示,将电熨斗底板视作截面相同的平板,建立平板物理模型A,平板物理模型A的参数包括平板的厚度d、底面积S、密度ρ、导热系数k和比热c、室温Tf以及平板对空气的对流换热系数hf,平板的厚度d、底面积S、密度ρ、导热系数k和比热c均与电熨斗底板的属性相同;
2)如图2所示,将电熨斗加热管视作U形管,建立U形管物理模型B,U形管为圆柱形管,U形管物理模型B的参数包括U形管中心线与U形管一侧直管部分的中轴线之间的距离r1、U形管两侧直管部分的半径r2和U形管两侧直管部分的长度h;
3)如图3所示,根据电熨斗加热管与电熨斗底板的位置,将U形管物理模型B平行嵌入到平板物理模型A上建立有限元模型C,有限元模型C中,U形管物理模型B和平板物理模型A的对称中心线重合,U形管物理模型B嵌入平板物理模型A的深度为U形管两侧直管部分的半径r2
4)根据有限元模型C通过有限元法计算得到平板下表面中心线的中点温度作为关键点温度,平板下表面等同为电熨斗底板的底面,计算各时刻下的平板温度分布作为电熨斗的底板温度分布。
所述的步骤4)中计算各时刻下的平板温度分布具体采用以下方法:
设定温度常数T,当该关键点温度大于温度常数T时,将近似热流密度qw代入所述的有限元模型C中计算当前时刻下的平板温度;当该关键点温度小于温度常数T时,将近似热流密度qw等于零代入所述的有限元模型C中计算当前时刻下的平板温度。
所述的温度常数T为电熨斗内的双金属片控温元件的断开温度。
所述的近似热流密度qw根据电熨斗的功率Q和U形管物理模型B通过以下公式计算得到:
q w = Q 2 π r 2 ( π r 1 + 2 h ) · R w
其中,Rw为电熨斗加热管的加热丝热量传递到电熨斗加热管的外部套管上热量的百分比,其数值跟内部氧化镁填充物的压实程度有关。
本发明的原理如下:
1)假定电熨斗初始时刻温度分布均匀,且为室温Tc=20℃,hf为平板对空气的对流换热系数,简称为换热系数。底板上表面,受到加热管的热流密度为qw的间断性加热;底板下表面处始终受到空气对流的冷却,得到以下边界条件:
其中,Tf为与底板下表面相交换的空气温度,hf为对应的换热系数。
2)假定底板各向材料性能相同,材料导热系数相同,热扩散率a=k/ρc=9.63×10-5m2/s。以底板上表面中心轴线与底板末端的交点为原点,中心轴线方向为y轴,沿底板方向为z轴,垂直上下表面方向为x轴,建立控制微分方程为以下公式:
∂ T ∂ t = a ( ∂ 2 T ∂ x 2 + ∂ 2 T ∂ y 2 + ∂ 2 T ∂ z 2 )
给定初始条件Tf=20℃和边界条件hf=15W/(m2·℃),即可得到当前时刻温度的唯一解。
本发明的实施例如下:
根据图1~图3所示的电熨斗底板模型,包括电熨斗底板、电熨斗加热管,所述电熨斗加热管嵌入所述电熨斗底板上。电熨斗底板温度分布计算方法具体步骤如下:
1)测量得电熨斗底板视为厚度d,其值为0.004m,底面积S,其值为0.115m2,密度ρ=2700kg/m3的平板;
其三维物理模型示意图如1所示,底板上表面受到加热管热源的加热,其热流密度为qw,下表面受到空气的对流换热,其表面传热系数为hf=15W/(m2·℃),空气温度为Tf=20℃;
2)设温控器接通后的加热管,电熨斗加热管的加热丝热量传递到电熨斗加热管的外部套管上热量的百分比Rw=0.7;
3)底板的导热系数k=237W/(m·℃)、比热c=900J/(kg·℃)为常数,不随温度变化而变化;
4)加热管为U型圆柱体,嵌于底板上表面,其发热量可用热流密度表示,加载在与底板接触的内表面上,发热管中心截面示意图如2所示,其热流密度qw,其值为:
q w = Q 2 π r 2 ( π r 1 + 2 h ) · R w = 1.337e5 W / m 2 ;
7)电熨斗的底板视为平板后,其下表面中心轴线中点为关键点:
在0~130s中,当此关键点的计算温度小于175℃时,始终在加热管加载qw,其值qw=1.337e5W/m2,当此关键点的计算温度大于230℃时,对加热管不加载热流密度。
在130~400s中,当此关键点的计算温度小于175℃时,始终在加热管加载qw,其值qw=1.337e5W/m2,当此关键点的计算温度大于220℃时,对加热管不加载热流密度。
然后实际中,用福禄克的红外热像仪Ti25测量电熨斗的实际底板温度,将实施例得到的第一次峰值温度电熨斗的温度极值(t=39s)并与实际温度进行比较,得到以下数据,如表1所示:
表1
将实施例得到的第一次峰值温度电熨斗的温度极值(t=149s)并与实际温度进行比较,得到以下数据,如表2所示:
表2
将实施例得到的第一次峰值温度电熨斗的温度极值(t=158s)并与实际温度进行比较,得到以下数据,如表3所示:
表3
这种用有限元法计算电熨斗底板温度分布的方法计算方便,精度高,避免了导热硅胶和热电偶对温度的影响。
上述具体实施方式用来解释说明本发明,而不是对本发明进行限制,在本发明的精神和权利要求的保护范围内,对本发明作出的任何修改和改变,都落入本发明的保护范围。

Claims (4)

1.一种电熨斗底板温度分布的计算方法,电熨斗底板内嵌有电熨斗加热管,其特征在于包括以下步骤:
1)将电熨斗底板视作截面相同的平板,建立平板物理模型A,平板物理模型A的参数包括平板的厚度d、底面积S、密度ρ、导热系数k和比热c、室温Tf以及平板对空气的对流换热系数hf,平板的厚度d、底面积S、密度ρ、导热系数k和比热c均与电熨斗底板的属性相同;
2)将电熨斗加热管视作U形管,建立U形管物理模型B,U形管物理模型B的参数包括U形管中心线与U形管一侧直管部分的中轴线之间的距离r1、U形管两侧直管部分的半径r2和U形管两侧直管部分的长度h;
3)根据电熨斗加热管与电熨斗底板的位置,将U形管物理模型B平行嵌入到平板物理模型A上建立有限元模型C,有限元模型C中,U形管物理模型B和平板物理模型A的对称中心线重合,U形管物理模型B嵌入平板物理模型A的深度为U形管两侧直管部分的半径r2
4)根据有限元模型C通过有限元法计算得到平板下表面中心线的中点温度作为关键点温度,计算各时刻下的平板温度分布作为电熨斗的底板温度分布。
2.根据权利要求1所述的一种电熨斗底板温度分布的计算方法,其特征在于:所述的步骤4)中计算各时刻下的平板温度分布具体采用以下方法:
设定温度常数T,当该关键点温度大于温度常数T时,将近似热流密度qw代入所述的有限元模型C中计算当前时刻下的平板温度;当该关键点温度小于温度常数T时,将近似热流密度qw等于零代入所述的有限元模型C中计算当前时刻下的平板温度。
3.根据权利要求2所述的一种电熨斗底板温度分布的计算方法,其特征在于:所述的温度常数T为电熨斗内的双金属片控温元件的断开温度。
4.根据权利要求2所述的一种电熨斗底板温度分布的计算方法,其特征在于:所述的近似热流密度qw根据电熨斗的功率Q和U形管物理模型B通过以下公式计算得到:
q w = Q 2 π r 2 ( π r 1 + 2 h ) · R w
其中,Rw为电熨斗加热管的加热丝热量传递到电熨斗加热管的外部套管上热量的百分比。
CN201410265935.2A 2014-06-14 2014-06-14 一种电熨斗底板温度分布的计算方法 Expired - Fee Related CN104050320B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410265935.2A CN104050320B (zh) 2014-06-14 2014-06-14 一种电熨斗底板温度分布的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410265935.2A CN104050320B (zh) 2014-06-14 2014-06-14 一种电熨斗底板温度分布的计算方法

Publications (2)

Publication Number Publication Date
CN104050320A true CN104050320A (zh) 2014-09-17
CN104050320B CN104050320B (zh) 2017-02-15

Family

ID=51503148

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410265935.2A Expired - Fee Related CN104050320B (zh) 2014-06-14 2014-06-14 一种电熨斗底板温度分布的计算方法

Country Status (1)

Country Link
CN (1) CN104050320B (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050056807A1 (en) * 2001-08-13 2005-03-17 New Qu Energy, Ltd. Medium having a high heat transfer rate
CN101290299A (zh) * 2007-04-16 2008-10-22 上海轮胎橡胶(集团)股份有限公司轮胎研究所 一种可变导热系数的测量装置及方法
CN102661969A (zh) * 2012-04-28 2012-09-12 宝山钢铁股份有限公司 一种测试不同水膜厚度钢板换热系数的方法和装置
CN103411996A (zh) * 2013-08-05 2013-11-27 电子科技大学 固体材料导热系数测量装置及测量方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050056807A1 (en) * 2001-08-13 2005-03-17 New Qu Energy, Ltd. Medium having a high heat transfer rate
CN101290299A (zh) * 2007-04-16 2008-10-22 上海轮胎橡胶(集团)股份有限公司轮胎研究所 一种可变导热系数的测量装置及方法
CN102661969A (zh) * 2012-04-28 2012-09-12 宝山钢铁股份有限公司 一种测试不同水膜厚度钢板换热系数的方法和装置
CN103411996A (zh) * 2013-08-05 2013-11-27 电子科技大学 固体材料导热系数测量装置及测量方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
郑坚璐: ""电热型家电产品热性能分析及测试方法研究"", 《中国优秀硕士学位论文全文数据库》 *
郑坚璐等: ""电熨斗在线温度检测系统的设计"", 《传感器与微系统》 *

Also Published As

Publication number Publication date
CN104050320B (zh) 2017-02-15

Similar Documents

Publication Publication Date Title
CN101113963B (zh) 一种测量液体导热系数的方法及其装置
Yu et al. Thermal response test and numerical analysis based on two models for ground-source heat pump system
Lu et al. Inverse estimation of the inner wall temperature fluctuations in a pipe elbow
CN108490024B (zh) 一种基于虚拟热源原理测量有限厚度材料异质含量的方法
CN103940847B (zh) 一种基于热流传感器各向异性薄膜热导率测试方法及装置
CN202814902U (zh) 一种沥青混合料导温系数测试装置
CN104062034A (zh) 一种基于管道外壁温度测量的非接触式管程流体温度测量方法
CN105004949A (zh) 一种在线运行耐张线夹最大载流量的测试方法和测试装置
CN105004749B (zh) 固‑液相变材料熔化传热性能参数测试系统及其方法
CN103196949A (zh) 热阻式热流计校准方法及其实施装置
CN102521439B (zh) 结合有限元法和反传热法计算淬火介质换热系数的方法
US20170184432A1 (en) Flow speed measurement method and flow speed measurement system
Taler Determination of local heat transfer coefficient from the solution of the inverse heat conduction problem
CN107643131A (zh) 一种热管式温度测定装置及方法
CN105675646B (zh) 基于本征光热信息同时测量高温半透明介质热导率及吸收系数的方法
CN106645284A (zh) 一种圆管材料导热系数测量系统及其测量方法
JP2017053758A (ja) 関係情報設定方法、流速計測方法、関係情報設定システム及び流速計測システム
CN104050320A (zh) 一种电熨斗底板温度分布的计算方法
CN103995017B (zh) 一种测定周期性换热系数的实验方法
CN105463142B (zh) 一种高炉炉缸内铁水温度测量的方法
CN108051475A (zh) 一种对流换热系数的快速测量方法
CN106885634A (zh) 基于红外热像测温技术的非定常壁面加热热流分布测量方法
CN107543841A (zh) 一种测量池中流体导热系数的新型便携式测量仪
CN103135631A (zh) 基于能量等量交换的软化水箱用软化水预热系统
CN207114051U (zh) 一种便携式热电偶/热电阻检测用精密恒温箱

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
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20170215