1. 引言
随着人们对浴室用品(如毛巾,浴巾等)的卫生程度要求越来越高,而浴室内常年潮湿且阴暗,因此长时间放置在浴室里潮湿的毛巾,浴巾以及毛巾架本身等很容易滋生细菌和螨虫 [1] [2] [3] 。干燥是保持毛巾等织物洁净的简单有效方式之一 [4] [5] ,如何对晾晒在浴室内潮湿的毛巾和浴巾进行烘干越来越被人们所重视。电热毛巾架随之进入人们视野,电热毛巾架是电暖的一种延伸,采用整架一体发热,能够快速烘干毛巾,高效灭杀99.99%的霉菌滋生 [6] [7] 。其实,电热毛巾架一开始属于高档卫浴用品的配套产品 [8] ,由于其出色的干燥及抑菌性能,逐渐进入普通家庭的卫浴用品中。
目前关于电热毛巾架的仿真分析比较少 [9] ,佛山科学技术学院工业设计与陶瓷艺术学院的李天赠 [10] 等提出了一种基于数值模拟技术的产品概念设计方法,阐述了如何将概念设计方法应用在电热毛巾架的设计当中。南通航运职业技术学院的杜虹 [11] 提出一种基于Arduino的电热毛巾架实现了恒温控制等技术。
本文通过对合作单位主流的电热毛巾架的结构进行提取,建立其主要结构的三维模型,通过ANSYS Fluent数值模拟软件对电热毛巾架的发热以及干燥过程进行分析,根据分析结果对毛巾架结构进行优化,最后根据优化后的结构再次进行仿真分析,通过优化前后干燥效果的对比,证明了结构优化的合理性。
2. 毛巾架模型的建立及传热问题
2.1. 电热毛巾架模型的建立
根据合作单位主流的电热毛巾架晾晒方式,并对结构特征进行提取仅保留电热毛巾架的发热杆件以及毛巾模型,得到两种类型的电热毛巾架干燥方式:一是毛巾直接与电热毛巾架发热杆接触,如图1所示,即毛巾悬挂在电热毛巾架发热杆件上;二是毛巾不与电热毛巾发热杆接触晾晒,如图2所示,即毛巾悬挂在与电热毛巾架距离一定的位置上。
由于电热毛巾架的安装是固定在墙面且与地面垂直,并由两种电热毛巾架的干燥方式可以知道,接触式的电热毛巾架由于在干燥过程中毛巾与发热杆件始终接触,杆件产生的热量大部分会被毛巾吸收,因此接触式电热毛巾架在其结构上相对合理。但由于发热杆件始终被毛巾覆盖,产生的热量会被聚集起来,在使用过程中会产生一定的风险。而分离式电热毛巾架由于发热杆件和毛巾具有一定的距离,热量不会聚集在毛巾内部从而降低了风险,但是由于折叠后的毛巾只有一半的面积可以直接接收到发热杆件产生的热量,因此分离式的电热毛巾架在干燥速率以及热量的利用率上大打折扣。
基于上述分析,本文对分离式电热毛巾架进行模拟分析。对电热毛巾架的结构进行提取,建立电热毛巾架、毛巾和壁面的三维模型如图3所示。
![](//html.hanspub.org/file/2-2990434x8_hanspub.png?20230609095519011)
Figure 1. Contact type electric towel rack
图1. 接触式电热毛巾架
![](//html.hanspub.org/file/2-2990434x9_hanspub.png?20230609095519011)
Figure 2. Separate electric towel rack
图2. 分离式电热毛巾架
![](//html.hanspub.org/file/2-2990434x10_hanspub.png?20230609095519011)
Figure 3. Sketch of structural analysis of electric towel rack
图3. 电热毛巾架结构分析简图
2.2. 电热毛巾架温度场分析的数学基础
电热毛巾架的传热问题可通过计算流体的方法进行求解,本文主要通过ANSYS Fluent作为数值求解平台对电热毛巾架的传热问题进行求解。将其传热过程视为自然对流模型,求解过程主要设计三大方程即能量方程、粘性方程以及辐射模型。
能量方程的求解原理如下:
(1)
式中:
为有效导热率;
为组分的j扩散通量;
包含化学反应放(吸)热以及任何其它的由用户定义的体积热源。
其中
(2)
式中:显焓
的定义(对理想气体)为
。
由于电热毛巾架的传热过程相对简单,可以将其看作为自然对流过程,因此将求解中选择层流模型。对于辐射模型的在本文中选择DO模型进行求解,其求解原理如式(3)的求解如下:
(3)
式中:
为位置向量;
为方向向量;
为散射方向;s为沿程长度;a为吸收系数;n为折射率;
为散射系数;I为辐射强度;
相位函数;
为空间立体角。
3. 热毛巾架的温度场仿真
3.1. 原电热毛巾架的温度场仿真
基于以上的数值分析原理,采用ANSYS Fluent对电热毛巾架干燥毛巾的过程进行模拟分析。电热毛巾架材料设置为铝合金,其发热功率设置为30,000 W/m3,密度及热导率等相关参数按照Fluent中铝的默认参数进行设置。毛巾相关参数的设置对分析过程以及结果的收敛程度比较重要,因此毛巾的模型密度、比热以及热导率等数值由材料决定,本文所选择毛巾为棉质毛巾由相关资料得到其主要热力学参数如表1。相关参数设置完成后,对毛巾架及毛巾进行网格划分,考虑到网格质量及数量对求解时长和精度等影响,将计算域的网格尺寸控制在50 mm,而接触面的网格进行加密设置为3 mm。最终整个计算模型的网格单元数量为658,114,节点数为898,672,其中电热毛巾架及毛巾模型的网格划分如图4。
![](Images/Table_Tmp.jpg)
Table 1. Towel model parameter settings
表1. 毛巾模型参数设置
![](//html.hanspub.org/file/2-2990434x25_hanspub.png?20230609095519011)
Figure 4. Grid division of electric towel rack
图4. 电热毛巾架网格划分
根据仿真分析结果,得到毛巾干燥过程温度分布云图,如图5所示;电热毛巾架烘干过程横截面温度分布云图,如图5(c)所示;毛巾温度分布区域占比图,如图6所示。从图5可以看出,由于毛巾折叠为均匀的两层放置在电热毛巾架一侧,由于空气被发热杆件加热后密度减小,被加热后的热气向上运动导致毛巾受热不均匀,最高温度在42℃左右,而毛巾底端温度只有37℃左右,且由于重力的作用,水分将汇集到毛巾最低端,这将导致毛巾最低端干燥时间大大增加。而背对电热毛巾架的一面,电热毛巾架产生的热量基本靠近电热毛巾架的一面吸收,其干燥过程中温度只能维持在36℃左右。
通过图5(c)截面温度分布云图可以看出,电热毛巾架产生的热量在放置毛巾的一侧,其热量可以很好的被毛巾吸收,而靠近墙面的一侧,其产生的热量基本传递到了空气当中,没有被充分利用。从图6可以看出,毛巾在干燥过程中其温度分布的百分比基本分布在两个温度区间,这与图5所分析的结果一致,即毛巾折叠放置在电热毛巾架的一侧进行干燥,一侧可以得到很好的干燥效果,而另一侧则不理想。
(a) 靠近毛巾架的一侧 (b) 靠近毛巾架的一侧 (c) 轴截面温度分布
Figure 5. Temperature distribution cloud of towel drying process
图5. 毛巾干燥过程温度分布云图
![](//html.hanspub.org/file/2-2990434x27_hanspub.png?20230609095519011)
Figure 6. Towel temperature distribution area share chart
图6. 毛巾温度分布区域占比图
3.2. 优化后电热毛巾架的温度场仿真
通过上述分析,上述电热毛巾架干燥毛巾的效果不理想,产生的大部分热量被浪费,不能很好地被毛巾吸收。因此根据上述分析结果,在不改变电热毛巾架单侧干燥毛巾的前提下,改进电热毛巾架结构,优化后的电热毛巾架结构,如图7所示。
![](//html.hanspub.org/file/2-2990434x28_hanspub.png?20230609095519011)
Figure 7. Sketch of the structure of the optimized electric towel rack
图7. 优化后电热毛巾架结构简图
根据3.1节数值模拟分析参数设置,对优化后的电热毛巾架干燥毛巾的过程再进行一次数值模拟分析。得到优化后毛巾温度分布云图,如图8所示;电热毛巾架横截面温度分布云图,如图8(c)所示;毛巾温度所占整个毛巾的百分比,如图9所示。通过图8可以看出,优化后的电热毛巾架干燥毛巾的效果相对比较理想,靠近电热毛巾架的一侧最高温度基本维持在43℃左右,毛巾最低端温度也能达到41℃左右,这将大大提高毛巾最低端的干燥速率,而另一侧的温度也维持在了较高的38℃附近,提高了毛巾的整体干燥速度。
通过优化后的电热毛巾架结构可以看出,由于底端的3~4根加热杆件向外延伸,而毛巾折叠后最低端位于加热杆件正上方,加快了毛巾最低端的干燥速度,同时当毛巾底部周围空气被加热后由于膨胀向上运动,干燥且高温的空气进入毛巾两夹层之间,充分利用了电热毛巾架的热效率。
(a) 靠近毛巾架的一侧 (b) 靠近毛巾架的一侧 (c) 轴截面温度分布
Figure 8. Cloud map of towel temperature distribution after optimization
图8. 优化后毛巾温度分布云图
![](//html.hanspub.org/file/2-2990434x30_hanspub.png?20230609095519011)
Figure 9. Percentage distribution of towel temperature after optimization
图9. 优化后毛巾温度分布百分比
从图9优化后毛巾温度所占百分比的直方图可以看出,在双层折叠毛巾的前提下,毛巾的温度分布主要集中在40℃和43℃左右,这是与优化前的毛巾温度分布在两个主要温度区间一致,但是相比优化前的数值模拟结果有所改善。背对电热毛巾架的毛巾一面温度由原来的36℃左右提升到了38℃左右,而另一面则有原来的42℃提升到了43℃左右,而毛巾最低端温度提高了4℃左右,毛巾整体温度得到了提升,加快了毛巾的干燥速度。
4. 结论
通过对结构优化前后的电热毛巾架数值模拟分析结果对比,发现毛巾架的设计结构对电热毛巾架的能量利用率影有一定的影响。在保持输入功率时间不变的条件下,优化后的毛巾温度比初始结构下的毛巾温度整体提升了3℃左右,同时大大提高了最难干燥的毛巾低端的温度,提高了毛巾干燥速率,增加用户体验的,同时也减少了能源的浪费。证明了结构优化设计在电热毛巾架热效率提升上的可行性,为后来者提供一份理论参考,同时也为我国在2060年前实现碳中和贡献一份力量。
基金项目
凤阳县科技计划项目:高效智能电热毛巾架研制与产业化开发(GY2022-05)。
NOTES
*通讯作者。