发明名称 与直线压缩机最优匹配的惯性管型脉冲管冷指的设计方法
摘要 本发明公开了一种与直线压缩机最优匹配的惯性管型脉冲管冷指的设计方法,分为八步:1)将惯性管型脉冲管制冷机冷指等效成为交流电路;2)根据压缩机参数计算最优匹配下脉冲管冷指阻抗;3)计算压缩机出口处体积流率;4)设置合理的制冷温度,制冷量以及制冷效率目标;5)赋予脉冲管冷指各个部件初始值;6)计算连管入口处阻抗值;7)计算制冷量及制冷效率;8)比较计算所得阻抗值,制冷量以及制冷效率是否与理论最优值相等,若相等,则设计完成,若不等,返回步骤五调整初始参数,重复步骤六到八。本发明所提出的设计方法对于高效的惯性管型脉冲管制冷机的实用化发展具有非常积极的意义。
申请公布号 CN105485954A 申请公布日期 2016.04.13
申请号 CN201510868765.1 申请日期 2015.12.01
申请人 中国科学院上海技术物理研究所 发明人 党海政;谭军
分类号 F25B9/14(2006.01)I 主分类号 F25B9/14(2006.01)I
代理机构 上海新天专利代理有限公司 31213 代理人 郭英
主权项 一种与直线压缩机最优匹配的惯性管型脉冲管冷指的设计方法,其特征在于,所述的设计方法包括以下步骤:步骤一:惯性管型高频脉冲管制冷机包括直线压缩机(1),连管(2),级后冷却器(3),蓄冷器(4),冷端换热器(5),脉冲管(6),热端换热器(7),惯性管(8),气库(9);其中连管(2),级后冷却器(3),蓄冷器(4),冷端换热器(5),脉冲管(6),热端换热器(7),惯性管(8)以及气库(9)组成了脉冲管冷指(10),直线压缩机(1)和脉冲管冷指(10)通过连管(2)相连接;根据电路类比模型,高频脉冲管制冷机中的压力被等效为电动势,体积流率被等效为电流,流阻,流容以及惯性被分别等效为电路中的电阻,电容和电感,整个高频脉冲管制冷机冷指可以等效成为交流电路;步骤二:测量给定直线压缩机(1)中磁体的磁场强度,活塞的面积,活塞的机械阻尼,线圈的长度,线圈的电阻,板弹簧的轴向刚度以及动子质量的大小,直线压缩机(1)与脉冲管冷指(10)匹配后的电机效率的表达式为:<img file="dest_path_FDA0000937280260000011.GIF" wi="1717" he="213" />表达式(1)中的η为直线压缩机(1)的输入电功转换为脉冲管冷指(10)入口处声功的转换效率;|Z<sub>a</sub>|为脉冲管冷指(10)阻抗的幅值,θ为脉冲管冷指(10)阻抗的相位角,B为直线压缩机(1)中磁体的磁场强度;L为线圈长度;A<sub>p</sub>为活塞面积;b为活塞机械阻尼;R<sub>e</sub>为线圈电阻;m为动子质量,ω为角频率;k<sub>x</sub>为板弹簧轴向刚度;基于压缩机电机效率的表达式(1),得到实现最优匹配 下的脉冲管冷指(10)阻抗的幅值,相位角以及运行频率;步骤三:根据直线压缩机(1)的可运行最大活塞行程,设定合适的压缩机活塞行程,并根据活塞表面体积流率的计算表达式(2)得出直线压缩机(1)出口处的体积流率大小:<img file="dest_path_FDA0000937280260000021.GIF" wi="1142" he="93" />表达式(2)中的<img file="dest_path_FDA0000937280260000022.GIF" wi="107" he="79" />为直线压缩机(1)出口处体积流率,A<sub>P</sub>为活塞面积,ω为角频率,X为活塞行程;步骤四:根据实际的应用需求,设置合理的脉冲管冷指(10)的目标制冷温度,制冷量以及制冷效率;步骤五:赋予初值给脉冲管冷指(10)的各个部件,包括连管(2)的横截面面积与长度,级后冷却器(3)的横截面面积、长度及孔隙率,蓄冷器(4)的横截面面积、长度、丝网直径及孔隙率,冷端换热器(5)的横截面面积、长度及孔隙率,脉冲管(6)的横截面面积与长度,热端换热器(7)的横截面面积、长度及孔隙率,惯性管(8)的横截面面积与长度,以及气库(9)的体积;步骤六:赋予充气压力以及气库(9)入口处的体积流率<img file="dest_path_FDA0000937280260000023.GIF" wi="62" he="79" />初始值,利用表达式(3)和表达式(4)计算气库(9)入口的动态压力以及阻抗值:<img file="dest_path_FDA0000937280260000024.GIF" wi="1190" he="126" /><img file="dest_path_FDA0000937280260000025.GIF" wi="1190" he="134" />表达式(3)中p<sub>9</sub>为气库(9)入口处动态压力,γ为绝热系数,P<sub>m</sub>为充气压力,ω为角频率,V<sub>9</sub>为气库(9)的体积,i为虚部,<img file="dest_path_FDA0000937280260000026.GIF" wi="63" he="78" />为气库(9)入口体积流率,表达式(4)中Z<sub>9</sub>为气库(9)的阻抗;利用表达式(5)、表达式(6)和表达 式(7)计算惯性管(8)入口处的动态压力,体积流率以及阻抗:<img file="dest_path_FDA0000937280260000031.GIF" wi="1412" he="150" /><img file="dest_path_FDA0000937280260000032.GIF" wi="1382" he="142" /><img file="dest_path_FDA0000937280260000033.GIF" wi="1204" he="79" />表达式(5)中的p<sub>8</sub>为惯性管(8)入口处动态压力,p<sub>9</sub>为气库(9)入口处动态压力,l<sub>8</sub>为惯性管(8)的长度,ω为角频率,ρ<sub>8</sub>为惯性管(8)中工质气体的密度,i为虚部,A<sub>8</sub>为惯性管(8)的横截面面积,μ为动粘性系数,S<sub>8</sub>为惯性管(8)的横截面周长,δ<sub>v</sub>为粘性渗透深度,<img file="dest_path_FDA0000937280260000034.GIF" wi="96" he="79" />为惯性管(8)中与气库(9)入口距离为x位置处的体积流率,表达式(6)中的<img file="dest_path_FDA0000937280260000035.GIF" wi="59" he="78" />为惯性管(8)入口处的体积流率,<img file="dest_path_FDA0000937280260000036.GIF" wi="70" he="78" />为气库(9)入口处体积流率,γ为绝热系数,P<sub>m</sub>为充气压力,p<sub>x‑8</sub>为惯性管(8)中与气库(9)入口距离为x位置处的动态压力,表达式(7)中的Z<sub>8</sub>为惯性管(8)入口处的阻抗;利用表达式(8)、表达式(9)和表达式(10)计算热端换热器(7)入口处的动态压力,体积流率以及阻抗:<img file="dest_path_FDA0000937280260000037.GIF" wi="1365" he="159" /><img file="dest_path_FDA0000937280260000038.GIF" wi="1350" he="151" /><img file="dest_path_FDA0000937280260000039.GIF" wi="1245" he="86" />表达式(8)中的p<sub>7</sub>为热端换热器(7)入口处动态压力,p<sub>8</sub>为惯性管(8)入口处动态压力,l<sub>7</sub>为热端换热器(7)的长度,ω为角频率,ρ<sub>7</sub>为热端换热器(7)中工质气体的密度,i为虚部,<img file="dest_path_FDA00009372802600000310.GIF" wi="65" he="62" />为热端换热器(7)的孔隙率,A<sub>7</sub>为热端换 热器(7)的横截面面积,r<sub>7</sub>为热端换热器(7)中流阻,<img file="dest_path_FDA0000937280260000041.GIF" wi="101" he="79" />为热端换热器(7)中与惯性管(8)入口距离为x位置处的体积流率,表达式(9)中的<img file="dest_path_FDA0000937280260000042.GIF" wi="61" he="79" />为热端换热器(7)入口处的体积流率,<img file="dest_path_FDA0000937280260000043.GIF" wi="65" he="79" />为惯性管(8)入口处体积流率,γ为绝热系数,P<sub>m</sub>为充气压力,p<sub>x‑7</sub>为热端换热器(7)中与惯性管(8)入口距离为x位置处的动态压力,表达式(10)中的Z<sub>7</sub>为热端换热器(7)入口处的阻抗;利用表达式(11)、表达式(12)以及表达式(13)计算脉冲管(6)入口处的动态压力,体积流率以及阻抗:p<sub>6</sub>=p<sub>7</sub>       (11) <img file="dest_path_FDA0000937280260000044.GIF" wi="1254" he="135" /><img file="dest_path_FDA0000937280260000045.GIF" wi="1101" he="71" />表达式(11)中的p<sub>6</sub>为脉冲管(6)入口处动态压力,p<sub>7</sub>为热端换热器(7)入口处动态压力,表达式(12)中的<img file="dest_path_FDA0000937280260000046.GIF" wi="62" he="78" />为脉冲管(6)入口处的体积流率,<img file="dest_path_FDA0000937280260000047.GIF" wi="65" he="78" />为热端换热器(7)入口处体积流率,l<sub>6</sub>为脉冲管(6)的长度,ω为角频率,A<sub>6</sub>为脉冲管(6)的横截面面积,i为虚部,γ为绝热系数,P<sub>m</sub>为充气压力,表达式(13)中的Z<sub>6</sub>为脉冲管(6)入口处的阻抗;利用表达式(14)、表达式(15)以及表达式(16)计算冷端换热器(5)入口处的动态压力,体积流率以及阻抗:<img file="dest_path_FDA0000937280260000048.GIF" wi="1332" he="151" /><img file="dest_path_FDA0000937280260000049.GIF" wi="1310" he="134" /><img file="dest_path_FDA00009372802600000410.GIF" wi="1220" he="85" />表达式(14)中的p<sub>5</sub>为冷端换热器(5)入口处动态压力,p<sub>6</sub>为脉冲管(6)入 口处动态压力,l<sub>5</sub>为冷端换热器(5)的长度,ω为角频率,ρ<sub>5</sub>为冷端换热器(5)中工质气体的密度,i为虚部,<img file="dest_path_FDA0000937280260000051.GIF" wi="64" he="62" />为冷端换热器(5)的孔隙率,A<sub>5</sub>为冷端换热器(5)的横截面面积,r<sub>5</sub>为冷端换热器(5)中流阻,<img file="dest_path_FDA0000937280260000052.GIF" wi="96" he="77" />为冷端换热器(5)中与脉冲管(6)入口距离为x位置处的体积流率,表达式(15)中的<img file="dest_path_FDA0000937280260000053.GIF" wi="62" he="78" />为冷端换热器(5)入口处的体积流率,<img file="dest_path_FDA0000937280260000054.GIF" wi="70" he="77" />为脉冲管(6)入口处体积流率,γ为绝热系数,P<sub>m</sub>为充气压力,p<sub>x‑5</sub>为冷端换热器(5)中与脉冲管(6)入口距离为x位置处的动态压力,表达式(16)中的Z<sub>5</sub>为冷端换热器(5)入口处的阻抗;利用表达式(17)、表达式(18)以及表达式(19)计算蓄冷器(4)入口处的动态压力,体积流率以及阻抗:<img file="dest_path_FDA0000937280260000055.GIF" wi="1372" he="150" /><img file="dest_path_FDA0000937280260000056.GIF" wi="1430" he="150" /><img file="dest_path_FDA0000937280260000057.GIF" wi="1286" he="86" />表达式(17)中的p<sub>4</sub>为蓄冷器(4)入口处动态压力,p<sub>5</sub>为冷端换热器(5)入口处动态压力,l<sub>4</sub>为蓄冷器(4)的长度,ω为角频率,ρ<sub>4</sub>为蓄冷器(4)中工质气体的密度,i为虚部,<img file="dest_path_FDA0000937280260000058.GIF" wi="68" he="63" />为蓄冷器(4)的孔隙率,A<sub>4</sub>为蓄冷器(4)的横截面面积,r<sub>4</sub>为蓄冷器(4)中流阻,<img file="dest_path_FDA0000937280260000059.GIF" wi="99" he="78" />为蓄冷器(4)中与冷端换热器(5)入口距离为x位置处的体积流率,表达式(18)中的<img file="dest_path_FDA00009372802600000510.GIF" wi="63" he="78" />为蓄冷器(4)入口处的体积流率,<img file="dest_path_FDA00009372802600000511.GIF" wi="70" he="78" />为冷端换热器(5)入口处体积流率,γ为绝热系数,P<sub>m</sub>为充气压力,p<sub>x‑4</sub>为蓄冷器(4)中与冷端换热器(5)入口距离为x位置处的动态压力,g为温度梯度引起的控制源项,表达式(19)中的Z<sub>4</sub>为蓄冷器(4)入口处的阻抗;利用表达式(20)、表达式(21)以及表达式(22)计算级后冷却器(3)入口处的动态压力,体积流率以及阻抗:<img file="dest_path_FDA0000937280260000061.GIF" wi="1277" he="151" /><img file="dest_path_FDA0000937280260000062.GIF" wi="1286" he="134" /><img file="dest_path_FDA0000937280260000063.GIF" wi="1165" he="87" />表达式(20)中的p<sub>3</sub>为级后冷却器(3)入口处动态压力,p<sub>4</sub>为蓄冷器(4)入口处动态压力,l<sub>3</sub>为级后冷却器(3)的长度,ω为角频率,ρ<sub>3</sub>为级后冷却器(3)中工质气体的密度,i为虚部,<img file="dest_path_FDA0000937280260000064.GIF" wi="69" he="61" />为级后冷却器(3)的孔隙率,A<sub>3</sub>为级后冷却器(3)的横截面面积,r<sub>3</sub>为级后冷却器(3)中流阻,<img file="dest_path_FDA0000937280260000065.GIF" wi="95" he="78" />为级后冷却器(3)中与蓄冷器(4)入口距离为x位置处的体积流率,表达式(21)中的<img file="dest_path_FDA0000937280260000066.GIF" wi="58" he="79" />为级后冷却器(3)入口处的体积流率,<img file="dest_path_FDA0000937280260000067.GIF" wi="68" he="78" />为蓄冷器(4)入口处体积流率,γ为绝热系数,P<sub>m</sub>为充气压力,p<sub>x‑3</sub>为级后冷却器(3)中与蓄冷器(4)入口距离为x位置处的动态压力,表达式(22)中的Z<sub>3</sub>为级后冷却器(3)入口处的阻抗;利用表达式(23)、表达式(24)以及表达式(25)计算连管(2)入口处的动态压力,体积流率以及阻抗:<img file="dest_path_FDA0000937280260000068.GIF" wi="1446" he="143" /><img file="dest_path_FDA0000937280260000069.GIF" wi="1414" he="142" /><img file="dest_path_FDA00009372802600000610.GIF" wi="1270" he="92" />表达式(23)中的p<sub>2</sub>为连管(2)入口处动态压力,p<sub>3</sub>为级后冷却器(3)入口处动态压力,l<sub>2</sub>为连管(2)的长度,ω为角频率,ρ<sub>2</sub>为连管(2)中工质气体的密度,i为虚部,A<sub>2</sub>为连管(2)的横截面面积,μ为动粘性系数,S<sub>2</sub>为连管(2)横截面周长,δ<sub>v</sub>为粘性渗透深度,<img file="dest_path_FDA00009372802600000611.GIF" wi="98" he="77" />为连管(2)中与级后冷却器(3)入口距离为x位置处的体积流率,表达式(24)中的<img file="dest_path_FDA00009372802600000612.GIF" wi="63" he="78" />为连管(2)入口处的 体积流率,<img file="dest_path_FDA0000937280260000071.GIF" wi="63" he="79" />为级后冷却器(3)入口处体积流率,γ为绝热系数,P<sub>m</sub>为充气压力,p<sub>x‑2</sub>为级后连管(2)中与级后冷却器(3)入口距离为x位置处的动态压力,表达式(25)中的Z<sub>2</sub>为连管(2)入口处的阻抗;步骤七:对比步骤六中计算得到的连管(2)入口的体积流率与步骤二计算得到的体积流率值是否相等,如果相等,利用表达式(26)和表达式(27)分别计算脉冲管冷指(10)的制冷量及制冷效率:<img file="dest_path_FDA0000937280260000072.GIF" wi="1564" he="103" /><img file="dest_path_FDA0000937280260000073.GIF" wi="1198" he="133" />表达式(26)中的<img file="dest_path_FDA0000937280260000074.GIF" wi="95" he="79" />为脉冲管冷指(10)的制冷量,(PV)<sub>5</sub>为冷端换热器(5)入口处的声功,ε为蓄冷器效率,C<sub>p</sub>等压比热容,P<sub>m</sub>为充气压力,<img file="dest_path_FDA0000937280260000076.GIF" wi="65" he="71" />为蓄冷器(4)入口处体积流率,<img file="dest_path_FDA0000937280260000077.GIF" wi="59" he="72" />为冷端换热器(5)入口处体积流率,π为圆周率,R<sub>g</sub>为气体常数,表达式(27)中的COP为制冷效率,(PV)<sub>2</sub>为连管(2)入口处的声功,然后进行步骤八;如果不相等,跳转至步骤五,调整气库(9)入口处的体积流率<img file="dest_path_FDA0000937280260000075.GIF" wi="59" he="78" />的初始值,并重复步骤六和步骤七;步骤八:比较步骤六中所得的连管(2)入口的阻抗值与步骤一计算所得的最优脉冲管冷指阻抗值,并同时比较步骤七中计算所得的制冷量和制冷效率与步骤三中的目标制冷量和制冷效率,如果均相等,那么设计完成,此时脉冲管冷指(10)各部件的尺寸能实现与直线压缩机(1)的最优匹配。如果有一项不相同或均不相等,返回步骤六,调整脉冲管冷指(10)各个部件的尺寸值,然后重复步骤六至步骤八。
地址 200083 上海市虹口区玉田路500号