适用于分析叶片-楔形缘板阻尼器的非线性接触模型构建方法及系统

技术特征：
1.一种适用于分析叶片-楔形缘板阻尼器的非线性接触模型构建方法，其特征在于，包括：步骤1、获取叶片-楔形缘板阻尼器网格数据结构对象及动力学参数；步骤2、基于所述网格数据结构和动力学参数构建叶片-楔形缘板阻尼器动力学模型；步骤3、基于叶片-楔形缘板阻尼器动力学模型，构建叶片-楔形缘板阻尼器非线性接触响应模型；步骤4、基于所述模型、网格、节点及所述动力学参数输出非线性分析文件，以供非线性接触分析求解器生成非线性接触分析模型并进行数值计算。2.根据权利要求1所述的一种适用于分析叶片-楔形缘板阻尼器的非线性接触模型构建方法，其特征在于，所述叶片-楔形缘板阻尼器网格数据结构对象及动力学参数，具体包括：楔形缘板阻尼器与两侧叶片接触面上的所有网格节点编号和所述接触面网格节点对应的三维坐标；所有所述叶片经过固定界面模态综合法缩减后的网格节点编号和所述叶片网格节点对应的三维坐标；所述阻尼器的几何参数及所述叶片有限元网格节点对应的质量矩阵、刚度矩阵；以及所述叶片-楔形缘板阻尼器有限元模型在压紧力作用下接触面网格节点对应的节点力矩阵。3.根据权利要求2所述的一种适用于分析叶片-楔形缘板阻尼器的非线性接触模型构建方法，其特征在于，所述步骤2具体包括：步骤201、根据叶片-楔形缘板阻尼器网格数据结构对象及动力学参数，进行离心力作用下所述阻尼器与叶片的静力结构计算，得到所述阻尼结构在柔性条件下的静变形和实际接触正压力；步骤202、根据前述阻尼器接触面网格节点的静变形和接触力计算结果，判定所述阻尼器与两侧叶片接触面上网格节点的接触状态，据此在所述接触面上选择恰当的接触点，建立简化的叶片-楔形缘板阻尼器接触运动模型；步骤203、建立所述叶片-楔形缘板阻尼器静平衡状态下广义位移边界条件，及静平衡状态下所述接触面的受力边界条件。4.根据权利要求3所述的一种适用于分析叶片-楔形缘板阻尼器的非线性接触模型构建方法，其特征在于，选择接触点的位置和数目，具体包括：步骤a、获取所述阻尼结构与两侧缘板接触面正压力分布，判定接触面的接触状态；步骤b、将所述接触面按接触状态及正压力分布划分区域，得到接触应力相近的接触区域；步骤c、对所述接触区域计算正压力合力位置坐标，在所述合力位置设定接触点；步骤d、在所述接触点建立粘滞-滑移单元表征接触对间作用力和力矩。5.根据权利要求4所述的一种适用于分析叶片-楔形缘板阻尼器的非线性接触模型构建方法，其特征在于，所述步骤3具体包括：步骤301、建立所述阻尼器两侧叶片的离散的非线性振动微分方程组，采用数值谐波平衡法将叶片离散的非线性振动微分方程组转化为非线性代数方程组，取无干摩擦情况下所
述阻尼器两侧叶片的响应或上一频率点的稳态响应作为初值；步骤302、将所述阻尼器两侧的缘板处接触点的振动响应转化到时域中，结合所述楔形缘板阻尼结构的动力学模型，应用变步长dormand-prince runge-kutta 45求解方法在时域内求解摩擦力，然后通过线性插值得到符合采样要求的等步长的时域信号傅里叶变换求解成本，再通过傅里叶变换将所述非线性摩擦力时域信号投影到频域中，重复循环此变步长时频域迭代求解过程直到满足收敛条件，得到振动位移与非线性摩擦力之间的显式关系。6.根据权利要求1至5之一所述的一种适用于分析叶片-楔形缘板阻尼器的非线性接触模型构建方法，其特征在于，所述模型、网格、节点及所述动力学参数输出非线性分析文件，以供非线性接触分析求解器生成非线性接触分析模型并进行数值计算。7.一种根据权利要求1-6之一所述的适用于分析叶片-楔形缘板阻尼器的非线性接触模型构建方法的构建系统，其特征在于，包括：获取模块，用于获取叶片-楔形缘板阻尼器网格数据结构对象及动力学参数；构建模块，用于基于所述网格数据结构和动力学参数构建叶片-楔形缘板阻尼器动力学模型；响应模块，用于基于叶片-楔形缘板阻尼器动力学模型，构建叶片-楔形缘板阻尼器非线性接触响应模型；输出模块，用于基于所述模型、网格、节点及所述动力学参数输出非线性分析文件，以供非线性接触分析求解器生成非线性接触分析模型并进行数值计算。8.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现如权利要求1至6任一项所述的一种适用于分析叶片-楔形缘板阻尼器的非线性接触模型构建方法的步骤。9.一种非暂态计算机可读存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现如权利要求1至6任一项的所述一种适用于分析叶片-楔形缘板阻尼器的非线性接触模型构建方法的步骤。

技术总结
本发明提供一种适用于分析叶片-楔形缘板阻尼器的非线性接触模型构建方法及系统。获取叶片-楔形缘板阻尼器网格数据结构对象及动力学参数；基于所述网格数据结构和动力学参数构建叶片-楔形缘板阻尼器动力学模型；基于叶片-楔形缘板阻尼器动力学模型，构建叶片-楔形缘板阻尼器非线性接触响应模型；基于所述模型、网格、节点及所述动力学参数输出非线性分析文件，以供非线性接触分析求解器生成非线性接触分析模型并进行数值计算。本发明通过建立叶片-楔形缘板阻尼器的接触运动模型、发展多重减缩技术和时频域求解技术，提出了能反映阻尼器滚转、惯性和柔性影响的求解效率更高、更稳健的叶片-楔形缘板阻尼器动力学模型的强迫振动响应求解方法。动响应求解方法。动响应求解方法。


技术研发人员：张大义 孙赫 佟若昀
受保护的技术使用者：北京航空航天大学
技术研发日：2022.09.14
技术公布日：2022/11/11