Arrayed high-energy laser counter-uav test method and apparatus

An arrayed high-energy laser counter-UAV (Unmanned Aerial Vehicle) test method, comprising the following steps:

• Step 1:
• preparing a test apparatus and a UAV, wherein the test apparatus comprises a high-energy laser array assembly, an automatic pan-tilt gimbal, a high-speed camera, an infrared camera, and a laser sight;
• the high-energy laser array assembly comprises a carrier box and a plurality of high-power laser devices, the plurality of high-power laser devices are arranged in the carrier box in parallel to form a laser array, the laser array is configured to emit a high-power laser beam to the UAV to ablate the UAV, the automatic pan-tilt gimbal is configured to drive the carrier box to rotate to adjust an emitting direction of the high-power laser beam, the high-speed camera is configured to take an image of a test process, the infrared camera is configured to take a heat distribution map of the UAV when ablated, and the laser sight is configured to emit aiming laser parallel to the high-power laser beam for detecting the UAV;
• Step 2:
• conducting a laser counter-UAV test;
• turning on the automatic pan-tilt gimbal, the high-speed camera, the infrared camera and the laser sight, remotely controlling the UAV to fly into a test region, and when the laser sight detects the UAV, remotely controlling the laser array to emit the high-power laser beam to the UAV to ablate the UAV;
• Step 3:
• observing an ablation process and analyzing an ablated morphology of the UAV to obtain a radius and ablation depth of a tested ablation area of the UAV;
• Step 4:
• performing finite element simulation on the test process to obtain a radius and ablation depth of a simulated ablation area of the UAV;
• and
• Step 5:
• comparing the radii and ablation depths of the simulated and tested ablation areas of the UAV, adjusting finite element simulation parameters to make a simulation result coincide with a test result, and deriving an ablation depth-time variation curve of a key point in the simulated ablation area;
• wherein Step 4 specifically comprises:
• establishing a two-dimensional geometric model of a UAV housing in COMSOL software and meshing the two-dimensional geometric model;
• setting the finite element simulation parameters and inputting the parameters into a Gaussian beam algorithm in a laser spectral resonator to perform finite element simulation analysis on the ablation process;
• and
• outputting a simulation result to obtain the radius and ablation depth of the simulated ablation area of the UAV, and deriving an ablation depth-time variation curve of a key point in the simulated ablation area;
• the Gaussian beam algorithm in the laser spectral resonator specifically comprises:
• enabling q-parameters of a Gaussian beam to pass through an optical system with a transformation matrix of
• M
• =
• [
• A
• B
• C
• D
• ]
• ,
• wherein a transformation of the q-parameter also follows ABCD law, that is:
• q
• 2
• =
• Aq
• 1
• +
• B
• Cq
• 1
• +
• D
• ;
• q
• 2
• -
• 1
• =
• C
• +
• Dq
• 1
• -
• 1
• A
• +
• Bq
• 1
• -
• 1
• ;
• in the formula:
• q1- and
• q
• 1
• -
• 1
• are q-parameters before transmission transformation, q2 and
• q
• 2
• -
• 1
• are q-parameters after transmission transformation;
• if the Gaussian beam passes through an optical system with a plurality of optical elements, the transformation matrix M is determined by a product of various optical element transformation matrices M1, M2, . . . , Mn;
• when q1- or
• q
• 1
• -
• 1
• -parameter and M1, M2, . . . , Mn are known, solving q- or q−1-parameter at any position z through the ABCD law, and then performing an operation of separating real and imaginary parts based on the q- or q−1-parameter to obtain an isophase surface curvature radius R(z), a spot radius w(z), a position from a beam waist, and a beam waist radius at the position;
• setting that the Gaussian beam in the laser spectral resonator is represented by the q-parameters, that a round-trip transmission matrix at a cavity mirror 1 is
• [
• a
• b
• c
• d
• ]
• ,
• and that the q- and q−1-parameter of a resonant beam satisfy a self-reproduction condition:
• q
• -
• 1
• =
• c
• +
• dq
• -
• 1
• a
• +
• bq
• -
• 1
• ;
• q
• =
• aq
• +
• b
• cq
• +
• d
• ;
• solving for the q−1-parameter to obtain:
• q
• -
• 1
• =
• d
• -
• a
• 2
• ⁢
• b
• ±
• i
• ⁢
• 1
• -
• (
• a
• +
• d
• 2
• )
• 2
• b
• ,
• b
• ≠
• 0
• ;
• q
• -
• 1
• =
• c
• a
• -
• d
• ,
• b
• =
• 0
• ;
• solving for the q-parameter to obtain:
• q
• =
• a
• -
• d
• 2
• ⁢
• c
• ±
• i
• ⁢
• 1
• -
• (
• a
• +
• d
• 2
• )
• 2
• c
• ,
• c
• ≠
• 0
• ;
• q
• =
• b
• d
• -
• a
• ,
• c
• =
• 0
• ;
• wherein i is a unit matrix;
• a plurality of solutions are present when solving for the q-parameter numerically, wherein a rule of choosing the solution is that a beam half-width square of the beam is positive;
• separating the real and imaginary parts after solving for the q-parameter, wherein the real part and the imaginary part correspond to a beam waist position and a beam waist radius of the beam, respectively;
• and solving for beam parameters of the Gaussian beam transformed by the optical system by using the ABCD law to obtain a beam waist position and a beam waist radius of the beam after transmission transformation;
• and
• separating the real and imaginary parts after solving for the q−1-parameter, wherein the real part and the imaginary part correspond to an isophase surface curvature radius and a spot radius of the beam, respectively;
• solving for beam parameters of the Gaussian beam transformed by the optical system by using the ABCD law to obtain an isophase surface curvature radius and a spot radius of the beam after transmission transformation.