CN103412984B - 一种规则波导端口微波部件电磁仿真的功率馈电方法 - Google Patents

一种规则波导端口微波部件电磁仿真的功率馈电方法 Download PDF

Info

Publication number
CN103412984B
CN103412984B CN201310309258.5A CN201310309258A CN103412984B CN 103412984 B CN103412984 B CN 103412984B CN 201310309258 A CN201310309258 A CN 201310309258A CN 103412984 B CN103412984 B CN 103412984B
Authority
CN
China
Prior art keywords
waveguide
field
holotype
power
section
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201310309258.5A
Other languages
English (en)
Other versions
CN103412984A (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.)
China Academy of Space Technology CAST
Original Assignee
China Academy of Space Technology CAST
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 China Academy of Space Technology CAST filed Critical China Academy of Space Technology CAST
Priority to CN201310309258.5A priority Critical patent/CN103412984B/zh
Publication of CN103412984A publication Critical patent/CN103412984A/zh
Application granted granted Critical
Publication of CN103412984B publication Critical patent/CN103412984B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Constitution Of High-Frequency Heating (AREA)

Abstract

本发明公开了一种规则波导端口微波部件电磁仿真的功率馈电方法,根据输入端口截面的形状建立以该输入端口为横截面的无限长波导,确定在所述无限长波导中传输的主模式场;根据激励信号确定所述主模式场平均传输功率,根据所述平均传输功率和输入功率确定所述主模式场的初始电场值功率系数,根据所述主模式场的初始电场值功率系数确定主模式场分布下的初始电场值;设置连接边界区域并将其分解成多个网格单元,在沿主模式场的电场线方向的每个网格单元的棱边中点上加入该点对应的主模式场的初始电场值,实现功率馈电。

Description

一种规则波导端口微波部件电磁仿真的功率馈电方法
技术领域
本发明涉及电磁仿真技术应用领域,尤其涉及一种规则波导端口微波部件电磁仿真的功率馈电方法。
背景技术
对于微波部件而言,馈入输入信号是进行微波部件产品设计与电磁仿真的关键性步骤。而规则波导是微波部件最为常见的输入输出端口形式之一。实际工程中测试信号常以功率作为单位,当结合电磁仿真方法与粒子模拟方法进行大功率微波部件等离子体仿真时,无法直接获得输入功率对应下的输出结果,会导致数值仿真结果无法直接进行实验验证。显然,为使电磁仿真方法更好地与其他仿真方法进行结合,满足更广泛的工程要求,尤其是大功率条件下的仿真要求,必须研究一种能够针对规则波导端口微波部件在其电磁仿真中实现功率馈电的方法,提供与工程相符的激励形式。现有技术对规则波导端口微波部件电磁仿真的输入信号馈入方法进行研究时,均通过假定规则波导端口输入截面的初始电磁场分布已知进行,都未涉及如何进行功率馈电的问题。
发明内容
本发明的目的在于,提供一种规则波导端口微波部件电磁仿真的功率馈电方法,该方法能够在规则波导端口微波部件电磁仿真中进行功率馈电,实现功率信号的加载,该方法易于实现,与其他仿真方法,例如粒子模拟方法兼容性高,具有广阔的应用前景。
本发明的上述目的是通过如下技术方案予以实现的:
一种规则波导端口微波部件电磁仿真的功率馈电方法,包括如下步骤:
a.建立规则波导端口微波部件的三维几何模型,并根据三维几何模型提取激励信号输入端口截面的形状参数;
b.确定规则波导端口微波部件的工作频率、激励信号和输入功率;
c.根据所述输入端口截面的形状建立以该输入端口为横截面的无限长波导,根据麦克斯韦方程组与所述无限长波导的边界条件确定所述无限长波导中波导横截面的各模式场分布和对应的截止频率;
d.根据所述各模式场分布和对应的截止频率,以及规则波导端口微波部件的工作频率确定在所述无限长波导中传输的主模式场;根据激励信号确定所述主模式场平均传输功率,根据所述平均传输功率和输入功率确定所述主模式场的初始电场值功率系数,根据所述主模式场的初始电场值功率系数确定主模式场分布下的初始电场值;
e.在所述无限长波导中与输入端口截面一定距离处设置连接边界区域,所述连接边界区域的尺寸和外形与所述输入端口截面的形状相同;将所述连接边界区域分解成多个网格单元,在沿主模式场的电场线方向的每个网格单元的棱边中点上加入该点对应的主模式场的初始电场值,从而在连接边界区域激励起主模式场。
所述步骤c中,所述无限长波导中波导横截面的各模式场分布通过分离变量法解析求解得到。
所述步骤d中,
当激励信号为时谐信号时,根据坡印亭定理得到所述主模式场平均传输功率。
所述步骤d中,
当激励信号为非时谐信号时,对非时谐信号进行傅立叶变换得到其谱域表达式,结合坡印亭定理得到所述主模式场平均传输功率。
本发明与现有技术相比具有如下有益效果:
本发明提出了一种基于电磁仿真进行微波部件功率馈电的方法,使得时域电磁算法能够应用于需要以功率作为输入单位的应用场合,实现与多种算法的无缝连接,例如粒子模拟算法或自主研发的其他算法。本方法拓展了时域电磁算法的应用领域,为多种实际应用问题的解决提供了有效途径,例如需要改变输入信号功率的微放电电磁粒子仿真平台的开发,具有广阔的应用前景。
附图说明
图1矩形波导端口激励源的引入示意图。
具体实施方式
下面结合附图对本发明具体实施方式作进一步详细的描述:
实施例一
本实施例,以输入信号为时谐信号时的矩形波导端口微波部件功率馈电为例,对本发明的规则波导端口微波部件电磁仿真的功率馈电方法进行说明,具体实施过程如下:
a.利用CAD软件建立待仿真的微波部件的三维几何模型,并根据三维几何模型提取激励信号输入端口截面1的形状参数,长度a为22.86mm,宽度b为10.16mm,如图1所示。
b.确定待仿真微波部件的工作频点f0为10GHz、激励信号为正弦波信号、输入功率为P0
c.如图1所示,根据输入端口截面1的形状建立以该输入端口为横截面的无限长波导2,建立以沿输入端口截面1的长边为x方向,宽边为y方向,沿无限长波导2的方向为z方向,以输入端口截面1的顶点为原点的直角坐标系,根据麦克斯韦方程组与所述无限长波导的边界条件,采取分离变量法获取无限长波导2中波导横截面的模式场分布:
E t ( x , y ) = E x ( x , y ) x → + E y ( x , y ) y → , - - - ( 1.1 )
其中: E x ( x , y ) = Σ m , n A 0 A mn ( f ) cos ( mπ a x ) sin ( nπ b y ) = A 0 Σ m , n E ~ x ( m , n , x , y ) , - - - ( 1.2 )
E y ( x , y ) = Σ m , n A 0 B mn ( f ) sin ( mπ a x ) cos ( nπ b y ) = A 0 Σ m , n E ~ y ( m , n , x , y ) - - - ( 1.3 )
为横向单位矢量,A0为以输入功率P0进行功率馈电时对应的初始电场值功率系数,Amn(f)、Bmn(f)为mn模式场幅度,是模式和频率的函数,m和n为0到无穷大的整数,进而有:
E t ( x , y ) = A 0 E ~ t ( x , y ) = A 0 [ x → Σ m , n E ~ x ( m , n , x , y ) + y → Σ m , n E ~ y ( m , n , x , y ) ] - - - ( 1.4 )
确定mn模式场对应的截止频率为
f mn = ( mπ a ) 2 + ( nπ b ) 2 2 π v - - - ( 1.5 )
其中:v为无限长波导2中的光速。
d.确定主模式场的初始电场值
根据步骤c所确定的模式场分布和规则波导端口微波部件的工作频率确定在所述无限长波导中传输的主模式场。根据待仿真规则波导端口微波部件的工作频点高于模式场的截止频率的要求,确定出满足上述要求的的m,n(本实施例中m=1,n=0);由此可见在所述无限长波导2中传输的主模式场为TE10模式场。
假设初始电场值功率系数A0为1时,根据激励信号利用坡印亭定理得到所述无限长波导2中TE10模式场平均传输功率
P ‾ ( f 0 ) = 1 2 | Z c ( f 0 ) | ∫ ∫ s | E ~ t ( x , y ) | 2 ds , - - - ( 1.6 )
其中Zc(f0)为频率为工作频点f0时TE10模式下矩形波导的波阻抗,s的积分范围为输入端口截面1对应区域。
根据输入功率P0与平均传输功率确定以输入功率P0进行功率馈电时对应的初始电场值功率系数结合步骤c所述mn模式场分布公式(1.1)~(1.4)得到每一个坐标点(x,y)处的TE10模式场的初始电场值Et(x,y)。
e.如图1所示,在所述无限长波导中与输入端口截面1一定距离d处(本实施例中为其中λ0=v/f0)设置连接边界区域3,所述连接边界区域3的尺寸和外形与所述输入端口截面1的形状相同,将所述连接边界区域3分解成许多小网格单元,在沿主模式电场线方向的每个小网格单元棱边中点(x1,y1)上加入该点对应的TE10模式场的初始电场值Et(x1,y1),从而在连接边界区域激励起TE10模式场。
实施例二
本实施例,以输入信号为非时谐高斯信号的矩形波导端口微波部件功率馈电为例,对本发明的规则波导端口微波部件电磁仿真的功率馈电方法进行说明,具体实施过程如下:
a.利用CAD软件建立待仿真的微波部件的三维几何模型,并根据三维几何模型提取激励信号输入端口截面1的形状参数,长度a为22.86mm,宽度b为10.16mm,如图1所示。
b.确定待仿真微波部件的工作频段为10GHz至12GHz、激励信号为非时谐高斯信号h(t)、输入功率为P0
c.如图1所示,根据输入端口截面1的形状建立以该输入端口为横截面的无限长波导2,建立以沿输入端口截面1的长边为x方向,宽边为y方向,沿无限长波导2的方向为z方向,以输入端口截面1的顶点为原点的直角坐标系,根据麦克斯韦方程组与所述无限长波导的边界条件,采取分离变量法获取无限长波导2中波导横截面的模式场分布
E t ( x , y ) = E x ( x , y ) x → + E y ( x , y ) y → , - - - ( 1.1 )
其中: E x ( x , y ) = Σ m , n A 0 A mn ( f ) cos ( mπ a x ) sin ( nπ b y ) = A 0 Σ m , n E ~ x ( m , n , x , y ) , - - - ( 1.2 )
E y ( x , y ) = Σ m , n A 0 B mn ( f ) sin ( mπ a x ) cos ( nπ b y ) = A 0 Σ m , n E ~ y ( m , n , x , y ) , - - - ( 1.3 )
为横向单位矢量,A0为以输入功率P0进行功率馈电时对应的初始电场值功率系数,Amn(f)、Bmn(f)为mn模式场幅度,是模式和频率的函数,m和n为0到无穷大的整数,进而有:
E t ( x , y ) = A 0 E ~ t ( x , y ) = A 0 [ x → Σ m , n E ~ x ( m , n , x , y ) + y → Σ m , n E ~ y ( m , n , x , y ) ] - - - ( 1.4 )
确定mn模式场对应的截止频率为
f mn = ( mπ a ) 2 + ( nπ b ) 2 2 π v , - - - ( 1.5 )
其中:v为无限长波导2中的光速。
d.确定主模式场的初始电场值
根据步骤c所确定的模式场分布和规则波导端口微波部件的工作频率确定在所述无限长波导中传输的主模式场。根据待仿真规则波导端口微波部件的工作频点高于主模式场的截止频率的要求,对应于10GHz至12GHz的工作频段确定出满足上述要求的m,n(本实施例中m=1,n=0);由此可见在所述无限长波导2中传输的主模式场为TE10模式场。
对激励信号非时谐信号进行傅立叶变换得到其谱域表达式h(f),其中f为工作频段10GHz至12GHz之间的任意一个工作频点,假设初始电场值功率系数A0为1时,根据激励信号与工作频段利用坡印亭定理得到所述无限长波导2中TE10模式场平均传输功率
P ‾ = ∫ f h ( f ) 2 | Z c ( f ) | df ∫ ∫ s | E ~ t ( x , y ) | 2 ds - - - ( 1.6 )
其中Zc(f)为频率为f时TE10模式下矩形波导的波阻抗,频率f的积分范围为工作频段10GHz至12GHz,s的积分范围为输入端口截面1对应区域。
根据输入功率P0与平均传输功率确定以输入功率P0进行功率馈电时对应的初始电场值功率系数结合步骤c所述mn模式场分布公式(1.1)~(1.4)得到每一个坐标点(x,y)处的TE10模式场的初始电场值Et(x,y)。
e.如图1所示,在所述无限长波导中与输入端口截面1一定距离d处(本实施例中为其中λ0=v/f0,f0为工作频段的中心频点,即11GHz)设置连接边界区域3,所述连接边界区域3的尺寸和外形与所述输入端口截面1的形状相同,将所述连接边界区域3分解成许多小网格单元,在沿主模式电场线方向的每个小网格单元棱边中点(x1,y1)上加入该点对应的TE10模式场的初始电场值Et(x1,y1),从而在连接边界区域激励起TE10模式场。
以上所述,仅为本发明最佳的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本领域的技术人员在本发明揭露的技术范围内,可轻易想到的变化或替换,都应涵盖在本发明的保护范围之内。
本发明说明书中未作详细描述的内容属于本领域技术人员的公知技术。

Claims (1)

1.一种规则波导端口微波部件电磁仿真的功率馈电方法,其特征在于,包括如下步骤:
a.建立规则波导端口微波部件的三维几何模型,并根据三维几何模型提取激励信号输入端口截面的形状参数;所述规则波导为矩形波导;
b.确定规则波导端口微波部件的工作频段、激励信号和输入功率P0;激励信号为非时谐信号h(t);
c.根据所述输入端口截面的形状建立以该输入端口为横截面的无限长波导,根据麦克斯韦方程组与所述无限长波导的边界条件确定所述无限长波导中波导横截面的各模式场分布和对应的截止频率;所述无限长波导中波导横截面的各模式场分布通过分离变量法解析求解得到;
d.根据所述各模式场分布和对应的截止频率,以及规则波导端口微波部件的工作频率确定在所述无限长波导中传输的主模式场;假设初始电场值功率系数A0为1时,根据激励信号确定所述主模式场平均传输功率根据所述主模式场平均传输功率和输入功率P0确定所述主模式场的初始电场值功率系数A0,其中对非时谐信号进行傅立叶变换得到其谱域表达式,结合坡印亭定理得到所述主模式场平均传输功率
e.在所述无限长波导中与输入端口截面一定距离处d设置连接边界区域,f0为工作频段的中心频点,υ为无限长波导中的光速;所述连接边界区域的尺寸和外形与所述输入端口截面的形状相同;将所述连接边界区域分解成多个网格单元,在沿主模式场的电场线方向的每个网格单元的棱边中点上加入该点对应的主模式场的初始电场值,从而在连接边界区域激励起主模式场。
CN201310309258.5A 2013-07-22 2013-07-22 一种规则波导端口微波部件电磁仿真的功率馈电方法 Active CN103412984B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310309258.5A CN103412984B (zh) 2013-07-22 2013-07-22 一种规则波导端口微波部件电磁仿真的功率馈电方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310309258.5A CN103412984B (zh) 2013-07-22 2013-07-22 一种规则波导端口微波部件电磁仿真的功率馈电方法

Publications (2)

Publication Number Publication Date
CN103412984A CN103412984A (zh) 2013-11-27
CN103412984B true CN103412984B (zh) 2014-10-08

Family

ID=49605996

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310309258.5A Active CN103412984B (zh) 2013-07-22 2013-07-22 一种规则波导端口微波部件电磁仿真的功率馈电方法

Country Status (1)

Country Link
CN (1) CN103412984B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113283135B (zh) * 2021-05-17 2023-02-21 核工业西南物理研究院 一种托卡马克真空室电磁仿真等离子体电流数据提取方法

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102655258A (zh) * 2011-03-01 2012-09-05 陈艳丽 输出端口隔离的h面波导t型结
CN102930071A (zh) * 2012-08-29 2013-02-13 电子科技大学 基于非匹配网格的周期结构的三维电磁场仿真模拟方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102655258A (zh) * 2011-03-01 2012-09-05 陈艳丽 输出端口隔离的h面波导t型结
CN102930071A (zh) * 2012-08-29 2013-02-13 电子科技大学 基于非匹配网格的周期结构的三维电磁场仿真模拟方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
《磁诱导三维拓扑绝缘体波导中的导模》;徐谊等;《光子学报》;20130531;第42卷(第5期);第565页到568页 *
徐谊等.《磁诱导三维拓扑绝缘体波导中的导模》.《光子学报》.2013,第42卷(第5期),第565页到568页.

Also Published As

Publication number Publication date
CN103412984A (zh) 2013-11-27

Similar Documents

Publication Publication Date Title
CN102930071B (zh) 基于非匹配网格的周期结构的三维电磁场仿真模拟方法
CN101944145B (zh) 可去除伪直流模式的微波管高频电路有限元仿真方法
Wazwaz New (3+ 1)-dimensional nonlinear evolution equation: multiple soliton solutions
Uçar et al. On q-Laplace type integral operators and their applications
Mohanty et al. Dynamics of exact closed-form solutions to the Schamel Burgers and Schamel equations with constant coefficients using a novel analytical approach
CN103412984B (zh) 一种规则波导端口微波部件电磁仿真的功率馈电方法
Kumar et al. Topological solitons and Lie symmetry analysis for the Kadomtsev-Petviashvili-Burgers equation with power law nonlinearity in dusty plasmas
Maksimenko et al. On the theory of longitudinally inhomogeneous waveguides with impedance walls
Wang et al. Implicit DG method for time domain Maxwell’s equations involving metamaterials
CN107609233A (zh) 一种行波管通用注波互作用模型的不连续场匹配方法
Asfar et al. Circular waveguide with sinusoidally perturbed walls
Seviour et al. Effects of high power on microwave metamaterials
Geng et al. Quasi-periodic solutions of the discrete mKdV hierarchy
CN107239586A (zh) 对无条件稳定时域有限差分法有效的区域分解并行方法
Kravchenko et al. Nonlinear theory of relativistic Cherenkov generators based on irregular waveguides with allowance for finite wall conductance.
Bai et al. Corrugated waveguide mode purifier for TEM output in a dual-mode operation overmoded coaxial millimeter-wave generator
CN108536929A (zh) 一种应用arpack求解波导结构色散特性的方法
Altınbasak et al. 2D Quasistatic TLM Field Solver for High Speed PCB Design
Xie et al. Design of Quasi-optical Mode Converter for 800-GHz TE 8, 5 Mode Gyrotron
Liu et al. Accelerated computing FDTD algorithm based on lossy dispersion medium in low frequency
Han et al. A study of waveguide loaded with periodical dielectrics
Zhang et al. The optimization of the parabolic reflector antenna layout with genetic algorithm and geometrical theory of diffraction
Mesa et al. On the Computation of the Dispersion Diagram of Symmetric Periodic One-dimensional Structures
Yu et al. Impedance matching design for capacitively coupled plasmas considering coaxial cables
Zharkova et al. Stages of particle acceleration in a 3D reconnecting current sheet

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant