CN109459206A - 地面试验非定常气动力加载方法 - Google Patents
地面试验非定常气动力加载方法 Download PDFInfo
- Publication number
- CN109459206A CN109459206A CN201811541338.2A CN201811541338A CN109459206A CN 109459206 A CN109459206 A CN 109459206A CN 201811541338 A CN201811541338 A CN 201811541338A CN 109459206 A CN109459206 A CN 109459206A
- Authority
- CN
- China
- Prior art keywords
- structural
- triangle
- aerodynamic force
- pneumatic
- unsteady aerodynamic
- 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
- 238000011068 loading method Methods 0.000 title claims abstract description 17
- 238000002474 experimental method Methods 0.000 title abstract 2
- 238000000034 method Methods 0.000 claims abstract description 57
- 238000012360 testing method Methods 0.000 claims abstract description 36
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 13
- 238000005259 measurement Methods 0.000 claims abstract description 11
- 239000012530 fluid Substances 0.000 claims abstract description 10
- 239000011159 matrix material Substances 0.000 claims description 27
- 238000006073 displacement reaction Methods 0.000 claims description 22
- 230000009467 reduction Effects 0.000 claims description 19
- 230000004044 response Effects 0.000 claims description 17
- 230000005284 excitation Effects 0.000 claims description 11
- 230000008569 process Effects 0.000 claims description 10
- 238000006243 chemical reaction Methods 0.000 claims description 9
- 238000012545 processing Methods 0.000 claims description 6
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 3
- 238000000354 decomposition reaction Methods 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 description 8
- 238000012512 characterization method Methods 0.000 description 3
- 238000004088 simulation Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- 229910000838 Al alloy Inorganic materials 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000001747 exhibiting effect Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M9/00—Aerodynamic testing; Arrangements in or on wind tunnels
- G01M9/06—Measuring arrangements specially adapted for aerodynamic testing
Landscapes
- Physics & Mathematics (AREA)
- Fluid Mechanics (AREA)
- General Physics & Mathematics (AREA)
- Aerodynamic Tests, Hydrodynamic Tests, Wind Tunnels, And Water Tanks (AREA)
Abstract
本发明公开了一种地面试验非定常气动力加载方法，用于解决现有非定常气动力加载方法实用性差的技术问题。技术方案是首先对试验结构测量点数据进行处理，其次采用常体积转换法实现从分布式非定常气动力到激振集中力的减缩，再建立基于计算流体力学的非定常气动力降阶模型，将基于Volterra级数的非定常气动力降阶模型转换成状态空间形式。由于采用常体积转换法实现从分布式非定常气动力到激振集中力的减缩，同时适用于二维和三维复杂表面；基于Volterra级数理论建立基于计算流体力学的非定常气动力降阶模型，计算精度高，刻画非线性系统能力强，使用系统最小特征实现算法将降阶模型转换成状态空间形式便于工程应用，实用性好。
Description
技术领域
本发明涉及一种非定常气动力加载方法，特别涉及一种地面试验非定常气动力加载方法。
背景技术
为了得到飞行器的结构性能，通常需要采用风洞试验的方法，但是风洞试验成本较高，而且对于一些特殊的试验，试验模型设计工作、飞行工况的真实模拟会变得非常复杂。因此，国内外研究人员尝试采用地面试验来代替常规风洞试验，即称之为干风洞试验。文献“GVT-BASED GROUND FLUTTER TEST WITHOUT WIND TUNNEL，52nd AIAA/ASME/ASCE/AHS/ASC Structures,Structural Dynamics and Materials Conference，AIAA Paper2011-1942,2011”公开了一种应用于地面干风洞试验的非定常气动力加载方法，该法采用无限平板样条插值法实现非定常气动力减缩，利用工程算法如活塞理论或者偶极子格网法建立非定常气动力降阶模型，对矩形铝合金平板进行了气动弹性分析，利用该法进行地面试验得到的平板颤振边界与参考值吻合的较好。文献公开的方法适用于干风洞试验，相比于风洞试验可以显著降低试验成本并提高试验效率。文献所述方法存在两个主要问题，一是非定常气动力减缩技术采用无限平板样条法，该法具有以下的局限性：1)在二维情况下要求所有的结构点不能共线，不能有重合的点，三维情况下要求所有的结构点不能共面，也不能有重合的点。2)该法是一种标量方法，对于给定了一个方向的位移或者力，就不能求解在其他方向上的位移。另一个问题在于采用工程算法建立非定常气动力降阶模型，其计算精度较低，对复杂结构或者非线性较强的情况表征能力不足。
发明内容
为了克服现有非定常气动力加载方法实用性差的不足，本发明提供一种地面试验非定常气动力加载方法。该方法首先对试验结构测量点数据进行处理，其次采用常体积转换法实现从分布式非定常气动力到激振集中力的减缩，再建立基于计算流体力学的非定常气动力降阶模型，使用系统最小特征实现算法将基于Volterra级数的非定常气动力降阶模型转换成状态空间形式。由于采用常体积转换法实现从分布式非定常气动力到激振集中力的减缩，同时适用于二维和三维复杂表面；基于Volterra级数理论建立基于计算流体力学的非定常气动力降阶模型，计算精度高，刻画非线性系统能力强，使用系统最小特征实现算法将降阶模型转换成状态空间形式便于工程应用，实用性好。
本发明解决其技术问题所采用的技术方案：一种地面试验非定常气动力加载方法，其特点是包括以下步骤：
步骤一、对试验结构测量点数据进行处理，在小变形情况下，结构变形表达为固有振型的叠加形式，即
其中，x为结构实际位移变形，φi为第i阶振型信息，ξi为第i阶广义位移。
记xc为测量点采集得到的位移变形信息，Φc为测量点对应的振型信息，ξc为广义位移，将式(1)改写为矩阵表达形式：
xc＝Φcξc (2)
根据式(2)得到由测量点变形表达的广义位移
采集到测量点变形信息之后通过式(3)计算得到广义位移。
步骤二、采用常体积转换法实现从分布式非定常气动力到激振集中力的减缩，在结构网格中为每一个气动网格节点寻找最小三角形。具体的搜索过程为，对每一个气动网格节点ga，寻找包含该气动网格节点的四边形结构网格单元，通过面积坐标法实现。气动网格投影点与四个结构节点组成四个三角形，面积分别为A1、A2、A3和A4。四边形网格的面积为A，该四边形结构网格单元包含该气动网格投影点时面积之比为 (A1+A2+A3+A4)/A＝1.0+ε，ε是给定误差。
包含气动网格投影点的结构单元四个顶点构成四个相邻的结构三角形，通过面积坐标的方法判断包含气动节点的三角形单元，气动网格点到结构三角形顶点的距离Li表示为
其中，(xa,ya,za)是气动网格点ga的坐标，是第m个结构三角形第i个顶点坐标。第m个三角形的最大距离为
根据式(5)，具有最小的三角形就是关于气动网格点ga的最小结构三角形。将常体积变换插值过程用标量方程组表述，在每一个ga所对应的最小结构三角形单元中，(x1,y1,z1)，(x2,y2,z2)和(x3,y3,z3)表示结构三个顶点坐标，(x,y,z)表示三角单元内某一任意点gs的坐标。结构三角单元组成的平面方程为
式中，α+β+γ＝1。ga与结构三角单元组成的四面体体积为V＝1/3SH，S是结构三角单元的面积，H是气动点到结构三角单元平面的距离，有
式中，l、m和n分别是叉乘积向量的三个分量。将上式代入体积表达式中，有
记ga在结构三角单元平面上的投影为坐标是(xp,yp,zp)。则气动点ga与其投影点组成的直线方程是
其中，μ是待求系数。联立方程(6)和(10)以及α+β+γ＝1，组成一个新的方程组
由上式求出α、β、γ和μ，再回代到式(10)中得到的值，之后用同样的方法寻找最小结构三角形。
记第j个气动网格点上的载荷为
其中，和是气动网格点在x、y和z方向的力。则每个最小结构三角形节点上的载荷根据下式计算
其中，F1 j、F2 j和F3 j是由Fj转换到结构三角形三个顶点上的载荷。α、β和γ是气动网格点j在结构三角形中的面积坐标。对于每个结构节点，作用在其上的总载荷等于所有Fj,j＝1,2,…,Na转换到其上的载荷之和，通过上述过程组合得到气动力减缩矩阵。
步骤三、建立基于计算流体力学的非定常气动力降阶模型。在小扰动下Euler方程和N-S方程具有弱非线性特性，因此非定常气动力能够以二阶Volterra级数的形式精确表示：
选用阶跃响应激励气弹系统，通过辨识近似一阶核函数建立基于Volterra级数的降阶模型。近似一阶Volterra核函数的定义如下：
式中，s(n)为阶跃响应，n为离散时间步，ξ0为阶跃响应的幅值。
步骤四、使用系统最小特征实现算法将基于Volterra级数的非定常气动力降阶模型转换成状态空间形式，系统最小特征实现算法得到的线性时不变离散状态空间形式如下：
Aa,Ba,Ca,Da分别对应于气动力系统的系统矩阵、输入矩阵、输出矩阵和前馈矩阵，xa为系统的状态变量，ξ为系统的输入，Fa为系统的非定常气动力输出，q为动压。
系统的脉冲响应输出为：
Fa(0)＝Da
记为M×N的矩阵，将Fa(n)用近似一阶核函数代替，构造Hankel矩阵：
其中，α,β为可调整的整数。对Hαβ(0)作奇异值分解
Hαβ(0)＝UΣVT (19)
得到矩阵Aa,Ba,Ca,Da的表达式如下
其中和的定义如下
本发明的有益效果是：该方法首先对试验结构测量点数据进行处理，其次采用常体积转换法实现从分布式非定常气动力到激振集中力的减缩，再建立基于计算流体力学的非定常气动力降阶模型，使用系统最小特征实现算法将基于Volterra级数的非定常气动力降阶模型转换成状态空间形式。由于采用常体积转换法实现从分布式非定常气动力到激振集中力的减缩，同时适用于二维和三维复杂表面；基于Volterra级数理论建立基于计算流体力学的非定常气动力降阶模型，计算精度高，刻画非线性系统能力强，使用系统最小特征实现算法将降阶模型转换成状态空间形式便于工程应用，实用性好。
下面结合附图和具体实施方式对本发明作详细说明。
附图说明
图1是本发明地面试验非定常气动力加载方法的流程图。
图2是本发明方法实施例面积坐标定义图。
图3是本发明方法实施例气动网格点相邻的结构三角形示意图。
图4是本发明方法实施例两种数值方法以及试验所得无量纲颤振速度结果对比图。
图5是本发明方法实施例两种数值方法以及试验所得无量纲颤振频率结果对比图。
图6是本发明方法实施例马赫数0.96下一阶气动力系数响应结果对比图。
具体实施方式
参照图1-6。
对国际标准气动弹性算例AGARD 445.6进行模拟仿真。AGARD 445.6机翼是 NASA兰利研究中心跨声速风洞中用于研究颤振特性的弹性模型，现已成为研究非线性非定常气动力的标准气动弹性算例。该机翼是45°后掠翼，翼型为NACA64A004，梢根比为0.6576，展现比为1.65。
本发明地面试验非定常气动力加载方法具体步骤如下：
1、对试验结构测量点数据进行处理，在建立非定常气动力降阶模型的过程中，为了方便统一处理，需要将试验结构测量点采集得到的变形信息转化为广义位移的形式，根据模态叠加原理，在小变形情况下，结构的变形表达为固有振型的叠加形式，即
其中x为结构实际位移变形，φi为第i阶振型信息，ξi为第i阶广义位移。
记xc为测量点采集得到的位移变形信息，Φc为测量点对应的振型信息，ξc为广义位移，将式(1)改写为矩阵表达形式：
xc＝Φcξc (2)
根据式(2)得到由测量点变形表达的广义位移
采集到测量点变形信息之后通过式(3)计算得到广义位移。
2、气动力减缩技术其功能是实现分布式气动力到激振点处气动力的减缩，将分布式气动力减缩为激振点处集中力的形式，采用常体积转换法实现，常体积转换方法属于表面跟踪法，它利用变形前后的体积守恒来实现插值计算。包括投影、展开和恢复表面三个过程。常体积转换法可以实现位移和载荷信息的传递，具有适用于三维复杂表面气动力减缩的特点。
首先在结构网格中为每一个气动网格节点寻找最小三角形。具体的搜索过程为，对每一个气动网格节点ga，寻找包含该气动网格节点的四边形结构网格单元，通过面积坐标法来实现，图2中白色点为四边形网格单元的顶点，黑色点表示气动网格点。气动网格投影点与四个结构节点组成四个三角形，面积分别为A1、A2、A3和A4。四边形网格的面积为A，该四边形结构网格单元包含该气动网格投影点时面积之比为 (A1+A2+A3+A4)/A＝1.0+ε，ε是给定误差。
包含气动网格投影点的结构单元四个顶点构成四个相邻的结构三角形(参见图3)，仍然通过面积坐标的方法判断包含气动节点的三角形单元，气动网格点到结构三角形顶点的距离Li表示为
其中(xa,ya,za)是气动网格点ga的坐标，是第m个结构三角形第i个顶点坐标。第m个三角形的最大距离为
根据式(5)，具有最小的三角形就是关于气动网格点ga的最小结构三角形。将常体积变换插值过程用标量方程组来表述，在每一个ga所对应的最小结构三角形单元中，(x1,y1,z1)，(x2,y2,z2)和(x3,y3,z3)表示结构三个顶点坐标，(x,y,z)表示三角单元内某一任意点gs的坐标。结构三角单元组成的平面方程为
式中α+β+γ＝1。ga与结构三角单元组成的四面体体积为V＝1/3SH，S是结构三角单元的面积，H是气动点到结构三角单元平面的距离，有
式中的l、m和n分别是叉乘积向量的三个分量。将上式代入体积表达式中，有
记ga在结构三角单元平面上的投影为坐标是(xp,yp,zp)。则气动点ga与其投影点组成的直线方程是
其中μ是待求系数。联立方程(6)和(10)以及α+β+γ＝1，组成一个新的方程组
由上式求出α、β、γ和μ，再回代到式(10)中得到的值，之后用同样的方法寻找最小结构三角形。
记第j个气动网格点上的载荷为
其中和是气动网格点在x、y和z方向的力。则每个最小结构三角形节点上的载荷根据下式计算
其中，F1 j、F2 j和F3 j是由Fj转换到结构三角形三个顶点上的载荷。α、β和γ是气动网格点j在结构三角形中的面积坐标。对于每个结构节点，作用在其上的总载荷等于所有Fj,j＝1,2,…,Na转换到其上的载荷之和，通过上述过程组合得到气动力减缩矩阵。
3、基于Volterra级数理论建立基于计算流体力学的非定常气动力降阶模型。在小扰动下Euler方程和N-S方程具有弱非线性特性，因此非定常气动力能够以二阶Volterra级数的形式来精确表示：
通过辨识气动弹性系统近似一阶Volterra核函数建立基于Volterra级数的降阶模型。由于阶跃响应的数值稳定性更好，且相比于脉冲响应能够更好的刻画系统特性，因此选用阶跃响应来激励气弹系统，本发明所述方法采用的阶跃响应幅值为ξ0＝10-5，选取结构前4阶模态进行激励，定义近似一阶Volterra核函数的表达形式如下：
式中，s(n)为阶跃响应。利用计算流体力学求解器进行计算并辨识出近似一阶核函数，进而建立基于Volterra级数的非定常气动力降阶模型；
4、使用系统最小特征实现算法将基于Volterra级数的非定常气动力降阶模型转换成状态空间形式，系统最小特征实现算法得到的线性时不变离散状态空间形式如下：
Aa,Ba,Ca,Da分别对应于气动力系统的系统矩阵、输入矩阵、输出矩阵和前馈矩阵，xa为系统的状态变量，ξ为系统的输入，Fa为系统的非定常气动力输出，q为动压。
系统的脉冲响应输出为：
记为M×N的矩阵，参数M和N的值均为4，将Fa(n)用近似一阶核函数代替，构造Hankel矩阵如下
其中α的值为800，β的值为60。令k＝1，对Hαβ(0)作奇异值分解
Hαβ(0)＝UΣVT (19)
定义矩阵Aa,Ba,Ca,Da的表达式如下
其中和的定义如下
至此即依据本发明所述方法建立得到了地面试验非定常气动力快速加载模型，接下来给出该模型预测得到的气动弹性响应结果以及与试验结果和CFD/CSD全阶模型仿真结果之间的对比。从图4、图5可以看出，在Ma<1时，CFD/CSD全阶模型计算结果，地面试验非定常气动力快速加载模型计算结果均与试验结果吻合的较好，在 Ma>1时CFD/CSD全阶模型结果高于试验值，这种现象同样出现在许多已有文献中，表明CFD/CSD全阶模型计算结果是正确的，此时可看到地面试验非定常气动力快速加载模型计算结果与CFD/CSD全阶模型计算结果仍然吻合的很好。
跨声速状态下系统的非线性程度较强，图6为分别使用CFD/CSD全阶模型和地面试验非定常气动力快速加载模型计算得到的马赫数为0.96，颤振边界附近的一阶气动力系数响应结果对比图，可以看出虽然两者之间有些许差异，但仍然可以满足工程高精度要求。
Claims (1)
1.一种地面试验非定常气动力加载方法，其特征在于包括以下步骤：
步骤一、对试验结构测量点数据进行处理，在小变形情况下，结构变形表达为固有振型的叠加形式，即
其中，x为结构实际位移变形，φi为第i阶振型信息，ξi为第i阶广义位移；
记xc为测量点采集得到的位移变形信息，Φc为测量点对应的振型信息，ξc为广义位移，将式(1)改写为矩阵表达形式：
xc＝Φcξc (2)
根据式(2)得到由测量点变形表达的广义位移
采集到测量点变形信息之后通过式(3)计算得到广义位移；
步骤二、采用常体积转换法实现从分布式非定常气动力到激振集中力的减缩，在结构网格中为每一个气动网格节点寻找最小三角形；具体的搜索过程为，对每一个气动网格节点ga，寻找包含该气动网格节点的四边形结构网格单元，通过面积坐标法实现；气动网格投影点与四个结构节点组成四个三角形，面积分别为A1、A2、A3和A4；四边形网格的面积为A，该四边形结构网格单元包含该气动网格投影点时面积之比为(A1+A2+A3+A4)/A＝1.0+ε，ε是给定误差；
包含气动网格投影点的结构单元四个顶点构成四个相邻的结构三角形，通过面积坐标的方法判断包含气动节点的三角形单元，气动网格点到结构三角形顶点的距离Li表示为
其中，(xa,ya,za)是气动网格点ga的坐标，是第m个结构三角形第i个顶点坐标；第m个三角形的最大距离为
根据式(5)，具有最小的三角形就是关于气动网格点ga的最小结构三角形；将常体积变换插值过程用标量方程组表述，在每一个ga所对应的最小结构三角形单元中，(x1,y1,z1)，(x2,y2,z2)和(x3,y3,z3)表示结构三个顶点坐标，(x,y,z)表示三角单元内某一任意点gs的坐标；结构三角单元组成的平面方程为
式中，α+β+γ＝1；ga与结构三角单元组成的四面体体积为V＝1/3SH，S是结构三角单元的面积，H是气动点到结构三角单元平面的距离，有
式中，l、m和n分别是叉乘积向量的三个分量；将上式代入体积表达式中，有
记ga在结构三角单元平面上的投影为坐标是(xp,yp,zp)；则气动点ga与其投影点组成的直线方程是
其中，μ是待求系数；联立方程(6)和(10)以及α+β+γ＝1，组成一个新的方程组
由上式求出α、β、γ和μ，再回代到式(10)中得到的值，之后用同样的方法寻找最小结构三角形；
记第j个气动网格点上的载荷为
其中，和是气动网格点在x、y和z方向的力；则每个最小结构三角形节点上的载荷根据下式计算
其中，F1 j、和是由Fj转换到结构三角形三个顶点上的载荷；α、β和γ是气动网格点j在结构三角形中的面积坐标；对于每个结构节点，作用在其上的总载荷等于所有Fj,j＝1,2,…,Na转换到其上的载荷之和，通过上述过程组合得到气动力减缩矩阵；
步骤三、建立基于计算流体力学的非定常气动力降阶模型；在小扰动下Euler方程和N-S方程具有弱非线性特性，因此非定常气动力能够以二阶Volterra级数的形式精确表示：
选用阶跃响应激励气弹系统，通过辨识近似一阶核函数建立基于Volterra级数的降阶模型；近似一阶Volterra核函数的定义如下：
式中，s(n)为阶跃响应，n为离散时间步，ξ0为阶跃响应的幅值；
步骤四、使用系统最小特征实现算法将基于Volterra级数的非定常气动力降阶模型转换成状态空间形式，系统最小特征实现算法得到的线性时不变离散状态空间形式如下：
Aa,Ba,Ca,Da分别对应于气动力系统的系统矩阵、输入矩阵、输出矩阵和前馈矩阵，xa为系统的状态变量，ξ为系统的输入，Fa为系统的非定常气动力输出，q为动压；系统的脉冲响应输出为：
Fa(0)＝Da
记为M×N的矩阵，将Fa(n)用近似一阶核函数代替，构造Hankel矩阵：
其中，α,β为可调整的整数；对Hαβ(0)作奇异值分解
Hαβ(0)＝UΣVT (19)
得到矩阵Aa,Ba,Ca,Da的表达式如下
其中和的定义如下
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811541338.2A CN109459206B (zh) | 2018-12-17 | 2018-12-17 | 地面试验非定常气动力加载方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811541338.2A CN109459206B (zh) | 2018-12-17 | 2018-12-17 | 地面试验非定常气动力加载方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109459206A true CN109459206A (zh) | 2019-03-12 |
| CN109459206B CN109459206B (zh) | 2020-10-27 |
Family
ID=65613486
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811541338.2A Expired - Fee Related CN109459206B (zh) | 2018-12-17 | 2018-12-17 | 地面试验非定常气动力加载方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109459206B (zh) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110162822A (zh) * | 2019-03-19 | 2019-08-23 | 北京机电工程研究所 | 耦合结构模态的时域快速非定常气动力计算方法 |
| CN110705154A (zh) * | 2019-09-24 | 2020-01-17 | 中国航空工业集团公司西安飞机设计研究所 | 航空器开环气动伺服弹性系统模型均衡降阶的优选方法 |
| CN114818550A (zh) * | 2022-06-30 | 2022-07-29 | 中国飞机强度研究所 | 一种飞机振动试验中时变气动载荷地面等效模拟方法 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050225747A1 (en) * | 2002-06-16 | 2005-10-13 | Rafael - Armament Development Authothority Ltd. | Method and system for evaluating optical disturbances occurring in a supersonic flow field |
| CN101532906A (zh) * | 2009-04-27 | 2009-09-16 | 东南大学 | 风力发电机叶片的流体动力学和结构力学分析方法 |
| CN103364170A (zh) * | 2013-07-05 | 2013-10-23 | 北京航空航天大学 | 气动弹性稳定性的地面模拟预测方法和系统 |
| CN104133933A (zh) * | 2014-05-29 | 2014-11-05 | 温州职业技术学院 | 一种高超声速飞行器热环境下气动弹性力学特性分析方法 |
| CN104317980A (zh) * | 2014-08-26 | 2015-01-28 | 中国直升机设计研究所 | 一种共轴刚性旋翼气动弹性响应方法 |
| CN105115692A (zh) * | 2015-08-12 | 2015-12-02 | 中国航天空气动力技术研究院 | 一种cfd数值模拟与风洞试验相结合的气动预测方法 |
-
2018
- 2018-12-17 CN CN201811541338.2A patent/CN109459206B/zh not_active Expired - Fee Related
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20050225747A1 (en) * | 2002-06-16 | 2005-10-13 | Rafael - Armament Development Authothority Ltd. | Method and system for evaluating optical disturbances occurring in a supersonic flow field |
| CN101532906A (zh) * | 2009-04-27 | 2009-09-16 | 东南大学 | 风力发电机叶片的流体动力学和结构力学分析方法 |
| CN103364170A (zh) * | 2013-07-05 | 2013-10-23 | 北京航空航天大学 | 气动弹性稳定性的地面模拟预测方法和系统 |
| CN104133933A (zh) * | 2014-05-29 | 2014-11-05 | 温州职业技术学院 | 一种高超声速飞行器热环境下气动弹性力学特性分析方法 |
| CN104317980A (zh) * | 2014-08-26 | 2015-01-28 | 中国直升机设计研究所 | 一种共轴刚性旋翼气动弹性响应方法 |
| CN105115692A (zh) * | 2015-08-12 | 2015-12-02 | 中国航天空气动力技术研究院 | 一种cfd数值模拟与风洞试验相结合的气动预测方法 |
Non-Patent Citations (2)
| Title |
|---|
| 李勇: "基于Volterra级数的非定常气动力和气动弹性分析", 《中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑（月刊）》 * |
| 陈华霆 等: "振型叠加法合理振型数目的确定", 《应用力学学报》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110162822A (zh) * | 2019-03-19 | 2019-08-23 | 北京机电工程研究所 | 耦合结构模态的时域快速非定常气动力计算方法 |
| CN110705154A (zh) * | 2019-09-24 | 2020-01-17 | 中国航空工业集团公司西安飞机设计研究所 | 航空器开环气动伺服弹性系统模型均衡降阶的优选方法 |
| CN114818550A (zh) * | 2022-06-30 | 2022-07-29 | 中国飞机强度研究所 | 一种飞机振动试验中时变气动载荷地面等效模拟方法 |
| CN114818550B (zh) * | 2022-06-30 | 2022-09-13 | 中国飞机强度研究所 | 一种飞机振动试验中时变气动载荷地面等效模拟方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109459206B (zh) | 2020-10-27 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Silva et al. | Development of reduced-order models for aeroelastic analysis and flutter prediction using the CFL3Dv6. 0 code | |
| CN109459206B (zh) | 地面试验非定常气动力加载方法 | |
| Raveh | CFD-based models of aerodynamic gust response | |
| CN109933876B (zh) | 一种基于广义气动力的非定常气动力降阶方法 | |
| CN109726465B (zh) | 基于非结构曲边网格的三维无粘低速绕流的数值模拟方法 | |
| CN110610065A (zh) | 基于混合动网格技术的飞行器多体分离cfd仿真方法及系统 | |
| Chen et al. | Parametric reduced-order modeling of unsteady aerodynamics for hypersonic vehicles | |
| CN108763683A (zh) | 一种三角函数框架下新weno格式构造方法 | |
| CN115422654B (zh) | 基于cfd/csd技术的跨/超声速飞行器高效高精度颤振时域分析方法 | |
| Valente et al. | A doublet-lattice method correction approach for high fidelity gust loads analysis | |
| Timme et al. | Linear reduced order modelling for gust response analysis using the DLR-TAU code | |
| CN109753690A (zh) | 基于计算流体力学的非线性非定常气动力降阶方法 | |
| CN105205293A (zh) | 用于获得飞机部件气动载荷的方法和系统 | |
| He et al. | A Coupled Newton--Krylov Time Spectral Solver for Flutter Prediction | |
| Ritter et al. | Comparison of nonlinear aeroelastic methods for maneuver simulation of very flexible aircraft | |
| CN115329643A (zh) | 适用于分析叶片-楔形缘板阻尼器的非线性接触模型构建方法及系统 | |
| Wu et al. | Computation of frequency responses and their sensitivities for undamped systems | |
| CN116151156A (zh) | 适用任意气动参数和模态振型的气动力降阶建模方法 | |
| Saldía et al. | Numerical simulations using tvd schemes of two-dimensional supersonic flow in chemical equilibrium | |
| CN104809301B (zh) | 一种反映非线性刚柔混合连接特性的时域子结构方法 | |
| Rozov et al. | Antisymmetric boundary condition for small disturbance CFD | |
| Huang et al. | Correlation studies of geometrically nonlinear aeroelastic formulations with beam and shell elements | |
| Li et al. | Wind-resistant analysis for large-span single-layer reticulated shells | |
| CN112100900B (zh) | 一种空间非合作目标点云初始姿态测量方法 | |
| CN111008492B (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 | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20201027 Termination date: 20211217 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |