2019, Vol. 40
2. 农业部南京农业机械化研究所, 江苏 南京 210014
2. Nanjing Research Institute for Agricultural Mechanization Ministry of Agriculture, Nanjing 210014, China
电动汽车的快速发展使得线控制动技术受到越来越多研究人员的重视.随着人们对汽车制动系统要求逐渐提高[1-3], 传统液压制动系统结构变得越发复杂, 而由于其自身的结构特点, 使得汽车在制动时始终存在制动延迟、脉动的特点.近些年研究人员开发的轮毂制动技术使得汽车的制动效果显著改善[4-5], 但其复杂的工艺结构以及较高的密封要求无疑会增加汽车制造成本, 且目前轮毂制动所能产生的制动力矩仍不能满足汽车所需.磁流变制动器结构简单紧凑, 有着制动连续可控、响应迅速的特点[6].本文主要基于结构较轻便的盘式磁流变制动器进行反应性能的相关研究[7-8].
1 建立制动力矩数学模型盘式磁流变制动器的工作原理如图 1所示.不需制动时, 励磁线圈不施加电流, 磁流变液呈牛顿流体状态, 制动盘此时仅受可以忽略的黏性力矩; 当汽车需要进行制动时, 励磁线圈中通入电流, 磁流变液中的磁性颗粒所形成的链状结构会与接触面产生一定的摩擦力, 阻碍接触面运动, 从而对制动盘产生剪切制动力矩.
|
图 1 线圈外置型磁流变制动器原理图 Fig.1 Schematic diagram of MRB with external coil |
研究表明, 磁流变液的理想特性可由Bingham塑性模型来加以描述[9].根据该模型, 磁流变液的剪切应力表示如下:
|
(1) |
式中:τ为磁流变液的剪切应力; τy为磁流变液的动态屈服应力; η为磁流变液的黏度; 为磁流变液的剪切应变率.
磁流变液的动态屈服应力与磁场强度的函数近似如下:
|
(2) |
式中, k与β为与磁流变液本身有关的定常数.
假设磁流变液是非流动状态, 则其剪切应变率为
|
(3) |
式中:ω为制动盘的转动角速度; δ为制动器工作间隙的大小.
在制动盘上, 取工作在半径r处宽度为dr圆环区域, 其产生的转矩为
|
(4) |
则制动力矩T1为
|
(5) |
由于制动盘制动时两端面产生制动力矩, 因此总的制动力矩T为
|
(6) |
根据现有加工及实验条件, 对制动器各参数进行设定.首先, 工作间隙δ是影响制动力矩大小的重要因素[10], 工作间隙越大, 制动器所产生的制动力矩越小, 而较小的工作间隙虽然能产生较大的制动力矩, 却对加工精度的要求更高, 也增大了装配的难度, 一般其取值范围为0.25~2 mm, 本文选取δ值为0.5 mm.
由于实验室实验台的环境条件等对制动器的径向大小有一定的限制, 因此需要根据实际情况确定制动器的径向尺寸, 设定制动器最外径R4=70 mm, 线圈宽度b=20 mm, 制动盘的厚度h=6 mm, 制动盘内径R1=18 mm, 设定线圈外径R3、制动盘外径R2以及壳体厚度l为优化变量.
理论上, 磁流变制动器的传动轴除需要足够的强度外, 还要有较低的导磁率, 从而使线圈产生的磁场更能集中穿过工作间隙的磁流变液, 1Cr18Ni9Ti合金钢是一种高强度的弱磁性材料, 因此选其作为实验环境的传动轴; 左右外壳、制动盘及外环均需具有较好导磁性能, 因此选用Q235碳钢, 制动器的励磁线圈不仅要满足低磁导率要求, 还应具有良好导热性, 铜的导热率为384 W/(m·℃), 因此采用漆包线最大外径为0.96 mm、铜芯的直径为0.9 mm型号规格铜漆包线.
2.2 制动器结构仿真优化利用ANSYS有限元分析软件进行制动器结构优化.仿真中, 定义相关材料单元类型PLANE13;边界条件定义为平行; 施加的载荷类型为电流密度.根据制动器结构及前文设置的参数值, 设计变量的设置范围如下:
|
(7) |
优化设计一般是使目标函数最小化, 且目标函数值需为正值, 而本文优化目标为max T, 即实现制动力矩最大化, 因此将目标函数变为-max T, 同时需使目标函数值加上一个较大的正数以防止目标函数值最终出现负值.优化方法选择一阶法.最终确定l=9 mm, R2=46 mm, R3=60 mm.整理优化结果得出制动器结构参数如表 1所示.
|
表 1 制动器结构参数 Table 1 Brake structural parameters |
制动器的示意图如图 2所示.
|
图 2 制动器示意图 Fig.2 Schematic photo of the brake |
为方便制动器进行相应实验, 将制动器固定于支架上.
3 磁流变制动控制器设计 3.1 控制器硬件设计本文使用如图 3所示的Arduino开发板作为控制器的核心部件.
|
图 3 Arduino实物图 Fig.3 Physical map of Arduino |
因磁流变制动器所需电流一般最大为3 A, 因此需要对电流进行放大, 本文选用L9349型号功率驱动芯片, 该芯片最大负载电流为5 A, 完全满足实验的需要.芯片的工作电压为4.5~32 V, 通过输入PWM控制信号, 就能方便地控制输出端以控制磁流变制动器进行制动.
在电路中, 利用Arduino开发板的两个模拟信号端口A0与A1分别接收来自信号发生器模拟的速度信号和转矩转速传感器产生的转速信号; 利用数字信号输出端口9作为PWM信号的输出端口, 并将其与L9349的信号输入端口2相连; 将数字信号端口10与L9349芯片使能端相连, 通过控制端口10输出高低电平来控制功率驱动芯片工作与否; 将功率驱动芯片的工作电压引脚VS通过导线与开关电源的正极相连, 将负极引脚PGND与开关电源的负极相连; 将芯片的输出引脚2与磁流变制动器的线圈一端连接, 并将制动器另一端的线圈与开关电源的正极相连; 将开发板的信号地引脚与芯片的信号地引脚相连.
3.2 控制器软件设计在汽车ABS控制系统中, PID控制是目前实际应用最为广泛的一种自动控制方式, 是按偏差的比例、积分和微分进行控制的, 其原理见图 4[11].
|
图 4 PID算法控制原理 Fig.4 PID algorithm control principle |
Arduino内部的PID库函数可免去PID控制编程, 直接调用即可.
PID的连续时间控制律可以表示为
|
(8) |
式中:u(t)为控制器的输出; e(t)为期望值与实际输入值的差值; KP, KI, KD分别为控制器的比例、积分、微分系数, 这3个系数的取值直接决定着控制系统的性能, 因此必须正确地设定3个系数使得控制系统的响应特性更迅速, 本文使用PID参数设定方法为Ziegler-Nichols方法, 其主要调整方式如下.
首先定义一个临界的P增益(KC).确定KC大小的方法如下:
1) 设置KI, KD为0;
2) KP值增加直到系统以一个稳定的速率振荡, 即为KC, 一旦KC被定义, 那么必须测量振荡的周期.根据周期PC, 可以用表 2计算KP, KI, KD参数的数值.
|
|
表 2 PID控制参数表 Table 2 PID control parameters |
首先需要找到临界的P增益, 设置KI与KD为0, 只利用比例调节.设置中逐渐增加参数P的值, 记录系统的响应.图 5是4种典型的系统特点, 图 5a表示P值设定过低, 需要继续调高; 图 5b表示P值仍然过低; 图 5c表示P值仍然过高, 导致系统不稳定; 图 5d是合适的P值调节出的效果, 这时系统出现稳定的振荡.这时P值即为KC的值, 系统振荡的周期即为参数PC, 则根据表 2得出系统的KP, KI, KD的数值.
|
图 5 PID参数调节方法 Fig.5 PID parameter adjustment method (a)—临界P增益过低；(b)—临界P增益较低；(c)—临界P增益过高；(d)—临界P增益适合. |
磁流变制动器线性可控的特点是其重要的优势之一, 制动器的制动力矩与线圈中通入的电流的线性关系可极大地方便制动力矩的控制.
绘制制动力矩控制实验的流程图如图 6所示.
|
图 6 制动力矩控制实验流程图 Fig.6 Experimental flow chart of braking torque controlling |
磁流变制动器之所以有替代传统液压制动系统的潜力, 不仅在于其仅仅依靠电流的改变即可实现制动力矩的改变即单一的控制参数, 还在于其快速的响应速度.但磁流变制动器制动时也存在制动动作的响应延迟以及时间的滞后, 为探索磁流变制动器的响应时间及滞后时间, 本文将制动力矩期望值首先设置为阶跃信号规律.
实验过程如下:由信号发生器产生阶跃规律的制动力矩信号-制动力矩期望值TE, 并将其通入Arduino开发板的模拟信号输入端口A1;同时, 将由转速转矩传感器测得的信号-制动力矩测量值TC通入Arduino开发板模拟信号输入端口A2;Arduino开发板内部将TE，TC模拟信号转换为数字信号, 并对二者进行比较, 其差值为误差eS, 作为PID控制的输入, 经过PID控制后输出具有适当的占空比N的PWM信号, 将PWM信号输入功率驱动芯片的信号输入引脚17, 进而控制磁流变制动器, 使其产生所需的制动力矩; 由转速转矩传感器测量制动器产生的制动力矩, 并将其输入Arduino开发板的模拟信号输入端口A2, 至此完成一个控制周期.
为进行制动力矩的相关实验, 根据制动力矩流程图搭建实验平台, 如图 7所示.
|
图 7 制动力矩控制实验 Fig.7 Braking torque control test |
实验台中, 选择型号为SK11E-90S/4的减速电机, 该电机的各个参数如表 3所示.
|
|
表 3 减速电机参数 Table 3 Gear motor parameters |
惯性飞轮作为一种吸收和释放能量的机构, 在实验中主要起到缓冲的作用, 防止制动器突然制动对电机造成伤害; 信号发生器负责产生具有阶跃规律的制动力矩期望值; 笔记本电脑作为上位机, 可利用Arduino特有的编程器直接通过USB接口将程序代码下载到Arduino内部的微控制器中.上位机除利用USB数据线对Arduino进行程序传输外, 还对开发板进行供电, 免去了供电电源和烧录器的使用; 此外Arduino开发板能在使用过程中对PID控制参数进行及时调整.
实验中, 首先给减速电机供电, 使其正常转动, 便于进行制动力矩的产生; 其次给开关电源供电使其对功率驱动芯片电源进行供电并将已写好的程序下载到开发板中, 进行制动力矩控制实验; 同时打开信号发生器, 使其产生预先设定的制动力矩的期望值; 最后将实验中产生的数据存储, 以便后续的处理及使用.
阶跃制动力矩实验结果如图 8所示.实验中, 设置了3种不同的阶跃制动力矩值, 制动力矩分别为4, 6, 8 N·m, 从前部分实验结果得知, 制动器的响应时间为40 ms左右; 设定制动力矩首次达到最大值时对应的时间与响应时间的差值为制动器的滞后时间, 由实验数据显示得制动器的滞后时间约为70 ms.总体响应加滞后时间约110 ms, 与液压制动系统(200~800 ms)相比, 极大地缩短了响应时间.
|
图 8 阶跃制动力矩测量值 Fig.8 Measurement of step braking torque |
为体现磁流变制动器可实现连续可变的制动力矩控制效果, 设定制动力矩期望值为正弦变化规律.Arduino开发板理论最大采样周期为10 kHz, 设置制动力矩期望值的信号周期为200 ms, 即制动力矩的期望值的频率为5 Hz, 将期望值通入开发板的端口进行制动力矩跟随实验, 实验结果如图 9所示.
|
图 9 制动力矩跟随实验 Fig.9 Test of braking torque following |
从正弦制动力矩实验图可知, 控制器控制磁流变制动器产生的制动力矩相比制动力矩期望值延后20 ms左右.以上实验结果表明制动器响应迅速, 能够产生较好的跟随效果.
4.4 制动减速度实验分析制动减速度也是衡量制动器制动性能的关键因素, 它的优劣直接影响到制动距离的长短, 因此, 对制动器做减速度实验.车载电流一般最大为3 A, 实验时用控制器给制动器施加3 A电流, 产生的制动力矩稳定在14 N·m, 在电机启动之后断开电磁离合器开始制动减速度过程, 得到转速与时间关系如图 10所示.从图 10中可以看出转速近似以恒定减速度减小, 制动时间在1 s之内, 说明制动器制动迅速, 性能稳定.
|
图 10 转速与时间关系 Fig.10 Relationship between rotational speed and time |
本文主要设计了一种盘式磁流变制动器, 利用ANSYS对制动器的结构进行仿真, 并进一步对制动器的结构进行优化, 对比优化前后的制动力矩, 体现了优化后的制动器的结构更加合理; 加工磁流变制动器样机; 利用Arduino开发板以及L9349芯片设计磁流变制动控制器, 并进行了制动力矩控制相关实验, 结果表明相比液压制动系统而言, 其优势明显.此研究对磁流变制动器后续的实际应用奠定了基础.
| [1] |
Clarke P, Muneer T, Cullinane K.
Cutting vehicle emissions with regenerative braking[J]. Transportation Research Part D, 2010, 15(3): 160–167.
DOI:10.1016/j.trd.2009.11.002 |
| [2] |
Gerla M, Lee E K, Pau G, et al.
Internet of vehicles:from intelligent grid to autonomous cars and vehicular clouds[J]. Internet of Things, IEEE, 2016, 16(1): 241–246.
|
| [3] |
Finn A, Scheding S.
Developments and challenges for autonomous unmanned vehicles[M]. Berlin: Intelligent Systems Reference Library, 2012: 56-77.
|
| [4] |
贺辉, 宋大凤, 杨南南, 等.
轮毂马达液驱系统控制与仿真[J]. 吉林大学学报(工学版), 2012, 42(1): 27–31.
( He Hui, Song Da-feng, Yang Nan-nan, et al. Control and simulation of hydraulic in-wheel motor propulsion system[J]. Journal of Jilin University(Engineering and Technology Edition), 2012, 42(1): 27–31. ) |
| [5] |
陈勇, 张建荣, 张大明.
电动轮技术在电动汽车中的应用及发展趋势[J]. 机械设计与制造, 2006(10): 169–171.
( Chen Yong, Zhang Jian-rong, Zhang Da-ming. State-of-art of application of hub-motor in electric vehicle[J]. Machinery Design and Manufacture, 2006(10): 169–171. DOI:10.3969/j.issn.1001-3997.2006.10.070 ) |
| [6] |
Gudmundsson K H, Jonsdottir F, Thorsteinsson F.
A geometrical optimization of a magneto-rheological rotary brake in a prosthetic knee[J]. Smart Materials and Structures, 2010, 19(19): 035023.
|
| [7] |
Park E J, Luz L F D, Suleman A.
Multidisciplinary design optimization of an automotive magnetorheological brake design[J]. Computers and Structures, 2008, 86(35): 207–216.
|
| [8] |
Farjoud A, Vahdati N, Fah Y F.
Mathematical model of drum-type MR brakes using herschel-bulkley shear model[J]. Journal of Intelligent Material Systems and Structures, 2008, 19(5): 565–572.
DOI:10.1177/1045389X07077851 |
| [9] |
Ashour O.
Magnetorheological fluids:materials, characterization, and devices[J]. Journal of Intelligent Material Systems and Structures, 2010, 7(2): 123–130.
|
| [10] |
Song W L, Li D H, Tao Y, et al.
Simulation and experimentation of a magnetorheological brake with adjustable gap[J]. Journal of Intelligent Material Systems and Structures, 2016, 28(12): 1614–1626.
|
| [11] |
宋宇.磁流变液ABS制动器结构研究与性能仿真[D].哈尔滨: 东北林业大学, 2011.
( Song Yu.Structure research and performance simulation of magnetorheological fluid ABS brake[D].Harbin: Northeast Forestry University, 2011. ) |