CN106469850B - A kind of Thickness Design Method of antenna house - Google Patents
A kind of Thickness Design Method of antenna house Download PDFInfo
- Publication number
- CN106469850B CN106469850B CN201610824698.8A CN201610824698A CN106469850B CN 106469850 B CN106469850 B CN 106469850B CN 201610824698 A CN201610824698 A CN 201610824698A CN 106469850 B CN106469850 B CN 106469850B
- Authority
- CN
- China
- Prior art keywords
- antenna
- thickness
- antenna house
- value
- house
- 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
Classifications
-
- H—ELECTRICITY
- H01—ELECTRIC ELEMENTS
- H01Q—ANTENNAS, i.e. RADIO AERIALS
- H01Q1/00—Details of, or arrangements associated with, antennas
- H01Q1/42—Housings not intimately mechanically associated with radiating elements, e.g. radome
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/30—Circuit design
- G06F30/36—Circuit design at the analogue level
- G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
Abstract
The present invention relates to a kind of Thickness Design Method of antenna house, include the following steps: that the thickness of antenna house is carried out discretization according to the short transverse of cover by (1);(2) it determines the value range of thickness value, and initial value is assigned to the thickness at discrete point;(3) according to the thickness value at discrete point, its statistical value is calculated, (4) calculate the far field of antenna, and therefrom extract gain G1With main beam position B1These electrical performance indexes;(5) transmission coefficient of cover is calculated with transmission line theory, (6) calculate the far field F ' (θ, φ) with cover antenna system, draw far-field pattern, (7) establish mathematical optimization models;(8) whether the electrical performance indexes and thickness distribution for judging the antenna house obtained after optimization meet preset requirement.It improves the electrical performance indexes of antenna house.
Description
Technical field
The invention belongs to Radar Antenna System field, the Thickness Design Method of specifically a kind of antenna house can be used for flying
The structure of device antenna house designs.
Background technique
Antenna house is the wave transparent shell for protecting an antenna from natural environment influence, is made of natural or artificial dielectric's material
Covering, or the bright window of electromagnetism of special shape being made of the dielectric housing of truss support.The antenna house of excellent in design,
Other than having the function of protectiveness, conductibility, reliability, concealment and dicoration etc., whole system each section can also be extended
Service life, reduce Life Cost and operating cost, simplify design, reduce maintenance cost, guaranteeing antenna surface and position
Accuracy creates good working environment to antenna operation personnel.But antenna house can also produce the electromagnetic radiation of ideal antenna
It is raw to influence, so that ideal antenna electric performance is decreased.
Aircraft antenna house uses streamlined cover, often in order to meet the requirement of pneumatic property so as to cause antenna house pair
The electrical property influence of internal antenna is excessive, needs to optimize antenna house to improve its electrical property.Traditional equal thickness
Optimisation technique improves limitation to the electrical property of antenna house, and discrete by carrying out antenna house thickness along cover short transverse
Change, carries out Varying-thickness optimization design, the electrical property of streamlined cover can be significantly improved.
FAN XUEPING utilizes heredity in 2011 paper " optimization based on genetic algorithm to two-dimensional antenna cover collimating fault "
The collimating fault of two dimensional phased array antenna cover is optimized in algorithm, by by antenna house wall thickness along short transverse equal length
Different zones are divided into, using the wall thickness value of each region as optimized variable, and are missed with maximum aim at during antenna scanning
Difference is used as optimization aim, optimal wall Thickness Distribution has been solved using genetic algorithm, to significantly reduce taking aim at for antenna house
Quasi- error.The deficiency of this method is: not accounting for constraining Thickness Distribution, causes the wall thickness change of gained antenna house acute
It is strong, it is difficult to which that processing is realized, this important electrical performance indexes of the gain loss of antenna house are not furthermore also accounted for.
Summary of the invention
It is an object of the invention in view of the above shortcomings of the prior art, provide a kind of electrical property that can improve antenna house
Index, and the Thickness Design Method of the antenna house of the manufacture difficulty of antenna house can be reduced.
To achieve the above object, it the technical scheme is that a kind of Thickness Design Method of antenna house, including walks as follows
It is rapid:
(1) thickness of antenna house is subjected to discretization according to the short transverse of cover;
It is that x/y plane establishes a coordinate system O- using the bottom center of antenna house as origin, bottom surface along the height of cover
Xyz, cover height in the z-direction, uniformly choose 10 discrete points in the direction, and coordinate is denoted as z1,z2,…,z10;
(2) according to the requirement of antenna cover structure performance and electrical property, determine the value range of thickness value, and to discrete point at
Thickness assign initial value;
(3) according to the thickness value at discrete point, its statistical value is calculated, the variance including variance, greatest gradient and its slope.
(4) far field of antenna is calculated, and therefrom extracts gain G1With main beam position B1These electrical performance indexes.
(5) according to the structural parameters of antenna house and material parameter, the transmission coefficient of cover is calculated with transmission line theory
And according to known antenna aperture field E (x, y), calculate through the aperture field after antenna house:
(6) according to through the aperture field E ' (x, y) after antenna house, the far field F ' (θ, φ) with cover antenna system is calculated, is drawn
Far-field pattern processed, and gain G is extracted from the far-field pattern2With main beam position B2These electrical performance indexes, and then determine
Gain loss TL caused by antenna house and collimating fault BSE.
(7) using the thickness value at the discrete point in step (1) as design variable, with variance, the maximum of thickness in step (3)
Electrical performance indexes in slope, the variance of slope and step (6) are design object, establish mathematical optimization models, utilize population
Optimization algorithm solves this model, obtains the thickness distribution of antenna house;
(8) whether the electrical performance indexes and thickness distribution for judging the antenna house obtained after optimization meet preset requirement, if
Meet, then antenna cover structure design scheme is qualified, otherwise, modifies the discretization method and optimization algorithm parameter of antenna house, lays equal stress on
Multiple step (1) to step (8), until obtaining electrical performance indexes and thickness distribution meets the design scheme of preset requirement.
The step (2) determines the value range of thickness value, presses according to the requirement of antenna cover structure performance and electrical property
Following process carries out:
(2a) integrally formed aircraft antenna house generallys use half-wavelength wall thickness, the calculation formula of half-wavelength wall thickness are as follows:
Wherein λ is wavelength, εrIt is the relative dielectric constant of radome material, α is incidence angle of the incident electromagnetic wave to cover wall.
(2b) determines n according to the structural performance requirements of cover, and n is bigger, and cover wall is thicker, and structural behaviour is better, while electrically
Can also can be poorer, therefore take the smallest n value for being able to satisfy structural performance requirements.
(2c) is calculated in the scanning process of internal antenna, the minimum of the electromagnetic wave of aerial radiation to the incidence angle of cover wall
Value αminWith maximum value αmax。
(2d) determines the minimum value d of Thickness as the following formulaminWith maximum value dmax:
Step (4) the specific implementation step is:
(4a) indicates the component in the x, y, z direction for the coordinate system established in step (1) with i, j, k respectively, according to known
Antenna aperture field distribution E (x, y), calculate the far field value F (θ, φ) of antenna:
Wherein, θ, φ are spherical coordinates angle of the point of observation in O-xyz,λ is the wavelength of antenna, according to antenna work
Working frequency f and light velocity c, passes through formulaIt is calculated, s is the area of integral unit;
(4b) draws Antenna Far Field directional diagram according to the far field F (θ, φ) of ideal antenna, and gain is extracted from directional diagram
G1With main beam position B1These electrical performance indexes.
Step (5) the specific implementation step is:
(5a) is directed to the coordinate system that step (1) is established, and the component in its x, y, z direction is indicated with i, j, k respectively, it is known that
Antenna aperture field distribution is denoted as E (x, y).
(5b) establishes the geometrical model of antenna house according to the structure type of antenna house in Mercantile Models analysis software, if
Setting side length of element is 0.2 λ, and wherein λ is the wavelength of antenna, carries out grid dividing to model;
(5c) calculates the transmission on covering at each point with transmission line theory according to the structural parameters and material parameter of antenna house
Coefficient
(5c1) according to geometry of antenna house and incident aperture field, find out incident angle α on antenna house at each point and
Normal angle at electromagnetic wave incident line and incidence point is denoted as incident angle α by polarizing angle β, by polarization of electromagnetic wave direction with
The angle of plane of incidence is denoted as polarizing angle β, and wherein plane of incidence is made of electromagnetic wave incident line and the normal at incidence point;
(5c2) is according to the thickness value d at discrete point1,d2,….,d10, obtained using cubic Bézier curves method smooth
The thickness profile of antenna house:
D (z)=fcubic-B-spline(d1,...,d10)
(5c3) is according to antenna house thickness d everywhere, relative dielectric constant εr, losstangenttanδ, calculate each on covering
Horizontal polarized components transmission coefficient at pointWith perpendicular polarisation components transmission coefficient
Wherein, ZH=cos α,These
Parameter is intermediate variable;TH、TVRespectivelyModulus value, ηH、ηVRespectivelyPhase;
(5c4) is according to horizontal polarized components transmission coefficientWith perpendicular polarisation components transmission coefficientObtain main polarization point
The transmission coefficient of amount:
Wherein,For intermediate variable;
(5d) multiplied by the transmission coefficient of its corresponding point position, calculates the aperture field being incident on antenna house through after covering
Aperture field:
Step (6) the specific implementation step is:
(6a) is calculated as follows according to being acquired in step (1) through the aperture field E ' (x, y) after antenna house through antenna
The far field F ' (θ, φ) that aperture field after cover generates:
Wherein, θ, φ are spherical coordinates angle of the point of observation in rectangular coordinate system O-xyz, k0For free-space propagation constant,
By formulaIt calculates, λ is that the wavelength of antenna according to operating frequency of antenna f and light velocity c passes through formulaIt calculates
It arrives, s is the area of integral unit;
(6b) draws the Antenna Far Field side after covering according to the far field F ' (θ, φ) generated through the aperture field after antenna house
Xiang Tu, and gain G is extracted from directional diagram2With main beam position B2These electrical performance indexes;
(6c) calculates gain loss TL caused by antenna house according to the electrical performance indexes calculated in step (4) and step (5)
With collimating fault BSE:
TL=G1- G2
BSE=| B1- B2|
Step (7) the specific implementation step is:
(7a) using the thickness value at the discrete point in step (1) as design variable, with the variance of thickness in step (3), most
Gain loss caused by antenna house and collimating fault are design object in big slope, the variance of slope and step (6), are established as follows
Consider the aircraft antenna house Varying-thickness Optimized model of thickness control:
Wherein BSEmaxIt is the maximum value of the collimating fault caused by antenna house under all operating conditions, TLmaxIt is in all operating conditions
The maximum value of gain loss caused by lower antenna house, vd are the variances of design variable, and mds is the maximum of the slope of design variable
Value, vds is the variance of the slope of design variable, BSE0,TL0,vd0,mds0And vds0It is the normalization provided according to particular problem
Coefficient, dminIt is the value lower limit of design variable, dmaxIt is the value upper limit of design variable;
(7b) is using particle swarm optimization algorithm (PSO) come to the aircraft antenna house Varying-thickness optimization mould for considering thickness control
Type is solved, and optimal antenna house thickness distribution is obtained, and the population scale of particle swarm optimization algorithm is taken as 50, and evolutionary generation takes
It is 200, for inertia weight as evolutionary generation is from 0.9 linear decrease to 0.4, aceleration pulse is taken as 2.
The present invention considers the gain loss of antenna house due to introducing thickness control factor, thus with existing change
Thickness Design Method is compared, and the electrical performance indexes of antenna house have not only been improved, but also reduces the manufacture difficulty of antenna house.
Detailed description of the invention
Fig. 1 is realization general flow chart of the invention;
Fig. 2 is the sub-process figure calculated in the present invention through the aperture field after antenna house;
Fig. 3 is the antenna and antenna house relation schematic diagram that the present invention uses;
Fig. 4 is the structural schematic diagram for certain aircraft antenna house that present invention emulation uses;
Fig. 5 is to optimize the resulting collimating fault comparison diagram in front and back to certain aircraft antenna house with the present invention;
Fig. 6 is to optimize the resulting gain loss comparison diagram in front and back to certain aircraft antenna house with the present invention;
Fig. 7 is to optimize the resulting thickness distribution comparison diagram in front and back to certain aircraft antenna house with the present invention.
Specific embodiment
Present invention is further described in detail referring to the drawings.
Referring to Fig.1, the specific steps of the present invention are as follows:
The thickness of antenna house is carried out discretization according to the short transverse of cover by step (1).
As shown in Fig. 2, being that x/y plane establishes one using the bottom center of antenna house as origin, bottom surface along the height of cover
Coordinate system O-xyz, cover height in the z-direction, uniformly choose 10 discrete points in the direction, and coordinate is denoted as z1,z2,…,z10,
θ in figuresIndicate the scan angle of antenna.
Step (2) determines the value range of Thickness, and assigns initial value to the thickness at discrete point.
(2a) integrally formed aircraft antenna house generallys use half-wavelength wall thickness, the calculation formula of half-wavelength wall thickness are as follows:
Wherein λ is wavelength, εrIt is the relative dielectric constant of radome material, α is incidence angle of the incident electromagnetic wave to cover wall.
(2b) determines n according to the structural performance requirements of cover, and n is bigger, and cover wall is thicker, and structural behaviour is better, while electrically
Can also can be poorer, therefore take the smallest n value for being able to satisfy structural performance requirements.
(2c) is calculated in the scanning process of internal antenna, the minimum of the electromagnetic wave of aerial radiation to the incidence angle of cover wall
Value αminWith maximum value αmax。
(2d) determines the minimum value d of Thickness as the following formulaminWith maximum value dmax:
(2e) is in section [dmin,dmax] in randomly select one group of number as the thickness value at discrete point in step (1).
Step (3) calculates the variance of its variance, greatest gradient and its slope according to the thickness value at discrete point;
According to the thickness value d at discrete point1,d2,….,d10, the variance vd of these numbers is calculated, further according to its coordinate
z1,z2,…,z10, calculate the slope value of adjacent discrete point:
Calculate the maximum value mds and its variance vds of these slope values.
Step (4), calculates the far field of antenna, and therefrom extracts gain G1With main beam position B1These electrical performance indexes;
(4a) indicates the component in the x, y, z direction for the coordinate system established in step (1) with i, j, k respectively, according to known
Antenna aperture field distribution E (x, y), calculate the far field value F (θ, φ) of antenna:
Wherein, θ, φ are spherical coordinates angle of the point of observation in O-xyz,λ is the wavelength of antenna, according to antenna work
Working frequency f and light velocity c, passes through formulaIt is calculated, s is the area of integral unit;
(4b) draws Antenna Far Field directional diagram according to the far field F (θ, φ) of ideal antenna, and gain is extracted from directional diagram
G1With main beam position B1These electrical performance indexes.
Step (5) calculates the transmission coefficient of cover with transmission line theory according to the structural parameters and material parameter of antenna houseAnd according to known antenna aperture field E (x, y), calculate through the aperture field after antenna house:
Referring to Fig. 3, this step is implemented as follows:
(5a) is directed to the coordinate system that step (1) is established, and the component in its x, y, z direction is indicated with i, j, k respectively, it is known that
Antenna aperture field distribution is denoted as E (x, y).
(5b) establishes the geometrical model of antenna house according to the structure type of antenna house in Mercantile Models analysis software, if
Setting side length of element is 0.2 λ, and wherein λ is the wavelength of antenna, carries out grid dividing to model;
(5c) calculates the transmission on covering at each point with transmission line theory according to the structural parameters and material parameter of antenna house
Coefficient
(5c1) according to geometry of antenna house and incident aperture field, find out incident angle α on antenna house at each point and
Normal angle at electromagnetic wave incident line and incidence point is denoted as incident angle α by polarizing angle β, by polarization of electromagnetic wave direction with
The angle of plane of incidence is denoted as polarizing angle β, and wherein plane of incidence is made of electromagnetic wave incident line and the normal at incidence point;
(5c2) is according to the thickness value d at discrete point1,d2,….,d10, obtained using cubic Bézier curves method smooth
The thickness profile of antenna house:
D (z)=fcubic-B-spline(d1,...,d10)
(5c3) is according to antenna house thickness d everywhere, relative dielectric constant εr, losstangenttanδ, calculate each on covering
Horizontal polarized components transmission coefficient at pointWith perpendicular polarisation components transmission coefficient
Wherein, ZH=cos α,These
Parameter is intermediate variable;TH、TVRespectivelyModulus value, ηH、ηVRespectivelyPhase;
(5c4) is according to horizontal polarized components transmission coefficientWith perpendicular polarisation components transmission coefficientObtain main polarization point
The transmission coefficient of amount:
Wherein,For intermediate variable;
(5d) multiplied by the transmission coefficient of its corresponding point position, calculates the aperture field being incident on antenna house through after covering
Aperture field:
Step (6) calculates the far field F ' (θ, φ) with cover antenna, draws according to through the aperture field E ' (x, y) after antenna house
Far-field pattern processed, and gain G is extracted from the far-field pattern2With main beam position B2These electrical performance indexes, and then determine
Gain loss TL caused by antenna house and collimating fault BSE.
(6a) is calculated as follows according to being acquired in step (1) through the aperture field E ' (x, y) after antenna house through antenna
The far field F ' (θ, φ) that aperture field after cover generates:
Wherein, θ, φ are spherical coordinates angle of the point of observation in rectangular coordinate system O-xyz, k0For free-space propagation constant,
By formulaIt calculates, λ is that the wavelength of antenna according to operating frequency of antenna f and light velocity c passes through formulaIt calculates
It arrives, s is the area of integral unit.
(6b) draws the Antenna Far Field side after covering according to the far field F ' (θ, φ) generated through the aperture field after antenna house
Xiang Tu, and gain G is extracted from directional diagram2With main beam position B2These electrical performance indexes.
(6c) calculates gain loss TL caused by antenna house according to the electrical performance indexes calculated in step (4) and step (5)
With collimating fault BSE:
TL=G1- G2
BSE=| B1- B2|
Step (7), is established and solving optimization designs a model.
(7a) using the thickness value at the discrete point in step (1) as design variable, with the variance of thickness in step (3), most
Gain loss caused by antenna house and collimating fault are design object in big slope, the variance of slope and step (6), are established as follows
Consider the aircraft antenna house Varying-thickness Optimized model of thickness control:
Wherein BSEmaxIt is the maximum value of the collimating fault caused by antenna house under all operating conditions, TLmaxIt is in all operating conditions
The maximum value of gain loss caused by lower antenna house, vd are the variances of design variable, and mds is the maximum of the slope of design variable
Value, vds is the variance of the slope of design variable, BSE0,TL0,vd0,mds0And vds0It is the normalization provided according to particular problem
Coefficient, dminIt is the value lower limit of design variable, dmaxIt is the value upper limit of design variable.
(7b) is using particle swarm optimization algorithm (PSO) come to the aircraft antenna house Varying-thickness optimization mould for considering thickness control
Type is solved, and optimal antenna house thickness distribution is obtained.The population scale of particle swarm optimization algorithm is taken as 50, and evolutionary generation takes
It is 200, for inertia weight as evolutionary generation is from 0.9 linear decrease to 0.4, aceleration pulse is taken as 2.
Step (8), judges whether the electrical performance indexes of the antenna house obtained after optimization and thickness distribution meet preset requirement.
According to the permitted electrical performance indexes knots modification of antenna, judges covering and rear system is optimized to antenna house
Electrical performance indexes knots modification and the thickness distribution of antenna house whether meet preset requirement, if it is satisfied, then antenna cover structure
Design scheme is qualified;Otherwise, the discretization method and Optimal Parameters of antenna house are modified, and repeats step (1) and arrives step (8), directly
It is met the requirements to result.
Advantages of the present invention can be further illustrated by following emulation experiment:
1. simulation parameter
Certain aircraft antenna house, shape as shown in figure 4, basal diameter be 0.5 meter, be highly 1 meter, radome material is glass
Glass Steel material, the relative dielectric constant of material are 4, and magnetic loss angle is just cut to 0.015, and cover internal antenna bore is 0.22 meter, work
Frequency is 9.4GHz, and aperture field is constant amplitude with being mutually distributed, and the scanning angular region of antenna is 0 °~90 °.
2. emulation content and result
Accounting for thickness control to above-mentioned aircraft antenna house using the present invention is that thickness optimization designs, simulation result
As shown in Fig. 5, Fig. 6 and Fig. 7, emulation data are as shown in table 1.
In Fig. 5, Fig. 6 and Fig. 7, design 1 represents equal thickness design scheme, and design 2, which represents, utilizes traditional thickness design side
The design scheme that method obtains, 3 representative of design are obtained using the thickness design method provided by the invention for considering thickness control
Design scheme.
The electrical performance indexes of 1 system of table
From above-mentioned data as it can be seen that the electrical property of antenna house has clear improvement, but its generation after using traditional thickness design
Valence is that the thickness change of antenna house is violent, causes very big difficulty to antenna house, use is provided by the invention
After design method, the improvement to electrical performance indexes not only ensure that, while the thickness change of antenna house is gentle, it is easy to process
Manufacture.
Above-mentioned emulation data experiment proves that the present invention can be effectively improved the electrical property of aircraft antenna house, while can reduce
Its manufacture difficulty.
Claims (6)
1. a kind of Thickness Design Method of antenna house, it is characterized in that: including the following steps:
(1) thickness of antenna house is subjected to discretization according to the short transverse of cover;
It is that x/y plane establishes a coordinate system O-xyz using the bottom center of antenna house as origin, bottom surface along the height of cover, cover
Body height in the z-direction, uniformly chooses 10 discrete points in the direction, and coordinate is denoted as z1,z2,…,z10;
(2) according to the requirement of antenna cover structure performance and electrical property, the value range of thickness value is determined, and to the thickness at discrete point
Degree assigns initial value;
(3) according to the thickness value at discrete point, its statistical value is calculated, the variance including variance, greatest gradient and its slope, according to
Thickness value d at discrete point1,d2,…,d10, the variance vd of these numbers is calculated, further according to its coordinate z1,z2,…,z10, meter
Calculate the slope value of adjacent discrete point:
Calculate the maximum value mds and its variance vds of these slope values;
(4) far field of antenna is calculated, and therefrom extracts gain G1With main beam position B1These electrical performance indexes;
(5) according to the structural parameters of antenna house and material parameter, the transmission coefficient of cover is calculated with transmission line theoryAnd according to
Known antenna aperture field E (x, y) calculates through the aperture field after antenna house:
(6) according to through the aperture field E ' (x, y) after antenna house, the far field F ' (θ, φ) with cover antenna system is calculated, is drawn remote
Field pattern, and gain G is extracted from the far-field pattern2With main beam position B2These electrical performance indexes, and then determine antenna
Gain loss TL and collimating fault BSE caused by covering;
(7) using the thickness value at the discrete point in step (1) as design variable, tiltedly with the variance of thickness, maximum in step (3)
Electrical performance indexes in rate, the variance of slope and step (6) are design object, establish mathematical optimization models, excellent using population
Change algorithm and solve this model, obtains the thickness distribution of antenna house;
(8) whether the electrical performance indexes and thickness distribution for judging the antenna house obtained after optimization meet preset requirement, if it is satisfied,
Then antenna cover structure design scheme is qualified, otherwise, modifies the discretization method and optimization algorithm parameter of antenna house, and repeat step
(1) to step (8), until obtaining electrical performance indexes and thickness distribution meets the design scheme of preset requirement.
2. the Thickness Design Method of a kind of antenna house according to claim 1, it is characterized in that: described step (2) basis
The requirement of antenna cover structure performance and electrical property determines the value range of thickness value, carries out according to the following procedure:
(2a) integrally formed aircraft antenna house uses half-wavelength wall thickness, the calculation formula of half-wavelength wall thickness are as follows:
Wherein λ is wavelength, εrIt is the relative dielectric constant of radome material, α is incidence angle of the incident electromagnetic wave to cover wall;
(2b) determines n according to the structural performance requirements of cover, and n is bigger, and cover wall is thicker, and structural behaviour is better, while electrical property
Can be poorer, therefore take the smallest n value for being able to satisfy structural performance requirements;
(2c) is calculated in the scanning process of internal antenna, the minimum value of the electromagnetic wave of aerial radiation to the incidence angle of cover wall
αminWith maximum value αmax,
(2d) determines the minimum value d of Thickness as the following formulaminWith maximum value dmax:
3. the Thickness Design Method of a kind of antenna house according to claim 1, it is characterized in that: the step (4) is specific
Realizing step is:
(4a) indicates the component in the x, y, z direction for the coordinate system established in step (1) with i, j, k respectively, according to known day
Line Aperture field distribution E (x, y) calculates the far field value F (θ, φ) of antenna:
Wherein, θ, φ are spherical coordinates angle of the point of observation in O-xyz,λ is the wavelength of antenna, according to Antenna Operation frequency
Rate f and light velocity c, passes through formulaIt is calculated, s is the area of integral unit;
(4b) draws Antenna Far Field directional diagram according to the far field F (θ, φ) of ideal antenna, and gain G is extracted from directional diagram1And master
Beam position B1These electrical performance indexes.
4. the Thickness Design Method of a kind of antenna house according to claim 1, it is characterized in that: the step (6) is specific
Realizing step is:
(5a) be directed to step (1) establish coordinate system, the component in its x, y, z direction is indicated with i, j, k respectively, it is known that antenna
Aperture field distribution is denoted as E (x, y);
(5b) establishes the geometrical model of antenna house according to the structure type of antenna house, net is arranged in Mercantile Models analysis software
Lattice side length is 0.2 λ, and wherein λ is the wavelength of antenna, carries out grid dividing to model;
(5c) calculates the transmission coefficient on covering at each point with transmission line theory according to the structural parameters and material parameter of antenna house
The geometry and incident aperture field of (5c1) according to antenna house, find out the incident angle α on antenna house at each point and polarization
Normal angle at electromagnetic wave incident line and incidence point is denoted as incident angle α by angle beta, by polarization of electromagnetic wave direction and incidence
The angle of plane is denoted as polarizing angle β, and wherein plane of incidence is made of electromagnetic wave incident line and the normal at incidence point;
(5c2) is according to the thickness value d at discrete point1,d2,…,d10, smooth antenna house is obtained using cubic Bézier curves method
Thickness profile:
D (z)=fcubic-B-spline(d1,...,d10)
(5c3) is according to antenna house thickness d everywhere, relative dielectric constant εr, losstangenttanδ, calculate on covering at each point
Horizontal polarized components transmission coefficientWith perpendicular polarisation components transmission coefficient
Wherein, ZH=cos α,These parameters are intermediate variable;TH、
TVRespectivelyModulus value, ηH、ηVRespectivelyPhase;
(5c4) is according to horizontal polarized components transmission coefficientWith perpendicular polarisation components transmission coefficientObtain main polarization component
Transmission coefficient:
Wherein,For intermediate variable;
(5d) multiplied by the transmission coefficient of its corresponding point position, calculates the aperture field being incident on antenna house through the bore after covering
:
5. the Thickness Design Method of a kind of antenna house according to claim 1, it is characterized in that: the step (6) is specific
Realizing step is:
(6a) is calculated as follows according to being acquired in step (1) through the aperture field E ' (x, y) after antenna house through after antenna house
Aperture field generate far field F ' (θ, φ):
Wherein, θ, φ are spherical coordinates angle of the point of observation in rectangular coordinate system O-xyz, k0For free-space propagation constant, by formulaIt calculates, λ is that the wavelength of antenna according to operating frequency of antenna f and light velocity c passes through formulaIt is calculated, s is
The area of integral unit;
(6b) draws the Antenna Far Field directional diagram after covering according to the far field F ' (θ, φ) generated through the aperture field after antenna house,
And gain G is extracted from directional diagram2With main beam position B2These electrical performance indexes;
(6c) calculates gain loss TL caused by antenna house and takes aim at according to the electrical performance indexes calculated in step (4) and step (5)
Quasi- error BSE:
TL=G1- G2
BSE=| B1- B2|。
6. the Thickness Design Method of a kind of antenna house according to claim 1, it is characterized in that: the step (7) is specific
Realizing step is:
(7a) is oblique with the variance of thickness, maximum in step (3) using the thickness value at the discrete point in step (1) as design variable
Gain loss caused by antenna house and collimating fault are design object in rate, the variance of slope and step (6), establish following consideration
The aircraft antenna house Varying-thickness Optimized model of thickness control:
Wherein, dmin≤di≤dmax, i=1 ..., 10
Wherein BSEmaxIt is the maximum value of the collimating fault caused by antenna house under all operating conditions, TLmaxIt is the day under all operating conditions
The maximum value of gain loss caused by irdome, vd are the variances of design variable, and mds is the maximum value of the slope of design variable, vds
It is the variance of the slope of design variable, BSE0,TL0,vd0,mds0And vds0It is the normalization coefficient provided according to particular problem,
dminIt is the value lower limit of design variable, dmaxIt is the value upper limit of design variable;
(7b) seeks the aircraft antenna house Varying-thickness Optimized model of consideration thickness control using particle swarm optimization algorithm
Solution, obtains optimal antenna house thickness distribution, and the population scale of particle swarm optimization algorithm is taken as 50, and evolutionary generation is taken as 200, is used to
For property weight as evolutionary generation is from 0.9 linear decrease to 0.4, aceleration pulse is taken as 2.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610824698.8A CN106469850B (en) | 2016-09-14 | 2016-09-14 | A kind of Thickness Design Method of antenna house |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610824698.8A CN106469850B (en) | 2016-09-14 | 2016-09-14 | A kind of Thickness Design Method of antenna house |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106469850A CN106469850A (en) | 2017-03-01 |
CN106469850B true CN106469850B (en) | 2019-03-15 |
Family
ID=58230535
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610824698.8A Active CN106469850B (en) | 2016-09-14 | 2016-09-14 | A kind of Thickness Design Method of antenna house |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106469850B (en) |
Families Citing this family (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108920831B (en) * | 2018-06-29 | 2020-11-03 | 西安电子科技大学 | Method for rapidly calculating influence of high-temperature ablation of high-speed aircraft hood on electrical performance of antenna |
CN110728012B (en) * | 2018-06-29 | 2022-11-04 | 中国石油天然气股份有限公司 | Method and device for acquiring well bore storage coefficient |
CN109490836A (en) * | 2018-07-13 | 2019-03-19 | 中国航空工业集团公司济南特种结构研究所 | A kind of honeycomb mixing sandwich radome |
CN108984902B (en) * | 2018-07-16 | 2020-11-03 | 西安电子科技大学 | Rapid thickness design method for aircraft radome based on phase adjustment |
CN110738002B (en) * | 2019-10-18 | 2023-07-04 | 上海无线电设备研究所 | Comprehensive performance design method for microwave/infrared composite radome |
CN110889216B (en) * | 2019-11-20 | 2023-11-14 | 上海无线电设备研究所 | Self-adaptive rapid design method for curved surface frequency selective surface radome |
CN111177890B (en) * | 2019-12-11 | 2022-09-16 | 西安电子科技大学 | Antenna housing electrical performance uncertainty analysis method based on polygonal interval |
CN111144048B (en) * | 2019-12-23 | 2023-04-07 | 北京机电工程研究所 | Design optimization method of wave-transparent, invisible and heat-proof radome |
CN111342208B (en) * | 2020-03-18 | 2022-01-18 | 维沃移动通信有限公司 | Antenna assembly and electronic equipment |
CN113625062B (en) * | 2021-07-29 | 2022-05-17 | 西安电子科技大学 | Antenna housing electrical property estimation method based on Taylor expansion method |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7894925B1 (en) * | 2008-02-27 | 2011-02-22 | Lockheed Martin Corporation | Method for making a seamed radome for an array antenna and radome with optimal seam locations |
CN102590656A (en) * | 2012-01-03 | 2012-07-18 | 西安电子科技大学 | Antenna cover electric property forecasting method based on distant field |
CN102799782A (en) * | 2012-07-17 | 2012-11-28 | 西安电子科技大学 | Electrical performance prediction method for metal truss-type antenna cover in consideration of scattering and transmission |
CN103401070A (en) * | 2013-07-13 | 2013-11-20 | 西安电子科技大学 | Far field-based wall thickness grinding method for thin shell-type radome |
JP2013244621A (en) * | 2012-05-24 | 2013-12-09 | Mitsubishi Electric Corp | Radome manufacturing method and radome |
CN104484537A (en) * | 2014-12-30 | 2015-04-01 | 西安电子科技大学 | Electromechanical integrated optimization method for sectional dimensions of ribs of metal truss type antenna housing |
-
2016
- 2016-09-14 CN CN201610824698.8A patent/CN106469850B/en active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7894925B1 (en) * | 2008-02-27 | 2011-02-22 | Lockheed Martin Corporation | Method for making a seamed radome for an array antenna and radome with optimal seam locations |
CN102590656A (en) * | 2012-01-03 | 2012-07-18 | 西安电子科技大学 | Antenna cover electric property forecasting method based on distant field |
JP2013244621A (en) * | 2012-05-24 | 2013-12-09 | Mitsubishi Electric Corp | Radome manufacturing method and radome |
CN102799782A (en) * | 2012-07-17 | 2012-11-28 | 西安电子科技大学 | Electrical performance prediction method for metal truss-type antenna cover in consideration of scattering and transmission |
CN103401070A (en) * | 2013-07-13 | 2013-11-20 | 西安电子科技大学 | Far field-based wall thickness grinding method for thin shell-type radome |
CN104484537A (en) * | 2014-12-30 | 2015-04-01 | 西安电子科技大学 | Electromechanical integrated optimization method for sectional dimensions of ribs of metal truss type antenna housing |
Non-Patent Citations (2)
Title |
---|
THICKNESS PROFILE DESIGN OF AIRBORNE RADOMES WITH THICKNESS CONTROL;W.Y.Xu;《Fifth Asia International Symposium on Mechatronics(AISM 2015)》;20151010;全文 |
天线罩机械结构因素对其电性能的影响研究;许万业;《中国博士学位论文全文数据库 信息科技辑》;20160315(第3期);全文 |
Also Published As
Publication number | Publication date |
---|---|
CN106469850A (en) | 2017-03-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106469850B (en) | A kind of Thickness Design Method of antenna house | |
CN106654566B (en) | A kind of quick Thickness Design Method of aircraft antenna house | |
CN106294970B (en) | A kind of novel netted antenna structure design method with shaped-beam | |
Yang et al. | Iterative hybrid method for electromagnetic scattering from a 3-D object above a 2-D random dielectric rough surface | |
CN104036078B (en) | Comprehensive design method of array antenna radiation and scattering performance based on installation height | |
CN104484537A (en) | Electromechanical integrated optimization method for sectional dimensions of ribs of metal truss type antenna housing | |
CN109472066B (en) | Reflecting surface antenna random error analysis method based on unit central point displacement | |
CN106252873B (en) | A kind of Novel Interval Methods of conformal bearer Antenna Far Field power radiation pattern | |
CN102590656A (en) | Antenna cover electric property forecasting method based on distant field | |
Zhu et al. | Multiplexing the aperture of a metasurface: inverse design via deep-learning-forward genetic algorithm | |
CN107404011A (en) | The airborne multiaerial system of low frequency and its design method of a kind of feature based theory of modules | |
CN103745060B (en) | Large-scale antenna dome electrical performance compensation method based on reflecting surface figuration | |
CN109616774A (en) | A kind of absorbing material based on spatial distribution PML model and miniature darkroom | |
CN108984902A (en) | The quick Thickness Design Method of aircraft antenna house based on phase adjustment | |
CN108009355B (en) | Darkroom spherical array compact field quiet zone characteristic spectrum analysis method | |
Wang et al. | An efficient integral equation/modified surface integration method for analysis of antenna-radome structures in receiving mode | |
CN111339606B (en) | Wing shielding effect calculation method based on diffraction principle | |
CN108446504A (en) | Near-field array Antenna measuring table method based on convex optimization | |
CN105740515B (en) | A kind of mixed topology optimization method based on sensitivity analysis | |
Aye et al. | Rectangular microstrip patch antenna array for RFID application using 2.45 GHz frequency range | |
Li et al. | Achievement of beam steering in terahertz band based on frequency-scanning grating-reflector antenna | |
CN202050050U (en) | Flat lens antenna based on electric resonant structure | |
Li et al. | Investigation of leaky-wave propagation and radiation in a metal cut-wire array | |
CN103928763A (en) | Multi-beam antenna based on transformation optics | |
Ataloglou et al. | Synthesis of modulated dielectric metasurfaces for precise antenna beamforming |
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 |