(19)中华人民共和国国家知识产权局
(12)发明专利申请
(10)申请公布号 (43)申请公布日 (21)申请号 201911097132.X
(22)申请日 2019.11.11
(71)申请人 北京空间机电研究所
地址 100076 北京市丰台区南大红门路1号
9201信箱5分箱
(72)发明人 张月　董士奎　高鹏　张学敏　
邓红艳　赵号　张琢　苏云　食品安全英语作文
郑国宪　
(74)专利代理机构 中国航天科技专利中心
梦见老鼠是什么意思
11009
代理人 张晓飞
(51)Int.Cl.
G06F  30/20(2020.01)
G06T  17/00(2006.01)
(54)发明名称一种动目标大气扰动特性仿真方法(57)摘要一种动目标大气扰动特性仿真方法，包括获取目标几何参数、建立目标三维模型、目标及流场求解区网格划分、设定流场计算条件、选定流场计算模型、获得目标大气扰动流场参数、计算流场折射率分布、根据Barron算子计算折射率梯度、利用
kevin costnerRunge -Kutta光线追迹计算扰动大气偏折传输路径、获得光线穿过扰动场的偏折角。本发明首次对地物光线穿过高速目标形成的大气扰动区，所产生的光线偏折进行仿真分析。本发明可用于大气层内飞行的任何目标的大气扰动仿真分析，可分析任意谱段的地物光线经过目标大气扰动区的偏折情况，也可分析不同谱段地物光线以不同入射角穿过目标的大气扰动区形成的偏折情况。本发明为目标大气扰动仿真分析提供一种新手段，
可应用于大气扰动场特性研究。权利要求书2页  说明书5页  附图2页CN 111125869 A 2020.05.08
C N  111125869
A
1.一种动目标大气扰动特性仿真方法，其特征在于步骤如下：
1)获取目标几何参数；
2)获取目标三维模型；
3)目标及流场求解区网格划分；
4)设定流场计算条件；
周秀娜男友5)选定流场计算模型；
6)获得目标大气扰动流场参数；
7)计算得到流场折射率分布；
8)根据Barron算子计算折射率梯度；
9)利用Runge-Kutta光线追迹计算扰动大气偏折传输路径；
10)获得光线穿过扰动场的偏折角。
2.根据权利要求1所述的一种动目标大气扰动特性仿真方法，其特征在于：所述步骤1)的具体过程为：根据目标实物进行三维扫描，以获取目标几何参数；或根据目标缩比实物模型进行三维扫描，以获取目标几何参数；或查文献、公开报道中目标几何参数或图片，进行目标几何参数推演。
3.根据权利要求1所述的一种动目标大气扰动特性仿真方法，其特征在于：所述步骤2)的具体过程为：根据步骤1)获取的目标几何参数构建目标三维模型；或直接获取目标三维模型；所述目标三维模型为目标的外观三维模型，构建重点为目标外形特性，不涉及目标内部结构。
4.根据权利要求1所述的一种动目标大气扰动特性仿真方法，其特征在于：所述步骤3)的具体过程为：利用流场网格制作软件对目标三维模型以及目标大气扰动场区域进行网格划分，并根据目标三维外形细节特征、非细节特征以及目标运动过程中大气扰动突变区域、缓变区域特性对网格进行紧密、稀疏以及结构化、非结构化区分。
5.根据权利要求1所述的一种动目标大气扰动特性仿真方法，其特征在于：所述步骤4)的具体过程为：设定目标飞行条件和环境大气条件；所述目标飞行条件包括在大气层内的目标飞行高度、飞行速度、飞行姿态；所述环境大气条件包括在大气层内以目标为中心四周一定求解区域内的标准大气环境、湍流大气环境。
6.根据权利要求1所述的一种动目标大气扰动特性仿真方法，其特征在于：所述步骤5)的具体过程为：将步骤3)完成的目标三维模型及大气扰动区网格读入到流场计算软件中，并选定流场计算模型、流场计算边界条件；其中湍流模型选择k-e模型、流场求解选择LHS和RHS，流场计算边界条件选择远场边界、气体和壁面边界。
7.根据权利要求1所述的一种动目标大气扰动特性仿真方法，其特征在于：所述步骤6)的具体过程为：对流场进行求解，获得目标大气扰动流场参数，包括流场内每一网格节点处的密度、压力、温度。
8.根据权利要求1所述的一种动目标大气扰动特性仿真方法，其特征在于：所述步骤7)的具体过程为：
利用步骤6)获得的目标大气扰动流场参数，通过格拉斯通-戴尔公式求解得到目标大气扰动场内每一网格节点处的折射率分布；
所述格拉斯通-戴尔公式
式中n为气体折射率，ρ为气体密度，K G-D＝f(λ,T,P)为格拉斯通-戴尔常数。K G-D可通过参考文献获得；
为提升折射率计算精度，将某一网格点P(x,y,z)周围M个点的折射率代入n p(x,y,z)，M 为正整数；求解该点更高精度的折射率梯度结果；重复步骤7)，以完成流场内所有网格点的求解，得到所有网格点更高精度的折射率梯度分布；
流场中任一点P(x,y,z)折射率计算公式
其中，n i(i＝1,2,…)为P点周围空间点上的折射率值，
为(x j,y j,z j)点与P点的距离。
9.根据权利要求1所述的一种动目标大气扰动特性仿真方法，其特征在于：所述步骤8)的具体过程为：利用步骤7)计算得到的流场折射率分布结果，求解某点P(x,y,z)的折射率梯度分布；将P(x,y,z)点周围N个点的折射率代入Barron算子法，求解该P(x,y,z)点的折射率梯度；重复步骤8)，以完成流场内所有网格点的求解，得到所有网格点的折射率梯度分布。
10.根据权利要求1所述的一种动目标大气扰动特性仿真方法，其特征在于：所述步骤9)的具体过程为：采用Runge-Kutta光线追迹方法，实时记录光线到达点坐标的方式，标记任意谱段地物光线以任意入射角进入大气扰动场的位置，利用步骤8)计算得到的折射率梯度以及光线偏折与折射率梯度关系式，获得光线在标记处偏折到下一点的位置；重复步骤9)，完成光线在整个流场内传输过程的求解，从而得到光线穿出整个大气扰动场时的出射角度。
11.根据权利要求1所述的一种动目标大气扰动特性仿真方法，其特征在于：所述步骤10)的具体过程为：求解步骤9)得到的光线穿出整个大气扰动场时的出射角度与光线进入大气扰动场的入射角度的差值，获得光线穿过扰动场的偏折角。
一种动目标大气扰动特性仿真方法
技术领域
[0001]本发明涉及一种动目标大气扰动特性仿真方法，是一种对大气层内目标高速运动引发的大气扰动特性仿真新方法，可仿真分析不同地物光线以不同入射角穿过目标大气扰动区形成的光线偏折情况，可为以光学手段实现大气扰动可视化观测提供设计输入，为大气扰动特性研究提供一种新的技术手段，可应用于高速运动目标气动外形优化、高速运动目标低分辨率大幅宽探测等，可适用于地物背景，突破气动特性必须在风洞内获取的局限性，大幅降低气动特性获取的难度及成本。
背景技术
[0002]大气层内高速运动目标高速飞行必将引起剧烈的大气扰动，大气扰动情况可直接反映高速运动目标的气动特性，但是大气扰动不可视。通常需在风洞中开展吹风试验以研究目标的气动特性，但是风洞尺寸有限，试验目标通常为缩比目标，不能完全真实反映目标全尺寸气动特性。而且，目标高速飞行引起的大气扰动范围很大，风洞出风口尺寸有限，不能有效模拟整个扰动区范围的大气来流。另外，目标与大气的相对运动状态不同，风洞中目标为静止状态、大气为运动状态，实际飞行中目标为运动状态、大气为近似静止状态。更为重要的是，地面风洞试验中流场可视化设备为纹影仪，受限于纹影仪口径问题，风洞试验仅能对部分区域的大气扰动进行可视化观测。所以，利用风洞试验对大气层内高速运动目标的大气扰动特性开展研究，具有很多局限性。
发明内容
[0003]本发明解决的技术问题是：克服现有技术的不足，提供一种动目标大气扰动特性仿真方法，解决现有风洞手段获取高速运动目标大气扰动特性过程中存在的问题。[0004]本发明的技术解决方案是：一种动目标大气扰动特性仿真方法，其特征在于步骤如下：
[0005]1)获取目标几何参数；
[0006]2)获取目标三维模型；
[0007]3)目标及流场求解区网格划分；
[0008]4)设定流场计算条件
[0009]5)选定流场计算模型；
[0010]6)获得目标大气扰动流场参数；
[0011]7)计算得到流场折射率分布；
[0012]8)根据Barron算子计算折射率梯度；
[0013]9)利用Runge-Kutta光线追迹计算扰动大气偏折传输路径；
[0014]10)获得光线穿过扰动场的偏折角。
[0015]所述步骤1)的具体过程为：根据目标实物进行三维扫描，以获取目标几何参数；或根据目标缩比实物模型进行三维扫描，以获取目标几何参数；或查文献、公开报道中目标
几何参数或图片，进行目标几何参数推演。
[0016]所述步骤2)的具体过程为：根据步骤1)获取的目标几何参数构建目标三维模型；或直接获取目标三维模型；所述目标三维模型为目标的外观三维模型，构建重点为目标外形特性，不涉及目标内部结构。
[0017]所述步骤3)的具体过程为：利用流场网格制作软件对目标三维模型以及目标大气扰动场区域进行网格划分，并根据目标三维外形细节特征、非细节特征以及目标运动过程中大气扰动突变区域、缓变区域特性对网格进行紧密、稀疏以及结构化、非结构化区分。[0018]所述步骤4)的具体过程为：设定目标飞行条件和环境大气条件；所述目标飞行条件包括在大气层内的目标飞行高度、飞行速度、飞行姿态；所述环境大气条件包括在大气层内以目标为中心四周一定求解区域内的标准大气环境、湍流大气环境。
[0019]所述步骤5)的具体过程为：将步骤3)完成的目标三维模型及大气扰动区网格读入到流场计算软件中，并选定流场计算模型、流场计算边界条件；其中湍流模型选择k-e模型、流场求解选择LHS和RHS，流场计算边界条件选择远场边界、气体和壁面边界。
[0020]所述步骤6)的具体过程为：对流场进行求解，获得目标大气扰动流场参数，包括：流场内每一网格节点处的密度、压力、温度。
[0021]所述步骤7)的具体过程为：利用步骤6)获得的目标大气扰动流场参数，通过格拉斯通-戴尔公式求解得到目标大气扰动场内每一网格节点处的折射率分布；
[0022]所述格拉斯通-戴尔公式
[0023]
[0024]式中n为气体折射率，ρ为气体密度，K G-D＝f(λ,T,P)为格拉斯通-戴尔常数。K G-D可通过参考文献获得；
[0025]为提升折射率计算精度，将某一网格点P(x,y,z)周围M个点的折射率代入n p(x,y, z)，M为正整数，求解该点更高精度的折射率梯度结果；重复步骤7)，以完成流场内所有网格点的求解，得到所有网格点更高精度的折射率梯度分布；汽车首付一般多少
[0026]流场中任一点P(x,y,z)折射率计算公式
[0027]
[0028]其中，n i(i＝1,2,…)为P点周围空间点上的折射率值，
为(x j,y j,z j)点与P点的距离。
八百到底为什么不能播了
[0029]所述步骤8)的具体过程为：利用步骤7)计算得到的流场折射率分布结果，求解某点P(x,y,z)的折射率梯度分布；将P(x,y,z)点周围N个点的折射率代入Barron算子法，求解该P(x,y,z)点的折射率梯度。重复步骤8)，以完成流场内所有网格点的求解，得到所有网格点的折射率梯度分布；