system: [SYSTEM]

user: 3. 系统计算机仿真和理论分析 3.1. 电动空调分析 电动空调主要作用是给液体加热或者制冷。根据电动大巴的工作状态分析，结合箱体各种材料的换热系数及电池本身的发热。 传热基本方程： Q = k A Δ t (1) 式(1)中：k：热阻，对稳态热传导：λ/δ (导热系数除以厚度)暂不考虑热对流、热辐射的传热；A：传热面积；Δt：温差；电池箱中使用了聚氨酯保温材料，其导热系数：0.025 W/(m*K)，燃烧等级：阻燃B2，保温材料厚度δ = 0.008m；其结构使用玻璃钢板材质，其导热系数约为0.23 W/m*K，玻璃钢壳体厚度δ = 0.005 m；电池箱散热总面积S = 25 m²。在本地如果电池箱系统与环境温差达到10摄氏度左右，基本上满足是设计要求。并结合计算电池本身发热量，可以得到设计功率。以上分析计算得，一辆大巴要使电池箱内温度在20度到30度之内的恒温状态，需要大概制冷功率在2800瓦左右，如果加上管道等器材的热量损失，设计功率在3000瓦。以上计算结果为静态的，得到水冷空调的最小功率。图2是设计结构图和仿真示意图。 图1. 系统示意图 图2. 设计结构图和仿真示意图 相关电动空调的仿真和计算，现阶段比较精确和方便的理论是稳态空调仿真。相关软件仿真非常多。这里使用了Matlab [ 1 ] 软件进行电动空调的稳态分析 [ 2 ] ，其精度满足应用需求。根据上面分析得到的功率等数据，按要求输入相关参数，就可以得出空调4大器材运行参数和配型要求，并结合水冷空调的特点，就能设计出比较合理的系统。这里需要考虑的是膨胀水箱的设计和水质软化稳定及结冰等情况。至此，本系统的核心设备设计基本设计完毕。 3.2. 管道及其控制系统分析 水管管道是控制水的流向和流速的，不仅如此，还具有保温功能。关于管道的选型要从安装是否方便，成本等因素考虑。处于设计阶段，管道参数的计算是比较麻烦的。基本上应用了流体动力学所有知识。 现阶段，系统只仿真9个电池箱体的管道参数。在设计之前，必须了解水的参数和纲量及一些弯管，水阀，膨胀阀，止回阀等特殊管材的等效长度(局部阻力)，分析水管沿程阻力和导热参数。以上参数均有相关表格给出数据。被研究的能量载体是水，属于不可压缩牛顿流体，仿真其运动是比较有规律的。整个管道系统中流体能量大致分为以温度为表现形式的水的内能，以流体运动为形式的水的动能和势能。这些能量符合守恒定律。水的动能来自于水泵和管壁沿程阻力及局部阻力，水的内能来自于电动空调，管道内壁散热，电池芯体热能等与之实现热交换的设备。掌握了流体运动和能量转化的规律，可以比较客观的仿真系统。建立分散式仿真模型，运用能量和动量守恒定律，输入系统已知参数(最小功率等)和管道的固有参数 [ 3 ] ，可以得到主水泵最小运行功率，管径，最大散热功率，每个分水管的水阀开度(流速)等数据。依据参数搭建实验实体。不考虑管路热交换和液体内能变化。下式(2)表示了能量守恒定律： W 势 能 + W 动 能 + W 内 能 + W 压 强 能 = 固 定 值 (2) 其中： W 势 能 = ρ g h ;   W 动 能 = 1 2 m υ 2 ；   W 压 强 能 = p V ；   V = ( π d 2 4 ) υ t ;   m = ρ V ; 这里p和v为某时间点的流体压强和流速，d是管径，h为管路取样位置， ρ 是流体密度。 W 内 能 这里只指示流体的分子平均动能(温度决定的)。联立公式建立p和v的偏微分方程，利用Matlab软件进行仿真。考虑管路的沿程损耗，引入管路系统的干扰因子，可以更精确的计算数值。 水循环恒温控制系统主要实现各个设备的开关机顺序及能量PID控制等。设计特点体现在使用了4路模拟温度探头和3路数字温度探头，可以将现场温度的突变和精准度充分的采集，更好控制系统的能量分配，节省运行成本。同时使用了大容量的存储设备进行关键数据保存和紧急意外处理机理，实现系统的增值服务，提高系统稳定性 [ 4 ] 。图3所示系统控制器简图。 图3. 系统控制器简图 图3中6路程序控制，按照控制策略进行控制，并且智能的处理设备之间的能量转化。控制策略的最优方法的选取，关键建立整个系统的数学模型。首先建立系统的受控状态方程： x ∘ = f ( x ， p ， v ， t ) ， (3) 其状态方程可以从式(2)中获取。寻求不受控制量u(与水冷空调输出功率成线性关系)，使系统从初始状态 x ( t 0 ) = x 0 ， (4) 在有限的时间间隔 [ t 0 ， t g ] 内，将系统状态转移到 x ( t g ) 下，并满足能量守恒定律，即 W ( x ( t g ) ， t g ) = 固 定 值 . (5) 这样系统建立泛函数模型，可以解出最优控制函数。控制系统中的CAN通讯与电池管理系统进行温度和充电数据的交互。同时232通讯与上位机或者与GPS进行通讯，传输有关系统的重要的温度和电池箱体内部照片等信息。图4为控制系统采集软件，通过与系统的232接口通讯，可以对7路温度进行采集分析和实时观测电池箱内的状态，在实体实验中，为系统的分析提供必要的数据支持。控制系统工作在电池箱内，工作环境比较恶劣，同时需要处理模拟数据，在设计电源时更加注意，图5所示电源设计思路。 图4. 控制系统采集软件 3.3. 电池箱体仿真和分析 电池箱内的温度控制是本文的重点又是难点。系统采用了电池箱体内相对封闭的结构，使用导热比较差的玻璃钢箱体，并在内部进行了隔热处理。将冷水流过的表冷器放入箱体内，使用涡流风扇，进行箱体内冷风循环，进而将电芯的热量带走。图6为电池箱仿真模型。 图5. 控制系统电源设计 图6. 电池箱仿真模型 上图使用了Matlab软件进行系统仿真 [ 2 ] 。这里介绍图中字母意义如下：空气中温度为X a ，电池的质量为m，电池的平均比热容为c，与空气的传热面积为A，电池与空气的当量传热系数为K，可以分析电池温度X随时间的关系可以建立微分方程。 K A ( X a − X ) = c m d X d t (6) 这里表冷器主要为能量输出设备，电池是被控制设备，空气等为媒介。V 1 为表冷器输出空气体积，假设表冷器输出空气为10摄氏度，V 2 为电池箱内的空气体积， ρ 为空气密度。并假设气体V 1 和气体V 2 有相互接触的边界，并且有一定的相互渗入，可以增加修正微分方程，建立电池箱空气温度X a 随时间变化的方程。 ∫ t t ξ c ρ V 1 ( X a − 10 ) d t = c ρ V 2 d X a d t + W 损 耗 (7) 联立式(6)和式(7)，建立状态函数 φ = φ ( X a ( t ) ， X ( t ) ， t ) ; (8) 进而进行变分法分析，约束相关变量，得电池温度X仿真图形。在Matlab软件里输入仿真参数， 假设电池箱体不蓄热，电池箱内空气体积为0.3 m 3 ，又假设箱内温度为45度，送风温度为10度，电池本身和风机等产生的热量为800瓦，空气热容为1000，单位时间内从表冷器输出的空气体积量为6 L/S*4 (4个涡流风扇)，由于电池箱存在热交换，则系统渗出热为100瓦。可以获得电芯发出不同热能时，随时间最终的箱体内温度平衡的曲线图。如图7所示的电池箱仿真曲线图。 图7. 电池箱仿真曲线图 这里以热负载为质点形式给出，没有考虑到热负载的散热面积等因素。并结合静态理论分析，可以估算实际情况。根据图7可以发现，当表冷器输出温度越低，平衡温度越低；输出风量越大，越接近表冷器翘片最低温度。这里仿真的原则是风道进行畅通循环状态下。如果进行双边强对流，中间自然扩散的设计原理。最终平衡的温度可以控制在30度以下。故在设计箱体系统时候，充分考虑以上理论原理。 3.4. 表冷器仿真和分析 电池箱体内的表冷器设备是系统能量交换中心点，是整个恒温系统的关键点。表冷器一方面要将冷水的能量充分释放，一方面通过风循环，将电芯的热量充分吸收。它是典型的气液热交换系统，但它又区别于其他的水冷系统。主要区别为芯体结构特殊，箱体空间比较小并且空间相对封闭，热量比较分散等。根据恒温系统要求，进行如图8所示的表冷器设计。 图8. 表冷器仿真 这里的表冷器为翘片冲缝多片并联换热器。保压测试大于1 Mpa；散热面积大于3平方米。仿真时 [ 2 ] ，忽略金属器壁的热容，物性参数不受温度影响。图示中字母说明如下，M a ：空气质量流量；M w ：冷水质量流量；T ai ：空气入口温度；T wi ：冷水入口温度；T ao ：空气出口温度；T wo ：冷水出口温度。 根据能量守恒得式(9)： c 水 M a ( T w o − T w i ) = c 气 M a ( T a o − T a i ) ; (9) 由于 M = ρ S υ t ; (10) 其中 ρ 物体密度，S流体窗口截面积，v为流体质量流速。理想情况下， T w i ≤ T a o ≤ T w o 。式(9)、式(10)结合式(2)、式(6)、式(7)。可以动态的分析表冷器状态。 如果输入表冷器固有参数(一般按照设计要求选取)，并假设风机吹出的风速各点一致，不考虑挡风板的因素，进入表冷器的各点风速流速一致，且冷水质量流速连续。可以仿真出冷水出口温度，空气出口温度和湿度。更改不同的表冷器固有参数，可以计算出符合我们要求的输出参数，对我们选型表冷器有重大的作用。经过仿真分析获得冷水质量流量值和空气质量流量值，均比较大的时候，输出制冷效果比较好。因为系统相对比较封闭(仿真的时候，假设空气含湿量为系统原则上最大值，即考虑了很少湿空气进入系统的情况)，其结果在出风口处含湿量比较少，冷凝水出现比较少，用在电池箱内，符合安全使用条件。 电芯温度主要靠箱体内风冷的原理进行热交换的，这样表冷器周围的风场分析是很重要的，然而风场属于可压缩牛顿流体，研究它比较困难，在工程中也只是估算风场状态的。为了降低分析难度，从而进行风静场研究。图9为电池箱内风场示意图。 图9. 电池箱内风场示意图 图9中，展示了风场中层流和湍流在电池箱体的独特特点。当冷风从表冷器中出来 [ 3 ] ，由于轴流风机风压高，通过电池箱内部较大风道，走到电池箱内部底面。其风速各点几乎相似，即形成围绕电池(除底面)的稳定低温场。低温区域，空气密度比较大，高温区域空气密度比较小，形成低温向高温自然对流扩散(温度差决定了扩散能力)。在风机入风口出，由于风阻比较大，导致在风机入口处形成了湍流的回风区域，其区域温度比较低。这是基于静态风场的分析。结合表冷器仿真结果，可以比较准确的仿真电池箱体内的温度分布和风场分布。