CN116415348A - 一种多孔层合结构构件的刚度等效模型构建方法 - Google Patents

一种多孔层合结构构件的刚度等效模型构建方法 Download PDF

Info

Publication number
CN116415348A
CN116415348A CN202310162626.1A CN202310162626A CN116415348A CN 116415348 A CN116415348 A CN 116415348A CN 202310162626 A CN202310162626 A CN 202310162626A CN 116415348 A CN116415348 A CN 116415348A
Authority
CN
China
Prior art keywords
unit cell
laminate
model
constructing
equivalent
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
CN202310162626.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.)
AECC Sichuan Gas Turbine Research Institute
Original Assignee
AECC Sichuan Gas Turbine 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 Sichuan Gas Turbine Research Institute filed Critical AECC Sichuan Gas Turbine Research Institute
Priority to CN202310162626.1A priority Critical patent/CN116415348A/zh
Publication of CN116415348A publication Critical patent/CN116415348A/zh
Pending legal-status Critical Current

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/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
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08—Thermal analysis or thermal optimisation
    • 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)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明提供了一种多孔层合结构构件的刚度等效模型构建方法，包括S1、获取层板上细孔的单胞特征尺寸参数；S2、构建层板的单胞实体模型；S3、构建层板的单胞实体有限元分析模型；S4、对单胞实体有限元分析模型进行前处理；S5、计算不同温度下层板的刚度等效模量；S6、重复步骤S1～S5，计算多孔层合结构构件中所有层板在不同温度下各方向的刚度等效模量；S7、对步骤S6中计算的刚度等效模量进行一体化边界处理，构建多孔层合结构构件的刚度等效模型。本发明的方法可以将多孔层合结构构件简化成无孔层合结构模型，可以很好地体现多孔层合结构构件的力学性质，同时能大幅提高计算效率和计算精度，具有广泛的工程应用前景。

Description

一种多孔层合结构构件的刚度等效模型构建方法
技术领域
本发明属于航空发动机领域，涉及带孔构件刚度等效模型设计技术，具体涉及一种多孔层合结构构件的刚度等效模型构建方法。
背景技术
航空发动机中部分构件，例如涡轮叶片、火焰筒、隔热屏等含有大量细孔，孔径不超过1mm，且孔总面积占构件总面积50％以上。对这些构件进行常规有限元分析时，若保留这些细孔的结构特征，构建构件计算模型时使得网格划分工作十分困难且非常耗时，且网格数量十分巨大，达到数百万量级，从而导致了计算效率低下、难以模拟真实使用场景；若不考虑细孔结构，虽然简化了构件计算模型，但是但是计算精度不能满足工程要求，不能准确的反映出构件在工作状态下的真实情况。
鉴于此，有必要对含有细孔构件的计算模型构建方法进行改进。
发明内容
为解决对航空发动机中带有大量细孔构件进行有限元分析时，采用现有方法构建的计算模型，难以模拟真实使用场景、计算结果准确度低、可信度低等问题，本发明公开了一种多孔层合结构构件的刚度等效模型构建方法，该方法是一种“刚度等效”方法，可以将多孔层合结构构件简化成无孔层合结构模型，可以很好地体现多孔层合结构构件的力学性质，同时能大幅提高计算效率和计算精度，具有广泛的工程应用前景。
实现发明目的的技术方案如下：一种多孔层合结构构件的刚度等效模型构建方法，多孔层合结构构件包括多个层板，多孔层合结构构件的刚度等效模型构建方法包括以下步骤：
S1、获取层板上细孔的单胞特征尺寸参数；
S2、单胞特征尺寸参数，构建层板的单胞实体模型；
S3、对单胞实体模型分网，构建层板的单胞实体有限元分析模型；
S4、对单胞实体有限元分析模型施加载性能参数、温度、约束进行前处理；
S5、计算不同温度下层板的刚度等效模量；
S6、重复步骤S1～S5，计算多孔层合结构构件中所有层板在不同温度下各方向的刚度等效模量；
S7、对步骤S6中计算的刚度等效模量进行一体化边界处理，构建多孔层合结构构件的刚度等效模型。
进一步地，步骤S1中，层板上细孔的单胞特征尺寸参数的获取方法，包括：
S11、构建层板的结构模型，依据细孔孔径确定该层板的细孔分布规律；
S12、基于细孔分布规律，获取层板上细孔的单胞特征尺寸参数。
更进一步地，步骤S12中，单胞特征尺寸参数包括细孔半径、相邻两个细孔孔心间横向距离、相邻两个细孔孔心间纵向距离、层板厚度；其中，相邻两个细孔是指孔径相同的两个细孔。
进一步地，步骤S3中，单胞实体模型分网时，网格类型为四面体网格单元或者六面体网格单元，且分网时对单胞实体模型中孔边位置进行加密划分。
进一步地，步骤S2中，若层板中有≥2个规格的单胞实体模型，则在步骤S3中对各规格的单胞实体模型分别分网建构单胞实体有限元分析模型，再组合形成该层板的单胞实体模型。
进一步地，步骤S4中，所述性能参数为材料性能参数，所述约束为边界约束和载荷约束。
进一步地，步骤S5中，不同温度下层板的刚度等效模量的计算方法，包括：
S51、在不同温度下，提取单胞实体有限元分析模型中各网格单元在各方向的应力及应变；
S52、计算每一个网格单元各方向的应力及应变与网格单元体积的乘积，求和后除网格单元总体积，获得单胞实体有限元分析模型的等效应力及等效应变；
S53、基于等效应力及等效应变，计算层板的刚度等效模量。
与现有技术相比，本发明的有益效果是：本发明提供的多孔层合结构构件的刚度等效模型构建方法，是一种“刚度等效”的方法，该方法简化了带有大量细孔的多孔层合结构构件的刚度等效模型，降低分析难度，提升计算效率，提高分析结果的准确性及可靠性，提供了理论支撑，为航空发动机中多孔层合结构的设计与分析提供基础支持，具有工程应用前景。
同时，采用本发明的方法，可比原模型降低有限元分析网格数量70％以上，节省了80％以上的计算分析时间；分析时由于施加了刚度等效模量的材料性能参数，进而能够很好地体现原多孔层合结构构件的力学性质，提高了计算精度。经验证，本发明的方法在某型发动机某多孔层合结构构件的有限元分析上已经成功应用，并在此基础上建立了完整的分析流程。
附图说明
为了更清楚地说明本发明实施例技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍。
图1为本发明多孔层合结构构件的刚度等效模型构建方法的流程图；
图2为具体实施方式中层板上细孔的分布示意图；
图3为具体实施方式中多孔层合结构构件一个层板的单胞实体模型网格划分后形成的单胞实体有限元分析模型的示意图。
具体实施方式
下面结合具体实施例来进一步描述本发明，本发明的优点和特点将会随着描述而更为清楚。但这些实施例仅是范例性的，并不对本发明的范围构成任何限制。本领域技术人员应该理解的是，在不偏离本发明的精神和范围下可以对本发明技术方案的细节和形式进行修改或替换，但这些修改和替换均落入本发明的保护范围内。
本具体实施方式针对包括多个层板的多孔层合结构构件构建刚度等效模型，用以对航空发动机中一些带有大量细孔的多层构件的分析提供基础支持，同时需要特别说明的是多个层板中具有相同细孔分布规律的层板可以视为同一个层板。
参见图1所示，本具体实施方式公开了一种多孔层合结构构件的刚度等效模型构建方法，包括以下步骤：
S1、获取层板上细孔的单胞特征尺寸参数；
S2、单胞特征尺寸参数，构建层板的单胞实体模型；
S3、对单胞实体模型分网，构建层板的单胞实体有限元分析模型；
S4、对单胞实体有限元分析模型施加载性能参数、温度、约束进行前处理；
S5、计算不同温度下层板的刚度等效模量；
S6、重复步骤S1～S5，计算多孔层合结构构件中所有层板在不同温度下各方向的刚度等效模量；
S7、对步骤S6中计算的刚度等效模量进行一体化边界处理，构建多孔层合结构构件的刚度等效模型。
以下通过具体的实施例对上述刚度等效模型构建方法中各步骤进行说明：
在一个实施例中，步骤S1中，层板上细孔的单胞特征尺寸参数的获取方法，包括：
S11、构建层板的结构模型，依据细孔孔径确定该层板的细孔分布规律；
S12、基于细孔分布规律，获取层板上细孔的单胞特征尺寸参数；
其中，单胞特征尺寸参数包括细孔半径、相邻两个细孔孔心间横向距离、相邻两个细孔孔心间纵向距离、层板厚度；其中，相邻两个细孔是指孔径相同的两个细孔。
在此需要说明的是，参见图2所示，同一个层板上，可能含有一种或多种尺寸的细孔，对这些种类的细孔，要分别计算单胞特征尺寸参数。
在一个实施例中，步骤S2中，单胞实体模型是用单胞特征尺寸参数，采用CAD、UG等几何建模软件中构建的。
需要说明的是，当层板中有≥2个规格的单胞实体模型，则在步骤S3中对各规格的单胞实体模型分别分网建构单胞实体有限元分析模型，再组合形成该层板的单胞实体模型。
在一个实施例中，步骤S3中，层板的单胞实体有限元分析模型是通过将单胞实体模型导入ANSA、Hypermesh等分网软件进行有限元网格划分，采用10节点四面体单元或者20节点六面体单元进行网格划分。
需要说明的是：为了保证求解的效率和计算结果的准确性，单胞实体模型网格划分时不应太多也不能太少，且在层板厚度方向上网格至少保证2层，同时对单胞实体模型中孔边位置进行加密，通常来说孔边位置的网格单元数量一般不低于20，本具体实施方式中，参见图3所示为多孔层合结构构件其中一个层板的单胞实体模型进行有限元网格划分后形成的单胞实体有限元分析模型的示意图。
在一个实施例中，需要对构建的单胞实体有限元分析模型进行预处理，即步骤S4中，对单胞实体有限元分析模型施加载性能参数、温度、约束进行前处理，其中所述性能参数为材料性能参数，所述约束为边界约束和载荷约束。
需要说明的是，由于在不同温度下，材料性能参数是不同的，因此在对单胞实体有限元分析模型进行前处理时，还确保材料性能参数与温度的对应关系。
在施加边界约束和载荷约束时，由于相邻两个细孔孔心间横向距离、相邻两个细孔孔心间纵向距离不一定相等，因此单胞实体有限元分析模型在xy平面内的内力学性能可能呈现出各项异性，即Ex≠Ey，就单胞细观结构而言，由于其具有三个正交的对称平面，该等效材料是正交各向异性材料，独立的工程材料变量有Exx、Eyy、Ezz、Gxy、Gxz、Gyz，其中Exx为六面体网格单元x方向的弹性模量，Eyy为六面体网格单元y方向的弹性模量，Ezz为六面体网格单元z方向的弹性模量，Gxy为六面体网格单元xy平面的剪切模量，Gxz为六面体网格单元xz平面的剪切模量，Gyz为六面体网格单元yz平面的剪切模量。
在求解Exx时，单胞约束边界条件为约束单胞x方向的一侧，另一侧施加0.1mm的强迫位移；y、z方向两侧面一侧的约束为0，另一侧保持变形一致，施加位移耦合约束。Ezz、Eyy、Gxy、Gxz、Gyz的约束边界与Exx的约束边界类似。需要说明的是，不同温度下的材料性能参数需要分别计算，模型一样，材料性能不一样，分别进行求解。
在一个实施例中，步骤S5中，不同温度下层板的刚度等效模量的计算方法，包括：
S51、在不同温度下，提取单胞实体有限元分析模型中各网格单元在各方向的应力σij(x)及应变εij(x)：
S52、计算每一个网格单元各方向的应力及应变与网格单元体积的乘积，求和后除网格单元总体积，获得单胞实体有限元分析模型的等效应力
Figure BDA0004094810700000071
及等效应变/>
Figure BDA0004094810700000072
S53、基于等效应力及等效应变，计算层板的刚度等效模量，即
Figure BDA0004094810700000073
其中，σ(x)为材料细观应力，ε(x)为材料细观应变；ΩX为x处体积单元所占区域，v为网格单元体积，x为细观坐标，
Figure BDA0004094810700000074
为多孔层合结构构件刚度等效模型的刚度等效模量。
通过步骤S1～S5，能够计算出多孔层合结构构件中所有层板在不同温度下各方向的刚度等效模量。在对整个多孔层合结构构件分析时，使用其各层板的刚度等效模量，对其进行一体化边界处理，即可构建出多孔层合结构构件的刚度等效模型。
本具体实施方式提供的多孔层合结构构件的刚度等效模型构建方法，是一种“刚度等效”的方法，该方法简化了带有大量细孔的多孔层合结构构件的刚度等效模型，降低分析难度，提升计算效率，提高分析结果的准确性及可靠性，提供了理论支撑，为航空发动机中多孔层合结构的设计与分析提供基础支持，具有工程应用前景。
同时，采用本发明的方法，可比原模型降低有限元分析网格数量70％以上，节省了80％以上的计算分析时间；分析时由于施加了刚度等效模量的材料性能参数，进而能够很好地体现原多孔层合结构构件的力学性质，提高了计算精度。经验证，本发明的方法在某型发动机某多孔层合结构构件的有限元分析上已经成功应用，并在此基础上建立了完整的分析流程。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。
此外，应当理解，虽然本说明书按照实施方式加以描述，但并非每个实施方式仅包含一个独立的技术方案，说明书的这种叙述方式仅仅是为清楚起见，本领域技术人员应当将说明书作为一个整体，各实施例中的技术方案也可以经适当组合，形成本领域技术人员可以理解的其他实施方式。

Claims (7)

1.一种多孔层合结构构件的刚度等效模型构建方法，其特征在于，多孔层合结构构件包括多个层板，多孔层合结构构件的刚度等效模型构建方法包括以下步骤：
S1、获取层板上细孔的单胞特征尺寸参数；
S2、单胞特征尺寸参数，构建层板的单胞实体模型；
S3、对单胞实体模型分网，构建层板的单胞实体有限元分析模型；
S4、对单胞实体有限元分析模型施加载性能参数、温度、约束进行前处理；
S5、计算不同温度下层板的刚度等效模量；
S6、重复步骤S1～S5，计算多孔层合结构构件中所有层板在不同温度下各方向的刚度等效模量；
S7、对步骤S6中计算的刚度等效模量进行一体化边界处理，构建多孔层合结构构件的刚度等效模型。
2.根据权利要求1所述的刚度等效模型构建方法，其特征在于，步骤S1中，层板上细孔的单胞特征尺寸参数的获取方法，包括：
S11、构建层板的结构模型，依据细孔孔径确定该层板的细孔分布规律；
S12、基于细孔分布规律，获取层板上细孔的单胞特征尺寸参数。
3.根据权利要求2所述的刚度等效模型构建方法，其特征在于，步骤S12中，单胞特征尺寸参数包括细孔半径、相邻两个细孔孔心间横向距离、相邻两个细孔孔心间纵向距离、层板厚度；其中，相邻两个细孔是指孔径相同的两个细孔。
4.根据权利要求1所述的刚度等效模型构建方法，其特征在于，步骤S3中，单胞实体模型分网时，网格类型为四面体网格单元或者六面体网格单元，且分网时对单胞实体模型中孔边位置进行加密划分。
5.根据权利要求1所述的刚度等效模型构建方法，其特征在于，步骤S2中，若层板中有≥2个规格的单胞实体模型，则在步骤S3中对各规格的单胞实体模型分别分网建构单胞实体有限元分析模型，再组合形成该层板的单胞实体模型。
6.根据权利要求1所述的刚度等效模型构建方法，其特征在于，步骤S4中，所述性能参数为材料性能参数，所述约束为边界约束和载荷约束。
7.根据权利要求1所述的刚度等效模型构建方法，其特征在于，步骤S5中，不同温度下层板的刚度等效模量的计算方法，包括：
S51、在不同温度下，提取单胞实体有限元分析模型中各网格单元在各方向的应力及应变；
S52、计算每一个网格单元各方向的应力及应变与网格单元体积的乘积，求和后除网格单元总体积，获得单胞实体有限元分析模型的等效应力及等效应变；
S53、基于等效应力及等效应变，计算层板的刚度等效模量。
CN202310162626.1A 2023-02-24 2023-02-24 一种多孔层合结构构件的刚度等效模型构建方法 Pending CN116415348A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310162626.1A CN116415348A (zh) 2023-02-24 2023-02-24 一种多孔层合结构构件的刚度等效模型构建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310162626.1A CN116415348A (zh) 2023-02-24 2023-02-24 一种多孔层合结构构件的刚度等效模型构建方法

Publications (1)

Publication Number Publication Date
CN116415348A true CN116415348A (zh) 2023-07-11

Family

ID=87057343

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310162626.1A Pending CN116415348A (zh) 2023-02-24 2023-02-24 一种多孔层合结构构件的刚度等效模型构建方法

Country Status (1)

Country Link
CN (1) CN116415348A (zh)

Similar Documents

Publication Publication Date Title
Nguyen et al. Heuristic optimization method for cellular structure design of light weight components
JP2022521907A (ja) ハイブリッド繊維複合材料の板巻きシェル構造に対する高速協調最適化方法
CN105488623B (zh) 一种飞机质量特性数据处理方法
Hadi et al. Design configurations and creation of lattice structures for metallic additive manufacturing
CN103020358A (zh) 针对机械装备的可适应性动态设计平台构建方法
CN112765732B (zh) 一种基于选区激光熔化工艺的航空叶片拓扑优化设计方法
CN111177861B (zh) 适用于增材制造成形技术的常平环结构轻量化设计方法
Petrolo et al. Best theory diagrams for multilayered structures via shell finite elements
Mazumder et al. A structured method to generate conformal FE mesh for realistic textile composite micro-geometry
Wu et al. Effects of higher-order global–local shear deformations on bending, vibration and buckling of multilayered plates
Yuan et al. Contact stress analysis and fatigue life prediction of turbine disc–blade attachment with fir‐tree tenon structure
CN113902872A (zh) 非结构基质网格与裂缝连接性的检测方法、装置及介质
CN112949004B (zh) 风力发电机组轴承座的轻量化设计方法及其轴承座
CN116415348A (zh) 一种多孔层合结构构件的刚度等效模型构建方法
CN112487673A (zh) 基于机床工作状态的关键主机组件结构优化设计方法
CN113806951A (zh) 一种基于半边数据结构的自然邻近点搜索的弹性仿真方法
CN108491654A (zh) 一种三维实体结构拓扑优化方法及系统
CN112307659A (zh) 一种双层柱形耐压壳极限承载能力数值的计算方法
CN114912306B (zh) 一种2.5d机织复合材料平板结构多尺度模型的建模方法
CN113076671B (zh) 有限元仿真分析中的阻尼器简化方法和计算机设备
Yan et al. Concurrent material and structural optimization of hollow plate with truss-like material
CN112989524B (zh) 一种球轴承滚动体建模方法、系统及存储介质
Boccini et al. Towards structural topology optimization of rotor blisks
CN112380624B (zh) 一种直升机座舱骨架刚度优化方法
CN108595724A (zh) 复合材料回转件设计方法

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