CN112800588B - 一种用于瞬态行驶工况下机舱换热器进风量的仿真计算方法 - Google Patents
一种用于瞬态行驶工况下机舱换热器进风量的仿真计算方法 Download PDFInfo
- Publication number
- CN112800588B CN112800588B CN202110018783.6A CN202110018783A CN112800588B CN 112800588 B CN112800588 B CN 112800588B CN 202110018783 A CN202110018783 A CN 202110018783A CN 112800588 B CN112800588 B CN 112800588B
- Authority
- CN
- China
- Prior art keywords
- vehicle speed
- heat exchanger
- working condition
- intake
- simulation calculation
- 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
Images
Classifications
-
- 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
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/08—Thermal analysis or thermal optimisation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Air-Conditioning For Vehicles (AREA)
Abstract
本发明公开了一种用于瞬态行驶工况下机舱换热器进风量的仿真计算方法,包括:确定采样点ai,对各采样点ai进行稳态仿真计算,输出各采样点ai仿真进风量的统计平均值Q1;将所述采样点ai与对应的进风量平均值Q1进行拟合,得到车速与进风量的关系式选取测试点bi,将测试点bi带入所述关系式中,得到测试点计算进风量的统计平均值Q2;对测试点bi进行稳态仿真计算,得到测试点仿真进风量的统计平均值Q2b;当时,确定为精准关系式;当时,重新规划采样点ai;将瞬态行驶工况下的车速导入确定的精准关系式,得到瞬态工况下机舱换热器进风量的关系曲线;其中,机舱换热器进风量包括:散热器的进风量和冷凝器的进风量。通过将瞬态过程离散成多个稳态计算过程从而预测瞬态行驶工况下。
Description
技术领域
本发明涉及一种用于瞬态行驶工况下机舱换热器进风量的仿真计算方法,属于汽车热管理领域。
背景技术
瞬态行驶工况下,车速随时间不断的变化带来的是换热器模块进风量的不断变化。常规的瞬态计算,需要将总时长离散成多个物理时长,每个物理时长内需要设置多个迭代步数。以工况总时长1800s为例,时间步长通常取3e-4s(按库朗特数为1进行计算),内迭代步数为50,采用256核并行计算,每一步长计算时间大概为6s,则总计算时长为1.8e9s,即20833天。很显然,这种规模的仿真计算无法实现。目前国内外关于机舱换热器模块进风量仿真计算的研究大多集中在稳态工况方面,对于瞬态行驶工况的仿真则较少,尚无明确的仿真方法。因此,本发明旨在解决瞬态行驶工况换热器模块进风量仿真计算量大、成本高、周期长、无法为热管理性能一维仿真提供换热器模块进气量输入这一问题。
发明内容
本发明设计开发了一种用于瞬态行驶工况下机舱换热器进风量的仿真计算方法,通过将瞬态过程离散成多个稳态计算过程,找到车速与换热器模块进风量的对应关系,从而预测瞬态行驶工况下,发动机舱换热器的进风量。
本发明提供的技术方案为:
一种用于瞬态行驶工况下机舱换热器进风量的仿真计算方法,包括:
在瞬态行驶工况下,确定采样点ai,对各采样点ai进行稳态仿真计算,输出各采样点ai仿真进风量平均值Q1;
将瞬态行驶工况下的车速导入确定的精准关系式,得到瞬态工况下机舱换热器进风量的关系曲线;
其中,机舱换热器进风量包括:散热器的进风量和冷凝器的进风量。
优选的是,所述采样点的确定过程包括:
在瞬态行驶工况下,识别出最高车速和最低车速,以最低车速作为第一个采样点,最高车速作为最后一个采样点,中间每隔20kph确定为一个采样点;
当重新规划样本点时,减小样本点间隔,每隔15kph确定为一个样本点。
优选的是,还包括:
在瞬态行驶工况时间为1800s,最低车速为0kph,最高车速为120kph条件下,选取采样点进行稳态仿真计算,得到采样点车速下散热器和冷凝器进风量平均值。
优选的是,
选用4次多项式进行拟合,得到车速与冷凝器的进风量关系式:
f(xal)=pl1x^4+pl2x^3+pl3x^2+pl4x+pl5,
其中,pl1=3.321e-09、pl2=-6.776e-07、pl3=-7.608e-07,pl4=0.001228、pl5=-0.516;
车速与散热器的进风量关系式为:
f(xas)=ps1x^4+ps2x^3+ps3x^2+ps4x+ps5,
其中,ps1=5.328e-10、ps2=-1.974e-07、ps3=3.351e-05,ps4=0.0002911、ps5=0.1636。
优选的是,
在最低车速为10kph,最高车速为110kph条件下,选取测试点带入上述关系式中计算测试点的进风量,对测试点进行稳态仿真计算,经过误差对比后,确定上述关系式为精准关系式。
本发明所述的有益效果:本发明提供的用于瞬态行驶工况下机舱换热器进风量的仿真计算方法,与特定车速稳态计算及瞬态计算相比,本发明即可实现瞬态行驶工况换热器模块进风量的预测,又可大大降低计算量及计算时长,可为一维热管理性能仿真提供瞬态行驶工况下换热器进风量输入,且满足工程开发时间节点的要求。
附图说明
图1为本发明所述的瞬态行驶工况换热器模块进风量计算技术路线图。
图2为本发明所述的车速与散热器进风量关系拟合图。
图3为本发明所述的车速与冷凝器进风量关系拟合图。
具体实施方式
下面结合附图对本发明做进一步的详细说明,以令本领域技术人员参照说明书文字能够据以实施。
如图1-3所示,本发明提供一种用于瞬态行驶工况下机舱换热器进风量的仿真计算方法,通过将瞬态过程离散成多个稳态计算过程,找到车速与换热器模块进风量的对应关系,从而预测瞬态行驶工况下,换热器模块的进风量,具体包括:
在瞬态行驶工况下,确定采样点ai,对各采样点ai进行稳态仿真计算,输出在各采样点ai进风量平均值Q1;
将瞬态行驶工况下的车速导入精准关系式,得到瞬态工况下机舱换热器进风量的关系曲线;
其中,机舱换热器进风量包括:散热器的进风量和冷凝器的进风量。
其中,进风量平均值为进风量计算稳定段的统计平均值。
对瞬态行驶工况进行最高车速及最低车速识别,以最低车速作为第一个样本点、最高车速作为最后一个样本点,中间每隔20kph确定为一个样本点。
对采样点ai应用starccm+软件进行流场稳态计算,入口条件选取各测试点车速,出口为压力,设置机舱换热器进风量统计监测曲线,仿真计算后输出各采样点车速下换热器模块进风量的统计平均值Q1,其中,starccm+软件为2020.1.1版本,在windows10***运行。
选取采样点以外的测试点车速带入公式进行换热器模块进风量计算,输出测试点公式计算换热器模进风量Q2。
选取样本点以外的测试点车速bi进行稳态仿真计算,入口条件选取各测试点车速bi,出口为压力,设置换热器模块进风量统计监测曲线,仿真计算后输出测试点稳态仿真计算换热器模块进风量的统计平均值Q2b。
如图1所示,将测试点公式计算得到的换热器模块进风量Q2与测试点稳态仿真计算换热器模块进风量Q2b进行对比,若误差小于1.5%,则输出精准车速与换热器模块进风量关系式。若误差大于1.5%,则重新规划样本点,中间每隔15kph确定为一个样本点,重复上述过程。直至,测试点公式计算换热器进风量与测试点稳态仿真计算换热器进风量对比,误差小于1.5%为止。
实施例
如图2-3所示,在瞬态行驶工况时间为1800s,最低车速为0kph,最高车速为120kph条件下,选取0kph/20kph/40kph/60kph/80kph/100kph/120kph作为采样点车速进行稳态仿真计算,其中,对采样点应用ai应用starccm+软件(2020.1.1版本)在Windows10***运行进行流场稳态计算,物理模型选取Three dimensional、Steady、Gas、Segragated Flow、Turbulent、Reynalds-Averaged Navies-Stokes、K-Epsilon Turbulence、Realizable K-Epsilon Turbulence、Two-Layer All y+wall Treatment。车身造型面及底部外表面边界类型为Wall,入口条件选取各测试点车速,出口压力设置为0pa,设置机舱换热器进风量统计监测曲线,仿真计算后输出各采样点车速下换热器模块进风量平均值,应用Matalab软件(2017b版本)在Windows10***运行,选用4次多项式对车速与换热器模块进风量关系进行数据拟合,得到车速与散热器及冷凝器进风量的关系式:
车速与冷凝器的进风量关系式:
f(xal)=pl1x^4+pl2x^3+pl3x^2+pl4x+pl5,
其中,pl1=3.321e-09、pl2=-6.776e-07、pl3=-7.608e-07,pl4=0.001228、pl5=-0.516;
车速与散热器的进风量关系式为:
f(xas)=ps1x^4+ps2x^3+ps3x^2+ps4x+ps5,
其中,ps1=5.328e-10、ps2=-1.974e-07、ps3=3.351e-05,ps4=0.0002911、ps5=0.1636。
选取样本点以外的测试点10kph/30kph/50kph/70kph/90kph/110kph车速进行测试,测试点公式计算结果与稳态仿真计算结果相比,散热器进风量误差为分别为0.11%/0.08%/0.105%0.13%/0.14%0.09%,冷凝器进风量误差为0.29%/0.36%/0.31%0.34%/0.28%/0.32%,满足误差小于1.5%的要求,可以用于瞬态行驶工况换热器模块进风量预测。
将瞬态行驶工况随时间变化的车速导入精准的车速与换热器模块进风量关系式,即可得到随时间变化的换热器模块进风量关系曲线。
尽管本发明的实施方案已公开如上,但其并不仅仅限于说明书和实施方式中所列运用,它完全可以被适用于各种适合本发明的领域,对于熟悉本领域的人员而言,可容易地实现另外的修改,因此在不背离权利要求及等同范围所限定的一般概念下,本发明并不限于特定的细节和这里示出与描述的图例。
Claims (5)
1.一种用于瞬态行驶工况下机舱换热器进风量的仿真计算方法,其特征在于,包括:
在瞬态行驶工况下,确定采样点ai,对各采样点ai进行稳态仿真计算,输出各采样点ai仿真进风量平均值Q1;
将瞬态行驶工况下的车速导入确定的精准关系式,得到瞬态工况下机舱换热器进风量的关系曲线;
其中,机舱换热器进风量包括:散热器的进风量和冷凝器的进风量。
2.根据权利要求1所述的用于瞬态行驶工况下机舱换热器进风量的仿真计算方法,其特征在于,所述采样点的确定过程包括:
在瞬态行驶工况下,识别出最高车速和最低车速,以最低车速作为第一个采样点,最高车速作为最后一个采样点,中间每隔20kph确定为一个采样点;
当重新规划样本点时,减小样本点间隔,每隔15kph确定为一个样本点。
3.根据权利要求2所述的用于瞬态行驶工况下机舱换热器进风量的仿真计算方法,其特征在于,还包括:
在瞬态行驶工况时间为1800s,最低车速为0kph,最高车速为120kph条件下,选取采样点进行稳态仿真计算,得到采样点车速下散热器和冷凝器进风量平均值。
4.根据权利要求3所述的用于瞬态行驶工况下机舱换热器进风量的仿真计算方法,其特征在于,
选用4次多项式进行拟合,得到车速与冷凝器的进风量关系式:
f(xal)=pl1x^4+pl2x^3+pl3x^2+pl4x+pl5,
其中,pl1=3.321e-09、pl2=-6.776e-07、pl3=-7.608e-07,pl4=0.001228、pl5=-0.516;
车速与散热器的进风量关系式为:
f(xas)=ps1x^4+ps2x^3+ps3x^2+ps4x+ps5,
其中,ps1=5.328e-10、ps2=-1.974e-07、ps3=3.351e-05,ps4=0.0002911、ps5=0.1636。
5.根据权利要求4所述的用于瞬态行驶工况下机舱换热器进风量的仿真计算方法,其特征在于,
在最低车速为10kph,最高车速为110kph条件下,选取测试点带入上述关系式中计算测试点的进风量,对测试点进行稳态仿真计算,经过误差对比后,确定上述关系式为精准关系式。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110018783.6A CN112800588B (zh) | 2021-01-07 | 2021-01-07 | 一种用于瞬态行驶工况下机舱换热器进风量的仿真计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110018783.6A CN112800588B (zh) | 2021-01-07 | 2021-01-07 | 一种用于瞬态行驶工况下机舱换热器进风量的仿真计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112800588A CN112800588A (zh) | 2021-05-14 |
CN112800588B true CN112800588B (zh) | 2022-04-01 |
Family
ID=75808988
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110018783.6A Active CN112800588B (zh) | 2021-01-07 | 2021-01-07 | 一种用于瞬态行驶工况下机舱换热器进风量的仿真计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112800588B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113588288B (zh) * | 2021-07-30 | 2022-06-14 | 吉林大学 | 一种用于汽车台架测试的瞬态供风模拟***及控制方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008008882A (ja) * | 2006-06-01 | 2008-01-17 | Yokohama Rubber Co Ltd:The | タイヤのコーナリング特性の評価方法および装置 |
CN102436185A (zh) * | 2011-11-08 | 2012-05-02 | 浙江理工大学 | 一种热定型机的换热器建模方法 |
CN105956316A (zh) * | 2016-05-18 | 2016-09-21 | 奇瑞汽车股份有限公司 | 一种汽车发动机舱前端模块空气进风量计算方法 |
CN107391807A (zh) * | 2017-06-28 | 2017-11-24 | 西安交通大学 | 基于瞬态技术的板翅式换热器传热流动性能数值模拟方法 |
CN109339933A (zh) * | 2018-10-31 | 2019-02-15 | 吉林大学 | 非最小相位发动机冷却***水温控制方法 |
CN109858073A (zh) * | 2018-12-12 | 2019-06-07 | 北京航空航天大学 | 针对板翅式换热器效率模型的基于瞬态试验数据的实时快速构建方法 |
CN111125858A (zh) * | 2018-11-01 | 2020-05-08 | 广州汽车集团股份有限公司 | 一种整车冷却***的仿真计算方法及*** |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20160018252A1 (en) * | 2013-03-04 | 2016-01-21 | International Road Dynamics, Inc. | Sensor including electrical transmission-line parameter that changes responsive to vehicular load |
US20170087497A1 (en) * | 2015-09-29 | 2017-03-30 | Chevron U.S.A. Inc. | Inlet distributor device and methods for use and design thereof |
CN106840685B (zh) * | 2017-01-22 | 2019-01-04 | 中国科学院工程热物理研究所 | 一种燃气轮机动态工况数据库样本点选取方法 |
US10817033B2 (en) * | 2017-12-14 | 2020-10-27 | Schneider Electric It Corporation | Method and system for predicting effect of a transient event on a data center |
CN109115516B (zh) * | 2018-08-02 | 2020-01-10 | 安徽江淮汽车集团股份有限公司 | 一种汽车乘员舱新风热负荷测量方法 |
-
2021
- 2021-01-07 CN CN202110018783.6A patent/CN112800588B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2008008882A (ja) * | 2006-06-01 | 2008-01-17 | Yokohama Rubber Co Ltd:The | タイヤのコーナリング特性の評価方法および装置 |
CN102436185A (zh) * | 2011-11-08 | 2012-05-02 | 浙江理工大学 | 一种热定型机的换热器建模方法 |
CN105956316A (zh) * | 2016-05-18 | 2016-09-21 | 奇瑞汽车股份有限公司 | 一种汽车发动机舱前端模块空气进风量计算方法 |
CN107391807A (zh) * | 2017-06-28 | 2017-11-24 | 西安交通大学 | 基于瞬态技术的板翅式换热器传热流动性能数值模拟方法 |
CN109339933A (zh) * | 2018-10-31 | 2019-02-15 | 吉林大学 | 非最小相位发动机冷却***水温控制方法 |
CN111125858A (zh) * | 2018-11-01 | 2020-05-08 | 广州汽车集团股份有限公司 | 一种整车冷却***的仿真计算方法及*** |
CN109858073A (zh) * | 2018-12-12 | 2019-06-07 | 北京航空航天大学 | 针对板翅式换热器效率模型的基于瞬态试验数据的实时快速构建方法 |
Non-Patent Citations (5)
Title |
---|
A numerical simulation of wheel spray for simplified vehicle model based discrete phase method;Xingjun Hu等;《Advances in Mechanical Engineering》;20150730;第7卷(第7期);第1-8页 * |
Cross Flow Heat Exchanger Modeling of Transient Temperature Input Coditions;Tianyi Gao等;《IEEE Transactions on Components,Packaging and Manufacturing Technology》;20140919;第4卷(第11期);第1796-1807页 * |
某车型发动机舱稳态和瞬态的CFD分析;郭健忠等;《科技通报》;20170930;第33卷(第09期);第199-202页 * |
汽车侧窗风振噪声特性研究;胡兴军等;《华南理工大学学报(自然科学版)》;20201115;第48卷(第11期);第64-71页 * |
汽车发动机冷却风扇对机舱热管理影响的研究;李飞等;《汽车实用技术》;20190730;第2019卷(第14期);第95-98页 * |
Also Published As
Publication number | Publication date |
---|---|
CN112800588A (zh) | 2021-05-14 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108304684B (zh) | 一种火箭发动机尾喷射流仿真方法及*** | |
CN107977513B (zh) | 一种基于路径搜索的天然气动态流动温度预测方法 | |
CN112800588B (zh) | 一种用于瞬态行驶工况下机舱换热器进风量的仿真计算方法 | |
CN112632719B (zh) | 一种基于一维平均流线法的多级轴流压气机特性校正方法 | |
CN114254460B (zh) | 基于图卷积神经网络的透平机械气动鲁棒性优化方法 | |
CN115292946B (zh) | 一种基于变比热计算的高压涡轮效率评估方法及装置 | |
Brown et al. | Surrogate modeling of manufacturing variation effects on unsteady interactions in a transonic turbine | |
Watanabe et al. | The CFD application for efficient designing in the automotive engineering | |
CN111859667B (zh) | 一种预测汽车空调冷凝器性能的建模方法 | |
Guardiola et al. | Representation limits of mean value engine models | |
Mehdipour et al. | Mathematical simulation of a vehicle radiator by genetic algorithm method and comparison with experimental data | |
CN113343607B (zh) | 一种基于管壁液膜运动规律的携液临界气流速计算方法 | |
Chen et al. | An effective turbine blade parameterization and aerodynamic optimization procedure using an improved response surface method | |
CN113722900B (zh) | 一种航空换热器非设计点性能设计与分析方法 | |
CN116822299B (zh) | 一种服役载荷历程下航空发动机火焰筒热应力快速计算方法 | |
Stephenson et al. | Multi-objective optimization of a charge air cooler using modeFRONTIER | |
Lafferty et al. | Under-hood Thermal Simulation of a Class 8 Truck | |
Bhatnagar et al. | Simulation of Class 8 Truck Cooling System: Comparison to Experiment under Different Engine Operation Conditions | |
CN118424035A (zh) | 用于多管程绕管式换热器的高效换热控制方法 | |
CN106096112A (zh) | 液力变矩器流场计算方法 | |
CN114781069A (zh) | 一种基于整车热平衡试验数据库开发冷却***性能分析工具的构建方法以及分析方法 | |
Hopf | CFD Topology and Shape Optimization of Ford Applications using Tosca Fluid | |
CN116976002A (zh) | 一种发动机机舱流场的评价方法及*** | |
Rauseo et al. | PHYSICS GUIDED MACHINE LEARNING MODELLING OF COMPRESSOR AEROELASTIC FLUTTER | |
CN118133604A (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 | ||
CB03 | Change of inventor or designer information | ||
CB03 | Change of inventor or designer information |
Inventor after: Fu Yu Inventor after: Min Haitao Inventor after: Hu Xingjun Inventor after: Gao Bingzhao Inventor after: Ma Xiaole Inventor before: Fu Yu Inventor before: Hu Xingjun Inventor before: Gao Bingzhao Inventor before: Ma Xiaole |
|
GR01 | Patent grant | ||
GR01 | Patent grant |