CN115203851A - 一种大尺寸风扇转子叶片扭转刚度确定方法 - Google Patents

一种大尺寸风扇转子叶片扭转刚度确定方法 Download PDF

Info

Publication number
CN115203851A
CN115203851A CN202210869786.5A CN202210869786A CN115203851A CN 115203851 A CN115203851 A CN 115203851A CN 202210869786 A CN202210869786 A CN 202210869786A CN 115203851 A CN115203851 A CN 115203851A
Authority
CN
China
Prior art keywords
rotor blade
fan rotor
torsional rigidity
torsional
shoulder
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
Application number
CN202210869786.5A
Other languages
English (en)
Other versions
CN115203851B (zh
Inventor
沈锡钢
苏守丽
刘一雄
张海洋
伊锋
郭劼
苏军
曹航
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
AECC Shenyang Engine Research Institute
Original Assignee
AECC Shenyang Engine Research Institute
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by AECC Shenyang Engine Research Institute filed Critical AECC Shenyang Engine Research Institute
Priority to CN202210869786.5A priority Critical patent/CN115203851B/zh
Publication of CN115203851A publication Critical patent/CN115203851A/zh
Application granted granted Critical
Publication of CN115203851B publication Critical patent/CN115203851B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • 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
    • 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
    • G06F2111/00—Details relating to CAD techniques
    • G06F2111/04—Constraint-based CAD
    • 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/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Geometry (AREA)
  • General Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Structures Of Non-Positive Displacement Pumps (AREA)

Abstract

本申请提供了一种大尺寸风扇转子叶片扭转刚度确定方法，该大尺寸风扇转子叶片具有凸肩，包括：获取大尺寸风扇转子叶片的结构模型参数及材料数据；构建大尺寸风扇转子叶片扭转刚度仿真分析模型，通过在扭转刚度仿真分析模型中施加不同扭矩获取不同的扭转变形角和扭转刚度；确定风扇转子叶片的扭转刚度试验参数；开展大尺寸风扇转子叶片的扭转刚度试验，测定大尺寸风扇转子叶片的扭转角及应力分布；将仿真分析结果与试验测定结果进行对比；根据试验测定结果修正扭转刚度仿真分析模型，修正内容包括载荷加载位置、力臂取值、凸肩变形取值以及榫头约束方法；结合发动机的工作状态，利用修正后的扭转刚度仿真模型，对风扇转子叶片进行扭转刚度分析。

Description

一种大尺寸风扇转子叶片扭转刚度确定方法
技术领域
本申请属于航空发动机技术领域，特别涉及一种大尺寸风扇转子叶片扭转刚度确定方法。
背景技术
叶片凸肩作为一种常用的减振结构，在航空发动机叶片设计中广泛应用。带凸肩的大尺寸风扇转子叶片的扭转刚度是重要的结构设计参数之一，风扇转子叶片需要有足够的刚度以维持凸肩的挤压力，以增强风扇转子叶片的干摩擦阻尼减振效果，避免因为刚度不足而出现失稳、颤振等灾难性事故。因此，在风扇转子叶片设计过程中，需要通过控制风扇转子叶片的扭转刚度防止发生故障。
然而，由于大尺寸风扇转子叶片尺寸相比于现有叶片来说，其展弦比小、扭转恢复明显，现有的分析方法无法满足大尺寸风扇转子叶片的扭转刚度精度设计要求，也未采用试验数据对仿真模型进行修正。因此目前工程上尚无有效可实施的大尺寸风扇转子叶片的扭转刚度设计方法。
发明内容
本申请的目的是提供了一种大尺寸风扇转子叶片扭转刚度确定方法，以解决或减轻背景技术中的至少一个问题。
本申请的技术方案是：一种大尺寸风扇转子叶片扭转刚度确定方法，大尺寸风扇转子叶片具有凸肩，包括：
获取所述大尺寸风扇转子叶片的结构模型参数及材料数据；
根据所述大尺寸风扇转子叶片的结构模型参数及材料数据构建所述大尺寸风扇转子叶片扭转刚度仿真分析模型，通过在所述扭转刚度仿真分析模型中施加不同扭矩获取不同的扭转变形角和扭转刚度，从而形成不同工况下的扭矩-扭转刚度拟合曲线；
根据试验设备能力及测量误差确定风扇转子叶片的扭转刚度试验参数；
开展所述大尺寸风扇转子叶片的扭转刚度试验，测定所述大尺寸风扇转子叶片的扭转角及应力分布；
将所述大尺寸风扇转子叶片的仿真分析结果与试验测定结果进行对比；
根据试验测定结果修正所述大尺寸风扇转子叶片的扭转刚度仿真分析模型，修正内容包括载荷加载位置、力臂取值、凸肩变形取值以及榫头约束方法；
结合发动机的工作状态，利用修正后的扭转刚度仿真模型，对风扇转子叶片进行扭转刚度分析，获取不同工作状态叶片的扭转刚度，评价凸肩设计的可行性。
进一步的，开展所述大尺寸风扇转子叶片的扭转刚度试验时，记录施加的扭矩并测量所述大尺寸风扇转子叶片的扭转角度，扭矩施加于凸肩工作面。
进一步的，在施加扭矩时扭矩具有约束上限，以防止试验过程中扭矩过大损坏所述大尺寸风扇转子叶片。
进一步的，对比的所述大尺寸风扇转子叶片的仿真分析结果与试验测定结果包括扭转刚度、应力最大点位置、应力分布规律。
进一步的，所述载荷加载位置修正为：
在凸肩工作面施加集中力，且该集中力作用于凸肩坐标原点积叠轴处。
进一步的，所述力臂取值修正为：
根据测量的凸肩两个工作面集中力作用点处的距离，并利用啮合角进行折算得到扭矩力臂。
进一步的，所述凸肩变形取值修正为：
在进行扭转角度计算时，凸肩变形取值在凸肩工作面处选取，盆、背侧位置在凸肩坐标原点对称位置。
进一步的，所述榫头约束方法修正为：
榫头约束条件与试验状态相同，采取固支夹持。
本申请提供的大尺寸风扇转子叶片扭转刚度设计方法可用于指导大尺寸带凸肩的风扇转子叶片扭转刚度设计，该方法基于试验结果修正扭转刚度仿真模型，从而可以提升扭转刚度仿真分析精度，提高风扇叶片设计水平，为大涵道比航空发动机研制提供坚实基础。
附图说明
为了更清楚地说明本申请提供的技术方案，下面将对附图作简单地介绍。显而易见地，下面描述的附图仅仅是本申请的一些实施例。
图1为本申请的大尺寸风扇转子叶片扭转刚度设计方法流程图。
图2为本申请中的扭转方向示意图。
具体实施方式
为使本申请实施的目的、技术方案和优点更加清楚，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行更加详细的描述。
为了满足带凸肩的大尺寸风扇转子叶片的扭转刚度设计，更好的支撑发动机设计，本申请中提出了一种工程上适用的带凸肩大尺寸风扇转子叶片扭转刚度设计方法。
如图1所示，本申请提供的大尺寸风扇转子叶片扭转刚度设计方法包括以下步骤：
S1、数据输入；
获取大尺寸风扇转子叶片(以下或简称风扇转子叶片)的结构模型参数或数据及材料数据。
S2、装配状态扭转刚度分析；
构建风扇转子叶片的扭转刚度仿真分析模型，通过有限元分析法进行装配状态下的风扇转子叶片扭转刚度分析，通过施加不同扭矩获取不同的扭转变形角和扭转刚度，从而形成不同工况下的扭矩-扭转刚度拟合曲线。
S3、确定风扇转子叶片的扭转刚度试验参数；
扭转刚度的测试需要施加扭矩并测量扭转角度，试验时以控制扭转角度为主，给定扭转角度时，需确保扭转方向要使凸肩工作面受压，如图2所示。
考虑到扭转刚度及测试扭转角度均存在一定的误差，试验扭转角度需综合考虑试验设备能力及测量误差给出。为防止试验过程中扭矩过大损坏叶片，还需明确扭矩约束上限。
S4、扭转刚度试验；
风扇转子叶片扭转刚度试验在室温条件下进行，试验时通过夹紧风扇转子叶片的榫头，在风扇转子叶片的凸肩工作面施加扭矩，从而测定风扇转子叶片的扭转角及应力分布。
S5、理论计算与试验结果对比分析；
将风扇转子叶片的理论计算结果或仿真计算结果与试验结果进行对比，对比内容包括扭转刚度、应力最大点位置、应力分布规律等，得到两者之间的差异。
在本申请中，可在试验的基础上通过记录不同扭转角度及应力部分值等数据得到风扇转子叶片实测的扭转刚度分散度，以了解风扇转子叶片的工艺分散性。
S6、根据试验测试结果修正风扇转子叶片的扭转刚度仿真分析模型；
扭转刚度试验为凸肩设计的重要参数之一，修正扭转刚度仿真分析模型根据仿真分析与试验测量的误差对凸肩接触应力的影响程度确定。
采用试验结果修正时，参考因素包括载荷的加载位置、力臂的取值、凸肩变形的取值以及榫头的约束方法等，具体包括：
1)在凸肩工作面尽量施加集中力，该集中力作用于凸肩坐标原点O’积叠轴处；
2)测量凸肩两个工作面集中力作用点处的距离，并利用啮合角进行折算得到扭矩力臂；
3)在进行扭转角度计算时，凸肩变形取值在凸肩工作面处选取，盆、背侧位置尽量在凸肩坐标原点O’对称位置；
4)扭转刚度为叶片的固有特性，不受发动机状态离心负荷影响，为保证计算精度，榫头约束条件与试验状态相同，通常情况下榫头采取固支夹持。
S7、采用修正后的模型，进行工作状态下的风扇转子叶片扭转刚度分析。
风扇转子叶片的扭转刚度选取影响凸肩工作面接触应力及紧度，从而影响凸肩的减振效果，因此，需结合发动机的工作状态，利用修正后的扭转刚度仿真模型，对风扇转子叶片进行扭转刚度分析，获取不同工作状态叶片的扭转刚度，评价凸肩设计的可行性。
本申请提供的大尺寸风扇转子叶片扭转刚度设计方法可用于指导大尺寸带凸肩的风扇转子叶片扭转刚度设计，该方法基于试验结果修正扭转刚度仿真模型，从而可以提升扭转刚度仿真分析精度，提高风扇叶片设计水平，为大涵道比航空发动机研制提供坚实基础。
以上所述，仅为本申请的具体实施方式，但本申请的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本申请揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本申请的保护范围之内。因此，本申请的保护范围应以所述权利要求的保护范围为准。

Claims (8)

1.一种大尺寸风扇转子叶片扭转刚度确定方法，大尺寸风扇转子叶片具有凸肩，其特征在于，包括：
获取所述大尺寸风扇转子叶片的结构模型参数及材料数据；
根据所述大尺寸风扇转子叶片的结构模型参数及材料数据构建所述大尺寸风扇转子叶片扭转刚度仿真分析模型，通过在所述扭转刚度仿真分析模型中施加不同扭矩获取不同的扭转变形角和扭转刚度，从而形成不同工况下的扭矩-扭转刚度拟合曲线；
根据试验设备能力及测量误差确定风扇转子叶片的扭转刚度试验参数；
开展所述大尺寸风扇转子叶片的扭转刚度试验，测定所述大尺寸风扇转子叶片的扭转角及应力分布；
将所述大尺寸风扇转子叶片的仿真分析结果与试验测定结果进行对比；
根据试验测定结果修正所述大尺寸风扇转子叶片的扭转刚度仿真分析模型，修正内容包括载荷加载位置、力臂取值、凸肩变形取值以及榫头约束方法；
结合发动机的工作状态，利用修正后的扭转刚度仿真模型，对风扇转子叶片进行扭转刚度分析，获取不同工作状态叶片的扭转刚度，评价凸肩设计的可行性。
2.如权利要求1所述的大尺寸风扇转子叶片扭转刚度确定方法，其特征在于，开展所述大尺寸风扇转子叶片的扭转刚度试验时，记录施加的扭矩并测量所述大尺寸风扇转子叶片的扭转角度，扭矩施加于凸肩工作面。
3.如权利要求2所述的大尺寸风扇转子叶片扭转刚度确定方法，其特征在于，在施加扭矩时扭矩具有约束上限，以防止试验过程中扭矩过大损坏所述大尺寸风扇转子叶片。
4.如权利要求1所述的大尺寸风扇转子叶片扭转刚度确定方法，其特征在于，对比的所述大尺寸风扇转子叶片的仿真分析结果与试验测定结果包括扭转刚度、应力最大点位置、应力分布规律。
5.如权利要求1所述的大尺寸风扇转子叶片扭转刚度确定方法，其特征在于，所述载荷加载位置修正为：
在凸肩工作面施加集中力，且该集中力作用于凸肩坐标原点积叠轴处。
6.如权利要求1所述的大尺寸风扇转子叶片扭转刚度确定方法，其特征在于，所述力臂取值修正为：
根据测量的凸肩两个工作面集中力作用点处的距离，并利用啮合角进行折算得到扭矩力臂。
7.如权利要求1所述的大尺寸风扇转子叶片扭转刚度确定方法，其特征在于，所述凸肩变形取值修正为：
在进行扭转角度计算时，凸肩变形取值在凸肩工作面处选取，盆、背侧位置在凸肩坐标原点对称位置。
8.如权利要求1所述的大尺寸风扇转子叶片扭转刚度确定方法，其特征在于，所述榫头约束方法修正为：
榫头约束条件与试验状态相同，采取固支夹持。
CN202210869786.5A 2022-07-22 2022-07-22 一种大尺寸风扇转子叶片扭转刚度确定方法 Active CN115203851B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210869786.5A CN115203851B (zh) 2022-07-22 2022-07-22 一种大尺寸风扇转子叶片扭转刚度确定方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210869786.5A CN115203851B (zh) 2022-07-22 2022-07-22 一种大尺寸风扇转子叶片扭转刚度确定方法

Publications (2)

Publication Number Publication Date
CN115203851A true CN115203851A (zh) 2022-10-18
CN115203851B CN115203851B (zh) 2023-04-07

Family

ID=83584813

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210869786.5A Active CN115203851B (zh) 2022-07-22 2022-07-22 一种大尺寸风扇转子叶片扭转刚度确定方法

Country Status (1)

Country Link
CN (1) CN115203851B (zh)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101510230A (zh) * 2009-03-11 2009-08-19 同济大学 车辆道路载荷的仿真方法
US10808681B2 (en) * 2018-01-23 2020-10-20 General Electric Company Twist correction factor for aerodynamic performance map used in wind turbine control
CN112182749A (zh) * 2020-09-23 2021-01-05 吉林大学 赛车车架性能分析方法、装置、设备及可存储介质
CN113009827A (zh) * 2021-02-19 2021-06-22 济南中科数控设备有限公司 一种基于虚拟调试的数控雕刻机动态性能优化方法
CN114528641A (zh) * 2022-02-25 2022-05-24 重庆长安汽车股份有限公司 一种平台化车身结构的优化方法及存储介质
CN114636623A (zh) * 2020-12-15 2022-06-17 中国航发商用航空发动机有限责任公司 航空发动机叶片抗冲击试验方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101510230A (zh) * 2009-03-11 2009-08-19 同济大学 车辆道路载荷的仿真方法
US10808681B2 (en) * 2018-01-23 2020-10-20 General Electric Company Twist correction factor for aerodynamic performance map used in wind turbine control
CN112182749A (zh) * 2020-09-23 2021-01-05 吉林大学 赛车车架性能分析方法、装置、设备及可存储介质
CN114636623A (zh) * 2020-12-15 2022-06-17 中国航发商用航空发动机有限责任公司 航空发动机叶片抗冲击试验方法
CN113009827A (zh) * 2021-02-19 2021-06-22 济南中科数控设备有限公司 一种基于虚拟调试的数控雕刻机动态性能优化方法
CN114528641A (zh) * 2022-02-25 2022-05-24 重庆长安汽车股份有限公司 一种平台化车身结构的优化方法及存储介质

Also Published As

Publication number Publication date
CN115203851B (zh) 2023-04-07

Similar Documents

Publication Publication Date Title
Rezaei et al. Modal-based damage identification for the nonlinear model of modern wind turbine blade
US20120303314A1 (en) Angle measurement device, angle measurement method and construction machinery
US20170284375A1 (en) Method and arrangement for continuous calibration of a wind direction measurement
JP2008102049A (ja) バランス修正装置
EP3225838A1 (en) Method and arrangement for performing a wind direction measurement
CN115203851B (zh) 一种大尺寸风扇转子叶片扭转刚度确定方法
CN111191391B (zh) 一种基于局部应力分布加权系数的模拟件设计方法
Kyoung et al. Time domain structural fatigue analysis of floating offshore platforms: A response based technique
CN115293061A (zh) 用于确定塔式起重机风载荷的方法、处理器及存储介质
CN104749959B (zh) 一种基于广义滑模估计器的机组变桨容错控制方法
KR101733361B1 (ko) 로터 해석 방법
CN111898238B (zh) 一种有约束的高速动平衡力学解算方法
CN112395754B (zh) 一种机械连接结构结合面刚度动态弱化预测方法
CN110489916B (zh) 基于损伤状态倾角影响线曲率的等截面梁损伤识别方法
CN111428413B (zh) 连续梁分段抗弯刚度识别方法
CN110501127B9 (zh) 一种基于损伤状态倾角斜率的等截面梁损伤识别方法
RU2410661C2 (ru) Способ контроля профиля зоны соединения между цилиндрической частью и затылованной поверхностью детали газотурбинного двигателя
CN112594015A (zh) 一种汽轮机寿命预测方法及系统
US20230008848A1 (en) Method and apparatus for determining upvalue factors for expansion measurements on machine elements
KR101699856B1 (ko) 말뚝 및 타워형 구조물의 최 상단 가속도와 각속도를 이용한 손상 추정 방법
Szwedowicz et al. Numerical and experimental damping assessment of a thin-walled friction damper in the rotating set-up with high pressure turbine blades
US20080049887A1 (en) Method for Monitoring the State of Turbines Based on their Coasting Time
CN112432746A (zh) 一种确定直升机桨叶挥摆载荷耦合程度的方法
KR101790733B1 (ko) 다단계 클러치 댐퍼 모델을 가지는 스텝 길이 제어 장치 및 그 방법
KR20210082948A (ko) 수치해석을 이용한 선박풍속 교정 장치 및 방법

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