CN113806866B - 基于光线偏折量的高速飞行器探测方法及仿真方法 - Google Patents
基于光线偏折量的高速飞行器探测方法及仿真方法 Download PDFInfo
- Publication number
- CN113806866B CN113806866B CN202111069676.2A CN202111069676A CN113806866B CN 113806866 B CN113806866 B CN 113806866B CN 202111069676 A CN202111069676 A CN 202111069676A CN 113806866 B CN113806866 B CN 113806866B
- Authority
- CN
- China
- Prior art keywords
- laser
- simulation
- aircraft
- speed
- refractive index
- 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.)
- Active
Links
- 238000004088 simulation Methods 0.000 title claims abstract description 34
- 238000000034 method Methods 0.000 title claims abstract description 20
- 238000001514 detection method Methods 0.000 title claims abstract description 19
- 238000012544 monitoring process Methods 0.000 claims abstract description 16
- 238000004364 calculation method Methods 0.000 claims description 12
- 230000005540 biological transmission Effects 0.000 claims description 7
- 238000013461 design Methods 0.000 description 3
- 238000004422 calculation algorithm Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 229920002430 Fibre-reinforced plastic Polymers 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000000576 coating method Methods 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000011151 fibre-reinforced plastic Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- 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/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Mathematical Optimization (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Analysis (AREA)
- Evolutionary Computation (AREA)
- Pure & Applied Mathematics (AREA)
- Computer Hardware Design (AREA)
- Aviation & Aerospace Engineering (AREA)
- Computational Mathematics (AREA)
- Automation & Control Theory (AREA)
- Algebra (AREA)
- Computing Systems (AREA)
- Fluid Mechanics (AREA)
- Mathematical Physics (AREA)
- Optical Radar Systems And Details Thereof (AREA)
Abstract
一种基于光线偏折量的高速飞行器探测方法及仿真方法,处于临近空间的飞行平台向目标识别区域发射激光，并记录其出射光线方向，发送至地面监测站。地面监测站接受激光信号，记录激光的入射方向，得到激光入射方向与出射方向的差值，得到激光光线偏折量；通过分析激光光线偏折量和当地区域的大气特性，得到目标识别区域是否存在高速飞行器的结论。本发明使用光线探测维持了原有的探测效率。
Description
技术领域
本发明属于高速飞行器的探测技术领域，具体是一种基于光线偏折量的高速飞行器探测方法及仿真方法。
背景技术
在高速飞行器不断发展下，新型高速武器常使用隐身外形设计及隐身涂层实现对传统雷达隐蔽性。
目前，通常使用新型算法或加大雷达信号来实现对高速飞行器的可探测性，但是效果十分有限。隐身外形设计及隐身材料的迅速发展，使得传统的雷达原理及探测方式对高速飞行器的探测效果不佳。
与专利CN 104199045 A相比，该方法不依赖于尾焰特性，对于不使用发动机的高动能飞行器，或在高速飞行器尾焰被隐藏时，依然可以对高速飞行器进行探测。
发明内容
为克服上述现有技术的不足，本发明利用高速飞行器飞行时产生的明显的大气扰动，引起大气折射率的变化，通过观测光线穿越目标区域的偏折大小，对目标区域的高速飞行器进行探测。
本发明的技术解决方案是：
一方面，本发明提供一种基于光线偏折量的高速飞行器探测方法，其特点在于，该方法包括如下步骤：
步骤S1：在临近空间区域设置飞行平台，在目标识别地面区域设置地面监测站；
步骤S2：所述的飞行平台向地面监测站发射指定波长λ的激光，并将激光的发射坐标r0及发射方向t0通过无线电信号传播给地面监测站；
步骤S3：地面监测站接收激光入射信号，包括其入射坐标r∞及入射方向t∞，同时接收激光于飞行平台的出射信息，无线电信号，并记录激光的入射方向t0；
步骤S4：将激光入射方向与出射方向进行差值，得到激光光线偏折量t∞-t0，分析该激光光线偏折量，当激光光线偏折角度大于0.0001rad时，则认为目标识别区域内存在高速飞行器，否则，则不存在高速飞行器。
所述的临近空间区域是指距离海平面100km-300km高。
另一方面，本发明还提供一种基于光线偏折量的高速飞行器探测仿真方法，其特点在于，该方法包括如下步骤：
步骤S1：选取一高超声速飞行器作为仿真对象，建立流场网格，并在高超声速飞行器的尾部保留大的流场区域；
步骤S2：根据预定工况确定高超声速飞行器仿真高度及仿真速度，将高超声速飞行器仿真高度输入至大气模型中，其特点是一大气数据库，通过输入当地高度索引当地大气环境特性，包括当地温度，当地压力，当地大气密度；
步骤S3：利用流体力学求解器，根据工况将仿真速度及当地位置的温度、压力和大气密度，获得的飞行器的速度场、压力场、密度场和温度场；
步骤S4：计算海平面高度下光线折射率，公式如下：
其中，λ代表仿真模拟的平台出射激光的波长，n0,760代表海平面高度的大气特性下该波长的光线折射率；
计算目标识别区域内的大气折射率nT,P，公式如下：
其中，T为目标识别区域内温度，P为目标识别区域内压力P；
步骤S5：通过梯度算子计算折射率梯度场，公式如下：
其中，g代表某一标量的值，i,j代表在正交的规则网格中，其在x，y方向上的坐标；
步骤S6：对光线传输方程进行离散，得到可迭代求解的光线传输方程：
其中rk代表在第k步时，激光光线在计算域中的坐标，tk代表在第k步时，激光光线的传播方向，n代表在rk坐标下，该区域的折射率大小；通过确定仿真中模拟的平台激光出射方向t0，出射坐标r0，其余变量为龙格库塔迭代中的过程变量；
通过确定仿真中模拟的平台激光出射方向t0，出射坐标r0。可对计算域内光线传播路径进行求解，当激光穿越整个计算域时，记录其出射方向t∞，即模拟地面监测站接收到的激光入射角度。其差值即tdiff＝|t∞-t0|，至此整个仿真过程结束。
与现有技术相比本发明有益效果为：
(1)本发明考虑了高超声速飞行器其对尾部空气产生显著扰动，导致尾部空气折射率的变化明显，因此激光穿越该流场会发生相对显著的偏折。这种扰动是由于激波的特性而纯在的，因此不易被常用隐身设计及隐身技术干扰。
(2)本发明考虑了光路的传播特性，使用光线探测保证了新的探测方式维持的原有的探测效率。
(3)本发明中给出探测手段及仿真方法，可以进行仿真测试验证该探测手段对不同型号的飞行器的有效性。
附图说明
图1为本发明方法流程图
图2为本发明实施过程的具体示意图
具体实施方式
下面结合附图及实例对本发明做详细说明。
一种光线偏折量的高速飞行器探测方式及仿真方法，该探测方式的具体实施步骤为：
步骤一：在临近空间区域(距离海平面高度100km-300km距离)设置飞行平台，其向目标区域发射波长为λ的激光光束，并将其出射坐标r0及出射方向t0信息通过无线电磁波，发送于地面监测站；
步骤二：在地面监测站接受激光入射的信息，包括其入射坐标r∞及入射方向t∞，同时地面监测站接收激光信号，记录激光的入射方向t0；
步骤三：通过对激光的光线偏折量t∞-t0进行分析，当其大小大于0.0001rad时确定目标区域是否存在高速飞行器。
其同时有仿真方法，具体实施步骤为
选取一高速飞行器作为仿真对象，建立其流场网格，其特征为在高速飞行器尾部保留较大的流场区域，以便进行后续计算。
根据预定工况确定飞行器仿真高度及仿真速度，将飞行器仿真高度输入至Nrlmsise-00模型中，其特点是一大气数据库，通过输入当地高度索引当地大气环境特性，包括当地温度，当地压力，当地大气密度。
为对光路通过飞行器的偏折特性进行仿真，需要对仿真计算域的折射率进行计算。通过海平面折射率计算公式：
其中λ代表仿真模拟的平台出射激光的波长大小，n0,760代表海平面高度的大气特性下其对该波长的光线折射率大小。因此，在确定激光波长的情况下，可以得到海平面高度下大气的折射率特性。
通过折射率换算公式：
其中nT,P代表在当地温度为T，当地压力为P的情况下，当地大气的折射率与海平面大气折射率的关系。通过以上公式，在已知当地温度及当地压力的情况下，当地大气折射率的大小。将计算域应用该公式进行计算，可计算出计算区域内的折射率场。对计算域内的折射率梯度场进行求解，由于在模拟光线穿越计算域时，需要使用光线传输方程，其要求对折射率梯度场进行计算。通过梯度算子：
其中，g代表某一标量的值，i,j代表在正交的规则网格中，其在x，y方向上的坐标。该公式可对任意物理标量的梯度进行求解，并取得良好精度。通过以上方法，可求解折射率场的梯度。
对光线传输方程进行离散，得到可迭代求解的光线传输方程：
其中rk代表在第k步时，激光光线在计算域中的坐标，tk代表在第k步时，激光光线的传播方向，n代表在rk坐标下，该区域的折射率大小，公式中的其余参数代表龙格库塔算法的过程变量，无实际物理意义。
通过确定仿真中模拟的平台激光出射方向t0，出射坐标r0。可对计算域内光线传播路径进行求解，当激光穿越整个计算域时，记录其出射方向t∞，即模拟地面监测站接收到的激光入射角度。其差值即tdiff＝|t∞-t0|。至此整个仿真过程结束。
本发明结果可以对飞行平台扫描区域内的空气介质折射率进行检测，进而发现目标区域内是否存在高速飞行器。其原理直接，特征明显，具有良好的可检特性，弥补了当前雷达对隐身设计的速飞行器的检测性不佳的不足。
本发明未详细说明部分属于本领域技术人员公知常识。
Claims (3)
1.一种基于光线偏折量的高速飞行器探测方法，其特征在于，该方法包括如下步骤：
步骤S1：在临近空间区域设置飞行平台，在目标识别地面区域设置地面监测站；
步骤S2：所述的飞行平台向地面监测站发射指定波长的激光，并将激光的发射坐标及发射方向通过无线电信号传播给地面监测站；
步骤S3：地面监测站接收激光信号，并记录激光的入射方向，同时接收激光于飞行平台的出射信息，无线电信号；
步骤S4：将激光入射方向与出射方向进行差值，得到激光光线偏折量，分析该激光光线偏折量，当激光光线偏折角度大于0.0001rad时，则认为目标识别区域内存在高速飞行器，否则，则不存在高速飞行器。
2.根据权利要求1所述的基于光线偏折量的高速飞行器探测方法，其特征在于，所述的临近空间区域是指距离海平面100km-300km高。
3.一种基于光线偏折量的高速飞行器探测仿真方法，其特征在于，该方法包括如下步骤：
步骤S1：选取一高超声速飞行器作为仿真对象，建立流场网格，并在高超声速飞行器的尾部保留大的流场区域；
步骤S2：根据预定工况确定高超声速飞行器仿真高度及仿真速度，将高超声速飞行器仿真高度输入至大气模型中，其特点是一大气数据库，通过输入当地高度索引当地大气环境特性，包括当地温度，当地压力，当地大气密度；
步骤S3：利用流体力学求解器，根据工况将仿真速度及当地位置的温度、压力和大气密度，获得的飞行器的速度场、压力场、密度场和温度场；
步骤S4：计算海平面折射率，公式如下：
其中，λ代表仿真模拟的平台出射激光的波长，n0,760代表海平面高度的大气特性下该波长的光线折射率；
计算目标识别区域内的大气折射率nT,P，公式如下：
其中，T为目标识别区域内温度，P为目标识别区域内压力P；
步骤S5：通过梯度算子计算折射率梯度场，公式如下：
其中，g代表某一标量的值，i,j代表在正交的规则网格中，其在x，y方向上的坐标；
步骤S6：对光线传输方程进行离散，得到可迭代求解的光线传输方程：
其中rk代表在第k步时，激光光线在计算域中的坐标，tk代表在第k步时，激光光线的传播方向，n代表在rk坐标下，该区域的折射率大小；通过确定仿真中模拟的平台激光出射方向t0，出射坐标r0，其余变量为龙格库塔迭代中的过程变量；当激光穿越整个计算域时，记录其出射方向t∞，即地面监测站接收到的激光入射角度；
计算激光入射方向与出射方向进行差值，即tdiff＝|t∞-t0|。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111069676.2A CN113806866B (zh) | 2021-09-13 | 2021-09-13 | 基于光线偏折量的高速飞行器探测方法及仿真方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111069676.2A CN113806866B (zh) | 2021-09-13 | 2021-09-13 | 基于光线偏折量的高速飞行器探测方法及仿真方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113806866A CN113806866A (zh) | 2021-12-17 |
| CN113806866B true CN113806866B (zh) | 2024-03-26 |
Family
ID=78941086
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111069676.2A Active CN113806866B (zh) | 2021-09-13 | 2021-09-13 | 基于光线偏折量的高速飞行器探测方法及仿真方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113806866B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117030013B (zh) * | 2023-10-09 | 2024-04-05 | 中国科学院西安光学精密机械研究所 | 一种光束大气直线传输补偿方法及系统 |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4050067A (en) * | 1976-04-21 | 1977-09-20 | Elmore Jr Ethelbert P | Airborne microwave path modeling system |
| KR102089539B1 (ko) * | 2019-06-24 | 2020-03-16 | 한화시스템 주식회사 | 고층 기상에 따른 장거리 탐지 레이다의 전파 특성 및 표적 탐지 성능 분석 방법 및 시스템 |
| CN111060899A (zh) * | 2019-12-31 | 2020-04-24 | 武汉大学 | 星地一体化激光雷达回波波形仿真方法及系统 |
| CN111125869A (zh) * | 2019-11-11 | 2020-05-08 | 北京空间机电研究所 | 一种动目标大气扰动特性仿真方法 |
-
2021
- 2021-09-13 CN CN202111069676.2A patent/CN113806866B/zh active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4050067A (en) * | 1976-04-21 | 1977-09-20 | Elmore Jr Ethelbert P | Airborne microwave path modeling system |
| KR102089539B1 (ko) * | 2019-06-24 | 2020-03-16 | 한화시스템 주식회사 | 고층 기상에 따른 장거리 탐지 레이다의 전파 특성 및 표적 탐지 성능 분석 방법 및 시스템 |
| CN111125869A (zh) * | 2019-11-11 | 2020-05-08 | 北京空间机电研究所 | 一种动目标大气扰动特性仿真方法 |
| CN111060899A (zh) * | 2019-12-31 | 2020-04-24 | 武汉大学 | 星地一体化激光雷达回波波形仿真方法及系统 |
Non-Patent Citations (1)
| Title |
|---|
| 航空平台间激光通信捕获链路功率分析与仿真;徐春凤;韩成;姜会林;兵工学报;20161231;第37卷(第11期);全文 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113806866A (zh) | 2021-12-17 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110095784B (zh) | 一种复杂环境影响下的海洋-低层大气激光传输建模方法 | |
| CN109900256B (zh) | 一种自适应海洋移动声层析系统和方法 | |
| Le Pichon et al. | Acoustic propagation and atmosphere characteristics derived from infrasonic waves generated by the Concorde | |
| CN107577241B (zh) | 一种基于避障系统的消防无人机航迹规划方法 | |
| CN110045341A (zh) | 雷达高度表低截获性能测试方法 | |
| Raffel et al. | Background-oriented schlieren imaging for full-scale and in-flight testing | |
| CN113806866B (zh) | 基于光线偏折量的高速飞行器探测方法及仿真方法 | |
| CN110880071A (zh) | 基于算子的无源雷达作战效能柔性评估建模方法 | |
| CN114966899B (zh) | 一种基于多源多元遥感技术协同的区域能见度预测方法 | |
| CN102768354A (zh) | 一种用于获得水下目标回波数据的方法及系统 | |
| CN106772386A (zh) | 一种利用lpso算法由雷达回波反演大气波导方法 | |
| CN110031801A (zh) | 一种高速临近空间飞行器的探测方法及系统 | |
| Greenwood II | Fundamental rotorcraft acoustic modeling from experiments (frame) | |
| Yunus et al. | Toward inclusion of atmospheric effects in the aircraft community noise predictions | |
| CN104101869B (zh) | 一种极坐标下的地波雷达运动目标仿真建模方法 | |
| CN209485353U (zh) | 一种自适应海洋移动声层析系统 | |
| Harris et al. | Aircraft wake vortices: a comparison of wind-tunnel data with field trial measurements by laser radar | |
| Smith et al. | Development of Methods for Top-Down Methane Emission Measurements of Oil and Gas Facilities in an Offshore Environment Using a Miniature Methane Spectrometer and Long-Endurance UAS | |
| CN116660831A (zh) | 一种基于海面杂波信号监测数据的大气波导反演方法 | |
| JP3696135B2 (ja) | 乱気流予測装置及び乱気流予測方法 | |
| Jensen et al. | Computational Evaluations of Experimental Data for Sonic Boom Models with Nozzle Jet Flow Interactions | |
| Alkandari et al. | Optimization of Visible Light Positioning in Industrial Applications using Machine Learning | |
| CN114136324A (zh) | 一种基于雷达探测概率云图的隐身飞机航迹规划方法 | |
| CN117574730B (zh) | 一种超声速民机机动飞行聚焦声爆的数值预测方法 | |
| Rainer | In-situ observation of undisturbed surface layer scaler profiles for characterizing evaporative duct properties |
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 |