1. 引言
新能源汽车电池包的隔热散热关系到电池安全及经历长期工况后的一致性问题,影响系统安全性和可靠性。电池热管理机组作为控制电池包温度的关键部件之一,机组内部布置了大量的管路和其他零部件,这些零部件通过螺栓或焊接的方式与机组钣金框架相连。整个机组受到车辆不同运行工况下的激励载荷,如果载荷频率与电池热管理机组某一阶固有频率接近或相同,可能会引起共振,导致局部应力过大,零部件及系统可靠性大大降低。GB/T 21361-2008对汽车用空调器振动试验的要求作了明确规定。
文中首先在前处理软件中建立整个电池热管理机组的有限元模型,利用 Abaqus 求解器进行系统的模态和谐响应分析;根据分析结果判定设计方案的可靠性和合理性。
2. 电池热管理机组的有限元模型建立
2.1 几何模型建立
几 何 模 型 是 建 立 CAE 有 限 元 模 型 的 基 础 , 在SolidWorks 中建立如图 1 所示的几何模型,为便于有限元前处理,将模型导出为.stp 格式。
图 1. 某新能源汽车电池热管理机组
图 1 中,电池热管理机组各个主要部件为:①安装脚;②钣金框架;③内部管路;④板式换热器;⑤压缩机;⑥冷凝器;⑦PTC 加热器;⑧风机。
2.2 有限元模型建立
针对电池热管理机组,采用专用前处理软件进行网格划分。在进行网格划分时,根据部件的几何特征,钣金结构①②、管路③使用壳单元,非考核部件④⑤⑦⑧使用壳单元+附加质量模拟,使其质量与重心位置与实际相同。单元类型为 S3 和 S4R,单元基本尺寸设为 10mm;
冷凝器⑥用六面体单元,单元类型为 C3D8R,单元基本尺寸设为 8mm;螺栓用 MPC 替代,由于模态分析及谐响应分析是线性分析技术,非线性的接触设置不起作用,因此部件间接触采用 Tie 模拟。最终有限元模型共有节点数62858,单元数 65696,电池热管理机组有限元模型如图2 所示,局部细节图如图 3 所示。
图 2. 电池热管理机组有限元模型
图 3. 电池热管理机组局部有限元模型
3. 电池热管理机组的模态分析
模态分析用于确定结构的固有频率和振型,固有频率和振型是承受动态载荷结构设计中的重要参数。模态分析也是响应谱分析和模态叠加法谐响应分析必需的前期分析过程。模态分析是研究结构动力特性的一种方法,是系统辨别方法在工程振动领域中的应用,是结构动态设计及设备故障诊断的重要方法[1-2]。本次分析采用Abaqus 求解器及分块 Lanczos 方法进行求解。
3.1 材料设定
在电池热管理机组中,参照图 1,其中①-⑤、⑦的材料为碳钢 Q235,⑥的材料为铝合金 3003,⑧的材料为ABS 塑料,计算中用到的材料属性如表 1 所示。
表 1 电池热管理机组材料属性
3.2 约束边界条件
电池热管理机组两侧有安装脚与框架相连,每侧安装脚各有两个螺栓孔。如图 4 红色虚线框所示,约束螺栓孔处节点的 1~6 自由度。
图 4. 热管理机组的安装脚约束示意图
3.3 模态分析
GBT 21361-2008 汽车用空调器振动试验方法中规定振动试验频率为 33Hz 或 67Hz。提取的模态频率范围应将其包含在内,设定提取模态数 0~20 阶,提交 Abaqus求解器进行求解,得到模态频率值如表 2 所示,若提取模态不足可重新设定更高阶数。图 5 为电池热管理机组前四阶模态振型图。
表 2 电池热管理机组的固有频率
(a)第一阶振型云图(振幅比/无量纲)
(b)第二阶振型云图(振幅比/无量纲)
(c)第三阶振型云图(振幅比/无量纲)
(d)第四阶振型云图(振幅比/无量纲)
图 5 电池热管理机组前四阶振型云图
电池热管理机组通过安装脚与振动试验台相连,从模态分析结果来看,第 13 阶模态频率和 GBT 21361-2008汽车用空调器振动试验方法中规定的试验频率67Hz一致,有可能发生共振现象。
4. 电池热管理机组的谐响应分析
谐响应是结构在周期载荷作用下产生的周期响应。谐响应分析能够预测结构的持续动力特性,验证结构设计是否能克服共振、疲劳及其他受迫振动引起的有害效果。谐响应分析计算结构在不同频率下的响应,可以观察到峰值响应频率对应的应力。谐响应分析是线性分析技术,任何非线性特性将被忽略[3]。本次分析采用 Abaqus求解器及模态叠加法进行求解。
4.1 载荷边界条件
按 GBT 21361-2008 汽车用空调器振动试验方法的规定,表 3 给出了 33Hz 和 67Hz 时电池热管理机组的载荷及加载方向。
表 3 电池热管理机组载荷及加载方向
4.2 谐响应分析
电池热管理机组在受到加速度简谐载荷时,因框架与压缩机、换热器等零部件的固定板均是钣金结构,这些钣金件受破坏风险最大,影响安装在上面部件的安全运行。因此,文章谐响应分析时,主要考察电池热管理机组的钣金件[4] 。
图 6~8 为电池热管理机组在激励频率 33Hz 时三种工况下的应力云图。从计算分析结果来看,横向 X 载荷工况下,其最大应力为 110.5Mpa,最大应力出现在背面钣金与底壳钣金的中间连接螺栓孔附近;纵向 Y 载荷工况下,其最大应力为 59.9Mpa,最大应力出现在板式换热器固定支架折弯附近;垂向 Z 载荷工况下,其最大应力为 86.2Mpa,最大应力出现在板式换热器固定支架折弯附近。
图 6 33 Hz 横向 X 工况应力云图(应力/MPa)
图 7 33 Hz 垂向 Y 工况应力云图(应力/MPa)
图 8 33 Hz 纵向 Z 工况应力云图(应力/MPa)
图 9 67 Hz 横向 X 工况应力云图(应力/MPa)
图 10 67 Hz 垂向 Y 工况应力云图(应力/MPa)
图 11 67 Hz 纵向 Z 工况应力云图(应力/MPa)
从计算分析结果来看,横向 X 载荷工况下,其最大应力为 53.7Mpa,最大应力出现在正面钣金与进出水口相连位置;纵向 Y 载荷工况下,其最大应力为 75.5Mpa,最大应力出现在板式换热器固定支架折弯附近;垂向 Z载荷工况下,其最大应力为 111.9Mpa,最大应力出现在出口管和高压加注口相连位置。
两种激励频率、三个载荷方向工况下,电池热管理机组的最大应力均未超过其材料的屈服强度 235Mpa。
5. 结论
文 章 以 某 电 池 热 管 理 机 组 为 研 究 对 象 , 运 用SolidWorks 建立了机组几何模型,利用专用前处理软件建立了有限元模型,在此基础上,用 Abaqus 求解器进行了模态分析和谐响应分析。分析表明,各工况下机组最大应力均未超过材料屈服强度,整体满足设计要求。
文章对某新能源轻卡电池热管理机组进行模态分析后,发现机组存在与试验频率 67Hz 相同的模态频率,后期还需进行优化设计,避开共振区域
[ABAQUS] Abaqus多物理场仿真如何解决复杂工程耦合问题
2025-06-17
[ABAQUS] 揭秘!Abaqus计算时间过长的5大原因及解决方案
2025-06-17
[ABAQUS] 别踩坑!Abaqus材料非线性设置的十大禁忌 - Abaq
2025-06-17
2025-06-16
2025-06-16
[行业资讯] Celanese与SIMULIA合作应对复杂的客户请求
2025-06-16
[ABAQUS] 为什么Abaqus网格划分失败显示边界三角形很差
2025-06-13
[CST] CST studio suite里的各种激励端口怎么用?
2025-06-13
2025-06-13
[行业资讯] 正版Abaqus软件价格及有限元分析软件Abaqus??榻?/a>
2025-06-12
2023-08-29
2023-08-24
[ABAQUS] ABAQUS如何建模?ABAQUS有限元分析教程
2023-07-07
[ABAQUS] 有限元分析软件abaqus单位在哪设置?【操作教程】
2023-09-05
[ABAQUS] ABAQUS软件中Interaction功能??橹械陌蠖ń?/a>
2023-07-19
[ABAQUS] Abaqus怎么撤回上一步操作?Abauqs教程
2024-05-01
[ABAQUS] abaqus里面s11、s12和u1、u2是什么意思?s和
2023-08-30
[ABAQUS] Abaqus单位对应关系及参数介绍-Abaqus软件
2023-11-20
[ABAQUS] ABAQUS软件教程|场变量输出历史变量输出
2023-07-18
[ABAQUS] Abaqus中的S、U、V、E、CF分别是什么意思?
2024-05-11
2025-06-16
2025-06-16
[行业资讯] Celanese与SIMULIA合作应对复杂的客户请求
2025-06-16
2025-06-12
[有限元知识] 如何通过达索代理商购买Abaqus有限元分析软件
2025-06-12
[有限元知识] 收藏!一文搞懂有限元分析的底层逻辑、操作要点和学习路径
2025-06-11
[有限元知识] 有限元分析时,网格尺寸通常如何确定?一般取多大?
2025-06-11
[行业资讯] SIMULIA电磁仿真线上技术汇|专题二:面向汽车电缆设计
2025-06-10
2025-06-10