CN101267062A - Method for predicting antenna electric performance based on simulated distortion reflective side - Google Patents

Method for predicting antenna electric performance based on simulated distortion reflective side Download PDF

Info

Publication number
CN101267062A
CN101267062A CNA2008100181069A CN200810018106A CN101267062A CN 101267062 A CN101267062 A CN 101267062A CN A2008100181069 A CNA2008100181069 A CN A2008100181069A CN 200810018106 A CN200810018106 A CN 200810018106A CN 101267062 A CN101267062 A CN 101267062A
Authority
CN
China
Prior art keywords
antenna
sigma
match
coordinate
node
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.)
Granted
Application number
CNA2008100181069A
Other languages
Chinese (zh)
Other versions
CN101267062B (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 CN200810018106A priority Critical patent/CN101267062B/en
Publication of CN101267062A publication Critical patent/CN101267062A/en
Application granted granted Critical
Publication of CN101267062B publication Critical patent/CN101267062B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The present invention provides a prediction method of the antenna electrical property based on the fitting deformation reflector, mainly resolving the problem of electro-mechanical separation in the antenna design. The method includes following steps: obtaining a node displacement after the deforming of the reflector based on the finite element analysis of the antenna structure; designing the space position relationship of the coordination and the deformed coordinate according to the theory of the reflector node, processing fit to the deformed reflector; calculating the normal error of the antenna surface and the phase error caused by the normal error by utilizing the fitting reflector; calculating antenna far fields electric field distribution according to the amplitude and the phase parameter of the antenna aperture; calculating antenna far fields electrical pattern according to the antenna far fields electric field distribution, obtaining the electrical parameter from the direction pattern, realizing the electrical property prediction of the reflector antenna structure displacement field and the electromagnetic field coupling. The invention can guide the structure design of the antenna, and comprehensively analyzes and assesses the antenna electromechanical properties under different conditions.

Description

Method for predicting antenna electric performance based on the match distortion reflective side
Technical field
The invention belongs to antenna technical field, specifically is a kind of method for predicting antenna electric performance based on the match distortion reflective side.Be used in reference to structural design, the performance simulation A+E of lead antenna.
Background technology
Along with the development of communication, radio astronomy cause, surface antenna develops to high band, bigbore direction.Large-scale surface antenna structure is typical electromechanical integrated structure, and its mechanical performance and electrical property influence each other, mutual restriction.In the engineering, the Electrical Engineer proposes the Design of Mechanical Structure requirement, and the structural engineer can only distribute the design accuracy of each building block by rule of thumb.Two kinds of situations appear in the result, and the one, the latter has used up all ways, has used best process equipment and means, still can't meet the demands; The 2nd, in actual production, what the accuracy of manufacture was high is not to satisfy electrical performance indexes, and some accuracy of manufacture and so high instead can not satisfy electrical performance indexes.The result causes antenna manufacturing cost height, and the lead time is long, and its performance can't fundamentally guarantee.
Because design, manufacturing and the testing expense of large-scale surface antenna are very high, so require its design should one-shot forming.But because of the bore of large-scale antenna reaches tens meters, nearly 100 tons of its weight are brought great difficulty to structural design; Deform owing to this large-scale antenna structure is vulnerable to the external environment condition effect very much simultaneously, antenna electric performance is affected.As malformation make antenna efficiency reduce, minor level uprises, directivity variation etc.When the operating frequency of antenna of high band reaches the Ka frequency range, antenna structure distortion to the influence of antenna electric performance with even more serious.Owing to can't determine the quantitative relationship between antenna structure distortion and the electrical property in the prior art, certainly exist dynamo-electric problem of separating when causing the antenna structure design.
At present, solving the most frequently used method of the dynamo-electric separate design problem of antenna both at home and abroad has following several:
(1) from comprehensive angle antenna is carried out integrated analysis, with the thought of optimizing modeling subject design such as each machinery, electromagnetism are required to unify to consider, this method has been considered the benefit of electrical and mechanical comprehensive design.As at J.S.Liu, L.Hollaway.Integrated structure-electromagnetic optimization of large reflector antenna systems.Structuraland Multidisciplinary Optimization, VOL.16, NO.1, the method that is adopted among the July 1998 is exactly this comprehensive optimization method.But it is how to influence antenna electric performance that this method is not fundamentally analyzed the antenna structure distortion, promptly can not satisfy the scheme that provides reduction structural design difficulty under the electrical performance indexes prerequisite.
(2) utilize the antenna reflective face warping function, obtain of the contribution of each node to the antenna electric field, thereby analyze the antenna electric performance situation of change under the different distortion situation, as at K.Bahadori, Y.Rahmat-samii.Characterization of effects ofperiodic and aperiodic surface distortions on membrane reflector antennas.IEEE Trans.Antennasand Propagation, VOL.53, NO.9, the method that is adopted among the September 2005 is exactly this method.This method only is that the hypothesis distorted reflector satisfies certain trigonometric function distribution, but malformation is difficult to provide with a certain concrete function in the reality.The electromechanical properties analysis-by-synthesis of this method is to be based upon on the basis of malformation shape hypothesis simultaneously, can not reflect the true influence relation between antenna structure distortion and the antenna electric performance.
(3) measurement point on the antenna deformation curved surface in the actual engineering of employing, and theoretical node simulation analysis distortion back is as calculating object, analyze the influence of antenna deformation, as (Hua Mulin) just adopting this method in the document in " modern radar " 1994 the 1st phases " a kind of approximating method of antenna deformation curved surface " to antenna electric performance.The method engineering using value is big, but key is the antenna that actual processing will be arranged, assemble, and needs to carry out on the antenna material object Measurement and analysis.General antenna structure designer need know antenna electric performance under the current structure in the design of Simulation stage, and need to judge whether change in view of the above or redesign antenna structure, and can not determine at antenna structure, reflecting surface shapes, assembles under the also completed situation, analyzes the electromechanical properties of antenna again.
The content of invention
The objective of the invention is to avoid the deficiency of above-mentioned art methods, propose a kind of method for predicting antenna electric performance, instruct the Electromechanical Design of antenna structure, to reduce design cost, to improve antenna electrical and mechanical comprehensive performance based on the match distortion reflective side.
The technical scheme that realizes the object of the invention is, based on the antenna structure finite element analysis, obtain the modal displacement after the distorted reflector, Design Theory coordinate and distortion recoil target spatial relation according to the reflecting surface node, distortion reflective side is carried out match, and calculate the phase error of the bore face that antenna surface normal error and this method phase error cause, calculate antenna far field Electric Field Distribution, draw the directional diagram of antenna far field electric field, carry out surface antenna electromechanical properties analysis-by-synthesis.Detailed process is as follows:
(1) according to Reflector Panel, back of the body frame, the centerbody structural parameters of antenna, determines the antenna structure FEM (finite element) model, obtain the theoretical coordinate P (x of reflecting surface sampling node i, y i, z i);
(2) utilize finite element software, the antenna reflective face structure is carried out finite element analysis, obtain displacement P (the Δ x of distorted reflector post-sampling node i, Δ y i, Δ z i), and carry out Coordinate Conversion by the structural model coordinate system, obtain the new coordinate system of antenna reflective face;
(3) in new coordinate system, utilize distorted reflector posterior nodal point coordinate and the minimum difference that is out of shape the front nodal point coordinate, construct the paraboloidal equation group of match, obtain the equation of match distortion reflective side;
(4) utilize match distortion reflective side equation, obtain distorted reflector posterior nodal point P (x i+ Δ x i, y i+ Δ y i, z i+ Δ z i) normal direction cosine and normal direction deviation ε i, calculate the phase error δ that each node causes at the bore face i
(5) the antenna opening diametric plane is divided into N ring territory, calculates n the K on the ring territory nThe bore face phase error of individual node correspondence obtains n phase error δ that encircles the territory correspondence n, n=1 ..., N;
(6) determine antenna opening diametric plane field distribution of amplitudes Q (ρ '), according to the phase error δ in each ring territory n,, calculate the electrical quantity of antenna by antenna far field Electric Field Distribution function;
(7), judge whether the antenna electrical quantity that calculates meets the demands, if meet the demands then the antenna structure design is qualified according to the requirement on electric performance of Antenna Design; Otherwise, revise parameter of structure design, and repeating step (1) is to step (6), until meeting the demands.
Described step (3), carry out according to the following procedure:
(3a) in new coordinate system, calculate the error of coordinate Δ r of practical distortion face and match distortion reflective side, Δ r=r (P 1)-r (P 0), in the formula, P 1Be the coordinate of antenna deformation face up-sampling point, P 0It is the coordinate of match distortion reflective side up-sampling point;
(3b) utilize practical distortion to face the error of coordinate Δ r of match distortion reflective side, according to the principle of least square and integration extreme value theorem, equationof structure group A β=H,
A is a coefficient,
A = Σ i = 1 num x i 2 2 f Σ x i y i 2 f - Σ x i - Σ x i y i Σ x i 2 Σ x i z i f Σ i = 1 num x i y i 2 f Σ y i 2 2 f - Σ y i - Σ y i 2 Σ x i y i Σ y i z i f Σ i = 1 nu m x i z i 2 f Σ y i z i 2 f - Σ z i - Σ y i z i Σ x i z i Σ z i 2 f Σ i = 1 num x i 2 f Σ y i 2 f - n - Σ y i Σ x i Σ z i f
β is the paraboloidal parameter of match, β=(Δ x Δ y Δ z φ xφ yΔ f) T,
H = ( Σ i = 1 num ( z i - z i ′ ) x i Σ i = 1 num ( z i - z i ′ ) y i Σ i = 1 num ( z i - z i ′ ) z i Σ i = 1 num ( z i - z i ′ ) ) T ;
(3c) find the solution above-mentioned equation group, obtain the paraboloidal parameter beta of match, be i.e. Δ x, Δ y, Δ z, φ x, φ yAnd Δ f, wherein Δ x, Δ y, Δ z are the parabolic displacement of summit in former coordinate system of match, φ x, φ yBe respectively the corner of the paraboloidal focal axis of match around former reference axis x, y, Δ f is the focal length variations amount;
(3d), obtain the equation of match distortion reflective side with in the paraboloidal parameter substitution of the match match reflecting surface formula:
z = ( x - Δx ) 2 + ( y - Δy ) 2 4 ( f + Δf ) + Δz + y φ x - x φ y
In the formula, f is the focal length of distortion front-reflection face.
Described step (4), carry out according to the following procedure:
(4a) the paraboloidal normal direction cosine of match is equaled the parabolic normal direction cosine of going up the relevant position of former design, obtain the node P (x after the distorted reflector according to the distortion posterior nodal point i+ Δ x i, y i+ Δ y i, z i+ Δ z i) normal direction cosine (l i, m i, n i):
l i = - x i 2 f ( f + z i ) , m i = - y i 2 f ( f + z i ) , n i = f f + z i ;
(4b) by distortion front nodal point P (x i, y i, z i) normal direction cosine, obtain through node P 1(x 0, y 0, z 0) normal equation:
x - x 0 l = y - y 0 m = z - z 0 n ;
(4c) find the solution normal equation, must be the normal direction straight line of actual node and the z coordinate of the parabolic intersection point of match, and the node P after the distortion 1(x 0, y 0, z 0) the intersection point P of the paraboloidal normal direction of corresponding match 0(x ' 0, y ' 0, z ' 0) coordinate;
(4d) by following formula, calculate the node P of distortion reflective side 1With respect to the parabolic node P of match 0Normal direction deviation ε i:
ϵ i = ( x 0 - x 0 ′ ) 2 + ( y 0 - y 0 ′ ) 2 + ( z 0 - z 0 ′ ) 2 ;
(4e) calculate new focus and focal axis, obtain new focal distance f;
(4f) according to the normal direction deviation ε of each node i, calculate the phase error δ that each node causes at the antenna opening diametric plane i:
δ i = 4 π λ ϵ i cos ξ i
Wherein: λ is the antenna operation wavelength, tan ζ ii'/2f, ρ i' be node P 1Respective radius at the bore face.
Described step (6), carry out according to the following procedure:
(6a) calculate antenna aperture field distribution of amplitudes Q (ρ) according to following formula:
Q ( ρ ) = B + C ( 1 - ρ 2 a 2 ) P
In the formula, B+C=1, P are the aperture field distributed constant, and a is the reflector antenna radius;
(6b) each piece territory phase error δ jBe used as a constant, as the phase error of ring domain center correspondence;
(6c) according to above-mentioned aperture field distribution of amplitudes and phase distribution parameters, by antenna far field electric field formula E = ∫ ∫ A Q ( ρ ′ ) e jδ ( ρ ′ ) e jk ρ → ′ · r ^ ds ′ , Calculate antenna far field Electric Field Distribution, ρ=ρ ' in the formula/a, ρ '=(a N-1+ a n)/2=(2n-1) a/2N, Be the distance vector of far field point of observation to initial point, k is a wave constant;
(6d) according to antenna far field Electric Field Distribution, calculate antenna far field direction of an electric field figure, from directional diagram, obtain the unit for electrical property parameters such as gain, minor level and lobe width of antenna.
The present invention compared with prior art has following advantage:
1. utilize distorted reflector posterior nodal point coordinate and the minimum difference that is out of shape the front nodal point coordinate, the paraboloidal equation of structure match simulates distortion reflective side, is easy to find normal direction of being out of shape node and the normal error that calculates this point by this match reflecting surface;
2. owing to be incorporated into the node normal error in the phase error of antenna opening diametric plane by optical path difference, antenna structure displacement field and electromagnetic field are closely connected, thereby avoided only using reflector precision to judge the deficiency of antenna performance, realize the electrical and mechanical comprehensive analysis of antenna;
3. by calculating antenna far field Electric Field Distribution, obtain the antenna electrical quantity, and with the electrical quantity and the antenna electric performance requirement contrast of calculating, predict by the electrical property in the structured design process, judge the reasonability of structural design, avoid carrying out by rule of thumb the antenna structure design, realized the Electromechanical Design of antenna.
Description of drawings
Fig. 1 is a reflector antenna electrical property prediction flow chart of the present invention;
Fig. 2 is a match distortion reflective side flow chart of the present invention;
Fig. 3 is that the present invention calculates bore face phase error flow chart;
Fig. 4 is that the present invention calculates reflector antenna electrical quantity flow chart;
Fig. 5 is the schematic diagram that reflector antenna is in any pose;
Fig. 6 is the schematic diagram that reflector antenna is in the state of looking up to heaven;
Fig. 7 is that the paraboloidal geometric parameter of match concerns schematic diagram;
Fig. 8 is the normal error schematic diagram of distortion reflective side antenna;
Fig. 9 is that reflector antenna bore face divides the ring schematic diagram;
Figure 10 is the used 7.3m antenna structure threedimensional model schematic diagram of emulation of the present invention.
Embodiment
Followingly the present invention is described in further detail with reference to accompanying drawing.
With reference to Fig. 1, concrete steps of the present invention are as follows:
Step 1 is set up the reflector antenna FEM (finite element) model.
According to given antenna aperture D, focal distance f, Reflector Panel, radiation beam, ring beam, centerbody basic parameter, determine the antenna structure FEM (finite element) model, obtain the theoretical coordinate P (x of reflecting surface sampling node i, y i, z i).Simultaneously, determine the coordinate h of reflecting surface summit in model coordinate systems, be generally Z to height.
Step 2 is analyzed FEM (finite element) model, the sampling node information after obtaining being out of shape.
Utilize finite element software, the antenna-reflected body structure is carried out finite element analysis under deadweight, wind lotus, vibration, the different operating modes with the ice and snow load of temperature load, obtain displacement P (the Δ x of reflecting surface sampling node i, Δ y i, Δ z i), and carry out Coordinate Conversion according to the coordinate system of structural model, and after obtaining being out of shape is initial point with the reflecting surface summit, and the antenna opening diametric plane is the XY plane, and focal axis is the new coordinate system of Z axle, and its process such as Fig. 5 and Fig. 6 institute are not.
Azimuth shown in Figure 5 and the angle of pitch are respectively the antenna of any pose of az and el, use the Coordinate Conversion formula: (x *y *z *1) R=(x y z 1) y(90 °-az) R x(el) R x(90 °) are converted to the state of looking up to heaven as shown in Figure 6.In the formula, (x *, y *, z *) for be in assigned address (az, antenna el) is a coordinate among O '-x ' y ' z ' at the state coordinate of looking up to heaven, (x, y is z) for being in assigned address (az, the coordinate of antenna el) in coordinate system O-xyz, R x(θ) and R yAll be the coordinate transformation matrix (φ), they are respectively:
R x ( θ ) = 1 0 0 0 0 cos θ sin θ 0 0 - sin θ cos θ 0 0 0 0 1 , R y ( φ ) = cos φ 0 - sin φ 0 0 1 0 0 sin φ 0 cos φ 0 0 0 0 1
Step 3, the match distortion reflective side.
With reference to Fig. 2, the concrete steps of distortion reflective side match are as follows:
1) in new coordinate system, calculates the error of coordinate Δ r of practical distortion face and match distortion reflective side, Δ r=r (P 1)-r (P 0), in the formula, P 1Be the coordinate of antenna deformation face up-sampling point, P 0Be the coordinate of match distortion reflective side up-sampling point, this practical distortion face and match reflection relation of plane are as shown in Figure 7;
2) the eyeball P of employing deformation surface 1With match face corresponding points P 0Axial error come the structural environment equation, according to the principle of least square and integration extreme value theorem, obtain following equation group:
A·β=H (1)
Wherein,
A = Σ i = 1 num x i 2 2 f Σ x i y i 2 f - Σ x i - Σ x i y i Σ x i 2 Σ x i z i f Σ i = 1 num x i y i 2 f Σ y i 2 2 f - Σ y i - Σ y i 2 Σ x i y i Σ y i z i f Σ i = 1 nu m x i z i 2 f Σ y i z i 2 f - Σ z i - Σ y i z i Σ x i z i Σ z i 2 f Σ i = 1 num x i 2 f Σ y i 2 f - n - Σ y i Σ x i Σ z i f
β=(ΔxΔyΔzφ xφ yΔf) T
H = ( Σ i = 1 num ( z i - z i ′ ) x i Σ i = 1 num ( z i - z i ′ ) y i Σ i = 1 num ( z i - z i ′ ) z i Σ i = 1 num ( z i - z i ′ ) ) T
Wherein: num is the sampling node sum, (x i, y i, z i') be node theoretical coordinate (x i, y i, z i) at the node coordinate of distortion reflective side correspondence, Δ x, Δ y, Δ z are the parabolic displacement of summit in former coordinate system of match, φ x, φ yBe respectively the corner of the paraboloidal focal axis of match around former reference axis x, y, Δ f is the focal length variations amount, and f is the focal length of distortion front-reflection face;
3) find the solution above-mentioned equation group (1), obtain paraboloidal 6 the geometric parameter β undetermined of match, i.e. Δ x, Δ y, Δ z, φ x, φ yAnd Δ f;
4) with in the paraboloidal parameter substitution of the match match reflecting surface formula, obtain the equation of match distortion reflective side:
z = ( x - Δx ) 2 + ( y - Δy ) 2 4 ( f + Δf ) + Δz + y φ x - x φ y . - - - ( 2 )
Step 4, the phase error of calculating bore face.
With reference to Fig. 3, the computational process of bore face phase error is as follows:
1) according to the distortion posterior nodal point the paraboloidal normal direction cosine of match is equaled the parabolic normal direction cosine of going up the relevant position of former design, obtain the node P (x after the distorted reflector i+ Δ x i, y i+ Δ y i, z i+ Δ z i) normal direction cosine be
l i = - x i 2 f ( f + z i ) , m i = - y i 2 f ( f + z i ) , n i = f f + z i - - - ( 3 )
Wherein: (l i, m i, n i) be the normal direction cosine of node P.
2) by the node P (x after the distortion i+ Δ x i, y i+ Δ y i, z i+ Δ z i) normal direction cosine, obtain through node P 1(x 0, y 0, z 0) normal equation:
x - x 0 l = y - y 0 m = z - z 0 n - - - ( 4 )
Point P (x P, y P, z P) be the position before the distortion of reflecting surface node, put P 1Be the physical location after the node P distortion, some P 0Be node P 1Along normal direction and the paraboloidal intersection point of match, as shown in Figure 8, Fig. 8 has described the geometry of the reflector antenna with surface error.
3) find the solution this equation (4), the node P after obtaining being out of shape 1(x 0, y 0, z 0) the intersection point P of the paraboloidal normal direction of corresponding match 0(x ' 0, y ' 0, z ' 0) coordinate is:
x 0 ′ = c 1 · z + c 2 y 0 ′ = c 3 · z + c 4 z 0 ′ = 1 - s 3 - s 4 - ( s 3 + s 4 - 1 ) 2 - 4 · s 2 s 5 - 4 · s 2 s 6 2 · s 2 - - - ( 5 )
In the formula:
c 1 = l n , c 2 = - l n z 0 + x 0 , c 3 = m n , c 4 = - m n z 0 + y 0 ,
s1=4(f+Δf), s 2 = c 1 2 + c 3 2 s 1 , s 3 = 2 [ c 1 ( c 2 - Δx ) + c 3 ( c 4 - Δy ) ] s 1 ,
s4=c3·φ x-c1·φ y s 5 = ( c 2 - Δx ) 2 + ( c 4 - Δy ) 2 s 1 , s6=Δz+c4·φ x-c2·φ y
4) the node P by distortion reflective side 1Calculating this point with respect to the paraboloidal normal direction deviation of match is:
ϵ i = ( x 0 - x 0 ′ ) 2 + ( y 0 - y 0 ′ ) 2 + ( z 0 - z 0 ′ ) 2 - - - ( 6 )
5) determine the position of focal axis
Theoretical paraboloidal coordinate is O-xyz, and the paraboloidal coordinate of match is O '-x ' y ' z ', and wherein O and O ' are respectively their summit, and Oz and O ' z ' are respectively their focal axis.
Theoretical paraboloidal equation is: x 2+ y 2=4fz
The paraboloidal equation of match is: x ' 2+ y ' 2=4 (f+ Δ f) z '
Make φ again zBe the corner of coordinate system O '-x ' y ' z ' around former reference axis Oz.Because of φ x, φ y, φ zBe trace, then can ignore their second order trace, so the Coordinate Conversion equation of two coordinate systems is
x 1 = ( x - Δx ) - z φ y + y φ z y 1 = ( y - Δy ) - x φ z + z φ x z 1 = ( z - Δz ) - y φ x + x φ y - - - ( 7 )
Through type (2) as can be known, the paraboloidal summit of match, the coordinate of initial point in former coordinate system O-xyz of promptly new coordinate system be (Δ x, Δ y, Δ z), newly the coordinate of the focal point F of focal axis ' in O-xyz is (x F ', y F ', z F '), the coordinate in new coordinate system O '-x ' y ' z ' is (0,0, f+ Δ f).Therefore, the coordinate of the initial point of former coordinate system O-xyz in new coordinate system O '-x ' y ' z ' is (Δ x ,-Δ y ,-Δ z), and its coordinate system is respectively-φ around the corner of new reference axis O ' x ', O ' y ', O ' z ' x,-φ yWith-φ z, then can get the coordinate of new focus in former coordinate system O-xyz and be
x F ′ = Δx + ( f + Δf ) · φ y y F ′ = Δy - ( f + Δf ) · φ x z F ′ = f + Δf + Δz . - - - ( 8 )
6) according to the normal direction deviation ε of each node i, calculate the phase error δ that each node causes at the bore face i
δ i = 4 π λ ϵ i cos ξ i - - - ( 9 )
Wherein: λ is the antenna operation wavelength, tan ζ ii'/2f, ρ i' be node P 1In the respective radius of bore face, f is new focal coordinates;
Step 5 is calculated the phase error in each ring territory.
As shown in Figure 9, the antenna actinal surface be divided into N the ring territory, n (n=1 ..., N) on the individual ring territory K is arranged nIndividual node.Suppose to have on the antenna reflective face num impact point, then num and K nSatisfy following relation:
num = Σ n = 1 N K n - - - ( 10 )
Calculate K on n the ring territory according to formula (9) nThe bore face phase error of individual node correspondence obtains n phase error δ that encircles the territory correspondence nFor:
δ n = 1 K n Σ j = 1 K n δ j 2 - - - ( 11 )
Step 6 is calculated the antenna electrical quantity.
With reference to Fig. 4, antenna computation of electrical step is as follows:
1) calculating aperture field distribution of amplitudes Q (ρ) is
Q ( ρ ) = B + C ( 1 - ρ 2 a 2 ) P - - - ( 12 )
Wherein, B+C=1 selects suitable B, just can produce different aperture field edge illumination level, and aperture field distributed constant P is used for controlling the shape that aperture field distributes;
2) each piece territory phase error δ jBe used as a constant, as the phase error of ring domain center correspondence;
3) according to above-mentioned aperture field distribution of amplitudes and phase distribution parameters, the transformation for mula that distributes by bore calculates Electric Field Distribution:
E = ∫ ∫ A Q ( ρ ′ ) e jδ ( ρ ′ ) e jk ρ → ′ · r ^ ds ′ - - - ( 13 )
Formula (13) is carried out discretization handle, to calculate the far field Electric Field Distribution of antenna:
E = 2 π a 2 Σ n = 1 N E n , n - 1 e j δ n - - - ( 14 )
Wherein
E n,n-1=E n-E n-1 E n = ∫ 0 n / N Q ( ρ ) J 0 ( uρ ) ρdρ , ρ '=(a N-1+ a n)/2=(2n-1) a/2N, ρ=ρ '/a, u=kasin θ, a are the actinal surface radius, and k is the free space wave constant, and θ is the azimuth of far field impact point, J 0(*) be first kind Bessel function.
4) according to antenna far field Electric Field Distribution, calculate antenna far field direction of an electric field figure, from directional diagram, obtain the unit for electrical property parameters such as gain, minor level and lobe width of antenna.
Step 7 according to the requirement on electric performance of Antenna Design, judges whether the antenna electrical quantity calculate meets the demands, if meet the demands then the antenna structure design is qualified; Otherwise, revise parameter of structure design, and repeating step (1) meets the demands until the antenna electrical quantity to step (6).
Advantage of the present invention can further specify by following emulation experiment:
1. simulated conditions
Antenna electromechanical properties comprehensive analysis method of the present invention is compiled to " reflector antenna electromechanical properties analysis-by-synthesis software ", on 7.3m bore boat-carrying parabolic antenna, carried out antenna reflector structural design and electrical property prognostic experiment.
This antenna reflective face is real template, and back of the body frame belongs to the truss class.The piece panel adopts the design of rigid panel structural shape, is riveted and is formed by longitudinal rib, the hoop muscle of stretching covering and drawn moulding.The material of monolith surface board member all adopts duralumin metal plate LY12-M, is considered as shell unit when finite element analysis, gets Shell 63 in ANSYS software.Entire antenna radially is divided into three circles, is evenly distributed with 16 radiation beams along the centerbody circumference, 48 annular girders.Working frequency range is Ku frequency range and C frequency range, this Ku band gain requirement 〉=57.2dB, this C band gain requirement 〉=50.2dB.The secondary lobe envelope satisfies the CCIR.580-2 requirement, and work wind speed is 20m/s, and the existence wind speed is 55m/s, and operating ambient temperature range is-45 ℃ to 60 ℃, and the surface normal precision index is 0.5mm, entire antenna reflector structure weight≤1.5 tons.Utilize antenna structure displacement field of the present invention and electromagnetic field couples model, 7.3m ship-board antenna reflector structure weight is reduced to 1.2 tons from 1.5 tons, and antenna gain is brought up to 58.76dB.Here getting the antenna irradiation taper is-10dB, and the aperture field distributed constant equals 1, and operating frequency is 14GHz.
2. simulation result
The antenna structure model that utilizes above condition emulation to be set up, as shown in figure 10.Emulation is under deadweight and the effect of wind lotus, and during the different operating elevation angle, the antenna electric performance parameter of eight kinds of working conditions is as shown in table 1 altogether.
The antenna electric performance major parameter of table 1 emulation under the different operating situation
Data can be found out from table 1, and the antenna of setting up based on dynamo-electric integrated approach of the present invention is that the antenna gain under 90 ° the deadweight operating mode is 58.761dB at the angle of pitch, satisfy 〉=designing requirement of 57.2dB.
By the experimental result of this software, prove and adopt method of the present invention to can be used for carrying out reflector antenna structural design and electrical property prediction.

Claims (4)

1. method for predicting antenna electric performance based on the match distortion reflective side comprises following process:
(1) according to Reflector Panel, back of the body frame, the centerbody structural parameters of antenna, determines the antenna structure FEM (finite element) model, obtain the theoretical coordinate P (x of reflecting surface sampling node i, y i, z i);
(2) utilize finite element software, the antenna reflective face structure is carried out finite element analysis, obtain displacement P (the Δ x of distorted reflector post-sampling node i, Δ y i, Δ z i), and carry out Coordinate Conversion by the structural model coordinate system, obtain the new coordinate system of antenna reflective face;
(3) in new coordinate system, utilize distorted reflector posterior nodal point coordinate and the minimum difference that is out of shape the front nodal point coordinate, construct the paraboloidal equation group of match, obtain the equation of match distortion reflective side;
(4) utilize match distortion reflective side equation, obtain distorted reflector posterior nodal point P (x i+ Δ x i, y i+ Δ y i, z i+ Δ z i) normal direction cosine and normal direction deviation ε i, calculate the phase error δ i that each node causes at the bore face;
(5) the antenna opening diametric plane is divided into N ring territory, calculates n the K on the ring territory nThe bore face phase error of individual node correspondence obtains n phase error δ that encircles the territory correspondence n, n=1 ..., N;
(6) determine antenna opening diametric plane field distribution of amplitudes Q (ρ '), according to the phase error δ in each ring territory n,, calculate the electrical quantity of antenna by antenna far field Electric Field Distribution function;
(7), judge whether the antenna electrical quantity that calculates meets the demands, if meet the demands then the antenna structure design is qualified according to the requirement on electric performance of Antenna Design; Otherwise, revise parameter of structure design, and repeating step (1) is to step (6), until meeting the demands.
2, method for predicting antenna electric performance according to claim 1 is characterized in that step (3) carries out according to the following procedure:
(3a) in new coordinate system, calculate the error of coordinate Δ r of practical distortion face and match distortion reflective side, Δ r=r (P 1)-r (P 0), in the formula, P 1Be the coordinate of antenna deformation face up-sampling point, P 0It is the coordinate of match distortion reflective side up-sampling point;
(3b) utilize practical distortion to face the error of coordinate Δ r of match distortion reflective side, according to the principle of least square and integration extreme value theorem, equationof structure group A β=H,
A is a coefficient,
A = Σ i = 1 num x i 2 2 f Σ x i y i 2 f - Σ x i - Σ x i y i Σ x i 2 Σ x i z i f Σ i = 1 num x i y i 2 f Σ y i 2 2 f - Σ y i - Σ y i 2 Σ x i y i Σ y i z i f Σ i = 1 nu m x i z i 2 f Σ y i z i 2 f - Σ z i - Σ y i z i Σ x i z i Σ z i 2 f Σ i = 1 num x i 2 f Σ y i 2 f - n - Σ y i Σ x i Σ z i f
β is the paraboloidal parameter of match, β=(Δ x Δ y Δ z φ xφ yΔ f) T
H = ( Σ i = 1 num ( z i - z i ′ ) x i Σ i = 1 num ( z i - z i ′ ) y i Σ i = 1 num ( z i - z i ′ ) z i Σ i = 1 num ( z i - z i ′ ) ) T ;
(3c) find the solution above-mentioned equation group, obtain the paraboloidal parameter beta of match, be i.e. Δ x, Δ y, Δ z, φ x, φ yAnd Δ f, wherein Δ x, Δ y, Δ z are the parabolic displacement of summit in former coordinate system of match, φ x, φ yBe respectively the corner of the paraboloidal focal axis of match around former reference axis x, y, Δ f is the focal length variations amount;
(3d), obtain the equation of match distortion reflective side with in the paraboloidal parameter substitution of the match match reflecting surface formula:
z = ( x - Δx ) 2 + ( y - Δy ) 2 4 ( f + Δf ) + Δz + y φ x - x φ y
In the formula, f is the focal length of distortion front-reflection face.
3, method for predicting antenna electric performance according to claim 1 is characterized in that step (4), carries out according to the following procedure:
(4a) the paraboloidal normal direction cosine of match is equaled the parabolic normal direction cosine of going up the relevant position of former design, obtain the node P (x after the distorted reflector according to the distortion posterior nodal point i+ Δ x i, y i+ Δ y i, z i+ Δ z i) normal direction cosine (l i, m i, n i):
l i = - x i 2 f ( f + z i ) , m i = - y i 2 f ( f + z i ) , n i = f f + z i ;
(4b) by distortion front nodal point P (x i, y i, z i) normal direction cosine, obtain through node P 1(x 0, y 0, z 0) normal equation:
x - x 0 l = y - y 0 m = z - z 0 n ;
(4c) find the solution normal equation, must be the normal direction straight line of actual node and the z coordinate of the parabolic intersection point of match, and the node P after the distortion 1(x 0, y 0, z 0) the intersection point P of the paraboloidal normal direction of corresponding match 0(x ' 0, y ' 0, z ' 0) coordinate;
(4d) by following formula, calculate the node P of distortion reflective side 1With respect to the parabolic node P of match 0Normal direction deviation ε i:
ϵ i = ( x 0 - x 0 ′ ) 2 + ( y 0 - y 0 ′ ) 2 + ( z 0 - z 0 ′ ) 2 ;
(4e) calculate new focus and focal axis, obtain new focal distance f;
(4f) according to the normal direction deviation ε of each node i, calculate the phase error δ that each node causes at the antenna opening diametric plane i:
δ i = 4 π λ ϵ i cos ξ i
Wherein: λ is the antenna operation wavelength, tan ζ ii'/2f, ρ i' be node P 1Respective radius at the bore face.
4, method for predicting antenna electric performance according to claim 1 is characterized in that step (6), carries out according to the following procedure:
(6a) calculate antenna aperture field distribution of amplitudes Q (ρ) according to following formula:
Q ( ρ ) = B + C ( 1 - ρ 2 a 2 ) P
In the formula, B+C=1, P are the aperture field distributed constant;
(6b) each piece territory phase error δ jBe used as a constant, as the phase error of ring domain center correspondence;
(6c) according to above-mentioned aperture field distribution of amplitudes and phase distribution parameters, by antenna far field electric field formula E = ∫ ∫ A Q ( ρ ′ ) e jδ ( ρ ′ ) e jk ρ → ′ · r ^ ds ′ , Calculate antenna far field Electric Field Distribution, ρ=ρ ' in the formula/a, ρ '=(a N-1+ a n)/2=(2n-1) a/2N, Be the distance vector of far field point of observation to initial point, k is a wave constant;
(6d) according to antenna far field Electric Field Distribution, calculate antenna far field direction of an electric field figure, from directional diagram, obtain the unit for electrical property parameters such as gain, minor level and lobe width of antenna.
CN200810018106A 2008-04-30 2008-04-30 Method for predicting antenna electric performance based on simulated distortion reflective side Expired - Fee Related CN101267062B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN200810018106A CN101267062B (en) 2008-04-30 2008-04-30 Method for predicting antenna electric performance based on simulated distortion reflective side

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN200810018106A CN101267062B (en) 2008-04-30 2008-04-30 Method for predicting antenna electric performance based on simulated distortion reflective side

Publications (2)

Publication Number Publication Date
CN101267062A true CN101267062A (en) 2008-09-17
CN101267062B CN101267062B (en) 2012-10-24

Family

ID=39989287

Family Applications (1)

Application Number Title Priority Date Filing Date
CN200810018106A Expired - Fee Related CN101267062B (en) 2008-04-30 2008-04-30 Method for predicting antenna electric performance based on simulated distortion reflective side

Country Status (1)

Country Link
CN (1) CN101267062B (en)

Cited By (23)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101989682A (en) * 2010-07-21 2011-03-23 中国科学院国家天文台南京天文光学技术研究所 Antenna reflection surface fitting method
CN102073754A (en) * 2010-09-20 2011-05-25 西安电子科技大学 Comprehensive electromechanical analysis method of reflector antenna based on error factor
CN102156764A (en) * 2010-12-27 2011-08-17 南京理工大学 Multi-resolution precondition method for analyzing aerial radiation and electromagnetic scattering
CN102708257A (en) * 2012-05-23 2012-10-03 西安电子科技大学 Method for fast determining tolerance of active phased array antenna structure
CN102788920A (en) * 2012-07-30 2012-11-21 西安电子科技大学 Electrical property prediction method of offset reflector antenna based on electromechanical coupling model
CN102968532A (en) * 2012-11-16 2013-03-13 西安电子科技大学 Electromechanical integration design method of 65m-aperture large-size reflector antenna structure
CN103020328A (en) * 2011-09-20 2013-04-03 深圳光启高等理工研究院 Optimum design method for small antenna
CN103473421A (en) * 2013-09-17 2013-12-25 西安电子科技大学 Foundation antenna temperature compensation method
CN104063426A (en) * 2014-06-06 2014-09-24 西安电子科技大学 Radiation-and-scattering-oriented quick determination method for active phased-array antenna structure tolerance
CN104111053A (en) * 2014-07-08 2014-10-22 哈尔滨工业大学 Method of analyzing reflection plane precision of parabolic antenna based on coefficient matrix QR decomposition calculation method
CN104157986A (en) * 2013-05-29 2014-11-19 广东通宇通讯股份有限公司 Microwave antenna suitable for quite-low-interference communication system, and optimization method thereof
WO2014190514A1 (en) * 2013-05-29 2014-12-04 Tongyu Communication Inc. Microwave antennas for extremely low interference communications systems
CN104537188A (en) * 2015-01-16 2015-04-22 西安电子科技大学 Reflector antenna panel functionality design method based on contour inversion technology
CN104682011A (en) * 2015-03-08 2015-06-03 西安电子科技大学 Ground reflector in cable network structure
CN104925563A (en) * 2015-05-20 2015-09-23 西安电子科技大学 Cutting method of spiral deployable antenna reflective surface and unfolding device thereof
CN105337045A (en) * 2015-10-12 2016-02-17 中国电子科技集团公司第五十四研究所 Optimal selection method of distribution of reflector antenna aperture field
CN105576385A (en) * 2016-02-02 2016-05-11 西安电子科技大学 Gain and direction oriented large-scale deformed paraboloid antenna panel fit rotation adjusting method
CN105718697A (en) * 2016-02-02 2016-06-29 西安电子科技大学 Antenna pointing oriented large deformed parabolic antenna panel movement fitting adjustment method
CN105740554A (en) * 2016-02-02 2016-07-06 西安电子科技大学 Gain-oriented large-scale deformed parabolic antenna panel precision adjusting method
CN107239602A (en) * 2017-05-22 2017-10-10 中国电子科技集团公司第四十研究所 A kind of probe antenna model quick calculation method based on curve matching
US9835664B2 (en) 2013-05-29 2017-12-05 Tongyu Communication Inc. Microwave antennas for extremely low interference communications systems
CN108256203A (en) * 2018-01-12 2018-07-06 西安电子科技大学 A kind of minor face deformation position posture quick calculation method based on the strain of antenna support leg
CN110375705A (en) * 2019-08-20 2019-10-25 大连理工大学 Antenna reflector and its type face deformation measurement method and measurement of comparison method

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN200950469Y (en) * 2006-05-25 2007-09-19 深圳市华达玻璃钢通信制品有限公司 Sheet molding plastic antenna reflector

Cited By (40)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101989682A (en) * 2010-07-21 2011-03-23 中国科学院国家天文台南京天文光学技术研究所 Antenna reflection surface fitting method
CN102073754A (en) * 2010-09-20 2011-05-25 西安电子科技大学 Comprehensive electromechanical analysis method of reflector antenna based on error factor
CN102073754B (en) * 2010-09-20 2012-06-27 西安电子科技大学 Comprehensive electromechanical analysis method of reflector antenna based on error factor
CN102156764A (en) * 2010-12-27 2011-08-17 南京理工大学 Multi-resolution precondition method for analyzing aerial radiation and electromagnetic scattering
CN102156764B (en) * 2010-12-27 2013-06-26 南京理工大学 Multi-resolution precondition method for analyzing aerial radiation and electromagnetic scattering
CN103020328B (en) * 2011-09-20 2015-02-25 深圳光启创新技术有限公司 Optimum design method for small antenna
CN103020328A (en) * 2011-09-20 2013-04-03 深圳光启高等理工研究院 Optimum design method for small antenna
CN102708257A (en) * 2012-05-23 2012-10-03 西安电子科技大学 Method for fast determining tolerance of active phased array antenna structure
CN102708257B (en) * 2012-05-23 2014-04-02 西安电子科技大学 Method for fast determining tolerance of active phased array antenna structure
CN102788920A (en) * 2012-07-30 2012-11-21 西安电子科技大学 Electrical property prediction method of offset reflector antenna based on electromechanical coupling model
CN102788920B (en) * 2012-07-30 2015-05-20 西安电子科技大学 Electrical property prediction method of offset reflector antenna based on electromechanical coupling model
CN102968532A (en) * 2012-11-16 2013-03-13 西安电子科技大学 Electromechanical integration design method of 65m-aperture large-size reflector antenna structure
CN102968532B (en) * 2012-11-16 2015-09-23 西安电子科技大学 The dynamo-electric integral design method of the large-scale reflector antenna structure of 65m bore
CN104157986A (en) * 2013-05-29 2014-11-19 广东通宇通讯股份有限公司 Microwave antenna suitable for quite-low-interference communication system, and optimization method thereof
WO2014190514A1 (en) * 2013-05-29 2014-12-04 Tongyu Communication Inc. Microwave antennas for extremely low interference communications systems
CN104157986B (en) * 2013-05-29 2017-03-22 广东通宇通讯股份有限公司 Microwave antenna suitable for quite-low-interference communication system, and optimization method thereof
US9835664B2 (en) 2013-05-29 2017-12-05 Tongyu Communication Inc. Microwave antennas for extremely low interference communications systems
CN103473421B (en) * 2013-09-17 2016-06-29 西安电子科技大学 foundation antenna temperature compensation method
CN103473421A (en) * 2013-09-17 2013-12-25 西安电子科技大学 Foundation antenna temperature compensation method
CN104063426A (en) * 2014-06-06 2014-09-24 西安电子科技大学 Radiation-and-scattering-oriented quick determination method for active phased-array antenna structure tolerance
CN104063426B (en) * 2014-06-06 2017-02-15 西安电子科技大学 Radiation-and-scattering-oriented quick determination method for active phased-array antenna structure tolerance
CN104111053A (en) * 2014-07-08 2014-10-22 哈尔滨工业大学 Method of analyzing reflection plane precision of parabolic antenna based on coefficient matrix QR decomposition calculation method
CN104537188A (en) * 2015-01-16 2015-04-22 西安电子科技大学 Reflector antenna panel functionality design method based on contour inversion technology
CN104682011A (en) * 2015-03-08 2015-06-03 西安电子科技大学 Ground reflector in cable network structure
CN104682011B (en) * 2015-03-08 2017-06-13 西安电子科技大学 Ground cable net structure reflector
CN104925563A (en) * 2015-05-20 2015-09-23 西安电子科技大学 Cutting method of spiral deployable antenna reflective surface and unfolding device thereof
CN104925563B (en) * 2015-05-20 2016-10-19 西安电子科技大学 The cutting method of a kind of spiral deployable antenna reflecting surface and expanding unit thereof
CN105337045A (en) * 2015-10-12 2016-02-17 中国电子科技集团公司第五十四研究所 Optimal selection method of distribution of reflector antenna aperture field
CN105337045B (en) * 2015-10-12 2018-01-30 中国电子科技集团公司第五十四研究所 A kind of method for optimizing of reflector antenna Aperture field distribution
CN105740554B (en) * 2016-02-02 2018-10-09 西安电子科技大学 A kind of large-scale deformation parabola antenna panel precision regulation method towards gain
CN105718697A (en) * 2016-02-02 2016-06-29 西安电子科技大学 Antenna pointing oriented large deformed parabolic antenna panel movement fitting adjustment method
CN105576385A (en) * 2016-02-02 2016-05-11 西安电子科技大学 Gain and direction oriented large-scale deformed paraboloid antenna panel fit rotation adjusting method
CN105718697B (en) * 2016-02-02 2019-02-19 西安电子科技大学 The large-scale mobile fitting method of adjustment of deformation parabola antenna panel is directed toward towards antenna
CN105576385B (en) * 2016-02-02 2018-05-25 西安电子科技大学 Towards gain method of adjustment is rotated with being directed toward large-scale deformation parabola antenna panel and coincideing
CN105740554A (en) * 2016-02-02 2016-07-06 西安电子科技大学 Gain-oriented large-scale deformed parabolic antenna panel precision adjusting method
CN107239602A (en) * 2017-05-22 2017-10-10 中国电子科技集团公司第四十研究所 A kind of probe antenna model quick calculation method based on curve matching
CN107239602B (en) * 2017-05-22 2020-07-03 中国电子科技集团公司第四十一研究所 Probe antenna model fast calculation method based on curve fitting
CN108256203A (en) * 2018-01-12 2018-07-06 西安电子科技大学 A kind of minor face deformation position posture quick calculation method based on the strain of antenna support leg
CN108256203B (en) * 2018-01-12 2020-06-23 西安电子科技大学 Rapid calculation method for secondary surface deformation position posture based on antenna supporting leg strain
CN110375705A (en) * 2019-08-20 2019-10-25 大连理工大学 Antenna reflector and its type face deformation measurement method and measurement of comparison method

Also Published As

Publication number Publication date
CN101267062B (en) 2012-10-24

Similar Documents

Publication Publication Date Title
CN101267062B (en) Method for predicting antenna electric performance based on simulated distortion reflective side
CN101308177B (en) Initiative reflecting plane antenna electrical behavior prediction method
CN101344564B (en) Active phase array antenna electrical property prediction method based on mechanical, electric and thermal three-field coupling
CN105977649A (en) Shaped surface-oriented quick determination method for adjustment quantity of active panel of large parabolic antenna
CN102968532B (en) The dynamo-electric integral design method of the large-scale reflector antenna structure of 65m bore
CN105740502A (en) Electric property prediction method of cylinder conformal array antenna under load distortion
CN102073754B (en) Comprehensive electromechanical analysis method of reflector antenna based on error factor
CN102590656B (en) Antenna cover electric property forecasting method based on distant field
CN104112051A (en) Element shape function based reflector antenna electromechanical integration design method
CN106981728A (en) A kind of bare cloth line array integrated approach of directional diagram reconstructable
CN105718697A (en) Antenna pointing oriented large deformed parabolic antenna panel movement fitting adjustment method
CN105302962A (en) Electromechanical integration optimization design method of reflector antenna on the basis of structure-electromagnetic hybrid unit
CN101877022B (en) Carrier antenna structure meshing method based on method of moment (MOM) for triangular vector basis function
CN106021766B (en) The high pointing accuracy design method of circular polarisation net-shape antenna based on electromechanical integrated optimization
CN105206941A (en) Electromechanical coupling-based direction adjustment method of large-scale forming double-reflection surface antenna
CN102788920B (en) Electrical property prediction method of offset reflector antenna based on electromechanical coupling model
CN106991210B (en) A kind of Shaped-beam reflector antenna electrical behavior prediction method based on electromechanical Coupling Model
CN101252226B (en) Method for positioning reflecting plane aerial feed source
CN105576385A (en) Gain and direction oriented large-scale deformed paraboloid antenna panel fit rotation adjusting method
Wang et al. Optimal surface adjustment by the error-transformation matrix for a segmented-reflector antenna
CN106599504B (en) Space mesh antenna dynamic load analysis method based on electromechanical coupling model
CN105740554A (en) Gain-oriented large-scale deformed parabolic antenna panel precision adjusting method
CN103490172B (en) Deformation reflector antenna phase center correction method based on far field
Tian et al. Geometry modeling of truss structure for space deployable truss antenna with multi-module
CN104063587B (en) The method that panel mismachining tolerance is affected on electrical property is calculated based on block form

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: 20121024

Termination date: 20150430

EXPY Termination of patent right or utility model