CN108052711A - 一种并联组合发动机安装结构的热力耦合分析方法 - Google Patents

一种并联组合发动机安装结构的热力耦合分析方法 Download PDF

Info

Publication number
CN108052711A
CN108052711A CN201711231186.1A CN201711231186A CN108052711A CN 108052711 A CN108052711 A CN 108052711A CN 201711231186 A CN201711231186 A CN 201711231186A CN 108052711 A CN108052711 A CN 108052711A
Authority
CN
China
Prior art keywords
analysis
installation structure
engine installation
parallel combination
mechanical method
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
Application number
CN201711231186.1A
Other languages
English (en)
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.)
Shenyang Aircraft Design and Research Institute Aviation Industry of China AVIC
Original Assignee
Shenyang Aircraft Design and Research Institute Aviation Industry of China AVIC
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 Shenyang Aircraft Design and Research Institute Aviation Industry of China AVIC filed Critical Shenyang Aircraft Design and Research Institute Aviation Industry of China AVIC
Priority to CN201711231186.1A priority Critical patent/CN108052711A/zh
Publication of CN108052711A publication Critical patent/CN108052711A/zh
Pending legal-status Critical Current

Links

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]

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)
  • Automation & Control Theory (AREA)
  • Aviation & Aerospace Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Testing Of Engines (AREA)

Abstract

本发明涉及飞机设计领域，特别涉及一种并联组合发动机安装结构的热力耦合分析方法，包括如下步骤：确定主承力构件，建立所述主承力构件的体单元细节有限元模型；对所述有限元模型进行静强度非线性分析，分析时对约束、接触、单元选择、载荷以及边界条件参数进行预定控制；对所述主承力构件进行热强度分析，分析时对材料和温度载荷参数进行预定控制；根据所述步骤二的静强度非线性分析结果和所述步骤三的热强度分析结果，对所述主承力构件进行热力耦合强度分析。本发明的并联组合发动机安装结构的热力耦合分析方法，与传统计算相比，真实模拟了并联组合发动机安装结构强度的非线性特性，因此计算分析结果更加准确。

Description

一种并联组合发动机安装结构的热力耦合分析方法
技术领域
本发明涉及飞机设计领域，特别涉及一种并联组合发动机安装结构的热力耦合分析方法。
背景技术
发动机为飞机内部最大的热源，发动机安装结构最高使用温度可达到500℃，同时发动机安装结构承受高量级的载荷，为保证飞行安全，需要进行发动机安装结构热力耦合强度分析。
发动机安装结构的使用温度与其距发动机的远近相关，靠近发动机处温度高，远离发动机处温度低，传统的工程计算方法无法计算由于温度梯度引起的热应力；发动机安装结构由不同材料的零件装配而成，由于不同材料的线膨胀系数不同，因此产生的热应力也不同，传统的工程计算方法无法计算由于材料线膨胀系数不同引起的热应力；发动机安装结构采用推力销与推力梁镶嵌式的装配关系，推力销与推力梁通过接触传力，传统的工程计算方法无法计算接触关系的边界非线性特征。
发明内容
本发明的目的是提供了一种并联组合发动机安装结构的热力耦合分析方法，以解决现有发动机安装结构热力耦合强度分析方法存在的至少一个问题。
本发明的技术方案是：
一种并联组合发动机安装结构的热力耦合分析方法，包括如下步骤：
步骤一、确定主承力构件，建立所述主承力构件的体单元细节有限元模型；
步骤二、对所述有限元模型进行静强度非线性分析，分析时对约束、接触、单元选择、载荷以及边界条件参数进行预定控制；
步骤三、对所述主承力构件进行热强度分析，分析时对材料和温度载荷参数进行预定控制；
步骤四、根据所述步骤二的静强度非线性分析结果和所述步骤三的热强度分析结果，对所述主承力构件进行热力耦合强度分析。
可选的，在所述步骤一中，所述主承力构件包括发动机的推力销、推力梁以及环形框。
可选的，在所述步骤一中，建立所述主承力构件的体单元细节有限元模型时，略去仅实现功能的零件，且对不影响零件静强度计算结果的细节部位的工艺孔进行填充，再去除倒角。
可选的，在所述步骤二中，约束、接触、单元选择、载荷以及边界条件参数的选取条件包括：
在所述有限元模型中，推力销与推力梁之间设定法向的接触和切向的摩擦接触，并通过接触算法模拟推力销与推力梁镶嵌的结构关系；
在所述有限元模型中，环形框与推力梁采用螺栓连接，在环形框与推力梁对接面处设置法向的接触和切向的摩擦接触。
可选的，在所述步骤二分析计算过程中，所述有限元模型中推力销、推力梁以及环形框均采用线性减缩积分C3D8R六面体单元建模。
可选的，在所述步骤二分析计算过程中，所述推力销在厚度方向具有四层以上的网格密度。
可选的，在所述步骤二分析计算过程中，在所述推力销端面处施加集中载荷，在环形框及推力梁上施加约束x、y、z方向的自由度。
可选的，在所述步骤二分析计算过程中，是按照载荷的10％为一步，分成10步进行计算。
可选的，在所述步骤三分析计算过程中，是设定常温为计算初始温度，来确定材料的线膨胀系数、弹性模量、泊松比随温度的变化关系。
可选的，在所述步骤三分析计算过程中，是在所述推力销、推力梁以及环形框的节点上分别施加温度载荷，以模拟温度场。
发明效果：
本发明的并联组合发动机安装结构的热力耦合分析方法，与传统计算相比，真实模拟了并联组合发动机安装结构强度的非线性特性，因此计算分析结果更加准确。
附图说明
图1是本发明并联组合发动机安装结构的热力耦合分析方法的流程图。
具体实施方式
为使本发明实施的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行更加详细的描述。在附图中，自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。所描述的实施例是本发明一部分实施例，而不是全部的实施例。下面通过参考附图描述的实施例是示例性的，旨在用于解释本发明，而不能理解为对本发明的限制。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。下面结合附图对本发明的实施例进行详细说明。
在本发明的描述中，需要理解的是，术语“中心”、“纵向”、“横向”、“前”、“后”、“左”、“右”、“竖直”、“水平”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本发明保护范围的限制。
下面结合附图1对本发明并联组合发动机安装结构的热力耦合分析方法做进一步详细说明。
本发明提供了一种并联组合发动机安装结构的热力耦合分析方法，包括如下步骤：
步骤一、确定主承力构件，建立所述主承力构件的体单元细节有限元模型。具体地，主承力构件包括发动机的推力销、推力梁以及环形框，体单元细节有限元模型是关于推力销、推力梁以及环形框的有限元模型。进一步，建立主承力构件的体单元细节有限元模型时，略去仅实现功能的零件，且对不影响零件静强度计算结果的细节部位的工艺孔进行填充，再去除倒角。
步骤二、对所述有限元模型进行静强度非线性分析，分析时对约束、接触、单元选择、载荷以及边界条件参数进行预定控制，提交运算并调试。
具体地，约束、接触、单元选择、载荷以及边界条件参数的选取条件包括：
在有限元模型中，推力销与推力梁之间设定法向的接触和切向的摩擦接触，并通过接触算法模拟推力销与推力梁镶嵌的结构关系；
在有限元模型中，环形框与推力梁采用螺栓连接，在环形框与推力梁对接面处设置法向的接触和切向的摩擦接触。
进一步，在步骤二分析计算过程中，有限元模型中推力销、推力梁以及环形框均采用线性减缩积分C3D8R六面体单元建模。
进一步，在步骤二分析计算过程中，所述推力销在厚度方向具有四层以上的网格密度。
进一步，在步骤二分析计算过程中，在推力销端面处施加集中载荷，在环形框及推力梁上施加约束x、y、z方向的自由度。
进一步，在步骤二分析计算过程中，是按照载荷的10％为一步，分成10步进行计算。
步骤三、对主承力构件进行热强度分析，分析时对材料和温度载荷参数进行预定控制。具体地，在步骤三分析计算过程中，是设定常温为计算初始温度，来确定材料的线膨胀系数、弹性模量、泊松比随温度的变化关系。
进一步，在步骤三分析计算过程中，是在推力销、推力梁以及环形框的节点上分别施加温度载荷，以模拟温度场。
步骤四、根据步骤二的静强度非线性分析结果和步骤三的热强度分析结果，提交ABAQUS/Standard模块对主承力构件进行热力耦合强度分析。
本发明的并联组合发动机安装结构的热力耦合分析方法，与传统计算相比，真实模拟了并联组合发动机安装结构强度的非线性特性，因此计算分析结果更加准确。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以所述权利要求的保护范围为准。

Claims (10)

1.一种并联组合发动机安装结构的热力耦合分析方法，其特征在于，包括如下步骤：
步骤一、确定主承力构件，建立所述主承力构件的体单元细节有限元模型；
步骤二、对所述有限元模型进行静强度非线性分析，分析时对约束、接触、单元选择、载荷以及边界条件参数进行预定控制；
步骤三、对所述主承力构件进行热强度分析，分析时对材料和温度载荷参数进行预定控制；
步骤四、根据所述步骤二的静强度非线性分析结果和所述步骤三的热强度分析结果，对所述主承力构件进行热力耦合强度分析。
2.根据权利要求1所述的并联组合发动机安装结构的热力耦合分析方法，其特征在于，在所述步骤一中，所述主承力构件包括发动机的推力销、推力梁以及环形框。
3.根据权利要求2所述的并联组合发动机安装结构的热力耦合分析方法，其特征在于，在所述步骤一中，建立所述主承力构件的体单元细节有限元模型时，略去仅实现功能的零件，且对不影响零件静强度计算结果的细节部位的工艺孔进行填充，再去除倒角。
4.根据权利要求2或3所述的并联组合发动机安装结构的热力耦合分析方法，其特征在于，在所述步骤二中，约束、接触、单元选择、载荷以及边界条件参数的选取条件包括：
在所述有限元模型中，推力销与推力梁之间设定法向的接触和切向的摩擦接触，并通过接触算法模拟推力销与推力梁镶嵌的结构关系；
在所述有限元模型中，环形框与推力梁采用螺栓连接，在环形框与推力梁对接面处设置法向的接触和切向的摩擦接触。
5.根据权利要求4所述的并联组合发动机安装结构的热力耦合分析方法，其特征在于，在所述步骤二分析计算过程中，所述有限元模型中推力销、推力梁以及环形框均采用线性减缩积分C3D8R六面体单元建模。
6.根据权利要求5所述的并联组合发动机安装结构的热力耦合分析方法，其特征在于，在所述步骤二分析计算过程中，所述推力销在厚度方向具有四层以上的网格密度。
7.根据权利要求6所述的并联组合发动机安装结构的热力耦合分析方法，其特征在于，在所述步骤二分析计算过程中，在所述推力销端面处施加集中载荷，在环形框及推力梁上施加约束x、y、z方向的自由度。
8.根据权利要求7所述的并联组合发动机安装结构的热力耦合分析方法，其特征在于，在所述步骤二分析计算过程中，是按照载荷的10％为一步，分成10步进行计算。
9.根据权利要求4所述的并联组合发动机安装结构的热力耦合分析方法，其特征在于，在所述步骤三分析计算过程中，是设定常温为计算初始温度，来确定材料的线膨胀系数、弹性模量、泊松比随温度的变化关系。
10.根据权利要求9所述的并联组合发动机安装结构的热力耦合分析方法，其特征在于，在所述步骤三分析计算过程中，是在所述推力销、推力梁以及环形框的节点上分别施加温度载荷，以模拟温度场。
CN201711231186.1A 2017-11-29 2017-11-29 一种并联组合发动机安装结构的热力耦合分析方法 Pending CN108052711A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711231186.1A CN108052711A (zh) 2017-11-29 2017-11-29 一种并联组合发动机安装结构的热力耦合分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711231186.1A CN108052711A (zh) 2017-11-29 2017-11-29 一种并联组合发动机安装结构的热力耦合分析方法

Publications (1)

Publication Number Publication Date
CN108052711A true CN108052711A (zh) 2018-05-18

Family

ID=62121237

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711231186.1A Pending CN108052711A (zh) 2017-11-29 2017-11-29 一种并联组合发动机安装结构的热力耦合分析方法

Country Status (1)

Country Link
CN (1) CN108052711A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108959822A (zh) * 2018-08-01 2018-12-07 中国航空工业集团公司沈阳飞机设计研究所 一种减小热结构内部热应力作用的设计方法
CN109540962A (zh) * 2018-11-30 2019-03-29 中国航空工业集团公司沈阳飞机设计研究所 一种隔热结构的隔热效能表征方法
CN111898209A (zh) * 2020-08-05 2020-11-06 中铁磁浮科技（成都）有限公司 一种对制动夹钳的静强度特性分析方法

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6301970B1 (en) * 1998-08-06 2001-10-16 The United States Of America The Secretary Of The Navy Cumulative damage model for structural analysis of filed polymeric materials
CN103984803A (zh) * 2014-04-17 2014-08-13 中国航空工业集团公司沈阳飞机设计研究所 一种有限元载荷施加方法
CN104408244A (zh) * 2014-11-19 2015-03-11 中国航空工业集团公司沈阳飞机设计研究所 一种适用于座舱盖结构的静、热强度耦合分析方法

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6301970B1 (en) * 1998-08-06 2001-10-16 The United States Of America The Secretary Of The Navy Cumulative damage model for structural analysis of filed polymeric materials
CN103984803A (zh) * 2014-04-17 2014-08-13 中国航空工业集团公司沈阳飞机设计研究所 一种有限元载荷施加方法
CN104408244A (zh) * 2014-11-19 2015-03-11 中国航空工业集团公司沈阳飞机设计研究所 一种适用于座舱盖结构的静、热强度耦合分析方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
罗金亮 等: "飞机发动机固定装置的载荷计算方法", 《洪都科技》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108959822A (zh) * 2018-08-01 2018-12-07 中国航空工业集团公司沈阳飞机设计研究所 一种减小热结构内部热应力作用的设计方法
CN109540962A (zh) * 2018-11-30 2019-03-29 中国航空工业集团公司沈阳飞机设计研究所 一种隔热结构的隔热效能表征方法
CN109540962B (zh) * 2018-11-30 2021-07-09 中国航空工业集团公司沈阳飞机设计研究所 一种隔热结构的隔热效能表征方法
CN111898209A (zh) * 2020-08-05 2020-11-06 中铁磁浮科技（成都）有限公司 一种对制动夹钳的静强度特性分析方法

Similar Documents

Publication Publication Date Title
CN108052711A (zh) 一种并联组合发动机安装结构的热力耦合分析方法
JP6862148B2 (ja) 可撓体のシム合わせのためのシステム及び方法
Rezaei et al. Modal-based damage identification for the nonlinear model of modern wind turbine blade
Benabou et al. A thermo-mechanical cohesive zone model for solder joint lifetime prediction
Botez et al. Laminar flow control on a research wing project presentation covering a three year period
CN103559390B (zh) 一种基于平均失效指数的复合材料π形胶接连接结构拉伸强度预测方法
CN107145676A (zh) 一种复杂载荷条件下的结构层级式有限元仿真分析方法
CN109388820A (zh) 确定不规则温循剖面下焊点损伤累积完整法则的有限元仿真方法
CN103324787B (zh) 批产卫星热设计验证方法
Keye et al. Validation of wing deformation simulations for the NASA CRM model using fluid-structure interaction computations
Sepe et al. Crack growth behavior of welded stiffened panel
Lin et al. Reliability analysis of aero-engine blades considering nonlinear strength degeneration
Liu et al. Crack growth analysis and fatigue life estimation in the piston rod of a Kaplan hydro turbine
KR101428720B1 (ko) 비행자료를 이용한 응력 스펙트럼 생성 방법
Yun et al. An improved crack tracking algorithm with self‐correction ability of the crack path and its application in a continuum damage model
Hahn et al. Design study of dovetail geometries of turbine blades using abaqus and isight
CN108090260A (zh) 一种接头约束载荷的分析方法
CN108562447A (zh) 一种民用飞机刹车系统综合环境试验箱的设计方法
Cartieri et al. Using RANS computations to calculate support interference effects on the Common Research Model
CN115186522A (zh) 一种考虑蠕变疲劳作用的飞行器头锥结构失效分析方法
Vrchota et al. Using wing model deformation for improvement of CFD results of ESWIRP project
WO2019078851A1 (en) SYSTEMS AND METHODS FOR ASSESSING STRUCTURAL HEALTH
CN111428398B (zh) 一种C/SiC舵面热强度计算方法
Elliott et al. Model calibration efforts for the international space station's solar array mast
Nielsen Adjoint-based aerodynamic design of complex aerospace configurations

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20180518

RJ01 Rejection of invention patent application after publication