CN114065559A - 一种基板集成波导天线的设计方法 - Google Patents

一种基板集成波导天线的设计方法 Download PDF

Info

Publication number
CN114065559A
CN114065559A CN202210053123.6A CN202210053123A CN114065559A CN 114065559 A CN114065559 A CN 114065559A CN 202210053123 A CN202210053123 A CN 202210053123A CN 114065559 A CN114065559 A CN 114065559A
Authority
CN
China
Prior art keywords
cylindrical conductor
cylindrical
basic solution
wave
waveguide
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
CN202210053123.6A
Other languages
English (en)
Other versions
CN114065559B (zh
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.)
Changzhou Xingyu Automotive Lighting Systems Co Ltd
Original Assignee
Changzhou Xingyu Automotive Lighting Systems Co Ltd
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 Changzhou Xingyu Automotive Lighting Systems Co Ltd filed Critical Changzhou Xingyu Automotive Lighting Systems Co Ltd
Priority to CN202210053123.6A priority Critical patent/CN114065559B/zh
Publication of CN114065559A publication Critical patent/CN114065559A/zh
Application granted granted Critical
Publication of CN114065559B publication Critical patent/CN114065559B/zh
Priority to JP2023517786A priority patent/JP2024504894A/ja
Priority to EP22865879.5A priority patent/EP4239514A4/en
Priority to PCT/CN2022/126192 priority patent/WO2023138127A1/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • H—ELECTRICITY
    • H01—ELECTRIC ELEMENTS
    • H01Q—ANTENNAS, i.e. RADIO AERIALS
    • H01Q13/00—Waveguide horns or mouths; Slot antennas; Leaky-waveguide antennas; Equivalent structures causing radiation along the transmission path of a guided wave
    • H01Q13/06—Waveguide mouths
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • G06F17/13—Differential equations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/30—Circuit design
    • G06F30/36—Circuit design at the analogue level
    • G06F30/367—Design verification, e.g. using simulation, simulation program with integrated circuit emphasis [SPICE], direct methods or relaxation methods
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/30—Circuit design
    • G06F30/39—Circuit design at the physical level
    • G06F30/398—Design verification or optimisation, e.g. using design rule check [DRC], layout versus schematics [LVS] or finite element methods [FEM]
    • H—ELECTRICITY
    • H01—ELECTRIC ELEMENTS
    • H01Q—ANTENNAS, i.e. RADIO AERIALS
    • H01Q1/00—Details of, or arrangements associated with, antennas
    • H01Q1/36—Structural form of radiating elements, e.g. cone, spiral, umbrella; Particular materials used therewith
    • H—ELECTRICITY
    • H01—ELECTRIC ELEMENTS
    • H01Q—ANTENNAS, i.e. RADIO AERIALS
    • H01Q1/00—Details of, or arrangements associated with, antennas
    • H01Q1/50—Structural association of antennas with earthing switches, lead-in devices or lightning protectors
    • H—ELECTRICITY
    • H01—ELECTRIC ELEMENTS
    • H01Q—ANTENNAS, i.e. RADIO AERIALS
    • H01Q21/00—Antenna arrays or systems
    • H01Q21/06—Arrays of individually energised antenna units similarly polarised and spaced apart
    • H01Q21/061—Two dimensional planar arrays
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2115/00—Details relating to the type of the circuit
    • G06F2115/12—Printed circuit boards [PCB] or multi-chip modules [MCM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2117/00—Details relating to the type or aim of the circuit design
    • G06F2117/12—Sizing, e.g. of transistors or gates
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Computer Hardware Design (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Pure & Applied Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Operations Research (AREA)
  • Microelectronics & Electronic Packaging (AREA)
  • Waveguides (AREA)
  • Variable-Direction Aerials And Aerial Arrays (AREA)
  • Waveguide Aerials (AREA)

Abstract

本发明公开了一种基板集成波导天线的设计方法，基板集成波导天线包括由若干阵列布置的圆柱导体组成的波导组件和连接于波导组件轴向两端的介质基板，其中一个介质基板连接一同轴电缆；所述方法包括以下步骤：S1：计算散射波的基本解；S2：利用多种边界条件将S1中散射波的基本解进一步变形；S3：采用全域基函数的矩量法，将电磁场的计算转化为关于平面内柱面坐标的计算问题。本发明在对电磁波辐射进行分析时，采用了基于全域基函数的矩量法，将天线设计的主要问题转化为关于柱面坐标的计算问题，大大减少了仿真计算的数据量。

Description

一种基板集成波导天线的设计方法
技术领域
本发明涉及天线技术领域，尤其涉及一种基板集成波导天线的设计方法。
背景技术
近年来，无线通信技术得到了快速发展并获得了广泛应用。在通信系统中，由于有源器件的非线性特性，在工作频带外会产生大量的杂波信号，一方面会恶化系统性能，另一方面会对其他系统产生干扰。因此，对天线和滤波器的频率选择性提出了很高的要求，从而达到抑制带外杂波和干扰的目的。
基板集成波导是近年来研究和应用较多的波导结构，它不仅保留了传统金属波导的辐射损耗小、插入损耗低、功率容量大等特点，同时具有成本低、尺寸小、重量轻、易集成等优点，因此具有很高的工程应用价值。
常用天线的频率选择特性很差，对系统的带外杂波和干扰抑制不明显，因此需要一种具有良好频率选择特性的天线，抑制带外杂波和干扰，从而降低对系统中其他器件设计指标的要求，降低系统的实现成本。而且目前由于针对天线的波导仿真计算数据很大，很难计算出合理的天线参数值，使天线的散射波强度达到最佳。
发明内容
本发明提出一种基板集成波导天线的设计方法，用于解决目前针对天线的波导仿真计算数据很大，很难计算出合理的天线参数值的技术问题。
一种基板集成波导天线，包括由若干阵列布置的圆柱导体组成的波导组件和连接于波导组件轴向两端的介质基板，其中介质基板连接一同轴电缆。
一种基板集成波导天线的设计方法，所述方法采用以上所述的基板集成波导天线，包括以下步骤：
S1：计算散射波的基本解；
S2：利用多种边界条件将S1中散射波的基本解进一步变形；
S3：采用全域基函数的矩量法，将电磁场水平和垂直方向的计算转化为关于柱坐标下水平方向的计算问题。
进一步的，步骤S1中散射波的基本解的计算过程如下：
S11：列出入射波的电场强度公式：
Figure 41527DEST_PATH_IMAGE001
，其中，
Figure 372014DEST_PATH_IMAGE002
（
Figure 364241DEST_PATH_IMAGE003
为时间因子，φ0为入射角度，k0为传播常数，μ0为真空中的透磁率，ε0为真空中的介电常数，ρ为圆柱导体1上任一点的柱面坐标，E0为入射波的电场强度）。
S12：由圆柱导体上的电磁感应电流产生的散射波满足下述的亥姆霍兹波动方程式：
Figure 121981DEST_PATH_IMAGE004
，其中，
Figure 566869DEST_PATH_IMAGE005
（η为传播导纳，
Figure 701047DEST_PATH_IMAGE006
为圆柱导体上的电磁感应电流，
Figure 10000254929754
为媒介的电导率）。
S13：将S12中的方程式转化为如下微分方程式：
Figure 813359DEST_PATH_IMAGE007
，并求其基本解为：
Figure 414105DEST_PATH_IMAGE008
。
S14：在S12中的参数代入S13的基本解中，得出散射波的基本解为：
Figure 736502DEST_PATH_IMAGE009
，其中，C=C1+C2+…+CN。
进一步的，所述步骤S2包括：
S21：根据圆柱导体上的全电场等于零得出方程式：
Figure 284158DEST_PATH_IMAGE010
；并将该方程式代入S14的基本解中，并得出基本解的第一转化形式：
Figure 110032DEST_PATH_IMAGE011
，其中，
Figure 350520DEST_PATH_IMAGE012
。
S22：将圆柱导体表面的点用极坐标轴表示，则圆柱导体上任一点的坐标函数为
Figure 894634DEST_PATH_IMAGE013
；将该坐标函数代入步骤S21中基本解的第一转化形式，得出基本解的第二转化形式：
Figure 245981DEST_PATH_IMAGE014
。
S23：将步骤S22中基本解的第二转化形式积分，得出基本解的第三转化形式：
Figure 191940DEST_PATH_IMAGE015
（av表示第v号圆柱导体的半径，au表示第u号圆柱导体的半径）。
进一步的，所述步骤S3包括：
S31：将步骤S23中基本解的第三转化形式转化为如下的联立线性方程式：
Figure 400067DEST_PATH_IMAGE016
，其中A表示圆柱导体的电流感应，Aii表示第i个圆柱导体的自感，Aij表示第i个圆柱导体对第j个圆柱导体的互感，Xi表示第i个圆柱导体产生的散射波强度，Bi表示第i个圆柱导体上的电流分布；
S32：将系数行列式对称化，通过矩阵变换得出如下散射波强度的公式：
Figure 306844DEST_PATH_IMAGE017
，其中，
Figure 320936DEST_PATH_IMAGE018
进一步的，令所述圆柱导体的直径为D，沿平行波导通道的方向上，相邻两个圆柱导体的中心距为P，则D/P=0.5。
优选的，所述圆柱导体为金属管。
进一步的，所述同轴电缆设置于所述波导组件的中心线上。
进一步的，输入电磁波的波长为λ；沿垂直波导通道的方向上，相邻两个圆柱导体的中心距为D1，则满足如下关系：D1=0.9λ。
进一步的，所述同轴电缆至介质基板闭口处的距离D2与输入电磁波的波长λ的关系为D2=0.25λ。
进一步的，所述波导组件的波导通道长度D3与输入电磁波的波长λ的关系为D3=2.4λ。
本发明的有益效果是：
（1）本发明所述的基板集成波导天线的设计方法，在对电磁波辐射进行分析时，采用了基于全域基函数的矩量法，将天线设计的主要问题转化为关于柱面坐标的计算问题，大大减少了仿真计算的数据量。
（2）本发明所述的基板集成波导天线的设计方法，通过圆柱导体的直径D和沿波导通道方向相邻圆柱导体的中心距P的合理配比，使天线结构简化的同时避免漏波现象发生，达到最优性价比的设计性能。
附图说明
下面结合附图和实施例对本发明进一步说明。
图1 是本发明所述的基板集成波导天线的具体实施方式的示意图；
图2 是图1所示基板集成波导天线的具体实施例的电场分布图；
图3 是均匀平面波垂直入射到两种不同媒质的分界面时的电场和磁场方向示意图。
图中，1、圆柱导体，2、介质基板，3、同轴电缆，4、波导通道。
具体实施方式
下面详细描述本发明的实施例，所述实施例的示例在附图中示出，其中自始至终相同或类似的标号表示相同或类似的元件或具有相同或类似功能的元件。下面通过参考附图描述的实施例是示例性的，仅用于解释本发明，而不能理解为对本发明的限制。
如图1所示，一种基板集成波导天线，包括由若干阵列布置的圆柱导体1组成的波导组件和连接于波导组件轴向两端的介质基板2，波导组件通常呈长方形布置，相邻两排圆柱导体1之间形成波导通道4，波导通道4沿波导组件的长度方向延伸，波导通道4的长度F即为波导组件的长度，两个介质基板2与每个圆柱导体1连接，其中介质基板2连接一同轴电缆3。
本发明假设圆柱导体1的轴向方向为Z轴方向，平行圆柱导体1的中心轴的平面内，与波导通道4平行的方向为X轴方向，与波导通道4垂直的方向为Y轴方向。
本实施例中，圆柱导体1为金属管，圆柱导体1共设有两排，每排设有五根圆柱导体1，两排圆柱导体1形成导波组件的两个侧壁，两个所述侧壁之间形成波导通道4，介质基板2采用单层PCB板工艺，降低成本，且容易批量生产。同轴电缆3通常由同轴由内至外依次套设的内导体、绝缘体和外导体组成，同轴电缆3的内导体插入到介质基板2内部，同轴电缆3的外导体与介质基板2的表面连接。同轴电缆3优选设置于波导组件的中心线上。
如图3所示，为均匀平面波垂直入射到两种不同媒质（煤质一和煤质二）的分界面时的电场和磁场方向示意图，图中Ei、Et、Er分别为入射波、透射波和反射波的电场方向，Hi为入射波的磁场方向，ki为入射方向介质的传播常数，Ht为透射波的磁场方向，kt为透射方向介质的传播常数，Hr为反射波的磁场方向，kr为反射方向介质的传播常数。
本发明中，煤质一为空气，煤质二为圆柱导体1，当观察点
Figure 262347DEST_PATH_IMAGE019
沿圆柱导体1的圆柱表面移动时，圆柱导体1上的电场强度总和等于0，由于入射波垂直入射圆柱导体1（在本发明的方向定义中，则入射波的入射方向与Z轴方向平行），故只产生圆柱导体1水平方向上的磁场和垂直方向上的电场（即磁场方向与XOY平面平行，电场方向与Z轴方向平行），而电场在圆柱导体1表面又呈均匀分布，故在仿真分析时，只需要考虑圆柱导体1断面的电场分布情况（即电场在XOY平面的分布情况）。故本发明提供一种基板集成波导天线的设计方法，该方法在对电磁波辐射进行分析时，采用了基于全域基函数的矩量法，电磁场的计算问题可以由传统的三个坐标轴方向的计算问题转化为关于平面内柱面坐标的计算问题，大大减少了仿真计算的数据量。
所述方法包括以下步骤：
S1：计算散射波的基本解；
S2：利用多种边界条件将S1中散射波的基本解进一步变形；
S3：采用全域基函数的矩量法，将电磁场水平和垂直方向的计算转化为关于柱坐标下水平方向的计算问题。此处所述水平方向是指水平平面，竖直方向是指与水平平面垂直的方向。
大多数的电磁场问题都是基于如下初始公式：
Figure 438114DEST_PATH_IMAGE020
，其中，L是线性算子，g是已知函数，h是未知函数。由于对该公式直接求严密解难度较大，所以根据线性算子L的线性特性，求解其数值解，其中使用了矩量法，具体方法如下：
首先，将未知函数h按如下方式展开成N个基底函数
Figure 832186DEST_PATH_IMAGE021
，其中，an是未知系数。
接着将关于h的函数代入初始公式中，通过对线性算子L进行线性变换，得到如下公式：
Figure 649969DEST_PATH_IMAGE022
。
最后，定义N个权重函数Wm，对线性变换后的公式进行积分，得到如下等式：
Figure 445887DEST_PATH_IMAGE023
，该等式也可通过如下行列式表示：
Figure 323713DEST_PATH_IMAGE024
，其中，
Figure 1819DEST_PATH_IMAGE025
将
Figure 233080DEST_PATH_IMAGE026
代入上述行列式中可得：
Figure 8138DEST_PATH_IMAGE027
，则未知函数h可通过如下公式表示：
Figure 932232DEST_PATH_IMAGE028
。
下面将上述的推演步骤应用至本发明所述的基板集成波导天线的电磁场分析中，具体计算过程如下：
步骤S1中散射波的基本解的计算过程如下：
S11：若时间因子为e-jωt（t为时间，ω=2πf，ω为震荡角频率，f为频率），基板集成波导天线的电磁场分量为时谐场，即随时间变化而变化，则入射波Z向传播的磁场强度公式为：
Figure 159951DEST_PATH_IMAGE001
，其中，
Figure 194903DEST_PATH_IMAGE002
，
Figure 824467DEST_PATH_IMAGE003
为时间因子，φ0为入射角度，k0为传播常数，μ0为真空中的透磁率，ε0为真空中的介电常数，ρ为圆柱导体1上任一点的柱面坐标，E0为入射波的电场强度。
S12：由圆柱导体1上的电磁感应电流产生的散射波满足下述的亥姆霍兹波动方程式：
Figure 919462DEST_PATH_IMAGE029
，其中，
Figure 900057DEST_PATH_IMAGE005
（η为传播导纳，
Figure 473120DEST_PATH_IMAGE006
为圆柱导体1上的电磁感应电流，
Figure 10000254943452
为媒介的电导率）。
S13：将S12中的方程式转化为如下微分方程式：
Figure 222771DEST_PATH_IMAGE007
，并求其基本解为：
Figure 488667DEST_PATH_IMAGE008
。
S14：将S12中的参数代入S13的基本解中，得出散射波的基本解为：
Figure 363082DEST_PATH_IMAGE009
，其中，C=C1+C2+…+CN。
接着通过多个边界条件对上述基本解进行转化，具体步骤如下：
S21：考察将观察点
Figure 130050DEST_PATH_IMAGE019
在圆柱导体1上移动，则圆柱导体1上的全电场等于0，即：
Figure 609573DEST_PATH_IMAGE010
；将该方程式代入S14的基本解中，并得出基本解的第一转化形式：
Figure 905425DEST_PATH_IMAGE011
，其中，
Figure 470398DEST_PATH_IMAGE030
为第v号的圆柱导体1上的电流分布。
又因为电流分布是沿着圆柱导体1的圆柱表面呈圆形分布的，即电流分布是周期2π的相关函数，所以，第v号的圆柱导体1上的电流分布可以由如下公式表示：
Figure 509898DEST_PATH_IMAGE012
。
S22：因为圆柱导体1的断面是圆形，所以通过坐标系变换，将圆柱导体1表面的点用极坐标轴来表示的话，可得圆柱导体1上任一点的坐标函数为
Figure 109507DEST_PATH_IMAGE013
；将该坐标函数代入步骤S21中基本解的第一转化形式，得出基本解的第二转化形式：
Figure 513944DEST_PATH_IMAGE031
。
S23：将步骤S22中基本解的第二转化形式左右两边都乘以
Figure 690847DEST_PATH_IMAGE032
，并沿Cu积分，即关于φ沿-π到π进行积分，可得基本解的第三转化形式：
Figure 409404DEST_PATH_IMAGE033
（av表示第v号圆柱导体1的半径，au表示第u号圆柱导体1的半径）。
为了计算未知系数
Figure 988153DEST_PATH_IMAGE034
，利用加法定理对上述公式的左右两边进行改写，首先左边改写成向量行列式：
Figure 766753DEST_PATH_IMAGE035
，其中，
Figure 165374DEST_PATH_IMAGE036
。所以，
Figure 687622DEST_PATH_IMAGE037
，接着对式子的右边改写，将关于
Figure 120877DEST_PATH_IMAGE034
的求和式提出，剩余的部分写成系数行列式
Figure 867116DEST_PATH_IMAGE038
，则：
Figure 628399DEST_PATH_IMAGE039
。
S31：将系数行列式
Figure 78972DEST_PATH_IMAGE038
和向量行列式
Figure 242100DEST_PATH_IMAGE040
以联立线性方程式重新整理：联立线性方程式为：
Figure 487136DEST_PATH_IMAGE016
，其中A表示圆柱导体1的电流感应，Aii表示第i个圆柱导体1的自感，Aij表示第i个圆柱导体1对第j个圆柱导体1的互感，比如：A11表示第一个圆柱导体1的自感，A12表示第一个圆柱导体1对第二个圆柱导体1的互感，以此类推。Xi表示第i个圆柱导体1产生的散射波强度，X1表示第一个圆柱导体1产生的散射波强度，以此类推。Bi表示第i个圆柱导体1上的电流分布，比如B1表示第一个圆柱导体1上的电流分布，以此类推。
S32：对权重函数使用全域基函数方式，将系数行列式对称化，则对角行列式
Figure 1294DEST_PATH_IMAGE038
可表示为
Figure 724400DEST_PATH_IMAGE041
，行列式Bu则表示为
Figure 742034DEST_PATH_IMAGE042
。
通过矩阵变换得出如下散射波强度的公式：
Figure 423551DEST_PATH_IMAGE043
，其中，
Figure 159426DEST_PATH_IMAGE044
通过上述公式可见，最终天线设计的主要问题转变为对向量
Figure 358326DEST_PATH_IMAGE045
的计算问题。即波导组件中圆柱导体1的直径和相邻圆柱导体1之间间隔的问题。
下表示意出上述基板集成波导天线实施例的电磁场的初始相关参数值：
Figure 620681DEST_PATH_IMAGE047
在该实施例中，经计算和数值模拟，在如下关系下电磁波聚集效果更好，电磁波不会漏出：令圆柱导体1的直径为D，沿平行波导通道4的方向上，相邻两个圆柱导体1的中心距为P，则D/P=0.5。设λ为输入电磁波的波长，沿垂直波导通道4的方向上，相邻两个圆柱导体1的中心距为D1，则D1=0.9λ。
作为优选的，同轴电缆3至介质基板2闭口处的距离为D2，则D2=0.25λ。
更进一步的，波导组件的波导通道4长度为D3，则D3=2.4λ。
图2所示为本实施例所述的基板集成波导天线结构下的电场分布图（横轴为波导通道4的宽度方向距离，左侧纵轴为波导通道的长度方向距离），图中的高亮区域即为波导通道4，波导通道4的左右两侧为圆柱导体1，从图中可以看出，电磁场均位于波导通道4内，没有穿过相邻圆柱导体1向外泄露，因此实现了电磁波的高效率传播。
本发明所述的基板集成波导天线结构简单，同时采用通过圆柱导体1的直径D和波导通道4方向相邻圆柱导体1的中心距P的合理配比，可以使天线结构避免漏波现象发生，达到最优性价比的设计性能。
参数目录：
L：线性算子，
g：已知函数，
h：未知函数；
an：未知系数；
hn：未知函数h的基底函数；
Wm：权重函数；
t：时间；
ω：震荡角频率；
f：频率；
φ0：入射角度，
k0：传播常数，
μ0：真空中的透磁率，
ε0：真空中的介电常数，
λ：输入电磁波的波长；
ρ’：从柱坐标原点到圆柱导体1圆周上任一点的向量
ρv：从柱坐标原点到第v号导体圆柱中心点的向量
E0：入射波电场强度；
Figure 348465DEST_PATH_IMAGE048
：Z轴方向上任意点
Figure 165111DEST_PATH_IMAGE049
的入射波电场强度，
η：传播导纳；
Figure 636544DEST_PATH_IMAGE050
：Z方向上电场在任意点
Figure 753405DEST_PATH_IMAGE051
的感应电流，
Figure 386511DEST_PATH_IMAGE052
：Z方向上电场在第v号的圆柱导体1圆周表面上ρ’位置的感应电流；
Figure 221612DEST_PATH_IMAGE053
：需求解的第v个圆柱导体1上的第n个划分网格点的行列式；
Figure 10000254965088
：媒介的电导率；
Figure 27894DEST_PATH_IMAGE054
：Z轴方向上任意点
Figure 874627DEST_PATH_IMAGE055
的散射电场强度；
Figure 92286DEST_PATH_IMAGE056
点的散乱波电场强度；
Jm(k0au)：第m个网格区域的感应电流；
av：第v号圆柱导体1的半径，
au：第u号圆柱导体1的半径；
Figure 24470DEST_PATH_IMAGE057
：系数行列式；
Figure 696759DEST_PATH_IMAGE058
：向量行列式； A：圆柱导体1的电流感应，A11：第一个圆柱导体1的自感，A12：第一个圆柱导体1对第二个圆柱导体1的互感，
X1：第一个圆柱导体1产生的散射波强度；
B1：第一个圆柱导体1上的电流分布；
D1：沿垂直波导通道4的方向上，相邻两个圆柱导体1的中心距；
D：圆柱导体1的直径；
P：沿平行波导通道4的方向上，相邻两个圆柱导体1的中心距；
D2：同轴电缆3至介质基板2闭口处的距离；
D3：波导通道4的长度。
在本说明书中，对所述术语的示意性表述不一定指的是相同的实施例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或多个实施例中以合适的方式结合。
以上述依据本发明的理想实施例为启示，通过上述的说明内容，相关工作人员完全可以在不偏离本项发明技术思想的范围内，进行多样的变更以及修改。本项发明的技术性范围并不局限于说明书上的内容，必须要根据权利要求范围来确定其技术性范围。

Claims (10)

1.一种基板集成波导天线的设计方法，其特征在于：所述方法采用如下的基板集成波导天线，所述基板集成波导天线包括由若干阵列布置的圆柱导体（1）组成的波导组件和连接于波导组件轴向两端的介质基板（2），其中介质基板（2）连接一同轴电缆（3）；所述方法包括以下步骤：
S1：计算散射波的基本解；
S2：利用多种边界条件将S1中散射波的基本解进一步变形；
S3：采用全域基函数的矩量法，将电磁场水平和垂直方向的计算转化为关于柱坐标下水平方向的计算问题。
2.根据权利要求1所述的基板集成波导天线的设计方法，其特征在于：步骤S1中散射波的基本解的计算过程如下：
S11：列出入射波的电场强度公式：
Figure 838978DEST_PATH_IMAGE001
，其中，
Figure 813887DEST_PATH_IMAGE002
（
Figure 764526DEST_PATH_IMAGE003
为时间因子，φ0为入射角度，k0为传播常数，μ0为真空中的透磁率，ε0为真空中的介电常数，ρ为圆柱导体1上任一点的柱面坐标，E0为入射波的电场强度）；
S12：由圆柱导体（1）上的电磁感应电流产生的散射波满足下述的亥姆霍兹波动方程式：
Figure 240507DEST_PATH_IMAGE004
，其中，
Figure 61832DEST_PATH_IMAGE005
（η为传播导纳，
Figure 332276DEST_PATH_IMAGE006
为圆柱导体（1）上的电磁感应电流，
Figure 10000154905051
为媒介的电导率）；
S13：将S12中的方程式转化为如下微分方程式：
Figure 239053DEST_PATH_IMAGE007
，并求其基本解为：
Figure 518724DEST_PATH_IMAGE008
；
S14：将S12中的参数代入S13的基本解中，得出散射波的基本解为：
Figure 194556DEST_PATH_IMAGE009
，其中，C=C1+C2+…+CN。
3.根据权利要求2所述的基板集成波导天线的设计方法，其特征在于：所述步骤S2包括：
S21：根据圆柱导体（1）上的全电场等于零得出方程式：
Figure 901481DEST_PATH_IMAGE010
；并将该方程式代入S14的基本解中，并得出基本解的第一转化形式：
Figure 295553DEST_PATH_IMAGE011
，其中，
Figure 113337DEST_PATH_IMAGE012
；
S22：将圆柱导体（1）表面的点用极坐标轴表示，则圆柱导体（1）上任一点的坐标函数为
Figure 174833DEST_PATH_IMAGE013
；将该坐标函数代入步骤S21中基本解的第一转化形式，得出基本解的第二转化形式：
Figure 52660DEST_PATH_IMAGE015
；
S23：将步骤S22中基本解的第二转化形式积分，得出基本解的第三转化形式：
Figure 668449DEST_PATH_IMAGE016
（av表示第v号圆柱导体（1）的半径，au表示第u号圆柱导体（1）的半径）。
4.根据权利要求3所述的基板集成波导天线的设计方法，其特征在于：所述步骤S3包括：
S31：将步骤S23中基本解的第三转化形式转化为如下的联立线性方程式：
Figure 289923DEST_PATH_IMAGE017
，其中A表示圆柱导体（1）的电流感应，Aii表示第i个圆柱导体（1）的自感，Aij表示第i个圆柱导体（1）对第j个圆柱导体（1）的互感，Xi表示第i个圆柱导体（1）产生的散射波强度，Bi表示第i个圆柱导体（1）上的电流分布；
S32：将系数行列式对称化，通过矩阵变换得出如下散射波强度的公式：
Figure 205926DEST_PATH_IMAGE019
，其中，
Figure DEST_PATH_IMAGE020
。
5.根据权利要求1所述的基板集成波导天线的设计方法，其特征在于：令所述圆柱导体（1）的直径为D，沿平行波导通道（4）的方向上，相邻两个圆柱导体（1）的中心距为P，则D/P=0.5。
6.根据权利要求1所述的基板集成波导天线的设计方法，其特征在于：所述圆柱导体（1）为金属管。
7.根据权利要求5所述的基板集成波导天线的设计方法，其特征在于：所述同轴电缆（3）设置于所述波导组件的中心线上。
8.根据权利要求5所述的基板集成波导天线的设计方法，其特征在于：输入电磁波的波长为λ；沿垂直波导通道（4）的方向上，相邻两个圆柱导体（1）的中心距为D1，则满足如下关系：D1=0.9λ。
9.根据权利要求5所述的基板集成波导天线的设计方法，其特征在于：所述同轴电缆（3）至介质基板（2）闭口处的距离D2与输入电磁波的波长λ的关系为D2=0.25λ。
10.根据权利要求5所述的基板集成波导天线的设计方法，其特征在于：所述波导组件的波导通道（4）长度D3与输入电磁波的波长λ的关系为D3=2.4λ。
CN202210053123.6A 2022-01-18 2022-01-18 一种基板集成波导天线的设计方法 Active CN114065559B (zh)

Priority Applications (4)

Application Number Priority Date Filing Date Title
CN202210053123.6A CN114065559B (zh) 2022-01-18 2022-01-18 一种基板集成波导天线的设计方法
JP2023517786A JP2024504894A (ja) 2022-01-18 2022-10-19 基板一体型導波管アンテナの設計方法
EP22865879.5A EP4239514A4 (en) 2022-01-18 2022-10-19 METHOD FOR DESIGNING A SUBSTRATE-INTEGRATED WAVEGUIDE ANTENNA
PCT/CN2022/126192 WO2023138127A1 (zh) 2022-01-18 2022-10-19 一种基板集成波导天线的设计方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210053123.6A CN114065559B (zh) 2022-01-18 2022-01-18 一种基板集成波导天线的设计方法

Publications (2)

Publication Number Publication Date
CN114065559A true CN114065559A (zh) 2022-02-18
CN114065559B CN114065559B (zh) 2022-04-15

Family

ID=80231266

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210053123.6A Active CN114065559B (zh) 2022-01-18 2022-01-18 一种基板集成波导天线的设计方法

Country Status (4)

Country Link
EP (1) EP4239514A4 (zh)
JP (1) JP2024504894A (zh)
CN (1) CN114065559B (zh)
WO (1) WO2023138127A1 (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023138127A1 (zh) * 2022-01-18 2023-07-27 常州星宇车灯股份有限公司 一种基板集成波导天线的设计方法

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118350255B (zh) * 2024-06-17 2024-08-30 西安中电科西电科大雷达技术协同创新研究院有限公司 一种基于有限元的圆波导器件电磁特性的计算方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103226644A (zh) * 2013-05-17 2013-07-31 南京理工大学 基于柱面等效源区域分解的电磁散射特性仿真方法
CN108763153A (zh) * 2018-04-11 2018-11-06 杭州电子科技大学 一种地下各向同性介质球电磁散射的计算方法
CN111668598A (zh) * 2019-03-06 2020-09-15 安波福技术有限公司 具有寄生特征的狭槽阵列天线
CN112287588A (zh) * 2020-11-09 2021-01-29 南京理工大学 一种基于特征模理论的天线阵列快速分析方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US10468738B2 (en) * 2016-03-15 2019-11-05 Aptiv Technologies Limited Signaling device including a substrate integrated waveguide coupled to a signal generator through a ball grid array
US11139581B2 (en) * 2019-03-07 2021-10-05 Aptiv Technologies Limited Signaling device including a slot transition between a substrate integrated waveguide and a signal generator
CN112688046B (zh) * 2020-12-04 2022-03-29 华南理工大学 一种近场聚焦全息阵列天线及调控方法
CN114065559B (zh) * 2022-01-18 2022-04-15 常州星宇车灯股份有限公司 一种基板集成波导天线的设计方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103226644A (zh) * 2013-05-17 2013-07-31 南京理工大学 基于柱面等效源区域分解的电磁散射特性仿真方法
CN108763153A (zh) * 2018-04-11 2018-11-06 杭州电子科技大学 一种地下各向同性介质球电磁散射的计算方法
CN111668598A (zh) * 2019-03-06 2020-09-15 安波福技术有限公司 具有寄生特征的狭槽阵列天线
CN112287588A (zh) * 2020-11-09 2021-01-29 南京理工大学 一种基于特征模理论的天线阵列快速分析方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
何小祥等: "大规模柱面频率选择表面电磁特性分析", 《南京理工大学学报》 *
王元源: "圆柱共形微带结构的快速算法研究", 《万方学位论文全文数据库》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2023138127A1 (zh) * 2022-01-18 2023-07-27 常州星宇车灯股份有限公司 一种基板集成波导天线的设计方法

Also Published As

Publication number Publication date
EP4239514A4 (en) 2024-06-19
EP4239514A1 (en) 2023-09-06
WO2023138127A1 (zh) 2023-07-27
JP2024504894A (ja) 2024-02-02
CN114065559B (zh) 2022-04-15

Similar Documents

Publication Publication Date Title
CN114065559B (zh) 一种基板集成波导天线的设计方法
Wang et al. Theory and analysis of leaky coaxial cables with periodic slots
CN113113278A (zh) 一种类梯形交错双栅慢波结构
Persson et al. Calculating the mutual coupling between apertures on a convex circular cylinder using a hybrid UTD-MoM method
CN113964549B (zh) 一种基于干扰对消的空间取样天线的设计方法及装置
Liu et al. Crosstalk reduction design and analysis of the planar meander transmission lines
Liu et al. Analysis of a quadruple corner-cut ridged/vane-loaded circular waveguide using scaled boundary finite element method
CN104636551B (zh) 一种网状反射面天线金属丝网等效电磁参数推演方法
Dong et al. Hybrid ray tracing method for microwave lens simulation
CN110323522A (zh) 一种基于H-T接头功分网络的TE10-TEn0的模式变换器
CN116151038A (zh) 一种圆极化mimo微带天线阵列自解耦技术的分析方法
Yu et al. Analysis of finite and curved frequency‐selective surfaces using the hybrid volume‐surface integral equation approach
Villegas et al. Leakage fields from planar semi-infinite transmission lines
CN111143991B (zh) 一种介质包裹导线的横磁波传输模型及其构建方法
Douglas Design and Analysis of microstrip antenna for 2.4 GHz applications
Addamo et al. Bloch wave analysis of long leaky coaxial cables
JP3978503B1 (ja) ストリップ線路型左手系線路でなるアンテナ
Lue et al. The equivalent parameters for the radiating slot on a sectoral waveguide
Berdnik et al. Waveguide T-junctions with resonant coupling between sections of different dimensions
Amendola et al. Numerical and analytical characteristic modes for conducting elliptic cylinders
Hongo et al. Field distribution in a flanged parallel-plate waveguide
CN109167143B (zh) 一种双层匹配隔板波导功分器结构及优化方法
da Silva et al. Analysis of the junction between smooth and corrugated cylindrical waveguides in mode converters
Sarnowski et al. Characterization of diffraction anomalies in 2-D photonic bandgap structures
Yang et al. Analysis of strip-grated conducting rod using asymptotic strip boundary conditions

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant