1. 引言
随着汽车工业的不断发展,制动系统的性能和可靠性要求日益提高。电子机械制动器作为一种新型的制动技术,EMB电驱动制动系统,作为新能源汽车领域的一项革命性技术,正逐步成为未来智能汽车的标配。相较于传统的液压制动,EMB以其独特的优势引领着制动技术的新一轮变革,其装配工艺和测试技术对于保障系统性能和安全性至关重要。具有响应速度快、控制精度高、易于集成等优点,汽车电子机械制动器是目前国内研究的热点[1]。汽车底盘本土龙头厂商伯特利在2023年报中表示已开展对EMB的研发工作,EMB的A轮首样已于2023年8月完成制作,并已进行冬季试验验证。对电子机械制动器进行准确地建模仿真,能够更好地模拟汽车制动时,EMB的制动情况,为产品的研发和优化提供技术性支持。本文对EMB进行设计建模、仿真,模拟真实汽车制动情形[2]。依据目标车型参数,对EMB各部件如电机、减速增矩装置、运动转换装置进行选型设计;利用CATIA绘制出其三维结构图,利用ANSYS/Workbench对其三维模型中的刹车片和刹车盘进行瞬态动力学的有限元分析。
2. 电子机械制动器的设计
2.1. EMB组成及工作原理
Figure 1. EMB working principle diagram
图1. EMB工作原理图
电子机械制动器主要由电机、行星齿轮减速器、滚珠丝杠、制动钳和制动盘等部件组成[3]。其工作原理是给电机供电,电机转动产生的动力经行星齿轮减速器进行减速增扭,然后通过滚珠丝杠将旋转运动转换为活塞的直线运动,从而推动制动钳对制动盘施加制动力,实现汽车减速。
EMB工作原理图如图1所示。
2.2. 制动效能校核
乘用车的液压盘式制动器制动管路最大油压为12 MPa,制动器活塞直径为50 mm [4] EMB最大制动力的计算如下式:
(1)
式中:Fmax——最大夹紧力;Pmax——制动最大油压;d——活塞直径。
汽车制动过程是一个高温、高扭矩、高压下的滑摩过程,忽略惯性力偶和阻力力偶,参考车辆车轮在良好的路面上的受力情况如图2所示。
Figure 2. Wheel stress situation
图2. 车轮受力情况
由力矩平衡原理可知:
(2)
式中:FXb——地面制动力;Tu——摩擦力矩;r——轮胎半径。
地面最大制动力是地面对车轮的法向作用力和摩擦系数的乘积,因此地面制动力要小于等于地面约束条件,即:
(3)
式中:Fφ——地面最大制动力;Fz——地面对车轮的法向作用力;φ——附着系数。
单个制动器的制动力:
(4)
2.3. 力矩电机选型
电子机械制动器的力矩电机选型要求力矩电机的转速要超过380 r/s,初选减速器的减速比为16,初选EMB的电机型号为80BF0140U30-24VD的无刷直流力矩电机。
2.4. 滚珠丝杠选型计算
滚珠丝杠的选型由平均载荷和平均转速决定,汽车行驶过程中一般不会出现猛踩刹车踏板的情况,整体处于中上等制动[5],因此取2/3的平均载荷即可。参考车辆的刹车盘和刹车片的间隙为2 mm,消除间隙的时间0.1 s。取平均转速的1/3用于丝杠额定动载荷的校核。校核公式如下:
(5)
2.5. 行星齿轮减速机构的设计
行星齿轮减速机构具有传动比大,效率高等优点,因此电子机械制动器采用行星齿轮减速机构[6] [7]。一般行星齿轮传动机构中,行星齿轮数目一般为3个,安装时应当满足以下条件:
a) 相邻两个行星齿轮的齿顶圆半径之和要小于它们的中心距。
b) 考虑到安装的需要,太阳轮齿数、内齿圈齿数、行星齿轮的齿数通常要满足太阳轮和内齿圈齿数之和是行星齿轮的整数倍的关系[8]。
根据上述条件取太阳轮齿数Z1 = 12,内齿圈齿数Z2 = 66,传动比i = 3则:
行星齿轮齿数:
(6)
得到行星齿轮Z3 = 26。
3. Catia三维建模
基于上述选型计算结果,使用Catia软件创建电子机械制动器各部件的三维模型,并进行精确装配,确保模型的准确性和完整性,为其相关零部件的有限元分析做铺垫。三维图形如图3所示。
Figure 3. Three dimensional model of electronic mechanical brake
图3. 电子机械制动器的三维模型
4. ANSYS瞬态动力学仿真
EMB产品研发过程中,利用ANSYS进行仿真,能够对其研发提供技术支持,能够规避其在实际使用工程中出现的问题[9],据此可进行改进和优化,可大大缩短研发周期,提高产品的可靠性。由于在刹车过程中只有刹车盘和刹车片进行接触摩擦。
1) 模型导入与处理
将Catia建立的刹车盘和刹车片模型导入ANSYS/Workbench,进行必要的简化和清理。
2) 材料属性定义
盘式制动器刹车片的材料为陶瓷,因此设置刹车片材料属性时,密度为3900 kg/m3,杨式模量为300 GPa,泊松比为0.23。刹车盘的材料一般为灰铸铁,密度为7300 kg/m3,杨氏模量为130 GPa,泊松比为0.25。
3) 网格划分
制动时,滚珠丝杠将行星齿轮减速机构传来的旋转运动,转换为直线运动,并且推动活塞前进,活塞推动刹车片夹紧刹车盘,依靠刹车片和刹车盘的摩擦作用实现汽车的减速因此,为了更加精确地模拟EMB的受力情况,在Geometry中对刹车盘、刹车片进行分切割及共节点操作,并且为了提高其计算精度,采用合适的网格类型和尺寸对模型进行划分如图4所示,以平衡计算精度和效率。
Figure 4. Grid division of brake pads and brake discs
图4. 刹车片和刹车盘的网格划分
网格划分后总共形成26,428个单元,42,567个节点,设置接触条件,将刹车盘和刹车片的接触设为摩擦的,并且设置摩擦系数为0.2.由于制动过程中,刹车盘随着车轮一起转动,因此还要设置旋转副,并且设置旋转速度为3.14 r/s,即每秒转半圈。
4) 边界条件与载荷施加
对模型施加约束,对刹车片施加分量类型的正压力,X方向和Z方向施加0 MPa的压力,Y方向施加−2 MPa的正压力。添加位移约束,并且设置在X和Z方向的分量上不允许模型移动,在Y方向上设置自由的分量,来观察刹车盘在正压力下的变形情况。
5) 求解设置
添加总变形、等效应力、接触工具来观察仿真结果,分析设置里面修改初始时部为0.05 s,最小时部为0.05 s,最大时部为0.1 s,启动瞬态动力学求解。
5. 仿真结果与分析
经过求解得到力收敛图如图5所示。
Figure 5. Force convergence diagram
图5. 力收敛图
由图5结果可知,随着刹车盘的转动,力收敛曲线呈现周期性收敛。
制动力响应曲线如图6所示。
Figure 6. Braking force response curve
图6. 制动力响应曲线
获得制动力随时间的变化曲线,由曲线可知制动力随着时间的增加在1.81 MPa处呈现周期性上下波动。
总变形、压力分布云图如图7、图8所示。
由总变形和压力分布云图可知[10],刹车盘的最大变形处发生在与刹车片接触的最外侧边缘处,最大应力发生在刹车片的下方左右两侧边缘,最大应力为118.39 MPa。陶瓷的主要成分是Al2O3,其抗压强度为850 Mpa。灰铸铁的抗压强度则在600~1000 MPa。取安全系数为2,因此根据第三强度校核,EMB制动时的最大应力仅为安全极限的27.85%,满足设计要求。
Figure 7. Cloud map of total deformation distribution
图7. 总变形分布云图
Figure 8. Cloud map of pressure distribution
图8. 压力分布云图
6. 结论
通过详细的选型计算、Catia建模和ANSYS/Workbench瞬态动力学仿真,对电子机械制动器的性能进行了全面评估。研究结论如下:
1) 最大应力为118.39 MPa,EMB制动时的最大应力仅为安全极限的27.85%,满足制动强度的要求。
2) 根据地面最大制动力的约束条件所选的电机、滚珠丝杠和行星齿轮减速器能够满足制动需求,各部件的应力和变形在合理范围内。
3) 仿真结果与实际结果吻合度在90%以上,可以有效提高车辆行驶的安全性,为EMB的研发提供理论基础。