CN109657297A - 一种气动载荷与结构有限元模型双向传递方法 - Google Patents
一种气动载荷与结构有限元模型双向传递方法 Download PDFInfo
- Publication number
- CN109657297A CN109657297A CN201811449284.7A CN201811449284A CN109657297A CN 109657297 A CN109657297 A CN 109657297A CN 201811449284 A CN201811449284 A CN 201811449284A CN 109657297 A CN109657297 A CN 109657297A
- Authority
- CN
- China
- Prior art keywords
- point
- aerodynamic
- displacement
- structural
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 39
- 238000006073 displacement reaction Methods 0.000 claims abstract description 57
- 238000005516 engineering process Methods 0.000 claims abstract description 5
- 239000011159 matrix material Substances 0.000 claims description 9
- 239000008186 active pharmaceutical agent Substances 0.000 claims description 3
- 238000005452 bending Methods 0.000 claims description 3
- 239000007787 solid Substances 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 3
- 230000005540 biological transmission Effects 0.000 description 2
- 230000036244 malformation Effects 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 241000201966 Goura Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000000547 structure data Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
Classifications
-
- 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)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
一种气动载荷与结构有限元模型双向传递方法,属于飞机结构强度技术领域,包括如下步骤:一、构建结构模型,划分气动网格;二、提取气动网格的网格点,构建气动表面点到结构单元的MPC;三、N‑S方程计算气动表面点的节点力;四、通过MPC,采用虚功原理将节点力加载到结构有限元模型,利用Nastran线性或非线性求解序列进行结构计算得到结构点位移;五、通过MPC计算得到气动表面点的位移,采用弹性变形技术生成新的气动网格;六、重复步骤三至五,直至收敛,得到结构模型的气动特性。该方法效率高,适合处理复杂的流/固体数据传递、大型数据插值问题。
Description
技术领域
本申请属于飞机结构强度技术领域,特别涉及一种基于多点约束的气动载荷与结构有限元模型双向传递方法。
背景技术
对于气动与结构间载荷位移传递方法,早在1970年Harder和Desmarais提出了无限平板样条(IPS)插值方法,该方法基于无限平板的载荷位移偏微分平衡方程和虚功原理推导得出,IPS方法计算简单且可同时实现载荷与位移的传递,但其应用范围仅限于二维问题。随后Appa将IPS插值方法改进为薄板样条插值(TPS)方法,TPS方法有效实现了三维梁静气弹计算中载荷位移传递。2001年Goura等人提出了一种常体积转换(CVT)方法,其基本假设为结构变形前后气动节点在结构三角形内的面积坐标保持不变,由结构节点和气动节点组成的四面体的体积在变形过程中也守恒,这样就可根据结构位移变形计算得出气动节点的位移,徐敏等通过面积因子限制提出了对该法的修正形式(ICVT),解决了变形畸形过渡时,CVT方法气动网格质量较差甚至出现负体积等问题。CVT方法只限于两种网格间的位移传递,对载荷传递并不适用,这限制了它的应用面。与此类似的还有2003年徐敏等受有限元四边形单元形函数启发提出的一种有限元四节点(FEFN)插值方法,该方法计算简单,计算量小,但也仅适合于气动与结构网格间的位移插值使用。
因此,希望有一种技术方案来克服或至少减轻现有技术的至少一个上述缺陷。
发明内容
本申请的目的是提供一种气动载荷与结构有限元模型双向传递方法,以解决上述问题,实现气动载荷、结构位移数据双向传递的同时,提高数据高精度传递。
本申请的技术方案是:
一种气动载荷与结构有限元模型双向传递方法,包括如下步骤:
一、构建结构模型,划分气动网格;
二、提取所述气动网格的网格点,构建气动表面点到结构单元的MPC,所述MPC指所述气动表面点与所述结构单元内的结构点的位移与力的关系函数;
三、N-S方程计算所述气动表面点的节点力;
四、利用MPC将所述节点力加载到结构有限元模型,利用Nastran线性或非线性求解序列进行结构计算得到结构点位移;
五、通过MPC计算得到所述气动表面点的位移,采用弹性变形技术生成新的气动网格;
六、重复步骤三至五,直至收敛,得到所述结构模型的气动特性。
根据本申请的至少一个实施方式,所述结构模型为考虑弯曲和扭转的梁模型。
根据本申请的至少一个实施方式,所述步骤三中,利用N-S方程通过有限体积法求解梁气动表面点的节点力:
其中,Ω是控制体;δΩ是Ω的边界;dS是δΩ的面积矢量;W,Fc,Fv分别为守恒量、对流项和黏性项。
根据本申请的至少一个实施方式,所述步骤四包括,将所述气动节点力以节点力的行形式通过MPC加载到有限元模型上。
根据本申请的至少一个实施方式,所述步骤四包括:
获得结构点的变形;
根据所述变形获得所述结构点的平动位移向量T及转动位移向量R。
根据本申请的至少一个实施方式,所述获得结构点的变形包括:
根据结构点的刚度矩阵求解结构点的变形;或者
根据结构点的柔度矩阵求解结构点的变形。
根据本申请的至少一个实施方式,若对所述气动表面点影响的结构点数量为两个,则所述步骤五包括:
确定第一结构点位移引起的所述气动表面点的第一位移;
确定第二结构点位移引起的所述气动表面点的第二位移;
获取所述气动表面点到由两个结构点构成的连线的垂点;
根据所述垂点到所述两个结构点的距离分配所述第一位移及所述第二位移的权重;
根据所述第一位移与其权重以及第二位移与其权重确定所述气动表面点的最终位移。
根据本申请的至少一个实施方式,确定第一结构点位移或所述第二结构点位移引起的所述气动表面点的第一位移包括:
获取气动表面点相对于所述第一结构点或所述第二结构点的初始位置向量r;
根据公式Ta=T+R×r计算气动表面点的位移Ta;
其中,T为所述第一结构点或所述第二结构点的平动位移向量,R为所述第一结构点或所述第二结构点的转动向量。
本申请至少存在以下有益技术效果:
本申请提供的气动载荷与结构有限元模型双向传递方法,相比传统的结构气动界面插值方法,具有:
1)可统一处理单梁、蒙皮、实体等不同形式的插值问题;
2)计算时间少、效率高,适合处理大型数据插值问题;
3)适合处理复杂、极度复杂外形的流/固体数据传递;
4)所加载的载荷可统一考虑压力、摩擦力等,加载精度高,无漏载等问题。
附图说明
图1是本申请提供的气动载荷与结构有限元模型双向传递方法流程图;
图2是本申请中构建的结构模型;
图3是本申请梁的有限元模型;
图4是本申请中梁单元的两点约束的示意图;
图5是本申请中四边形单元的多点约束示意图。
具体实施方式
为使本申请实施的目的、技术方案和优点更加清楚,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行更加详细的描述。在附图中,自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。所描述的实施例是本申请一部分实施例,而不是全部的实施例。下面通过参考附图描述的实施例是示例性的,旨在用于解释本申请,而不能理解为对本申请的限制。基于本申请中的实施例,本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例,都属于本申请保护的范围。下面结合附图对本申请的实施例进行详细说明。
下面结合附图1至图5对本申请做进一步详细说明。
本申请提供的一种气动载荷与结构有限元模型双向传递方法,包括如下步骤:
一、构建结构模型,划分气动网格;
二、提取所述气动网格的网格点,构建气动表面点到结构单元的MPC,所述MPC指所述气动表面点与所述结构单元内的结构点的位移与力的关系函数;
三、N-S方程计算所述气动表面点的节点力;
四、通过MPC,采用虚功原理将所述节点力加载到结构有限元模型,利用Nastran线性或非线性求解序列进行结构计算得到结构点位移;
五、通过MPC计算得到所述气动表面点的位移,采用弹性变形技术生成新的气动网格;
六、重复步骤三至五,直至收敛,得到所述气动表面点的气动特性。
在本实施例中,步骤一、如图2与图3所示,构建一种考虑弯曲和扭转的梁(或机翼)模型作为结构模型,将梁模型上的气动表面点应用Nastran进行结构有限元建模,并对结构模型进行划分成气动网格;所述气动网格划分后,提取气动网格上的网格点,气动网格形成一个表面,网格点为气动网格表面上的点,简称为气动表面点。
在本实施例中,步骤二、构建的结构模型是梁模型,气动网格为梁模型上的蒙皮,即代表一种流体、气动表面,所述梁代表结构单元,构建所述气动表面点到所述梁的结构单元的MPC,构建的MPC表示所述气动表面点与所述结构单元内的结构点的位移与力的关系函数,构建的MPC能够实现流场(气动表面)与结构场(结构单元)之间双向数据传递,以便于结构变形的计算。
步骤三,在结构模型产生变形时,利用N-S方程通过有限体积法求解梁气动表面点的节点力:
其中,Ω是控制体;δΩ是Ω的边界;dS是δΩ的面积矢量;W,Fc,Fv分别为守恒量、对流项和黏性项。
步骤四、采用虚功原理将所述气动节点力以节点力的行形式通过MPC加载到有限元模型上,利用Nastran线性或非线性求解序列进行结构计算得到结构点(结构单元上的一点)位移;
获得结构点的变形;
采用Nastran线性求解线性静力学方程:
KX=Fg+Faero
对于柔度矩阵,求解刚度矩阵的逆矩阵,采用IPS插值方法得到梁表面气动网格单元的柔度矩阵C,求解静力学方程:
X=C(Fg+Faero)
其中,K和X分别为刚度矩阵和结构点的变形,Fg和Faero分别为重力和气动表面力。
根据所述变形获得所述结构点的平动位移向量T及转动位移向量R。
步骤五、如图4所示,若对所述气动表面点影响的结构点数量为两个,则所述步骤五包括:
确定第一结构点位移N1引起的所述气动表面点的第一位移Ta1;
Ta1=T1+R1×r1
确定第二结构点位移N2引起的所述气动表面点的第二位移Ta2;
Ta2=T2+R2×r2
其中,气动表面点相对于第一结构点位移N1和第二结构点位移N2的初始位置向量为r1和r2,气动表面点相对于第一结构点位移N1和第二结构点位移N2的初始位置转动位移向量为R1和R2。
获取所述气动表面点到由两个结构点构成的连线的垂点P;
根据所述垂点P到所述两个结构点的距离分配所述第一位移N1及所述第二位移N2的权重;根据所述第一位移N1与其权重以及第二位移N2与其权重确定所述气动表面点的最终位移Ta;如下式:
将上述得到的气动表面点的位移,采用弹性变形技术生成新的气动网格。
步骤六、重复步骤三至步骤五,直至收敛,得到变形后的气动网格,即利用MPC方法完成了气动载荷与结构有限元模型双向传递,得到所述结构模型的气动特性。
在一些可选的实施方式中,若对所述气动表面点影响的结构点数量为多个点,即为多点约束,在该实施例中,以四边形为例,搜索距离所述气动表面点最近的四边形单元,构建气动表面点到结构单元的映射关系,如图5所示,则所述步骤五为:
投影气动点(气动表面点)p到结构单元内于q点,则q点在结构单元内的二维参数坐标(ε,η)。通过以下公式中四边形等参单元的形函数可以计算出q点的四边形内的插值参数N1至N4:
假设在p点的变形过程中保持相对于结构单元的位置不变,即参数坐标不变,将插值系数作为MPC的权系数,构建气动表面点到每个结构点的MPC。通过以下公式计算出每个结构点引起的气动点变形,采用插值系数加权平均计算出单元变形引起的结构变形,如以下公式所示:
rp=N1r1+N2r2+N3r3+N4r4
其中,r1、r2、r3以及r4分别为所述气动点p到N1至N4的初始位置向量。
以上所述,仅为本申请的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应以所述权利要求的保护范围为准。
Claims (8)
1.一种气动载荷与结构有限元模型双向传递方法,其特征在于,包括如下步骤:
一、构建结构模型,划分气动网格;
二、提取所述气动网格的网格点,构建气动表面点到结构单元的MPC,所述MPC指所述气动表面点与所述结构单元内的结构点的位移与力的关系函数;
三、获取所述气动表面点的节点力;
四、利用MPC将所述节点力加载到结构有限元模型,利用Nastran线性或非线性求解序列进行结构计算得到结构点位移;
五、通过MPC计算得到所述气动表面点的位移,采用弹性变形技术生成新的气动网格;
六、重复步骤三至步骤五,直至收敛,得到所述结构模型的气动特性。
2.根据权利要求1所述的气动载荷与结构有限元模型双向传递方法,其特征在于,所述结构模型为考虑弯曲和扭转的梁模型。
3.根据权利要求1所述的气动载荷与结构有限元模型双向传递方法,其特征在于,所述步骤三中,利用如下N-S方程求解梁气动表面点的节点力:
其中,Ω是控制体;δΩ是Ω的边界;dS是δΩ的面积矢量;W,Fc,Fv分别为守恒量、对流项和黏性项。
4.根据权利要求1所述的气动载荷与结构有限元模型双向传递方法,其特征在于,所述步骤四包括:
将所述气动节点力以节点力的行形式通过MPC加载到有限元模型上。
5.根据权利要求1所述的气动载荷与结构有限元模型双向传递方法,其特征在于,所述步骤四包括:
获得结构点的变形;
根据所述变形获得所述结构点的平动位移向量T及转动位移向量R。
6.根据权利要求5所述的气动载荷与结构有限元模型双向传递方法,其特征在于,所述获得结构点的变形包括:
根据结构点的刚度矩阵求解结构点的变形;或者
根据结构点的柔度矩阵求解结构点的变形。
7.根据权利要求5所述的气动载荷与结构有限元模型双向传递方法,其特征在于,若对所述气动表面点影响的结构点数量为两个,则所述步骤五包括:
确定第一结构点位移引起的所述气动表面点的第一位移;
确定第二结构点位移引起的所述气动表面点的第二位移;
获取所述气动表面点到由两个结构点构成的连线的垂点;
根据所述垂点到所述两个结构点的距离分配所述第一位移及所述第二位移的权重;
根据所述第一位移与其权重以及第二位移与其权重确定所述气动表面点的最终位移。
8.根据权利要求7所述的气动载荷与结构有限元模型双向传递方法,其特征在于,确定第一结构点位移或所述第二结构点位移引起的所述气动表面点的第一位移包括:
获取气动表面点相对于所述第一结构点或所述第二结构点的初始位置向量r;
根据公式Ta=T+R×r计算气动表面点的位移Ta;
其中,T为所述第一结构点或所述第二结构点的平动位移向量,R为所述第一结构点或所述第二结构点的转动向量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811449284.7A CN109657297A (zh) | 2018-11-30 | 2018-11-30 | 一种气动载荷与结构有限元模型双向传递方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201811449284.7A CN109657297A (zh) | 2018-11-30 | 2018-11-30 | 一种气动载荷与结构有限元模型双向传递方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109657297A true CN109657297A (zh) | 2019-04-19 |
Family
ID=66111949
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201811449284.7A Pending CN109657297A (zh) | 2018-11-30 | 2018-11-30 | 一种气动载荷与结构有限元模型双向传递方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109657297A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112052525A (zh) * | 2020-09-25 | 2020-12-08 | 中国直升机设计研究所 | 一种气动载荷分配方法 |
CN112131665A (zh) * | 2020-09-23 | 2020-12-25 | 上海交通大学 | 一种基于网格分类的火箭流固载荷混合转换方法 |
CN113792374A (zh) * | 2021-07-19 | 2021-12-14 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种飞机飞行载荷确定方法 |
CN117540613A (zh) * | 2024-01-10 | 2024-02-09 | 西安羚控电子科技有限公司 | 一种气动载荷二次插值处理方法及装置 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108052772A (zh) * | 2017-12-30 | 2018-05-18 | 北京航空航天大学 | 一种基于结构降阶模型的几何非线性静气动弹性分析方法 |
-
2018
- 2018-11-30 CN CN201811449284.7A patent/CN109657297A/zh active Pending
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108052772A (zh) * | 2017-12-30 | 2018-05-18 | 北京航空航天大学 | 一种基于结构降阶模型的几何非线性静气动弹性分析方法 |
Non-Patent Citations (1)
Title |
---|
黄炜等: "基于多点约束的大展弦比机翼静气动弹性计算", 《北京航空航天大学学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112131665A (zh) * | 2020-09-23 | 2020-12-25 | 上海交通大学 | 一种基于网格分类的火箭流固载荷混合转换方法 |
CN112052525A (zh) * | 2020-09-25 | 2020-12-08 | 中国直升机设计研究所 | 一种气动载荷分配方法 |
CN112052525B (zh) * | 2020-09-25 | 2023-03-17 | 中国直升机设计研究所 | 一种气动载荷分配方法 |
CN113792374A (zh) * | 2021-07-19 | 2021-12-14 | 中国航空工业集团公司沈阳飞机设计研究所 | 一种飞机飞行载荷确定方法 |
CN117540613A (zh) * | 2024-01-10 | 2024-02-09 | 西安羚控电子科技有限公司 | 一种气动载荷二次插值处理方法及装置 |
CN117540613B (zh) * | 2024-01-10 | 2024-04-05 | 西安羚控电子科技有限公司 | 一种气动载荷二次插值处理方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109657297A (zh) | 一种气动载荷与结构有限元模型双向传递方法 | |
CN113111430B (zh) | 基于非线性气动力降阶的弹性飞机飞行动力学建模方法 | |
Smith et al. | CFD-based analysis of nonlinear aeroelastic behavior of high-aspect ratio wings | |
CN108427322B (zh) | 一种大柔性飞行器基于在线辨识的建模方法 | |
CN109726437B (zh) | 一种舱门气动载荷等效节点力处理方法 | |
Dang et al. | Accelerated loosely-coupled CFD/CSD method for nonlinear static aeroelasticity analysis | |
CN108363843A (zh) | 一种基于结构降阶模型的几何非线性静气动弹性全机配平方法 | |
CN112580241B (zh) | 一种基于结构降阶模型的非线性气动弹性动响应分析方法 | |
CN109740182A (zh) | 一种基于再生核粒子的无网格物理变形仿真方法 | |
CN109446471B (zh) | 一种考虑载荷不确定性的流固耦合界面数据传递方法 | |
CN111460622A (zh) | 一种基于参数化彩色水平集的多组件布局拓扑优化方法 | |
CN109726465A (zh) | 基于非结构曲边网格的三维无粘低速绕流的数值模拟方法 | |
CN112001109A (zh) | 再生核粒子算法实现结构冲击动力学仿真方法 | |
CN110717216A (zh) | 不规则波下带柔性气囊直升机横摇响应预报方法 | |
Xie et al. | Geometrically nonlinear aeroelastic stability analysis and wind tunnel test validation of a very flexible wing | |
Ritter et al. | Comparison of nonlinear aeroelastic methods for maneuver simulation of very flexible aircraft | |
CN109459206B (zh) | 地面试验非定常气动力加载方法 | |
CN111341449A (zh) | 一种虚拟血管介入手术训练的模拟方法 | |
CN113128085B (zh) | 一种基于状态空间形式涡格法的阵风减缓分析方法 | |
McDaniel et al. | Efficient mesh deformation for computational stability and control analyses on unstructured viscous meshes | |
CN116050303B (zh) | 一种cfd并行计算下的周期性边界条件施加方法 | |
BATOZ et al. | Development of two simple shell elements. | |
CN110147571B (zh) | 一种组件结构的拓扑优化方法及装置 | |
CN109635370B (zh) | 开裂式阻力方向舵静气动弹性特性分析方法 | |
Guo et al. | Nonlinear gust response analysis of free flexible aircraft |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190419 |