CN204086184U - 一种可同时测固体材料导热系数及热扩散率的系统 - Google Patents
一种可同时测固体材料导热系数及热扩散率的系统 Download PDFInfo
- Publication number
- CN204086184U CN204086184U CN201420475647.5U CN201420475647U CN204086184U CN 204086184 U CN204086184 U CN 204086184U CN 201420475647 U CN201420475647 U CN 201420475647U CN 204086184 U CN204086184 U CN 204086184U
- Authority
- CN
- China
- Prior art keywords
- water bath
- heat conductivity
- solid material
- thermal diffusivity
- temperature
- 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.)
- Expired - Fee Related
Links
Landscapes
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
Abstract
本实用新型涉及一种可同时测固体材料导热系数及热扩散率的系统,即利用超级恒温水浴作为热源,与试样盒底部的黄铜板形成恒温边界,温度均匀性和稳定性好,且易于实现,直接基于恒温边界一维非稳态传热模型,利用参数估计法同时测得固体材料的导热系数和热扩散率,且测得的结果可靠、准确、精度高。
Description
技术领域
本实用新型属于固体材料热物性参数测试技术领域,具体涉及一种可同时测固体材料导热系数及热扩散率的系统。
背景技术
随着我国对节能及能源利用效率的日益重视,材料热物性及其测试方法方面的研究方兴未艾。固体材料热物性测试总体上分为稳态法和非稳态法,稳态法测试由于所需时间长、测试环境要求高,而不能适应现代科技迅速发展的需要,非稳态法应运而生并因为其具有快速、准确、可实现多参数测量的优势而得到了更为广泛的应用。然而,非稳态法在实际应用中也存在一些问题,例如目前多采用加热板、加热丝、加热棒等固态电热器作为加热热源,存在着热源热流密度的准确给定、加热探头热容及时间延迟性、热丝电阻温度依变性和功率波动等问题,这些问题在一定程度上影响了热物性参数测试的精度。近年来,发展了一种陶瓷基底刻蚀双螺旋铂金属丝加热膜,热稳定性好且厚度极薄,可有效提高材料热物性测试精度,但成本不菲且制作困难,难以推广应用。因此,目前尚无一种精度理想的可同时测固体材料导热系数及热扩散率的系统。
实用新型内容
本实用新型要解决的技术问题是提供一种同时测固体材料导热系数及热扩散率的系统,以克服不能同时测试固体材料导热系数及热扩散率以及测 试精度不理想的缺陷。
为解决上述技术问题,本实用新型采用以下技术方案:
一种可同时测固体材料导热系数及热扩散率的系统,其特征在于:包括超级恒温水浴系统、试样盒和温度采集与分析系统,所述超级恒温水浴系统包括水介质(1)、加热元件(2)、温度传感器(3)、温控单元(4)、搅拌器(5)、绝热水浴箱(6)组成,所述水介质(1)装在所述绝热水浴箱(6)中,所述加热元件(2)安装在所述绝热水浴箱(6)的内侧下方,与所述温控单元(4)电连接,所述温控单元(4)安装在所述绝热水浴箱(6)的顶部,所述温控单元(4)输出端与所述温度传感器(3)电连接,所述温度传感器(3)输出端进入所述绝热水浴箱(6),所述搅拌器(5)向内安装在所述绝热水浴箱(6)的底部;所述试样盒安装在所述绝热水浴箱(6)的顶部,由底部黄铜板(7)、侧绝热板(8)和上绝热板(9)组成,所述温度采集与分析系统包括热电偶(11)、信号采集与传输单元(12)和微型计算机(13),所述热电偶(11)一端通过所述上绝热板(9)进入所述试样盒,另一端通过所述采集与传输单元(12)与微型计算机(13)电连接。
优选的,所述试样盒的内腔尺寸满足长与厚的比值大于8-10,宽与厚度的比值大于8-10。
本实用新型与现有技术相比所取得的有益效果为:利用超级恒温水浴作为热源,与试样盒底部的黄铜板形成恒温边界,温度均匀性和稳定性好,且易于实现,达到了同时测试固体材料导热系数及热扩散率的目的,且测得的结果可靠、精度高。
附图说明
1.图1为本实用新型的硬件结构原理图;
2.图2为本实用新型计算原理模型;
图1和图2中,1、水介质 2、加热元件 3、温度传感器 4、温控单元 5、搅拌器 6、绝热水浴箱 7、底部黄铜板 8、侧绝热板 9、上绝热板 10、试样 11、热电偶 12、信号采集与传输单元 13、微型计算机。
3.图3和图4为本发明待估参数灵敏度系数计算曲线图。
具体实施方式
下面结合附图和具体实施方式对本实用新型做进一步详细的说明。
如图1所示,水介质1温度由超级恒温水浴系统控制,超级恒温水浴系统包括加热元件2、温度传感器3、温控单元4、搅拌器5、绝热水浴箱6等组成,其中加热元件2用于加热水介质1,温度传感器3用于检测水介质1的实时温度,并反馈给温控单元4,控制加热元件2的启闭,以保证水介质1维持恒定温度。搅拌器5用以搅拌水介质1,以使其温度均匀,绝热层6用于使水介质1与周围环境绝热,降低热量散发因素的影响。试样盒由底部黄铜板7、侧绝热板8和上绝热板9组成,当其底部黄铜板7与恒温的水介质1接触时,将形成恒温边界。如在某一时刻τ突然将厚度为δ,初始温度为t0的试样10与黄铜板7接触,试样10中将被加热,而产生恒温边界条件下的一维热传导过程。当试样盒的内腔尺寸满足长和宽与厚度的比值大于8-10时,能更好地形成一维非稳态传热条件。试样10被加热过程中,上绝热面处的温度变化通过热电偶11检测,并通过信号采集与传输单元12传送给微型计算机13,进而利用数学模型计算得 到试样热物性参数。
如图2所示,设试样导热系数为λ,密度ρ、比热Cp等均为常数,其初始温度为t0,坐标原点取于试样上绝热面,则试样10与恒温T=tw的黄铜板7接触后,其内的非稳态导热微分方程及其初始、边界条件如下:
τ=0,t=t0 (2a)
x=0,
x=δ,t=tw (2c)
式(1)的解析解如下:
式中θ(x,τ)为过余温度,θ=t-tw;热扩散率a=λ/(ρCp)。
式(3)为无穷级数,无法直接求解。而由于恒温边界条件等价于毕渥数Bi→∞时的热对流边界条件,从而式(3)可改写为:
式中βn为超越方程ctgβ=β/Bi的系列解,对应不同的Bi数,β值也不同(具体可查表),而在Bi→∞时β1=1.5708,βn=β1+(n-1)×π(n=1,2,…),试算表明, 针对固体材料热物性测试需求取前15项的计算误差即可控制在1%以内,显然在试样10的密度ρ、导热系数λ、比热容Cp及厚度δ等物理参数已知的情况下,可计算得到试样10内任意位置x任意时刻τ的温度分布情况,此即恒温边界条件下的非稳态传热数学模型。
在式(4)的基础上,且试样热物性参数未知的情况下,由测得的x=0处的温度数据进行反演推算,也称参数估计。如前所述,从某一时刻开始试样下表面被温度为tw的热板加热,该温度扰动为输入信号,记为u(τ),试样内任意一点处的温升记为Y(τ,η),其不但与试样导热系数、比热容等热物性参数有关,还取决于热板温度tw、时间、空间位置等,所有这些参数构成一个向量:
η=(η1,η2,η3,......,ηm)T (5)
其中m为参数个数。
此处试样导热系数λ、比热容Cp未知,密度ρ可实测获得,首先获得x=0表面处k(k>2)个离散测点处的温度测量值Yi(τ,η)(i=1…k),然后基于估计的参数值,通过式(4)计算各离散测点的温度值θi(τ,η)(i=1…k),并进行对比:
由于测量误差以及参数估计误差的存在,ε(η)>0,估计试样导热系数及热扩散率的反问题可归结为对式(4)进行迭代计算,最终使目标函数ε(η)→min。
参数向量η中的j个参数要能同时以足够的精度估计出来,必须满足在最小二乘估计意义下,测量时间范围内(参数估计的时间区间内),参数的灵敏度线性无关。灵敏度系数反映了参数βj的微小扰动对输出的影响。
式(4)中待估计参数为:η1=λ,η2=Cp
为检验对不同导热系数大小的材料参数估计可行性,分别选取导热系数较大的玻璃和导热系数较小的松散煤体进行参数灵敏度分析。设玻璃物理参数为:厚度5mm,密度2530kg/m3,导热系数0.701W/m.℃,比热容750J/kg.℃。松散煤体物理参数为:密度1010kg/m3,试样厚度60mm,粒度范围0-6mm,导热系数0.14W/m.℃,比热容840J/kg.℃。
灵敏度系数计算结果分别如图3和图4所示。从图3和图4可以看出,玻璃和松散煤体的导热系数和比热容参数均灵敏度线性无关,但比热容Cp灵敏度极小(仅10e-3数量级),意味着直接对其进行估计误差将比较大。实际上在测得温度θ(x,τ)的情况下,式(4)中的λ/(ρCp)项即热扩散率a,将只有此一个未知数,可首先将热扩散率a求解出来,然后利用其对估计的比热容Cp进行修正。也即,本系统首先采用随机共轭梯度法估计热扩散率,然后再对比热容估计结果进行修正,计算过程如下:
1)估计热扩散率a
(1)根据先验信息确定待反演参数的初始猜测值η0(此时参数向量η=(a)T)、设定一较小的大于零的数ζ作为退出控制误差。
(2)修改模型中的参数向量ηj:
式中ηj new和ηj old分别为修改前后的第j个参数;dηj表示第j个参数的基本偏离值,通过手工设定,用于控制模型参数的修改量;sign为符号函数,表示在相同概率下(0.5的概率)取值为1和-1;Rand()产生[0,1]范围内的随机数。
(3)将参数修改前后的目标函数分别记为εold(η)和εnew(η),其差值记为: △ε=εold(η)-εnew(η),当△ε<0时,p=1;△ε>0时,p=0。
(4)当p=1时,接受修改后的模型,计算共轭系数、梯度、搜索方向及搜索步长等参量;共轭梯度法收敛后,保存反演估计值
(5)充分步骤(2)、(3)和(4),比较每次接受的反演估计值,直到收敛到允许误差范围。
2)将待估参数向量更改为η=(λ,Cp)T,然后按照步骤(1)至(5)估计导热系数λ及比热容Cp。
3)根据热扩散率a及导热系数λ估计值修正比热容Cp估计值。
理论上每个特定时刻τ均可得到一组数值λ和Cp,对计算出的一系列值求加权平均作为最终参数估计值,同时可消除测试过程中因干扰因素造成的个别数据波动影响。
以上实例所假设的参数和尺寸都是帮助理解本实用新型的方法和核心思想,针对不同的固体材料应该查相关手册来确定其参数以及根据实际情况来测量其尺寸大小。综上所述,本说明书参数和尺寸不应理解为对本实用新型的限制,在不脱离本实用新型的精神和范围情况下,还可以做出各种变化,因此所有等同的技术方案也属于本实用新型的范畴。
Claims (2)
1.一种可同时测固体材料导热系数及热扩散率的系统,其特征在于:包括超级恒温水浴系统、试样盒和温度采集与分析系统,所述超级恒温水浴系统包括水介质(1)、加热元件(2)、温度传感器(3)、温控单元(4)、搅拌器(5)、绝热水浴箱(6),所述水介质(1)装在所述绝热水浴箱(6)中,所述加热元件(2)安装在所述绝热水浴箱(6)的内侧下方,与所述温控单元(4)电连接,所述温控单元(4)安装在所述绝热水浴箱(6)的顶部,所述温控单元(4)输出端与所述温度传感器(3)电连接,所述温度传感器(3)输出端进入所述绝热水浴箱(6),所述搅拌器(5)向内安装在所述绝热水浴箱(6)的底部;所述试样盒安装在所述绝热水浴箱(6)的顶部,由底部黄铜板(7)、侧绝热板(8)和上绝热板(9)组成,所述温度采集与分析系统包括热电偶(11)、信号采集与传输单元(12)和微型计算机(13),所述热电偶11一端通过所述上绝热板(9)进入所述试样盒,另一端通过所述采集与传输单元(12)与微型计算机(13)电连接。
2.根据权利要求1所述的可同时测固体材料导热系数及热扩散率的系统,其特征在于:所述试样盒的内腔尺寸满足长与厚的比值大于8-10,宽与厚度的比值大于8-10。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201420475647.5U CN204086184U (zh) | 2014-08-21 | 2014-08-21 | 一种可同时测固体材料导热系数及热扩散率的系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201420475647.5U CN204086184U (zh) | 2014-08-21 | 2014-08-21 | 一种可同时测固体材料导热系数及热扩散率的系统 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN204086184U true CN204086184U (zh) | 2015-01-07 |
Family
ID=52178885
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201420475647.5U Expired - Fee Related CN204086184U (zh) | 2014-08-21 | 2014-08-21 | 一种可同时测固体材料导热系数及热扩散率的系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN204086184U (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106198354A (zh) * | 2016-08-19 | 2016-12-07 | 中国华电科工集团有限公司 | 一种渗流、应力、温度耦合试验机 |
CN113138207A (zh) * | 2021-04-22 | 2021-07-20 | 安徽理工大学 | 一种正交各向异性固体材料热扩散系数测试系统及方法 |
-
2014
- 2014-08-21 CN CN201420475647.5U patent/CN204086184U/zh not_active Expired - Fee Related
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106198354A (zh) * | 2016-08-19 | 2016-12-07 | 中国华电科工集团有限公司 | 一种渗流、应力、温度耦合试验机 |
CN106198354B (zh) * | 2016-08-19 | 2023-12-05 | 中国华电科工集团有限公司 | 一种渗流、应力、温度耦合试验机 |
CN113138207A (zh) * | 2021-04-22 | 2021-07-20 | 安徽理工大学 | 一种正交各向异性固体材料热扩散系数测试系统及方法 |
WO2022222776A1 (zh) * | 2021-04-22 | 2022-10-27 | 安徽理工大学 | 一种正交各向异性固体材料热扩散系数测试系统及方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104215660B (zh) | 一种可同时测固体材料导热系数及热扩散率的方法及系统 | |
CN103411996B (zh) | 固体材料导热系数测量装置及测量方法 | |
CN101126729B (zh) | 双热流计稳态法材料热导率测量方法 | |
CN105548606B (zh) | 基于mems的柔性流速传感器的流速测量方法 | |
CN109752113B (zh) | 薄片温度传感器及其应用中位置确定方法和电路设计方法 | |
CN101113963A (zh) | 一种测量液体导热系数的方法及其装置 | |
EP2126533A1 (en) | Apparatuses and methods for measuring and controlling thermal insulation | |
CN106338323B (zh) | 一种流量计现场检定方法及检定系统 | |
CN103954650B (zh) | 一种固体材料热扩散系数测试方法与系统 | |
CN109443601A (zh) | 基于热电阻的高精度多点温度测量系统及温度校准方法 | |
CN107037079B (zh) | 一种支撑梁式mems复合传感器及其制备和测试方法 | |
CN102879421A (zh) | 导热系数的测量装置及其测量方法 | |
CN204086184U (zh) | 一种可同时测固体材料导热系数及热扩散率的系统 | |
CN109738092B (zh) | 一种双向热电堆式薄膜热流计及热流测量方法 | |
CN101871900A (zh) | 一类用于导热系数测量的传感器 | |
CN203502367U (zh) | 一种瞬态平面热源法测试材料导热系数的装置 | |
CN103257154B (zh) | 大温差样品的半球向全发射率与导热系数的测量方法 | |
CN108051475A (zh) | 一种对流换热系数的快速测量方法 | |
CN108956688B (zh) | 一种建筑材料表面综合辐射吸收系数的测量系统及方法 | |
CN107966472B (zh) | 一种高温接触热阻的无损快速测量方法 | |
CN204128607U (zh) | 一种热膜式流量传感芯片 | |
CN101285786B (zh) | 谐波探测技术用于微槽道局部对流换热系数测定的方法 | |
CN108776152A (zh) | 一种建筑材料表面对流热交换系数的测量系统及方法 | |
Cortellessa et al. | Experimental and numerical analysis in heat flow sensors calibration | |
CN103245692A (zh) | 基于稳态分析的半球向全发射率与导热系数的测量方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
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: 20150107 Termination date: 20150821 |
|
EXPY | Termination of patent right or utility model |