CN106599422A - 叶片泵转子系统的振动仿真分析方法及其装置 - Google Patents
叶片泵转子系统的振动仿真分析方法及其装置 Download PDFInfo
- Publication number
- CN106599422A CN106599422A CN201611099497.2A CN201611099497A CN106599422A CN 106599422 A CN106599422 A CN 106599422A CN 201611099497 A CN201611099497 A CN 201611099497A CN 106599422 A CN106599422 A CN 106599422A
- Authority
- CN
- China
- Prior art keywords
- impeller
- transient
- water body
- stress
- transient state
- 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
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
- 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
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
Abstract
本发明提供一种叶片泵转子系统的振动仿真分析方法及装置。该方法包括如下步骤:利用所述三维水体网格建立计算流体力学模型;以稳态计算结果为初始值,根据预设参数获取每个叶片网格单元和每个轮毂网格单元瞬态计算结果,并通过积分的方法获得叶轮的瞬态受力;以叶轮的瞬态受力作为叶轮质点的载荷条件,采用有限元分析根据预设的转子系统信息对由主轴、轴承和叶轮质点组成的转子系统进行瞬态动力学分析,并获得瞬态动力学分析结果;根据瞬态动力学分析结果获得转子系统的振动情况。该方法将计算流体力学与转子动力学有限元分析两大模型结合起来,可以快速、准确地获得水力诱导下叶片泵转子系统的振动情况,满足叶片泵子系统设计需求。
Description
【技术领域】
本发明涉及一种叶片泵仿真分析技术领域,具体地,涉及一种叶片泵转子系统的振动仿真分析方法及装置。
【背景技术】
叶片泵是应用最为广泛的流体输送管网系统核心装备和能量转换装置,对流体输送系统的安全稳定运行起到决定性作用。一般情况下,叶片泵的运行状况与其结构振动密切相关,尤其是水力作用下转子系统的振动特征。如何在各种可能的运行工况下,控制转系系统的振动烈度,并合理地安排装配工艺,是叶片泵设计制造最重要的内容之一。目前主要针对设计方案试制叶片泵样机,开展水力性能试验并使用振动传感采集设备获得轴承等转子系统部件的振动情况。尽管计算流体力学技术已经被广泛地用于叶片泵领域的水力性能预测,各种结构分析有限元技术也普及至叶片泵领域的结构力学仿真领域,也有将计算流体力学与有限元分析耦合进行叶轮强度校核的应用案例;但是,对于如何采用仿真分析的方法,准确地获得流体诱导下转子部件的振动信息,尚没有公知的信息披露和现有的应用案例。
【发明内容】
本发明提供一种叶片泵转子系统的振动仿真分析方法及其装置,通过将计算流体力学与转子系统动力学有限元分析结合起来,快速、准确获得不同运行工况点下叶片泵转子系统的振动情况。
为实现上述目的,本发明提供一种叶片泵转子系统的振动仿真分析方法,包括如下步骤:
对叶片泵的三维水体模型进行网格划分得到三维水体网格,并利用所述三维水体网格建立计算流体力学模型,所述三维水体网格包括多个叶片网格单元和多个轮毂网格单元;
以稳态计算结果为初始值,根据预设参数获取每个叶片网格单元和每个轮毂网格单元瞬态计算结果,并通过积分的方法获得叶轮的瞬态受力,所述预设参数包括时间步长和第一计算总时长;
以叶轮的瞬态受力作为叶轮质点的载荷条件,采用有限元分析根据预设的转子系统信息对由主轴、轴承和叶轮质点组成的转子系统进行瞬态动力学分析,并获得瞬态动力学分析结果,其中,有限元分析的计算总时长与第一计算总时长相同;
根据瞬态动力学分析结果获得转子系统的振动情况。
优选的,所述对叶片泵的三维水体模型进行网格划分得到三维水体网格,并利用所述三维水体网格建立计算流体力学模型包括:
分别单独绘制叶片泵的叶轮水体、吸水室水体和压水室水体,再将三者进行装配组合为一个整体进行网格划分得到三维水体网格;
利用所述三维水体网格建立计算流体力学模型;
确定计算流体力学模型的进口边界条件和出口边界条件。
优选的,所述确定计算流体力学模型的进口边界条件和出口边界条件包括:
分别建立所述吸水室水体与所述叶轮水体以及所述叶轮水体与所述压水室水体的交界面联系,并以所述吸水室水体与所述叶轮水体的交界面为进口,所述叶轮水体与所述压水室水体的交界面为出口;
根据设计参数的工况点确定计算流体力学模型的进口边界条件和出口边界条件,所述设计参数包括流体介质的密度和粘度以及叶轮水体的旋转速度。
优选的,所述以稳态计算结果为初始值,根据预设参数获取每个叶片网格单元和每个轮毂网格单元瞬态计算结果,并通过积分的方法获得叶轮的瞬态受力,所述预设参数包括时间步长和第一计算总时长包括:
根据计算流体力学模型进行稳态计算,获得稳态计算结果;
以稳态计算结果为初始值,设定计算时间步长和第一计算总时长,根据计算流体力学模型计算每个叶片网格单元和每个轮毂网格单元表面水压分布,获得每个叶片网格单元水平面和竖直面轴向瞬态受力和每个轮毂网格单元水平面和竖直面轴向瞬态受力;
针对每个时间步长,分别对所述叶片网格单元和所述轮毂网格单元的水平面以及所述叶片网格单元和所述轮毂网格单元的竖直面轴向瞬态受力进行积分,获得水平面和竖直面上叶轮的瞬态受力随时间的分布图。
优选的,所述以叶轮的瞬态受力作为叶轮质点的载荷条件,采用有限元分析根据预设的转子系统信息对由主轴、轴承和叶轮质点组成的转子系统进行瞬态动力学分析,并获得瞬态动力学分析结果,其中,有限元分析的计算总时长与第一计算总时长相同包括:
对由主轴、轴承和叶轮质点组成的转子系统的三维有限元瞬态分析模型进行有限元网格划分;
分别将水平面上叶轮的瞬态受力随时间的分布和竖直面上叶轮的瞬态受力随时间的分布导入所述三维有限元瞬态分析模型,作为叶轮质点受力的载荷条件;
设定转子系统信息和计算总时长,对所述转子系统进行瞬态动力学分析并获得瞬态动力学分析结果,所述转子系统信息包括主轴转速和轴承约束条件,其计算总时长与所述第一计算总时长相同。
此外,为实现上述目的,本发明同时提供一种叶片泵转子系统的振动仿真分析装置,所述叶片泵转子系统的振动仿真分析装置包括:
计算流体力学模型建模模块,用于对叶片泵的三维水体模型进行网格划分得到三维水体网格并利用所述三维水体网格建立计算流体力学模型,所述三维水体网格包括多个叶片网格单元和多个轮毂网格单元;
叶轮的瞬态受力获取模块,用于以稳态计算结果为初始值,根据预设参数获取每个叶片网格单元和每个轮毂网格单元瞬态计算结果,并通过积分的方法得到叶轮的瞬态受力;
转子系统瞬态动力学分析模块,用于以叶轮的瞬态受力作为叶轮质点的载荷条件,采用有限元分析根据预设的转子系统信息对由主轴、轴承和叶轮质点组成的转子系统进行瞬态动力学分析并获得瞬态动力学分析结果;
仿真分析结果获取模块,用于根据瞬态动力学分析结果获得转子系统的振动情况。
优选的,所述计算流体力学模型建模模块包括:
三维建模单元,用于分别单独绘制叶片泵的叶轮水体、吸水室水体和压水室水体,再将三者进行装配组合为一个整体;
网格划分单元,用于进行网格划分得到三维水体网格;
力学模型建立单元,用于利用所述三维水体网格建立计算流体力学模型;
边界条件确定单元,用于确定计算流体力学模型的进口边界条件和出口边界条件。
优选的,所述边界条件确定单元包括:
交界面建立子单元,用于建立所述吸水室水体与所述叶轮水体以及所述叶轮水体与所述压水室水体的交界面联系,并以所述吸水室水体与所述叶轮水体的交界面为进口,所述叶轮水体与所述压水室水体的交界面为出口;
条件确定子单元,用于根据设计参数的工况点确定计算流体力学模型的进口边界条件和出口边界条件,所述设计参数包括流体介质的密度和粘度以及叶轮水体的旋转速度。
优选的,所述叶轮的瞬态受力获取模块包括:
稳态计算单元,用于根据计算流体力学模型进行稳态计算,获得稳态计算结果;
瞬态计算单元,用于根据计算流体力学模型计算每个叶片网格单元和每个轮毂网格单元表面水压分布,获得每个叶片网格单元水平面和竖直面轴向瞬态受力和每个轮毂网格单元水平面和竖直面轴向瞬态受力;
积分单元,用于针对每个时间步长,分别对所述叶片网格单元和所述轮毂网格单元的水平面以及所述叶片网格单元和所述轮毂网格单元的竖直面轴向瞬态受力进行积分,获得水平面和竖直面上叶轮的瞬态受力随时间的分布。
优选的,所述转子系统瞬态动力学分析模块包括:
有限元网格划分单元,用于对由主轴、轴承和叶轮质点组成的转子系统的三维有限元瞬态分析模型进行有限元网格划分;
叶轮质点受力导入单元,用于将水平面上叶轮的瞬态受力随时间的分布和竖直面上叶轮的瞬态受力随时间的分布导入所述三维有限元瞬态分析模型,作为叶轮质点受力的载荷条件;
瞬态分析单元,用于对所述转子系统进行瞬态动力学分析并获得瞬态动力学分析结果。
相较于现有技术,本发明提供的叶片泵转子系统的振动仿真分析方法及其装置具有以下有益效果:本发明将计算流体力学与转子动力学有限元分析两大模型结合起来,首先通过对叶片泵进行网格划分建立计算流体力学模型,再通过计算流体力学模型得到每个叶片网格单元和每个轮毂网格单元的瞬态计算结果,通过积分的方法可获得一个周期内叶轮的瞬态受力情况,将叶轮的瞬态受力情况作为叶轮质点的载荷条件导入转子动力学有限元分析模型,最后可以获得转子振动系统的振动情况,作为叶片泵设计制造的依据。该方法不但可以快速、准确地获得水力诱导下叶片泵转子系统的振动情况,而且振动仿真分析预测精度高,与实测值的偏差在15%以内,满足叶片泵子系统设计需求。
【附图说明】
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其它的附图,其中:
图1是本发明提供的叶片泵转子系统的振动仿真分析方法第一实施例的步骤流程图;
图2是本发明提供的叶片泵转子系统的振动仿真分析方法第一实施例所使用的叶片泵的三维水体模型示意图;
图3是本发明提供的叶片泵转子系统的振动仿真分析方法第一实施例所得到的叶片泵的三维水体网格示意图;
图4是本发明提供的叶片泵转子系统的振动仿真分析方法第一实施例一个转动周期叶轮的瞬态受力情况示意图;
图5是本发明提供的叶片泵转子系统的振动仿真分析方法第一实施例建立的转子系统有限元模型示意图;
图6是本发明提供的叶片泵转子系统的振动仿真分析方法第一实施例所建立的主轴三维模型示意图;
图7是本发明提供的叶片泵转子系统的振动仿真分析方法第一实施主轴有限元网格示意图;
图8是本发明提供的叶片泵转子系统的振动仿真分析方法第一实施例主轴变形情况示意图;
图9是本发明提供的叶片泵转子系统的振动仿真分析方法第一实施例一个转动周期内叶轮质点处的主轴速度变化示意图;
图10是本发明提供的叶片泵转子系统的振动仿真分析方法第二实施例的步骤流程图;
图11是本发明提供的叶片泵转子系统的振动仿真分析方法第三实施例的步骤流程图;
图12是本发明提供的叶片泵转子系统的振动仿真分析方法第四实施例的步骤流程图;
图13是本发明提供的叶片泵转子系统的振动仿真分析方法第五实施例的步骤流程图;
图14是本发明提供的叶片泵转子系统的振动仿真分析装置第一实施例的功能模块图;
图15是本发明提供的叶片泵转子系统的振动仿真分析装置第二实施例的功能模块图;
图16是本发明提供的叶片泵转子系统的振动仿真分析装置第三实施例的功能模块图;
图17是本发明提供的叶片泵转子系统的振动仿真分析装置第四实施例的功能模块图;
图18是本发明提供的叶片泵转子系统的振动仿真分析装置第五实施例的功能模块图。
【具体实施方式】
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅是本发明的一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其它实施例,都属于本发明保护的范围。
本发明实施例一至六的执行主体可为叶片泵转子系统的振动仿真分析装置,例如,叶片泵转子系统的振动仿真分析设备或安装在计算机上的具有相应功能的模块组件,本发明对此不进行限定。
本发明提供一种叶片泵转子系统的振动仿真分析方法,请参阅图1,为本发明提供的叶片泵转子系统的振动仿真分析方法第一实施例的步骤流程图。在本实施例中,所述叶片泵转子系统的振动仿真分析方法包括如下步骤:
步骤S1、对叶片泵的三维水体模型进行网格划分得到三维水体网格,并利用所述三维水体网格建立计算流体力学模型,所述三维水体网格包括多个叶片网格单元和多个轮毂网格单元;
具体地,应用计算机辅助制图以及数学辅助放算软件绘制三维水体模型并进行网格划分,检查网格质量,然后利用所述三维水体网格建立计算流体力学模型。本实施例的叶片泵选择350S16型单级双吸中开式离心泵,该泵的进出口直径均为0.35m,转速为1480rpm,额定流量为1260m3/h,该叶片泵的三维水体模型示意图见图2,三维水体网格示意图见图3。
步骤S2、以稳态计算结果为初始值,根据预设参数获取每个叶片网格单元和每个轮毂网格单元瞬态计算结果,并通过积分的方法获得叶轮的瞬态受力,所述预设参数包括时间步长和第一计算总时长;
具体地,为了获得叶轮瞬态受力,首先根据步骤S1建立的计算流体力学模型进行稳态计算,将获得的稳态计算结果作为初始值,再进行非稳态计算,获得每个叶片网格单元和每个轮毂网格单元瞬态计算结果,然后通过积分的方法获得叶轮的瞬态受力情况。额定工况下,根据计算流体力学模型计算求得的一个转动周期内叶轮瞬态受力情况见图4。
步骤S3、以叶轮的瞬态受力作为叶轮质点的载荷条件,采用有限元分析根据预设的转子系统信息对由主轴、轴承和叶轮质点组成的转子系统进行瞬态动力学分析,并获得瞬态动力学分析结果,其中,有限元分析的计算总时长与第一计算总时长相同;
具体地,转子系统由主轴、两个轴承和叶轮质点组成,其有限元模型请参见图5,在有限元分析中,将主要考虑主轴和轴承的振动情况,叶轮可简化了一个点,称为叶轮质点。应用计算机辅助制图以及数学辅助放算软件绘制主轴模型并导入到有限元模拟软件进行网格划分,主轴三维模型示意图见图6,该主轴的材质为45号钢,轴的总长度为1.3m,最大直径为0.05mm。主轴划分的有限元网格请参见图7。将所述步骤S2获取的叶轮的瞬态受力作为叶轮质点的载荷条件导入,并采用有限元分析法对转子系统进行瞬态动力学分析,并获得瞬态动力学分析结果,其中,有限元分析的计算总时长与步骤S2中的第一计算总时长相同。
步骤S4、根据瞬态动力学分析结果获得转子系统的振动情况。
计算得到的某时刻主轴变形情况请见图8,一个转动周期内叶轮质点处的主轴速度变化情况请见图9。
实验测试证明,本发明提供的叶片泵转子系统的振动仿真分析预测精度高,与实测值的偏差在15%以内,满足叶片泵子系统设计需求。
进一步地,请参照图10,基于本发明提供的叶片泵转子系统的振动仿真分析方法第一实施例,在本发明提供的叶片泵转子系统的振动仿真分析方法第二实施例中,步骤S1包括:
步骤S11、分别单独绘制叶片泵的叶轮水体、吸水室水体和压水室水体,再将三者进行装配组合为一个整体进行网格划分得到三维水体网格;
步骤S12、利用所述三维水体网格建立计算流体力学模型;
步骤S13、确定计算流体力学模型的进口边界条件和出口边界条件。
进一步地,请参照图11,基于本发明提供的叶片泵转子系统的振动仿真分析方法第二实施例,在本发明提供的叶片泵转子系统的振动仿真分析方法第三实施例中,步骤S13包括:
步骤S131、分别建立所述吸水室水体与所述叶轮水体以及所述叶轮水体与所述压水室水体的交界面联系,并以所述吸水室水体与所述叶轮水体的交界面为进口,所述叶轮水体与所述压水室水体的交界面为出口;
步骤S132、根据设计参数的工况点确定计算流体力学模型的进口边界条件和出口边界条件,所述设计参数包括流体介质的密度和粘度以及叶轮水体的旋转速度。
进一步地,请参照图12,基于本发明提供的叶片泵转子系统的振动仿真分析方法第一实施例,在本发明提供的叶片泵转子系统的振动仿真分析方法第四实施例中,步骤S2包括:
S21、根据计算流体力学模型进行稳态计算,获得稳态计算结果;
S22、以稳态计算结果为初始值,设定计算时间步长和第一计算总时长,根据计算流体力学模型计算每个叶片网格单元和每个轮毂网格单元表面水压分布,获得每个叶片网格单元水平面和竖直面轴向瞬态受力和每个轮毂网格单元水平面和竖直面轴向瞬态受力;
S23、针对每个时间步长,分别对所述叶片网格单元和所述轮毂网格单元的水平面以及所述叶片网格单元和所述轮毂网格单元的竖直面轴向瞬态受力进行积分,获得水平面和竖直面上叶轮的瞬态受力随时间的分布。
进一步地,请参照图13,基于本发明提供的叶片泵转子系统的振动仿真分析方法第一实施例,在本发明提供的叶片泵转子系统的振动仿真分析方法第五实施例中,步骤S3包括:
S31、对由主轴、轴承和叶轮质点组成的转子系统的三维有限元瞬态分析模型进行有限元网格划分;
S32、分别将水平面上叶轮的瞬态受力随时间的分布和竖直面上叶轮的瞬态受力随时间的分布导入所述三维有限元瞬态分析模型,作为叶轮质点受力的载荷条件;
S33、设定转子系统信息和计算总时长,对所述转子系统进行瞬态动力学分析并获得瞬态动力学分析结果,所述转子系统信息包括主轴转速和轴承约束条件,所述计算总时长与所述第一计算总时长相同。
本发明还提供一种叶片泵转子系统的振动仿真分析装置,请参阅图14,为本发明提供的叶片泵转子系统的振动仿真分析装置第一实施例的功能模块图。在本实施例中,所述叶片泵转子系统的振动仿真分析装置100包括:
计算流体力学模型建模模块10,用于对叶片泵的三维水体模型进行网格划分得到三维水体网格并利用所述三维水体网格建立计算流体力学模型,所述三维水体网格包括多个叶片网格单元和多个轮毂网格单元;
叶轮的瞬态受力获取模块30,用于以稳态计算结果为初始值,根据预设参数获取每个叶片网格单元和每个轮毂网格单元瞬态计算结果,并通过积分的方法得到叶轮的瞬态受力;
转子系统瞬态动力学分析模块50,用于以叶轮的瞬态受力作为叶轮质点的载荷条件,采用有限元分析根据预设的转子系统信息对由主轴、轴承和叶轮质点组成的转子系统进行瞬态动力学分析并获得瞬态动力学分析结果;
仿真分析结果获取模块70,用于根据瞬态动力学分析结果获得转子系统的振动情况。
进一步地,请参照图15,为本发明提供的叶片泵转子系统的振动仿真分析装置第二实施例的功能模块图。基于上述实施例在本实施例中,所述计算流体力学模型建模模块10包括:
三维建模单元11,用于分别单独绘制叶片泵的叶轮水体、吸水室水体和压水室水体,再将三者进行装配组合为一个整体;
网格划分单元13,用于进行网格划分得到三维水体网格;
力学模型建立单元15,用于利用所述三维水体网格建立计算流体力学模型;
边界条件确定单元17,用于确定计算流体力学模型的进口边界条件和出口边界条件。
进一步地,请参照图16,为本发明提供的叶片泵转子系统的振动仿真分析装置第三实施例的功能模块图。基于本发明叶片泵转子系统的振动仿真分析装置第二实施例,所述边界条件确定单元17包括:
交界面建立子单元171,用于建立所述吸水室水体与所述叶轮水体以及所述叶轮水体与所述压水室水体的交界面联系,并以所述吸水室水体与所述叶轮水体的交界面为进口,所述叶轮水体与所述压水室水体的交界面为出口;
条件确定子单元173,用于根据设计参数的工况点确定计算流体力学模型的进口边界条件和出口边界条件,所述设计参数包括流体介质的密度和粘度以及叶轮水体的旋转速度。
进一步地,请参照图17,为本发明提供的叶片泵转子系统的振动仿真分析装置第四实施例的功能模块图。基于本发明叶片泵转子系统的振动仿真分析装置第一实施例,所述叶轮的瞬态受力获取模块30包括:
稳态计算单元31,用于根据计算流体力学模型进行稳态计算,获得稳态计算结果;
瞬态计算单元33,用于根据计算流体力学模型计算每个叶片网格单元和每个轮毂网格单元表面水压分布,获得每个叶片网格单元水平面和竖直面轴向瞬态受力和每个轮毂网格单元水平面和竖直面轴向瞬态受力;
积分单元35,用于针对每个时间步长,分别对所述叶片网格单元和所述轮毂网格单元的水平面以及所述叶片网格单元和所述轮毂网格单元的竖直面轴向瞬态受力进行积分,获得水平面和竖直面上叶轮的瞬态受力随时间的分布。
进一步地,请参照图18,为本发明提供的叶片泵转子系统的振动仿真分析装置第五实施例的功能模块图。基于本发明叶片泵转子系统的振动仿真分析装置第一实施例,所述转子系统瞬态动力学分析模块50包括:
有限元网格划分单元51,用于对由主轴、轴承和叶轮质点组成的转子系统的三维有限元瞬态分析模型进行有限元网格划分;
叶轮质点受力导入单元53,用于将水平面上叶轮的瞬态受力随时间的分布和竖直面上叶轮的瞬态受力随时间的分布导入所述三维有限元瞬态分析模型,作为叶轮质点受力的载荷条件;
瞬态分析单元55,用于对所述转子系统进行瞬态动力学分析并获得瞬态动力学分析结果。
本发明提供的叶片泵转子系统的振动仿真分析方法及装置具有以下有益效果:本发明将计算流体力学与转子动力学有限元分析两大模型结合起来,首先通过对叶片泵进行网格划分建立计算流体力学模型,再通过计算流体力学模型得到每个叶片网格单元和每个轮毂网格单元的瞬态计算结果,通过积分的方法可获得一个周期内叶轮的瞬态受力情况,将叶轮的瞬态受力情况作为叶轮质点的载荷条件导入转子动力学有限元分析模型,最后可以获得转子振动系统的振动情况,作为叶片泵设计制造的依据。该方法不但可以快速、准确地获得水力诱导下叶片泵转子系统的振动情况,而且振动仿真分析预测精度高,与实测值的偏差在15%以内,满足叶片泵子系统设计需求。
以上所述仅为本发明的优选实施例而已,并不用于限制本发明,对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (10)
1.一种叶片泵转子系统的振动仿真分析方法,其特征在于,该方法包括如下步骤:
对叶片泵的三维水体模型进行网格划分得到三维水体网格,并利用所述三维水体网格建立计算流体力学模型,所述三维水体网格包括多个叶片网格单元和多个轮毂网格单元;
以稳态计算结果为初始值,根据预设参数获取每个叶片网格单元和每个轮毂网格单元瞬态计算结果,并通过积分的方法获得叶轮的瞬态受力,所述预设参数包括时间步长和第一计算总时长;
以叶轮的瞬态受力作为叶轮质点的载荷条件,采用有限元分析根据预设的转子系统信息对由主轴、轴承和叶轮质点组成的转子系统进行瞬态动力学分析,并获得瞬态动力学分析结果,其中,有限元分析的计算总时长与第一计算总时长相同;
根据瞬态动力学分析结果获得转子系统的振动情况。
2.根据权利要求1所述的叶片泵转子系统的振动仿真分析方法,其特征在于,所述对叶片泵的三维水体模型进行网格划分得到三维水体网格,并利用所述三维水体网格建立计算流体力学模型包括:
分别单独绘制叶片泵的叶轮水体、吸水室水体和压水室水体,再将三者进行装配组合为一个整体进行网格划分得到三维水体网格;
利用所述三维水体网格建立计算流体力学模型;
确定计算流体力学模型的进口边界条件和出口边界条件。
3.根据权利要求2所述的叶片泵转子系统的振动仿真分析方法,其特征在于,所述确定计算流体力学模型的进口边界条件和出口边界条件包括:
分别建立所述吸水室水体与所述叶轮水体以及所述叶轮水体与所述压水室水体的交界面联系,并以所述吸水室水体与所述叶轮水体的交界面为进口,所述叶轮水体与所述压水室水体的交界面为出口;
根据设计参数的工况点确定计算流体力学模型的进口边界条件和出口边界条件,所述设计参数包括流体介质的密度和粘度以及叶轮水体的旋转速度。
4.根据权利要求1所述的叶片泵转子系统的振动仿真分析方法,其特征在于,所述以稳态计算结果为初始值,根据预设参数获取每个叶片网格单元和每个轮毂网格单元瞬态计算结果,并通过积分的方法获得叶轮的瞬态受力,所述预设参数包括时间步长和第一计算总时长包括:
根据计算流体力学模型进行稳态计算,获得稳态计算结果;
以稳态计算结果为初始值,设定计算时间步长和第一计算总时长,根据计算流体力学模型计算每个叶片网格单元和每个轮毂网格单元表面水压分布,获得每个叶片网格单元水平面和竖直面轴向瞬态受力和每个轮毂网格单元水平面和竖直面轴向瞬态受力;
针对每个时间步长,分别对所述叶片网格单元和所述轮毂网格单元的水平面以及所述叶片网格单元和所述轮毂网格单元的竖直面轴向瞬态受力进行积分,获得水平面和竖直面上叶轮的瞬态受力随时间的分布图。
5.根据权利要求4所述的叶片泵转子系统的振动仿真分析方法,其特征在于,所述以叶轮的瞬态受力作为叶轮质点的载荷条件,采用有限元分析根据预设的转子系统信息对由主轴、轴承和叶轮质点组成的转子系统进行瞬态动力学分析,并获得瞬态动力学分析结果,其中,有限元分析的计算总时长与第一计算总时长相同包括:
对由主轴、轴承和叶轮质点组成的转子系统的三维有限元瞬态分析模型进行有限元网格划分;
分别将水平面上叶轮的瞬态受力随时间的分布和竖直面上叶轮的瞬态受力随时间的分布导入所述三维有限元瞬态分析模型,作为叶轮质点受力的载荷条件;
设定转子系统信息和计算总时长,对所述转子系统进行瞬态动力学分析并获得瞬态动力学分析结果,所述转子系统信息包括主轴转速和轴承约束条件,其计算总时长与所述第一计算总时长相同。
6.一种叶片泵转子系统的振动仿真分析装置,其特征在于,所述叶片泵转子系统的振动仿真分析装置包括:
计算流体力学模型建模模块,用于对叶片泵的三维水体模型进行网格划分得到三维水体网格并利用所述三维水体网格建立计算流体力学模型,所述三维水体网格包括多个叶片网格单元和多个轮毂网格单元;
叶轮的瞬态受力获取模块,用于以稳态计算结果为初始值,根据预设参数获取每个叶片网格单元和每个轮毂网格单元瞬态计算结果,并通过积分的方法得到叶轮的瞬态受力;
转子系统瞬态动力学分析模块,用于以叶轮的瞬态受力作为叶轮质点的载荷条件,采用有限元分析根据预设的转子系统信息对由主轴、轴承和叶轮质点组成的转子系统进行瞬态动力学分析并获得瞬态动力学分析结果;
仿真分析结果获取模块,用于根据瞬态动力学分析结果获得转子系统的振动情况。
7.根据权利要求6所述的叶片泵转子系统的振动仿真分析装置,其特征在于,所述计算流体力学模型建模模块包括:
三维建模单元,用于分别单独绘制叶片泵的叶轮水体、吸水室水体和压水室水体,再将三者进行装配组合为一个整体;
网格划分单元,用于进行网格划分得到三维水体网格;
力学模型建立单元,用于利用所述三维水体网格建立计算流体力学模型;
边界条件确定单元,用于确定计算流体力学模型的进口边界条件和出口边界条件。
8.根据权利要求7所述的叶片泵转子系统的振动仿真分析装置,其特征在于,所述边界条件确定单元包括:
交界面建立子单元,用于建立所述吸水室水体与所述叶轮水体以及所述叶轮水体与所述压水室水体的交界面联系,并以所述吸水室水体与所述叶轮水体的交界面为进口,所述叶轮水体与所述压水室水体的交界面为出口;
条件确定子单元,用于根据设计参数的工况点确定计算流体力学模型的进口边界条件和出口边界条件,所述设计参数包括流体介质的密度和粘度以及叶轮水体的旋转速度。
9.根据权利要求6所述的叶片泵转子系统的振动仿真分析装置,其特征在于,所述叶轮的瞬态受力获取模块包括:
稳态计算单元,用于根据计算流体力学模型进行稳态计算,获得稳态计算结果;
瞬态计算单元,用于根据计算流体力学模型计算每个叶片网格单元和每个轮毂网格单元表面水压分布,获得每个叶片网格单元水平面和竖直面轴向瞬态受力和每个轮毂网格单元水平面和竖直面轴向瞬态受力;
积分单元,用于针对每个时间步长,分别对所述叶片网格单元和所述轮毂网格单元的水平面以及所述叶片网格单元和所述轮毂网格单元的竖直面轴向瞬态受力进行积分,获得水平面和竖直面上叶轮的瞬态受力随时间的分布。
10.根据权利要求9所述的叶片泵转子系统的振动仿真分析装置,其特征在于,所述转子系统瞬态动力学分析模块包括:
有限元网格划分单元,用于对由主轴、轴承和叶轮质点组成的转子系统的三维有限元瞬态分析模型进行有限元网格划分;
叶轮质点受力导入单元,用于将水平面上叶轮的瞬态受力随时间的分布和竖直面上叶轮的瞬态受力随时间的分布导入所述三维有限元瞬态分析模型,作为叶轮质点受力的载荷条件;
瞬态分析单元,用于对所述转子系统进行瞬态动力学分析并获得瞬态动力学分析结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611099497.2A CN106599422B (zh) | 2016-12-02 | 2016-12-02 | 叶片泵转子系统的振动仿真分析方法及其装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611099497.2A CN106599422B (zh) | 2016-12-02 | 2016-12-02 | 叶片泵转子系统的振动仿真分析方法及其装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106599422A true CN106599422A (zh) | 2017-04-26 |
CN106599422B CN106599422B (zh) | 2020-05-26 |
Family
ID=58596500
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611099497.2A Active CN106599422B (zh) | 2016-12-02 | 2016-12-02 | 叶片泵转子系统的振动仿真分析方法及其装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106599422B (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108875235A (zh) * | 2018-06-28 | 2018-11-23 | 奇瑞汽车股份有限公司 | 汽车发动机的水套的三维稳态传热性能分析方法及装置 |
CN112257309A (zh) * | 2020-09-29 | 2021-01-22 | 北京科技大学 | 一种水下履带装备流场及力学响应预测方法及系统 |
CN112329357A (zh) * | 2020-10-29 | 2021-02-05 | 湘潭大学 | 一种清水离心泵振动故障的简易诊断方法及诊断系统 |
CN113297809A (zh) * | 2021-05-06 | 2021-08-24 | 湘潭大学 | 一种离心泵叶轮口环间隙泄漏量的仿真估算方法 |
CN116305641A (zh) * | 2023-03-15 | 2023-06-23 | 浙江大学 | 一种柱销式双作用叶片泵的结构参数设计方法和装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101908088A (zh) * | 2010-07-22 | 2010-12-08 | 北京航空航天大学 | 一种基于时域双向迭代的叶轮机叶片颤振应力预测方法 |
CN105956254A (zh) * | 2016-04-27 | 2016-09-21 | 西安交通大学 | 一种用于转子系统动力学建模的刚体单元法 |
US20160324452A1 (en) * | 2014-01-31 | 2016-11-10 | The University Of Toledo | Integrated Magneto-Optic Modulator/Compensator System, Methods of Making, and Methods of Using the Same |
-
2016
- 2016-12-02 CN CN201611099497.2A patent/CN106599422B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101908088A (zh) * | 2010-07-22 | 2010-12-08 | 北京航空航天大学 | 一种基于时域双向迭代的叶轮机叶片颤振应力预测方法 |
US20160324452A1 (en) * | 2014-01-31 | 2016-11-10 | The University Of Toledo | Integrated Magneto-Optic Modulator/Compensator System, Methods of Making, and Methods of Using the Same |
CN105956254A (zh) * | 2016-04-27 | 2016-09-21 | 西安交通大学 | 一种用于转子系统动力学建模的刚体单元法 |
Non-Patent Citations (2)
Title |
---|
XIN WANG 等: "Vibration Analysis of Large Bulb Tubular Pumping Station Considering Flow-Structure Interaction", 《2010 ASIA-PACIFIC POWER AND ENERGY ENGINEERING CONFERENCE》 * |
胡凤兰 等: "基于ANSYS的叶片泵有限元分析", 《CAD/CAM与制造业信息化》 * |
Cited By (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108875235A (zh) * | 2018-06-28 | 2018-11-23 | 奇瑞汽车股份有限公司 | 汽车发动机的水套的三维稳态传热性能分析方法及装置 |
CN112257309A (zh) * | 2020-09-29 | 2021-01-22 | 北京科技大学 | 一种水下履带装备流场及力学响应预测方法及系统 |
CN112257309B (zh) * | 2020-09-29 | 2023-08-25 | 北京科技大学 | 一种水下履带装备流场及力学响应预测方法及系统 |
CN112329357A (zh) * | 2020-10-29 | 2021-02-05 | 湘潭大学 | 一种清水离心泵振动故障的简易诊断方法及诊断系统 |
CN113297809A (zh) * | 2021-05-06 | 2021-08-24 | 湘潭大学 | 一种离心泵叶轮口环间隙泄漏量的仿真估算方法 |
CN113297809B (zh) * | 2021-05-06 | 2022-04-26 | 湘潭大学 | 一种离心泵叶轮口环间隙泄漏量的仿真估算方法 |
CN116305641A (zh) * | 2023-03-15 | 2023-06-23 | 浙江大学 | 一种柱销式双作用叶片泵的结构参数设计方法和装置 |
CN116305641B (zh) * | 2023-03-15 | 2024-02-13 | 浙江大学 | 一种柱销式双作用叶片泵的结构参数设计方法和装置 |
Also Published As
Publication number | Publication date |
---|---|
CN106599422B (zh) | 2020-05-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106599422A (zh) | 叶片泵转子系统的振动仿真分析方法及其装置 | |
CN103226634B (zh) | 基于三维动网格的旋喷泵非定常流场的计算方法 | |
CN105201901B (zh) | 基于流固耦合计算的离心泵叶轮水力及结构设计方法 | |
Liu et al. | Three-dimensional flow simulation of transient power interruption process of a prototype pump-turbine at pump mode | |
CN103226635B (zh) | 基于三维动网格的旋转叶轮机械非定常流场的计算方法 | |
CN104166752B (zh) | 液力变矩器全流道瞬态数值模拟计算方法 | |
CN109033488B (zh) | 一种基于水机电耦合模型的抽水蓄能机组寿命分析方法 | |
CN107122512B (zh) | 液环泵非稳态气体流场及吸入压缩性能的简化计算方法 | |
CN113048086B (zh) | 一种基于径向基神经网络模型的低噪声不等距离心风扇优化设计方法 | |
CN102141064A (zh) | 空间过滤法建立湍流模型的构建方法 | |
CN108506223A (zh) | 基于弱可压流动分析的大功率离心泵水力性能预测方法 | |
CN107273565A (zh) | 一种基于les/lbm的离心泵内流动诱导噪声数值计算方法 | |
CN105889136B (zh) | 一种用来降低离心泵内部压力脉动的离心泵改进方法 | |
Khare et al. | CFD approach for flow characteristics of hydraulic francis turbine | |
Li et al. | Stability optimization and analysis of a bidirectional shaft extension pump | |
Jemal et al. | COMPREHENSIVE REVIEW OF PUMP AS TURBINE (PAT). | |
CN112412829A (zh) | 基于物质平均大涡模拟的离心泵水力性能测量方法 | |
Zhou et al. | Internal flow numerical simulation of double-suction centrifugal pump using DES model | |
Kaniecki et al. | Computational fluid dynamics simulations of small capacity Kaplan turbines | |
Hermez et al. | Comparison of inlet curved disk arrangements for suppression of recirculation in centrifugal pump impellers | |
Chen et al. | 3-D Simulation of a prototype pump-turbine during starting period in turbine model | |
Anagnostopoulos | Numerical calculation of the flow in a centrifugal pump impeller using Cartesian grid | |
Krzemianowski et al. | Experimental validation of numerical model within a flow configuration of the model Kaplan turbine | |
El-Hennawi et al. | NUMERICAL AND EXPERIMENTAL STUDY OF A SMALL-SCALE SLIDING VANE PUMP | |
CN117648780B (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 | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20220317 Address after: 411100 office building, No. 30, bonded Road, Heping Street, Xiangtan Economic Development Zone, Xiangtan City, Hunan Province Patentee after: HUNAN M&W PUMPS CO.,LTD. Address before: 410528 No. 14-e043, 14th floor, Yannong complex building, intersection of lushong road and Luquan Road, high tech Development Zone, Changsha, Hunan Patentee before: CHANGSHA MOUNTAIN WATER ENERGY-SAVING RESEARCH INSTITUTE Co.,Ltd. |