CN109858054B - 一种飞行器前缘结构有限元模型建立方法 - Google Patents
一种飞行器前缘结构有限元模型建立方法 Download PDFInfo
- Publication number
- CN109858054B CN109858054B CN201711246410.4A CN201711246410A CN109858054B CN 109858054 B CN109858054 B CN 109858054B CN 201711246410 A CN201711246410 A CN 201711246410A CN 109858054 B CN109858054 B CN 109858054B
- Authority
- CN
- China
- Prior art keywords
- shell
- aircraft
- unit
- front edge
- finite element
- 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.)
- Active
Links
Images
Classifications
-
- 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
Abstract
本发明公开了一种飞行器前缘结构有限元模型建立方法,在飞行器前缘的对称面建立壳单元,选取前缘结构中心对称面的第一个圆孔区域进行壳单元网格划分,对划分的壳单元网格进行右手螺旋定律检查,使其严格满足右手螺旋定律,单孔区域壳单元沿飞行器前缘的对称面平移8次建立完整的壳单元,完整的壳单元沿飞行器前缘母线以30等分拉伸生成单层体单元,这样,能够方便地对拉伸和镜像生成的全结构体单元孔内表面施加对流边界条件以及前缘外表面的热流施加。生成的体单元再通过镜像生成前缘全结构有限元模型。本发明在国内首次提出并成功实现了一种飞行器前缘结构有限元模型建立方法,该方法物理概念明晰、操作方法简单、实施效果良好。
Description
技术领域
本发明属于飞行器前缘结构有限元建模领域,具体涉及一种飞行器前缘结构有限元模型建立方法。
背景技术
航空航天飞行器进行高超音速飞行的情况下,其表面受到高速气流的气动加热和对流换热联合作用,飞行器前缘最高表面温度甚至超高1800℃,会对飞行器结构产生不良影响,因此必须对飞行器前缘结构进行数值仿真,通过数值仿真确定出强迫对流换热系数。
由于飞行器前缘承受气动加热,气动加热作为施加热流与结构表面温度存在耦合关系,因此,建立有限元模型时,必须能够快速提取飞行器前缘体单元表面编号和飞行器前缘体内部圆孔体单元表面编号。
发明内容
本发明的目的是:提供一种飞行器前缘结构有限元模型建立方法。
本发明的技术方案如下:
一种飞行器前缘结构有限元模型建立方法,包括以下步骤:
步骤一、选取结构中心面的第一个圆孔区域进行平面网格划分形成壳单元;
步骤二、对当前壳单元沿着面法线方向拉伸为单层体单元,拉伸控制参数取为1,对体单元孔内表面施加对流边界条件,得知单元表面有两种编号,把两种编号改变为一种编号,则在全结构体单元孔内表面施加对流边界条件才能方便地选择;
步骤三、修改壳单元编号,使之右手螺旋定律,并删除单层体单元;
步骤四、对单孔区域的壳单元平移形成对称面上的全部壳单元;
步骤五、将对称面上的全部壳单元作为基单元沿两条曲线滑动产生新单元;
步骤六、镜像生成另一部分体单元网格;
步骤七、删除壳单元网格,合并重合点;
步骤八、在飞行器前缘外表面施加热流载荷;
步骤九、在飞行器前缘的九个孔的内表面施加强迫对流边界条件;
进一步地,所述步骤五中产生新单元的拉伸控制参数取为30。
进一步地,步骤一中,对所述划分形成的壳单元网格进行右手螺旋定律检查,使其严格满足右手螺旋定律。
进一步地,步骤四中,所述的单孔区域壳单元沿飞行器前缘的对称面平移8次建立完整的壳单元。
附图说明
图1是飞行器前缘结构示意图;
图2是单孔壳单元连接图;
图3是完整壳单元有限元网格;
图4是沿飞行器前缘母线拉伸生成的体单元有限元网格;
图5是沿通过镜像生成的完整体单元有限元网格;
图6是飞行器前缘外表面施加的热流载荷;
图7是在飞行器前缘的九个孔的内表面施加强迫对流边界条件;
图8是飞行器前缘某时刻的温度云图显示。
具体实施方式
发明原理:
一个完整的有限元模型包含网格划分、载荷、约束的施加以及边界条件施加,对于结构比较复杂的模型,根据载荷、约束、边界条件的不同,在建立模型时需要通盘考虑。
对于给定的飞行器前缘,其母线、施加对流的内部通道尺寸和数量已经完全确定,有限元数值分析时要求在飞行器前缘外表面施加热流,在飞行器前缘内部圆孔表面施加对流边界条件,由于结构形式的复杂性,在有限元模型建立时必须采取一些技巧才能快速、准确地对家里的模型施加载荷与边界条件。
下面结合附图和实施例对本发明做进一步说明:
步骤一、对飞行器前缘在对称中面实施切分,得到结构中心面(图1),对结构中心面的第一个圆孔区域进行平面网格划分,划分结果见图2,在有限元建模软件平台自动形成的单元连接方式见表1;
表1Quad4网格连接方式
步骤二、表1满足右手螺旋定律模型建立规则,对当前壳单元沿着面法线方向拉伸为单层体单元(拉伸控制参数取为1),对体单元孔内表面施加对流边界条件,得知单元表面有两种编号,必须把两种编号改变为统一编号,则在全结构体单元孔内表面施加对流边界条件才能方便地选择;
步骤三、手动修改表1的壳单元CQUAD4编号为表2,由于单层体单元的建立仅仅是为了显示圆孔的单元内表面的编号,圆孔的单元内表面的编号确定后,删除单层体单元;
步骤四、对图2在当前坐标系进行平移复制,形成对称面上的全部壳单元(图3);
表2修改后的Quad4网格连接方式
步骤五、控制参数取为30,将对称面上的全部壳单元作为基单元沿两条曲线滑动生成体单元,如图4所示。
步骤六、镜像(M i rror)生成另一部分体单元网格,如图5所示;
步骤七、删除壳单元网格,合并重合点;
步骤八、在飞行器前缘外表面施加热流载荷,如图6所示;
步骤九、在飞行器前缘的九个孔的内表面施加强迫对流边界条件,如图7所示。
Claims (4)
1.一种飞行器前缘结构有限元模型建立方法,其特征在于,包括以下步骤:
步骤一、选取结构中心面的第一个圆孔区域进行平面网格划分形成壳单元;
步骤二、对当前壳单元沿着面法线方向拉伸为单层体单元,拉伸控制参数取为1,对体单元孔内表面施加对流边界条件,得知单元表面有两种编号,把两种编号改变为一种编号,则在全结构体单元孔内表面施加对流边界条件才能方便地选择;
步骤三、修改壳单元编号,使之满足右手螺旋定律,并删除单层体单元;
步骤四、对单孔区域的壳单元平移形成对称面上的全部壳单元;
步骤五、将对称面上的全部壳单元作为基单元沿两条曲线滑动产生新单元;
步骤六、镜像生成另一部分体单元网格;
步骤七、删除壳单元网格,合并重合点;
步骤八、在飞行器前缘外表面施加热流载荷;
步骤九、在飞行器前缘的九个孔的内表面施加强迫对流边界条件。
2.根据权利要求1所述的一种飞行器前缘结构有限元模型建立方法,其特征在于,所述步骤五中产生新单元的拉伸控制参数取为30。
3.根据权利要求1所述的一种飞行器前缘结构有限元模型建立方法,其特征在于,步骤一中,对所述划分形成的壳单元网格进行右手螺旋定律检查,使其严格满足右手螺旋定律。
4.根据权利要求1所述的一种飞行器前缘结构有限元模型建立方法,其特征在于,步骤四中,所述的单孔区域壳单元沿飞行器前缘的对称面平移8次建立完整的壳单元。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711246410.4A CN109858054B (zh) | 2017-11-30 | 2017-11-30 | 一种飞行器前缘结构有限元模型建立方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711246410.4A CN109858054B (zh) | 2017-11-30 | 2017-11-30 | 一种飞行器前缘结构有限元模型建立方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109858054A CN109858054A (zh) | 2019-06-07 |
CN109858054B true CN109858054B (zh) | 2023-05-26 |
Family
ID=66889453
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711246410.4A Active CN109858054B (zh) | 2017-11-30 | 2017-11-30 | 一种飞行器前缘结构有限元模型建立方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109858054B (zh) |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104038295B (zh) * | 2014-06-06 | 2016-03-02 | 西安电子科技大学 | 一种基于机电耦合的变形阵列天线散射性能分析方法 |
CN105205209B (zh) * | 2015-08-20 | 2018-04-10 | 哈尔滨工业大学 | 飞行器平板侧窗的瞬态热响应分析方法 |
CN106570204B (zh) * | 2016-09-23 | 2019-07-23 | 西安交通大学 | 一种基于cpu+gpu异构并行计算的透平机械叶片静强度特性分析方法 |
-
2017
- 2017-11-30 CN CN201711246410.4A patent/CN109858054B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN109858054A (zh) | 2019-06-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11003807B2 (en) | Techniques for generating materials to satisfy design criteria | |
Li et al. | Aerodynamic design optimization of nacelle/pylon position on an aircraft | |
CN108090275B (zh) | 一种面向参数化模型的涡轮气冷叶片构型方法 | |
US20150193559A1 (en) | Method for creating three dimensional lattice structures in computer-aided design models for additive manufacturing | |
CN108009336A (zh) | 一种微桁架结构承载和热防护结构的多目标优化方法 | |
Azamatov et al. | Comprehensive aircraft configuration design tool for integrated product and process development | |
CN106599405B (zh) | 一种计算主减速器与机身任意连接点位置实际载荷的方法 | |
Shimada | Current issues and trends in meshing and geometric processing for computational engineering analyses | |
CN112257321A (zh) | 基于无监督学习的卫星组件温度场确定方法 | |
Mieloszyk et al. | Rapid geometry definition for multidisciplinary design and analysis of an aircraft | |
CN106529005A (zh) | 一种基于有限元仿真的超高强度钢螺纹类工件局部回火感应线圈的设计方法 | |
CN107527385B (zh) | 一种网格自动投影方法 | |
Andrejašič et al. | A mesh morphing based FSI method used in aeronautical optimization applications | |
Sclafani et al. | Drag prediction for the common research model using CFL3D and OVERFLOW | |
CN109858054B (zh) | 一种飞行器前缘结构有限元模型建立方法 | |
CN108319775B (zh) | 一种热传导问题中复合材料的近场动力学建模方法 | |
US20110191068A1 (en) | Multiscale substructures in finite element analysis | |
Ferrero et al. | A numerical method for conjugate heat transfer problems in hypersonic flows | |
Sclafani et al. | Analysis of the Common Research Model Using Structured and Unstructured Meshes | |
CN109766641A (zh) | 基于知识重用的多物理场cae建模方法 | |
CN108804791A (zh) | 一种适用于埋入式进气道布局的飞行器参数化方法 | |
US20220405441A1 (en) | Designing a structural product including structural analysis and post-processing | |
Matsushima et al. | Wing design for supersonic transports using integral equation method | |
Zhu et al. | An efficient method for helicopter fuselage shape optimization | |
CN109241572B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |