跨超声速大迎角配平翼的抖振响应获取方法及相关装置与流程

1.本技术涉及数据处理技术领域，具体涉及一种跨超声速大迎角配平翼的抖振响应获取方法及相关装置。


背景技术：

2.跨超声速飞行中，由非定常流动不稳定所引起的振动，形成跨超声速抖振现象。跨超声速抖振具有复杂的非定常、非线性特性，其引起的脉动载荷容易造成结构疲劳，从而引发飞行事故，极大地制约了飞行器的飞行包线。因而跨超声速抖振成为飞行器研发领域的难点和热点。跨超声速抖振的研究内容主要是：抖振的机理研究，抖振边界的判定，抖振载荷分析和抖振抑制等。由于跨超声速抖振的复杂性，早期的研究工作主要依赖风洞试验，且随着现代试验测量技术的不断进步和完善，试验研究方法在揭示抖振机理等领域起到了重要的作用。同时，随着计算流体力学和计算机技术的发展，数值模拟方法在抖振研究中得到了广泛的应用，且相对于风洞试验，数值模拟方法的成本微乎其微，因而，数值方法正逐步成为跨超声速抖振研究的主要手段。
3.抖振最早是由humphreys
1.于1951年的试验中发现的，而数值方法应用于抖振研究主要集中在近二十年。brunet
2.和grossi
3.等采用zone-des方法模拟了cat3d翼身组合体的抖振流动，并成功捕捉到了大部分的失稳区和分离区。sartor和timme
[4-5]
开展的研究表明，urans结合适当的湍流模型也能合理的复现抖振的主要流动特征，并指出s-a模型和剪切应力输运k-w模型预测较准确。ohmichi
[6]
、kenway
[7]
和ishida
[8]
等也基于urans方法开展了crm机翼的抖振仿真，对抖振机理进行了探讨。国内蔻家庆
[9]
通过本征正交分解和动态模态分解对ota15a翼型的跨声速抖振进行了分析。高传强
[10]
通过urans和动态模态分解方法研究了crm等典型机翼的抖振流动特性和主要失稳模态。
[0004]
目前研究存在的主要问题：目前的研究认为抖振问题的流过耦合特性较弱，气动载荷和结构运动无关，因此都是基于刚体模型进行的研究，导致了获取的抖振响应的准确性较低。


技术实现要素：

[0005]
本技术实施例提供一种跨超声速大迎角配平翼的抖振响应获取方法及相关装置，能够提升抖振响应获取时的准确性。
[0006]
本技术实施例的第一方面提供了一种跨超声速大迎角配平翼的抖振响应获取方法，所述方法包括：
[0007]
获取跨超声速大迎角配平翼的初始条件参数和边界条件参数；
[0008]
根据所述初始条件参数和边界条件参数，确定跨超声速大迎角配平翼的流场网格点的非定常气动载荷；
[0009]
将流场网格点的非定常气动载荷传递至结构网格点，以得到第一传递后的参数；
[0010]
根据跨超声速大迎角配平翼的结构动力学方程和所述第一传递后的参数，确定跨
超声速大迎角配平翼的结构网格点的位移信息；
[0011]
将所述结构网格点的位移信息传递至所述流场网格点，以得到第二传递后的参数；
[0012]
根据所述第二传递后的参数对流场网格进行调整，以得到调整后的流场网格，根据所述调整后的流场网格按照时间推进重复上述过程，以确定跨超声速大迎角配平翼的抖振响应。
[0013]
结合第一方面，在一个可能的实现方式中，所述根据所述初始条件参数和边界条件参数，确定跨超声速大迎角配平翼的流场网格点的非定常气动载荷，包括：
[0014]
通过数值方法求解如下公式，确定所述非定常气动载荷：
[0015][0016]
其中q＝[ρ，ρu，ρv，ρw，e]
t
是守恒流动变量，ρ是密度，u，v，w分别为x，y，z方向的速度分量，e为总能。fc(q)是对流通量，fv(q)是粘性通量，代表单元体积的ω边界的外法向单位矢量。
[0017]
结合第一方面，在一个可能的实现方式中，所述根据跨超声速大迎角配平翼的结构动力学方程和所述第一传递后的参数，以确定跨超声速大迎角配平翼的结构网格点的位移信息，包括：
[0018]
所述位移信息包括广义位移，通过数值方法求解如下公式，确定结构网格点的位移信息：
[0019][0020][0021][0022]
其中，{q(t)}为广义位移，w(x，y，z；t)为真实位移，{f(t)}为广义的气动力，δfi(x，y，z；t)为真实载荷矢量，为结构模态。[m]、[d]、[k]分别为结构的广义质量、阻尼和刚度矩阵。
[0023]
结合第一方面，在一个可能的实现方式中，所述将流场网格点的非定常气动载荷传递至结构网格点，以得到第一传递后的参数，包括：
[0024]
获取插值矩阵；
[0025]
根据所述插值矩阵将流场网格点的非定常气动载荷传递至结构网格点，以得到第一传递后的参数。
[0026]
结合第一方面，在一个可能的实现方式中，所述将结构网格点的位移传递至流场网格点，以得到第二传递后的参数，包括：
[0027]
获取插值矩阵；
[0028]
根据所述插值矩阵将结构网格点的位移传递至流场网格点，以得到第二传递后的参数。
[0029]
本技术实施例的第二方面提供了一种跨超声速大迎角配平翼的抖振响应获取装
置，所述装置包括：
[0030]
获取单元，用于获取跨超声速大迎角配平翼的初始条件参数和边界条件参数；
[0031]
第一确定单元，用于根据所述初始条件参数和边界条件参数，确定跨超声速大迎角配平翼的流场网格点的非定常气动载荷；
[0032]
第一传递单元，用于将流场网格点的非定常气动载荷传递至结构网格点，以得到第一传递后的参数；
[0033]
第二确定单元，用于根据跨超声速大迎角配平翼的结构动力学方程和所述第一传递后的参数，确定跨超声速大迎角配平翼的结构网格点的位移信息；
[0034]
第二传递单元，用于将所述结构网格点的位移信息传递至所述流场网格点，以得到第二传递后的参数；
[0035]
第三确定单元，用于根据所述第二传递后的参数对流场网格进行调整，以得到调整后的流场网格，根据所述调整后的流场网格按照时间推进重复上述过程，以确定跨超声速大迎角配平翼的抖振响应。
[0036]
结合第二方面，在一个可能的实现方式中，所述第一确定单元用于：
[0037]
通过数值方法求解如下公式，确定所述非定常气动载荷：
[0038][0039]
其中q＝[ρ，ρu，ρv，ρw，e]
t
是守恒流动变量，ρ是密度，u，v，w分别为x，y，z方向的速度分量，e为总能。fc(q)是对流通量，fv(q)是粘性通量，代表单元体积的ω边界的外法向单位矢量。
[0040]
结合第二方面，在一个可能的实现方式中，第二确定单元用于：
[0041]
所述位移信息包括广义位移，通过数值方法求解如下公式，确定结构网格点的位移信息：
[0042][0043][0044][0045]
其中，{q(t)}为广义位移，w(x，y，z；t)为真实位移，{f(t)}为广义的气动力，δfi(x，y，z；t)为真实载荷矢量，为结构模态。[m]、[d]、[k]分别为结构的广义质量、阻尼和刚度矩阵。
[0046]
结合第二方面，在一个可能的实现方式中，所述第一传递单元用于：
[0047]
获取插值矩阵；
[0048]
根据所述插值矩阵将流场网格点的非定常气动载荷传递至结构网格点，以得到第一传递后的参数。
[0049]
结合第二方面，在一个可能的实现方式中，所述第二传递单元用于：
[0050]
获取插值矩阵；
[0051]
根据所述插值矩阵将结构网格点的位移传递至流场网格点，以得到第二传递后的
参数。
[0052]
本技术实施例的第三方面提供一种终端，包括处理器、输入设备、输出设备和存储器，所述处理器、输入设备、输出设备和存储器相互连接，其中，所述存储器用于存储计算机程序，所述计算机程序包括程序指令，所述处理器被配置用于调用所述程序指令，执行如本技术实施例第一方面中的步骤指令。
[0053]
本技术实施例的第四方面提供了一种计算机可读存储介质，其中，上述计算机可读存储介质存储用于电子数据交换的计算机程序，其中，上述计算机程序使得计算机执行如本技术实施例第一方面中所描述的部分或全部步骤。
[0054]
本技术实施例的第五方面提供了一种计算机程序产品，其中，上述计算机程序产品包括存储了计算机程序的非瞬时性计算机可读存储介质，上述计算机程序可操作来使计算机执行如本技术实施例第一方面中所描述的部分或全部步骤。该计算机程序产品可以为一个软件安装包。
[0055]
实施本技术实施例，至少具有如下有益效果：
[0056]
通过获取跨超声速大迎角配平翼的初始条件参数和边界条件参数，根据所述初始条件参数和边界条件参数，确定跨超声速大迎角配平翼的流场网格点的非定常气动载荷，将流场网格点的非定常气动载荷传递至结构网格点，以得到第一传递后的参数，根据跨超声速大迎角配平翼的结构动力学方程和所述第一传递后的参数，确定跨超声速大迎角配平翼的结构网格点的位移信息，将所述结构网格点的位移信息传递至所述流场网格点，以得到第二传递后的参数，根据所述第二传递后的参数对流场网格进行调整，以得到调整后的流场网格，根据所述调整后的流场网格按照时间推进重复上述过程，以确定跨超声速大迎角配平翼的抖振响应。
附图说明
[0057]
为了更清楚地说明本技术实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
[0058]
图1a为本技术中用于计算气动载荷所用的cfd方法中的节点中心型有限控制体示意图；
[0059]
图1b为本技术实施例提供了一种跨超声速大迎角配平翼的抖振响应获取方法的流程示意图；
[0060]
图2为本技术中用于计算跨超声速抖振的流程图；
[0061]
图3为本技术中用于抖振计算的着陆巡视器基本外形参数；
[0062]
图4为本技术中配平翼在着陆巡视器上展开状态的安装示意图；
[0063]
图5为本技术用于抖振计算的cfd网格示意图；
[0064]
图6为本技术中配平翼后加密区网格；
[0065]
图7至图11为本技术中用于抖振计算模型的前五阶模态及响应的固有频率；
[0066]
图12为马赫数1.2时抖振响应分析中翼尖位移曲线；
[0067]
图13为马赫数1.5时抖振响应分析中翼尖位移曲线；
[0068]
图14为马赫数2.0时抖振响应分析中翼尖位移曲线；
[0069]
图15为马赫数2.5时抖振响应分析中翼尖位移曲线；
[0070]
图16为马赫数3.0时抖振响应分析中翼尖位移曲线
[0071]
图17为本技术实施例提供的一种终端的结构示意图；
[0072]
图18为本技术实施例提供了一种跨超声速大迎角配平翼的抖振响应获取装置的结构示意图。
具体实施方式
[0073]
下面将结合本技术实施例中的附图，对本技术实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本技术一部分实施例，而不是全部的实施例。基于本技术中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本技术保护的范围。
[0074]
本技术的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别不同对象，而不是用于描述特定顺序。此外，术语“包括”和“具有”以及它们任何变形，意图在于覆盖不排他的包含。例如包含了一系列步骤或单元的过程、方法、系统、产品或设备没有限定于已列出的步骤或单元，而是可选地还包括没有列出的步骤或单元，或可选地还包括对于这些过程、方法、产品或设备固有的其他步骤或单元。
[0075]
在本技术中提及“实施例”意味着，结合实施例描述的特定特征、结构或特性可以包含在本技术的至少一个实施例中。在说明书中的各个位置出现该短语并不一定均是指相同的实施例，也不是与其它实施例互斥的独立的或备选的实施例。本领域技术人员显式地和隐式地理解的是，本技术所描述的实施例可以与其它实施例相结合。
[0076]
本技术实施例中抖振响应分析中的非定常气动力采用基于混合网格的navier-stokes方程求解器，湍流采用基于rans方法的k-w sst模型。控制方程在非结构网格上基于节点中心型有限体积法进行离散，控制体如图1a所示，空间离散格式采用hlle格式。时间推进格式选择使用lu-sgs隐式方法。空间离散和时间推进精度均达到二阶精度。
[0077]
请参阅图1b，图1b为本技术实施例提供了一种跨超声速大迎角配平翼的抖振响应获取方法的流程示意图。如图1b所示，
[0078]
101、获取跨超声速大迎角配平翼的初始条件参数和边界条件参数。
[0079]
所述初始条件参数包括速度分布、压强分布、温度分布等，所述边界条件参数包括边界面类型和数值赋值等。
[0080]
102、根据所述初始条件参数和边界条件参数，确定跨超声速大迎角配平翼的流场网格点的非定常气动载荷。
[0081]
根据计算流体力学来确定出非定常气动载荷。其中，初始条件参数只用于首个物理时间步。
[0082]
103、将流场网格点的非定常气动载荷传递至结构网格点，以得到第一传递后的参数。
[0083]
可以通过插值矩阵，基于力和力矩守恒的原则，来将流场网格点的非定常气动载荷传递至结构网格点，以得到第一传递后的参数。
[0084]
104、根据跨超声速大迎角配平翼的结构动力学方程和所述第一传递后的参数，确
定跨超声速大迎角配平翼的结构网格点的位移信息。
[0085]
可以根据结构网格点的结构模态等参量，来确定出位移信息。
[0086]
105、将所述结构网格点的位移信息传递至所述流场网格点，以得到第二传递后的参数。
[0087]
可以参照前述步骤103中传递的方法，以得到第二传递后的参数，此处不再赘述。
[0088]
106、根据所述第二传递后的参数对流场网格进行调整，以得到调整后的流场网格。
[0089]
基于调整后的流场网格再次重复步骤102-106的运算，直至确定出满足要求的广义位移时间曲线，从而得到跨超声速大迎角配平翼的抖振响应。
[0090]
本技术在确定抖振响应时，(1)更加关注涡的形成和涡的发展传递过程，所以抖振不但需要较高的边界层网格，还需要在飞行器附近较大空间内的网格同样保持较高的质量；(2)抖振响应的计算中的时间步长的选取和总的迭代步数，既要考虑流场本身大涡和小涡的频域特性，还要考虑结构本身的频域特性，所以需要更小的时间步长和更多的时间推进步数，也就决定了更大的总计算量；(3)抖振计算时需要关注结构在整个时间历程里的动态载荷情况以及响应情况。
[0091]
本示例中，通过获取跨超声速大迎角配平翼的初始条件参数和边界条件参数，根据所述初始条件参数和边界条件参数，确定跨超声速大迎角配平翼的流场网格点的非定常气动载荷，将流场网格点的非定常气动载荷传递至结构网格点，以得到第一传递后的参数，根据跨超声速大迎角配平翼的结构动力学方程和所述第一传递后的参数，确定跨超声速大迎角配平翼的结构网格点的位移信息，将所述结构网格点的位移信息传递至所述流场网格点，以得到第二传递后的参数，根据所述第二传递后的参数对流场网格进行调整，以得到调整后的流场网格，根据所述调整后的流场网格按照时间推进重复上述过程，以确定跨超声速大迎角配平翼的抖振响应，从而提升了抖振响应确定时的准确性。
[0092]
在一个可能的实现方式中，所述根据所述初始条件参数和边界条件参数，确定跨超声速大迎角配平翼的流场网格点的非定常气动载荷的方法，包括：
[0093]
通过数值方法求解如下公式，确定所述非定常气动载荷：
[0094][0095]
其中q＝[ρ,ρu,ρv,ρw,e]
t
是守恒流动变量，ρ是密度，u，v，w分别为x，y，z方向的速度分量，e为总能。fc(q)是对流通量，fv(q)是粘性通量，代表单元体积的ω边界的外法向单位矢量。
[0096]
在一个可能的实现方式中，一种可能的述根据跨超声速大迎角配平翼的结构动力学方程和所述第一传递后的参数，以确定跨超声速大迎角配平翼的结构网格点的位移信息的方法包括：
[0097]
通过数值方法求解如下公式，确定结构网格点的位移信息：
[0098][0099]
[0100][0101]
其中，{q(t)}为广义位移，w(x，y，z；t)为真实位移，{f(t)}为广义的气动力，δfi(x，y，z；t)为真实载荷矢量，为结构模态。[m]、[d]、[k]分别为结构的广义质量、阻尼和刚度矩阵。
[0102]
在一个可能的实现方式中，一种可能的将流场网格点的非定常气动载荷传递至结构网格点，以得到第一传递后的参数的方法，包括：
[0103]
a1、获取插值矩阵；
[0104]
a2、根据所述插值矩阵将流场网格点的非定常气动载荷传递至结构网格点，以得到第一传递后的参数。
[0105]
在一个可能的实现方式中，一种可能的将结构网格点的位移传递至结构网格点，以得到第二传递后的参数的方法，包括：
[0106]
b1、获取插值矩阵；
[0107]
b2、根据所述插值矩阵将结构网格点的位移传递至结构网格点，以得到第二传递后的参数。
[0108]
其中，插值矩阵可以通过经验值或历史数据确定。
[0109]
可以通过如下公式所示的方法，将流场网格点的非定常气动载荷传递至结构网格点，以得到第一传递后的参数：
[0110][0111][0112]
其中，设为结构网格点上的力，为引起的结构点位移，为气动网格点上的载荷,为引起的流场网格点位移。
[0113]
在一个具体的实施例中，本技术实施例还提供一种具体的抖振响应获取方法，具体如下：
[0114]
本技术实例中采用了基于紧耦合技术的cfd-csd分析方法，分析历程如图2所示，具体步骤如下：
[0115]
(a)将初始流场输入计算流体力学求解器；
[0116]
(b)采用计算流体力学求解器(cfd)得出非定常气动载荷，区别去颤振需要更高质量的cfd网格，同时需要采用des脱体涡模型模拟湍流；
[0117]
(c)将流体网格的气动载荷通过数据交换方法传递到结构网格点；
[0118]
(d)求解结构动力学方程(csd)计算结构点的位移，输出结构位移，区别于颤振，抖振需要整个时间历程的动态载荷来确定动态响应；
[0119]
(e)结构点位移通过数据交换方法传递到流体网格节点；
[0120]
(f)对气动计算域网格进行重整；
[0121]
(g)返回b)直至满足计算要求退出。
[0122]
上述过程是求解抖振过程的一个时间步，区别于颤振，抖振的时间步长更小。
[0123]
本技术实施例以图1、图所示的火星着陆器及安装在其上的展开状态的配平翼模型进行抖振响应的计算，cfd网格采用gridgen生成的混合网格，着陆巡视器物面为三角形单元，配平翼物面为四面形单元，如图所示。边界层设置20层网格，首层网格厚度为0.03，如图所示，着陆巡视器的边界层内为三棱柱单元，配平翼边界层内为六面体单元。边界层与远场间填充四面体单元，并通过金字塔单元过渡。配平翼后流场区域建立加密区进行加密，如图所示。总网格量约为1200万。
[0124]
本技术实例采用模态叠加法计算抖振响应，弹性结构为配平翼，并选择前五阶模态作为抖振响应分析的结构输入，包括两个一弯模态、一扭模态、二弯模态和二扭模态，如图至图所示，其中，图7示出了第一阶模态，面内一弯，25.22hz；图8示出了第二阶模态，面外一弯，28.44hz；图9示出了第三阶模态，一阶扭转，84.74hz；图10示出了第四阶模态，面外二弯，123.69hz；图11示出了第五阶模态，二阶扭转，199.89hz；
[0125]
本技术实例所进行的抖振响应分析针对真实的飞行状态，包括五个超声速飞行工况，包括马赫数1.2、1.5、2.0、2.5、3.0，每个马赫数条件下分别包括三个飞行攻角状态，包括-10
°
、0
°
、10
°
，具体飞行参数在表1、
[0126]
表2中列出。
[0127]
表1进入弹道状态及相应大气参数
[0128]
马赫数速度(m/s)密度(kg/m3)压力(pa)温度(k)动压(pa)1.2277.9738.47e-03352.225217.212327.161.5347.4667.25e-03301.468217.212437.542.0463.5675.31e-03221.062217.212570.332.5575.0883.93e-03161.039213.991649.223.0686.2123.19e-03129.184211.587749.89
[0129]
表2真实弹道计算状态
[0130]
研究内容马赫数攻角α(
°
)抖振响应分析1.2、1.5、2.0、2.5、3.0-10、0、 10
[0131]
本技术实例的抖振响应分析方法用于装有配平翼的火星着陆器抖振问题，并对上述各飞行工况进行抖振分析，并获得了各工况下翼尖位移曲线，如图至图所示，其中，图12示出了马赫1.2，翼尖位移；图13示出了马赫1.5，翼尖位移；图14示出了马赫2.0，翼尖位移；图15示出了马赫2.5，翼尖位移；图16示出了马赫3.0，翼尖位移。
[0132]
对于所研究的五个马赫数，-10
°
与10
°
之间的变形大小比例约为1.5。抖振计算采用era方法对气动阻尼进行辨识；对于衰减的响应，era方法可以辨识出正的气动阻尼，而对于负气动阻尼情况，如果仿真时间足够长，可以得到极限环振荡(lco)。
[0133]
与上述实施例一致的，请参阅图4，图4为本技术实施例提供的一种终端的结构示意图，如图所示，包括处理器、输入设备、输出设备和存储器，处理器、输入设备、输出设备和存储器相互连接，其中，所述存储器用于存储计算机程序，所述计算机程序包括程序指令，所述处理器被配置用于调用所述程序指令，上述程序包括用于执行以下步骤的指令；
[0134]
获取跨超声速大迎角配平翼的初始条件参数和边界条件参数；
[0135]
根据所述初始条件参数和边界条件参数，确定跨超声速大迎角配平翼的流场网格点的非定常气动载荷；
[0136]
将流场网格点的非定常气动载荷传递至结构网格点，以得到第一传递后的参数；
[0137]
根据跨超声速大迎角配平翼的结构动力学方程和所述第一传递后的参数，确定跨超声速大迎角配平翼的结构网格点的位移信息；
[0138]
将所述结构网格点的位移信息传递至所述流场网格点，以得到第二传递后的参数；
[0139]
根据所述第二传递后的参数对流场网格进行调整，以得到调整后的流场网格，根据所述流场网格按照时间推进重复上述过程，以确定跨超声速大迎角配平翼的抖振响应。
[0140]
上述主要从方法侧执行过程的角度对本技术实施例的方案进行了介绍。可以理解的是，终端为了实现上述功能，其包含了执行各个功能相应的硬件结构和/或软件模块。本领域技术人员应该很容易意识到，结合本文中所提供的实施例描述的各示例的单元及算法步骤，本技术能够以硬件或硬件和计算机软件的结合形式来实现。某个功能究竟以硬件还是计算机软件驱动硬件的方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用使用不同方法来实现所描述的功能，但是这种实现不应认为超出本技术的范围。
[0141]
本技术实施例可以根据上述方法示例对终端进行功能单元的划分，例如，可以对应各个功能划分各个功能单元，也可以将两个或两个以上的功能集成在一个处理单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件功能单元的形式实现。需要说明的是，本技术实施例中对单元的划分是示意性的，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式。
[0142]
与上述一致的，请参阅图18，图18为本技术实施例提供了一种跨超声速大迎角配平翼的抖振响应获取装置的结构示意图。如图18所示，所述装置包括：
[0143]
获取单元301，用于获取跨超声速大迎角配平翼的初始条件参数和边界条件参数；
[0144]
第一确定单元302，用于根据所述初始条件参数和边界条件参数，确定跨超声速大迎角配平翼的流场网格点的非定常气动载荷；
[0145]
第一传递单元303，用于将流场网格点的非定常气动载荷传递至结构网格点，以得到第一传递后的参数；
[0146]
第二确定单元304，用于根据跨超声速大迎角配平翼的结构动力学方程和所述第一传递后的参数，确定跨超声速大迎角配平翼的结构网格点的位移信息；
[0147]
第二传递单元305，用于将所述结构网格点的位移信息传递至所述流场网格点，以得到第二传递后的参数；
[0148]
第三确定单元306，用于根据所述第二传递后的参数对流场网格进行调整，以得到调整后的流场网格，根据所述调整后的流场网格按照时间推进重复上述过程，以确定跨超声速大迎角配平翼的抖振响应。
[0149]
一个可能的实现方式中，所述第一确定单元302用于：
[0150]
通过数值方法求解如下公式，确定所述非定常气动载荷：
[0151][0152]
其中，其中q＝[ρ,ρu,ρv,ρw,e]
t
是守恒流动变量，ρ是密度，u，v，w分别为x，y，z方向的速度分量，e为总能。fc(q)是对流通量，fv(q)是粘性通量，代表单元体积的ω边界
的外法向单位矢量。
[0153]
在一个可能的实现方式中，第二确定单元304用于：
[0154]
所述位移信息包括广义位移，通过数值方法求解如下公式，确定结构网格点的位移信息：
[0155][0156][0157][0158]
其中，{q(t)}为广义位移，w(x，y，z；t)为真实位移，{f(t)}为广义的气动力，δfi(x，y，z；t)为真实载荷矢量，为结构模态。[m]、[d]、[k]分别为结构的广义质量、阻尼和刚度矩阵。
[0159]
在一个可能的实现方式中，所述第一传递单元303用于：
[0160]
获取插值矩阵；
[0161]
根据所述插值矩阵将流场网格点的非定常气动载荷传递至结构网格点，以得到第一传递后的参数。
[0162]
在一个可能的实现方式中，所述第一传递单元305用于：
[0163]
获取插值矩阵；
[0164]
根据所述插值矩阵将结构网格点的位移信息至流场网格点，以得到第二传递后的参数。
[0165]
本技术实施例还提供一种计算机存储介质，其中，该计算机存储介质存储用于电子数据交换的计算机程序，该计算机程序使得计算机执行如上述方法实施例中记载的任何一种跨超声速大迎角配平翼的抖振响应获方法的部分或全部步骤。
[0166]
本技术实施例还提供一种计算机程序产品，所述计算机程序产品包括存储了计算机程序的非瞬时性计算机可读存储介质，该计算机程序使得计算机执行如上述方法实施例中记载的任何一种跨超声速大迎角配平翼的抖振响应获方法的部分或全部步骤。
[0167]
需要说明的是，对于前述的各方法实施例，为了简单描述，故将其都表述为一系列的动作组合，但是本领域技术人员应该知悉，本技术并不受所描述的动作顺序的限制，因为依据本技术，某些步骤可以采用其他顺序或者同时进行。其次，本领域技术人员也应该知悉，说明书中所描述的实施例均属于优选实施例，所涉及的动作和模块并不一定是本技术所必须的。
[0168]
在上述实施例中，对各个实施例的描述都各有侧重，某个实施例中没有详述的部分，可以参见其他实施例的相关描述。
[0169]
在本技术所提供的几个实施例中，应该理解到，所揭露的装置，可通过其它的方式实现。例如，以上所描述的装置实施例仅仅是示意性的，例如所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，例如多个单元或组件可以结合或者可以集成到另一个系统，或一些特征可以忽略，或不执行。另一点，所显示或讨论的相互之间的耦合或直接耦合或通信连接可以是通过一些接口，装置或单元的间接耦合或通信连接，可以是电性或其它的形式。
[0170]
所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部单元来实现本实施例方案的目的。
[0171]
另外，在申请明各个实施例中的各功能单元可以集成在一个处理单元中，也可以是各个单元单独物理存在，也可以两个或两个以上单元集成在一个单元中。上述集成的单元既可以采用硬件的形式实现，也可以采用软件程序模块的形式实现。
[0172]
所述集成的单元如果以软件程序模块的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储器中。基于这样的理解，本技术的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的全部或部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储器中，包括若干指令用以使得一台计算机设备(可为个人计算机、服务器或者网络设备等)执行本技术各个实施例所述方法的全部或部分步骤。而前述的存储器包括：u盘、只读存储器(read-only memory，rom)、随机存取存储器(random access memory，ram)、移动硬盘、磁碟或者光盘等各种可以存储程序代码的介质。
[0173]
本领域普通技术人员可以理解上述实施例的各种方法中的全部或部分步骤是可以通过程序来指令相关的硬件来完成，该程序可以存储于一计算机可读存储器中，存储器可以包括：闪存盘、只读存储器、随机存取器、磁盘或光盘等。
[0174]
以上对本技术实施例进行了详细介绍，本文中应用了具体个例对本技术的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本技术的方法及其核心思想；同时，对于本领域的一般技术人员，依据本技术的思想，在具体实施方式及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本技术的限制。