CN108682046A - 带有柔性基材的悬浮喷嘴计算方法 - Google Patents

带有柔性基材的悬浮喷嘴计算方法 Download PDF

Info

Publication number
CN108682046A
CN108682046A CN201810471784.4A CN201810471784A CN108682046A CN 108682046 A CN108682046 A CN 108682046A CN 201810471784 A CN201810471784 A CN 201810471784A CN 108682046 A CN108682046 A CN 108682046A
Authority
CN
China
Prior art keywords
parent metal
flexible parent
suspension nozzle
flow field
suspension
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
CN201810471784.4A
Other languages
English (en)
Other versions
CN108682046B (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.)
Shenzhen Xinyuren Polytron Technologies Inc
Shenzhen Xinyuren Technology Co Ltd
Original Assignee
Shenzhen Xinyuren Polytron Technologies Inc
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 Shenzhen Xinyuren Polytron Technologies Inc filed Critical Shenzhen Xinyuren Polytron Technologies Inc
Priority to CN201810471784.4A priority Critical patent/CN108682046B/zh
Publication of CN108682046A publication Critical patent/CN108682046A/zh
Application granted granted Critical
Publication of CN108682046B publication Critical patent/CN108682046B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/10—Constructive solid geometry [CSG] using solid primitives, e.g. cylinders, cubes
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
    • G06T17/205—Re-meshing

Landscapes

  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Application Of Or Painting With Fluid Materials (AREA)
  • Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)

Abstract

本发明属于R2R过程中悬浮干燥相关领域，并公开了一种带有柔性基材的悬浮喷嘴设计计算方法，该方法包括以下步骤：(1)建立悬浮喷嘴系统流体域和柔性基材的几何模型；(2)对计算域进行网格划分；(3)给定边界条件和初始条件；(4)利用数值模型求解柔性基材的变形与换热；(5)根据计算结果确定悬浮喷嘴的设计参数对变形和换热的影响，知道悬浮喷嘴的设计。本发明以流固耦合模拟方法作为基础，通过计算来指导悬浮喷嘴的设计，可有效降低悬浮喷嘴的制造成本并提高悬浮喷嘴的换热能力。本方法充分考虑了柔性基材变形对基材表面压力、换热特性的影响，通过对FLUENT进行二次开发，嵌入基材变形程序，提高了数值预测精度与可靠性。

Description

带有柔性基材的悬浮喷嘴计算方法
技术领域
本发明涉及计算机数值模拟及CAE领域，特别涉及带有柔性基材的悬浮喷嘴设计计算方法。
背景技术
制造业越来越多地采用R2R方法，悬浮干燥工艺是影响R2R制造速度的决定性工艺，基材经双面涂布之后，需要在悬浮烘箱中干燥。悬浮烘箱需要将基材稳定的悬浮于上下喷嘴之间，悬浮喷嘴的悬浮效果及换热能力直接影响干燥效率和干燥品质。现有的悬浮喷嘴计算绝大多数将基材视为刚体，不考虑基材的变形，这种计算方法无法准确考虑基材的变形以及变形对悬浮喷嘴换热能力的影响，造成设计的悬浮喷嘴难以将基材稳定的悬浮于上下喷嘴之间，造成基材与悬浮烘箱内的零部件发生刮擦，导致设计成本较高。
目前悬浮喷嘴的主流计算方法是将基材视为刚性，即不考虑基材的变形。例如文献(刘金金.悬浮式烘箱数值模拟分析及结构优化[D].西安理工大学,2015.)采用CFD方法计算了带有刚性基材的悬浮喷嘴的流场与热场分布。又如文献(卫延斌,侯和平,张海燕.悬浮烘干系统烘箱体结构优化研究[J].西安理工大学学报,2015(2):171-175.张海燕,刘金金,薛志成,等.涂布机悬浮烘箱的结构参数优化[J].机械设计与研究,2015(4):142-146.李徐佳,高殿荣,邸立明,等.基于流场数值模拟的锂电池极片干燥箱结构改进[J].中国机械工程,2010,21(18):2183-2187.李徐佳.锂电池极片干燥箱流场特性研究[D].燕山大学,2012.)等都介绍了悬浮喷嘴及悬浮烘箱的计算方法。然而，如锂离子电池极片、纸张等柔性基材在悬浮干燥烘箱干燥过程中，其会发生明显的变形，因此如果不考虑基材的变形，由计算模型所得到的计算结果可能偏离实际，导致计算结果不能准确地指导悬浮喷嘴的结构设计。
发明内容
本发明主要针对目前悬浮喷嘴模拟计算方法中存在的不足，提供一种带有柔性基材的悬浮喷嘴设计计算方法，该方法能够准确地分析出悬浮烘箱中柔性基材的变形以及变形对悬浮喷嘴换热能力的影响。基于此，该方法包括以下步骤：
步骤一、建立悬浮喷嘴及柔性基材几何模型：悬浮喷嘴在悬浮烘箱中呈周期性分布，且上下悬浮喷嘴之间的间距以及悬浮喷嘴的宽度远大于悬浮喷嘴出风口尺寸，对几何模型做如下简化：柔性基材两端简支梁支撑，仅考虑柔性基材在空气流场作用下垂直方向的变形。
步骤二、流场网格划分：对步骤一建立的流场几何模型进行网格划分，在柔性基材壁面增加网格密度，以便更好地捕捉柔性基材表面的换热特性。
步骤三、建立计算流体力学模型和柔性基材变形模型。计算流体力学模型包括连续性方程、动量方程和能量方程，分别为：
式中u和v分别表示流体流动速度张量与动网格速度张量，ρ,μ以及Pr表示空气的密度、动力粘度以及普朗特数，P为流场压力，T为空气温度，t为时间。采用SST k-w湍流模型封闭方程中的雷诺应力与雷诺热通量。
柔性基材控制方程为：
柔性基材边界条件和初始条件为：
y|x＝-L/2＝y|x＝L/2＝0 (5)
y(t,0)＝y(t,L)＝0 (7)
m为柔性基材单位面积的质量，y和Vw表示柔性基材的垂直位移与柔性基材的轴向运动速度，F表示柔性基材中的张力。
步骤四、进行初始定常流场数值计算：在FLUENT中，对计算参数进行初始化。流场入口给定速度和高恒温边界条件，出口给定平均静压以及出口低恒温边界条件，柔性基材表面给定无滑移、光滑低恒温壁面边界条件。基于上述边界条件和初始条件，不考虑流场参数(速度、压力以及温度)随时间的变化，不考虑柔性基材的变形，利用FLUENT中的SIMPLEC求解器进行订场流场数值计算，可以得到刚性壁面的流场区域的速度、压力以及温度分布。
步骤五、进行流固耦合数值计算：将步骤四中的定常流场数值计算结果作为初始条件进行流固耦合的数值计算。流固耦合数值计算方法为：
步骤5.1、读取柔性基材壁面的压力，通过方程(4)计算出柔性基材的变形。
步骤5.2、根据步骤5.1所计算的变形结果，通过FLUENT中的扩散光顺方法对网格进行更新，得到新的柔性基材表面位置。
步骤5.3、根据新的柔性基材表面位置，计算流体域的速度、压力以及温度分布。
步骤5.4、重复步骤5.1至步骤5.3的计算过程，直至达到预定求解时间或者数值计算结果收敛，得到最终柔性基材的变形量以及最终的流场与热场分布。
附图说明
图1悬浮喷嘴与柔性基材相对位置示意图。
图2是图1的截面结构示意图。
图3带有柔性基材的悬浮喷嘴几何模型示意图。
图4是带有柔性基材的悬浮喷嘴的流场控制域示意图。
图5是带有柔性基材的悬浮喷嘴的计算模型边框条件示意图。
图6是带有柔性基材的悬浮喷嘴的计算网格示意图。
图7是带有柔性基材的悬浮喷嘴的流固耦合计算流程示意图。
图8不考虑柔性基材变形的速度场与考虑柔性基材变形的速度场分布示意图。
图9不考虑柔性基材变形的温度场与考虑柔性基材变形的温度场分布示意图。
具体实施方式
实施例：
下面结合附图，对本发明的具体实施方式进行详细说明。
步骤一：建立悬浮喷嘴及柔性基材几何模型，见图1、图2和图3。
针对给定的悬浮喷嘴，在CATIA中建立流场与柔性基材的几何模型，如图4所示。柔性基材材料为铜箔，主要材料属性为：密度ρ＝8900kg/m3，弹性模量E＝121GPa，厚度t＝10μm。对柔性基材的变形做如下简化：柔性基材两端简支梁支撑，仅考虑柔性基材在垂直方向的变形。流场介质为空气：密度ρ＝1.225kg/m3，动力粘度μ＝1.7895e-5kg/(m·s)。
步骤二：流场网格划分。
对流场区域进行网格划分，如图5和图6所示，整个流场采用结构化网格，上下靠近壁面位置的网格进行加密，悬浮喷嘴出风口位置网格也进行加密，将建立的网格文件保存输出为*.msh格式。
步骤三：建立计算流体力学模型和柔性基材变形模型(见图7)。
计算流体力学模型包括连续性方程、动量方程和能量方程，分别为：
式中u和v分别表示流体流动速度张量与动网格速度张量，ρ,μ以及Pr表示空气的密度、动力粘度以及普朗特数，P为流场压力，T为空气温度，t为时间。采用SST k-w湍流模型封闭方程中的雷诺应力与雷诺热通量。
柔性基材控制方程为：
柔性基材边界条件和初始条件为：
y|x＝-L/2＝y|x＝L/2＝0 (5)
y(t,0)＝y(t,L)＝0 (7)
m为柔性基材单位面积的质量，y和Vw表示柔性基材的垂直位移与柔性基材的轴向运动速度，F表示柔性基材中的张力。
步骤四、进行初始定常流场数值计算：在FLUENT中，对计算参数进行初始化。流场入口给定速度和高恒温边界条件，出口给定平均静压以及出口低恒温边界条件，柔性基材表面给定无滑移、光滑低恒温壁面边界条件。基于上述边界条件和初始条件，不考虑流场参数(速度、压力以及温度)随时间的变化，不考虑柔性基材的变形，利用FLUENT中的SIMPLEC求解器进行订场流场数值计算，其中压力采用PRESTO离散格式，扩散项采用中心差分离散格式，对流上采用二阶迎风格式，可以得到刚性壁面的流场区域的速度、压力以及温度分布。
步骤五、进行流固耦合数值计算：将步骤四中的定常流场数值计算结果作为初始条件进行流固耦合的数值计算。流固耦合数值计算方法为：
步骤5.1、读取柔性基材壁面的压力，通过方程(4)计算出柔性基材的变形。采用如下的方法计算柔性基材的变形，计算流程如图3所示：
为求解方程(4)，基于柔性基材的变形特点，采用余弦级数近似柔性基材的变形：
变形方程能够自动满足边界条件和初始条件，代入到控制方程(4)，可得控制方程的残差：
并采用Galerkin最小二乘法求解余弦级数中的系数，其系数的求解方程为：
采用中心差分法可求出每个时间步的系数：
步骤5.2、根据步骤5.1所计算的变形结果，利用FLUENT自带User-Defined-Function功能，通过FLUENT中的扩散光顺方法对网格进行更新，得到新的柔性基材表面位置。
步骤5.3、根据新的柔性基材表面位置，计算流体域的速度、压力以及温度分布。
步骤5.4、重复步骤5.1至步骤5.3的计算过程，直至达到预定求解时间或者数值计算结果收敛，得到最终柔性基材的变形量以及最终的流场与热场分布。
图8不考虑柔性基材变形的速度场与考虑柔性基材变形的速度场分布示意图。
图8中，图(a)-(c)为雷诺数Re＝4000～6000的带有刚性基材的悬浮喷嘴速度场分布；
图(d)-(f)为雷诺数Re＝4000～6000的带有柔性基材的悬浮喷嘴速度场分布；
由图8可知，带有刚性基材的悬浮喷嘴速度场速度最大点位于基材位置附近，带有柔性基材的悬浮喷嘴速度场速度最大点位于入口位置附近。
图9不考虑柔性基材变形的温度场与考虑柔性基材变形的温度场分布示意图。
图9中，图(a)-(c)为雷诺数Re＝4000～6000的带有刚性基材的悬浮喷嘴温度场分布；
(d)-(f)为雷诺数Re＝4000～6000的带有柔性基材的悬浮喷嘴温度场分布；
由图9可知，带有刚性基材的悬浮喷嘴与带有柔性基材的悬浮喷嘴的温度场在入口与对称面位置之间明显不同。
因此，如果不考虑柔性基材的变形，悬浮喷嘴设计的结果将偏离实际情形，造成返工成本高，设计效率低下。

Claims (3)

1.一种带有柔性基材的悬浮喷嘴计算方法，其特征在于包括以下步骤：
步骤一、建立悬浮喷嘴及柔性基材几何模型：悬浮喷嘴在悬浮烘箱中呈周期性分布，且上下悬浮喷嘴之间的间距以及悬浮喷嘴的宽度远大于悬浮喷嘴出风口尺寸，对几何模型做如下简化：柔性基材两端简支梁支撑，仅考虑柔性基材在空气流场作用下垂直方向的变形；
步骤二、流场网格划分：对步骤一建立的流场几何模型进行网格划分，在柔性基材壁面增加网格密度，以便更好地捕捉柔性基材表面的换热特性；
步骤三、建立计算流体力学模型和柔性基材变形模型；计算流体力学模型包括连续性方程、动量方程和能量方程，分别为：
式中u和v分别表示流体流动速度张量与动网格速度张量，ρ,μ以及Pr表示空气的密度、动力粘度以及普朗特数，P为流场压力，T为空气温度，t为时间；采用SST k-w湍流模型封闭方程中的雷诺应力与雷诺热通量；
柔性基材控制方程为：
柔性基材边界条件和初始条件为：
y|x＝-L/2＝y|x＝L/2＝0 (5)
y(t,0)＝y(t,L)＝0 (7)
m为柔性基材单位面积的质量，y和Vw表示柔性基材的垂直位移与柔性基材的轴向运动速度，F表示柔性基材中的张力；
步骤四、进行初始定常流场数值计算：在FLUENT中，对计算参数进行初始化；流场入口给定速度和高恒温边界条件，出口给定平均静压以及出口低恒温边界条件，柔性基材表面给定无滑移、光滑低恒温壁面边界条件；基于上述边界条件和初始条件，不考虑流场参数随时间的变化，不考虑柔性基材的变形，利用FLUENT中的SIMPLEC求解器进行订场流场数值计算，可以得到刚性壁面的流场区域的速度、压力以及温度分布；
步骤五、进行流固耦合数值计算：将步骤四中的定常流场数值计算结果作为初始条件进行流固耦合的数值计算。
2.根据权利要求1所述的带有柔性基材的悬浮喷嘴计算方法，其特征在于：所述流固耦合数值计算方法为：
步骤5.1、读取柔性基材壁面的压力，通过方程(4)计算出柔性基材的变形；
步骤5.2、根据步骤5.1所计算的变形结果，通过FLUENT中的扩散光顺方法对网格进行更新，得到新的柔性基材表面位置；
步骤5.3、根据新的柔性基材表面位置，计算流体域的速度、压力以及温度分布；
步骤5.4、重复步骤5.1至步骤5.3的计算过程，直至达到预定求解时间或者数值计算结果收敛，得到最终柔性基材的变形量以及最终的流场与热场分布。
3.根据权利要求1或2所述的带有柔性基材的悬浮喷嘴计算方法，其特征在于：提出一种基于余弦级数的柔性基材变形数值解法，在FLUENT中直接快速求解柔性基材的变形。
CN201810471784.4A 2018-05-17 2018-05-17 带有柔性基材的悬浮喷嘴计算方法 Active CN108682046B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810471784.4A CN108682046B (zh) 2018-05-17 2018-05-17 带有柔性基材的悬浮喷嘴计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810471784.4A CN108682046B (zh) 2018-05-17 2018-05-17 带有柔性基材的悬浮喷嘴计算方法

Publications (2)

Publication Number Publication Date
CN108682046A true CN108682046A (zh) 2018-10-19
CN108682046B CN108682046B (zh) 2020-08-11

Family

ID=63806444

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810471784.4A Active CN108682046B (zh) 2018-05-17 2018-05-17 带有柔性基材的悬浮喷嘴计算方法

Country Status (1)

Country Link
CN (1) CN108682046B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113723018A (zh) * 2021-07-28 2021-11-30 广东省科学院智能制造研究所 一种用于气浮风刀设计的流固耦合计算方法及装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5923329A (en) * 1996-06-24 1999-07-13 National Research Council Of Canada Method of grid generation about or within a 3 dimensional object
CN103514344A (zh) * 2012-06-21 2014-01-15 中南大学 一种基于谱方法的铝合金热轧板带横向厚度分布建模方法
CN104298869A (zh) * 2014-10-07 2015-01-21 北京理工大学 一种弹性水翼的流固耦合特性数值预测方法
CN107832260A (zh) * 2017-12-07 2018-03-23 大连理工大学 一种平板冲击射流传热问题的数值模拟方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5923329A (en) * 1996-06-24 1999-07-13 National Research Council Of Canada Method of grid generation about or within a 3 dimensional object
CN103514344A (zh) * 2012-06-21 2014-01-15 中南大学 一种基于谱方法的铝合金热轧板带横向厚度分布建模方法
CN104298869A (zh) * 2014-10-07 2015-01-21 北京理工大学 一种弹性水翼的流固耦合特性数值预测方法
CN107832260A (zh) * 2017-12-07 2018-03-23 大连理工大学 一种平板冲击射流传热问题的数值模拟方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
TAREK JAMALEDDINE AND MADHUMITA B. RAY: "Application of Computational Fluid Dynamics for Simulation", 《DRYING TECHNOLOGY》 *
张海燕: "涂布机悬浮烘箱的结构参数优化", 《机械设计与研究》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113723018A (zh) * 2021-07-28 2021-11-30 广东省科学院智能制造研究所 一种用于气浮风刀设计的流固耦合计算方法及装置

Also Published As

Publication number Publication date
CN108682046B (zh) 2020-08-11

Similar Documents

Publication Publication Date Title
Wang et al. Flow control on the NREL S809 wind turbine airfoil using vortex generators
Afshari et al. On numerical methods; optimization of CFD solution to evaluate fluid flow around a sample object at low Re numbers
Mahon et al. Computational analysis of pressure and wake characteristics of an aerofoil in ground effect
CN109530648A (zh) 一种预测连铸生产中结晶器钢渣界面波动的方法
CN108280259B (zh) 一种风电叶片涡流发生器安装位置及其尺寸的确定方法
CN109444350A (zh) 一种基于无人机的大气污染物监测传感器的布局方法
CN109766627A (zh) 一种基于滑板间距的受电弓非定常特性的分析方法
CN107832260B (zh) 一种平板冲击射流传热问题的数值模拟方法
CN111209703A (zh) 一种计及延迟的区域蒸汽热网拓扑结构优化方法及系统
CN112883667A (zh) 一种基于大涡模拟在竖向脉动风荷载作用下的膜结构抗风优化方法
Zou et al. Large-eddy simulation of flow around cylinder arrays at a subcritical Reynolds number
CN110543677B (zh) 一种涡特征驱动的旋转湍流pans模型的构建方法
Wang Effects of Reynolds number on vortex structure behind a surface-mounted finite square cylinder with AR= 7
CN110175377A (zh) 一种关于水翼云状空化的数值预报方法
CN112163381A (zh) 一种适用于复杂地形风场流动数值模拟的侧向边界条件设置方法
Zhao et al. A case study of hydraulic collecting a single spherical particle
CN108682046A (zh) 带有柔性基材的悬浮喷嘴计算方法
Song et al. Skin friction reduction characteristics of nonsmooth surfaces inspired by the shapes of barchan dunes
Yue et al. Experimental and numerical methods for obtaining flow field formed by hydraulic nodule pick-up devices
CN109492280A (zh) 一种用于双面狭缝涂布接触式模头的模拟计算方法
Ma et al. A Partially-Averaged Navier-Stokes model for hill and curved duct flow
Zhang et al. Experimental and numerical investigations of aerodynamic loads and 3D flow over non‐rotating MEXICO blades
CN115455688A (zh) 一种感应加热下中间包内夹杂物粒子去除效果预报方法
CN108052692A (zh) 垃圾在气力输运管线中运移规律的计算方法
CN108536896A (zh) 一种光伏逆变器内温度场计算方法

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
GR01 Patent grant
GR01 Patent grant