CN116432311A - 一种不确定参数复合材料整体叶环振动特性计算方法 - Google Patents
一种不确定参数复合材料整体叶环振动特性计算方法 Download PDFInfo
- Publication number
- CN116432311A CN116432311A CN202310306832.5A CN202310306832A CN116432311A CN 116432311 A CN116432311 A CN 116432311A CN 202310306832 A CN202310306832 A CN 202310306832A CN 116432311 A CN116432311 A CN 116432311A
- Authority
- CN
- China
- Prior art keywords
- composite material
- blisk
- calculating
- vibration characteristics
- composite
- 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.)
- Pending
Links
- 239000002131 composite material Substances 0.000 title claims abstract description 125
- 238000000034 method Methods 0.000 title claims abstract description 42
- 239000006185 dispersion Substances 0.000 claims abstract description 11
- 238000013461 design Methods 0.000 claims abstract description 10
- 238000005259 measurement Methods 0.000 claims abstract description 7
- 238000004364 calculation method Methods 0.000 claims description 15
- 238000012360 testing method Methods 0.000 claims description 13
- 239000000463 material Substances 0.000 claims description 7
- 238000010586 diagram Methods 0.000 claims description 4
- 238000009864 tensile test Methods 0.000 claims description 3
- 230000002265 prevention Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 3
- 238000000342 Monte Carlo simulation Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000009795 derivation Methods 0.000 description 2
- 239000000835 fiber Substances 0.000 description 2
- HBMJWWWQQXIZIP-UHFFFAOYSA-N silicon carbide Chemical compound [Si+]#[C-] HBMJWWWQQXIZIP-UHFFFAOYSA-N 0.000 description 2
- 229910010271 silicon carbide Inorganic materials 0.000 description 2
- 229910001069 Ti alloy Inorganic materials 0.000 description 1
- RTAQQCXQSZGOHL-UHFFFAOYSA-N Titanium Chemical compound [Ti] RTAQQCXQSZGOHL-UHFFFAOYSA-N 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 239000011156 metal matrix composite Substances 0.000 description 1
- 230000003014 reinforcing effect Effects 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 239000010936 titanium Substances 0.000 description 1
- 229910052719 titanium Inorganic materials 0.000 description 1
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/15—Vehicle, aircraft or watercraft 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
- G06F2113/00—Details relating to the application field
- G06F2113/26—Composites
-
- 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
-
- 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)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Aviation & Aerospace Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Automation & Control Theory (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
本申请属于航空发动机整体叶环领域,特别涉及一种不确定参数复合材料整体叶环振动特性计算方法。步骤一、确定复合材料整体叶环的复合材料体积占比;步骤二、对所述复合材料整体叶环进行复合材料性能实测,得到复合材料弹性模量数据;步骤三、根据所述复合材料弹性模量数据,确定复合材料变异系数;步骤四、根据所述复合材料变异系数计算出复合材料整体叶环振动特性。本申请的不确定参数复合材料整体叶环振动特性计算方法,能够高效计算具有不确定性复合材料整体叶环频率的概率分布,进而精确开展整体叶环的避振设计,从而制定复合材料整体叶环结构频率分散度标准,为复合材料整体叶环结构的避振提供支撑。
Description
技术领域
本申请属于航空发动机整体叶环领域,特别涉及一种不确定参数复合材料整体叶环振动特性计算方法。
背景技术
复合材料整体叶环是使用连续单根碳化硅长丝增强的钛合金金属基复合材料制造的,单长丝增强的复合材料在高温下在沿纤维方向具有很高的比强度和比刚度。于是为适应高性能、高推重比发动机,复合材料整体叶环用于新一代压气机结构之中。通常在复合材料整体叶环结构的强度设计之中使用确定的结构参数、载荷得到确定的分析结果。然而在实际的工程中,由于连续纤维增强的钛基复合材料的制造、工艺等各方面因素的影响,使得材料的实测性能与预期设计性能具有一定的偏差,从而进一步导致整体叶环结构与所分析模型之间不可避免地存在一定的差异,这些不确定性可能会对发动机结构的安全性、完整性产生影响。
对于不确定性问题,国内外都开展了相应的研究,主要分为模糊方法、区间方法和概率方法。模糊方法采用隶属度函数来表征差异,其缺点是对于计算量较大的动力学问题很难获取不确定参数的表征函数。区间方法可以计算不确定参数区间范围,其计算时需要给出不确定参数的上下界但存在放大区间的风险。概率方法主要包括蒙特卡洛法和摄动法,蒙特卡洛法需要较大的计算成本但稳定性好,通常用于验证其他方法。摄动法是通过泰勒展开构造递推方程、求解构造的递推方程得到不确定参数的概率分布。对于不确定性较低的结构使用摄动法具有较高的计算效率和精度。此外,现有的不确定方法通常需要较为复杂的算法研究和理论公式推导,研究周期长且存在较大的计算量,在实际使用时存在一些限制,不利于实际工程设计中的使用。
因此,希望有一种技术方案来克服或至少减轻现有技术的至少一个上述缺陷。
发明内容
本申请的目的是提供了一种不确定参数复合材料整体叶环振动特性计算方法,以解决现有技术存在的至少一个问题。
本申请的技术方案是:
一种不确定参数复合材料整体叶环振动特性计算方法,包括:
步骤一、确定复合材料整体叶环的复合材料体积占比;
步骤二、对所述复合材料整体叶环进行复合材料性能实测,得到复合材料弹性模量数据;
步骤三、根据所述复合材料弹性模量数据,确定复合材料变异系数;
步骤四、根据所述复合材料变异系数计算出复合材料整体叶环振动特性。
在本申请的至少一个实施例中,步骤二中,通过拉伸试验进行复合材料性能实测。
在本申请的至少一个实施例中,步骤三中,所述根据所述复合材料弹性模量数据,确定复合材料变异系数,包括:
计算所述复合材料弹性模量数据的均值μ以及方差σ,所述复合材料变异系数为:
Cv=μ/σ
其中,Cv为复合材料变异系数。
在本申请的至少一个实施例中,步骤四中,所述根据所述复合材料变异系数计算出复合材料整体叶环振动特性,包括:
基于所述复合材料变异系数构造摄动项,对所述复合材料整体叶环采用有限元法进行3次有限元确定性计算,得到复合材料整体叶环频率的概率分布。
在本申请的至少一个实施例中,步骤四中,所述根据所述复合材料变异系数计算出复合材料整体叶环振动特性,还包括:
获取3组复合材料整体叶环振动特性的试验结果;
将复合材料整体叶环振动特性的计算结果与3组试验结果进行对比,若试验结果在计算结果的±3σ概率范围内,则验证得到具有不确定材料参数的复合材料整体叶环频率的概率分布。
在本申请的至少一个实施例中,还包括步骤五、通过频率分散度绘制坎贝尔图进行所述复合材料整体叶环的避振设计。
在本申请的至少一个实施例中,还包括步骤六、根据所述复合材料整体叶环振动特性的计算结果以及试验结果,形成复合材料整体叶环频率分散度标准。
发明至少存在以下有益技术效果:
本申请的不确定参数复合材料整体叶环振动特性计算方法,能够高效计算具有不确定性复合材料整体叶环频率的概率分布,进而精确开展整体叶环的避振设计,从而制定复合材料整体叶环结构频率分散度标准,为复合材料整体叶环结构的避振提供支撑。
附图说明
图1是本申请一个实施方式的不确定参数复合材料整体叶环振动特性计算方法流程图;
图2是本申请一个实施方式的复合材料整体叶环示意图;
图3是本申请一个实施方式的纤维加强环结构剖视图。
具体实施方式
为使本申请实施的目的、技术方案和优点更加清楚,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行更加详细的描述。在附图中,自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。所描述的实施例是本申请一部分实施例,而不是全部的实施例。下面通过参考附图描述的实施例是示例性的,旨在用于解释本申请,而不能理解为对本申请的限制。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。下面结合附图对本申请的实施例进行详细说明。
在本申请的描述中,需要理解的是,术语“中心”、“纵向”、“横向”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系,仅是为了便于描述本申请和简化描述,而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作,因此不能理解为对本申请保护范围的限制。
下面结合附图1至图3对本申请做进一步详细说明。
本申请提供了一种不确定参数复合材料整体叶环振动特性计算方法,包括以下步骤:
步骤一、确定复合材料整体叶环的复合材料体积占比;
步骤二、对复合材料整体叶环进行复合材料性能实测,得到复合材料弹性模量数据;
步骤三、根据复合材料弹性模量数据,确定复合材料变异系数;
步骤四、根据复合材料变异系数计算出复合材料整体叶环振动特性。
本申请的不确定参数复合材料整体叶环振动特性计算方法,可以根据现有的理论公式计算得到材料性能参数。在确定复合材料整体叶环的复合材料体积占比后,通过复合材料整体叶环的拉伸试验进行复合材料性能实测,测量出复合材料弹性模量数据。
在本申请的优选实施方式中,根据实测数据确定复合材料变异系数,复合材料变异系数表征参数的离散程度,复合材料的性能数据为不确定参数,根据复合材料弹性模量数据的均值μ以及方差σ,得到不确定参数即复合材料变异系数为:
Cv=μ/σ
其中,Cv为复合材料变异系数。
在本申请的优选实施方式中,基于改进的摄动法计算复合材料整体叶环振动特性,改进的摄动法计算不确定问题更加快速,不需要通过控制方程建立递推方程、避免了对随机参数矩阵的求导,将不确定问题计算转化成为确定性问题计算,具有较高的效率。本实施例中,根据复合材料变异系数计算出复合材料整体叶环振动特性,包括:
基于复合材料变异系数构造摄动项,对复合材料整体叶环采用有限元法进行3次有限元确定性计算,得到复合材料整体叶环频率的概率分布。
进一步,本实施例中,还包括:
获取3组复合材料整体叶环振动特性的试验结果;
将复合材料整体叶环振动特性的计算结果与3组试验结果进行对比,若试验结果在计算结果的±3σ概率范围内,则验证得到具有不确定材料参数的复合材料整体叶环频率的概率分布。
本申请的不确定参数复合材料整体叶环振动特性计算方法,还包括步骤五、通过频率分散度绘制坎贝尔图进行复合材料整体叶环的避振设计。还包括步骤六、根据复合材料整体叶环振动特性的计算结果以及试验结果,形成复合材料整体叶环频率分散度标准。本实施例中,通过频率分散度绘制坎贝尔图进行整体叶环的避振设计,准确给出共振区间进而进行避振,如图2-3所示,采用碳化硅限位加强环进行避振。最后,通过计算得到的频率高斯分布结果与试验频率结果,形成复合材料整体叶环频率分散度标准。
本申请的不确定参数复合材料整体叶环振动特性计算方法,仅通过较少的计算次数就能够得到具有不确定材料参数的复合材料整体叶环的振动特性的频率分布,实用性强,效率高,计算成本低。以此制定复合材料叶环结构频率分散度标准,可以为含复合材料整体叶环结构发动机的避振设计提供支撑。
以上所述,仅为本申请的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应以所述权利要求的保护范围为准。
Claims (7)
1.一种不确定参数复合材料整体叶环振动特性计算方法,其特征在于,包括:
步骤一、确定复合材料整体叶环的复合材料体积占比;
步骤二、对所述复合材料整体叶环进行复合材料性能实测,得到复合材料弹性模量数据;
步骤三、根据所述复合材料弹性模量数据,确定复合材料变异系数;
步骤四、根据所述复合材料变异系数计算出复合材料整体叶环振动特性。
2.根据权利要求1所述的不确定参数复合材料整体叶环振动特性计算方法,其特征在于,步骤二中,通过拉伸试验进行复合材料性能实测。
3.根据权利要求2所述的不确定参数复合材料整体叶环振动特性计算方法,其特征在于,步骤三中,所述根据所述复合材料弹性模量数据,确定复合材料变异系数,包括:
计算所述复合材料弹性模量数据的均值μ以及方差σ,所述复合材料变异系数为:
Cv=μ/σ
其中,Cv为复合材料变异系数。
4.根据权利要求3所述的不确定参数复合材料整体叶环振动特性计算方法,其特征在于,步骤四中,所述根据所述复合材料变异系数计算出复合材料整体叶环振动特性,包括:
基于所述复合材料变异系数构造摄动项,对所述复合材料整体叶环采用有限元法进行3次有限元确定性计算,得到复合材料整体叶环频率的概率分布。
5.根据权利要求4所述的不确定参数复合材料整体叶环振动特性计算方法,其特征在于,步骤四中,所述根据所述复合材料变异系数计算出复合材料整体叶环振动特性,还包括:
获取3组复合材料整体叶环振动特性的试验结果;
将复合材料整体叶环振动特性的计算结果与3组试验结果进行对比,若试验结果在计算结果的±3σ概率范围内,则验证得到具有不确定材料参数的复合材料整体叶环频率的概率分布。
6.根据权利要求5所述的不确定参数复合材料整体叶环振动特性计算方法,其特征在于,还包括步骤五、通过频率分散度绘制坎贝尔图进行所述复合材料整体叶环的避振设计。
7.根据权利要求6所述的不确定参数复合材料整体叶环振动特性计算方法,其特征在于,还包括步骤六、根据所述复合材料整体叶环振动特性的计算结果以及试验结果,形成复合材料整体叶环频率分散度标准。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310306832.5A CN116432311A (zh) | 2023-03-27 | 2023-03-27 | 一种不确定参数复合材料整体叶环振动特性计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310306832.5A CN116432311A (zh) | 2023-03-27 | 2023-03-27 | 一种不确定参数复合材料整体叶环振动特性计算方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN116432311A true CN116432311A (zh) | 2023-07-14 |
Family
ID=87080797
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310306832.5A Pending CN116432311A (zh) | 2023-03-27 | 2023-03-27 | 一种不确定参数复合材料整体叶环振动特性计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116432311A (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116663200A (zh) * | 2023-07-25 | 2023-08-29 | 中国航发四川燃气涡轮研究院 | 频率分散性可控的压气机整体叶盘叶片筛选方法及装置 |
-
2023
- 2023-03-27 CN CN202310306832.5A patent/CN116432311A/zh active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN116663200A (zh) * | 2023-07-25 | 2023-08-29 | 中国航发四川燃气涡轮研究院 | 频率分散性可控的压气机整体叶盘叶片筛选方法及装置 |
CN116663200B (zh) * | 2023-07-25 | 2023-10-20 | 中国航发四川燃气涡轮研究院 | 频率分散性可控的压气机整体叶盘叶片筛选方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Gao et al. | Probabilistic-based combined high and low cycle fatigue assessment for turbine blades using a substructure-based kriging surrogate model | |
CN116432311A (zh) | 一种不确定参数复合材料整体叶环振动特性计算方法 | |
CN110532723A (zh) | 一种基于egra的涡轮盘多失效模式可靠性优化方法 | |
US20210383035A1 (en) | Computer-implemented method for the probabilistic estimation of a probability of failure of a component, a data processing system, a computer program product and a computer-readable storage medium | |
CN110411850B (zh) | 一种高温合金涡轮叶片服役条件的评估方法 | |
CN103970965A (zh) | 燃气涡轮发动机加速寿命试验试车方法 | |
AU2019436655B2 (en) | Dynamic reliability evaluation method for coupling faults of middle trough of scraper conveyor | |
CN112989659B (zh) | 一种基于点权函数法的表面裂纹强度因子数据库建立方法 | |
CN108334704B (zh) | 基于自适应密度峰值-混合概率建模的结构损伤监测方法 | |
CN109086507B (zh) | 一种评价含复合点缺陷镍基合金变形能力的方法 | |
CN110411851A (zh) | 一种高温合金涡轮叶片服役损伤评价及蠕变寿命预测方法 | |
Gerberich et al. | Experimental and theoretical strain distributions for stationary and growing cracks | |
CN112100766A (zh) | 一种冷挤压强化孔结构的疲劳寿命预测方法 | |
Molaei et al. | Notched fatigue of additive manufactured metals under axial and multiaxial loadings, part II: Data correlations and life estimations | |
CN111366450A (zh) | 一种镍基高温合金涡轮叶片服役组织损伤评价方法 | |
Davidson et al. | Fatigue crack tip plastic strain in high‐strength aluminum alloys | |
CN113609584B (zh) | 重复使用火箭发动机推力室可靠性优化设计方法 | |
Balmes et al. | Constrained viscoelastic damping, test/analysis correlation on an aircraft engine | |
CN109522592A (zh) | 一种航空发动机叶片fod允修极限确定方法 | |
CN116541970B (zh) | 一种基于能量法的压气机叶片减振优化设计评估方法 | |
CN113008694A (zh) | 基于错配度的镍基高温合金涡轮叶片服役损伤评价方法 | |
CN115979804A (zh) | 一种增材制造材料寿命预测方法 | |
CN111521495A (zh) | 一种双相高强钢焊接粗晶区韧性优化方法 | |
CN105651603A (zh) | 一种tc18钛合金网篮组织断裂韧性预测方法 | |
Narasimhachary et al. | Life assessment of large gas turbine blades |
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 |