CN116110525B - 一种纸尿裤内部微气候热湿环境计算方法和系统 - Google Patents
一种纸尿裤内部微气候热湿环境计算方法和系统 Download PDFInfo
- Publication number
- CN116110525B CN116110525B CN202310393847.XA CN202310393847A CN116110525B CN 116110525 B CN116110525 B CN 116110525B CN 202310393847 A CN202310393847 A CN 202310393847A CN 116110525 B CN116110525 B CN 116110525B
- Authority
- CN
- China
- Prior art keywords
- paper diaper
- grid
- diaper
- heat
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 89
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 134
- 238000004364 calculation method Methods 0.000 claims abstract description 49
- 230000027939 micturition Effects 0.000 claims abstract description 19
- 239000010410 layer Substances 0.000 claims description 66
- 238000012546 transfer Methods 0.000 claims description 38
- 230000008569 process Effects 0.000 claims description 26
- 210000002700 urine Anatomy 0.000 claims description 26
- 239000012792 core layer Substances 0.000 claims description 17
- 238000001704 evaporation Methods 0.000 claims description 16
- 230000008020 evaporation Effects 0.000 claims description 16
- 230000008859 change Effects 0.000 claims description 15
- 239000000463 material Substances 0.000 claims description 14
- 238000010521 absorption reaction Methods 0.000 claims description 12
- 230000005540 biological transmission Effects 0.000 claims description 10
- 239000007788 liquid Substances 0.000 claims description 9
- 239000002344 surface layer Substances 0.000 claims description 9
- 238000009834 vaporization Methods 0.000 claims description 8
- 230000008016 vaporization Effects 0.000 claims description 8
- 238000010276 construction Methods 0.000 claims description 4
- 238000012545 processing Methods 0.000 claims description 4
- 229920006395 saturated elastomer Polymers 0.000 claims description 4
- 238000012887 quadratic function Methods 0.000 claims description 3
- 230000036541 health Effects 0.000 abstract description 5
- 230000036559 skin health Effects 0.000 description 6
- 238000013461 design Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000012356 Product development Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 238000012938 design process Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 239000002657 fibrous material Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 239000012782 phase change material Substances 0.000 description 1
- 230000036620 skin dryness Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C60/00—Computational materials science, i.e. ICT specially adapted for investigating the physical or chemical properties of materials or phenomena associated with their design, synthesis, processing, characterisation or utilisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computing Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Computational Biology (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Investigating Or Analyzing Materials Using Thermal Means (AREA)
- Absorbent Articles And Supports Therefor (AREA)
Abstract
本发明提供一种纸尿裤内部微气候热湿环境计算方法和系统,涉及婴幼儿健康护理技术领域,该方法包括获取纸尿裤的热力学参数,建立纸尿裤内部微气候模型并对纸尿裤部分进行网格划分;对建立的纸尿裤内部微气候模型进行初始化;利用有限体积法计算在一定排尿量的条件下,不同时刻纸尿裤每个网格的温度和含水量;根据计算结果获得皮肤损失的热量、不同时刻纸尿裤的温度分布和湿度分布,完成纸尿裤内部热湿环境的计算;本发明基于有限体积法进行计算,运算精度高、计算周期短,能快速精准地反映纸尿裤内微气候复杂的温湿度场及流场的分布情况,能够系统高效地预测纸尿裤微气候环境的湿热舒适性。
Description
技术领域
本发明涉及婴幼儿健康护理技术领域,更具体地,涉及一种纸尿裤内部微气候热湿环境计算方法和系统。
背景技术
纸尿裤对婴幼儿穿着舒适性和皮肤健康具有重要影响。纸尿裤通过多层结构有机组合实现以下功能:包括吸收尿液、保持透气性、防止尿液侧漏回渗等。在不同环境条件下,纸尿裤需要确保婴幼儿保排尿后,皮肤附近的微气候区域维持湿热平衡,处于一个让皮肤舒适的温湿度范围。
随着人们对健康护理理念的深刻认识,健康舒适纸尿裤产品的需求日益增加。为了把微气候维持在舒适区,纸尿裤产品的选择需要充分考虑婴幼儿的尿液量、日用夜用、冬季夏季和环境潮湿程度等因素。因此,有必要开发出适应各种要求的健康舒适产品。
然而,精准开发适应各种需求的健康舒适纸尿裤的产品有一定难度。首先,纸尿裤微气候环境中的湿热传递过程十分复杂,包括了液态水输运和吸收、湿气传递、多种形式的热传递以及湿热耦合作用。此外,排尿后的芯体吸湿和放湿过程、婴幼儿活动等因素大大增加了湿热舒适性预测的难度。同时,纸尿裤针对结构设计的性能要求也有一定矛盾性,例如增加厚度可以提高吸水量或促进皮肤干燥,但也会大幅增加热阻而影响热舒适性。目前的产品开发方法主要处于实验试错、评价测试和用户反馈阶段,缺乏深入系统的预测模型。
目前的现有技术公开了一种多尺度框架下的热物理传输数值仿真方法,提出了一组多尺度非线性模型用来描述人体,织物,纤维材料和相变材料颗粒之间的混合型耦合热湿传输;所述多尺度模型整合了上述实体之间的动态热边界条件,采用了耦合偏微分方程描述,所述偏微分方程通过有限体积方法进行离散化;同时,考虑到不同的穿戴情景,提出了智能保暖服装数字化仿真的方案;现有技术中的方法主要利用在服装设计领域,能够在智能保暖服装的工程设计过程中提前评估服装的热智能质量,对于湿热传递过程十分复杂的个性化纸尿裤设计领域,现有技术存在空间分辨率低、设计效率较低和热湿环境预测精度不高的缺陷。
发明内容
本发明为克服上述现有技术在针对个性化纸尿裤设计时效率较低且内部热湿环境计算精度不高的缺陷,提供一种纸尿裤内部微气候热湿环境计算方法和系统,能够系统高效地预测纸尿裤微气候环境的湿热舒适性,揭示不同结构材料对湿热传递性能的影响,对个性化产品快速优化设计和保护婴幼儿皮肤健康,具有实际应用价值和重要现实意义。
为解决上述技术问题,本发明的技术方案如下:
一种纸尿裤内部微气候热湿环境计算方法,包括以下步骤:
S1:获取纸尿裤的热力学参数,建立纸尿裤内部微气候模型,具体为:
所述纸尿裤内部微气候模型包括从内到外的三部分,具体为皮肤、纸尿裤和外界环境;
所述纸尿裤部分包括从内到外的面层、覆盖层、吸水芯体层和底膜层;
所述纸尿裤各层的厚度量级为毫米级;
所述皮肤和外界环境均为无厚度边界,皮肤边界为恒定温度;
基于以下假设:
假设一:纸尿裤内起主要锁水作用的吸水芯体层保水性强,吸水芯体层内部水分在短时间内传递缓慢,因此忽略纸尿裤内部的对流换热,将传热过程视为导热过程;
假设二:尿液在纸尿裤内扩散的区域有限,尿液所在区域的换热量远大于其他区域,因此将计算区域设置为尿液存在的区域;
假设三:纸尿裤内部各层之间存在间隙,在尿液存在的区域内含水量很高,因此本模型假设纸尿裤各层间隙内的空气湿度一直处于该位置当前温度下的饱和湿度,同时忽略空气的对流传热和纸尿裤侧边界的传热;
假设四:尿液存在的区域所有边界均存在蒸发相变过程和热量传递过程,对于同种材料来说侧边界水分的流失速度大于中心区域,因此本模型考虑水分质量流失速度在皮肤切向上的差异;
对纸尿裤内部微气候模型中的纸尿裤部分进行网格划分;
S2:对建立的纸尿裤内部微气候模型进行初始化;
所述初始化包括设置含水量关系拟合函数、相变类型、初始条件和边界条件;
所述含水量关系拟合函数包括对面层、覆盖层、吸水芯体层和底膜层均设置导热系数-含水量关系拟合函数、密度-含水量拟合函数、比热容-含水量拟合函数、水分传输比例-时间关系拟合函数和水分传输比例-水平分布拟合函数;
设置纸尿裤内部微气候模型的初始条件具体为:
设置纸尿裤内部微气候模型各部分的初始温度、纸尿裤的相对湿度、外界环境的初始温度和相对湿度,以及排尿次数;
设置纸尿裤内部微气候模型的边界条件具体为:
设置皮肤表面的边界条件为预设温度的恒温;
设置纸尿裤的边界条件为纸尿裤最外侧与外界环境的热交换过程包括蒸发和对流;
设置对流换热系数和水蒸气的蒸发速率;
S3:利用有限体积法计算在一定排尿量的条件下,不同时刻纸尿裤每个网格的温度和含水量,具体方法为:
所述纸尿裤的能量变化包括皮肤输入的热量、失去水分的内能、纸尿裤最外层对流换热和液态水汽化潜热;
纸尿裤每个网格的能量控制方程为
其中,为本时间步密度,/>为上一时间步密度,/>为本时间步比热容,/>为上一时间步比热容,/>为上一时间步导热系数,/>为网格宽度,/>为网格高度,/>为网格中心与右侧相邻网格中心距离,/>为网格中心与左侧相邻网格中心距离,/>为网格中心与上侧相邻网格中心距离,/>为网格中心与下侧相邻网格中心距离,/>为本时间步网格中心温度,/>为上一时间步网格中心温度,/>为右侧相邻网格中心点温度,/>为左侧相邻网格中心点温度,/>为上侧相邻网格中心点温度,/>为下侧相邻网格中心点温度;
、/>、/>、/>、/>和/>的系数分别为:
其中,为上一时间步温度系数,/>为本时间步温度系数,/>为右侧网格温度系数,/>为左侧网格温度系数,/>为上侧网格温度系数,/>为下侧网格温度系数;
根据每个网格的能量控制方程计算不同时刻纸尿裤各层的能量变化以及每个网格的温度;
纸尿裤每个网格的含水量计算方程为:
其中,为纸尿裤第i层的水分传递比例,/>为水平方向上的分配比例,为时间当前时间下的水分蒸发量,/>为新时间步网格含水量,为上一时间步网格含水量;
纸尿裤每个网格的总含水量为各层含水量之和;
根据每个网格的含水量计算方程计算不同时刻纸尿裤各层的含水量,根据每个网格不同时刻纸尿裤各层的含水量计算每个网格不同时刻纸尿裤的总含水量;
S4:根据不同时刻纸尿裤每个网格的温度和含水量获得皮肤损失的热量、不同时刻纸尿裤的温度分布和湿度分布,完成纸尿裤内部热湿环境的计算。
优选地,所述步骤S1中所述的纸尿裤的热力学参数包括:纸尿裤各层材料的导热系数、密度和孔隙率。
优选地,所述纸尿裤各层含水量的切向分布呈二次函数分布,在法向方向上具有线性关系。
优选地,所述皮肤输入的热量具体为:
其中,为皮肤温度,/>为皮肤相邻网格温度,/>为导热面积,具体为单个网格的面积,/>为单个网格的导热系数。
优选地,所述纸尿裤最外层对流换热具体为:
其中,h为对流换热系数,表示对流强度,为外界环境温度,/>为时间步长。
优选地,所述液态水汽化潜热具体为:
其中,为水蒸气质量。
本发明还提供一种纸尿裤内部微气候热湿环境计算系统,应用上述的一种纸尿裤内部微气候热湿环境计算方法,包括:
模型构建和网格划分单元:用于获取纸尿裤的热力学参数,建立纸尿裤内部微气候模型;
所述纸尿裤内部微气候模型包括从内到外的三部分,具体为皮肤、纸尿裤和外界环境;
并对纸尿裤内部微气候模型中的纸尿裤部分进行网格划分;
初始化单元:用于对建立的纸尿裤内部微气候模型进行初始化;
计算单元:用于利用有限体积法计算在一定排尿量的条件下,不同时刻纸尿裤每个网格的温度和含水量;
计算结果处理单元:用于根据不同时刻纸尿裤每个网格的温度和含水量获得皮肤损失的热量、不同时刻纸尿裤的温度分布和湿度分布,完成纸尿裤内部热湿环境的计算。
与现有技术相比,本发明技术方案的有益效果是:
本发明提供一种纸尿裤内部微气候热湿环境计算方法和系统,该方法通过获取纸尿裤的热力学参数,建立纸尿裤内部微气候模型并对纸尿裤部分进行网格划分;对建立的纸尿裤内部微气候模型进行初始化;利用有限体积法计算在一定排尿量的条件下,不同时刻纸尿裤每个网格的温度和含水量;根据不同时刻纸尿裤每个网格的温度和含水量获得皮肤损失的热量、不同时刻纸尿裤的温度分布和湿度分布,完成纸尿裤内部热湿环境的计算;
本发明基于有限体积法进行计算,运算精度高、计算周期短,能快速精准地反映纸尿裤内微气候复杂的温湿度场及流场的分布情况,能够系统高效地预测纸尿裤微气候环境的湿热舒适性,揭示不同结构材料对湿热传递性能的影响,对个性化产品快速优化设计和保护婴幼儿皮肤健康,具有实际应用价值和重要现实意义。
附图说明
图1为实施例1所提供的一种纸尿裤内部微气候热湿环境计算方法流程图。
图2为实施例2所提供的纸尿裤内部微气候模型结构图。
图3为实施例2所提供的纸尿裤各层含水量的分布图。
图4为实施例2所提供的能量变化示意图。
图5为实施例2所提供的一种纸尿裤内部微气候热湿环境计算方法计算结果图。
图6为实施例3所提供的一种纸尿裤内部微气候热湿环境计算系统结构图。
具体实施方式
附图仅用于示例性说明,不能理解为对本专利的限制;
为了更好说明本实施例,附图某些部件会有省略、放大或缩小,并不代表实际产品的尺寸;
对于本领域技术人员来说,附图中某些公知结构及其说明可能省略是可以理解的。
下面结合附图和实施例对本发明的技术方案做进一步的说明。
实施例1
如图1所示,本实施例提供一种纸尿裤内部微气候热湿环境计算方法,包括以下步骤:
S1:获取纸尿裤的热力学参数,建立纸尿裤内部微气候模型;
所述纸尿裤内部微气候模型包括从内到外的三部分,具体为皮肤、纸尿裤和外界环境;
对纸尿裤内部微气候模型中的纸尿裤部分进行网格划分;
S2:对建立的纸尿裤内部微气候模型进行初始化;
S3:利用有限体积法计算在一定排尿量的条件下,不同时刻纸尿裤每个网格的温度和含水量;
S4:根据不同时刻纸尿裤每个网格的温度和含水量获得皮肤损失的热量、不同时刻纸尿裤的温度分布和湿度分布,完成纸尿裤内部热湿环境的计算。
在具体实施过程中,首先获取纸尿裤的热力学参数,建立纸尿裤内部微气候模型并对纸尿裤部分进行网格划分;对建立的纸尿裤内部微气候模型进行初始化;利用有限体积法计算在一定排尿量的条件下,不同时刻纸尿裤每个网格的温度和含水量;根据不同时刻纸尿裤每个网格的温度和含水量获得皮肤损失的热量、不同时刻纸尿裤的温度分布和湿度分布,完成纸尿裤内部热湿环境的计算;
本发明基于有限体积法进行计算,运算精度高、计算周期短,能快速精准地反映纸尿裤内微气候复杂的温湿度场及流场的分布情况,能够系统高效地预测纸尿裤微气候环境的湿热舒适性,揭示不同结构材料对湿热传递性能的影响,对个性化产品快速优化设计和保护婴幼儿皮肤健康,具有实际应用价值和重要现实意义。
实施例2
本实施例提供一种纸尿裤内部微气候热湿环境计算方法,包括以下步骤:
S1:获取纸尿裤的热力学参数,建立纸尿裤内部微气候模型;
所述纸尿裤内部微气候模型包括从内到外的三部分,具体为皮肤、纸尿裤和外界环境;
对纸尿裤内部微气候模型中的纸尿裤部分进行网格划分;
S2:对建立的纸尿裤内部微气候模型进行初始化;
S3:利用有限体积法计算在一定排尿量的条件下,不同时刻纸尿裤每个网格的温度和含水量;
S4:根据不同时刻纸尿裤每个网格的温度和含水量获得皮肤损失的热量、不同时刻纸尿裤的温度分布和湿度分布,完成纸尿裤内部热湿环境的计算。
所述步骤S1中所述的纸尿裤的热力学参数包括:纸尿裤各层材料的导热系数、密度和孔隙率;
所述步骤S1中所建立的纸尿裤内部微气候模型具体为:
如图2所示,所述纸尿裤内部微气候模型包括从内到外的三部分,具体为皮肤、纸尿裤和外界环境;
所述纸尿裤部分包括从内到外的面层、覆盖层、吸水芯体层和底膜层;
所述纸尿裤各层的厚度量级为毫米级;
所述皮肤和外界环境均为无厚度边界,皮肤边界为恒定温度;
基于以下假设:
假设一:纸尿裤内起主要锁水作用的吸水芯体层保水性强,吸水芯体层内部水分在短时间内传递缓慢,因此忽略纸尿裤内部的对流换热,将传热过程视为导热过程;
假设二:尿液在纸尿裤内扩散的区域有限,尿液所在区域的换热量远大于其他区域,因此将计算区域设置为尿液存在的区域;
假设三:纸尿裤内部各层之间存在间隙,在尿液存在的区域内含水量很高,因此本模型假设纸尿裤各层间隙内的空气湿度一直处于该位置当前温度下的饱和湿度,同时忽略空气的对流传热和纸尿裤侧边界的传热;
假设四:尿液存在的区域所有边界均存在蒸发相变过程和热量传递过程,对于同种材料来说侧边界水分的流失速度大于中心区域,因此本模型考虑水分质量流失速度在皮肤切向上的差异。
如图3所示,所述纸尿裤各层含水量的切向分布呈二次函数分布,在法向方向上具有线性关系;
所述步骤S2中,对建立的纸尿裤内部微气候模型进行初始化的具体方法为:
所述初始化包括设置含水量关系拟合函数、相变类型、初始条件和边界条件;
所述含水量关系拟合函数包括对面层、覆盖层、吸水芯体层和底膜层均设置导热系数-含水量关系拟合函数、密度-含水量拟合函数、比热容-含水量拟合函数、水分传输比例-时间关系拟合函数和水分传输比例-水平分布拟合函数;
设置纸尿裤内部微气候模型的初始条件具体为:
设置纸尿裤内部微气候模型各部分的初始温度、纸尿裤的相对湿度、外界环境的初始温度和相对湿度,以及排尿次数;
设置纸尿裤内部微气候模型的边界条件具体为:
设置皮肤表面的边界条件为预设温度的恒温;
设置纸尿裤的边界条件为纸尿裤最外侧与外界环境的热交换过程包括蒸发和对流;
设置对流换热系数和水蒸气的蒸发速率;
所述步骤S3中,对初始化后的纸尿裤内部微气候模型中的纸尿裤部分进行网格划分,利用有限体积法计算在一定排尿量的条件下,不同时刻纸尿裤每个网格的温度和含水量,具体方法为:
如图4所示,所述纸尿裤的能量变化包括皮肤输入的热量、失去水分的内能、纸尿裤最外层对流换热和液态水汽化潜热;
纸尿裤每个网格的能量控制方程为
其中,为本时间步密度,/>为上一时间步密度,/>为本时间步比热容,/>为上一时间步比热容,/>为上一时间步导热系数,/>为网格宽度,/>为网格高度,/>为网格中心与右侧相邻网格中心距离,/>为网格中心与左侧相邻网格中心距离,/>为网格中心与上侧相邻网格中心距离,/>为网格中心与下侧相邻网格中心距离,/>为本时间步网格中心温度,/>为上一时间步网格中心温度,/>为右侧相邻网格中心点温度,/>为左侧相邻网格中心点温度,/>为上侧相邻网格中心点温度,/>为下侧相邻网格中心点温度;
、/>、/>、/>、/>和/>的系数分别为:
其中,为上一时间步温度系数,/>为本时间步温度系数,/>为右侧网格温度系数,/>为左侧网格温度系数,/>为上侧网格温度系数,/>为下侧网格温度系数;
根据每个网格的能量控制方程计算不同时刻纸尿裤各层的能量变化以及每个网格的温度;
纸尿裤每个网格的含水量计算方程为:
其中,为纸尿裤第i层的水分传递比例,/>为水平方向上的分配比例,为时间当前时间下的水分蒸发量,/>为新时间步网格含水量,为上一时间步网格含水量;
纸尿裤每个网格的总含水量为各层含水量之和;
根据每个网格的含水量计算方程计算不同时刻纸尿裤各层的含水量,根据每个网格不同时刻纸尿裤各层的含水量计算每个网格不同时刻纸尿裤的总含水量;
所述皮肤输入的热量具体为:
其中,为皮肤温度,/>为皮肤相邻网格温度,/>为导热面积,具体为单个网格的面积,/>为单个网格的导热系数;
所述纸尿裤最外层对流换热具体为:
其中,h为对流换热系数,表示对流强度,为外界环境温度,/>为时间步长;
所述液态水汽化潜热具体为:
其中,为水蒸气质量。
在具体实施过程中,首先获取纸尿裤的热力学参数,在本实施例中,首先通过实验法获取纸尿裤各层材料的导热系数、密度和孔隙率;
基于以下假设,建立纸尿裤内部微气候模型:
假设一:纸尿裤内起主要锁水作用的吸水芯体层保水性强,吸水芯体层内部水分在短时间内传递缓慢,因此忽略纸尿裤内部的对流换热,将传热过程视为导热过程;
假设二:尿液在纸尿裤内扩散的区域有限,尿液所在区域的换热量远大于其他区域,因此将计算区域设置为尿液存在的区域;
假设三:纸尿裤内部各层之间存在间隙,在尿液存在的区域内含水量很高,因此本模型假设纸尿裤各层间隙内的空气湿度一直处于该位置当前温度下的饱和湿度,同时忽略空气的对流传热和纸尿裤侧边界的传热;
假设四:尿液存在的区域所有边界均存在蒸发相变过程和热量传递过程,对于同种材料来说侧边界水分的流失速度大于中心区域,因此本模型考虑水分质量流失速度在皮肤切向上的差异;
所建立的纸尿裤内部微气候模型包括从内到外的三部分,具体为皮肤、纸尿裤和外界环境;
所述纸尿裤部分包括从内到外的面层、覆盖层、吸水芯体层和底膜层;
所述纸尿裤各层的厚度量级均为毫米级,在本实施例中,面层厚度为0.4mm,覆盖层厚度为0.2mm,吸水芯体层厚度为3mm,底膜层厚度为0.1mm;
并对纸尿裤部分进行网格划分,将纸尿裤部分划分为n×m的网格;
对建立的纸尿裤内部微气候模型进行初始化,具体为:设置含水量关系拟合函数、相变类型、初始条件和边界条件;
所述含水量关系拟合函数包括对面层、覆盖层、吸水芯体层和底膜层均设置导热系数-含水量关系拟合函数、密度-含水量拟合函数、比热容-含水量拟合函数、水分传输比例-时间关系拟合函数和水分传输比例-水平分布拟合函数;
设置纸尿裤内部微气候模型的初始条件具体为:
设置纸尿裤内部微气候模型各部分的初始温度、纸尿裤的相对湿度、外界环境的初始温度和相对湿度;
在本实施例中,纸尿裤内部微气候模型各部分的初始温度为36℃,纸尿裤的相对湿度为50%,外界环境的初始温度为25℃、相对湿度为50%,排尿次数设置为3次,每次排尿量为80mL;
设置纸尿裤内部微气候模型的边界条件具体为:
设置皮肤表面的边界条件为36℃的恒温;
设置纸尿裤的边界条件为纸尿裤最外侧与外界环境的热交换过程包括蒸发和对流;
设置对流换热系数和水蒸气的蒸发速率,在本实施例中,对流换热系数为8 W/(m2·K);
利用有限体积法计算在一定排尿量的条件下,不同时刻纸尿裤每个网格的温度和含水量;
所述纸尿裤的能量变化包括皮肤输入的热量、失去水分的内能、纸尿裤最外层对流换热和液态水汽化潜热;
纸尿裤每个网格的能量控制方程为
其中,为本时间步密度,/>为上一时间步密度,/>为本时间步比热容,/>为上一时间步比热容,/>为上一时间步导热系数,/>为网格宽度,/>为网格高度,/>为网格中心与右侧相邻网格中心距离,/>为网格中心与左侧相邻网格中心距离,/>为网格中心与上侧相邻网格中心距离,/>为网格中心与下侧相邻网格中心距离,/>为本时间步网格中心温度,/>为上一时间步网格中心温度,/>为右侧相邻网格中心点温度,/>为左侧相邻网格中心点温度,/>为上侧相邻网格中心点温度,/>为下侧相邻网格中心点温度;
、/>、/>、/>、/>和/>的系数分别为:
其中,为上一时间步温度系数,/>为本时间步温度系数,/>为右侧网格温度系数,/>为左侧网格温度系数,/>为上侧网格温度系数,/>为下侧网格温度系数;
根据每个网格的能量控制方程计算不同时刻纸尿裤各层的能量变化以及每个网格的温度;
纸尿裤每个网格的含水量计算方程为:
其中,为纸尿裤第i层的水分传递比例,/>为水平方向上的分配比例,为时间当前时间下的水分蒸发量,/>为新时间步网格含水量,为上一时间步网格含水量;
纸尿裤每个网格的总含水量为各层含水量之和;
根据每个网格的含水量计算方程计算不同时刻纸尿裤各层的含水量,根据每个网格不同时刻纸尿裤各层的含水量计算每个网格不同时刻纸尿裤的总含水量;
根据不同时刻纸尿裤每个网格的温度和含水量获得皮肤损失的热量、不同时刻纸尿裤的温度分布和湿度分布,完成纸尿裤内部热湿环境的计算;
所述皮肤输入的热量具体为:
其中,为皮肤温度,/>为皮肤相邻网格温度,/>为导热面积,具体为单个网格的面积,/>为单个网格的导热系数;
所述纸尿裤最外层对流换热具体为:
其中,h为对流换热系数,表示对流强度,为外界环境温度,/>为时间步长;
所述液态水汽化潜热具体为:
其中,为水蒸气质量;
本实施例中纸尿裤内部热湿环境的计算结果如图5所示,图5精确地展示出在14s和86s两个时刻纸尿裤内部的温度分布和湿度分布,计算精度较高,效果较好;
利用本实施例中的方法可以根据尿液所在区域大小不同计算婴儿不同尿液量时纸尿裤的热量损失;还可以预测不同材料搭配和厚度比例时纸尿裤的热损失;也可以通过设置户外室内不同的对流换热系数,或冬天夏天不同的环境温度来模拟在不同的外界环境下,纸尿裤内部的热湿环境,以指导开发出适应各种要求的健康舒适产品;
本发明基于有限体积法进行计算,运算精度高、计算周期短,能快速精准地反映纸尿裤内微气候复杂的温湿度场及流场的分布情况,能够系统高效地预测纸尿裤微气候环境的湿热舒适性,揭示不同结构材料对湿热传递性能的影响,对个性化产品快速优化设计和保护婴幼儿皮肤健康,具有实际应用价值和重要现实意义。
实施例3
如图6所示,本实施例提供一种纸尿裤内部微气候热湿环境计算系统,应用实施例1或2所述的一种纸尿裤内部微气候热湿环境计算方法,包括:
模型构建和网格划分单元301:用于获取纸尿裤的热力学参数,建立纸尿裤内部微气候模型;
所述纸尿裤内部微气候模型包括从内到外的三部分,具体为皮肤、纸尿裤和外界环境;
并对纸尿裤内部微气候模型中的纸尿裤部分进行网格划分;
初始化单元302:用于对建立的纸尿裤内部微气候模型进行初始化;
计算单元303:用于利用有限体积法计算在一定排尿量的条件下,不同时刻纸尿裤每个网格的温度和含水量;
计算结果处理单元304:用于根据不同时刻纸尿裤每个网格的温度和含水量获得皮肤损失的热量、不同时刻纸尿裤的温度分布和湿度分布,完成纸尿裤内部热湿环境的计算。
在具体实施过程中,首先模型构建和网格划分单元301获取纸尿裤的热力学参数,建立纸尿裤内部微气候模型,并对纸尿裤内部微气候模型中的纸尿裤部分进行网格划分;之后初始化单元302对建立的纸尿裤内部微气候模型进行初始化;计算单元303利用有限体积法计算在一定排尿量的条件下,不同时刻纸尿裤每个网格的温度和含水量;最后计算结果处理单元304根据不同时刻纸尿裤每个网格的温度和含水量获得皮肤损失的热量、不同时刻纸尿裤的温度分布和湿度分布,完成纸尿裤内部热湿环境的计算;
本发明基于有限体积法进行计算,运算精度高、计算周期短,能快速精准地反映纸尿裤内微气候复杂的温湿度场及流场的分布情况,能够系统高效地预测纸尿裤微气候环境的湿热舒适性,揭示不同结构材料对湿热传递性能的影响,对个性化产品快速优化设计和保护婴幼儿皮肤健康,具有实际应用价值和重要现实意义。
相同或相似的标号对应相同或相似的部件;
附图中描述位置关系的用语仅用于示例性说明,不能理解为对本专利的限制;
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明权利要求的保护范围之内。
Claims (7)
1.一种纸尿裤内部微气候热湿环境计算方法,其特征在于,包括以下步骤:
S1:获取纸尿裤的热力学参数,建立纸尿裤内部微气候模型,具体为:
所述纸尿裤内部微气候模型包括从内到外的三部分,具体为皮肤、纸尿裤和外界环境;
所述纸尿裤部分包括从内到外的面层、覆盖层、吸水芯体层和底膜层;
所述纸尿裤各层的厚度量级为毫米级;
所述皮肤和外界环境均为无厚度边界,皮肤边界为恒定温度;
基于以下假设:
假设一:纸尿裤内起主要锁水作用的吸水芯体层保水性强,吸水芯体层内部水分在短时间内传递缓慢,因此忽略纸尿裤内部的对流换热,将传热过程视为导热过程;
假设二:尿液在纸尿裤内扩散的区域有限,尿液所在区域的换热量远大于其他区域,因此将计算区域设置为尿液存在的区域;
假设三:纸尿裤内部各层之间存在间隙,在尿液存在的区域内含水量很高,因此本模型假设纸尿裤各层间隙内的空气湿度一直处于尿液存在的区域当前温度下的饱和湿度,同时忽略空气的对流传热和纸尿裤侧边界的传热;
假设四:尿液存在的区域所有边界均存在蒸发相变过程和热量传递过程,对于同种材料来说侧边界水分的流失速度大于中心区域,因此本模型考虑水分质量流失速度在皮肤切向上的差异;
对纸尿裤内部微气候模型中的纸尿裤部分进行网格划分;
S2:对建立的纸尿裤内部微气候模型进行初始化;
所述初始化包括设置含水量关系拟合函数、相变类型、初始条件和边界条件;
所述含水量关系拟合函数包括对面层、覆盖层、吸水芯体层和底膜层均设置导热系数-含水量关系拟合函数、密度-含水量拟合函数、比热容-含水量拟合函数、水分传输比例-时间关系拟合函数和水分传输比例-水平分布拟合函数;
设置纸尿裤内部微气候模型的初始条件具体为:
设置纸尿裤内部微气候模型各部分的初始温度、纸尿裤的相对湿度、外界环境的初始温度和相对湿度,以及排尿次数;
设置纸尿裤内部微气候模型的边界条件具体为:
设置皮肤表面的边界条件为预设温度的恒温;
设置纸尿裤的边界条件为纸尿裤最外侧与外界环境的热交换过程包括蒸发和对流;
设置对流换热系数和水蒸气的蒸发速率;
S3:利用有限体积法计算在一定排尿量的条件下,不同时刻纸尿裤每个网格的温度和含水量,具体方法为:
所述纸尿裤的能量变化包括皮肤输入的热量、失去水分的内能、纸尿裤最外层对流换热和液态水汽化潜热;
纸尿裤每个网格的能量控制方程为
其中,为本时间步密度,/>为上一时间步密度,/>为本时间步比热容,/>为上一时间步比热容,/>为上一时间步导热系数,/>为网格宽度,/>为网格高度,/>为网格中心与右侧相邻网格中心距离,/>为网格中心与左侧相邻网格中心距离,/>为网格中心与上侧相邻网格中心距离,/>为网格中心与下侧相邻网格中心距离,/>为本时间步网格中心温度,/>为上一时间步网格中心温度,/>为右侧相邻网格中心点温度,/>为左侧相邻网格中心点温度,/>为上侧相邻网格中心点温度,/>为下侧相邻网格中心点温度;
、/>、/>、/>、/>和/>的系数分别为:
其中,为上一时间步温度系数,/>为本时间步温度系数,/>为右侧网格温度系数,为左侧网格温度系数,/>为上侧网格温度系数,/>为下侧网格温度系数;
根据每个网格的能量控制方程计算不同时刻纸尿裤各层的能量变化以及每个网格的温度;
纸尿裤每个网格的含水量计算方程为:
其中,为纸尿裤第i层的水分传递比例,/>为水平方向上的分配比例,为时间当前时间下的水分蒸发量,/>为新时间步网格含水量,为上一时间步网格含水量;
纸尿裤每个网格的总含水量为各层含水量之和;
根据每个网格的含水量计算方程计算不同时刻纸尿裤各层的含水量,根据每个网格不同时刻纸尿裤各层的含水量计算每个网格不同时刻纸尿裤的总含水量;
S4:根据不同时刻纸尿裤每个网格的温度和含水量获得皮肤损失的热量、不同时刻纸尿裤的温度分布和湿度分布,完成纸尿裤内部热湿环境的计算。
2.根据权利要求1所述的一种纸尿裤内部微气候热湿环境计算方法,其特征在于,所述步骤S1中所述的纸尿裤的热力学参数包括:纸尿裤各层材料的导热系数、密度和孔隙率。
3.根据权利要求2所述的一种纸尿裤内部微气候热湿环境计算方法,其特征在于,所述纸尿裤各层含水量的切向分布呈二次函数分布,在法向方向上具有线性关系。
4.根据权利要求1所述的一种纸尿裤内部微气候热湿环境计算方法,其特征在于,所述皮肤输入的热量具体为:
其中,为皮肤温度,/>为皮肤相邻网格温度,/>为导热面积,具体为单个网格的面积,/>为单个网格的导热系数。
5.根据权利要求1所述的一种纸尿裤内部微气候热湿环境计算方法,其特征在于,所述纸尿裤最外层对流换热具体为:
其中,h为对流换热系数,表示对流强度,为外界环境温度,/>为时间步长。
6.根据权利要求1所述的一种纸尿裤内部微气候热湿环境计算方法,其特征在于,所述液态水汽化潜热具体为:
其中,为水蒸气质量。
7.一种纸尿裤内部微气候热湿环境计算系统,应用权利要求1~6任意一项中所述的一种纸尿裤内部微气候热湿环境计算方法,其特征在于,包括:
模型构建和网格划分单元:用于获取纸尿裤的热力学参数,建立纸尿裤内部微气候模型;
所述纸尿裤内部微气候模型包括从内到外的三部分,具体为皮肤、纸尿裤和外界环境;
并对纸尿裤内部微气候模型中的纸尿裤部分进行网格划分;
初始化单元:用于对建立的纸尿裤内部微气候模型进行初始化;
计算单元:用于利用有限体积法计算在一定排尿量的条件下,不同时刻纸尿裤每个网格的温度和含水量;
计算结果处理单元:用于根据不同时刻纸尿裤每个网格的温度和含水量获得皮肤损失的热量、不同时刻纸尿裤的温度分布和湿度分布,完成纸尿裤内部热湿环境的计算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310393847.XA CN116110525B (zh) | 2023-04-13 | 2023-04-13 | 一种纸尿裤内部微气候热湿环境计算方法和系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310393847.XA CN116110525B (zh) | 2023-04-13 | 2023-04-13 | 一种纸尿裤内部微气候热湿环境计算方法和系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116110525A CN116110525A (zh) | 2023-05-12 |
CN116110525B true CN116110525B (zh) | 2023-09-08 |
Family
ID=86260148
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310393847.XA Active CN116110525B (zh) | 2023-04-13 | 2023-04-13 | 一种纸尿裤内部微气候热湿环境计算方法和系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116110525B (zh) |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105928566A (zh) * | 2016-07-13 | 2016-09-07 | 湖南康程护理用品有限公司 | 一种纸尿裤温湿度检测系统及方法 |
CN106407532A (zh) * | 2016-10-18 | 2017-02-15 | 中山大学 | 一种着装人体热质传递并行仿真方法及其系统 |
CN108647400A (zh) * | 2018-04-18 | 2018-10-12 | 华南理工大学 | 一种多尺度框架下的热物理传输数值仿真方法 |
CN110907613A (zh) * | 2019-12-31 | 2020-03-24 | 露乐健康科技(广州)有限公司 | 婴幼儿纸尿裤“微气候”的测量评价装置及评价方法 |
CN111444594A (zh) * | 2020-03-04 | 2020-07-24 | 湖南科技大学 | 一种用于优化气体降温服结构的模拟方法 |
CN114324768A (zh) * | 2021-12-31 | 2022-04-12 | 露乐健康科技股份有限公司 | 纸尿裤微气候自动化动态监测系统及方法 |
CN216791978U (zh) * | 2021-12-31 | 2022-06-21 | 露乐健康科技股份有限公司 | 用于纸尿裤内部环境监测的仿真人偶 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7099734B2 (en) * | 2003-05-22 | 2006-08-29 | Kimberly-Clark Worldwide, Inc. | Method of evaluating the performance of a product using a virtual environment |
US7373284B2 (en) * | 2004-05-11 | 2008-05-13 | Kimberly-Clark Worldwide, Inc. | Method of evaluating the performance of a product using a virtual environment |
-
2023
- 2023-04-13 CN CN202310393847.XA patent/CN116110525B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105928566A (zh) * | 2016-07-13 | 2016-09-07 | 湖南康程护理用品有限公司 | 一种纸尿裤温湿度检测系统及方法 |
CN106407532A (zh) * | 2016-10-18 | 2017-02-15 | 中山大学 | 一种着装人体热质传递并行仿真方法及其系统 |
CN108647400A (zh) * | 2018-04-18 | 2018-10-12 | 华南理工大学 | 一种多尺度框架下的热物理传输数值仿真方法 |
CN110907613A (zh) * | 2019-12-31 | 2020-03-24 | 露乐健康科技(广州)有限公司 | 婴幼儿纸尿裤“微气候”的测量评价装置及评价方法 |
CN111444594A (zh) * | 2020-03-04 | 2020-07-24 | 湖南科技大学 | 一种用于优化气体降温服结构的模拟方法 |
CN114324768A (zh) * | 2021-12-31 | 2022-04-12 | 露乐健康科技股份有限公司 | 纸尿裤微气候自动化动态监测系统及方法 |
CN216791978U (zh) * | 2021-12-31 | 2022-06-21 | 露乐健康科技股份有限公司 | 用于纸尿裤内部环境监测的仿真人偶 |
Non-Patent Citations (1)
Title |
---|
Topology optimization of unsaturated flows in multi-material porous media:application to a simple diaper model;Fabian Wein等;《Communications in Nonlinear Science and Numerical Simulation》;第第78卷卷;第1-18页 * |
Also Published As
Publication number | Publication date |
---|---|
CN116110525A (zh) | 2023-05-12 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Tanabe et al. | Evaluation of thermal comfort using combined multi-node thermoregulation (65MN) and radiation models and computational fluid dynamics (CFD) | |
Atthajariyakul et al. | Real-time determination of optimal indoor-air condition for thermal comfort, air quality and efficient energy usage | |
Błażejczyk et al. | Principles of the new Universal Thermal Climate Index (UTCI) and its application to bioclimatic research in European scale | |
Catalina et al. | Evaluation of thermal comfort using combined CFD and experimentation study in a test room equipped with a cooling ceiling | |
Tokunaga et al. | Human-body exergy balance calculation under un-steady state conditions | |
Zhai et al. | On approaches to couple energy simulation and computational fluid dynamics programs | |
Prek | Thermodynamical analysis of human thermal comfort | |
Nagano et al. | Effects of ambient temperature steps on thermal comfort requirements | |
Zolfaghari et al. | A new simplified model for evaluating non-uniform thermal sensation caused by wearing clothing | |
Mao et al. | Comparative studies on using RSM and TOPSIS methods to optimize residential air conditioning systems | |
Norton et al. | Improving the representation of thermal boundary conditions of livestock during CFD modelling of the indoor environment | |
Stavrakakis et al. | Development of a computational tool to quantify architectural-design effects on thermal comfort in naturally ventilated rural houses | |
Wan et al. | A transient thermal model of the human body–clothing–environment system | |
Mao et al. | PMV-based dynamic optimization of energy consumption for a residential task/ambient air conditioning system in different climate zones | |
Vasilyev et al. | Simulation of heat and moisture transfer in a multiplex structure | |
Zhang et al. | Exploring buildings’ secrets: The ideal thermophysical properties of a building’s wall for energy conservation | |
ketan Nayak et al. | A simplified numerical model for evaporative cooling by water spray over roof surfaces | |
Lebon et al. | Numerical analysis and field measurements of the airflow patterns and thermal comfort in an indoor swimming pool: a case study | |
He et al. | Measurement and evaluation of the summer microclimate in the semi-enclosed space under a membrane structure | |
Nagano et al. | Development of outdoor thermal index indicating universal and separate effects on human thermal comfort | |
Caliskan | Energetic and exergetic comparison of the human body for the summer season | |
Zolfaghari et al. | A new predictive index for evaluating both thermal sensation and thermal response of the human body | |
Schweiker et al. | Unsteady-state human-body exergy consumption rate and its relation to subjective assessment of dynamic thermal environments | |
Ghaddar et al. | Integrated human-clothing system model for estimating the effect of walking on clothing insulation | |
Xiao et al. | Proposal and application of comprehensive thermal comfort evaluation model in heating seasons for buildings with solar Trombe wall |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB02 | Change of applicant information | ||
CB02 | Change of applicant information |
Address after: 11 Yucai Road, Hung Hom, Hong Kong, China Applicant after: The Hong Kong Polytechnic University Applicant after: Lule Health Technology Co.,Ltd. Address before: Room 188, room 406, 1 Yichuang street, Zhongxin Guangzhou Knowledge City, Huangpu District, Guangzhou City, Guangdong Province, 510555 Applicant before: Lule Health Technology Co.,Ltd. Applicant before: The Hong Kong Polytechnic University |
|
GR01 | Patent grant | ||
GR01 | Patent grant |