CN115048789A - 涡旋压缩机管路振动数值仿真计算方法 - Google Patents
涡旋压缩机管路振动数值仿真计算方法 Download PDFInfo
- Publication number
- CN115048789A CN115048789A CN202210670272.7A CN202210670272A CN115048789A CN 115048789 A CN115048789 A CN 115048789A CN 202210670272 A CN202210670272 A CN 202210670272A CN 115048789 A CN115048789 A CN 115048789A
- Authority
- CN
- China
- Prior art keywords
- compressor
- pipeline
- vibration
- load
- excitation load
- 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.)
- Granted
Links
Images
Classifications
-
- 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
-
- 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/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Geometry (AREA)
- Evolutionary Computation (AREA)
- Fluid Mechanics (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及压缩机管路振动领域,具体涉及一种涡旋压缩机管路振动数值仿真计算方法,能够更贴近实际的模拟压缩机管路振动情况,极大地缩短了仿真计算周期,同时能够对管路涉及方案进行优化。本发明涡旋压缩机管路振动数值仿真计算方法,包括:建立涡旋压缩机数值仿真计算模型,计算模型包括涡旋压缩机本体以及补气管路;对计算模型进行简化处理,添加约束条件,进行振动响应仿真数值计算,所述约束条件包括压缩机机械激励载荷及流体脉动激励载荷。本发明适用于涡旋压缩机管路振动数值仿真的计算。
Description
技术领域
本发明涉及压缩机管路振动领域,具体涉及一种涡旋压缩机管路振动数值仿真计算方法。
背景技术
涡旋压缩机普遍应用于大型中央空调系统,相比转子式压缩机,涡旋压缩机运行平稳、振动小、噪音小。一般为了在低温环境中实现较大的压缩比,涡旋压缩机采用了喷气增焓技术,在涡旋压缩机涡旋盘上开设补气孔,增加一个中间压力的补气回路,增加系统制冷剂质量流量,提升涡旋压缩机的压缩效率和能力。涡旋压缩机激励载荷复杂,涡盘运动过程中涡盘受压缩腔气体力作用产生周向力,同时受离心力和竖直方向的垂向力作用,对于补气管路来说,除受涡旋压缩机产生的机械振动载荷外,还受到涡旋压缩机喷气产生的周期性压力脉动,涡旋压缩机运行频率范围宽,一般运行频率越高振动越大。大型中央空调涡旋压缩机管路系统结构较为复杂,四通阀组件包含吸气管路、排气管路、补气管路等连接管,管路系统固有频率多,容易与涡旋压缩机激励载荷频率重合,形成共振及流固耦合振动。管路振动会引起管路疲劳断裂,导致系统失效,一般企业里通过振动应力试验测试管路应力,根据管路振动应力值评估管路设计方案是否可行。这种仅依靠试验测试的评估手段周期长、成本高、效率低,且无法在概念设计阶段事先对方案优良进行预判。
发明内容
本发明的目的是提供一种涡旋压缩机管路振动数值仿真计算方法,能够更贴近实际的模拟压缩机管路振动情况,极大地缩短了仿真计算周期,同时能够对管路涉及方案进行优化。
本发明采取如下技术方案实现上述目的,涡旋压缩机管路振动数值仿真计算方法,包括:
建立涡旋压缩机数值仿真计算模型,计算模型包括涡旋压缩机本体以及补气管路;
对计算模型进行简化处理,添加约束条件,并进行振动响应仿真数值计算,所述约束条件包括压缩机机械激励载荷及流体脉动激励载荷。
进一步的是,为了提高机械激励载荷计算的准确性,压缩机机械激励载荷确定的具体方法包括:
压缩机机械激励载荷F=[MZ DX DY DZ],MZ为压缩机本体竖直方向力矩载荷,DX为压缩机本体在X方向的位移载荷,DY为压缩机本体在Y方向的位移载荷,DZ为压缩机本体在Z方向的位移载荷。
位移载荷与压缩机运行频率关系为:
DX=KX×f,DY=KY×f,DZ=KZ×f,其中KX、KY、KZ为常数,f为压缩机运行频率,三个方向的位移载荷的相位相差n度,n≥90度。
进一步的是,为了提高流体脉动激励载荷的计算的准确性,流体脉动激励载荷确定的具体方法包括:
通过压力脉动试验直接测试补气管路内部流体压力脉动数据,通过公式拟合流体脉动激励载荷P,P=A sin(wt+φ),A为压力脉动载荷幅值,φ为相位。
流体脉动激励载荷加载的具体方法包括:
在补气管路表面加载流体脉动激励载荷P,载荷作用方向为流体流动方向,加载位置为补气管路各个弯位处。
进一步的是,压缩机机械激励载荷F优化的具体方法包括:
通过试验测试数据与仿真数据对压缩机机械激励载荷F进行优化,试验测试数据选择振动加速度、振动位移、振动速度、振动应力以及振动应变中任意一种测量值,测量点至少包括压缩机在圆周方向上相互垂直的两个点,其余点为压缩机表面其他位置的测量点;仿真数据中以相同位置的测量点为监测点,加载初始压缩机机械激励载荷F以及补气管路流体脉动激励载荷P,从仿真结果中获取与监测点相同的测量值,以载荷F中MZ、KX、KY、KZ为变量,以监测点的仿真值与实测值的最小差值为目标,通过迭代仿真计算,获得优化的目标值,确定优化的压缩机机械激励载荷F。
进一步的是,为了对管路结构进行优化,仿真计算方法还包括:通过多次迭代计算,获取优化后的压缩机激励载荷数据,基于优化后的压缩机激励载荷数据对管路系统进行数值仿真计算,优化管路结构。
本发明在进行振动响应仿真数值计算时,增加了约束条件:压缩机机械激励载荷及流体脉动激励载荷,这样能更接近实际的模拟涡旋压缩机管路振动特征;同时简化了模型结构,避免了复杂的流固耦合计算,极大地提高了仿真计算效率,缩短了仿真计算周期;并且提供了优化方法,可对管路结构进行优化,提高最终产品的可靠性。
附图说明
图1是本发明实施例提供的仿真优化计算的流程图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述。显然,所描述的实施例是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
下面结合附图对本发明的具体实施方式进行详细说明。
本发明涡旋压缩机管路振动数值仿真计算方法,包括:
建立涡旋压缩机数值仿真计算模型,在本发明的一种实施例中,计算模型包括涡旋压缩机本体、气液分离器、油液分离器、吸气管路、排气管路以及补气管路;
对计算模型进行简化处理,添加约束条件,并进行振动响应仿真数值计算,约束条件包括压缩机机械激励载荷及流体脉动激励载荷。
在本发明的一种实施例中,机械激励载荷主要由压缩机本体机械振动产生,流体脉动激励载荷主要由补气管路内流体产生,因此本发明压缩机机械激励载荷确定的具体方法包括:
压缩机机械激励载荷F=[MZ DX DY DZ],MZ为压缩机本体竖直方向力矩载荷,DX为压缩机本体在X方向的位移载荷,DY为压缩机本体在Y方向的位移载荷,DZ为压缩机本体在Z方向的位移载荷。
位移载荷与压缩机运行频率关系为:
DX=KX×f,DY=KY×f,DZ=KZ×f,其中KX、KY、KZ为常数,f为压缩机运行频率,三个方向的位移载荷的相位相差n度,n≥90度。提高了机械激励载荷计算的准确性。
流体脉动激励载荷确定的具体方法包括:
通过压力脉动试验直接测试补气管路内部流体压力脉动数据,通过公式拟合流体脉动激励载荷P,P=A sin(wt+φ),A为压力脉动载荷幅值,φ为相位。提高了流体脉动激励载荷的计算的准确性。
流体脉动激励载荷加载的具体方法包括:
在补气管路表面加载流体脉动激励载荷P,载荷作用方向为流体流动方向,加载位置为补气管路各个弯位处。
本发明压缩机机械激励载荷F优化的具体方法包括:
通过试验测试数据与仿真数据对压缩机机械激励载荷F进行优化,试验测试数据选择振动加速度、振动位移、振动速度、振动应力以及振动应变中任意一种测量值,测量点至少包括压缩机在圆周方向上相互垂直的两个点,其余点为压缩机表面其他位置的测量点;仿真数据中以相同位置的测量点为监测点,加载初始压缩机机械激励载荷F以及补气管路流体脉动激励载荷P,从仿真结果中获取与监测点相同的测量值,以载荷F中MZ、KX、KY、KZ为变量,以监测点的仿真值与实测值的最小差值为目标,通过迭代仿真计算,获得优化的目标值,确定优化的压缩机机械激励载荷F。
仿真计算方法还包括:通过多次迭代计算,获取优化后的压缩机激励载荷数据,基于优化后的压缩机激励载荷数据对管路系统进行数值仿真计算,优化管路结构。极大地提高了对管路的优化能力。
其中,优化管路结构具体包括对涡旋压缩机本体、气液分离器、油液分离器、吸气管路、排气管路以及补气管路整个的管路结构进行优化。
本发明实施例提供的仿真优化计算的流程图如图1所示,建立涡旋压缩机数值仿真计算模型,然后加载仿真载荷初始值,具体可以加载初始压缩机机械激励载荷F以及补气管路流体脉动激励载荷P,载荷初始值加载完成后进行仿真计算,判断仿真结果是否满足误差要求,若满足要求则进行仿真优化计算,具体通过多次迭代计算,获取优化后的压缩机激励载荷数据,基于优化后的压缩机激励载荷数据对管路系统进行数值仿真计算,优化管路结构。若不满足要求,则优化载荷后,再次进行加载仿真。
综上所述,本发明能更接近实际的模拟涡旋压缩机管路振动特征;同时简化了模型结构,避免了复杂的流固耦合计算,极大地提高了仿真计算效率,缩短了仿真计算周期;并且提供了优化方法,可对管路结构进行优化,提高最终产品的可靠性。
Claims (7)
1.涡旋压缩机管路振动数值仿真计算方法,其特征在于,包括:
建立涡旋压缩机数值仿真计算模型,计算模型包括涡旋压缩机本体以及补气管路;
对计算模型进行简化处理,添加约束条件,并进行振动响应仿真数值计算,所述约束条件包括压缩机机械激励载荷及流体脉动激励载荷。
2.根据权利要求1所述的涡旋压缩机管路振动数值仿真计算方法,其特征在于,压缩机机械激励载荷确定的具体方法包括:
压缩机机械激励载荷F=[MZ DX DY DZ],MZ为压缩机本体竖直方向力矩载荷,DX为压缩机本体在X方向的位移载荷,DY为压缩机本体在Y方向的位移载荷,DZ为压缩机本体在Z方向的位移载荷。
3.根据权利要求2所述的涡旋压缩机管路振动数值仿真计算方法,其特征在于,位移载荷与压缩机运行频率关系为:
DX=KX×f,DY=KY×f,DZ=KZ×f,其中KX、KY、KZ为常数,f为压缩机运行频率,三个方向的位移载荷的相位相差n度,n≥90度。
4.根据权利要求2所述的涡旋压缩机管路振动数值仿真计算方法,其特征在于,流体脉动激励载荷确定的具体方法包括:
通过压力脉动试验直接测试补气管路内部流体压力脉动数据,通过公式拟合流体脉动激励载荷P,P=Asin(wt+φ),A为压力脉动载荷幅值,φ为相位。
5.根据权利要求4所述的涡旋压缩机管路振动数值仿真计算方法,其特征在于,流体脉动激励载荷加载的具体方法包括:
在补气管路表面加载流体脉动激励载荷P,载荷作用方向为流体流动方向,加载位置为补气管路各个弯位处。
6.根据权利要求4所述的涡旋压缩机管路振动数值仿真计算方法,其特征在于,压缩机机械激励载荷F优化的具体方法包括:
通过试验测试数据与仿真数据对压缩机机械激励载荷F进行优化,试验测试数据选择振动加速度、振动位移、振动速度、振动应力以及振动应变中任意一种测量值,测量点至少包括压缩机在圆周方向上相互垂直的两个点,其余点为压缩机表面其他位置的测量点;仿真数据中以相同位置的测量点为监测点,加载初始压缩机机械激励载荷F以及补气管路流体脉动激励载荷P,从仿真结果中获取与监测点相同的测量值,以载荷F中MZ、KX、KY、KZ为变量,以监测点的仿真值与实测值的最小差值为目标,通过迭代仿真计算,获得优化的目标值,确定优化的压缩机机械激励载荷F。
7.根据权利要求6所述的涡旋压缩机管路振动数值仿真计算方法,其特征在于,仿真计算方法还包括:通过多次迭代计算,获取优化后的压缩机激励载荷数据,基于优化后的压缩机激励载荷数据对管路系统进行数值仿真计算,优化管路结构。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210670272.7A CN115048789B (zh) | 2022-06-14 | 2022-06-14 | 涡旋压缩机管路振动数值仿真计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210670272.7A CN115048789B (zh) | 2022-06-14 | 2022-06-14 | 涡旋压缩机管路振动数值仿真计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN115048789A true CN115048789A (zh) | 2022-09-13 |
CN115048789B CN115048789B (zh) | 2025-05-16 |
Family
ID=83161827
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210670272.7A Active CN115048789B (zh) | 2022-06-14 | 2022-06-14 | 涡旋压缩机管路振动数值仿真计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115048789B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115982863A (zh) * | 2023-03-16 | 2023-04-18 | 西安航天动力研究所 | 一种控制响应的冲击试验仿真方法、装置、设备及介质 |
CN119720706A (zh) * | 2025-03-03 | 2025-03-28 | 浙江三田汽车空调压缩机有限公司 | 一种基于有限元模型的涡旋盘应力应变分析方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108896258A (zh) * | 2018-03-30 | 2018-11-27 | 四川长虹空调有限公司 | 一种变频压缩机振动载荷计算方法 |
CN109063312A (zh) * | 2018-07-26 | 2018-12-21 | 四川长虹空调有限公司 | 变频空调器双转子压缩机配管系统振动仿真方法 |
CN110287569A (zh) * | 2019-06-18 | 2019-09-27 | 珠海格力电器股份有限公司 | 空调管路设计的分析方法和装置 |
-
2022
- 2022-06-14 CN CN202210670272.7A patent/CN115048789B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108896258A (zh) * | 2018-03-30 | 2018-11-27 | 四川长虹空调有限公司 | 一种变频压缩机振动载荷计算方法 |
CN109063312A (zh) * | 2018-07-26 | 2018-12-21 | 四川长虹空调有限公司 | 变频空调器双转子压缩机配管系统振动仿真方法 |
CN110287569A (zh) * | 2019-06-18 | 2019-09-27 | 珠海格力电器股份有限公司 | 空调管路设计的分析方法和装置 |
Non-Patent Citations (5)
Title |
---|
余洋: "涡旋压缩机动力特性及仿真模拟研究", 中国优秀博士学位论文全文数据库(电子期刊), 15 September 2014 (2014-09-15), pages 029 - 14 * |
秦振振 等: "基于虚拟迭代的空调管路疲劳寿命预估方法", 合肥工业大学学报(自然科学版), no. 05, 28 May 2019 (2019-05-28), pages 595 - 600 * |
谷欢欢 等: "涡旋压缩机激励载荷反求分析", 噪声与振动控制, vol. 41, no. 01, 18 February 2021 (2021-02-18), pages 1 - 2 * |
邓培生 等: "基于变频涡旋压缩机的管路振动仿真技术研究及应用", 2023年中国家用电器技术大会论文集, 18 May 2023 (2023-05-18), pages 709 - 713 * |
马堃: "空调室外机载流管道振动和应力分析研究", 中国优秀硕士学位论文全文数据库(电子期刊), 15 February 2017 (2017-02-15), pages 2 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115982863A (zh) * | 2023-03-16 | 2023-04-18 | 西安航天动力研究所 | 一种控制响应的冲击试验仿真方法、装置、设备及介质 |
CN119720706A (zh) * | 2025-03-03 | 2025-03-28 | 浙江三田汽车空调压缩机有限公司 | 一种基于有限元模型的涡旋盘应力应变分析方法 |
CN119720706B (zh) * | 2025-03-03 | 2025-05-30 | 浙江三田汽车空调压缩机有限公司 | 一种基于有限元模型的涡旋盘应力应变分析方法 |
Also Published As
Publication number | Publication date |
---|---|
CN115048789B (zh) | 2025-05-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109779894B (zh) | 一种基于神经网络算法的往复式压缩机故障诊断系统及方法 | |
CN115048789A (zh) | 涡旋压缩机管路振动数值仿真计算方法 | |
CN109063312B (zh) | 变频空调器双转子压缩机配管系统振动仿真方法 | |
CN101080700A (zh) | 用于真空泵的故障保护和预知维护的精确诊断方法及精确诊断系统 | |
CN114912371B (zh) | 一种贯流泵运行稳定检测的数字孪生体建模系统及方法 | |
CN108920804B (zh) | 制冷设备变频压缩机激励载荷仿真计算方法 | |
CN109063246A (zh) | 变频压缩机配管振动评价方法 | |
CN113688510B (zh) | 一种燃料电池车用离心压缩系统喘振预测方法 | |
CN114818205B (zh) | 一种航空发动机全生命周期叶尖间隙在线感知方法 | |
Manwaring et al. | Structures and Dynamics Committee Best Paper of 1996 Award: Inlet Distortion Generated Forced Response of a Low-Aspect-Ratio Transonic Fan | |
CN115859536B (zh) | 压气机转子叶片非同步振动锁频数值模拟方法 | |
CN114330083A (zh) | 一种轴承参数不确定的离心泵转子动力学稳健预测方法 | |
CN115438548A (zh) | 一种飞机液压管路系统有限元模型修正方法 | |
CN114117862A (zh) | 一种发动机整机振动测点选取方法和系统 | |
CN115435894B (zh) | 基于模拟退火算法的叶尖定时振动应力反演方法 | |
Lerche et al. | Dynamic stress prediction in centrifugal compressor blades using fluid structure interaction | |
CN116754134A (zh) | 基于试验与仿真数据融合的转子不平衡状态识别方法 | |
CN113218049B (zh) | 一种变频空调器压缩机与管路快速匹配方法 | |
Gan et al. | Numerical examination of lock-in hypothesis of non-synchronous vibration in an axial compressor | |
de Cherisey et al. | Optimization of a turbomachinery blade with regards to tip-rub events | |
Seinturier et al. | Turbine Mistuned Forced Response Prediction: Comparison with Experimental Results | |
CN111222277A (zh) | 输气站增压泵进出口管道的振动评价方法 | |
CN117405396A (zh) | 轴承座参振动力学参数附加质量试验辨识方法和装置 | |
Harrison et al. | From CAD to 1D: a direct approach to modeling scroll compressors with multi-physics simulation | |
CN112781722B (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 |