时变工况下制冷设备转子压缩机载荷激励辨识方法
专利汇可以提供时变工况下制冷设备转子压缩机载荷激励辨识方法专利检索，专利查询，专利分析的服务。并且本 发明 公开了一种时变工况下 制冷设备 转子 压缩机 载荷 激励辨识方法，其特征是对时变工况下制冷设备转子压缩机排回气管口进行振动测试，获得转子压缩机排回气管口的实测 加速 度响应时间历程 信号 ，构建加速度期望响应矩阵；利用有限元方法构建制冷设备转子压缩机虚拟样机模型，应用虚拟样机模型分析制冷设备转子压缩机在白噪声载荷作用下排回气管口测点的加速度响应，求出传递函数矩阵。由时变工况下转子压缩机排回气管口测点的振动测试结果反算时变工况下制冷设备转子压缩机的载荷激励，利用 迭代 方法消除误差。本发明实现了时变工况下制冷设备转子压缩机载荷激励的辨识分析，为压缩机‑管路系统振动分析和管路结构的可靠性设计及评估等提供依据。，下面是时变工况下制冷设备转子压缩机载荷激励辨识方法专利的具体信息内容。
1.一种时变工况下制冷设备转子压缩机载荷激励辨识方法，其特征是：对时变工况下制冷设备转子压缩机排回气管口进行振动测试，获得制冷设备转子压缩机排回气管口的实测加速度响应时间历程信号，构建加速度期望响应矩阵；利用有限元方法构建制冷设备转子压缩机的虚拟样机模型，并应用虚拟样机模型分析制冷设备转子压缩机在白噪声载荷作用下排回气管口测点的加速度响应，求出传递函数矩阵；由时变工况下转子压缩机排回气管口测点的振动测试结果反算时变工况下制冷设备转子压缩机的载荷激励，并利用迭代方法消除误差；
所述时变工况下制冷设备转子压缩机载荷激励辨识方法是按如下步骤进行：
步骤1、针对采用转子压缩机的制冷设备，在压缩机单体试验台架上或者在空调整机系统上进行时变工况下的排回气管口振动测试，所述转子压缩机包括压缩机本体(1)和储液桶(2)，排气口(3)位于压缩机本体(1)的顶部，回气口(4)位于储液桶(2)的顶部；
步骤2、分别取所述排气口(3)为测点A、取所述回气口(4)为测点B，并以测点A为坐标原点O，以过原点O的竖直方向为Z轴方向，以垂直于Z轴的平面为XOY平面，测点A和测点B的连线在XOY平面上的投影为X轴，Y轴垂直于XOZ平面，确立坐标系OXYZ；
步骤3、经实测并通过低通滤波信号处理获得时变工况下制冷设备转子压缩机测点A在X、Y及Z方向上的加速度响应时间历程信号为ax、ay、az，以及测点B在X、Y和Z方向上的加速度响应时间历程信号分别为bx、by、bz，构建加速度期望响应矩阵Y为：
Y＝[ax ay az bx by bz] (1)
步骤4、利用有限元方法构建制冷设备转子压缩机的虚拟样机模型，按如下过程分析制冷设备转子压缩机在白噪声激励载荷作用下测点A和测点B的试验加速度响应时间历程信号：
步骤a、对制冷设备转子压缩机的虚拟样机模型，在加载点P上分别施加不同方向上的力和力矩，分别是：
在X方向施加均值为零、方差为1的高斯白噪声载荷序列Fx0(k),将Fx0(k)也记为Fx0；
在Y方向施加均值为零、方差为1的高斯白噪声载荷序列Fy0(k),将Fy0(k)也记为Fy0；
在Z方向施加均值为零、方差为1的高斯白噪声载荷序列Fz0(k),将Fz0(k)也记为Fz0；
绕X方向施加均值为零、方差为1的高斯白噪声载荷序列Mx0(k),将Mx0(k)也记为Mx0；
绕Y方向施加均值为零、方差为1的高斯白噪声载荷序列My0(k),将My0(k)也记为My0；
绕Z方向施加均值为零、方差为1的高斯白噪声载荷序列Mz0(k),将Mz0(k)也记为Mz0；
通过时间响应历程分析得到：
在Fx0的作用下测点A的加速度响应的时间历程信号ax1、ay1、az1，以及测点B的加速度响应的时间历程信号：bx1、by1、bz1；
在Fy0的作用下测点A的加速度响应的时间历程信号ax2、ay2、az2，以及测点B的加速度响应的时间历程信号：bx2、by2、bz2；
在Fz0的作用下测点A的加速度响应的时间历程信号ax3、ay3、az3，以及测点B的加速度响应的时间历程信号：bx3、by3、bz3；
在Mx0的作用下测点A的加速度响应的时间历程信号ax4、ay4、az4，以及测点B的加速度响应的时间历程信号：bx4、by4、bz4；
在My0的作用下测点A的加速度响应的时间历程信号ax5、ay5、az5，以及测点B的加速度响应的时间历程信号：bx5、by5、bz5；
在Mz0的作用下测点A的加速度响应的时间历程信号ax6、ay6、az6，以及测点B的加速度响应的时间历程信号：bx6、by6、bz6；
每个载荷序列对应的六组传递函数为：
H11(w)＝FFT(ax1)/FFT(Fx0),H12(w)＝FFT(ay1)/FFT(Fx0).....H16(w)＝FFT(bz1)/FFT(Fx0)；
H21(w)＝FFT(ax2)/FFT(Fy0),H22(w)＝FFT(ax2)/FFT(Fy0)......H26(w)＝FFT(bz2)/FFT(Fy0)；
…….
H51(w)＝FFT(ax5)/FFT(My0),H52(w)＝FFT(ay5)/FFT(My0)......H56(w)＝FFT(bz5)/FFT(My0)；
H61(w)＝FFT(ax6)/FFT(Mz0),H62(w)＝FFT(ay6)/FFT(Mz0)......H66(w)＝FFT(bz6)/FFT(Mz0)；
其中：FFT为快速傅里叶变换函数，利用六组载荷的6×6频响函数构建传递函数矩阵G为：
所述加载点是指转子压缩机的载荷作用点，取所述加载点为压缩机本体1/3高度处圆截面与压缩机轴线的交叉点；
步骤b、利用式(3)得到时变工况下制冷设备转子压缩机的载荷激励：
-1
F＝W·G (3)
其中：
F＝(FFT(Fx),FFT(Fy),FFT(Fz),FFT(Mx),FFT(My),FFT(Mz))
W＝(FFT(ax),FFT(ay),FFT(az),FFT(bx),FFT(by),FFT(bz))
则：
(Fx,Fy,Fz,Mx,My,Mz)＝IFFT(F) (4)
其中：IFFT为快速傅里叶逆变换函数，Fx和Mx为加载点P上X方向的载荷激励，Fy和My为加载点P上Y方向的载荷激励，Fz和Mz为加载点P上Z方向的载荷激励；
步骤c、将所述Fx、Fy、Fz、Mx、My、Mz加载到制冷设备转子压缩机虚拟样机模型中的加载点，通过时间响应历程分析，得到虚拟样机模型中测点A的X、Y和Z方向加速度响应的时间历程分别为ax10、ay10、az10，测点B的X、Y和Z方向加速度响应的时间历程分别为bx20、by20、bz20，构建如式(5)所示的加速度响应矩阵Y0(k)，也记为Y0，
Y0＝[ax10 ay10 az10 bx20 by20 bz20] (5)
步骤d、将所述加速度期望响应矩阵Y和加速度响应矩阵Y0按式(6)进行比较：
其中：δ为设定的允许误差,k＝0,1,2……
若满足式(6)，步骤(b)中所获得的载荷激励Fx、Fy、Fz、Mx、My和Mz即为所要求解的时变工况下的制冷设备转子压缩机载荷激励；
若不能满足式(6)，则继续进入步骤(e)；
步骤e、利用式(7)获得更新的制冷设备转子压缩机载荷激励F1：
-1
F1＝IIFT(F)+IFFT(G ×FFT(Y-Y0)) (7)
步骤f、将式(7)得到的新的制冷设备转子压缩机载荷激励F1重复步骤c-步骤d，利用迭代的方法消除误差，直到满足式(6)的要求，获得时变工况下制冷设备转子压缩机的载荷激励。
时变工况下制冷设备转子压缩机载荷激励辨识方法
技术领域
背景技术
发明内容
具体实施方式
| 标题 | 发布/更新时间 | 阅读量 |
|---|---|---|
| 一种冰箱瞬冻储存控制方法及冰箱 | 2020-07-21 | 1 |
| 用于控制可变几何结构压缩机的方法和系统 | 2020-10-02 | 0 |
| 火焰传播管、燃烧器及包括其的燃气轮机 | 2020-09-11 | 1 |
| 一种测量任意角度下非球形颗粒曳力的设备和方法 | 2021-02-09 | 0 |
| 压缩空气储能发电装置 | 2021-10-21 | 0 |
| 一种瞬冻储存的控制方法和冰箱 | 2020-07-22 | 1 |
| 空气调节器 | 2020-06-19 | 2 |
| 一种分子筛制氧机及其分子筛装置 | 2021-01-22 | 1 |
| 一种水加热冷暖一体空调系统 | 2022-07-30 | 1 |
| 一种带水泵系统顶置式新能源H型空调机组 | 2020-07-19 | 1 |
高效检索全球专利专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台，是提供专利分析，专利查询，专利检索等数据服务功能的知识产权数据服务商。
我们的产品包含105个国家的1.26亿组数据，免费查、免费专利分析。
专利汇分析报告产品可以对行业情报数据进行梳理分析，涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度，系统通过AI智能系统对图表进行解读，只需1分钟，一键生成行业专利分析报告。
