Effects of different ambient temperatures on performance of electric vehicles' heat pump air conditioning
-
摘要: 针对电动汽车冬天取暖能耗较高,设计了基于蒸气压缩循环的热泵空调试验系统,研究了制热模式下不同的环境温度和压缩机转速对车室内平均温度、高压管路内部工质的温度和压力、系统能效比(COP, Coefficient of Performance)等参数的影响.试验结果表明,当环境温度和压缩机的转速较高时,压缩机出口、车室内换热器出口处工质的温度和压力值较大,车室内平均温度的上升速度越快,达到舒适温度所需的时间越短.在压缩机转速相同时,环境温度越高,系统的COP值越大.Abstract: A heat pump air conditioning experimental system for electric vehicles based on the vapor-compression cycle was designed because much energy was consumed for the electric vehicle air conditioning under the heating mode when the ambient temperature was low. The refrigerant temperature and pressure in the high-pressure pipeline, the rate of temperature rising in the cab and the coefficient of performance (COP) for the air conditioning system were studied in the heating mode under different ambient temperatures and different compressor speeds. The test results indicate that, when the ambient temperature and the compressor speeds are higher, the temperature and pressure of the refrigerant at the compressor outlet and the heat exchanger outlet are higher, and the average temperature increases faster and the time reached to the comfort temperature is shorter in the cab. When the ambient temperature is higher in case that the compressor kept the same speed,the COP is greater.
-
绳索具有线密度低、强度高、柔性好等特点，因此广泛应用于各种工程领域，包括斜拉桥[1]、船舶起重[2]、绳系卫星[3]、大型射电望远镜[4]、风洞试验[5-7]等。绳索的阻尼特性是影响机构功能的重要因素之一，也是设计中需要考虑的重要参数。绳索的材料、长度、直径、预紧力不同，都会导致绳索阻尼有较大差异，而目前尚缺乏可靠、完整的绳索阻尼数据库可供查询。
阻尼的理论计算方法有复数特征值法、模态应变能法、有限单元法和Rayleigh-Ritz法等[8-12]，但是这些方法基于线性黏弹性振动理论，即假定振动过程中张力、几何形状等不变，忽略了阻尼的非线性特性，且计算方法复杂，无法高效准确地对阻尼进行评价。现有文献对绳阻尼的研究均针对粗绳(直径＞2 mm)[13-16]，由于阻尼的非线性，文献[17-19]中提到的针对粗绳阻尼的理论计算方法不一定能适用于细绳(直径≤2 mm)，因此本文采用实验方法研究细绳的阻尼特性。
目前，绳系并联机构具有工作空间大、负载能力强、刚度高等特点，已成为国内外的研究热点[4-7]。特别在风洞试验领域，绳系并联机构作为一种新型支撑方式，与传统的硬式支撑相比，具有刚度大、对流场干扰小、易实现高速及复杂规律的运动等优点[20]，近年来引起国内外学者广泛关注。然而，在现有文献中，大部分未考虑绳阻尼的影响，个别虽考虑了绳阻尼，但绳阻尼系数是虚构值。汪选要等[21-22]将柔索简化成刚度为常值的弹簧，建立了并联柔索机构的动力学模型，进行了轨迹跟踪控制的研究，但没有考虑绳阻尼的影响。刘欣等[23]考虑了绳的弹性变形，进行了绳牵引并联机器人的运动控制研究，但没有考虑绳阻尼的影响。文献[24]的研究表明，绳阻尼对末端执行器的位姿具有显著的影响，但绳阻尼系数是虚构值。文献[25]考虑绳索弹性、阻尼和拉伸效应，给出了绳牵引并联支撑系统的振动特性，但绳阻尼比是虚构值。高估或低估绳索阻尼，将给绳系并联机构的设计带来较大误差，甚至影响机构整体的安全性。因此，研究绳索的阻尼特性，并分析绳阻尼对绳系并联机构的动力学特性的影响具有重要意义。
不失一般性，本文以应用于风洞试验新型支撑的绳系并联机器人(Wire-driven Parallel Robot，WDPR)为例，研究绳阻尼对WDPR在风洞试验中的动力学特性的影响。本文提出了测量绳阻尼的实验方法，解决了细绳阻尼的测量困难；通过采用绳阻尼比的实验测量值，建立WDPR的有阻尼动力学方程，分析风洞来流冲击作用下WDPR的动力学响应，给出了绳阻尼对WDPR动力学特性的影响判据。
1. 绳阻尼特性测试
阻尼比的定义是阻尼系数与临界阻尼系数之比[26]，用于表达结构体标准化的阻尼大小。根据阻尼的产生机制[27]，本文研究绳的内部阻尼，即来自于材料内部各种微观和宏观过程的机械能耗散，以下简称阻尼。
1.1 实验原理
目前，对绳阻尼的测量一般采用激光位移计或应变片测量的方式[15, 18-19]，而这2种方法均不适用于细绳(直径≤2 mm)的阻尼测量。因为绳太细，激光位移计的测量变得十分困难甚至不可用；粘贴应变片，则会改变细绳的动特性，影响阻尼测量结果的准确性。本文采用高速相机测量的方式，属于非接触式的测量方式，由于无需在绳索上附加任何质量元件，绳索本身的特性不被改变，特别适用于细绳(直径≤2 mm)阻尼的测量，同时也适用于粗绳(直径＞2 mm)。
实验原理如图 1(a)所示。实验时先将绳索张紧，在绳索中央悬挂重物施以载荷；待重物稳定后突然释放重物，即对绳索中央施加瞬间激励，使其做自由衰减振动；采用高速相机拍摄绳索中央部位的振动过程；从所拍绳索振动视频中提取绳索振幅最大处的位移信息进行数据处理和分析，得出与绳索阻尼特性相关的参数。通过该实验装置可测量不同长度、直径、预紧力和材料的绳索的阻尼参数。
绳阻尼实验装置的实物图如图 1(b)所示。实验台架采用欧标6060L双槽工业铝型材搭建，以保证实验台架具有良好的强度、刚度和稳定性；并以水平仪校核其水平度和垂直度。采用型号为CHC-S的拉力传感器(量程为200 N，精度为0.1%)测量绳张力。采用高速相机(型号：FASTCAM Mini AX200)获取绳索振动位移。标定板用于在数据后处理中确定绳索在振动中的实际位移。
如图 2所示，选择3种不同直径的Kevlar绳作为实验样本，其物性参数如表 1所示。通过设置不同的绳索参数(直径d、预紧力作用下的初始绳长L0和预紧力T0)进行实验，得到不同绳索参数下的阻尼比。
表 1 绳样本的物性参数Table 1. Physical parameters of wire samples样本序号 绳直径d/mm 材料 体密度ρw/(kg·m-3) 初始绳长L0/m 预紧力T0/N 1 0.6 Kevlar 1 440 0.6~1.4 20~150 2 1.3 3 2 1.2 数据处理方法
采用对数衰减法进行数据后处理。记绳中点振动位移响应为dy。图 3为实验数据dy随时间衰减的曲线。
图 3 绳在竖直面内的振动位移响应(d=0.6 mm，L0=1 m，T0=30 N)Figure 3. Vibration displacement response of wire in vertical plane(d=0.6 mm, L0=1 m, T0=30 N)如图 3所示，Ai代表绳振动衰减曲线dy的一个峰值，Ai+r代表第r个周期后的峰值，可得[27]
(1) 则对数衰减率δ为
(2) 由式(2)可得绳阻尼比ζ的表达式为
(3) 通过实验测得绳索的振动位移响应曲线，利用式(2)和式(3)进行数据处理，即可得到绳阻尼比。
1.3 实验结果
不同直径、不同初始绳长的Kevlar绳在不同预紧力作用下的阻尼实验结果如图 4~图 6所示。
图 4 绳阻尼比与预紧力的关系(d=0.6 mm，L0=1 m)Figure 4. Relationship between wire damping ratio and preload(d=0.6 mm, L0=1 m)
图 5 绳阻尼比与初始绳长的关系(d=0.6 mm)Figure 5. Relationship between wire damping ratio and wire length (d=0.6 mm)
图 6 绳阻尼比与绳直径的关系(L0=1 m)Figure 6. Relationship between wire damping ratio and wire diameter (L0=1 m)从图 4可见，绳阻尼比随预紧力的增加而减小。T0∈[20, 80] N时，绳阻尼比的下降速率很快；T0∈[90, 145] N时，绳阻尼比的下降速率趋于平缓，且绳阻尼比随预紧力的变化呈现非线性。
图 5表明，绳阻尼比随初始绳长的增加而增大，但曲线变化较为平坦。在不同预紧力作用下，绳阻尼比随初始绳长的变化速率基本相同。预紧力为80 N时，绳阻尼比的分散度更小。
从图 6可见，绳阻尼比随绳直径的增加而增大。当预紧力为30 N时，随着绳直径的增加，绳阻尼比明显增大；但当预紧力为80 N时，随着绳直径的增加，绳阻尼比的增大很不明显。因此，预紧力越小，绳阻尼比随绳直径的变化将不能忽略。
2. 动力学建模
2.1 WDPR动力学建模
图 7(a)为八绳牵引的六自由度WDPR原理样机示意图[28]。WDPR原理样机采用8根Kevlar绳做牵引绳，将飞机模型悬挂在空中；改变牵引绳的长度可对模型六自由度的运动进行控制。所建造的WDPR原理样机如图 7(b)所示。
如图 8所示，以静坐标系OXYZ为参考系，飞机模型的位姿记为X =[XP, YP, ZP, ф, θ, ψ]T，(XP, YP, ZP)为沿3个坐标轴的平动，(ф, θ, ψ)为绕3个坐标轴的转动。
在静坐标系OXYZ下，第i根绳的绳长矢量定义为
(4) 式中：Bi为第i根绳与滑轮的连接点，Pi为第i根绳与飞机模型的牵引点，
记Li为第i根绳的实时绳长，则有
(5) 根据Newton-Euler法，飞机模型的动力学方程如下：
(6) 式中：
在风洞试验中，作用在飞机模型上的气动力和力矩有升力、阻力、横向力、俯仰力矩、偏航力矩和滚转力矩，其表达式如下[29]：
(7) 式中：负号表示与静坐标系坐标轴方向相反；q为风洞试验段来流动压，对于低速风洞，
2.2 绳张力建模
考虑绳阻尼，采用弹簧模型对绳索进行建模，得到绳张力的表达式如下：
(8) 式中：T为绳张力；ku为未变形绳的刚度；c为绳阻尼系数；L为实时绳长；Lu为未变形的绳长。
根据绳阻尼系数和绳阻尼比的定义有
(9) 式中：ζ为绳阻尼比，其数值可通过第1节的实验测得；ms为绳单位长度质量。
(10) 
(11) 式中：E为绳弹性模量；A为绳未变形时的横截面积；ρw为绳的体密度。
再根据应变公式有
(12) 联立式(9)~式(12)，可得绳阻尼系数为
(13) 绳长变化速率
(14) 式中：
所以，绳张力矢量表达式为
(15) 式中：Ku=diag(ku1, ku2, …, ku8)；Cw=diag(c1, c2, …, c8)；L =[L1, L2, …，L8]T；Lu=[Lu1, Lu2, …，Lu8]T。
2.3 WDPR的有阻尼动力学方程
联立式(6)、式(14)和式(15)，可得WDPR的有阻尼动力学方程为
(16) 式中：
根据式(4)可知，绳长向量L随飞机模型位姿X的变化而变化；根据式(6)可知，哥氏离心力N与飞机模型绕质心转动的角速度矢量ω呈非线性关系；Jacobi矩阵JAT也会随飞机模型位姿X的变化而变化。因此，式(16)是一个非线性动力学方程。
3. 绳阻尼对WDPR动力学特性的影响
本节基于式(16)给出的WDPR的有阻尼动力学方程，研究风洞来流作用下，WDPR的动力学响应受绳阻尼的影响。
3.1 仿真条件
仿真中所用的绳参数及飞机模型参数如下：
1) 绳参数。如表 1所示，采用3种不同直径的Kevlar绳做牵引绳，为了便于比较，弹性模量统一取E=21.9 GPa。
2) 飞机模型采用SDM标模，模型质量m=1.093 kg，机翼参考面积S=0.026 594 m2，平均气动弦长cA= 0.092 m，模型关于质心的惯量矩阵为
3) 来流条件。以纵向测力试验[29](模型侧滑角为0°，在一系列攻角下进行测量)为例，来流速度V=17 m/s，空气密度ρ=1.29 kg/m3，飞机模型攻角12°，升力系数CL、阻力系数CD、俯仰力矩系数Cm见参考文献[28]。
为了克服非线性动力学方程的求解困难，本文采用隐式的变阶Runge-Kutta数值积分方法，对式(16)进行求解。
3.2 结果与分析
通过仿真分析，得到WDPR的动力学响应曲线，如图 9~图 12所示。其中，飞机模型的位姿响应以飞机模型沿OX方向的位姿变化和俯仰角变化为例。
如图 9所示，当绳直径为0.6 mm时，不考虑阻尼的情况下，飞机模型沿OX方向的位姿变化的峰-峰值为0.1 mm；考虑阻尼的情况下，飞机模型沿OX方向的位姿变化的初始峰-峰值也为0.1 mm，且随时间缓慢变小，但变化幅度不大。当绳直径为2 mm时，不考虑阻尼的情况下，飞机模型沿OX方向的位姿变化的峰-峰值仅为0.01 mm；考虑阻尼的情况下，飞机模型沿OX方向的位姿变化的初始峰-峰值也为0.01 mm，且随时间快速变小，在t=3 s后趋于稳定。
如图 10所示，当绳直径为0.6 mm时，不考虑阻尼的情况下，飞机模型俯仰角变化的峰-峰值为0.1°；考虑阻尼的情况下，飞机模型俯仰角变化的初始峰-峰值也为0.1°，且随时间缓慢变小，但变化幅度不大。当绳直径为2 mm时，不考虑阻尼的情况下，飞机模型俯仰角变化的峰-峰值仅为0.01°；考虑阻尼的情况下，飞机模型俯仰角变化的初始峰-峰值为0.006°，且随时间快速变小，在t=2 s后趋于稳定。
在来流作用下，飞机模型位姿变化的同时，绳张力也随之变化。令ΔT=T-T0，ΔT为绳张力的变化量，T为飞机模型位姿变化过程中的实时绳张力。以绳5为例，如图 11所示，对于不同粗细的绳，不考虑阻尼的情况下，绳张力变化的峰-峰值为2 N；考虑阻尼的情况下，绳张力变化的初始峰-峰值也为2 N。在考虑阻尼的情况下，直径为0.6 mm的绳的绳张力随时间缓慢变小；而直径为2 mm的绳的绳张力随时间快速变小，在t=2 s后趋于稳定。
综上所述，在来流作用下，采用直径较大的绳，飞机模型位姿变化和绳张力变化在短时间内能趋于稳定。因此，在绳对流场的影响可以接受的情况下，采用直径较大的绳，WDPR的稳定性更好。
图 12为飞机模型位姿响应和绳张力响应曲线的局部放大图。从图 12可见，无阻尼情况下，飞机模型的位姿和绳张力呈周期性振荡；考虑阻尼的情况下，飞机模型的位姿和绳张力呈衰减振荡，且绳直径较大时，振荡衰减得越快。说明绳阻尼对WDPR的动力学响应具有减振作用，且绳直径越大，绳阻尼的减振作用越明显。特别地，对于直径小于2 mm的细绳来说，相比无阻尼的情况，绳阻尼对飞机模型位姿及绳张力的振荡频率和相位的影响很小，基本可以忽略。
3.3 绳阻尼的影响判据
根据上述分析，直径不同的绳，对WDPR动力学响应的影响明显不同。以绳5为例，根据式(13)，不同直径绳的阻尼系数如表 2所示。由于阻尼比和绳直径不同，导致阻尼系数有很大差别，从而导致WDPR的动力学响应明显不同。
表 2 不同直径绳的阻尼系数(绳5)Table 2. Damping coefficient of wire with different diameters(Wire No.5)绳直径d/mm 初始绳长L0/m 预紧力T0/N 阻尼比ζ 阻尼系数c/(N·s·m-1) 0.6 0.685 30 0.034 0.131 2 0.685 30 0.075 3.198 直径不同的绳，其阻尼比不同，并且初始绳长、初始绳张力的不同，都会影响绳的阻尼系数。因此，以绳阻尼系数作为因变量，取2种不同直径的绳，分析其阻尼系数的变化对WDPR动力学响应的影响。定义当量系数ηX、ηP、ηT, 对飞机模型位姿和绳张力进行无量纲化处理：
(17) 式中：pX0为无阻尼时飞机模型沿OX方向位姿变化的峰-峰值；pXc为有阻尼时飞机模型沿OX方向位姿变化5 s后的峰-峰值；pP0为无阻尼时飞机模型俯仰角变化的峰-峰值；pPc为有阻尼时飞机模型俯仰角变化5 s后的峰-峰值；pT0为无阻尼时绳张力变化的峰-峰值；pTc为有阻尼时绳张力变化5 s后的峰-峰值。
如图 13所示，当绳阻尼系数c＞0.6 N·s/m时，无阻尼与有阻尼的动力学响应峰-峰值之比将超过10，即绳阻尼的减振作用变得很显著，特别是直径粗的绳比直径细的绳，其阻尼的作用更明显。因此，在WDPR的设计中，可根据式(13)和实测的绳阻尼比计算得到，当绳阻尼系数c满足以下条件：c＞0.6 N·s/m时，不论绳直径粗细如何，其阻尼对WDPR动力学特性的影响不能忽略。
图 13 绳阻尼系数变化对WDPR动力学响应的影响Figure 13. Influence of wire damping coefficient variation on dynamic response of WDPR4. 结论
本文以应用于风洞试验模型支撑的绳系并联机器人WDPR为例，研究绳阻尼对WDPR动力学特性的影响，得出以下结论：
1) 考虑绳阻尼的情况下，当绳直径较小(d＜1 mm)时，飞机模型位姿变化和绳张力变化的峰-峰值随时间缓慢变小，但变化幅度不大；当绳直径较大(d≥2 mm)时，飞机模型位姿变化和绳张力变化的峰-峰值随时间快速变小，并在3 s后趋于稳定。因此，在绳对流场的影响可以接受的情况下，采用直径较大的绳，WDPR的稳定性更好。
2) 绳阻尼对WDPR动力学响应的影响主要体现在响应幅值上，绳直径越大，绳阻尼对WDPR动力学响应的减振作用越明显。特别地，对于直径小于2 mm的细绳来说，相比无阻尼的情况，绳阻尼对飞机模型位姿及绳张力的振荡频率和相位的影响很小，基本可以忽略。
3) 当绳阻尼系数大于0.6 N·s/m时，不论绳直径粗细如何，其阻尼对WDPR动力学特性的影响不能忽略。
综上所述，本文提出的建模方法和研究结果可为WDPR的设计提供指导。此外，绳阻尼的非线性特征及其随不同参数的变化规律，将在今后进一步开展研究。
-
[1] 魏名山.汽车与环境[M].北京:化学工业出版社,2005:6-7 Wei Mingshan.Automobile and the environment[M].Beijing:Chemical Industry Press,2005:6-7(in Chinese) [2] 赵清,徐衍亮,安忠良,等.电动汽车的发展与环境保护[J].沈阳工业大学学报,2000,22(5):430-432 Zhao Qing,Xu Yanliang,An Zhongliang,et al.Development of electric vehicles and environmental protection[J].Journal of Shenyang University of Technology,2000,22(5):430-432(in Chinese) [3] 马国远,史保新,陈观生,等.电动汽车热泵空调系统的试验研究[J].低温工程,2000(4):40-44 Ma Guoyuan,Shi Baoxin,Chen Guansheng,et al.Experimental study of heat pump system for electric vehicle air conditioning[J].Cryogenics,2000(4):40-44(in Chinese) [4] 张优,张海东,黄绍军,等.纯电动汽车空调制热方案探讨[C]//SAE-C2010160.北京:机械工业出版社,2010 Zhang You,Zhang Haidong,Huang Shaojun,et al.Discussing the heating schemes of air-conditioning in pure electric vehicle[C]//SAE-C2010160.Beijing:China Machine Press,2010(in Chinese) [5] 王晓丹.纯电动客车空调系统参数匹配与设计研究[D].长春:吉林大学,2009 Wang Xiaodan.Study on design and parameter matching for air conditioning system of electric bus[D].Changchun:Jilin University,2009(in Chinese) [6] 宁秋宇.电动车空调系统的匹配设计[D].长春:吉林大学,2011 Ning Qiuyu.Matching design for air conditioning system of electric vehicle[D].Changchun:Jilin University,2011(in Chinese) [7] 陈清泉,孙逢春,祝嘉光.现代电动汽车技术[M].北京:北京理工大学出版社,2002:249 Chen Qingquan,Sun Fengchun,Zhu Jiaguang.Modern electric vehicle technology[M].Beijing:Beijing Institute of Technology Press,2002:249(in Chinese) [8] Tadashi O,Imanishi Y,Sekiya S.Thermal management system for electric vehicles[J].SAE International Journal of Materials and Manufacturing,2011,4(1)：1277-1285 [9] Rugh J P.Proposal for a vehicle level test procedure to measure air conditioning fuel use[C]//SAE Technical Papers.Warrendale,PA:SAE,2010 [10] 姚杨,马最良.浅议热泵定义[J].暖通空调,2002,32(3):33 Yao Yang,Ma Zuiliang.Notion of heat pump definition[J].Heating,Ventilating & Air Conditioning,2002,32(3):33(in Chinese) [11] 谢卓,陈江平,陈芝久.电动车热泵空调系统的设计分析[J].汽车工程,2006,28(8):763-765 Xie Zhuo,Chen Jiangping,Chen Zhijiu.On the design of heat pump air conditioning system for electric vehicles[J].Automotive Engineering,2006,28(8):763-765(in Chinese) [12] Hosoz M,Direk M.Performance evaluation of an integrated automotive air conditioning and heat pump system[J].Energy Conversion and Management,2006,47(5):545-559 [13] 史保新,马国远,陈观生.电动车用空调装置的研究[J].流体机械,2002,30(4):48-50 Shi Baoxin,Ma Guoyuan,Chen Guansheng.Research on heat pump system for electric vehicle air conditioning[J].Fluid Machinery,2002,30(4):48-50(in Chinese) [14] 魏名山,彭发展,李丽,等.一种变频电机驱动的电动汽车热泵空调系统:中国,201220659508.9[P].2013-05-15 Wei Mingshan,Peng Fazhan,Li Li, et al.A heat pump air conditioner drived by a variable-frequency motor for electric vehicle:China,201220659508.9[P].2013-05-15(in Chinese) [15] 陈万仁,王保东.热泵与中央空调节能技术[M].北京:化学工业出版社,2010:16 Chen Wanren,Wang Baodong.Energy-saving technology of heat pump[M].Beijing:Chemical Industry Press,2010:16(in Chinese) 期刊类型引用(4)
1. 王晓光，林麒. 基于绳牵引并联支撑的风洞复杂动态试验技术研究进展. 厦门大学学报(自然科学版). 2025(02): 213-224+276 . 
百度学术2. 吴彦宝，乔立春，陈楠，郭阳. CTU机器人装载轨迹自适应模糊滑模控制方法. 自动化与仪表. 2023(06): 51-56 . 百度学术3. 吴惠松，林麒，柳汀，刘震，师璐，王晓光. 风洞虚拟飞行试验模型绳系并联支撑机构. 航空学报. 2022(08): 314-326 . 百度学术4. 乔栋，谢亚龙，李博文，姚涛. 基于Matlab的并联机器人机构动力分析与动态仿真. 现代电子技术. 2021(24): 182-186 . 百度学术其他类型引用(3)
-
下载:
百度学术
点击查看大图
计量
- 文章访问数: 1726
- HTML全文浏览量: 131
- PDF下载量: 843
- 被引次数: 7
下载:
