CN104063587B - The method that panel mismachining tolerance is affected on electrical property is calculated based on block form - Google Patents

The method that panel mismachining tolerance is affected on electrical property is calculated based on block form Download PDF

Info

Publication number
CN104063587B
CN104063587B CN201410257555.4A CN201410257555A CN104063587B CN 104063587 B CN104063587 B CN 104063587B CN 201410257555 A CN201410257555 A CN 201410257555A CN 104063587 B CN104063587 B CN 104063587B
Authority
CN
China
Prior art keywords
sigma
panel
error
rho
delta
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
CN201410257555.4A
Other languages
Chinese (zh)
Other versions
CN104063587A (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 University
Original Assignee
Xidian University
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 University filed Critical Xidian University
Priority to CN201410257555.4A priority Critical patent/CN104063587B/en
Publication of CN104063587A publication Critical patent/CN104063587A/en
Application granted granted Critical
Publication of CN104063587B publication Critical patent/CN104063587B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Aerials With Secondary Devices (AREA)
  • Variable-Direction Aerials And Aerial Arrays (AREA)

Abstract

本发明涉及基于分块形式计算面板加工误差对电性能影响的方法,用于指导面板表面的加工和天线电性能仿真分析与评价,本发明把受表面误差影响的远场表示为理想场与误差场之和,在求天线平均功率时可将不同面板之间的误差功率简化为零处理,减少了计算量;根据天线的实际分块形式计算表面加工误差对电性能的影响,更接近于工程实际;利用精确的四重数值积分求解单块面板的平均误差功率,避免了传统积分时间隔矢量Δ无法完整表示的缺陷;通过仿真分析得出表面误差均方根的上限值,避免了凭经验提出过高的精度要求,降低了加工难度同时又满足电性能指标。

The invention relates to a method for calculating the influence of panel processing errors on electrical performance based on the block form, which is used to guide the processing of the panel surface and the simulation analysis and evaluation of antenna electrical performance. The invention expresses the far field affected by the surface error as an ideal field and error The sum of fields can simplify the error power between different panels to zero when calculating the average power of the antenna, which reduces the amount of calculation; calculate the influence of surface processing errors on electrical performance according to the actual block form of the antenna, which is closer to engineering Practical; use accurate quadruple numerical integration to solve the average error power of a single panel, avoiding the defect that the traditional integral time interval vector Δ cannot be fully represented; through simulation analysis, the upper limit of the root mean square of the surface error is obtained, avoiding the need to rely on Experience puts forward high precision requirements, which reduces the difficulty of processing and meets the electrical performance indicators at the same time.

Description

基于分块形式计算面板加工误差对电性能影响的方法Method of Calculating the Influence of Panel Processing Errors on Electrical Properties Based on Block Form

技术领域technical field

本发明属于天线技术领域,具体是基于分块形式计算面板加工误差对电性能影响的方法,用于指导面板表面的加工和天线电性能仿真分析与评价。The invention belongs to the technical field of antennas, and specifically relates to a method for calculating the influence of panel processing errors on electrical performance based on block form, which is used to guide panel surface processing and antenna electrical performance simulation analysis and evaluation.

背景技术Background technique

反射面天线作为最常用的微波毫米波高增益天线,广泛应用于雷达测控、卫星通信、射电天文等领域,将来还有可能作为微波武器以及能量传输的工具。由于其工作在较高频段,电磁性能易受到各种误差的影响。反射面天线表面误差主要有两种:系统误差和随机误差。系统误差主要包括天线的自重变形、太阳不均匀照射引起的热变形以及其它环境载荷引起的变形;随机误差则来自于面板表面的加工误差以及拼接时的安装不确定性,对天线性能有比较明显的影响。长期以来,Ruze公式在某种程度上有效地指导了天线的结构设计,其本人也被公认为系统研究微小相位误差对反射面天线增益损失以及副瓣电平影响的第一人。但是Ruze理论假设整个口径面的相位误差呈相关半径为c的高斯帽状分布,并没有考虑由面板拼接而成的反射面中的高斯帽跨界计算问题。其之后的学者尽管进一步完善了Ruze公式,但同样没有按照面板的拼接形式来计算表面制造误差对天线性能的影响。As the most commonly used microwave and millimeter wave high-gain antenna, the reflector antenna is widely used in radar measurement and control, satellite communication, radio astronomy and other fields. It may also be used as a microwave weapon and a tool for energy transmission in the future. Since it works in a higher frequency band, the electromagnetic performance is susceptible to various errors. There are two main types of surface errors for reflector antennas: systematic errors and random errors. The system error mainly includes the deformation of the antenna due to its own weight, the thermal deformation caused by uneven sunlight, and the deformation caused by other environmental loads; the random error comes from the processing error of the panel surface and the installation uncertainty during splicing, which has a significant impact on the performance of the antenna. Impact. For a long time, the Ruze formula has effectively guided the structural design of the antenna to a certain extent, and he himself is recognized as the first person to systematically study the influence of small phase errors on the gain loss and sidelobe level of reflector antennas. However, the Ruze theory assumes that the phase error of the entire aperture surface is distributed in a Gaussian hat shape with a relative radius c, and does not consider the Gaussian hat cross-border calculation problem in the reflective surface formed by splicing panels. Although subsequent scholars further improved the Ruze formula, they also did not calculate the impact of surface manufacturing errors on antenna performance according to the panel splicing form.

随着反射面向着大口径、高频段、高增益的方向发展,仍用Ruze公式作为指导将会对表面精度提出极为苛刻的要求,使结构工程人员在背架保型、面板加工及装配方面将面临巨大挑战。然而实际上也许并不需要达到如此高的精度就可以满足电性能的要求,因此,更为准确而详细的性能预测显得尤为重要。With the development of reflectors in the direction of large aperture, high frequency band, and high gain, still using the Ruze formula as a guide will put forward extremely stringent requirements on surface accuracy, so that structural engineers will be able to maintain the shape of the back frame, panel processing and assembly. Facing huge challenges. However, in fact, it may not be necessary to achieve such a high precision to meet the electrical performance requirements. Therefore, more accurate and detailed performance prediction is particularly important.

发明内容Contents of the invention

本发明的目的是提出一种基于分块形式计算面板加工误差对电性能影响的方法,以便通过准确的仿真分析给出面板表面加工精度的合理要求,从而指导面板加工。The purpose of the present invention is to propose a method for calculating the influence of panel processing errors on electrical properties based on block form, so as to provide reasonable requirements for panel surface processing accuracy through accurate simulation analysis, thereby guiding panel processing.

实现本发明目的的技术方案是,基于分块形式计算面板加工误差对电性能影响的方法,其特征是:至少包括如下步骤:The technical solution for realizing the purpose of the present invention is a method for calculating the influence of panel processing errors on electrical performance based on block form, which is characterized in that it includes at least the following steps:

(1)把单块面板远场积分表达式中的相位影响项进行一阶泰勒展开,将远场表示为理想场与误差场之和;(1) The phase influence term in the far-field integral expression of a single panel is first-order Taylor expanded, and the far field is expressed as the sum of the ideal field and the error field;

(2)根据加工工艺确定每块面板的随机误差相关半径,然后得到同块面板不同位置处的误差相关函数,最后求得每块面板任意两点法向随机误差γn,i1)与γn,i2)乘积的均值<γn,i1n,i2)>;(2) Determine the random error correlation radius of each panel according to the processing technology, then obtain the error correlation function at different positions of the same panel, and finally obtain the normal random error γ n,i1 ) of any two points on each panel The mean value of the product with γ n,i2 )<γ n,i1n,i2 )>;

(3)利用法向误差与相位误差的关系求得每块面板任意两点相位误差δn,i1)与δn,i2)乘积的均值<δn,i1n,i2)>;(3) Use the relationship between the normal error and the phase error to obtain the mean value of the product of the phase error δ n,i1 ) and δ n,i2 ) at any two points on each panel<δ n,i1n,i2 )>;

(4)将面板在远场的误差场强δEn,i与其共轭相乘并求平均,推导得出单块面板在远场的平均误差功率表达式 (4) Conjugate the error field strength δE n,i of the panel in the far field to its conjugate Multiply and average to derive the average error power expression of a single panel in the far field

(5)确定每块面板的积分区间参数,包括径向参数,即:单块面板的内、外径:a11,a12,a21,a22;周向参数,即:单块面板两侧的角度坐标:α11,α12,α21,α22(5) Determine the integral interval parameters of each panel, including radial parameters, namely: inner and outer diameters of a single panel: a 11 , a 12 , a 21 , a 22 ; circumferential parameters, namely: two Angular coordinates of sides: α 11 , α 12 , α 21 , α 22 ;

(6)将步骤(3)得出的<δn,i1n,i2)>代入平均误差功率表达式中,并将被积函数用积分变量表达,然后通过四重数值积分得到单块面板平均误差功率;(6) Substitute the <δ n,i1n,i2 )> obtained in step (3) into the average error power expression, and express the integrand as an integral variable, and then pass four The average error power of a single panel is obtained by multiple numerical integration;

(7)将所有面板在远场的平均误差功率叠加,所得结果与天线在理想情况下的功率求和,最终得到含有表面加工误差的整个天线在远场的平均功率方向图。(7) The average error power of all panels in the far field is superimposed, and the result is summed with the power of the antenna under ideal conditions, and finally the average power pattern of the entire antenna in the far field including surface processing errors is obtained.

步骤(1)所述的将每块面板的远场表示为理想场与误差场之和,按如下过程进行:The far field of each panel described in step (1) is expressed as the sum of the ideal field and the error field, as follows:

(1a)由面板拼接而成的反射面共分为N环,第n环(n=1,2,…,N)又可分为Kn个子块;每块面板的表面都会存在制造误差,对第(n,i)块面板而言,设其表面某点处的随机误差引起的相位误差为δ(ρ),则其远场方向图可表示为(1a) The reflective surface spliced by panels is divided into N rings, and the nth ring (n=1,2,...,N) can be divided into K n sub-blocks; there will be manufacturing errors on the surface of each panel, For the (n,i)th panel, assuming that the phase error caused by a random error at a certain point on its surface is δ(ρ), its far-field pattern can be expressed as

式中,为天线口径场幅度分布函数,其中B+C=1,通过改变P值可控制口径场幅度分布函数的形状,通常1≤P≤2,ρ为口径面上天线某点的半径,a为天线的半径。k为电磁波传播波数,为目标方向的单位矢量;In the formula, is the antenna aperture field amplitude distribution function, where B+C=1, the shape of the aperture field amplitude distribution function can be controlled by changing the P value, usually 1≤P≤2, ρ is the radius of a certain point of the antenna on the aperture surface, and a is the antenna of the radius. k is the wave number of electromagnetic wave propagation, is the unit vector of the target direction;

(1b)对上式中的ejδ(ρ)作一阶泰勒展开,有ejδ(ρ)≈1+jδ(ρ),然后得到(1b) Perform first-order Taylor expansion on e jδ(ρ) in the above formula, e jδ(ρ) ≈1 +jδ(ρ), and then get

式中,为理想面板在远场的场强方向图,δEn,i为表面随机制造误差引起的误差场强。In the formula, is the field strength pattern of the ideal panel in the far field, and δE n,i is the error field strength caused by random manufacturing errors on the surface.

步骤(2)按如下过程进行:Step (2) is carried out as follows:

(2a)根据加工工艺确定每块面板表面随机误差的相关半径Lp(2a) Determine the relative radius L p of the random error on the surface of each panel according to the processing technology;

(2b)设单块面板上间隔为Δ的任意两点的表面误差分别为γn,i1),γn,i2),得到两者之间的误差相关函数(2b) Suppose the surface errors of any two points with an interval of Δ on a single panel are γ n,i1 ) and γ n,i2 ), respectively, and the error correlation function between them can be obtained

式中,为第(n,i)块面板的表面均方误差,Δ为矢量Δ的模值。In the formula, is the surface mean square error of the (n,i)th panel, and Δ is the modulus of the vector Δ.

步骤(6),按如下过程进行:Step (6), proceed as follows:

(6a)设定目标相对于坐标系O-xyz所在的方向(θ,φ)以方向余弦表示为(cosαx,cosαy,cosαz),并根据目标与坐标系的空间几何关系,得到目标相对于坐标轴的夹角与方向余弦的关系为:(6a) Set the direction (θ, φ) of the target relative to the coordinate system O-xyz to be expressed as (cosα x , cosα y , cosα z ) by the direction cosine, and obtain the target according to the spatial geometric relationship between the target and the coordinate system The relationship between the angle relative to the coordinate axis and the direction cosine is:

cosαx=sinθcosφ,cosαy=sinθsinφ,cosαz=cosθ;cosα x = sinθ cosφ, cosα y = sinθ sinφ, cosα z = cosθ;

式中,各个参数含义是:θ,φ是目标点的俯仰角与方位角;cosαx,cosαy,cosαz是目标方向单位矢量的方向余弦。In the formula, the meaning of each parameter is: θ, φ are the pitch angle and azimuth angle of the target point; cosα x , cosα y , cosα z are the direction cosines of the target direction unit vector.

(6b)得到口径面上任意两点的矢量与观察方向单位矢量的点乘(6b) Obtain the point product of the vector of any two points on the aperture surface and the unit vector of the viewing direction

式中,φ′1,φ′2分别为ρ1,ρ2的角度坐标。然后有In the formula, φ′ 1 and φ′ 2 are the angular coordinates of ρ 1 and ρ 2 respectively. then there is

(6c)每块面板上任意两点间的间距Δ用积分变量表示为(6c) The distance Δ between any two points on each panel is expressed by the integral variable as

(6d)单块面板的平均误差功率表达式为(6d) The average error power expression of a single panel is

上式是一个四重积分,可利用数值积分进行求解。The above equation is a quadruple integral, which can be solved by numerical integration.

步骤(7)按如下过程进行:Step (7) is carried out as follows:

(7a)整个反射面天线在远场的功率方向图为(7a) The power pattern of the whole reflector antenna in the far field is

上式的平均值就是天线的平均功率方向图,即The average value of the above formula is the average power pattern of the antenna, that is,

式中,E0为无误差时整个反射面的场强方向图;假设各块面板上的表面随机误差在相关半径内服从均值为0,方差为σ2的正态分布,则对误差场强的均值而言有 In the formula, E 0 is the field strength pattern of the entire reflector without error; assuming that the surface random error on each panel obeys a normal distribution with a mean value of 0 and a variance of σ 2 within the relevant radius, then the error field strength the mean value of

(7b)由于各分块面板独立制造,不同面板之间的表面随机误差并不具有相关性,因此在对所有误差功率的求和,即上式第四项时,不同面板之间的误差场强相乘的均值等于各自均值的乘积,同样为0,即(7b) Since each block panel is manufactured independently, the surface random errors between different panels are not correlated, so when summing all error powers, that is, the fourth term of the above formula, the error field between different panels The mean of the strong multiplication is equal to the product of the respective means, which is also 0, that is

故只剩下同块面板自身的功率在远场的叠加,因此整个天线在远场的平均辐射功率方向图为Therefore, only the superposition of the power of the same panel in the far field remains, so the average radiation power pattern of the entire antenna in the far field is

本发明与现有技术相比,具有如下优点:Compared with the prior art, the present invention has the following advantages:

1.把受表面误差影响的远场表示为理想场与误差场之和,在求天线平均功率时可将不同面板之间的误差功率简化为零处理,减少了计算量;1. The far field affected by the surface error is expressed as the sum of the ideal field and the error field, and the error power between different panels can be simplified to zero when calculating the average power of the antenna, which reduces the amount of calculation;

2.根据天线的实际分块形式计算表面加工误差对电性能的影响,更接近于工程实际;2. According to the actual block form of the antenna, the influence of the surface processing error on the electrical performance is calculated, which is closer to the actual engineering;

3.利用精确的四重数值积分求解单块面板的平均误差功率,避免了传统积分时间隔矢量Δ无法完整表示的缺陷;3. Use accurate quadruple numerical integration to solve the average error power of a single panel, avoiding the defect that the traditional integration time interval vector Δ cannot be fully represented;

4.通过仿真分析得出表面误差均方根的上限值,避免了凭经验提出过高的精度要求,降低了加工难度同时又满足电性能指标。4. Through the simulation analysis, the upper limit of the root mean square of the surface error is obtained, which avoids the excessive precision requirements raised by experience, reduces the processing difficulty and meets the electrical performance indicators at the same time.

附图说明Description of drawings

图1是本发明基于面板分块形式计算加工误差对电性能影响的流程图;Fig. 1 is the flow chart of calculating the influence of processing error on electrical performance based on panel block form in the present invention;

图2是本发明所用的面板拼接组成的反射面示意图;Fig. 2 is a schematic diagram of a reflective surface formed by panel splicing used in the present invention;

图3是本发明所用的单块面板的积分计算参数示意图;Fig. 3 is the integral calculation parameter schematic diagram of the single panel used in the present invention;

图4是目标与坐标系的空间几何关系的示意图;Fig. 4 is a schematic diagram of the spatial geometric relationship between the target and the coordinate system;

图5是本发明所采用的仿真示例图;Fig. 5 is the emulation example figure that the present invention adopts;

图6是本发明仿真结果的存在误差与理想情况下的电性能对比图。Fig. 6 is a comparison chart of the existing error of the simulation result of the present invention and the electrical performance under ideal conditions.

具体实施方式detailed description

以下参照附图对本发明作进一步详细描述。The present invention will be described in further detail below with reference to the accompanying drawings.

参照图1,本发明的具体步骤如下:With reference to Fig. 1, concrete steps of the present invention are as follows:

步骤一,将单块面板的远场表示为理想场与误差场之和;(1a)参照图2,由面板拼接而成的反射面共分为N环,第n环(n=1,2,…,N)又可分为Kn个子块。每块面板的表面都会存在制造误差,对第(n,i)块面板而言,设其表面某点处的随机误差引起的相位误差为δ(ρ),则其远场方向图可表示为Step 1, the far field of a single panel is expressed as the sum of the ideal field and the error field; (1a) Referring to Figure 2, the reflective surface formed by splicing panels is divided into N rings, the nth ring (n=1,2 ,…,N) can be divided into K n sub-blocks. There will be manufacturing errors on the surface of each panel. For the (n,i)th panel, assuming that the phase error caused by a random error at a certain point on the surface is δ(ρ), its far-field pattern can be expressed as

式中,为天线口径场幅度分布函数,其中B+C=1,通过改变P值可控制口径场幅度分布函数的形状,通常1≤P≤2,ρ为口径面上天线某点的半径,a为天线的半径。k为电磁波传播波数,为目标方向的单位矢量;In the formula, is the antenna aperture field amplitude distribution function, where B+C=1, the shape of the aperture field amplitude distribution function can be controlled by changing the P value, usually 1≤P≤2, ρ is the radius of a certain point of the antenna on the aperture surface, and a is the antenna of the radius. k is the wave number of electromagnetic wave propagation, is the unit vector of the target direction;

(1b)对上式中的ejδ(ρ)作一阶泰勒展开,有ejδ(ρ)≈1+jδ(ρ),然后得到(1b) Perform first-order Taylor expansion on e jδ(ρ) in the above formula, e jδ(ρ) ≈1 +jδ(ρ), and then get

式中,为理想面板在远场的场强方向图,δEn,i为表面随机制造误差引起的误差场强。In the formula, is the field strength pattern of the ideal panel in the far field, and δE n,i is the error field strength caused by random manufacturing errors on the surface.

步骤二,计算单块面板任意两点法向随机误差γn,i1)与γn,i2)乘积的均值<γn,i1n,i2)>。Step 2: Calculate the mean value of the product of random error γ n,i1 ) and γ n,i2 ) of any two points on a single panel<γ n,i1n,i2 )>.

(2a)根据加工工艺确定每块面板表面随机误差的相关半径Lp(2a) Determine the relative radius L p of the random error on the surface of each panel according to the processing technology;

(2b)参照图3,设单块面板上间隔为Δ的任意两点的表面误差分别为γn,i1),γn,i2),得到两者之间的误差相关函数(2b) Referring to Fig. 3, suppose the surface errors of any two points with an interval of Δ on a single panel are γ n,i1 ) and γ n,i2 ) respectively, and the error correlation between them can be obtained function

式中,为第(n,i)块面板的表面均方误差,Δ为矢量Δ的模值。In the formula, is the surface mean square error of the (n,i)th panel, and Δ is the modulus of the vector Δ.

步骤三,计算单块面板任意两点随机相位误差δn,i1)与δn,i2)乘积的均值<δn,i1n,i2)>。Step 3: Calculate the mean value of the product of random phase error δ n,i1 ) and δ n,i2 ) at any two points on a single panel<δ n,i1n,i2 ) >.

结合表面误差与空间相位误差的关系可得Combined with the relationship between the surface error and the spatial phase error, we can get

式中,λ为天线工作波长。In the formula, λ is the working wavelength of the antenna.

步骤四,推导得出单块面板在远场的平均误差功率表达式。Step 4, derive the expression of the average error power of a single panel in the far field.

每块面板的误差场强为The error field strength of each panel is

故平均误差功率为So the average error power is

步骤五,确定单块面板的积分区间。Step five, determine the integration interval of the single panel.

参照图3,单块面板误差场强的积分区间为:径向[a11,a12],周向[α1112];误差场强共轭值的积分区间为:径向[a21,a22],周向[α21,α22]。Referring to Figure 3, the integral interval of the error field strength of a single panel is: radial [a 11 , a 12 ], circumferential [α 11 , α 12 ]; the integral interval of the conjugate value of the error field strength is: radial [a 21 ,a 22 ], circumferential direction [ α2 1,α 22 ].

步骤六,积分计算单块面板平均误差功率。Step six, integrally calculate the average error power of a single panel.

(6a)参照图4,设定目标相对于坐标系O-xyz所在的方向(θ,φ)以方向余弦表示为(cosαx,cosαy,cosαz),并根据目标与坐标系的空间几何关系,得到目标相对于坐标轴的夹角与方向余弦的关系为:(6a) Referring to Figure 4, set the direction (θ, φ) of the target relative to the coordinate system O-xyz to be expressed as (cosα x , cosα y , cosα z ) by the direction cosine, and according to the spatial geometry of the target and the coordinate system The relationship between the angle of the target relative to the coordinate axis and the direction cosine is obtained as:

cosαx=sinθcosφ,cosαy=sinθsinφ,cosαz=cosθ;cosα x = sinθ cosφ, cosα y = sinθ sinφ, cosα z = cosθ;

(6b)参照图3以及图4,得到口径面上某点的矢量与观察方向单位矢量的点乘(6b) Referring to Figure 3 and Figure 4, the point product of the vector of a certain point on the aperture surface and the unit vector of the viewing direction is obtained

式中,φ′1,φ′2分别为ρ1,ρ2的角度坐标。然后有In the formula, φ′ 1 and φ′ 2 are the angular coordinates of ρ 1 and ρ 2 respectively. then there is

(6c)每块面板上任意两点间的间距Δ可用积分变量表示为(6c) The distance Δ between any two points on each panel can be expressed as

(6d)将步骤五、六中各式代入单块面板平均误差功率表达式中,得到(6d) Substitute the formulas in steps 5 and 6 into the expression of the average error power of a single panel, and get

上式是一个四重积分,可利用数值积分进行求解。The above equation is a quadruple integral, which can be solved by numerical integration.

步骤七,计算整个反射面天线的平均辐射功率。Step 7, calculate the average radiated power of the entire reflector antenna.

(7a)整个反射面天线在远场的功率方向图为(7a) The power pattern of the whole reflector antenna in the far field is

上式的平均值就是天线的平均功率方向图,即The average value of the above formula is the average power pattern of the antenna, that is,

式中,E0为无误差时整个反射面的场强方向图。假设各块面板上的表面随机误差在相关半径内服从均值为0,方差为σ2的正态分布,则对误差场强的均值而言有 In the formula, E 0 is the field strength pattern of the entire reflector without error. Assuming that the surface random error on each panel obeys a normal distribution with a mean of 0 and a variance of σ2 within the relevant radius, then for the mean value of the error field strength,

(7b)由于各分块面板独立制造,不同面板之间的表面随机误差并不具有相关性,因此在对所有误差功率的求和(即上式第四项)时,不同面板之间的误差场强相乘的均值等于各自均值的乘积,同样为0,即(7b) Since each block panel is manufactured independently, the surface random errors between different panels are not correlated, so when summing all error powers (that is, the fourth term of the above formula), the error between different panels The mean value of the field strength multiplication is equal to the product of the respective mean values, which is also 0, that is

故只剩下同块面板自身的功率在远场的叠加,因此整个天线在远场的平均辐射功率方向图为Therefore, only the superposition of the power of the same panel in the far field remains, so the average radiation power pattern of the entire antenna in the far field is

步骤八,把不同的面板表面加工误差均方根代入上式,仿真分析其对电性能的影响,进而提出合理的表面加工精度要求。Step 8: Substitute the root mean square error of different panel surface processing into the above formula, simulate and analyze its impact on electrical performance, and then put forward reasonable surface processing accuracy requirements.

本发明的优点可通过以下仿真实验进一步说明:Advantages of the present invention can be further illustrated by following simulation experiments:

1.仿真条件1. Simulation conditions

将本发明的基于分块形式计算面板加工误差对电性能影响的方法编制为Matlab程序,在一35米分块反射面天线上进行电性能仿真分析。The method for calculating the impact of panel processing errors on electrical performance based on the block form of the present invention is compiled into a Matlab program, and the electrical performance simulation analysis is performed on a 35-meter block reflector antenna.

该反射面天线的分块结构示意图如图5所示,其反射面共分为9环,1~4环每环均由16块面板构成,第5环分为32块,6~9环每环又分为96块。天线工作频率为10GHz,焦径比为F/D=0.35,边缘照射锥削为ET=-10dB,口径场照射函数为表面加工误差相关半径Lp=0.4m。The block structure schematic diagram of the reflecting surface antenna is shown in Figure 5. The reflecting surface is divided into 9 rings. Each ring of the 1st to 4th rings is composed of 16 panels. The fifth ring is divided into 32 panels. The ring is further divided into 96 pieces. The working frequency of the antenna is 10GHz, the focal diameter ratio is F/D=0.35, the edge irradiation taper is ET=-10dB, and the aperture field irradiation function is Surface processing error-related radius L p =0.4m.

2.仿真结果2. Simulation results

设每块面板的表面误差均方根均为ε=λ/30,利用编制的Matlab程序对以上模型进行仿真,得到天线在存在表面制造误差时和理想情况下的方向图对比,如图6所示。从图中可以看出,在存在较小表面制造误差时,天线的主瓣展宽,近副瓣影响较小,但远区副瓣抬高明显。从仿真分析结果可知,采用本发明的方法可用于指导大型反射面天线的面板加工制造以及电性能的分析与评价。Assuming that the root mean square of the surface error of each panel is ε=λ/30, the above model is simulated by using the compiled Matlab program, and the comparison of the antenna pattern when there is surface manufacturing error and the ideal situation is obtained, as shown in Figure 6 Show. It can be seen from the figure that when there is a small surface manufacturing error, the main lobe of the antenna is broadened, and the influence of the near side lobe is small, but the side lobe in the far area is significantly raised. It can be seen from the simulation analysis results that the method of the invention can be used to guide the panel processing and manufacturing of large reflector antennas and the analysis and evaluation of electrical performance.

Claims (5)

1. the method that panel mismachining tolerance is affected on electrical property is calculated based on block form, be it is characterized in that:At least include following step Suddenly:
(1) the phase effect item in the integral expression of monolithic panel far field is carried out first order Taylor expansion, far field is expressed as into reason Think field and error field sum;
(2) the random error correlation radius of every piece of panel is determined according to processing technique, same block panel various location is then obtained Error correlation function, finally try to achieve every piece of panel any two points normal direction random error γn,i1) and γn,i2) product it is equal Value<γn,i1n,i2)>;
(3) relation using normal error with phase error tries to achieve every piece of panel any two points phase error δn,i1) and δn,i2) product average<δn,i1n,i2)>;
(4) by panel far field error field intensity δ En,iIt is conjugated with whichIt is multiplied and is averaging, is derived from monolithic panel and exists The average error power expression formula in far field
(5) the integrating range parameter of every piece of panel, including radial direction parameter are determined, i.e.,:The inside and outside footpath of monolithic panel:a11, a12, a21, a22;Circumferential parameter, i.e.,:The angle coordinate of monolithic panel two ends string:α11, α12, α21, α22
(6) step (3) is drawn<δn,i1n,i2)>Substitute in average error power expression formula, and integrand is used Integration variable is expressed, and then obtains monolithic panel average error power by quadruple numerical integration;
(7) average error power by all panels in far field is superimposed, and power of the acquired results with antenna in the ideal case is asked With finally give average power pattern of the whole antenna containing Surface Machining error in far field.
2. according to claim 1 that the method that panel mismachining tolerance is affected on electrical property is calculated based on block form, which is special Levy and be:The phase effect item in the integral expression of monolithic panel far field is carried out first order Taylor expansion described in step (1), will Far field is expressed as ideal field and error field sum, carries out according to the following procedure:
(1a) N rings are divided into by the reflecting surface of panel splicing, the n-th ring can be divided into K againnIndividual sub-block, wherein, n=1,2 ..., N;Can all there is foozle in the surface of every piece of panel, for (n, i) block panel, if the random error at its surface point The phase error for causing is δ (ρ), then its far-field pattern is represented by
E n , i = &Integral; &Integral; ( n , i ) z o n e Q ( &rho; ) e j &delta; ( &rho; ) e j k &rho; , r ^ dS n , i
In formula,For antenna aperture field amplitude distribution function, wherein B+C=1, ρ go up to the sky for bore face The radius of line point, radiuses of a for antenna, k are electromagnetic wave propagation wave number,For the unit vector of target direction;
(1b) to the e in above formulajδ(ρ)Make first order Taylor expansion, have ejδ(ρ)≈ 1+j δ (ρ), then obtains
E n , i &ap; &Integral; &Integral; ( n , i ) z o n e Q ( &rho; ) &lsqb; 1 + j &delta; ( &rho; ) &rsqb; e j k &rho; &RightArrow; &CenterDot; r ^ dS n , i = E n , i 0 + &delta;E n , i
In formula,For ideal panel far field field strength pattern, δ En,iFor the error field intensity that surface random manufacturing error causes.
3. according to claim 1 that the method that panel mismachining tolerance is affected on electrical property is calculated based on block form, which is special Levy and be:Step (2) is carried out according to the following procedure:
(2a) the correlation radius L of every piece of panel surface random error is determined according to processing techniquep
(2b) surface error for setting any two points on monolithic panel at intervals of Δ is respectively γn,i1), γn,i2), obtain two Error correlation function between person
C n , i ( &Delta; ) = < &gamma; n , i ( &rho; 1 ) &gamma; n , i ( &rho; 2 ) > &epsiv; n , i 2 = exp &lsqb; - ( &Delta; / L p ) 2 &rsqb;
In formula,For the surface mean square error of (n, i) block panel, modulus value of the Δ for vector Δ.
4. according to claim 1 that the method that panel mismachining tolerance is affected on electrical property is calculated based on block form, which is special Levy and be:Step (6), is carried out according to the following procedure:
(6a) direction (θ, φ) that target setting is located relative to coordinate system O-xyz is expressed as (cos α with direction cosinesx,cosαy, cosαz), and according to target and the space geometry relation of coordinate system, obtain angle and direction cosines of the target relative to coordinate axess Relation be:
cosαx=sin θ cos φ, cos αy=sin θ sin φ, cos αz=cos θ;
In formula, parameters are meant that:θ, φ are the pitching of impact point and azimuth;cosαx,cosαy,cosαzIt is object vector Direction cosines;
(6b) obtain the relation of the vector and direction of observation unit vector of certain point on bore face
&rho; 1 &CenterDot; r ^ = &rho; 1 sin&theta;cos&phi;cos&phi; 1 &prime; + &rho; 1 sin&theta;sin&phi;sin&phi; 1 &prime; &rho; 2 &CenterDot; r ^ = &rho; 2 sin&theta;cos&phi;cos&phi; 2 &prime; + &rho; 2 sin&theta;sin&phi;sin&phi; 2 &prime;
In formula, φ '1, φ '2Respectively ρ1, ρ2Angle coordinate;Then have
( &rho; 1 - &rho; 2 ) &CenterDot; r ^ = &rho; 1 s i n &theta; c o s ( &phi; - &phi; 1 &prime; ) - &rho; 2 s i n &theta; c o s ( &phi; - &phi; 2 &prime; )
(6c) on every piece of panel, the separation delta of any point-to-point transmission is expressed as with integration variable
&Delta; = | &Delta; | = ( &rho; 1 cos&phi; 1 &prime; - &rho; 2 cos&phi; 2 &prime; ) 2 + ( &rho; 1 sin&phi; 1 &prime; - &rho; 2 sin&phi; 2 &prime; ) 2 = &rho; 1 2 + &rho; 2 2 - 2 &rho; 1 &rho; 2 cos ( &phi; 1 &prime; - &phi; 2 &prime; )
(6d) by various substitution monolithic panel mean power expression formula in step 2, three, obtain
&delta;E n , i &CenterDot; &delta;E n , i * &OverBar; = &Integral; &alpha; 11 &alpha; 12 &Integral; &alpha; 11 &alpha; 12 &Integral; &alpha; 21 &alpha; 22 &Integral; &alpha; 21 &alpha; 22 Q ( &rho; 1 ) Q * ( &rho; 2 ) &CenterDot; ( 4 &pi; &lambda; ) 2 &epsiv; n , i 2 exp &lsqb; - ( &Delta; / L p ) 2 &rsqb; &CenterDot; e j k &lsqb; &rho; 1 sin &theta; cos ( &phi; - &phi; 1 &prime; ) - &rho; 2 sin &theta; cos ( &phi; - &phi; 2 &prime; ) &rsqb; &rho; 1 &rho; 2 d&rho; 2 d&phi; 2 &prime; d&rho; 1 d&phi; 1 &prime;
Above formula is a quad-slope integration, is solved using numerical integration.
5. according to claim 1 that the method that panel mismachining tolerance is affected on electrical property is calculated based on block form, which is special Levy is that step (7) is carried out according to the following procedure:
(7a) whole power radiation pattern of the reflector antenna in far field is
P = E &CenterDot; E * = &lsqb; &Sigma; n = 1 N &Sigma; i = 1 K n ( E n , i 0 + &delta;E n , i ) &rsqb; &CenterDot; &lsqb; &Sigma; m = 1 N &Sigma; l = 1 K m ( E m , l 0 * + &delta;E m , l * ) &rsqb; = &Sigma; n = 1 N &Sigma; m = 1 N &Sigma; i = 1 K n &Sigma; l = 1 K m ( E n , i 0 E m , l 0 * + E n , i 0 &delta;E m , l * + &delta;E n , i E m , l 0 * + &delta;E n , i &delta;E m , l * )
The meansigma methodss of above formula are exactly the average power pattern of antenna, i.e.,
P &OverBar; = E 0 &CenterDot; E 0 * + &Sigma; n = 1 N &Sigma; m = 1 N &Sigma; i = 1 K n &Sigma; l = 1 K m E n , i 0 &CenterDot; &delta;E m , l * &OverBar; + &Sigma; n = 1 N &Sigma; m = 1 N &Sigma; i = 1 K n &Sigma; l = 1 K m E m , l 0 * &CenterDot; &delta;E n , i &OverBar; + &Sigma; n = 1 N &Sigma; m = 1 N &Sigma; i = 1 K n &Sigma; l = 1 K m &delta;E n , i &CenterDot; &delta;E m , l * &OverBar;
In formula, E0For it is error free when whole reflecting surface field strength pattern;Assume the surface random error on each piece of panel in correlation half It is 0 that average is obeyed in footpath, and variance is σ2Normal distribution, then have for the average of error field intensity
(7b) as each piecemeal panel is independently manufactured, the surface random error between different panels does not have dependency, therefore In the summation to all error powers, i.e. above formula Section 4, the average that the error field intensity between different panels is multiplied is equal to each From the product of average, 0 is similarly, i.e.,
&Sigma; n = 1 N &Sigma; m = 1 N &Sigma; i = 1 K n &Sigma; l = 1 K m &delta;E n , i &CenterDot; &delta;E m , l * &OverBar; = &Sigma; n = 1 N &Sigma; m = 1 , m &NotEqual; n N &Sigma; i = 1 K n &Sigma; l = 1 K m &delta;E n , i &CenterDot; &delta;E m , l * &OverBar; + &Sigma; n = 1 N &Sigma; i = 1 K n &Sigma; l = 1 , l &NotEqual; i K m &delta;E n , i &CenterDot; &delta;E m , l * &OverBar; + &Sigma; n = 1 N &Sigma; i = 1 K n &delta;E n , i &CenterDot; &delta;E n , i * &OverBar; = &Sigma; n = 1 N &Sigma; i = 1 K n &delta;E n , i &CenterDot; &delta;E n , i * &OverBar;
Therefore only it is left the superposition with the block panel power of itself in far field, therefore whole antenna is in the average radiating power side in far field To figure it is
P &OverBar; = E 0 &CenterDot; E 0 * + &Sigma; n = 1 N &Sigma; i = 1 K n &delta;E n , i &CenterDot; &delta;E n , i * &OverBar; .
CN201410257555.4A 2014-06-11 2014-06-11 The method that panel mismachining tolerance is affected on electrical property is calculated based on block form Expired - Fee Related CN104063587B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410257555.4A CN104063587B (en) 2014-06-11 2014-06-11 The method that panel mismachining tolerance is affected on electrical property is calculated based on block form

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410257555.4A CN104063587B (en) 2014-06-11 2014-06-11 The method that panel mismachining tolerance is affected on electrical property is calculated based on block form

Publications (2)

Publication Number Publication Date
CN104063587A CN104063587A (en) 2014-09-24
CN104063587B true CN104063587B (en) 2017-04-05

Family

ID=51551297

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410257555.4A Expired - Fee Related CN104063587B (en) 2014-06-11 2014-06-11 The method that panel mismachining tolerance is affected on electrical property is calculated based on block form

Country Status (1)

Country Link
CN (1) CN104063587B (en)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107679336B (en) * 2017-10-20 2021-02-05 西安电子科技大学 Reflecting surface antenna surface random error analysis method based on second order approximation formula
CN109408958B (en) * 2018-10-23 2022-12-23 西安电子科技大学 Modeling Method of Average Power Pattern of Umbrella Antenna Considering Patch Composition Error
CN109408957B (en) * 2018-10-23 2022-12-23 西安电子科技大学 Umbrella-shaped antenna random error analysis method considering deterministic errors

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102253290A (en) * 2011-03-29 2011-11-23 王从思 Method for predicting electrical properties of deformed log-periodic antennae based on electromechanical coupling model
CN102890741A (en) * 2012-10-25 2013-01-23 西安电子科技大学 Method for predicting influence of vibration deformation on electric performance of array antenna

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102253290A (en) * 2011-03-29 2011-11-23 王从思 Method for predicting electrical properties of deformed log-periodic antennae based on electromechanical coupling model
CN102890741A (en) * 2012-10-25 2013-01-23 西安电子科技大学 Method for predicting influence of vibration deformation on electric performance of array antenna

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
Research on electromechanical field coupling model and solution strategy of large reflector antennas;Meng Wang等;《COUPLED 2013》;20130630;全文 *
结构误差对阵列天线极化特性的影响分析;王猛等;《系统工程与电子技术》;20121130;第34卷(第11期);第2193-2197页 *
面板位置误差对反射面天线功率方向图的影响机理;王伟等;《西安电子科技大学学报(自然科学版)》;20090831;第36卷(第4期);第708-713页 *

Also Published As

Publication number Publication date
CN104063587A (en) 2014-09-24

Similar Documents

Publication Publication Date Title
CN101267062B (en) Method for predicting antenna electric performance based on simulated distortion reflective side
Wang et al. A novel integrated simulation approach couples MCRT and Gebhart methods to simulate solar radiation transfer in a solar power tower system with a cavity receiver
CN102226840B (en) Radar cross-section layered calculation method of ship target within atmospheric duct range
CN101308177B (en) Initiative reflecting plane antenna electrical behavior prediction method
CN102243276B (en) Electric coupling analysis method of temperature loader of large-scale foundation surface antenna
CN102968532B (en) The dynamo-electric integral design method of the large-scale reflector antenna structure of 65m bore
CN104484537B (en) Electromechanical integrated optimization method for sectional dimensions of ribs of metal truss type antenna housing
CN106886021B (en) High Resolution Spaceborne SAR image quality method for improving
CN105430740A (en) Indoor wireless positioning method based on WiFi signal strength simulation and location fingerprint algorithm
CN104063587B (en) The method that panel mismachining tolerance is affected on electrical property is calculated based on block form
CN104898264B (en) A kind of rotationally symmetrical non-imaged free-form surface mirror of heavy caliber and its design method
CN106209269A (en) The calibration steps of sphere composite array near-field effect in a kind of radio freqency simulation system
CN101252226B (en) Method for positioning reflecting plane aerial feed source
CN106597364A (en) Target radiation source initial position estimation method for single-antenna single-station passive positioning
CN103037389A (en) Cable net reflective surface antenna molded surface shaping method based on sensitivity information
CN103604944A (en) Surface flow measurement method based on monostation shipborne high-frequency ground wave radar
He et al. Monte Carlo calculation of view factors between some complex surfaces: Rectangular plane and parallel cylinder, rectangular plane and torus, especially cold-rolled strip and W-shaped radiant tube in continuous annealing furnace
CN104021241B (en) A kind of earth light optical pressure perturbation modeling method of navigation satellite complex model
CN104063544A (en) Near Field Effect Error Analysis Method
CN105334739A (en) FAST whole network control method based on P type learning law of iterative learning
CN100574007C (en) Adjust the large-scale antenna reflecting plane gravity preset method at angle based on the best
CN105372633A (en) Phased array radar dimension reduction four-channel mainlobe sidelobe interference-resisting method
Wang et al. Edge sensor based real-time calculation method for active panel position of QTT
Rotgerink et al. Modelling of effects of nose radomes on radar antenna performance
CN114707380B (en) A method for adjusting and calculating accuracy of the main reflecting surface actuator of a submillimeter wave telescope

Legal Events

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

Granted publication date: 20170405