Multi-source cross-scale high-power microwave unmanned aerial vehicle damage simulation method and simulation system

Abstract

This invention discloses a multi-source, multi-scale, high-power microwave unmanned aerial vehicle (UAV) damage simulation method: 1. Constructing the angle-related transfer function of the key UAV system; 2. Constructing a multi-source synthetic equivalent source and selecting the equivalent source with the strongest coupling; 3. Using the equivalent source as excitation, combining the angle-related transfer function and the slot-cable path coupling function, calculating the induced voltage and current at the cable terminal, and simulating the back-door coupling effect; 4. Using the equivalent source as excitation, combining the angle-related transfer function, calculating the antenna port coupling voltage, and simulating the front-door coupling effect; 5. Performing a comprehensive damage assessment of the response based on preset damage criteria. This invention achieves efficient simulation of multi-source collaborative scenarios through a multi-source synthetic framework combining the equivalent source method and transfer function, while considering the influence of the UAV's flight state and supporting comparison and verification with local irradiation tests, ensuring the accuracy of cross-scale mapping and achieving efficient cross-scale coupling.

Description

Technical Field

[0001] This invention belongs to the field of electromagnetic compatibility and high-power microwave technology, specifically relating to a multi-source, multi-scale, high-power microwave UAV damage simulation method, and also to a multi-source, multi-scale, high-power microwave UAV damage simulation system that implements the above simulation method. Background Technology

[0002] High-power microwave sources can instantly interfere with, suppress, or even destroy the electronic systems of drones by directionally radiating extremely high-energy electromagnetic pulses.

[0003] However, the destructive effects of high-power microwaves on UAVs involve cross-scale physical processes ranging from macroscopic scales (meter-scale platforms) to microscopic scales (nanoscale devices), and existing simulation methods mainly suffer from the following technical problems:

[0004] (1) High complexity of multi-source collaborative simulation: In a scenario where multiple microwave sources irradiate simultaneously from different directions, the accurate evaluation of field superposition effect and angle diversity gain requires modeling each source separately. The computational load increases linearly with the number of sources, making it difficult to efficiently handle large-scale multi-source collaborative scenarios.

[0005] (2) The dynamic changes in flight status are difficult to characterize: The attitude of the UAV changes continuously during flight, which leads to a dynamic change in its relative positional relationship with each interference source. Existing methods are difficult to consider the influence of flight status on multi-source synthesis effect.

[0006] (3) Insufficient accuracy of cross-scale coupling modeling: The complete action chain from the external field to the internal device involves multiple spatial scales. Full-scale fine modeling is computationally infeasible, and existing simplification methods often lose key physical information.

[0007] (4) Disconnect between simulation and experimental verification: Traditional whole-machine irradiation test is costly and time-consuming, making it difficult to fully verify in the simulation stage, and lacking simulation methods that can be compared with local tests.

[0008] Therefore, there is an urgent need for a method to simulate the damage effects of unmanned aerial vehicles (UAVs) to solve the above problems. Summary of the Invention

[0009] The first objective of this invention is to provide a multi-source, multi-scale, high-power microwave UAV damage simulation method. By combining the equivalent source method with the transfer function into a multi-source synthesis framework, it achieves efficient simulation of multi-source collaborative scenarios. At the same time, it considers the influence of the UAV's flight state and supports comparison and verification with local irradiation tests, ensuring the accuracy of cross-scale mapping and achieving efficient cross-scale coupling.

[0010] The second objective of this invention is to provide a multi-source, multi-scale, high-power microwave UAV damage simulation system that performs the above-described simulation method.

[0011] The first technical solution adopted in this invention is: a multi-source, multi-scale, high-power microwave UAV damage simulation method, the specific method of which is as follows:

[0012] S1. Construct the angle-related transfer function of the key UAV system using the frequency sweeping method;

[0013] S2. Construct multi-source synthetic equivalent sources based on the UAV flight status, and select the equivalent sources with the strongest coupling;

[0014] S3. Using the selected equivalent source as excitation, construct the slot-cable path coupling function, combine the angle-related transfer function and the slot-cable path coupling function to calculate the induced voltage and current at the cable terminal, and perform backdoor coupling effect simulation.

[0015] S4. Using the selected equivalent source as excitation and combining the angle-related transfer function, calculate the antenna port coupling voltage and perform front door coupling effect simulation.

[0016] S5. Based on the preset damage criteria, perform a comprehensive damage assessment on the obtained response.

[0017] The invention is further characterized by:

[0018] A multi-source, multi-scale, high-power microwave UAV damage simulation method is presented below:

[0019] The angle-related transfer function for constructing the key UAV system in S1 specifically includes the following aspects:

[0020] S1.1 Constructing the transfer function of the aperture structure, the specific process is as follows: Establish a local geometric model containing typical apertures, use plane wave excitation, perform frequency sweep simulation within a preset frequency range and incident angle range, record the electric field response of key points inside the aperture, and obtain the transfer function of the aperture structure. :

[0021] (1);

[0022] in, For frequency, The incident elevation angle, It is the azimuth angle. For a strong incident field, The electric field response at key points inside the aperture;

[0023] The transfer function of the slot structure was constructed using the finite-difference time-domain method or the finite element method, with a frequency scanning range of 0.1~30 GHz and an elevation and azimuth scanning range of 0°~90°.

[0024] S1.2 Constructing the frequency-selective surface transfer function, the specific process is as follows: Establish a frequency-selective surface, obtain the variation of the transmission coefficient with frequency and incident angle through frequency sweep simulation, and obtain the frequency-selective surface transfer function. ;

[0025] S1.3 Constructing the antenna coupling transfer function, the specific process is as follows: Establish a single-element antenna model, obtain the relationship between the open-circuit voltage at the antenna port and the incident field strength through frequency sweep simulation, and obtain the antenna coupling transfer function. ;

[0026] S1.4 Identify the maximum coupling condition: Identify the frequency-angle combination with the strongest coupling from the transfer function data of the frequency-selective surface and antenna coupling, and use it as the focus of subsequent fine simulation.

[0027] The specific method for S2 is as follows:

[0028] S2.1, Set the drone's time. The position is Attitude angles include pitch, roll, and yaw; the positions of each high-power microwave source are set. And parameters such as peak power, operating frequency, pulse width, pulse repetition frequency, and transmit antenna gain;

[0029] S2.2 Identify key components: Based on the UAV structural analysis, identify at least one key component affected by the coupling effect, including heat dissipation holes, structural mounting seams, sensor windows, and antenna covers.

[0030] S2.3 Constructing a Local Equivalent Surface: Select a closed surface surrounding the critical part as a local equivalent surface. The local equivalent surface is a cuboid or cylindrical surface, with a size 3 to 5 times the feature size of the critical part, and a distance from the critical part. ,in The wavelength corresponding to the operating frequency of high-power microwaves;

[0031] S2.4 Calculate the total field distribution on the local equivalent surface: For each flight state, based on the geometric relationship between each high-power microwave source and the UAV, calculate the incident field strength of each equivalent source on the local equivalent surface, and perform vector superposition, as shown in the following formula:

[0032] (2);

[0033] in Overall strength It is the field strength of each high-power microwave source at the local equivalent surface, which is determined by the source's peak power, radiation pattern, propagation distance, transmitting antenna gain, and atmospheric attenuation; It includes the source initial phase and the propagation phase, and it is related to its operating frequency and initial phase; For the drone's position parameters, Position parameters for each high-power microwave source; i It is the serial number of the high-power microwave source. j It is the mathematical symbol for imaginary numbers;

[0034] S2.5 Constructing Local Equivalent Sources: Based on the total field distribution on the local equivalent surface, the equivalent current or magnetic current source distribution on the equivalent surface is obtained by solving the inverse problem, forming a multi-source synthetic equivalent source corresponding to the flight state;

[0035] S2.6 Filtering the maximum coupling equivalent source: For the multi-source synthetic equivalent sources corresponding to each flight state constructed in S2.5, the angle correlation transfer function constructed in S1 is used to quickly evaluate the internal coupling response of each equivalent source in the key parts, and filter out the equivalent source with the maximum coupling.

[0036] S2.7 Constructing an equivalent source at 1 / 4 wavelength: Shift the selected equivalent source outwards along the normal direction to a distance from the critical location. At this point, a new equivalent source is constructed as the input for subsequent detailed simulation.

[0037] The specific method for S3 is as follows:

[0038] S3.1 Calculation of the internal field of the pore: Using the equivalent source constructed in S2.7 as the excitation, and utilizing the pore structure transfer function constructed in S1. Calculate the internal field of the pore. Alternatively, a local full-wave simulation can be performed directly;

[0039] S3.2 Constructing the Through-Cable Path Coupling Function: Constructing the coupling function from the field behind the through-hole to the cable port through parametric simulation. ,in The distance from the center of the slot to the cable. The angle between the cable and the normal of the hole / slot. For cable length, For frequency;

[0040] The pinhole-cable path coupling function is constructed using a parametric scanning method, with a distance... The scanning range is 1 mm to 50 cm, with an included angle. The scanning range is 0°~90°, and the cable length is... The scanning range is 10cm to 1m. The simulation results are organized into a multidimensional interpolation table or fitted using a neural network.

[0041] S3.3 Cable Termination Response Calculation: Based on the actual cable layout parameters inside the UAV, the induced voltage or current at the cable termination is obtained by querying the slot-cable path coupling function constructed in S3.2.

[0042] (3);

[0043] in For the field inside the pore;

[0044] S3.4 Sensitive Device Injection Analysis: The induced voltage or current at the cable terminal obtained in S3.3 is used as an excitation source and injected into the subsequent circuit simulation model to analyze the response state of the sensitive device.

[0045] The specific method for S4 is as follows:

[0046] S4.1 Calculation of Frequency-Selective Surface Transmission Field: Using the equivalent source constructed in S2.7 as the excitation, and utilizing the frequency-selective surface transfer function constructed in S1... Calculate the synthesized field after passing through the frequency-selective surface. :

[0047] (4);

[0048] in For the first i The incident field vector of a high-power microwave source on the outer surface of a frequency-selective surface. For the first i A high-power microwave source polarization unit vector; i For drones and the i The azimuth angle formed between the high-power microwave sources, i.e. the rotation angle about the reference axis; For drones and the i The pitch angle formed between the high-power microwave sources, that is, the angle between the incident direction and the reference axis; The mathematical symbol for imaginary numbers;

[0049] S4.2 Antenna Port Coupling Voltage Calculation: Using the transmitted field inside the radome obtained in S4.1 as the excitation, the antenna coupling transfer function constructed in S1 is used... Calculate the antenna port coupling voltage :

[0050] (5);

[0051] The local elevation angle of the incident wave relative to the antenna mounting plane; The local azimuth angle of the incident wave relative to the antenna mounting plane;

[0052] S4.3 Simulation of the receiver front-end circuit: The antenna port coupling voltage obtained in S4.2 is... As an excitation source, the simulation model of the receiving front-end circuit is imported to analyze the response of the limiter and the low-noise amplifier.

[0053] S4.4 Device Temperature Rise and Damage Analysis: Based on the electrothermal analogy method, a thermal equivalent model of the low-noise amplifier is established. According to the pulse width and repetition frequency parameters, the device temperature rise is calculated to determine whether it exceeds the damage threshold.

[0054] Frequency-selective surface transfer function in S4.1 It can be replaced by existing radomes or frequency-selective surface pattern data;

[0055] If the surface pattern data only contains amplitude information, the phase information can be inverted using minimum phase approximation algorithms such as Hilbert transform and homomorphic filtering, or the phase information can be reconstructed based on a physical model that shows the transmission phase changing linearly with the angle.

[0056] The damage criteria in S5 include the following aspects:

[0057] Back door coupling damage judgment: The induced voltage and current of the cable terminal obtained in S3.4 are compared with the preset device voltage and current damage thresholds. If the thresholds are exceeded, permanent damage is determined to have occurred.

[0058] Front door coupling suppression judgment: The antenna port coupling voltage obtained from S4.2 is converted into the power density at the antenna and compared with the preset receiving density suppression threshold. If it exceeds the threshold, it is determined that functional suppression has occurred.

[0059] Front door coupling damage judgment: Compare the device temperature rise obtained in S4.4 with the preset device thermal damage threshold. If it exceeds the threshold, it is determined that permanent damage has occurred.

[0060] Multi-source synergistic cumulative effect assessment: For scenarios involving simultaneous or time-sharing irradiation from multiple sources, considering power superposition and energy accumulation effects, the cumulative criterion is applied to assess the overall damage probability.

[0061] (6);

[0062] in: This is a power superposition function, describing the time-varying behavior of multi-source pulses. t Total power; Let be the thermal response function, and be the thermal time constant of the device. t Exponentially decaying kernel; The damage energy threshold is expressed in joules.

[0063] The second technical solution adopted in this invention is: a multi-source, multi-scale, high-power microwave UAV damage simulation system, including a transfer function library construction module, a multi-source equivalent source construction module, a back-door coupling simulation module, a front-door coupling simulation module, and a damage assessment module.

[0064] The invention is further characterized by:

[0065] The transfer function library construction module is used to construct the angle-related transfer functions of key UAV systems using the frequency sweep method.

[0066] The multi-source equivalent source construction module is used to consider the flight state of the UAV, construct local equivalent surfaces near key parts, construct multi-source synthetic equivalent sources based on the total field distribution of each microwave source on the local equivalent surface, and select the equivalent source with the largest coupling.

[0067] The back-door coupling simulation module is used to calculate the induced voltage and current at the cable terminal by using the selected equivalent source as excitation and combining the transfer function of the slot structure and the slot-cable path coupling function.

[0068] The front-door coupling simulation module is used to calculate the antenna port coupling voltage and device response by using a selected equivalent source as excitation and combining the frequency-selective surface transfer function and the antenna coupling transfer function.

[0069] The damage assessment module is used to comprehensively assess the damage effect based on preset damage criteria.

[0070] The beneficial effects of this invention are:

[0071] (1) The multi-source cross-scale high-power microwave UAV damage simulation method of the present invention dynamically constructs multi-source synthetic equivalent sources based on the real-time flight state of the UAV, accurately characterizing the impact of the change in the relative position relationship between the UAV and each interference source on the damage effect during the combat process; by introducing parameterized modeling of flight state, a refined description of the UAV motion process is realized; by analyzing the damage probability under different flight attitudes, the UAV can be guided to avoid high-threat attitudes, or the deployment of interference sources can be optimized to cover the entire flight profile of the target, which can truly reflect the impact of flight state.

[0072] (2) The multi-source cross-scale high-power microwave UAV damage simulation method of the present invention realizes the rapid mapping of the complete transfer link from the external field to the internal cable and then to the device through the pre-construction of the transfer function and the hole-cable coupling function, avoiding multi-scale repeated simulation. Once the transfer function and coupling function are constructed, they can handle any combination of incident direction and frequency, greatly reducing the number of repeated simulations and improving the computational efficiency. The transfer function retains the resonance characteristics and phase information of the key structure, and the coupling function considers the geometric sensitivity of the cable layout, ensuring the accuracy of cross-scale mapping and realizing efficient cross-scale coupling.

[0073] (3) The multi-source, multi-scale, high-power microwave UAV damage simulation method of the present invention adopts a local equivalent surface strategy to enable direct comparison between the simulation and the local irradiation test, realize closed-loop verification between simulation and test, and improve the credibility of the simulation; compare the results of the local irradiation test with the simulation results based on the local equivalent source to verify the accuracy of the equivalent source construction and the effectiveness of the transfer function and coupling function; if the deviation is within an acceptable range (e.g., 3 dB), the simulation method is considered credible; if the deviation is large, the model parameters or the equivalent source construction method are corrected until they match.

[0074] (4) The multi-source cross-scale high-power microwave UAV damage simulation method of the present invention simplifies the multi-source collaborative problem into a single-source problem through the local equivalent surface strategy. The computational load is reduced inversely proportional to the cube of the number of sources. At the same time, only the local model needs to be processed instead of the whole model, which greatly improves the simulation efficiency. Attached Figure Description

[0075] Figure 1 This is an overall flowchart of the multi-source, multi-scale, high-power microwave UAV damage simulation method of the present invention;

[0076] Figure 2 This is a schematic diagram of the angle-related transfer function construction process in the simulation method of this invention;

[0077] Figure 3 This is a schematic diagram of the multi-source synthesis equivalent source construction process considering flight state in the simulation method of this invention;

[0078] Figure 4 This is a schematic diagram of the construction result of the multi-source synthesis equivalent source considering flight state in the simulation method of this invention;

[0079] Figure 5 This is a model diagram of the local equivalent source irradiation aperture in the simulation method of this invention;

[0080] Figure 6 This is a schematic diagram illustrating the construction of the aperture-cable path coupling function in the simulation method of this invention;

[0081] Figure 7 This is a schematic diagram of the comprehensive evaluation process of damage effects in the simulation method of this invention;

[0082] Figure 8 This is a schematic diagram of dual-antenna source synthesis in Embodiment 1 of the simulation method of the present invention;

[0083] Figure 9 This is a schematic diagram of the four-antenna source synthesis in Embodiment 1 of the simulation method of the present invention;

[0084] Figure 10 This is a schematic diagram of the equivalent dual-source result of a typical equivalent source in Embodiment 1 of the simulation method of the present invention;

[0085] Figure 11 This is a graph showing the evaluation results of electric field shielding effectiveness in Embodiment 1 of the simulation method of the present invention. Detailed Implementation

[0086] The present invention will now be described in detail with reference to the accompanying drawings and specific embodiments.

[0087] This invention relates to a multi-source, multi-scale, high-power microwave UAV damage simulation method, such as... Figure 1 As shown, the specific method is as follows:

[0088] S1. Construct a key system transfer function library: Construct angle-related transfer functions for key UAV systems using the frequency sweeping method;

[0089] Among them, such as Figure 2 As shown, the angle-dependent transfer function for constructing a key UAV system specifically includes the following aspects:

[0090] S1.1 Constructing the transfer function of the aperture structure, the specific process is as follows: Establish a local geometric model containing typical apertures, use plane wave excitation, perform frequency sweep simulation within a preset frequency range and incident angle range, record the electric field response of key points inside the aperture, and obtain the transfer function of the aperture structure. :

[0091] (1);

[0092] in, For frequency, The incident elevation angle, It is the azimuth angle. For a strong incident field, The electric field response at key points inside the aperture;

[0093] The transfer function of the slot structure was constructed using the finite-difference time-domain method or the finite element method, with a frequency scanning range of 0.1~30 GHz and an elevation and azimuth scanning range of 0°~90°.

[0094] S1.2 Constructing the frequency-selective surface transfer function, the specific process is as follows: Establish a frequency-selective surface, obtain the variation of the transmission coefficient with frequency and incident angle through frequency sweep simulation, and obtain the frequency-selective surface transfer function. ;

[0095] S1.3 Constructing the antenna coupling transfer function, the specific process is as follows: Establish a single-element antenna model, obtain the relationship between the open-circuit voltage at the antenna port and the incident field strength through frequency sweep simulation, and obtain the antenna coupling transfer function. ;

[0096] S1.4 Identify the maximum coupling condition: Identify the frequency-angle combination with the strongest coupling from the transfer function data of the frequency-selective surface and antenna coupling, and use it as the focus of subsequent fine simulation.

[0097] S2. Construct multi-source synthetic equivalent sources based on the UAV flight status, and select the equivalent sources with the strongest coupling;

[0098] like Figure 3 As shown, the specific method is as follows:

[0099] S2.1, Set the drone's time. The position is Attitude angles include pitch, roll, and yaw; the positions of each high-power microwave source are set. And parameters such as peak power, operating frequency, pulse width, pulse repetition frequency, and transmit antenna gain;

[0100] S2.2 Identify key components: Based on the UAV structural analysis, identify at least one key component affected by the coupling effect, including heat dissipation holes, structural mounting seams, sensor windows, and antenna covers.

[0101] S2.3 Constructing a Local Equivalent Surface: Select a closed surface surrounding the critical part as a local equivalent surface. The local equivalent surface is a cuboid or cylindrical surface, with a size 3 to 5 times the feature size of the critical part, and a distance from the critical part. ,in The wavelength corresponding to the operating frequency of high-power microwaves;

[0102] S2.4 Calculate the total field distribution on the local equivalent surface: For each flight state, based on the geometric relationship between each high-power microwave source and the UAV, calculate the incident field strength of each equivalent source on the local equivalent surface, and perform vector superposition, as shown in the following formula:

[0103] (2);

[0104] in Overall strength It is the field strength of each high-power microwave source at the local equivalent surface, which is determined by the source's peak power, radiation pattern, propagation distance, transmitting antenna gain, and atmospheric attenuation; It includes the source initial phase and the propagation phase, and it is related to its operating frequency and initial phase; For the drone's position parameters, Position parameters for each high-power microwave source; i It is the serial number of the high-power microwave source. j It is the mathematical symbol for imaginary numbers;

[0105] S2.5 Constructing Local Equivalent Sources: Based on the total field distribution on the local equivalent surface, the equivalent current or magnetic current source distribution on the equivalent surface is obtained by solving the inverse problem, forming a multi-source synthetic equivalent source corresponding to this flight state, such as... Figure 4 As shown;

[0106] S2.6 Filtering the maximum coupling equivalent source: For the multi-source synthetic equivalent sources corresponding to each flight state constructed in S2.5, the angle correlation transfer function constructed in S1 is used to quickly evaluate the internal coupling response of each equivalent source in the key parts, and filter out the equivalent source with the maximum coupling.

[0107] S2.7 Constructing an equivalent source (radiation source) at 1 / 4 wavelength: Shift the selected equivalent source outward along the normal direction to a distance from the critical location. At this point, a new equivalent source (radiation source) is constructed as the input for subsequent detailed simulation.

[0108] S3. Using the selected equivalent source as excitation, construct the slot-cable path coupling function. Combine the angle-related transfer function and the slot-cable path coupling function to calculate the induced voltage and current at the cable terminal, and perform backdoor coupling effect simulation. The specific method is as follows:

[0109] S3.1 Calculation of the field inside the aperture: Using the equivalent source (radiation source) constructed in S2.7 as the excitation, and utilizing the transfer function of the aperture structure constructed in S1. Calculate the internal field of the pore. Alternatively, a local full-wave simulation can be performed directly, such as a full-wave simulation model of a microwave equivalent source irradiating a hole in a cylinder, as shown in the example. Figure 5 As shown;

[0110] S3.2 Constructing the Through-Cable Path Coupling Function: Constructing the coupling function from the field behind the through-hole to the cable port through parametric simulation. ,in The distance from the center of the slot to the cable. The angle between the cable and the normal to the hole / slot. For cable length, For frequency;

[0111] like Figure 6 As shown, the pinhole-cable path coupling function is constructed using a parametric scanning method, with a distance... The scanning range is 1mm to 50cm, with an included angle. The scanning range is 0°~90°, and the cable length is... The scanning range is 10cm to 1m. The simulation results are organized into a multidimensional interpolation table or fitted using a neural network.

[0112] S3.3 Cable Termination Response Calculation: Based on the actual cable layout parameters inside the UAV, the induced voltage or current at the cable termination is obtained by querying the slot-cable path coupling function constructed in S3.2.

[0113] (3);

[0114] in For the field inside the pore;

[0115] S3.4 Sensitive Device Injection Analysis: The induced voltage or current at the cable terminal obtained in S3.3 is used as an excitation source and injected into the subsequent circuit simulation model to analyze the response state of the sensitive device.

[0116] S4. Using the selected equivalent source as excitation and combining the angle-dependent transfer function, calculate the antenna port coupling voltage and perform a front-door coupling effect simulation; the specific method is as follows:

[0117] S4.1 Calculation of Frequency-Selective Surface Transmission Field: Using the equivalent source constructed in S2.7 as the excitation, and utilizing the frequency-selective surface transfer function constructed in S1... Calculate the synthesized field after passing through the frequency-selective surface. :

[0118] (4);

[0119] in For the first i The incident field vector of a high-power microwave source on the outer surface of a frequency-selective surface. For the first i A high-power microwave source polarization unit vector; i For drones and the i The azimuth angle formed between the high-power microwave sources, i.e. the rotation angle about the reference axis; For drones and the i The pitch angle formed between the high-power microwave sources, that is, the angle between the incident direction and the reference axis; The mathematical symbol for imaginary numbers;

[0120] The frequency-selective surface transfer function It can be replaced by existing radomes or frequency-selective surface pattern data;

[0121] If the surface pattern data only contains amplitude information, the phase information can be inverted using minimum phase approximation algorithms such as Hilbert transform and homomorphic filtering, or the phase information can be reconstructed based on a physical model that shows the transmission phase changing linearly with the angle.

[0122] S4.2 Antenna Port Coupling Voltage Calculation: Using the radome transmission field obtained in S4.1 as the excitation, the antenna coupling transfer function constructed in S1 is used... Calculate the antenna port coupling voltage :

[0123] (5);

[0124] The local elevation angle of the incident wave relative to the antenna mounting plane; The local azimuth angle of the incident wave relative to the antenna mounting plane;

[0125] S4.3 Simulation of the receiver front-end circuit: The antenna port coupling voltage obtained in S4.2 is... As an excitation source, the simulation model of the receiving front-end circuit is imported to analyze the response of the limiter and the low-noise amplifier.

[0126] S4.4 Device Temperature Rise and Damage Analysis: Based on the electrothermal analogy method, a thermal equivalent model of the low-noise amplifier is established. According to the pulse width and repetition frequency parameters, the device temperature rise is calculated to determine whether it exceeds the damage threshold.

[0127] S5. Based on the preset damage criteria, perform a comprehensive damage assessment on the responses obtained from S3 and S4.

[0128] like Figure 7 As shown, the damage criteria include the following aspects:

[0129] Back door coupling damage judgment: The induced voltage and current of the cable terminal obtained in S3.4 are compared with the preset device voltage and current damage thresholds. If the thresholds are exceeded, permanent damage is determined to have occurred.

[0130] Front door coupling suppression judgment: The antenna port coupling voltage obtained from S4.2 is converted into the power density at the antenna and compared with the preset receiving density suppression threshold. If it exceeds the threshold, it is determined that functional suppression has occurred.

[0131] Front door coupling damage judgment: Compare the device temperature rise obtained in S4.4 with the preset device thermal damage threshold. If it exceeds the threshold, it is determined that permanent damage has occurred.

[0132] Multi-source synergistic cumulative effect assessment: For scenarios involving simultaneous or time-sharing irradiation from multiple sources, considering power superposition and energy accumulation effects, the cumulative criterion is applied to assess the overall damage probability.

[0133] (6);

[0134] in: It is an aggregated power function that describes the total power of multi-source pulses over time (including simultaneous superposition or time-division accumulation effects). The thermal response function is denoted by , and the thermal time constant of the device is denoted by . t Exponentially decaying kernels, for example ,in It is a constant, which is related to the heat accumulation and heat dissipation characteristics of the device, and can be obtained through actual measurement; The damage energy threshold is expressed in joules.

[0135] The method of the present invention also includes: local irradiation test verification, in which a horn antenna is used to irradiate a local structural component containing key parts in a microwave anechoic chamber, the internal response is measured, and the result is compared with the simulation result based on the local equivalent source to verify the accuracy of the equivalent source construction.

[0136] The present invention provides a multi-source, multi-scale, high-power microwave UAV damage simulation system for the above-mentioned simulation method, comprising a transfer function library construction module, a multi-source equivalent source construction module, a back-door coupled simulation module, a front-door coupled simulation module, and a damage assessment module.

[0137] The transfer function library construction module is used to construct the angle-related transfer functions of key UAV systems using the frequency sweep method.

[0138] The multi-source equivalent source construction module is used to consider the flight state of the UAV, construct local equivalent surfaces near key parts, construct multi-source synthetic equivalent sources based on the total field distribution of each microwave source on the local equivalent surface, and screen out the equivalent source with the strongest coupling.

[0139] The back-door coupling simulation module is used to calculate the induced voltage and current at the cable terminal by using the selected equivalent source as excitation and combining the transfer function of the slot structure and the slot-cable path coupling function.

[0140] The front-door coupling simulation module is used to calculate the antenna port coupling voltage and device response by using a selected equivalent source as excitation and combining the frequency-selective surface transfer function and the antenna coupling transfer function.

[0141] The damage assessment module is used to comprehensively assess the damage effect based on preset damage criteria.

[0142] The present invention provides a multi-source, multi-scale, high-power microwave UAV damage simulation method that can realistically reflect the impact of flight status.

[0143] The simulation method of this invention dynamically constructs a multi-source synthetic equivalent source based on the real-time flight status of the UAV, accurately characterizing the impact of changes in the relative positional relationship between the UAV and various interference sources on the damage effect during combat.

[0144] In real combat scenarios, UAVs are in constant motion, and their position and attitude (pitch, roll, yaw) change dynamically over time. This change causes the geometric relationship (distance, incident angle, polarization direction) between the UAV and various interference sources to constantly change, which in turn affects the superposition field distribution of multi-source microwaves on the surface of the UAV and the response of each coupling channel. Traditional simulation methods usually assume that the target is stationary or only consider a few typical attitudes, and cannot accurately characterize the dynamic electromagnetic environment in actual combat.

[0145] This invention achieves a refined description of the UAV's motion process by introducing parametric modeling of flight state. Specifically, the UAV's flight trajectory is decomposed into a series of discrete moments or typical states, each state represented by a position vector. and attitude angle ( , , The location is uniquely determined. Simultaneously, the fixed positions of each interference source are also determined. The radiation parameters (power, radiation pattern, initial phase, frequency) are used as known inputs.

[0146] By analyzing the probability of damage under different flight attitudes, we can guide UAVs to avoid high-threat attitudes or optimize the deployment of interference sources to cover the entire flight profile of the target.

[0147] This invention provides a multi-source, cross-scale, high-power microwave UAV damage simulation method that achieves efficient cross-scale coupling.

[0148] The destructive effects of high-power microwaves on UAVs involve multi-scale physical processes, from the external radiation field (meter-scale) to the internal sensitive devices (micrometer-scale). Traditional methods, if employing full-scale detailed modeling, suffer from enormous computational costs and difficulty in convergence; while stepwise simplification may result in the loss of physical information about critical coupling paths. By pre-constructing transfer functions and using aperture-cable coupling functions, a rapid mapping of the complete transfer link from the external field to the internal cable and then to the device is achieved, avoiding repetitive multi-scale simulations.

[0149] Once the transfer function and coupling function are constructed, they can handle any combination of incident directions and frequencies, significantly reducing the number of repeated simulations and improving computational efficiency.

[0150] The transfer function preserves the resonance characteristics and phase information of the key structure, while the coupling function takes into account the geometric sensitivity of the cable layout, ensuring the accuracy of cross-scale mapping.

[0151] New structures (such as different hole shapes and cable types) can be seamlessly integrated into existing frameworks simply by building their transfer functions separately, demonstrating strong modularity and scalability.

[0152] This invention provides a multi-source, multi-scale, high-power microwave UAV damage simulation method, which is easy to verify experimentally.

[0153] The reliability of simulation methods depends on comparison and verification with experiments. Traditional whole-system irradiation tests are costly, time-consuming, and have poor repeatability, and are difficult to fully verify during the simulation stage. The local equivalent surface strategy allows simulations to be directly compared with local irradiation tests, achieving closed-loop verification between simulation and experiment, and improving the reliability of simulations.

[0154] A test specimen with dimensions comparable to the local equivalent surface was fabricated, including real key structures (such as a metal plate with perforations and an antenna module with a radome). In a microwave anechoic chamber, the test specimen was illuminated using a standard gain horn antenna, with parameters such as incident angle, frequency, and polarization consistent with the simulation settings. Key responses (such as the field strength behind the perforations, cable termination voltage, and antenna port voltage) were measured.

[0155] The results of local irradiation experiments are compared with the simulation results based on local equivalent sources to verify the accuracy of the equivalent source construction and the effectiveness of the transfer function / coupling function. If the deviation is within an acceptable range (e.g., 3 dB), the simulation method is considered reliable. If the deviation is large, the model parameters or the equivalent source construction method are corrected until they match.

[0156] Improving the accuracy of local models can, in turn, improve the overall simulation accuracy.

[0157] This invention provides a multi-source, multi-scale, high-power microwave UAV damage simulation method that significantly reduces the computational load of multi-source collaborative simulation.

[0158] In multi-source collaborative scenarios, if each source is modeled separately and their interactions are considered, the computational cost increases linearly with the number of sources N, and the complex structure of the entire system must be handled, leading to a significant increase in simulation time. By using the local equivalent surface strategy, the multi-source collaborative problem is simplified into a single-source problem. The computational cost is reduced inversely proportional to the cube of the number of sources, and only the local model needs to be processed instead of the entire system model, which greatly improves simulation efficiency.

[0159] By constructing local equivalent surfaces near critical areas, the size of these surfaces only needs to encompass the critical structure (such as holes and slots), and is much smaller than the overall machine size. Through dimensionality reduction using local equivalent surfaces, the number of mesh elements can be reduced by one to two orders of magnitude.

[0160] Multi-source field synthesis is performed beforehand by pre-calculating the vector superposition of all source contributions on a local equivalent surface to obtain the total field distribution. This calculation only involves analytical formulas or fast field superposition and does not require full-wave simulation.

[0161] Based on the total field distribution, an equivalent source is constructed, and subsequent simulations such as aperture coupling and cable coupling are all excited by this equivalent source. In this way, regardless of the actual number of microwave sources, only one coupling simulation is required.

[0162] Example 1

[0163] This embodiment takes a certain type of fixed-wing UAV as the target object and describes in detail the specific implementation process of the method of the present invention, as follows:

[0164] Step 1: Pre-construction of critical system transfer function library:

[0165] The heat dissipation hole on the left side of the UAV fuselage was selected as a typical slot structure, with a slot size of 15cm × 5cm. A local model of the slot was established in electromagnetic simulation software, and frequency sweep simulation was performed using the finite-difference time-domain (FDTD) method. The frequency sweep range was 1~18 GHz, with a step of 0.1 GHz; the incident angle sweep range was θ=0°~90°, φ=0°~360°, with a step of 5°; the electric field response 1cm behind the slot was recorded to obtain the transfer function of the slot structure. Simulations revealed that resonance enhancement occurred near 14 GHz, with the transfer function amplitude reaching -8 dB.

[0166] The GPS navigation antenna (1.575 GHz) and data link antenna (2.4 GHz) were selected as the front-door coupling research objects; a single-element antenna model was established, and the antenna coupling transfer function was obtained through frequency sweep simulation. .

[0167] Step 2: Constructing equivalent sources for multi-source synthesis considering flight conditions:

[0168] Two high-power microwave sources, each with a frequency of 12 GHz and a radiated power of 1 GW, are set up at azimuth angles of 45° and -45° respectively, with pitch angles of 30° each, at a distance of 500m from the target. Typical flight states for the UAV are set, including level flight (yaw angles of 0°, 30°, and 60°), climb (pitch angle of 15°), and hover (roll angle of 30°), etc. Figure 8 As shown.

[0169] Four high-power microwave sources are configured, each with a frequency of 12 GHz and a radiated power of 1 GW. Their deployment locations are at azimuth angles of +45°, -45°, 135°, and 225°, with an elevation angle of 30° for each source. Other conditions remain unchanged. Figure 9 As shown.

[0170] The key part is the heat dissipation vent on the left side of the chassis, in front of the vent opening. A local equivalent surface (a cuboid with dimensions of 20cm × 10cm × 5cm) is constructed. For each flight state, the total field distribution of the source on the local equivalent surface is calculated, and an equivalent source is constructed. The equivalent results for dual sources are as follows: Figure 10 As shown.

[0171] The internal field of the aperture structure was rapidly evaluated under various flight conditions using the transfer function of the aperture structure. It was found that the field strength inside the aperture was the highest when the UAV was in level flight with a yaw angle of 30°, approximately 2.5 times higher than that of a single source. The equivalent source under this condition was then translated outwards along the normal direction to a distance from the aperture. This serves as the input for subsequent detailed simulations.

[0172] Step 3, Simulation of backdoor coupling effect:

[0173] Using the equivalent source constructed in step 2 as the excitation, the field distribution inside the pore is calculated using the finite-difference time-domain (FDTD) method, and the results are obtained. ,like Figure 11 As shown.

[0174] Constructing the orifice-cable path coupling function: In the simulation software, a local model containing the orifice and a section of cable is established. The cable length is L=50 cm and terminated with a 50Ω load. Parameter sweep: distance d from 1 cm to 30 cm, in 1 cm increments; included angle α from 0° to 90°, in 5° increments; frequency f from 13 GHz to 15 GHz, in 0.1 GHz increments; The simulation results are then compiled into a three-dimensional interpolation table.

[0175] Based on the internal cable layout of the UAV, the critical cable for the flight control system is located 15cm behind the slot, at a 45° angle to the normal of the slot, and is 50cm long. Querying the coupling function yields an induced voltage of 45V at the cable termination, exceeding the typical damage threshold of 20V for CMOS circuit I / O ports.

[0176] Step 4: Simulation of front door coupling effect:

[0177] Using the equivalent source constructed in step 2 as excitation, the antenna port coupling voltage was calculated to be 15V using the GPS antenna pattern. The coupling voltage was then fed into the receiver front-end circuit simulation. With the limiter threshold set to 10 dBm, 5V still entered the low-noise amplifier after limiting. Based on the electrothermal analogy model, the calculated channel temperature rise of the low-noise amplifier exceeded 200℃, approaching the damage threshold.

[0178] Step 5: Comprehensive assessment of damage effects:

[0179] Rear door coupling path: The induced voltage at the cable terminal is 45V>20V, indicating that the flight control computer interface may be permanently damaged.

[0180] Front door coupling path: The calculated power density at the antenna is -20dBm / m², which exceeds the 30dB suppression threshold compared to the typical receiving density of -90dBm, indicating that the GPS receiver is experiencing functional suppression; at the same time, the temperature rise of the low-noise amplifier is close to the damage threshold, posing a risk of damage.

[0181] Evaluation of the shielding effectiveness of the structural design.

[0182] The overall damage probability was calculated using the Monte Carlo method, taking into account factors such as transmission power fluctuation ±1dB, attitude change ±10°, and threshold distribution. After 10,000 simulations, the damage probability was approximately 65%.

[0183] Example 2

[0184] This embodiment is basically the same as embodiment 1, except that in step 4, antenna pattern data is used for front door coupling analysis.

[0185] The data used is a table of the transmission coefficients of a certain type of radome at a frequency of 2.4 GHz, which includes transmission amplitude data (without phase information) for incident angles θ=0°~60° and φ=0°~360°; the phase information is estimated by using the minimum phase approximation, assuming that the transmission phase changes linearly with the incident angle.

[0186] Import the directional pattern data into the system to replace... Multi-source transmission field synthesis is performed to calculate the synthesized field after passing through the radome, and then the data link antenna port coupling voltage is obtained; the subsequent steps are the same as in Example 1.

[0187] Example 3

[0188] This embodiment provides a specific implementation method for local irradiation test verification:

[0189] A partial structural component, including heat dissipation holes and a section of cable, was fabricated with dimensions equivalent to the local equivalent surface in step 2 of Example 1 (20cm × 10cm × 5cm). In a microwave anechoic chamber, the partial structural component was illuminated with a standard gain horn antenna at an incident angle of 45° and a frequency of 14GHz, and the voltage at the cable termination was measured.

[0190] Meanwhile, a local model identical to the test piece was constructed in the simulation software, and the cable terminal voltage was calculated by exciting it with a plane wave (corresponding to the horn antenna illumination in the test).

[0191] The simulation results are compared with the experimental results to verify the accuracy of the equivalent source construction and coupling function calculation; if the deviation is within 3dB, the method is considered effective.

[0192] Example 4

[0193] This embodiment presents a multi-source, multi-scale, high-power microwave UAV damage simulation method, the specific method of which is as follows:

[0194] S1. Construct a key system transfer function library: Construct angle-related transfer functions for key UAV systems using the frequency sweeping method;

[0195] Specifically, the angle-related transfer function for constructing the key UAV system includes the following aspects:

[0196] S1.1 Constructing the transfer function of the aperture structure, the specific process is as follows: Establish a local geometric model containing typical apertures, use plane wave excitation, perform frequency sweep simulation within a preset frequency range and incident angle range, record the electric field response of key points inside the aperture, and obtain the transfer function of the aperture structure. :

[0197] (1);

[0198] in, For frequency, The incident elevation angle, It is the azimuth angle. For a strong incident field, The electric field response at key points inside the aperture;

[0199] The transfer function of the slot structure was constructed using the finite-difference time-domain method or the finite element method, with a frequency scanning range of 0.1~30 GHz and an elevation and azimuth scanning range of 0°~90°.

[0200] S1.2 Constructing the frequency-selective surface transfer function, the specific process is as follows: Establish a frequency-selective surface, obtain the variation of the transmission coefficient with frequency and incident angle through frequency sweep simulation, and obtain the frequency-selective surface transfer function. ;

[0201] S1.3 Constructing the antenna coupling transfer function, the specific process is as follows: Establish a single-element antenna model, obtain the relationship between the open-circuit voltage at the antenna port and the incident field strength through frequency sweep simulation, and obtain the antenna coupling transfer function. ;

[0202] S1.4 Identify the maximum coupling condition: Identify the frequency-angle combination with the strongest coupling from the transfer function data of the frequency-selective surface and antenna coupling, and use it as the focus of subsequent fine simulation.

[0203] S2. Construct multi-source synthetic equivalent sources based on the UAV flight status, and select the equivalent sources with the strongest coupling;

[0204] S3. Using the selected equivalent source as excitation, construct the slot-cable path coupling function, combine the angle-related transfer function and the slot-cable path coupling function to calculate the induced voltage and current at the cable terminal, and perform backdoor coupling effect simulation.

[0205] S4. Using the selected equivalent source as excitation and combining the angle-related transfer function, calculate the antenna port coupling voltage and perform front door coupling effect simulation.

[0206] S5. Based on the preset damage criteria, perform a comprehensive damage assessment on the responses obtained from S3 and S4.

[0207] Example 5

[0208] The multi-source, multi-scale, high-power microwave UAV damage simulation method in this embodiment, based on embodiment 4, has the following specific method in S2:

[0209] S2.1, Set the drone's time. The position is Attitude angles include pitch, roll, and yaw; the positions of each high-power microwave source are set. And parameters such as peak power, operating frequency, pulse width, pulse repetition frequency, and transmit antenna gain;

[0210] S2.2 Identify key components: Based on the UAV structural analysis, identify at least one key component affected by the coupling effect, including heat dissipation holes, structural mounting seams, sensor windows, and antenna covers.

[0211] S2.3 Constructing a Local Equivalent Surface: Select a closed surface surrounding the critical part as a local equivalent surface. The local equivalent surface is a cuboid or cylindrical surface, with a size 3 to 5 times the feature size of the critical part, and a distance from the critical part. ,in The wavelength corresponding to the operating frequency of high-power microwaves;

[0212] S2.4 Calculate the total field distribution on the local equivalent surface: For each flight state, based on the geometric relationship between each high-power microwave source and the UAV, calculate the incident field strength of each equivalent source on the local equivalent surface, and perform vector superposition, as shown in the following formula:

[0213] (2);

[0214] in Overall strength It is the field strength of each high-power microwave source at the local equivalent surface, which is determined by the source's peak power, radiation pattern, propagation distance, transmitting antenna gain, and atmospheric attenuation; It includes the source initial phase and the propagation phase, and it is related to its operating frequency and initial phase; For the drone's position parameters, Position parameters for each high-power microwave source; i It is the serial number of the high-power microwave source. j It is the mathematical symbol for imaginary numbers;

[0215] S2.5 Constructing Local Equivalent Sources: Based on the total field distribution on the local equivalent surface, the equivalent current or magnetic current source distribution on the equivalent surface is obtained by solving the inverse problem, forming a multi-source synthetic equivalent source corresponding to the flight state;

[0216] S2.6 Filtering the maximum coupling equivalent source: For the multi-source synthetic equivalent sources corresponding to each flight state constructed in S2.5, the angle correlation transfer function constructed in S1 is used to quickly evaluate the internal coupling response of each equivalent source in the key parts, and filter out the equivalent source with the maximum coupling.

[0217] S2.7 Constructing an equivalent source (radiation source) at 1 / 4 wavelength: Shift the selected equivalent source outward along the normal direction to a distance from the critical location. At this point, a new equivalent source (radiation source) is constructed as the input for subsequent detailed simulation.

[0218] Example 6

[0219] The multi-source, multi-scale, high-power microwave UAV damage simulation method in this embodiment, based on embodiment 4, has the following specific method in S3:

[0220] S3.1 Calculation of the field inside the aperture: Using the equivalent source (radiation source) constructed in S2.7 as the excitation, and utilizing the transfer function of the aperture structure constructed in S1. Calculate the internal field of the pore. Alternatively, a local full-wave simulation can be performed directly, such as a full-wave simulation model of a microwave equivalent source irradiating a hole in a cylinder, as shown in the example. Figure 5 As shown;

[0221] S3.2 Constructing the Through-Cable Path Coupling Function: Constructing the coupling function from the field behind the through-hole to the cable port through parametric simulation. ,in The distance from the center of the slot to the cable. The angle between the cable and the normal to the hole / slot. For cable length, For frequency;

[0222] The pinhole-cable path coupling function is constructed using a parametric scanning method, with a distance... The scanning range is 1mm to 50cm, with an included angle. The scanning range is 0°~90°, and the cable length is... The scanning range is 10cm to 1m. The simulation results are organized into a multidimensional interpolation table or fitted using a neural network.

[0223] S3.3 Cable Termination Response Calculation: Based on the actual cable layout parameters inside the UAV, the induced voltage or current at the cable termination is obtained by querying the slot-cable path coupling function constructed in S3.2.

[0224] (3);

[0225] in For the field inside the pore;

[0226] S3.4 Sensitive Device Injection Analysis: The induced voltage or current at the cable terminal obtained in S3.3 is used as an excitation source and injected into the subsequent circuit simulation model to analyze the response state of the sensitive device.

[0227] Example 7

[0228] The multi-source, multi-scale, high-power microwave UAV damage simulation method in this embodiment, based on embodiment 4, has the following specific method in S4:

[0229] S4.1 Calculation of Frequency-Selective Surface Transmission Field: Using the equivalent source constructed in S2.7 as the excitation, and utilizing the frequency-selective surface transfer function constructed in S1... Calculate the synthesized field after passing through the frequency-selective surface. :

[0230] (4);

[0231] in For the first i The incident field vector of a high-power microwave source on the outer surface of a frequency-selective surface. For the first i A high-power microwave source polarization unit vector; i For drones and the i The azimuth angle formed between the high-power microwave sources, i.e. the rotation angle about the reference axis; For drones and the i The pitch angle formed between the high-power microwave sources, that is, the angle between the incident direction and the reference axis; The mathematical symbol for imaginary numbers;

[0232] The frequency-selective surface transfer function It can be replaced by existing radomes or frequency-selective surface pattern data;

[0233] If the surface pattern data only contains amplitude information, the phase information can be inverted using minimum phase approximation algorithms such as Hilbert transform and homomorphic filtering, or the phase information can be reconstructed based on a physical model that shows the transmission phase changing linearly with the angle.

[0234] S4.2 Antenna Port Coupling Voltage Calculation: Using the radome transmission field obtained in S4.1 as the excitation, the antenna coupling transfer function constructed in S1 is used... Calculate the antenna port coupling voltage :

[0235] (5);

[0236] The local elevation angle of the incident wave relative to the antenna mounting plane; The local azimuth angle of the incident wave relative to the antenna mounting plane;

[0237] S4.3 Simulation of the receiver front-end circuit: The antenna port coupling voltage obtained in S4.2 is... As an excitation source, the simulation model of the receiving front-end circuit is imported to analyze the response of the limiter and the low-noise amplifier.

[0238] S4.4 Device Temperature Rise and Damage Analysis: Based on the electrothermal analogy method, a thermal equivalent model of the low-noise amplifier is established. According to the pulse width and repetition frequency parameters, the device temperature rise is calculated to determine whether it exceeds the damage threshold.

[0239] Example 8

[0240] The damage simulation method for multi-source, multi-scale, high-power microwave UAVs in this embodiment, based on Embodiment 4, includes the following damage criteria:

[0241] Back door coupling damage judgment: The induced voltage and current of the cable terminal obtained in S3.4 are compared with the preset device voltage and current damage thresholds. If the thresholds are exceeded, permanent damage is determined to have occurred.

[0242] Front door coupling suppression judgment: The antenna port coupling voltage obtained from S4.2 is converted into the power density at the antenna and compared with the preset receiving density suppression threshold. If it exceeds the threshold, it is determined that functional suppression has occurred.

[0243] Front door coupling damage judgment: Compare the device temperature rise obtained in S4.4 with the preset device thermal damage threshold. If it exceeds the threshold, it is determined that permanent damage has occurred.

[0244] Multi-source synergistic cumulative effect assessment: For scenarios involving simultaneous or time-sharing irradiation from multiple sources, considering power superposition and energy accumulation effects, the cumulative criterion is applied to assess the overall damage probability.

[0245] (6);

[0246] in: It is an aggregated power function that describes the total power of multi-source pulses over time (including simultaneous superposition or time-division accumulation effects). The thermal response function is denoted by , and the thermal time constant of the device is denoted by . t Exponentially decaying kernels, for example ,in It is a constant, which is related to the heat accumulation and heat dissipation characteristics of the device, and can be obtained through actual measurement; The damage energy threshold is expressed in joules.

Claims

1. A multi-source, multi-scale, high-power microwave UAV damage simulation method, characterized in that, The specific method is as follows: S1. Construct the angle-related transfer function of the key UAV system using the frequency sweeping method; S2. Construct multi-source synthetic equivalent sources based on the UAV flight status, and select the equivalent sources with the strongest coupling; S3. Using the selected equivalent source as excitation, construct the slot-cable path coupling function, combine the angle-related transfer function and the slot-cable path coupling function to calculate the induced voltage and current at the cable terminal, and perform backdoor coupling effect simulation. S4. Using the selected equivalent source as excitation and combining the angle-related transfer function, calculate the antenna port coupling voltage and perform front door coupling effect simulation. S5. Based on the preset damage criteria, conduct a comprehensive damage assessment on the obtained response; The angle-related transfer function for constructing the key UAV system in S1 specifically includes the following aspects: S1.1 Constructing the transfer function of the aperture structure, the specific process is as follows: Establish a local geometric model containing typical apertures, use plane wave excitation, perform frequency sweep simulation within a preset frequency range and incident angle range, record the electric field response of key points inside the aperture, and obtain the transfer function of the aperture structure. : （1）； in, For frequency, The incident elevation angle, It is the azimuth angle. For a strong incident field, The electric field response at key points inside the aperture; The transfer function of the pore structure is constructed using the finite-difference time-domain method or the finite element method, with a frequency scanning range of 0.1~30 GHz and an elevation and azimuth scanning range of 0°~90°. S1.2 Constructing the frequency-selective surface transfer function, the specific process is as follows: Establish a frequency-selective surface, obtain the variation of the transmission coefficient with frequency and incident angle through frequency sweep simulation, and obtain the frequency-selective surface transfer function. ; S1.3 Constructing the antenna coupling transfer function, the specific process is as follows: Establish a single-element antenna model, obtain the relationship between the open-circuit voltage at the antenna port and the incident field strength through frequency sweep simulation, and obtain the antenna coupling transfer function. ; S1.4 Identify the maximum coupling condition: Identify the frequency-angle combination with the strongest coupling from the transfer function data of the frequency selective surface and antenna coupling, and use it as the focus of subsequent fine simulation; Constructing the aperture-cable path coupling function: A coupling function from the field behind the aperture to the cable port is constructed through parametric simulation. ,in The distance from the center of the slot to the cable. The angle between the cable and the normal to the hole / slot. For cable length, For frequency; The pinhole-cable path coupling function is constructed using a parametric scanning method, with a distance... The scanning range is 1 mm to 50 cm, with an included angle. The scanning range is 0°~90°, and the cable length is... The scanning range is 10cm to 1m. The simulation results are organized into a multidimensional interpolation table or fitted using a neural network.

2. The multi-source, multi-scale, high-power microwave UAV damage simulation method according to claim 1, characterized in that, The specific method of S2 is as follows: S2.1, Set the drone's time. The position is Attitude angles include pitch, roll, and yaw; the positions of each high-power microwave source are set. And peak power, operating frequency, pulse width, pulse repetition frequency, and transmit antenna gain; S2.2 Identify key components: Based on the UAV structural analysis, identify at least one key component affected by the coupling effect, including heat dissipation holes, structural mounting seams, sensor windows, and antenna covers. S2.3 Constructing a Local Equivalent Surface: Select a closed surface surrounding the critical part as a local equivalent surface. This local equivalent surface is a cuboid or cylindrical surface, with dimensions 3-5 times the characteristic dimension of the critical part, and a distance from the critical part. ,in The wavelength corresponding to the operating frequency of high-power microwaves; S2.4 Calculate the total field distribution on the local equivalent surface: For each flight state, based on the geometric relationship between each high-power microwave source and the UAV, calculate the incident field strength of each equivalent source on the local equivalent surface, and perform vector superposition, as shown in the following formula: （2）； in Overall strength It is the field strength of each high-power microwave source at the local equivalent surface, which is determined by the source's peak power, radiation pattern, propagation distance, transmitting antenna gain, and atmospheric attenuation; It includes the source initial phase and the propagation phase, and it is related to its operating frequency and initial phase; For the drone's position parameters, Position parameters for each high-power microwave source; i It is the serial number of the high-power microwave source. j It is the mathematical symbol for imaginary numbers; S2.5 Constructing Local Equivalent Sources: Based on the total field distribution on the local equivalent surface, the equivalent current or magnetic current source distribution on the equivalent surface is obtained by solving the inverse problem, forming a multi-source synthetic equivalent source corresponding to the flight state; S2.6 Filtering the maximum coupling equivalent source: For the multi-source synthetic equivalent sources corresponding to each flight state constructed in S2.5, the angle correlation transfer function constructed in S1 is used to quickly evaluate the internal coupling response of each equivalent source in the key parts, and filter out the equivalent source with the maximum coupling. S2.7 Constructing an equivalent source at 1 / 4 wavelength: Shift the selected equivalent source outwards along the normal direction to a distance from the critical location. At this point, a new equivalent source is constructed as the input for subsequent detailed simulation.

3. The multi-source, multi-scale, high-power microwave UAV damage simulation method according to claim 2, characterized in that, The specific method of S3 is as follows: S3.1 Calculation of the internal field of the pore: Using the equivalent source constructed in S2.7 as the excitation, and utilizing the pore structure transfer function constructed in S1. Calculate the internal field of the pore. Alternatively, a local full-wave simulation can be performed directly; S3.2 Construct the hole-cable path coupling function; S3.3 Cable Termination Response Calculation: Based on the actual cable layout parameters inside the UAV, the induced voltage or current at the cable termination is obtained by querying the slot-cable path coupling function constructed in S3.

2. （3）； in For the field inside the pore; S3.4 Sensitive Device Injection Analysis: The induced voltage or current at the cable terminal obtained in S3.3 is used as an excitation source and injected into the subsequent circuit simulation model to analyze the response state of the sensitive device.

4. The multi-source, multi-scale, high-power microwave UAV damage simulation method according to claim 3, characterized in that, The specific method of S4 is as follows: S4.1 Calculation of Frequency-Selective Surface Transmission Field: Using the equivalent source constructed in S2.7 as the excitation, and utilizing the frequency-selective surface transfer function constructed in S1... Calculate the synthesized field after passing through the frequency-selective surface. : （4）； in For the first i The incident field vector of a high-power microwave source on the outer surface of a frequency-selective surface. For the first i A high-power microwave source polarization unit vector; i For drones and the i The azimuth angle formed between the high-power microwave sources, i.e. the rotation angle about the reference axis; For drones and the i The pitch angle formed between the high-power microwave sources, that is, the angle between the incident direction and the reference axis; The mathematical symbol for imaginary numbers; S4.2 Antenna Port Coupling Voltage Calculation: Using the transmitted field inside the radome obtained in S4.1 as the excitation, the antenna coupling transfer function constructed in S1 is used... Calculate the antenna port coupling voltage : （5）； The local elevation angle of the incident wave relative to the antenna mounting plane; The local azimuth angle of the incident wave relative to the antenna mounting plane; S4.3 Simulation of the receiver front-end circuit: The antenna port coupling voltage obtained in S4.2 is... As an excitation source, the simulation model of the receiving front-end circuit is imported to analyze the response of the limiter and the low-noise amplifier. S4.4 Device Temperature Rise and Damage Analysis: Based on the electrothermal analogy method, a thermal equivalent model of the low-noise amplifier is established. According to the pulse width and repetition frequency parameters, the device temperature rise is calculated to determine whether it exceeds the damage threshold.

5. The multi-source, multi-scale, high-power microwave UAV damage simulation method according to claim 4, characterized in that, The frequency-selective surface transfer function in S4.1 It can be replaced by existing radomes or frequency-selective surface pattern data; If the surface pattern data only contains amplitude information, the phase information can be inverted using minimum phase approximation algorithms such as Hilbert transform and homomorphic filtering, or the phase information can be reconstructed based on a physical model that shows the transmission phase changing linearly with the angle.

6. The multi-source, multi-scale, high-power microwave UAV damage simulation method according to claim 5, characterized in that, The damage criteria in S5 include the following aspects: Back door coupling damage judgment: The induced voltage and current of the cable terminal obtained in S3.4 are compared with the preset device voltage and current damage thresholds. If the thresholds are exceeded, permanent damage is determined to have occurred. Front door coupling suppression judgment: The antenna port coupling voltage obtained from S4.2 is converted into the power density at the antenna and compared with the preset receiving density suppression threshold. If it exceeds the threshold, it is determined that functional suppression has occurred. Front door coupling damage judgment: Compare the device temperature rise obtained in S4.4 with the preset device thermal damage threshold. If it exceeds the threshold, it is determined that permanent damage has occurred. Multi-source synergistic cumulative effect assessment: For scenarios involving simultaneous or time-sharing irradiation from multiple sources, considering power superposition and energy accumulation effects, the cumulative criterion is applied to assess the overall damage probability. （6）； in: This is a power superposition function, describing the time-varying behavior of multi-source pulses. t Total power; Let be the thermal response function, and be the thermal time constant of the device. t Exponentially decaying kernel; The damage energy threshold is expressed in joules.

7. A multi-source, multi-scale, high-power microwave UAV damage simulation system for use in the simulation method described in any one of claims 1 to 6, characterized in that, It includes a transfer function library construction module, a multi-source equivalent source construction module, a backdoor coupling simulation module, a frontdoor coupling simulation module, and a damage assessment module.

8. The multi-source, multi-scale, high-power microwave UAV damage simulation system according to claim 7, characterized in that, The transfer function library construction module is used to construct the angle-related transfer functions of key UAV systems using the frequency sweeping method. The multi-source equivalent source construction module is used to consider the flight state of the UAV, construct a local equivalent surface near the key parts, construct a multi-source synthetic equivalent source based on the total field distribution of each microwave source on the local equivalent surface, and select the equivalent source with the largest coupling. The backdoor coupling simulation module is used to calculate the induced voltage and current at the cable terminal by using the selected equivalent source as excitation and combining the transfer function of the slot structure and the slot-cable path coupling function. The front-door coupling simulation module is used to calculate the antenna port coupling voltage and device response by using the selected equivalent source as excitation and combining the frequency-selective surface transfer function and the antenna coupling transfer function. The damage assessment module is used to comprehensively assess the damage effect based on preset damage criteria.

Images