申请/专利权人：北京航空航天大学

申请日：2021-11-15

公开（公告）日：2024-06-21

公开（公告）号：CN114186449B

主分类号：G06F30/23

分类号：G06F30/23;G06F30/10;G06F17/11;G06F17/16;G06F119/14

优先权：

专利状态码：有效-授权

法律状态：2024.06.21#授权;2022.04.01#实质审查的生效;2022.03.15#公开

摘要：本发明公开了一种基于梁弯曲理论的叶盘扇区间刚度评估方法，属于叶轮机械叶盘结构设计技术领域，所述方法包括：根据真实叶片三维几何模型，划分叶片网格模型并计算模态频率，同时提取真实轮盘三维几何模型的特征参数，建立简化的轮盘参数模型，并构建原始叶盘集总参数模型；结合梁弯曲理论计算得到不同模态振型的各个扇区之间轴向等效连接刚度；对原始叶盘集总参数模型的扇区间连接刚度进行更新，计算得到新刚度下的叶盘模态频率和振型，并重新基于梁弯曲理论计算得到新的等效刚度值，并与上一步结果对比。本发明在保证集总参数模型原有计算效率的前提下，提高了叶盘模态频率和振型的计算精度，为叶盘模态特性分析提供了一个新途径。

主权项：1.一种基于梁弯曲理论的叶盘扇区间刚度评估方法，其特征在于，包括：步骤1：根据真实叶片三维几何模型划分叶片网格模型，计算叶片的一阶弯曲模态频率；步骤2：提取真实轮盘三维几何模型的特征参数，建立简化的轮盘参数模型；步骤3：基于简化的轮盘参数模型和叶片一阶弯曲模态频率，构建原始叶盘集总参数模型；步骤4：基于原始叶盘集总参数模型计算叶盘的模态特性，包括模态频率和振型；步骤5：基于梁弯曲理论和步骤4得到的模态振型，计算不同模态振型下各个扇区之间的轴向等效连接刚度；步骤6：对原始叶盘集总参数模型的扇区间连接刚度进行更新，重复步骤4计算得到新的模态频率和振型，重复步骤5得到新的等效连接刚度值，并与上一步结果对比，若两次结果的差异小于所要求的量级，将当前的等效连接刚度值作为结果输出；否则，再次执行步骤6，直至结果满足要求；其中，所述步骤5包括：假设各个叶片由于变形而作用在轮盘轴向方向上的力视为集中力，则将轮盘外缘在这些力作用下的轴向变形简化为梁受等间隔分布的集中力作用下的弯曲变形；下面对梁的受力状态进行分析，根据梁弯曲变形时的内力F、弯矩W、转角θ和挠度s之间的关系，得到：内力：假设任意两个相邻叶片之间的梁段的内力F＝C1；弯矩：对内力在长度l上积分，并加上该梁段起始位置处的弯矩C2后得到该梁段的弯矩表达式W＝C1l+C2；转角：令弯矩除以弹性模量E和惯性矩I，然后在长度l上积分并加上该梁段起始位置处的转角C3后得到该梁段的转角表达式挠度：令转角在长度l上积分并加上该梁段起始位置处的挠度C4后得到该梁段的挠度表达式因此，任一梁段的挠度表达式都有C1、C2、C3和C4共4个未知参数，N个叶片的集中力将梁分为N个梁段，则共有4×N个未知参数，则需要构建4×N个方程来进行求解；对于一个给定的模态振型，带入2类已知条件：1第n个梁段介于第n和第n+1个叶片之间，其在第n-1个叶片和在第n个叶片处的挠度与其模态位移相等，得到2×N个方程；2第n-1个梁段和第n个梁段在第n个叶片位置处的转角和弯矩相等，得到2×N个方程；联立上述4×N个方程得到具有4×N个未知数的方程组，求解该方程组得到各个梁段的4个未知参数，其中C1为该梁段的内力，然后从第n个梁段到第n+1个梁段的内力变化主要由两个梁段连接处，即第n个叶片位置，对梁段所作用的集中力导致，故第n个叶片位置处的集中力由两侧的内力求差得到： 最终，第n个叶片位置处的受力和挠度关系由下式表示： 其中，为第n-1和第n个叶片之间的轴向等效连接刚度，为轮盘在第n-1和第n个叶片位置处的挠度的差，为第n和第n+1个叶片之间的轴向等效连接刚度，为轮盘在第n和第n+1个叶片位置处的挠度的差；联立N个叶片位置处的表达式得到由N个方程组成的具有N个未知数的线性方程组，但其中相互独立的方程只有N-1个，因此需要补充一个方程，本方法通过遍历当前模态下各个节点的模态位移，找到某个与其相邻两个节点的模态位移的差和最接近的节点p，令该节点与相邻两个节点的轴向等效连接刚度和相等，则上述方程组可解，通过数学方法求解方程组得到各个叶片之间的轴向等效连接刚度。

全文数据：

权利要求：

百度查询： 北京航空航天大学 基于梁弯曲理论的叶盘扇区间刚度评估方法

免责声明
1、本报告根据公开、合法渠道获得相关数据和信息，力求客观、公正，但并不保证数据的最终完整性和准确性。
2、报告中的分析和结论仅反映本公司于发布本报告当日的职业理解，仅供参考使用，不能作为本公司承担任何法律责任的依据或者凭证。

阅读全文 双屏查看 官方信息 专利公告 收藏专利 下载PDF 下载WORD