洪昊 赵志康 陈曦 何韩军

摘要：直线压缩机具有体积小、摩擦力小、效率高以及无油润滑等优点，其小型化设计在微電子冷却及微型空调领域具有重要的发展前景。利用Ansoft Maxwell 有限元分析软件建立了动磁式直线电机模型，确定了谐振工况下直线电机结构及性能参数，通过数值积分分析了空载和负载工况下的活塞位移。根据理论分析和模拟，研制了对置式直线压缩机实验样机，压缩机总质量为1.2 kg，总长度为118 mm。实验结果显示：以 R134a 为制冷剂时，在冷凝压力为0.79 MPa、蒸发压力为0.39 MPa 下，制冷量为57.3 W。模拟及实验结果对直线压缩机小型化设计具有参考价值。

关键词：直线压缩机;动磁式直线电机;有限元分析; R134a

中图分类号：
TF457    文献标志码：
A

Simulation and experimental research on small opposed linear compressor

HONG Hao，ZHAO Zhikang，CHEN Xi，HE Hanjun

（School of Energy and Power Engineering， University of Shanghai for Science and Technology， Shanghai 200093， China）

Abstract：Linear compressor has the advantages of small volume， small friction， high efficiency， and oil-free lubrication. Its miniaturization design is promising in the fields of microelectronic cooling and micro air-conditioning. A moving-magnet motor model was developed by the finite element analysis software of Ansoft Maxwell， by which the structure of the linear motor and performance under resonance conditions were determined. And the piston displacement under no- load and load conditions was analyzed by numerical integration. The experimental prototype of opposed linear compressor was developed based on the theoretical simulation. The compressor mass  was 1.2 kg  and  its  length  was 118 mm. The  experimental  results  showed  that  when evaporation pressure and condensation pressure were 0.39 MPa and 0.79 MPa， respectively with R134a  as  refrigerant，  the  cooling  capacity  of this  system  was 57.3 W. The  simulation  and experimental results can provide references for the miniaturization design of linear compressor.

Keywords：linear compressor; moving-magnet linear motor; finite element analysis; R134a

压缩机是蒸气压缩制冷系统的核心部件，它决定整个制冷系统的制冷量及性能优劣。目前，市场上的制冷设备大多采用传统往复式压缩机。该类型压缩机技术相对成熟，加工成本较低，但其内部的曲柄连杆机构导致压缩机整体系统结构庞大，同时在动力传递过程中，各部件间相互摩擦，造成传动效率降低、可靠性下降以及噪音与振动较大等不良影响。直线压缩机由直线电机直接驱动活塞运动，结构紧凑，无曲柄连杆机构，采用间隙密封技术实现了无油润滑，可替代往复式压缩机以提高制冷系统制冷性能和效率[1–3]。

直线电机小型化及其模拟设计是目前直线压缩机一个重要研究方向，关系到直线压缩机在航空航天、微电子、民用制冷设备等领域的应用与发展。

Bradshaw 等[4]为了测试线性压缩机适用于电子冷却的可行性，研发了一台小型线性压缩机，探讨了压缩机几何形状变化的敏感性，指出泄漏间隙和活塞偏心为线性压缩机设计时须考虑的重要参数。毕研强等[5]通过数值分析求解了动磁式直线压缩机的数学模型，得到压缩机启动特性以及制冷量和制冷系数随线圈电阻、动子质量、摩擦阻尼等参数的变化情况。

Oliveira 等[6]以Embraco公司设计的“智驱”小型直线压缩机为模型建立无油直线压缩机的热力学模型，通过有限体积法计算固体部件的导热以及压缩机壳体内的气体流动状况，通过模拟得到的吸气温度分布同实验结果基本吻合，但一些固体部件和排气温度分布与实验相比有较大差异。唐明生等[7]通过对比测量值与实验测试值，对压缩机动力学参数计算方法进行校正，给出了直线压缩机动力学参数的测量方法，并在线测量动力学参数以实现直线压缩机运行状态的故障诊断及在线监测。

为研究小型动磁式直线压缩机运行特性，本文对系统运行方程组进行数值求解，分析压缩机活塞位移响应。采用Ansoft Maxwell 电磁场有限元二维瞬态场分析软件建立动磁式直线电机模型，依据电机性能模拟结果研发对置式直线压缩机样机，搭建性能测试实验台并进行空载及负载实验。

1理论模型及位移响应分析

动磁式直线压缩机动子部分由永磁体、永磁体支架、板弹簧、连接件及活塞组成，其轴向剖面结构如图1所示。板弹簧为压缩机活塞提供径向支撑力和轴向往复力，定子部分由内轭铁、外轭铁、气缸以及线圈组成。压缩机的工作磁场主要有两个：一个是永磁体产生的恒定磁场，另一个是通电线圈产生的交变磁场，在两个磁场相互作用下，电磁驱动力推动压缩机活塞在轴向上作往复直线运动[8]。

根据法拉第电磁感应定律，在直线电机电路中产生的感应电动势，一部分是由于磁场强度变化产生涡旋电场导致电流定向运动的感生电动势，另一部分是永磁体在内、外轭铁间作直线往复运动产生的动生电动势。在线圈两端施加交变

电压u（t）时，根据基尔霍夫定律以及电压平衡关系，可得到电机驱动电压控制方程[9]，即

u（t）= iRe+ Le + K0（1）

式中：i为线圈电流， A ;Re 为线圈电阻，Ω; Le为线圈等效电感，H;K0为比推力系數，N · A-1; t为时间， s;x为活塞位移， m。

动磁式直线压缩机机械弹簧共振系统可简。在压缩机运行过程中，活塞受到电机电磁力、机械阻尼力、机械弹簧力及气体力作用，外力间的矢量关系如图2所示，图中：FMT（t）为电机电磁力， N ;FML（t）为机械弹簧力， N ;FAM（t）为机械阻尼力， N ; FG（t）为气体力， N ;x.为活塞速度， m · s?1;为活塞加速度，m · s?2;m为动子质量，kg。

根据力平衡关系可得到活塞控制方程[11]，即

电机驱动电压控制方程与活塞控制方程一同构成动磁式直线压缩机系统运行方程组，其机械弹簧共振系统中的势能和动能分别与动子的位移和速度相关，电磁系统中的磁能与励磁线圈中的电流相关，故选取活塞控制方程中的i、x（t）和速度v（t）作为系统状态变量，可得到适用于数值积分计算的系统状态方程组[12]，即

式中：c为等效阻尼系数，包括活塞运动过程中的线性阻尼系数，N ·（m · s）-1;k为等效弹簧刚度，包括板弹簧机械弹簧刚度和气体弹簧刚度，N · m-1。

通过求解系统状态方程组，可得到在一定输入电压条件下压缩机的活塞位移响应。方程组中已知常量参数可由式（4）得到，即式中：
Wi为压缩机指示功， J ;X为活塞振幅， m;u0为真空磁导率， H · A-1;D为永磁体与轭铁间隙的平均直径， m; L为轭铁轴向长度， m; n为永磁体轴向长度， m; N 为线圈匝数;Hc为矫顽力，A · m-1;gn为内轭铁与永磁体间隙宽度， m;gw为外轭铁与永磁体间隙宽度， m;ω为活塞运动的角频率，rad · s-1;w 为永磁体的宽度， m;f 为共振频率， Hz。

求解常微分方程一般采用龙格?库塔法（Runge?Kutta）。该方法是一种广泛应用于工程计算的数值积分算法，其优点在于计算精度、效率较高，稳定性好。在Matlab软件中采用龙格?库塔法进行数值分析求解，计算参数如表1所示。

压缩机空载时活塞不受气体力作用，运动平衡位置不发生偏移。图3（a）为输入电压有效值为6 V 时压缩机活塞位移响应，在电机开始运行后数个周期内活塞位移逐渐增加，运动中心并未偏移，之后达到行程为9 mm 的稳定运行工况。图3（b）为负载工况下输入电压有效值为23 V 时的活塞位移响应。从图中可以看出，当压缩机稳定运行时，活塞运动平衡位置发生偏移。分析其原因是因为活塞两侧压差产生的额外偏置力与机械弹簧力相平衡，使活塞往复位移中点偏离活塞启动位置。活塞偏置会造成余隙容积增大，压缩机效率降低。为消除活塞偏置，可采用在电源输入交变电压上增加定额电压的活塞行程控制方法[13]，附加电动势产生的电磁力抵消了活塞偏置力，使活塞位移中心重新回到初始启动位置。

2电磁场有限元分析

在Ansoft Maxwell 有限元分析软件中建立二维模型，模拟结构参数与实验样机一致。电机固有频率与驱动电源频率相等时即为谐振，此时在最小的驱动力下活塞能够达到最大位移。图4为电机在100 Hz 的谐振状态下活塞位移、电流和电磁力随时间的变化。当电流达到峰值，电磁力也达到最大值，活塞处于平衡位置，故电流与电磁力同相位，电流与活塞位移间的相位角为90°。

动子在位移平衡位置附近所受电磁力较大，瞬时磁感应强度大，轭铁易达到磁饱和点，从而影响压缩机稳定运行。图5为电机稳定运行中永磁体接近平衡位置时磁感应强度云图，内轭铁和外轭铁内部最大磁感应强度低于1.5 T。模拟中选用电工纯铁作为轭铁材料，其磁通密度在1.9 T 左右达到饱和，故整个工况下轭铁磁通密度未饱和，磁路稳定。

3性能测试实验与结果分析

3.1实验系统

根据理论分析及电磁场有限元模拟结果开发了小型双活塞对置式制冷压缩机样机，实物如图6所示。该压缩机整机总质量为1.2 kg，最大直径为62 mm，长度为118 mm。实验需测试样机在制冷系统中的运行工况参数，包括压缩机输入电压、驱动电源频率、压缩活塞行程、吸气压力、排气压力及制冷量等关键性能参数，故搭建了较完备的性能测试与参数调节系统。直线压缩机制冷性能测试系统如图7所示。

3.2空载特性实验

压缩机空载运行时理论上无气体力作用，系统刚度仅为板弹簧机械刚度，系统共振频率发生改变。共振频率计算式为

图8为空载工况下输入电压不变时单个直线电机活塞振幅、电机效率和输入功率随驱动电源频率的变化。随着驱动电源频率的增加，输入功率先减小后增大，电机效率和活塞振幅先增大后减小，在频率为74 Hz 时，输入功率最低，为17 W，电机效率也达到最大值，为0.82。在75 Hz 左右活塞振幅达到最大值，该频率与电机效率最大和输入功率最小时的频率基本吻合，可认为压缩机达到谐振状态。根据共振频率计算式得到的频率为71 Hz，稍低于实验共振频率。其原因是：由于压缩过程中吸气、排气阀片开启与关闭造成活塞两侧压力不相等，气体弹簧发挥作用，系统刚度大于板弹簧刚度，故实验共振频率略大于理论计算值。

3.3负载特性实验

图9～11分别为制冷量，压缩机吸气和排气压力，压缩机吸气和排气温度随运行时间的变化。由于压缩机吸气管道较短，过热度小，故压缩机吸气温度近似为蒸发温度。实验中施加了4.3 V 的偏置电压以抵消活塞偏置造成的不利影响，并在冷凝压力为0.79 MPa、蒸发压力为0.39 MPa 时测得最大制冷量为57.3 W，压缩机吸气温度降至7.3 C，此时的 COP 为1.15。实验中线性压缩机的排气温度持续上升导致制冷效果变差。其原因在于：电机绕组线圈电流负荷较高，发热量较大，热量传递给压缩腔导致排气温度上升。通过为直线电机配置主动冷却设备可有效控制排气温度。实验结果与理论分析及模拟结果有一定差距。其原因在于：压缩机实际机械阻尼大于电机模拟中的机械阻尼，活塞与气缸间摩擦较严重，另外，压缩腔内高压气体向背压腔泄漏导致压比降低。该不良影响主要由两方面原因造成：吸气、排气阀片延迟开启和闭合，导致气体回流，排气压力减小，背压腔压力升高;活塞磨损与偏置导致活塞与气缸之间間隙变大，间隙密封效果变差。

4结论

为研究小型动磁式直线压缩机的制冷性能，开发并搭建了直线压缩机制冷性能测试系统。测试了小型对置式直线压缩机空载和负载工况下的直线电机性能，电机达到谐振状态时活塞振幅与电机效率最大，输入功率最小。负载实验结果表明，在冷凝压力为0.79 MPa、蒸发压力为0.39 MPa 时，制冷温度为7.3 C，COP 为1.15，系统制冷量为57.3 W。

参考文献：

[1]邰晓亮.动磁式直线电机驱动微小型活塞压缩机理论分析及实验研究[D].上海：上海交通大学， 2009.

[2] LIANG K， STONE R， DADD M， et al. A novel linearelectromagnetic-drive  oil-free  refrigeration  compressorusing R134a[J]. International  Journal  of  Refrigeration，2014， 40：450-459.

[3] 陈曦， 张华， 吴亦农， 等.200 K/40 W 自由活塞斯特林制冷机的实验研究[J].上海理工大学学报 ， 2007， 29（4）：337-344.

[4] BRADSHAW C R， GROLL E A， GARIMELLA S V.Acomprehensive  model  of  a  miniature-scale  linear compressor  for  electronics  cooling[J]. International Journal of Refrigeration， 2011， 34（1）：63-73.

[5] 毕研强， 徐向华， 梁新刚.动磁式直线压缩机对制冷量与 COP 影响分析[J].清华大学学报：自然科学版 ， 2014， 54（7）：884-890.

[6] OLIVEIRA  M  J，  DINIZ  M  C，  DESCHAMPS  C  J.Thermal modelling and analysis of an oil-free linear compressor[C]//Proceedings  of  the 9th  International Conference on Compressors and their Systems. London， UK：
IOP， 2015：012016.

[7] 唐明生， 邹慧明， 王敏， 等.动力学参数在直线压缩机故障诊断中的应用[J].压缩机技术， 2019（3）：21-25，32.

[8] 谢洁飞， 金涛， 童水光.直线压缩机的研究现状与发展趋势[J].流体机械， 2004， 32（12）：31-35.

[9] 雷美珍， 王立强， 夏永明.动磁式双定子直线压缩机动态特性仿真研究[J].低温与超导， 2018， 46（1）：7-11，35.

[10]张博.基于动磁式直线振荡电机的直线缩机技术的研究[D].南京：南京理工大学， 2015.

[11]邓伟峰.大功率动磁式线性压缩机关键技术与实验研究[D].上海：中国科学院研究生院（上海技术物理研究所）， 2016.

[12]谢洁飞.动磁式直线压缩机理论与试验研究[D].杭州：浙江大学， 2005.

[13]崔丽娜.直线压缩机控制系统的开发[D].杭州：浙江大学， 2005.

猜你喜欢 有限元分析 对于Pro/mechanica的数控车刀的有限元分析中国科技纵横(2016年20期)2016-12-28S型压力传感器形变的有限元分析企业技术开发·下旬刊(2016年11期)2016-12-27横向稳定杆的侧倾角刚度仿真与试验研究现代商贸工业(2016年25期)2016-12-26多自由度升降翻转舞台的运动分析演艺科技(2016年11期)2016-12-24高精度细深孔的镗孔加工工艺分析山东工业技术(2016年23期)2016-12-23基于LS—DYNA的某汽车乘客安全气囊支架焊接强度的分析山东工业技术(2016年23期)2016-12-23自锚式悬索桥钢箱梁顶推施工阶段结构分析文艺生活·中旬刊(2016年11期)2016-12-13随机振动载荷下发射装置尾罩疲劳寿命分析航空兵器(2016年4期)2016-11-28有限元分析带沟槽平封头的应力集中科学与财富(2016年15期)2016-11-24飞机起落架支撑杆强度有限元分析科技视界(2016年18期)2016-11-03