CN108896258B - 一种变频压缩机振动载荷计算方法 - Google Patents
一种变频压缩机振动载荷计算方法 Download PDFInfo
- Publication number
- CN108896258B CN108896258B CN201810290263.9A CN201810290263A CN108896258B CN 108896258 B CN108896258 B CN 108896258B CN 201810290263 A CN201810290263 A CN 201810290263A CN 108896258 B CN108896258 B CN 108896258B
- Authority
- CN
- China
- Prior art keywords
- compressor
- vibration
- point
- frequency
- exhaust port
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M7/00—Vibration-testing of structures; Shock-testing of structures
- G01M7/02—Vibration-testing by means of a shake table
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M99/00—Subject matter not provided for in other groups of this subclass
- G01M99/005—Testing of complete machines, e.g. washing-machines or mobile phones
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/06—Power analysis or power optimisation
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Control Of Positive-Displacement Pumps (AREA)
Abstract
本发明公开了一种变频压缩机振动载荷计算方法,包括如下步骤:S1:通过振动测试***对变频压缩机进行振动试验测试;S2:以压缩机吸排气口的连线为x方向,在水平面内垂直于所述吸排气管口连线的方向为y方向,竖直方向为z方向,分别测试压缩机吸气口A点处与排气口B点处在各频率点下的x和y向振动位移;S3:利用有限元仿真技术计算压缩机载荷与吸气口、排气口振动位移之间的映射关系;S4:根据压缩机振动试验测试获得的吸气口A点与排气口B点在各频率点下的x和y方向振动位移,得到压缩机振动载荷。本发明可准确评估压缩机管路振动可靠性,计算简便,能够减少实验次数,降低研发成本,提高生产效率。
Description
技术领域
本发明属于变频压缩机技术领域,具体涉及一种变频压缩机振动载荷计算方法。
背景技术
压缩机结构主要包括电机转子和定子、汽缸、滚动活塞、轴承、储液器及压缩机壳体等部件,影响压缩机振动的主要激励载荷有电机驱动力、气体阻力、不平衡惯性力、电磁转矩负载、压力脉动等,其中不平衡惯性力是由压缩机转子不平衡结构在高速旋转时产生的惯性力,压力脉动是由压缩机周期性吸气和排气引起的脉动压力,压缩机振动主要是由驱动力矩与转子阻力矩的不平衡导致速度波动引起转子***振动,转子***的振动通过各部件间的传递,表现为压缩机壳体的振动,进而通过压缩机吸、排气管口将振动传递给管路,以及通过压缩机底板直接传递给室外机。由于压缩机为封闭的内部结构,其载荷识别难度高,技术复杂,导致压缩机管路***有限元仿真时输入载荷不确定,其结果难以与试验对标,本领域的技术人员通常给定一个估计的压缩机载荷,仅用于方案对比分析,无法定量评估管路振动可靠性。
现有技术中公开了一种制冷设备用转子压缩机载荷测试分析方法(公开号为CN102562568A),然而其仅适用于单个频率点的载荷计算,且计算效率低,准确性较差,无法描述变频压缩机各频率点的振动载荷与吸排气口振动位移的数学关系,难以适应企业大批量生产的需求。
发明内容
针对现有技术中所存在的不足,本发明提供了一种可用于准确评估压缩机管路振动可靠性,减少实验次数,降低研发成本,提高效率的变频压缩机振动载荷计算方法。
一种变频压缩机振动载荷计算方法,包括如下步骤:
S1:对压缩机单体进行振动试验测试,所述压缩机单体包括储液罐和压缩机刚体,所述压缩机刚体顶部设有排气口,所述储液罐顶部设有吸气口;
S2:以压缩机吸气口和排气口的连线为x方向,在水平面内垂直于所述吸气口和排气口连线的方向为y方向,竖直方向为z方向,分别测试压缩机吸气口A点处与排气口B点处在各频率点fi下的x和y向振动位移并分别记为A(xa1i,ya1i)和B(xb1i,yb1i);
S3:利用有限元仿真技术计算压缩机载荷与吸气口、排气口振动位移之间的映射关系:
S31:在压缩机载荷作用点施加z方向的转动力矩Mi,计算压缩机频率f1在不同载荷力矩下的振动位移A(xa2i,ya2i)和B(xb2i,yb2i),通过数据拟合得到频率f1下,力矩M与吸排气口振动位移之间的映射关系,M=f(xa2,ya2,xb2,yb2);
S32:利用有限元仿真技术计算压缩机振动频率响应,分别得到A、B两点在x和y方向的频率响应曲线,通过数据拟合得到振幅与频率的关系,A点x方向记为Axi=gx(fi),A点y方向记为Ayi=gy(fi),B点x方向记为Bxi=hx(fi),B点y方向记为Byi=hy(fi);
S33:根据步骤S32得到各频率点在x和y方向振动位移幅值,如频率f1的振动幅值为Ax1=gx(f1),Ay1=gy(f1),Bx1=hx(f1),By1=hy(f1),则记系数A1i=gx(fi)/gx(f1),A2i=gy(fi)/gy(f1),B1i=hx(fi)/hx(f1),B2i=hy(fi)/hy(f1),那么频率点fi的吸排气口振动位移与力矩载荷关系可表示为:
Mi=f(A1i·xa2,A2i·ya2,B1i·xb2,B2i·yb2) ①;
S4:根据压缩机试验测试获得的吸气口A点与排气口B点在各频率点fi下的x和y方向振动位移A(xa1i,ya1i)和B(xb1i,yb1i),由公式①可得压缩机振动载荷为:
Mi=f(A1i·xa1i,A2i·ya1i,B1i·xb1i,B2i·yb1i)。
相比于现有技术,本发明具有如下有益效果:
本发明通过一种数值计算的方法获得变频压缩机各频率点的振动载荷与吸排气口振动位移的数学关系,再输入实验测试数据可得到较为准确的压缩机载荷值,本发明方法采用的算法计算准确性高,能快速计算变频压缩机全频率段的振动载荷,用于定量评估管路振动可靠性,提高有限元技术的准确性,明显减少实验测试的工作量。
附图说明
图1为变频压缩机结构示意图;
图2为变频压缩机结构俯视图。
其中,1储液罐,2压缩机刚体,3吸气口,4排气口,A和B为振动测试点。
具体实施方式
为了使发明实现的技术手段、创作特征、达成目的与功效易于明白了解,下面结合具体图示,进一步阐述本发明。
一种变频压缩机振动载荷计算方法,包括如下步骤:
S1:对压缩机单体进行振动试验测试,如图1~2所示,所述压缩机单体包括储液罐1和压缩机刚体2,所述压缩机刚体2顶部设有排气口4,所述储液罐1顶部设有吸气口3。
S2:以压缩机吸气口3和排气口4的连线为x方向,在水平面内垂直于吸、排气管口连线为y方向,竖直方向为z方向,分别测试压缩机吸气口3A点处与排气口4B点处在各频率点fi下的x和y向振动位移并分别记为A(xa1i,ya1i)和B(xb1i,yb1i),z方向振动较小可忽略。
S3:利用有限元仿真技术计算压缩机载荷与吸气口3、排气口4振动位移之间的映射关系,包括如下步骤:
S31:在压缩机载荷作用点施加z方向的转动力矩Mi,计算压缩机频率f1在不同载荷力矩下的振动位移A(xa2i,ya2i)和B(xb2i,yb2i),通过数据拟合得到频率f1下,力矩M与吸气口3、排气口4振动位移之间的映射关系,M=f(xa2,ya2,xb2,yb2)。
S32:利用有限元仿真技术计算压缩机振动频率响应,分别得到A、B两点在x和y方向的频率响应曲线,通过数据拟合得到振幅与频率的关系,A点x方向记为Axi=gx(fi),A点y方向记为Ayi=gy(fi),B点x方向记为Bxi=hx(fi),B点y方向记为Byi=hy(fi)。
S33:根据步骤S32可得到各频率点在x和y方向振动位移幅值,如频率f1的振动幅值为Ax1=gx(f1),Ay1=gy(f1),Bx1=hx(f1),By1=hy(f1),则记系数A1i=gx(fi)/gx(f1),A2i=gy(fi)/gy(f1),B1i=hx(fi)/hx(f1),B2i=hy(fi)/hy(f1),那么频率点fi的吸气口3、排气口4振动位移与力矩载荷关系可表示为:
Mi=f(A1i·xa2,A2i·ya2,B1i·xb2,B2i·yb2) ①。
S4:根据压缩机试验测试获得的吸气管口A与排气管口B在各频率点fi下的x和y方向振动位移A(xa1i,ya1i)和B(xb1i,yb1i),由公式①可得压缩机振动载荷为:
Mi=f(A1i·xa1i,A2i·ya1i,B1i·xb1i,B2i·yb1i)。
在现有技术中已存在用于对变频压缩机进行振动试验测试的振动测试***。具体来说,在本申请中,除压缩机单体外,振动测试***可以采用制冷设备、变频控制器、振动传感器以及振动信号采集仪等设备。所述制冷设备用于控制压缩机单体吸排气的压力工况,所述变频控制器用于驱动压缩机逐点扫频运转,所述振动传感器用于检测测试点处(本申请中为A、B两点)的振动信号(比如加速度或振动位移等),所述振动信号采集仪与振动传感器连接,用于将通过振动传感器采集的振动信号汇集起来进行模数转换后再进行分析。
所述fi为变频压缩机工作的一系列频率,频率点f1是fi中任意一个确定的频率点。
本发明中的变频压缩机振动载荷计算方法,先通过测试压缩机吸、排气口X方向和Y方向的振动信号,再采用有限元技术得到某个频率点力矩M与吸排气口振动位移之间的映射关系以及结构的频响函数,然后得出变频压缩机各频率点的振动载荷与吸排气口振动位移的数学关系,最后通过试验数据求得压缩机实际的振动载荷。
这样一来,本发明方法采用的算法计算准确性高,能快速计算变频压缩机全频率段的振动载荷,可实现定量评估管路振动可靠性,提高了有限元技术的准确性,明显减少了实验测试的工作量。
以上所述仅为本发明的优选实施方式,本发明的保护范围并不仅限于上述实施方式,凡是属于本发明原理的技术方案均属于本发明的保护范围。对于本领域的技术人员而言,在不脱离本发明的原理的前提下进行的若干改进,这些改进也应视为本发明的保护范围。
Claims (1)
1.一种变频压缩机振动载荷计算方法,其特征在于包括如下步骤:
S1:对压缩机单体进行振动试验测试,所述压缩机单体包括储液罐和压缩机刚体,所述压缩机刚体顶部设有排气口,所述储液罐顶部设有吸气口;
S2:以压缩机吸气口和排气口的连线为x方向,在水平面内垂直于所述吸气口和排气管口连线的方向为y方向,竖直方向为z方向,分别测试压缩机吸气口A点处与排气口B点处在各频率点fi下的x和y向振动位移并分别记为A(xa1i,ya1i)和B(xb1i,yb1i),所述fi为变频压缩机工作的一系列频率;
S3:利用有限元仿真技术计算压缩机载荷与吸气口、排气口振动位移之间的映射关系:
S31:在压缩机载荷作用点施加z方向的转动力矩Mi,计算压缩机频率f1在不同载荷力矩下的振动位移A(xa2i,ya2i)和B(xb2i,yb2i),通过数据拟合得到频率f1下,力矩M与吸排气口振动位移之间的映射关系,M=f(xa2,ya2,xb2,yb2),所述频率点f1是fi中任意一个确定的频率点;
S32:利用有限元仿真技术计算压缩机振动频率响应,分别得到A、B两点在x和y方向的频率响应曲线,通过数据拟合得到振幅与频率的关系,A点x方向记为Axi=gx(fi),A点y方向记为Ayi=gy(fi),B点x方向记为Bxi=hx(fi),B点y方向记为Byi=hy(fi);
S33:根据步骤S32得到各频率点在x和y方向振动位移幅值,所述频率f1的振动幅值为Ax1=gx(f1),Ay1=gy(f1),Bx1=hx(f1),By1=hy(f1),则记系数A1i=gx(fi)/gx(f1),A2i=gy(fi)/gy(f1),B1i=hx(fi)/hx(f1),B2i=hy(fi)/hy(f1),那么频率点fi的吸排气口振动位移与力矩载荷关系可表示为:
Mi=f(A1i·xa2,A2i·ya2,B1i·xb2,B2i·yb2) ①;
S4:根据压缩机试验测试获得的吸气口A点与排气口B点在各频率点fi下的x和y方向振动位移A(xa1i,ya1i)和B(xb1i,yb1i),由公式①可得压缩机振动载荷为:
Mi=f(A1i·xa1i,A2i·ya1i,B1i·xb1i,B2i·yb1i)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810290263.9A CN108896258B (zh) | 2018-03-30 | 2018-03-30 | 一种变频压缩机振动载荷计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810290263.9A CN108896258B (zh) | 2018-03-30 | 2018-03-30 | 一种变频压缩机振动载荷计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108896258A CN108896258A (zh) | 2018-11-27 |
CN108896258B true CN108896258B (zh) | 2020-05-26 |
Family
ID=64342272
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810290263.9A Active CN108896258B (zh) | 2018-03-30 | 2018-03-30 | 一种变频压缩机振动载荷计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108896258B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111651925B (zh) * | 2020-06-05 | 2022-09-16 | 四川长虹空调有限公司 | 空调管路低频噪声声源识别方法及*** |
CN113239489B (zh) * | 2021-05-12 | 2022-12-09 | 四川长虹空调有限公司 | 变频单转子压缩机振动的仿真评价方法 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB0019434D0 (en) * | 2000-08-09 | 2000-09-27 | Rolls Royce Plc | A device and method for fatigue testing of materials |
KR101837984B1 (ko) * | 2011-10-06 | 2018-03-15 | 주식회사 대유위니아 | 압축기의 상선 오조립 검출방법 |
CN102562568B (zh) * | 2012-02-14 | 2014-08-06 | 合肥工业大学 | 制冷设备用转子压缩机载荷测试分析方法 |
CN104791234B (zh) * | 2015-05-04 | 2016-10-26 | 合肥工业大学 | 制冷设备转子压缩机启动工况下载荷激励测试分析方法 |
CN106529048B (zh) * | 2016-11-11 | 2019-04-09 | 四川长虹空调有限公司 | 空调压缩机配管振动仿真方法 |
CN106815446B (zh) * | 2017-01-24 | 2019-12-27 | 合肥工业大学 | 时变工况下制冷设备转子压缩机载荷激励辨识方法 |
CN107808056A (zh) * | 2017-10-30 | 2018-03-16 | 四川长虹空调有限公司 | 空调压缩机配管振动可靠性评估方法 |
-
2018
- 2018-03-30 CN CN201810290263.9A patent/CN108896258B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN108896258A (zh) | 2018-11-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Elhaj et al. | Numerical simulation and experimental study of a two-stage reciprocating compressor for condition monitoring | |
CN107808056A (zh) | 空调压缩机配管振动可靠性评估方法 | |
CN108896258B (zh) | 一种变频压缩机振动载荷计算方法 | |
CN108920804B (zh) | 制冷设备变频压缩机激励载荷仿真计算方法 | |
CN207906130U (zh) | 一种轴流风机测试*** | |
CN107167223A (zh) | 变频空调器压缩机与管路振动测试***及方法 | |
CN111413404A (zh) | 基于叶尖定时和支持向量机原理的叶片裂纹在线测量方法 | |
CN103742277B (zh) | 一种发动机摩擦扭矩计算方法及装置 | |
CN106762592A (zh) | 一种舰船用直驱式柱塞泵谐振频率测试方法 | |
CN205384141U (zh) | 工程机械变速箱总成性能检测分析试验台 | |
CN102128674A (zh) | 测量水轮机过流部件卡门涡振动的方法 | |
CN107748033A (zh) | 一种筒型径向力传感器标定装置及其标定方法 | |
CN106815446B (zh) | 时变工况下制冷设备转子压缩机载荷激励辨识方法 | |
CN104005975A (zh) | 一种轴流式通风机失速和喘振的诊断方法 | |
CN117028299B (zh) | 一种磁悬浮分子泵涡轮蠕胀测试工装 | |
CN103048153A (zh) | 一种用于圆筒形结构声学性能测试的试验*** | |
CN110702314B (zh) | 基于柔性转子加速响应信息的无试重瞬态高速动平衡方法 | |
CN203441721U (zh) | 一种压缩机试验装置 | |
Cui et al. | Instantaneous frequency estimation-based order tracking for bearing fault diagnosis under strong noise | |
CN113720369B (zh) | 仪器仪表加速失效激发智能运维试验平台 | |
CN207960910U (zh) | 空调器压缩机振动测试*** | |
KR20150119873A (ko) | 출력 결정 방법 및 터보 기계 | |
CN113139317B (zh) | 空调压缩机配管振动仿真方法 | |
Han et al. | Recursive engine in-cylinder pressure estimation merely using structural vibration signal | |
CN103234605B (zh) | 汽车流量计双向气流及反应速度的对比测试方法和测试仪 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |