CN102708257B - Method for fast determining tolerance of active phased array antenna structure - Google Patents

Method for fast determining tolerance of active phased array antenna structure Download PDF

Info

Publication number
CN102708257B
CN102708257B CN201210162446.5A CN201210162446A CN102708257B CN 102708257 B CN102708257 B CN 102708257B CN 201210162446 A CN201210162446 A CN 201210162446A CN 102708257 B CN102708257 B CN 102708257B
Authority
CN
China
Prior art keywords
antenna
radiating
tolerance
error
cos
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.)
Expired - Fee Related
Application number
CN201210162446.5A
Other languages
Chinese (zh)
Other versions
CN102708257A (en
Inventor
王从思
康明魁
段宝岩
王伟
徐慧娟
王猛
黄进
保宏
朱敏波
陈光达
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xidian Univ
Original Assignee
Xidian Univ
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Xidian Univ filed Critical Xidian Univ
Priority to CN201210162446.5A priority Critical patent/CN102708257B/en
Publication of CN102708257A publication Critical patent/CN102708257A/en
Application granted granted Critical
Publication of CN102708257B publication Critical patent/CN102708257B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02EREDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation

Abstract

The invention belongs to the technical field of radar antennas, and specifically relates to a method for fast determining tolerance of an active phased array antenna structure. The method comprises the steps of: 1, determining information of an antenna array surface radiation unit; 2, obtaining a height error of the array surface radiation unit; 3, obtaining a position mounting error of the array surface radiation unit; 4, calculating a phase error of an antenna port; 5, calculating electric field distribution of a far field of the antenna; 6, calculating an electric performance parameter related to the antenna; 7, judging and calculating whether the electric performance of the antenna meets the requirements under the current structure tolerance condition according to the requirements of tactical and technical targets, if not, carrying out the step 8 and the step 9, and if so, directly turning to the step 10, wherein the tolerance of the position and the height of the current radiation unit is the fast determined tolerance of the antenna structure. The problem of difficulty in fast determining the tolerance on the distribution antenna structure is effectively solved.

Description

A kind of fast determination method of active phase array antenna structure tolerance
Technical field
The invention belongs to Radar Antenna System field, be specifically related to a kind of fast determination method of active phase array antenna structure tolerance, can be used for instructing determining fast and the evaluation of organization plan of active phase array antenna structure tolerance.
Background technology
Active phase array antenna is in recent years just in the new technology of accelerated development, can meet antenna high-performance, failure-survival capability requirement simultaneously, is also the gordian technique that meets modern radar growth requirement.China develops and formally equipped in 2007 alert-2000 early warning planes of sky of China air force voluntarily as the No.1 military priority project of country second to none of air force, and the breakthrough of its most critical is exactly to have equipped Active Phased Array Radar.Along with radar is day by day urgent to increasing the demands such as operating distance, raising Stealth Fighter, active phase array antenna electrical property is just suggested higher requirement.The state-of-the-art fourth generation fighter that China is developing at present---destroy 20 for meeting the tactical and technical norms such as large detection range, high stealth, also will equip high performance Active Phased Array Radar.
And active phase array antenna electrical property not only depends on the amplitude-phase error of feed system, be more vulnerable to the site error impact of structure radiating element, and this is limited by structural design.As GBR ground radar in American National missile defense systems (NMD) has adopted large-scale flat surface active phased array antenna physique just; its radar weighs 1000 tons; and U.S. SBX sea base radar is as the current unique platform that can bring into play multiple important detecting function in ballistic missile defense system; its radar weighs 2000 tons; for supporting so heavy radar equipment; and meet radar antenna electrical performance indexes, and guarantee antenna array precision, Antenna Construction Design has been proposed to extremely strict requirement.
Meanwhile, because the athletic meeting of antenna mounting platform brings vibrational excitation, airborne radar for example, thereby cause, active phase array antenna is lower in working order there is malformation, causes structure radiating element error.Known, the front malformation that the manufacture of active phase array antenna structure, alignment error and environmental load cause, all will cause the position of distance between antenna array radiating element and each radiating element that relative variation occurs, thereby cause that radiating element exciting current and aperture field PHASE DISTRIBUTION etc. change, so that cause that antenna gain declines, minor level raises and beam position inaccurate etc., finally can reduce the electrical property of active phase array antenna.For this reason, how according to antenna electric performance index request, determine structure tolerance, carry out structural design, and evaluation structure scheme, be a difficult problem that will inevitably run in the development of high-performance active phase array antenna.
For solving this difficult problem, common way is the installation the way of restraint according to active phase array antenna front, suppose that the known distortion of certain rule occurs antenna array structure, typical form of distortion has flexural deformation, bowl-shape distortion, saddle facial disfigurement etc., and then according to the situation of change of antenna electric performance, provide the maximum deformation quantity of front, and then during definite structure tolerance, concrete visible Analysis of performance of active phased array antennas with distorted plane error.International Journal of Electronics, 2009, 96 volumes, 5 phases, 549-559..This method exists main problem to be, constrained when front require to be installed or in front center line, or front surrounding, or four jiaos of fronts, otherwise front deformation rule is difficult to adopt clear and definite mathematical formulae to be described.This causes its range of application to be very limited, in engineering, is difficult to promote, because antenna deformation form is unknown in reality, be not only three kinds of deformation patterns above.Other way, make exactly active phase array antenna exemplar, according to actual measurement electrical property result, through too much wheel scheme, repeatedly carry out to determine structure tolerance, this brings maximum problem is exactly that the whole solution formulation cycle is long, design cost is high, can not meet the demand of current China radar rapid development.
Summary of the invention
The object of the invention is to propose a kind of fast determination method of active phase array antenna structure tolerance, to effectively solve, when active phase array antenna Scheme design, be difficult to determine fast and the problem of distributing antenna structure tolerance, can instruct formulation and the distribution of active phase array antenna structure processing technology and installation accuracy, and the evaluation of antenna structure scheme.
The technical solution that realizes the object of the invention is that a kind of fast determination method of active phase array antenna structure tolerance, is characterized in that: at least comprise the steps:
Step 1, determines antenna array radiating element information;
Step 2, obtains the height error of front radiating element;
Step 3, obtains the position alignment error of front radiating element;
Step 4, calculates Antenna aperture phase error;
Step 5, calculates antenna far field Electric Field Distribution;
Step 6, calculates the relevant unit for electrical property parameters of antenna;
Step 7, according to tactical and technical norms requirement, whether the antenna electric performance that judgement calculates under current structure tolerance conditions meets the demands, if do not meet the demands, carries out step 8 and step 9; If meet the demands, directly forward step 10 to;
Step 8, revises front flatness;
Step 9, revises front installation accuracy, and re-starts step 2 to step 7;
Step 10, the position of current radiating element and height tolerance be exactly definite antenna structure tolerance fast.
Described step 1 determines that antenna array radiating element information comprises the steps:
1.1. obtain antenna length L x, width L y, frequency of operation f and wavelength X information, obtain line number M, columns N and the spacing d of radiating element in x, y direction of front radiating element x, d yinformation; Being numbered of radiating element (i, j) in front;
1.2. the information of above-mentioned antenna array radiating element is formed to the data file of set form according to the order of row, column numbering.
Described step 2 is obtained the height error of front radiating element;
2.1. making the front flatness average that distributing is 0, and standard deviation is σ znormal distribution stochastic error, choose prompt radiation cell height error, i.e. σ zit is 1 times of wavelength X;
2.2. making the design coordinate of (i, j) individual radiating element is (id x, jd y, 0), when front flatness z is when existing stochastic error Δ z, radiating element coordinate becomes:
(i·d x,j·d y,Δz ij) (1)。
Described step 3 is obtained the position alignment error of front radiating element;
3.1. making the front installation accuracy average that distributing is 0, and standard deviation is σ xynormal distribution stochastic error, choose prompt radiation cell position alignment error, i.e. σ x, σ yit is 1 times of wavelength X;
3.2. when radiating element (i, j) position alignment error exists after stochastic error Δ x and Δ y in x, y direction, the corresponding 0xy plane of radiating element internal coordinate becomes
(i·d x+Δx ij,j·d y+Δy ij) (2);
When the site error amount of radiating element is (Δ x ij, Δ y ij, Δ z ij), the new coordinate of radiating element is:
(i·d x+Δx ij,j·d y+Δy ij,Δz ij) (3)。
Described step 4 is the results according to step 2 and step 3, calculates Antenna aperture phase error:
4.1. target setting is expressed as (cos α with respect to the direction (θ, φ) at coordinate system O-xyz place with direction cosine x, cos α y, cos α z), target with respect to the angle of coordinate axis and the pass of direction cosine is
cos α x = sin θ cos φ cos α y = sin θ sin φ cos α z = cos θ - - - ( 4 ) ;
4.2. between adjacent two radiating elements of antenna, at target place (θ, φ), along the space quadrature of x-axis, y-axis and z-axis, be respectively:
Wherein, wave constant k=2 π/λ
Therefore, (m, n) individual radiating element (0≤m≤M-1,0≤n≤N-1) with respect to the phase differential of (0,0) individual radiating element is:
Wherein, (Δ x 00, Δ y 00, Δ z 00) be the displacement of (0,0) individual radiating element, β mnbe phase differential in (m, n) individual radiating element battle array, by phase shifter, control and realize beam scanning;
4.3. the phase differential with respect to reference radiation unit by each radiating element, becomes the form of matrix according to the sequential storage of radiating element Position Number, this matrix represents the phase differential of active phase array antenna actinal surface.
Described step 5 is the result of calculation according to step 4, calculates antenna far field Electric Field Distribution:
Draw the electric field value of Antenna Far Field region point (θ, φ); Change (θ, φ) numerical value, double counting process, can draw in certain concrete scope of far-field region electric field value a little, field value is taken the logarithm, and is plotted in coordinate system, can obtain the directional diagram of far field regional extent.
Described step 6 is the result of calculation according to step 5, calculates the relevant unit for electrical property parameters of antenna and comprises the steps:
6.1. from directional diagram, obtain the electrical properties such as minor level SLL, beam position BS and beam angle HPBW;
6.2. according to Electric Field Distribution E (θ, φ), calculate antenna gain
G = 10 log 10 4 π ∫ 0 2 π ∫ 0 π [ E ( θ , φ ) ] 2 sin θdθdφ - - - ( 8 ) ;
6.3. based on Antenna Design index, calculate the reducing amount of the electrical properties such as antenna gain G, minor level SLL, beam position BS.
Described step 8: revise front flatness, make standard deviation sigma zmultiple for wavelength X is distinguished value successively in the following order
[ 1 2 , 1 4 , 1 5 , 1 8 , 1 10 , 1 15 , 1 20 , 1 25 , 1 30 , 1 40 , 1 60 , 1 80 , 1 100 ] - - - ( 9 ) .
Described step 9: revise front installation accuracy, make standard deviation 0 x, σ yxy) for the multiple of wavelength X, distinguish in the following order successively value
[ 1 2 , 1 4 , 1 5 , 1 8 , 1 10 , 1 15 , 1 20 , 1 25 , 1 30 , 1 40 , 1 60 , 1 80 , 1 100 ] - - - ( 10 ) .
The present invention is due to the site error of front radiating element is incorporated in the phase error of antenna opening diametric plane by phase differential, antenna structure displacement field and electromagnetic field are closely connected, thereby avoided only by structure tolerance, judging the deficiency of antenna performance, thereby realized the electrical and mechanical comprehensive analysis of antenna; Antenna electric performance computing method owing to having built that front radiating element exists arbitrary height error, position alignment error, thereby can realize by random generation radiating element site error the evaluation of electrical property, avoided prior art only can analyze the problem that is difficult to engineering application that three kinds of front constraint types (flexural deformation, bowl-shape distortion, saddle facial disfigurement) cause, precision and speed that antenna electric performance calculates have been improved simultaneously, and the active phase array antenna structure that can analyze different frequency range, has good applicability.
Simulation result shows, structure tolerance when method of the present invention can realize the active phase array antenna solution formulation of different frequency range is determined fast, instruct formulation and the distribution of active phase array antenna structure processing technology and installation accuracy, and the evaluation of antenna structure scheme, make determining of antenna array structure tolerance more have explanation of force and feasibility.
Accompanying drawing explanation
Referring to accompanying drawing, the present invention is described in further detail:
Fig. 1 is that active phase array antenna structure tolerance of the present invention is determined process flow diagram;
Fig. 2 is that the radiating element of planar rectangular active phase array antenna is arranged schematic diagram;
Fig. 3 is radiating element height error schematic diagram;
Fig. 4 is radiating element position alignment error schematic diagram;
Fig. 5 is the space geometry graph of a relation of target;
Fig. 6 is flat surface active phased array antenna front schematic diagram;
Fig. 7 is the antenna radiation pattern that Shi Butong front flatness is corresponding;
Fig. 8 is the antenna radiation pattern that Shi Butong front flatness is corresponding;
Fig. 9 is the antenna radiation pattern that Shi Butong front installation accuracy is corresponding;
Figure 10 is the antenna radiation pattern that Shi Butong front installation accuracy is corresponding.
In figure, sequence number represents respectively: 1, antenna array radiating element; 2, antenna length; 3, day line width; 4, line number; 5, columns; 6, the spacing of radiating element in x direction; 7, the spacing of radiating element in y direction.
Embodiment
(with reference to Fig. 1) is as follows for concrete steps of the present invention:
Step 1, determines antenna array radiating element information.
1.1. obtain antenna length L x, width L y, frequency of operation f and wavelength X information, obtain line number M, columns N and the spacing d of radiating element in x, y direction of front radiating element x, d yinformation (see figure 2); Being numbered of radiating element (i, j) in front; Wherein i, j are the radiating element numbering in x, y direction respectively, and front lower right-hand corner is Base Serial Number, and the radiating element of front lower right-hand corner is numbered (0,0), and this is also the true origin that is positioned at the coordinate system Oxy of front simultaneously; Front normal direction is exactly the z axle of coordinate system Oxyz;
1.2. the information of above-mentioned antenna array radiating element is formed to the data file of set form according to the order of row, column numbering.
Step 2, obtains the height error of front radiating element.
2.1. active phase array antenna structure tolerance shows as front flatness and front installation accuracy at front, and front flatness is definite to height error by the z of radiating element;
2.2. making the front flatness average that distributing is 0, and standard deviation is σ znormal distribution stochastic error, choose prompt radiation cell height error, i.e. σ zit is 1 times of wavelength X;
2.3. making the design coordinate of (i, j) individual radiating element is (id x, jd y, 0), when front flatness z (see figure 3) when there is stochastic error Δ z, radiating element coordinate becomes:
(i·d x,j·d y,Δz ij) (11)。
Step 3, obtains the position alignment error of front radiating element.
3.1. front installation accuracy is determined by radiating element position alignment error in x, y direction, and the mounting means of radiating element is identical in x, y direction, therefore its installation accuracy is identical;
3.2. making the front installation accuracy average that distributing is 0, and standard deviation is σ xynormal distribution stochastic error, choose prompt radiation cell position alignment error, i.e. σ x, σ yit is 1 times of wavelength X;
3.3. (see figure 4) after radiating element (i, j) position alignment error exists stochastic error Δ x and Δ y in x, y direction, the corresponding 0xy plane of radiating element internal coordinate becomes
(i·d x+Δx ij,j·d y+Δy ij) (12);
When the site error amount of radiating element is (Δ x ij, Δ y ij, Δ x ij), the new coordinate of radiating element is:
(i·d x+Δx ij,j·d y+Δy ij,Δz ij) (13)。
Step 4, according to the result of step 2 and step 3, calculates Antenna aperture phase error.
4.1. target setting is expressed as (cos α with respect to the direction (θ, φ) at coordinate system O-xyz place with direction cosine x, cos α y, cos α z), according to the known target of Fig. 5, with respect to the angle of coordinate axis and the pass of direction cosine, be
cos α x = sin θ cos φ cos α y = sin θ sin φ cos α z = cos θ - - - ( 14 ) ;
4.2. between adjacent two radiating elements of antenna, at target place (θ, φ), along the space quadrature of x-axis, y-axis and z-axis, be respectively:
Wherein, wave constant k=2 π/λ
Therefore, (m, n) individual radiating element (0≤m≤M-1,0≤n≤N-1) with respect to the phase differential of (0,0) individual radiating element is:
Wherein, (Δ x 00, Δ y 00, Δ z 00) be the displacement of (0,0) individual radiating element, β mnbe phase differential in (m, n) individual radiating element battle array, by phase shifter, control and realize beam scanning;
4.3. the phase differential with respect to reference radiation unit by each radiating element, becomes the form of matrix according to the sequential storage of radiating element Position Number, this matrix represents the phase differential of active phase array antenna actinal surface.
Step 5, according to the result of calculation of step 4, calculates antenna far field Electric Field Distribution.
5.1. because changing, radiating element spacing will make radiating element the mutual coupling coefficient S mn, pqchange into S mn, pq(Δ x, Δ y), therefore (m, n) individual radiating element directional diagram is
f mn ( θ , φ ) = f e ( θ , φ ) · Σ p = 0 M - 1 Σ q = 0 N - 1 S mn , pq ′ ( Δx , Δy ) · expjk { [ ( p - m ) dx + Δx pq - Δ x mn ] · cos α x + - - - ( 17 ) ;
[ ( q - n ) dy + Δy pq - Δy mn ] · cos α y }
Wherein, f e(θ, φ) is radiating element ideal orientation figure (all radiating elements are identical here), S ' mn, pq(Δ x, Δ y) is the mutual coupling coefficient of (p, q) individual radiating element to (m, n) individual radiating element, and its value is
Here radiating element the mutual coupling coefficient can utilize the Electromagnetic Simulation softwares such as HFSS to calculate.
The exciting current amplitude of 5.2. applying to each radiating element in feeding network is phase place is β mn, therefore, consider that the radiating element exciting current amplitude after mutual coupling is
I mn = I mn e + Σ p = 0 M - 1 Σ q = 0 N - 1 S mn , pq ( Δx , Δy ) · I pq - - - ( 19 ) ;
5.3. therefore according to formula (17), (18), (19), calculate radiating element and exist the far field Electric Field Distribution of the active phase array antenna of height error and position alignment error to be
5.4. according to calculating antenna far field Electric Field Distribution, draw the far field direction of an electric field figure of antenna.Draw the electric field value of Antenna Far Field region point (θ, φ); Change (θ, φ) numerical value, double counting process, can draw in certain concrete scope of far-field region electric field value a little, field value is taken the logarithm, and is plotted in coordinate system, can obtain the directional diagram of far field regional extent.
Step 6, according to the result of calculation of step 5, calculates the relevant unit for electrical property parameters of antenna.
6.1. from directional diagram, obtain the electrical properties such as minor level SLL, beam position BS and beam angle HPBW;
From directional diagram, can directly read out the unit for electrical property parameters such as controlling antenna wave beam to point BS and beam angle HPBW, to former directional diagram normalization, can obtain normalized radiation pattern, therefrom can read minor level SLL.As for other electrical property, as efficiency, on the basis of existing electric field value, by corresponding calculating, can obtain.
6.2. according to Electric Field Distribution E (θ, φ), calculate antenna gain G
G = 10 log 10 4 π ∫ 0 2 π ∫ 0 π [ E ( θ , φ ) ] 2 sin θdθdφ - - - ( 21 ) ;
6.3. based on Antenna Design index, calculate reducing amount or the deterioration degree of the electrical properties such as antenna gain G, minor level SLL, beam position BS.
Step 8, revises front flatness, makes standard deviation sigma zmultiple for wavelength X is distinguished value successively in the following order
[ 1 2 , 1 4 , 1 5 , 1 8 , 1 10 , 1 15 , 1 20 , 1 25 , 1 30 , 1 40 , 1 60 , 1 80 , 1 100 ] - - - ( 22 ) .
Step 9, revises front installation accuracy, makes standard deviation sigma x, σ yxy) for the multiple of wavelength X, distinguish in the following order successively value
[ 1 2 , 1 4 , 1 5 , 1 8 , 1 10 , 1 15 , 1 20 , 1 25 , 1 30 , 1 40 , 1 60 , 1 80 , 1 100 ] - - - ( 23 ) .
Advantage of the present invention can further illustrate by following emulation experiment:
1. simulated conditions:
Adopt the X frequency range planar rectangular active phase array antenna schematic diagram providing as Fig. 6, the front radiating element 1 that distributing on the antenna of length 2 and width 3, its line number 4 is that M, columns 5 are N, the spacing 6 of radiating element in x direction is d x; The spacing 7 of radiating element in y direction is d y.
Whole antenna package is containing 1024 radiating elements, and design parameter is as shown in table 1, and supposes the even weighting of the exciting current employing constant amplitude homophase of Antenna aperture.
Table 1 antenna parameter
There is the antenna electric performance in height error and two kinds of structure tolerance situations of position alignment error in emulation front radiating element, a kind of is z impact on electrical property to different front flatnesses, and another kind is x, the y impact on electrical property to different front installation accuracies.
2. simulation result:
Because antenna is operated in high-frequency band, so according to actual processing and manufacturing ability, the standard deviation Initial value choice of stochastic error is 3.2mm,
(1) in front z direction, add standard deviation sigma zbe respectively normal distribution stochastic error Δ z after, emulation obtains with as shown in Figure 7 and Figure 8, corresponding unit for electrical property parameters is as shown in table 2 for the antenna radiation pattern of plane.
(2) in front x, y direction, add identical standard poor normal distribution stochastic error Δ x, after Δ y, emulation obtains with as shown in Figure 9 and Figure 10, corresponding unit for electrical property parameters is as shown in table 3 for the antenna radiation pattern of plane.
Antenna electric performance parameter under table 2 Different Plane degree tolerance
Antenna electric performance parameter under table 3 diverse location location tolerance
For this 32x32 active phase array antenna that works in 9.375GHz: from Fig. 7, Fig. 8 and table 2, flatness tolerance is larger, antenna gain loss is larger, and when flatness tolerance is λ/15(2.13mm) time, antenna gain loss is 0.454dB, minor level is raised as 0.447dB, meets the requirement that is less than 0.5dB in engineering for gain loss; From Fig. 9, Figure 10 and table 3, position location tolerance is larger, and antenna gain loss is larger, and when position location tolerance is λ/20(1.6mm) time, antenna gain loss is 0.471dB, and minor level is raised as 0.114dB, meets the requirement that is less than 0.5dB in engineering for gain loss.Data in table, the minor level deterioration degree under known Different Plane degree tolerance will be apparently higher than corresponding position location tolerance, and flatness and the installation accuracy of therefore tackling active phase array antenna structure tolerance propose respectively rational requirement.
Above-mentioned simulation numerical evidence, adopts the present invention can effectively determine fast the structure tolerance of active phase array antenna.

Claims (8)

1. a fast determination method for active phase array antenna structure tolerance, is characterized in that: at least comprise the steps:
Step 1, determines antenna array radiating element information;
Step 2, obtains the height error of front radiating element;
Step 3, obtains the position alignment error of front radiating element;
Step 4, calculates Antenna aperture phase error;
Step 5, calculates antenna far field Electric Field Distribution;
Step 6, calculates the relevant unit for electrical property parameters of antenna;
Step 7, according to tactical and technical norms requirement, whether the antenna electric performance that judgement calculates under current structure tolerance conditions meets the demands, if do not meet the demands, carries out step 8 and step 9; If meet the demands, directly forward step 10 to;
Step 8, revises front flatness;
Step 9, revises front installation accuracy, and re-starts step 2 to step 7;
Step 10, the position of current radiating element and height tolerance be exactly definite antenna structure tolerance fast;
In described step 1, determine that antenna array radiating element information comprises the steps:
1.1. obtain antenna length L x, width L y, frequency of operation f and wavelength X information, obtain line number M, columns N and the spacing d of radiating element in x, y direction of front radiating element x, d yinformation; Being numbered of radiating element (i, j) in front;
1.2. the information of above-mentioned antenna array radiating element is formed to the data file of set form according to the order of row, column numbering.
2. the fast determination method of a kind of active phase array antenna structure tolerance according to claim 1, is characterized in that: described step 2 is obtained the height error of front radiating element;
2.1. making the front flatness average that distributing is 0, and standard deviation is σ znormal distribution stochastic error, choose prompt radiation cell height error, i.e. σ zit is 1 times of wavelength X;
2.2. making the design coordinate of (i, j) individual radiating element is (id x, jd y, 0), when front flatness z is when existing stochastic error Δ z, radiating element coordinate becomes:
(i·d x,j·d y,Δz ij) (1)。
3. the fast determination method of a kind of active phase array antenna structure tolerance according to claim 1, is characterized in that: described step 3 is obtained the position alignment error of front radiating element;
3.1. making the front installation accuracy average that distributing is 0, and standard deviation is σ xynormal distribution stochastic error, choose prompt radiation cell position alignment error, i.e. σ x, σ yit is 1 times of wavelength X;
3.2. when radiating element (i, j) position alignment error exists after stochastic error Δ x and Δ y in x, y direction, the corresponding Oxy plane of radiating element internal coordinate becomes
(i·d x+Δx ij,j·d y+Δy ij) (2);
When the site error amount of radiating element is (Δ x ij, Δ y ij, Δ z ij), the new coordinate of radiating element is:
(i·d x+Δx ij,j·d y+Δy ij,Δz ij) (3)。
4. the fast determination method of a kind of active phase array antenna structure tolerance according to claim 1, is characterized in that: described step 4 is the results according to step 2 and step 3, calculates Antenna aperture phase error:
4.1. target setting is expressed as (cos α with respect to the direction (θ, φ) at coordinate system O-xyz place with direction cosine x, cos α y, cos α z), target with respect to the angle of coordinate axis and the pass of direction cosine is
cos α x = sin θ cos φ cos α y = sin θ sin φ cos α z = cos θ - - - ( 4 ) ;
4.2. between adjacent two radiating elements of antenna, at target place (θ, φ), along the space quadrature of x-axis, y-axis and z-axis, be respectively:
Wherein, wave constant k=2 π λ
Therefore, (m, n) individual radiating element (0≤m≤M-1,0≤n≤N-1) with respect to the phase differential of (0,0) individual radiating element is:
Wherein, (Δ x 00, Δ y 00, Δ z 00) be the displacement of (0,0) individual radiating element, β mnbe phase differential in (m, n) individual radiating element battle array, by phase shifter, control and realize beam scanning;
4.3. the phase differential with respect to reference radiation unit by each radiating element, becomes the form of matrix according to the sequential storage of radiating element Position Number, this matrix represents the phase differential of active phase array antenna actinal surface.
5. the fast determination method of a kind of active phase array antenna structure tolerance according to claim 1, is characterized in that: described step 5 is the result of calculation according to step 4, calculates antenna far field Electric Field Distribution:
Draw the electric field value of Antenna Far Field region point (θ, φ); Change (θ, φ) numerical value, double counting process, can draw in certain concrete scope of far-field region electric field value a little, field value is taken the logarithm, and is plotted in coordinate system, can obtain the directional diagram of far field regional extent; Wherein, I mnfor the radiating element exciting current amplitude after mutual coupling; f mn(θ, φ) is (m, n) individual radiating element direction.
6. the fast determination method of a kind of active phase array antenna structure tolerance according to claim 1, is characterized in that: described step 6 is the result of calculation according to step 5, calculates the relevant unit for electrical property parameters of antenna and comprises the steps:
6.1. from directional diagram, obtain the electrical properties such as minor level SLL, beam position BS and beam angle HPBW;
6.2. according to Electric Field Distribution E (θ, φ), calculate antenna gain
G = 10 log 10 4 π ∫ 0 2 π ∫ 0 π E ( θ , φ ) ] 2 sin θdθdφ - - - ( 8 ) ;
6.3. based on Antenna Design index, calculate reducing amount or the deterioration degree of the electrical properties such as antenna gain G, minor level SLL, beam position BS.
7. the fast determination method of a kind of active phase array antenna structure tolerance according to claim 1, is characterized in that: described step 8: revise front flatness, make standard deviation sigma zmultiple for wavelength X is distinguished value successively in the following order
[ 1 2 , 1 4 , 1 5 , 1 8 , 1 10 , 1 15 , 1 20 , 1 25 , 1 30 , 1 40 , 1 60 , 1 80 , 1 100 ] - - - ( 9 ) .
8. the fast determination method of a kind of active phase array antenna structure tolerance according to claim 1, is characterized in that: described step 9: revise front installation accuracy, make standard deviation sigma x, σ yxy) for the multiple of wavelength X, distinguish in the following order successively value
[ 1 2 , 1 4 , 1 5 , 1 8 , 1 10 , 1 15 , 1 20 , 1 25 , 1 30 , 1 40 , 1 60 , 1 80 , 1 100 ] - - - ( 10 ) .
CN201210162446.5A 2012-05-23 2012-05-23 Method for fast determining tolerance of active phased array antenna structure Expired - Fee Related CN102708257B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210162446.5A CN102708257B (en) 2012-05-23 2012-05-23 Method for fast determining tolerance of active phased array antenna structure

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210162446.5A CN102708257B (en) 2012-05-23 2012-05-23 Method for fast determining tolerance of active phased array antenna structure

Publications (2)

Publication Number Publication Date
CN102708257A CN102708257A (en) 2012-10-03
CN102708257B true CN102708257B (en) 2014-04-02

Family

ID=46901020

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210162446.5A Expired - Fee Related CN102708257B (en) 2012-05-23 2012-05-23 Method for fast determining tolerance of active phased array antenna structure

Country Status (1)

Country Link
CN (1) CN102708257B (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103116677B (en) * 2013-02-22 2015-08-12 西安电子科技大学 A kind of microchannel cold plates method for designing towards active phase array antenna
CN104063426B (en) * 2014-06-06 2017-02-15 西安电子科技大学 Radiation-and-scattering-oriented quick determination method for active phased-array antenna structure tolerance
CN104036093B (en) * 2014-06-26 2017-01-18 西安电子科技大学 Large deformation array antenna sidelobe performance prediction method based on array element mutual coupling
CN105718662B (en) * 2016-01-21 2018-10-16 西安电子科技大学 Characteristics of conformal array antenna structure tolerance fast determination method based on mechanical-electric coupling
CN109031226B (en) * 2018-08-23 2020-11-03 西安电子科技大学 Electromechanical coupling-based method for rapidly evaluating detection performance of active phased array radar
CN110532631B (en) * 2019-08-01 2021-01-05 西安电子科技大学 6G communication antenna array element position tolerance determination method based on channel capacity sensitivity
CN110488094B (en) * 2019-08-28 2021-08-17 西安电子科技大学 Interval analysis method for influence of manufacturing tolerance of microstrip patch antenna on electrical performance

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101267062A (en) * 2008-04-30 2008-09-17 西安电子科技大学 Method for predicting antenna electric performance based on simulated distortion reflective side

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP4190929B2 (en) * 2003-03-28 2008-12-03 三菱レイヨン株式会社 False twisted yarn and method for producing the same, and twisted yarn and woven / knitted fabric of the false twisted yarn

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101267062A (en) * 2008-04-30 2008-09-17 西安电子科技大学 Method for predicting antenna electric performance based on simulated distortion reflective side

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
JP特开2004-293006A 2004.10.21
切比雪夫带通滤波器的设计;王利众;《凯里学院学报》;20100630;第23-24页 *
王从思等.基于相位差的平面相控阵天线阵元精度分析方法.《2009年全国天线年会论文集(上)》.2009,第755-758页. *
王利众.切比雪夫带通滤波器的设计.《凯里学院学报》.2010,第23-24页.

Also Published As

Publication number Publication date
CN102708257A (en) 2012-10-03

Similar Documents

Publication Publication Date Title
CN102708257B (en) Method for fast determining tolerance of active phased array antenna structure
CN104036114B (en) A kind of fast determination method of the hexagon active phase array antenna structure tolerance based on mechanical-electric coupling
CN104038295B (en) A kind of distortion array antenna scattering property analytical method based on mechanical-electric coupling
CN101344564B (en) Active phase array antenna electrical property prediction method based on mechanical, electric and thermal three-field coupling
CN105740502B (en) A kind of characteristics of conformal array antenna electrical behavior prediction method under deformation under load
Wang et al. Coupled structural-electromagnetic-thermal modelling and analysis of active phased array antennas
CN104182636A (en) Fast implementation method for comprehensive minor lobes of array antenna radiation field and scattered field
CN104701637B (en) Anamorphic array antenna electric performance compensation method based on mechanical-electric coupling Yu least square method
CN104615836A (en) Rapid prediction method for impact of spaceborne microstrip antenna array thermal deformation on electrical performance
CN105718662B (en) Characteristics of conformal array antenna structure tolerance fast determination method based on mechanical-electric coupling
CN104036093B (en) Large deformation array antenna sidelobe performance prediction method based on array element mutual coupling
CN102253290B (en) Method for predicting electrical properties of deformed log-periodic antennae based on electromechanical coupling model
CN106021764B (en) The calculation method of active phase array antenna performance simulation confidence level towards mechanical-electric coupling
CN108110419A (en) A kind of method for correcting radome boresight error
CN106407723B (en) The determination method of sparse arrangement array antenna exciting current amplitude towards Sidelobe
CN104933213A (en) Large-scale phased antenna array wide-angle scanning optimization method based on space mapping
CN101931124A (en) Embattling method of logarithmic spiral array antennas
CN104036078A (en) Comprehensive design method of array antenna radiation and scattering performance based on installation height
CN105161860B (en) Deformation planar array electrical performance compensation method based on mechanical-electric coupling and Fourier transformation
CN104409859A (en) High-gain low-side lobe half-height waveguide wide side slot antenna and design method thereof
CN105787160B (en) The structure lightened design method of spaceborne active phase array antenna based on mechanical-electric coupling
CN106354909B (en) Sparse arrangement array antenna structure public affairs method for determining difference towards radiation and scattering
CN105760600A (en) Method for determining heat power consumption of satellite-borne active phased-array antenna assembly based on electromechanical coupling
CN104063426A (en) Radiation-and-scattering-oriented quick determination method for active phased-array antenna structure tolerance
CN109031226B (en) Electromechanical coupling-based method for rapidly evaluating detection performance of active phased array radar

Legal Events

Date Code Title Description
PB01 Publication
C06 Publication
SE01 Entry into force of request for substantive examination
C10 Entry into substantive examination
GR01 Patent grant
C14 Grant of patent or utility model
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20140402

Termination date: 20150523

EXPY Termination of patent right or utility model