J1-08 座舱一维多节点热模型与参数辨识(含 CFD→1D 分区映射)
课程代码 J1-08 | 板块 J 仿真分析 CAE / J1 一维系统仿真 时长 约 4 学时 | 难度 L3 专家 适合对象 座舱 / 整车热管理仿真工程师(必修);控制与标定工程师了解模型接口、状态向量与可信边界(选修);试验工程师了解辨识所需的激励设计与取数纪律(选修) 前置 J1-01 1D 热管理仿真导论(KULI/AMESim/GT/Flowmaster) 1D 热管理仿真导论;A1-06 热阻网络法与集总参数建模入门 热阻网络法与集总参数建模入门;A2-04 湿空气物理与焓湿图:空调负荷计算基石 湿空气物理与焓湿图;B1-04 目标设定:降温速率、冬季续航热损、快充温升 目标设定:降温速率、冬季续航热损、快充温升;C1-06 乘用车座舱制冷负荷计算书:逐项取值、UA 标定与环境舱对账 乘用车座舱制冷负荷计算书:逐项取值、UA 标定与环境舱对账 本讲义定位 讲清座舱这一个对象的多节点 1D 模型怎么定结构、怎么把参数从实测里解出来、怎么把 CFD 的分区结果喂回来、以及交付时必须随行的可信边界。它不是建模基础课(热阻网络、集总参数判据、RC 拓扑在 A1-06 热阻网络法与集总参数建模入门),不是负荷计算课(分项取值与来源列在 C1-06 乘用车座舱制冷负荷计算书:逐项取值、UA 标定与环境舱对账),不是三维课(座舱 CFD 建模与提取在 J4-01 乘员舱气流组织与温度场 CFD),也不交任何控制器实现。
引言:模型交出去的那一刻,误差就不再由你控制了
座舱是乘用车最大的单一热负荷对象,而在一维系统仿真里,它长期是被简化得最狠的那一块——制冷剂回路有专门的 1D 建模课,冷却液回路有,电池有,唯独座舱常常在整车模型里只剩一个「等效热容 + 等效 UA」的方块。这个方块在它被标定的那一条曲线上通常表现很好,问题出在它被拿去做别的事的时候。
本课要立的判断力只有一条主线:一个座舱模型,凭什么可以被下游拿去用。
它要过四道关,正好是本课的四块内容。
第一道关:结构本身要可辨识。节点划分不是「物理上分得越细越像」的问题。多加一个节点就多几个参数,而参数能不能被解出来,上界由你有几路独立测点、几段独立激励决定,不由物理细节决定。这一层的常见错误不是分得太粗,而是分得太细却没有对应的观测量——模型在标定工况上照样能拟合,误差只是在参数之间乱窜,一换工况就崩。第 1 讲处理这一层。
第二道关:边界项要写成能随工况变的形式。座舱的三条主要外部输入——玻璃太阳得热、不透明车身得热、渗漏换气——每一条都有一个「取常数」的诱人写法,而每一条取常数的偏差方向都是偏乐观。第 2 讲把这三条各自写成 1D 能吃、又能随辐照 / 车速变的形式。
第三道关:参数要真的被解出来,不是被优化器分配出来。这是本课最专属的一层。等效 UA 与有效太阳开口在暴晒升温段同向抬温,一条升温曲线只能定出它们的组合量;此时继续优化,得到的是一组数值上收敛、物理上互相抵消的参数。第 3 讲给出的处置不是「换更好的优化器」,是换激励,而且要在动手拟合之前就用灵敏度矩阵判断这组数据够不够。
第四道关:交付时必须带边界。没有可信边界与失效工况声明的模型,下游必然会拿它去跑不对称受晒、强直吹与开门瞬态——这三类 1D 本就不成立,而且偏差方向偏乐观,最容易被当成结论用。第 6 讲把这一页写成交付物的一部分。
中间还夹着一层本课独有的活:CFD 的结果怎么反向喂回 1D(第 4 讲),以及湿度这条状态怎么加进来而不把起雾掩盖掉(第 5 讲)。
本课不给什么,先说清楚。① 建模基础——热阻网络与集总参数拓扑(Cauer / Foster)、电-热类比、Bi 判据本体、有效参与质量,以及最小二乘反问题辨识方法学本体,全部承 A1-06 热阻网络法与集总参数建模入门;1D 建模范式与工具口径承 J1-01 1D 热管理仿真导论(KULI/AMESim/GT/Flowmaster)。② 单节点一阶模型与 pull-down 目标可行性预检承 B1-04 目标设定:降温速率、冬季续航热损、快充温升。③ 负荷分项口径、分部位 U 值取值来源列、内饰有效热容清单法、逐面式太阳得热、乘用车渗漏量级与来源、以及「先对边界、再对分项、最后才动参数」的残差归因顺序,全部承 C1-06 乘用车座舱制冷负荷计算书:逐项取值、UA 标定与环境舱对账。④ 座舱 CFD 的建模有效性(浮力判据、DO / S2S 取舍、y⁺ 与网格无关性)、sol-air 式本体、换气次数的两种口径与分区量的提取入口,全部承 J4-01 乘员舱气流组织与温度场 CFD;本课出的是映射方法学、加权口径与回验。⑤ 1D-3D 在线双向联算、宏时间步与耦合发散承 J2-02 机舱热管理与前端模块 CFD,本课是离线单向提参数。⑥ 人体侧的分区换热系数与等效温度承 A4-03 等效温度 EHT 与动态/瞬态热舒适评价,其对象是人体表面,与本课的舱内壁面换热系数不可互换。⑦ 降阶与相关性验证承 J7-02 模型降阶(ROM)与控制导向建模 与 J7-01 仿真—试验相关性验证(Correlation)。⑧ 控制器侧全部不在本课:玻璃温度估算器的输入裁剪与冷点偏置见 C5-03 低温起雾机理与主动防雾控制,产湿率在线推断的可辨识性边界见 C4-04 座舱湿度管理与舒适加湿,观测器整定与退出条件见 K3-10 热管理在线状态与负荷估计:观测器整定、软测量置信度与退出条件,分区送风策略见 A4-08 局部微气候与分区送风的舒适策略 与 C6-01 双区/四区自动空调控制逻辑。⑨ 试验规程本身承 L3-03 整车降温/升温性能主客观评价 与 L4-01 高温标定试验(吐鲁番等)方法。
数字上分三类,全课一致。物性与可复算的常数给准确值并写明来源级别;行业典型区间只给区间并标「典型」;与平台相关的量(等效 UA、有效热容、有效太阳开口、泄漏系数、舱容、分区换热系数)本课给出典型参考带与常见值锚点,用途是让你判断「自己辨识出来的数落在不在合理量级」,⚠ 它们是典型值、⛔ 不是验收依据,必须以本平台实测确认;拿到它们的路径在各自小节里写明。凡带代入值的算例,输入一律是教学假设值,⛔ 不对应任何平台、⛔ 禁止照抄取用。所引标准本课未逐字核对原文,故只写标准号与它管什么对象,⛔ 不写条款内容与限值;版次以现行有效目录为准,引用前须核。
符号约定,先立后用。本课横跨传热、光学、湿空气与辨识四套记号,而这四套在本课程体系里已各有属主课钉过约定。为了让读者能把本课的数直接搬到那几门课里去用,本课不另起一套,逐条沿用并在下表标明来源;三处必须换符号的地方也在表里说明理由。
| 符号 | 含义 | 单位 | 口径来源与说明 |
|---|---|---|---|
| T_i、C_i、R_ij | 节点温度、节点有效热容、节点间热阻 | ℃、J/K、K/W | RC 网络记法承 A1-06 热阻网络法与集总参数建模入门 |
| UA_id | 辨识出的座舱等效传热量 | W/K | 与 C1-06 乘用车座舱制冷负荷计算书:逐项取值、UA 标定与环境舱对账 的 UA_eff 是同一个物理量、不是同一个数,两者之别见 3.6 |
| (m·c_p)_id | 辨识出的座舱有效热容 | kJ/K | 与 C1-06 乘用车座舱制冷负荷计算书:逐项取值、UA 标定与环境舱对账 清单法给的 (m·c_p)_eff 同上,见 3.6 |
| τ_cab | 座舱时间常数 = (m·c_p)/UA | s | 承 C1-06 乘用车座舱制冷负荷计算书:逐项取值、UA 标定与环境舱对账 与 B1-04 目标设定:降温速率、冬季续航热损、快充温升;⛔ 全课不裸写 τ |
| (τA)_eff | 有效太阳开口 | m² | 本课首定,见 2.3 |
| τ_e | 玻璃太阳直接透射比 | 无量纲 | 承 J4-01 乘员舱气流组织与温度场 CFD;⚠ C1-06 乘用车座舱制冷负荷计算书:逐项取值、UA 标定与环境舱对账 写作 τ_g,是同一个量,见 2.2 |
| α_s | 短波(太阳)吸收比 | 无量纲 | 承 J4-01 乘员舱气流组织与温度场 CFD / C1-06 乘用车座舱制冷负荷计算书:逐项取值、UA 标定与环境舱对账;⛔ 不裸写 α(裸写的 α 是热扩散率) |
| ε_IR | 长波红外发射率 | 无量纲 | 承 J4-01 乘员舱气流组织与温度场 CFD;⛔ 不裸写 ε |
| h_o、h_i | 外 / 内表面综合换热系数(对流 + 已线性化的长波) | W/(m²·K) | 承 C1-06 乘用车座舱制冷负荷计算书:逐项取值、UA 标定与环境舱对账 1.7 与 J4-01 乘员舱气流组织与温度场 CFD 1.4;⛔ 不裸写 h |
| h_zone,k | 第 k 个分区的舱内壁面换热系数 | W/(m²·K) | 本课首定;⛔ 与 A4-03 等效温度 EHT 与动态/瞬态热舒适评价 的人体表面分区值不可互换,见 4.6 |
| ṁ_f、ṁ_inf、ṁ_air | 新风、渗漏、送风的质量流量 | kg/s | 承 C4-04 座舱湿度管理与舒适加湿;⚠ C5-03 低温起雾机理与主动防雾控制 把新风写作 ṁ_fresh,同一个量 |
| ACH_inf | 渗漏换气次数(分母是舱容 V_cab) | 次/h | 承 C1-06 乘用车座舱制冷负荷计算书:逐项取值、UA 标定与环境舱对账 4.3 |
| V_cab | 舱容 | m³ | 承 C1-06 乘用车座舱制冷负荷计算书:逐项取值、UA 标定与环境舱对账;⚠ 必须与实测 ACH_inf 时用的舱容同源同口径 |
| d、d_amb、d_sat | 舱内 / 外气 / 饱和含湿量(分母是干空气) | g/kg 干空气 | 承 A2-04 湿空气物理与焓湿图:空调负荷计算基石 与 C4-04 座舱湿度管理与舒适加湿 |
| m_da | 舱内干空气质量 = V_cab·ρ_da | kg 干空气 | 承 C4-04 座舱湿度管理与舒适加湿 2.1 |
| ṁ_gen、ṁ_dehum、ṁ_cond | 产湿、蒸发器除湿、冷表面析水 | g/h | 承 C4-04 座舱湿度管理与舒适加湿;⛔ 后两者不得合并,见 5.2 |
| T_g,i | 玻璃内表面温度 | ℃ | 承 C5-03 低温起雾机理与主动防雾控制 1.2;⚠ 口语里的 T_glass 是它的简写,本课一律带下标 |
| T_dew | 舱内空气露点 | ℃ | 承 C5-03 低温起雾机理与主动防雾控制;⚠ C4-04 座舱湿度管理与舒适加湿 写作 T_dp,同一个量 |
| ΔT_margin | 露点裕度 = T_g,i − T_dew | K | 承 C5-03 低温起雾机理与主动防雾控制,本课只输出这条曲线、⛔ 不定阈值 |
| S、r_jk | 灵敏度矩阵、第 j 与第 k 个参数的灵敏度相关系数 | — | 本课首定,见 3.4 |
⚠ 三处刻意与常见写法不同,理由必须先读一遍,否则第 3、5 讲会与兄弟课读串:
- 参数相关系数写 r_jk,⛔ 不写 ρ。统计学里相关系数惯用 ρ,但在本课程体系里裸写的 ρ 是材料 / 空气密度(C1-06 乘用车座舱制冷负荷计算书:逐项取值、UA 标定与环境舱对账 1.7 已钉死),而本课第 3 讲会在同一页里同时出现空气密度与相关系数。⇒ 相关系数一律写 r_jk。
- 条件数写全称 cond(SᵀS),⛔ 不引入单字母 κ。C4-04 座舱湿度管理与舒适加湿 3.5 已把 κ 用作「含内饰湿容的真实时间常数对空气时间常数的放大倍数」,而本课第 5 讲正要引用那个量。⇒ 本课不给条件数一个单字母符号。
- 玻璃节点温度写 T_g,i。C5-03 低温起雾机理与主动防雾控制 已把内 / 外表面分开钉死,而本课的玻璃节点温度正是它的输入——若本课写成裸的 T_glass,读者搬过去时无从判断这是内表面还是外表面,而两者在暴晒工况下可以差出十几开尔文。
第 1 讲 状态向量与节点划分:几个节点才够
这一讲不写任何一条方程的解,它只做一件事:把状态向量定下来。定完之后,你才知道要辨识哪几个参数、需要哪几路观测量、以及数据够不够。顺序反过来做——先建一个「看起来物理上很全」的模型再去凑数据——是本课要拆掉的第一个习惯。
1.1 节点表是交付物的第一页,三栏缺一即不可辨识
座舱多节点模型的状态向量,典型起步是这样一组:
- 空气节点:前排(后排是否单独分开,见 1.5);
- 固体节点:玻璃、顶棚、仪表台、座椅、地毯,必要时加门内板与前围。
但是「列出这几个名字」不叫定状态向量。每一个节点必须同时写清三栏:
| 栏 | 它回答什么 | 缺了会怎样 |
|---|---|---|
| 热容从哪来 | 这个 C_i 是清单法算的、手册物性算的,还是待辨识的 | 参数与常数混在一起,辨识时不知道该动哪个 |
| 边界是什么 | 它的外侧接谁(环境 / sol-air 等效环境 / 另一个固体节点),受不受晒、受不受直吹 | 两个边界不同的表面被合成一个节点,太阳项与围护项互相污染 |
| 对应哪个观测量 | 实测里哪一路热电偶 / 哪一个测点对得上它 | 这个节点的参数不可辨识——它可以取任何值,残差都差不多 |
第三栏是最常被空着的,也是最致命的。⇒ 本课的第一条可执行动作:节点表里第三栏为空的行,要么补测点,要么把这个节点合并掉。把它留在模型里而不去管,等于往参数空间里加了一个自由维度,优化器会把别处的误差塞进去。
1.2 Bi < 0.1 回答不了座舱的分节点问题
A1-06 热阻网络法与集总参数建模入门 给的集总参数判据是 Bi < 0.1。这条判据在本课仍然成立,但它回答的是另一个问题:这一块固体要不要沿厚度方向再分层。它管的是固体内部的导热梯度。
座舱要回答的是另一个问题:这两块东西能不能合成一个节点。而这个问题的判据在边界上,不在 Bi 上。
最干净的反例就是玻璃。停在阳光下的车,侧窗一块受晒、一块背阴:
- 两块玻璃材质、厚度完全相同,Bi 都远小于 0.1,各自内部确实近似等温;
- 但它们的外侧边界一个有 α_s·I 这一项、一个没有,稳态温差可以到十几开尔文量级;
- 合成一个节点之后,模型给出的是两者的平均值——而防雾判据要的是最冷的那一块(C5-03 低温起雾机理与主动防雾控制 3.x 的冷点问题),暴晒热浸要的是最热的那一块。两个下游要的都不是平均值。
⇒ 座舱的分节点判据是「边界不同就要分」,Bi 只管「要不要沿厚度分层」。两条判据是串联关系,不是替代关系:先按边界分节点,再对每个节点问一次 Bi 要不要分层。
1.3 座舱另立的三条节点数判据
除了「边界不同就要分」,座舱还需要三条判据来回答「拆到几个节点够用」。这三条按顺序问,任何一条答「是」就保留独立节点:
判据 ①:这个节点是否承担独立输出? 模型的每一路输出都必须有一个节点为它负责。玻璃内表面温度 T_g,i 要交给露点裕度(C5-03 低温起雾机理与主动防雾控制),内饰表面温度要交给触摸温度限,前排空气温度要交给舒适评价。⇒ 有独立输出的,节点必须独立存在,哪怕它的热容很小。
判据 ②:它的时间常数是否与已有节点相差一个数量级以上? 时间常数相近的两个节点,在任何一条实测曲线上都表现为同一个指数段,无论怎么拟合都分不开;此时合并不损失信息。相差一个数量级以上的必须分开——典型的例子是玻璃(薄、热容小、跟随快,分钟量级)与座椅泡沫(厚、热容大、跟随慢,几十分钟量级):合成一个节点之后,模型的降温曲线会既跟不上前段的快降,也跟不上后段的长尾。
判据 ③:它是否接受独立边界? 即 1.2 那一条的推广:受晒 / 背阴、直吹 / 非直吹、贴外板 / 贴内饰。直吹与非直吹这一对在 pull-down 前段尤其明显——直吹面的表面换热系数可以比非直吹面高数倍。
1.4 节点数与可辨识性是一对对偶关系
这是本讲最反直觉的一条,也是全课的地基:
节点数的上界由独立测点数与独立激励段数决定,⛔ 不由物理细节决定。
理由很直接。每加一个节点,至少多出一个未知的热容与一个未知的热阻;而你能从数据里提取的独立信息量,由「独立观测通道数 × 独立激励段数」封顶。物理上再合理的细分,一旦越过这条线,多出来的参数就进入了一个零空间:它们的任意组合都给出几乎相同的残差。
这件事在标定工况上看不出来——残差甚至会更小,因为多出来的自由度总能把噪声也拟合一点。它在换工况的时候才爆发,而那时模型已经交出去了。
⇒ 处置有三条,按优先级排:
- 先减节点——把没有独立输出、时间常数相近、边界相同的节点合并;
- 再加激励——多取一段辐照条件不同的数据,比多贴两个热电偶便宜(见第 3 讲);
- 最后才加测点——加在那个「有独立输出但没有测点」的节点上。
⚠ 三条之外还有一条不能做的:人为固定某个参数以「减少未知数」。它不解决问题,只是把不可辨识性藏起来——固定值若不对,误差会全部转移到其余参数上,而残差看起来正常。第 3 讲 3.4 会给出为什么它一定会失败。
1.5 前后排要不要分区,看两件事而不是一件
「座舱模型要不要分前后排」在项目上常常被当成一个精度问题来争,其实它是两个独立条件的合取:
- 条件 A(出口需求):模型的输出里要不要出现前后排温差,或者要不要支撑分区控制的评估?
- 条件 B(观测条件):后排有没有独立测点?
两条都成立,才做双空气节点。只有 A 没有 B 时,正确做法是单空气节点 + 后排偏置——即用一条固定或随工况查表的偏置量表示后排与前排的差,偏置量由实测统计得到而不是由模型解出。这样做诚实地承认了「这个差值不是本模型算出来的」,下游也就不会拿它去做分区控制的定量结论。
⚠ 边界要说清楚:分区送风策略本身、以及双区 / 四区控制逻辑,不在本课——去向 A4-08 局部微气候与分区送风的舒适策略 与 C6-01 双区/四区自动空调控制逻辑。本课只回答「模型的状态向量里要不要有第二个空气节点」。
1.6 本讲不重讲的部分与它的出口
不重讲:单节点一阶模型 τ_cab = (m·c_p)_eff/(UA + ṁ_air·c_p) 与 pull-down 目标可行性预检,承 B1-04 目标设定:降温速率、冬季续航热损、快充温升 第 2 讲;RC 网络拓扑(Cauer / Foster)、电-热类比与「有效参与质量」的判法,承 A1-06 热阻网络法与集总参数建模入门;负荷分项的口径与取值来源,承 B1-01 整车热平衡与全工况热负荷谱分析 / C1-06 乘用车座舱制冷负荷计算书:逐项取值、UA 标定与环境舱对账 / C6-03 大空间/商用车/客车空调系统设计 / A2-04 湿空气物理与焓湿图:空调负荷计算基石 / C3-01 PTC 采暖系统设计与功率匹配 / B4-02 冬季续航掉电与热管理能耗拆解 / A4-07 多感官舒适:湿度、气味、噪声与光的协同。
出口:这一讲定下的状态向量,降阶之后就是车载在线估计器的模型骨架。⚠ 但降阶形式与实时性核算见 J7-02 模型降阶(ROM)与控制导向建模,输入可得性裁剪见 C5-03 低温起雾机理与主动防雾控制,观测器整定与退出条件见 K3-10 热管理在线状态与负荷估计:观测器整定、软测量置信度与退出条件。本课交模型与辨识,⛔ 不交控制器实现。
典型参考值(本讲用到的两个平台相关量)。⚠ 以下是典型参考带,用途是判断自己的模型量级对不对,⛔ 不是验收依据,须以本平台实测确认。
| 量 | 典型带 | 常见值锚点 | 怎么拿到它 |
|---|---|---|---|
| 舱容 V_cab(M1 类五座车) | 2.5 ~ 4.0 m³ | ≈ 3.0 m³ | 整车几何模型扣除座椅与内饰件占体积;⚠ 必须与实测 ACH_inf 时用的舱容同源(C1-06 乘用车座舱制冷负荷计算书:逐项取值、UA 标定与环境舱对账 4.3) |
| 节点数(含空气节点) | 5 ~ 9 个 | 起步 6 个(空气、玻璃、顶棚、仪表台、座椅、地毯),玻璃按受晒 / 背阴拆开后 7 个 | 由 1.3 三条判据逐条判出,再用 1.4 的上界收敛 |
本讲小结:状态向量必须先定后建模,交付形态是一张四列的节点表(节点 / 热容来源 / 边界 / 观测量),第三栏为空的行不可辨识、只有补测点或合并两条路。A1-06 热阻网络法与集总参数建模入门 的 Bi < 0.1 回答的是「沿厚度要不要分层」,而座舱要回答的是「两块能不能合并」——判据在边界不同,两块 Bi 都极小的玻璃只要一受晒一背阴就必须拆开,因为防雾要最冷那块、热浸要最热那块,两个下游要的都不是平均值。座舱另立三条判据:承担独立输出、时间常数差一个数量级以上、接受独立边界,任一成立即保留。全讲的地基是那条对偶关系:节点数的上界由独立测点数与激励段数决定,越线之后多出来的参数进入零空间,在标定工况上看不出来、换工况才爆发;处置按「先减节点、再加激励、最后加测点」的顺序,⛔ 而人为固定参数不在处置之列。前后排分区是「出口需求」与「有无测点」两条的合取,只有前者时按单空气节点加后排偏置处理。
后面还有 5 讲正文 · 关键公式 · 案例拆解 · 常见误区 · 动手做