CN105512352A - 一种面向气动载荷的分析方法 - Google Patents
一种面向气动载荷的分析方法 Download PDFInfo
- Publication number
- CN105512352A CN105512352A CN201410535096.1A CN201410535096A CN105512352A CN 105512352 A CN105512352 A CN 105512352A CN 201410535096 A CN201410535096 A CN 201410535096A CN 105512352 A CN105512352 A CN 105512352A
- Authority
- CN
- China
- Prior art keywords
- point
- grid
- bin
- interpolation
- pneumatic
- 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
Links
Landscapes
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明提供了一种面向气动载荷的分析方法,其特征在于,包括以下步骤:1)根据计算模型划分气动计算网格,进行流场数值计算,获得气动计算网格点上的压力系数Cp值;2)根据载荷的给定要求对各部件进行网格划分,形成载荷插值网格;3)选择气动计算所使用软件,自动读入气动计算网格的网格点和压力分布数据;4)根据气动计算网格的网格点上压力分布数据对载荷插值网格的网格点上压力系数Cp值进行数据插值。
Description
技术领域:
本发明设计航空气动力设计实验领域。
背景技术:
在型号任务研制中,需要大批量向载荷专业发送气动载荷数据。目前对载荷的插值分析缺乏高精度的方法。传统方法采用四叉树、八叉树等算法,通过寻找离载荷插值网格目标点距离最近的若干气动计算网格上的边界点来进行加权平均。此种方法的缺点是插值精度低,当气动计算网格长细比过大并且流动变化剧烈时,不能精确捕捉到连续的压力系数Cp分布,给计算带来较大的误差。根据气动计算网格计算得到的等压力系数Cp线如图1所示,采用传统方法进行插值计算时,会出现计算结果如图2所示,由图示可见,误差比较大。
发明内容
本发明的任务和目的就是,提供一种面向气动载荷的分析方法。
气动计算网格和载荷插值网格的定义
气动计算网格:根据气动计算软件的要求,用于计算飞机模型的压力系数Cp分布而划分的网格;
载荷插值网格:由于载荷专业通常需要给定位置分布的压力系数Cp,根据载荷给定的位置要求而划分的网格;
面元:对于四边形计算网格来说,由四个顶点组成的一个网格单元定义为面元。
本发明提供了一种面向气动载荷的分析方法,其特征在于,包括以下步骤:
1)根据计算模型划分气动计算网格,进行流场数值计算,获得气动计算网格点上的压力系数Cp值;
2)根据载荷的给定要求对各部件进行网格划分,形成载荷插值网格;
3)选择气动计算所使用软件,自动读入气动计算网格的网格点和压力分布数据;
4)根据气动计算网格的网格点上压力分布数据对载荷插值网格的网格点上压力系数Cp值进行数据插值:
a)采用基于局部加密的方法:对载荷插值网格的网格点上每一个目标插值点,搜索其在气动计算网格中距离最小点所在面元的边界点,并记录最小距离R1,利用所在气动计算网格的面元边界点的Cp值进行局部加密,其中加密后产生的面元网格的边界点用相连面元边界点的距离加权平均计算得到;
b)计算加密后产生的面元网格的边界点与目标插值点的最小距离R2,当R2小于设定的容差R时,将距离R2所在的面元网格点的Cp值赋给目标插值点;
c)如果R2等于设定的容差R,则进行判断:如果R2小于上一轮距离最小值R1,则说明目标插值点在当前面元内,将当前面元进行进一步加密,直到距离最小值R2小于设定的容差R为止;如果R2等于上一轮最小距离值R1,则说明目标插值点不在当前面元内,那么退回到次小距离点所在面元,进行局部加密,重复第b)步的过程,直至当R2小于设定的容差R为止。
本发明提出的载荷数据分析技术,采用基于局部网格加密的方法,对每一个目标插值点,搜索其在源网格中所在面元,利用面元所在的顶点的Cp值,对面元进行局部加密,计算加密后面元网格点与目标插值点的距离,进行二次搜索插值。该方法突破了传统方法对气动网格长细比和流动变化程度敏感的缺点,提高了载荷计算的精度。同时,减少了人工干预,提高了工作效率,缩短飞机设计周期,对飞机设计流程具有一定实用意义。
附图说明
图1为气动计算网格等压力系数Cp线计算结果示意图;
图2为传统方法载荷插值网格等压力系数Cp线插值结果;
图3为本发明载荷插值网格等压力系数Cp线插值结果。
具体实施方式
一种面向气动载荷的分析方法,其特征在于,包括以下步骤:
1)根据计算模型划分气动计算网格,进行流场数值计算,获得气动计算网格点上的压力系数Cp值;
2)根据载荷的给定要求对各部件进行网格划分,形成载荷插值网格;
3)选择气动计算所使用软件,自动读入气动计算网格的网格点和压力分布数据;
4)根据气动计算网格的网格点上压力分布数据对载荷插值网格的网格点上压力系数Cp值进行数据插值:
a)采用基于局部加密的方法:对载荷插值网格的网格点上每一个目标插值点,搜索其在气动计算网格中距离最小点所在面元的边界点,并记录最小距离R1,利用所在气动计算网格的面元边界点的Cp值进行局部加密,其中加密后产生的面元网格的边界点用相连面元边界点的距离加权平均计算得到;
b)计算加密后产生的面元网格的边界点与目标插值点的最小距离R2,当R2小于设定的容差R时,将距离R2所在的面元网格点的Cp值赋给目标插值点;
c)如果R2等于设定的容差R,则进行判断:如果R2小于上一轮距离最小值R1,则说明目标插值点在当前面元内,将当前面元进行进一步加密,直到距离最小值R2小于设定的容差R为止;如果R2等于上一轮最小距离值R1,则说明目标插值点不在当前面元内,那么退回到次小距离点所在面元,进行局部加密,重复第b)步的过程,直至当R2小于设定的容差R为止。
得到结果如图3所示。可以看出,采用该插值方法能够捕捉到连续的压力分布,使精度满足工程需求。
本发明提出的载荷数据分析技术,采用基于局部网格加密的方法,对每一个目标插值点,搜索其在源网格中所在面元,利用面元所在的顶点的Cp值,对面元进行局部加密,计算加密后面元网格点与目标插值点的距离,进行二次搜索插值。该方法突破了传统方法对气动网格长细比和流动变化程度敏感的缺点,提高了载荷计算的精度。同时,减少了人工干预,提高了工作效率,缩短飞机设计周期,对飞机设计流程具有一定实用意义。
Claims (1)
1.一种面向气动载荷的分析方法,其特征在于,包括以下步骤:
1)根据计算模型划分气动计算网格,进行流场数值计算,获得气动计算网格点上的压力系数Cp值;
2)根据载荷的给定要求对各部件进行网格划分,形成载荷插值网格;
3)选择气动计算所使用软件,自动读入气动计算网格的网格点和压力分布数据;
4)根据气动计算网格的网格点上压力分布数据对载荷插值网格的网格点上压力系数Cp值进行数据插值:
a)采用基于局部加密的方法:对载荷插值网格的网格点上每一个目标插值点,搜索其在气动计算网格中距离最小点所在面元的边界点,并记录最小距离R1,利用所在气动计算网格的面元边界点的Cp值进行局部加密,其中加密后产生的面元网格的边界点用相连面元边界点的距离加权平均计算得到;
b)计算加密后产生的面元网格的边界点与目标插值点的最小距离R2,当R2小于设定的容差R时,将距离R2所在的面元网格点的Cp值赋给目标插值点;
c)如果R2等于设定的容差R,则进行判断:如果R2小于上一轮距离最小值R1,则说明目标插值点在当前面元内,将当前面元进行进一步加密,直到距离最小值R2小于设定的容差R为止;如果R2等于上一轮最小距离值R1,则说明目标插值点不在当前面元内,那么退回到次小距离点所在面元,进行局部加密,重复第b)步的过程,直至当R2小于设定的容差R为止。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410535096.1A CN105512352B (zh) | 2014-10-11 | 2014-10-11 | 一种面向气动载荷的分析方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410535096.1A CN105512352B (zh) | 2014-10-11 | 2014-10-11 | 一种面向气动载荷的分析方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105512352A true CN105512352A (zh) | 2016-04-20 |
CN105512352B CN105512352B (zh) | 2018-05-29 |
Family
ID=55720332
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410535096.1A Active CN105512352B (zh) | 2014-10-11 | 2014-10-11 | 一种面向气动载荷的分析方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105512352B (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106156444A (zh) * | 2016-07-27 | 2016-11-23 | 中国航空工业集团公司西安飞机设计研究所 | 一种飞机气动载荷处理方法及飞机强度计算方法 |
CN107194079A (zh) * | 2017-05-25 | 2017-09-22 | 中国商用飞机有限责任公司 | 像素细分载荷传递方法和系统 |
CN111008481A (zh) * | 2019-12-18 | 2020-04-14 | 北京空天技术研究所 | 航天器的气动分析方法及装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101599104B (zh) * | 2009-07-16 | 2011-06-22 | 北京航空航天大学 | 一种航空涡轮发动机叶片颤振边界的模拟方法 |
CN102012953B (zh) * | 2010-11-04 | 2013-05-08 | 西北工业大学 | Cfd/csd耦合求解非线性气动弹性仿真方法 |
-
2014
- 2014-10-11 CN CN201410535096.1A patent/CN105512352B/zh active Active
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106156444A (zh) * | 2016-07-27 | 2016-11-23 | 中国航空工业集团公司西安飞机设计研究所 | 一种飞机气动载荷处理方法及飞机强度计算方法 |
CN107194079A (zh) * | 2017-05-25 | 2017-09-22 | 中国商用飞机有限责任公司 | 像素细分载荷传递方法和系统 |
CN107194079B (zh) * | 2017-05-25 | 2020-12-15 | 中国商用飞机有限责任公司 | 像素细分载荷传递方法和系统 |
CN111008481A (zh) * | 2019-12-18 | 2020-04-14 | 北京空天技术研究所 | 航天器的气动分析方法及装置 |
CN111008481B (zh) * | 2019-12-18 | 2024-01-12 | 北京空天技术研究所 | 航天器的气动分析方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN105512352B (zh) | 2018-05-29 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105631125A (zh) | 基于降阶模型的气动-热-结构耦合分析方法 | |
KR101612506B1 (ko) | 전산 유체 역학을 이용한 항공기 공력해석 시스템 및 방법 | |
CN106777482A (zh) | 一种基于网格参数化的结构多学科设计优化方法 | |
CN104112167B (zh) | 可发电风资源分布的获取方法 | |
CN105512352A (zh) | 一种面向气动载荷的分析方法 | |
CN102682348A (zh) | 复杂装备部件维修级别优化系统及其建立方法 | |
CN114676522B (zh) | 融合gan和迁移学习的气动形状优化设计方法及系统及设备 | |
CN107871181B (zh) | 一种热电联产机组发电计划制定方法、装置及系统 | |
CN115688276A (zh) | 一种基于离散伴随方法的飞行器外形自动化优化方法、系统、设备、介质 | |
CN106709179B (zh) | 一种叶片气动力分布数据处理方法 | |
CN103699756A (zh) | 松耦合模拟系统及方法 | |
US11699013B2 (en) | Computer-implemented method for simulation of an electrical circuit | |
CN104657554A (zh) | 一种火箭三维模型设计信息快速提取系统 | |
CN105093932B (zh) | 一种确定lpv变增益控制器的鲁棒性的方法 | |
Liu et al. | A generalized Griffith importance measure for components with multiple state transitions | |
CN103218493A (zh) | 一种基于多重网格的快速等几何分析数值模拟方法 | |
CN104573331B (zh) | 一种基于MapReduce的K近邻数据预测方法 | |
GB2549028A (en) | Integrated a priori uncertainty parameter architecture in simulation model creation | |
CN105138766A (zh) | 基于模糊聚类的高超声速气动热降阶模型的加点方法 | |
US20120316842A1 (en) | Automatic extremum detection on a surface mesh of a component | |
CN103942368A (zh) | 一种激光切割机床的结构设计方法 | |
CN107194079B (zh) | 像素细分载荷传递方法和系统 | |
CN116224070A (zh) | 一种基于ocv-soc模型的全温度锂离子电池ocv评估方法及系统 | |
CN104850914B (zh) | 一种基于特征建模的新能源发电量预测方法 | |
CN105550424A (zh) | 一种基于rbf网格变形插值序列的筛选方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |