《基于时域谱元法的微波无源电路电磁热一体化分析方法.pdf》由会员分享，可在线阅读，更多相关《基于时域谱元法的微波无源电路电磁热一体化分析方法.pdf（9页珍藏版）》请在专利查询网上搜索。
1、10申请公布号CN104050307A43申请公布日20140917CN104050307A21申请号201310718765422申请日20131223G06F17/5020060171申请人南京理工大学地址210094江苏省南京市孝陵卫200号72发明人陈如山丁大志樊振宏盛亦军陈桂莲74专利代理机构南京理工大学专利中心32203代理人朱显国54发明名称基于时域谱元法的微波无源电路电磁热一体化分析方法57摘要本发明公开了一种基于谱元法的微波无源电路电磁热一体化分析方法该方法首先采用时域谱元法求解麦克斯韦方程组，求出电路在高功率脉冲作用下瞬时的电场和磁场分布，得出当前时刻的电磁损耗。假设模型内。
2、部电磁损耗全部转化为热能，将所得热能带入热传导方程，得到当前时刻各点的温度分布情况。利用介电参数随温度变化的关系式得到下一时刻材料的电特性参数并再次计算电磁场方程，得出电磁损耗。如此反复循环，直到完成预定加热时间。电磁热一体化分析可以清楚的得到滤波器在不同脉冲的作用下，模型内部温度随时间变化的分布情况，建模灵活，剖分方便，形成的矩阵具有良好的稀疏性，求解效率较高。51INTCL权利要求书1页说明书6页附图1页19中华人民共和国国家知识产权局12发明专利申请权利要求书1页说明书6页附图1页10申请公布号CN104050307ACN104050307A1/1页21一种基于时域谱元法的微波无源电路电。
3、磁热一体化分析方法，其特征在于步骤如下第一步，建立微波无源电路的求解模型，并采用曲六面体对模型进行剖分，得到模型的结构信息，包括六面体的单元信息及节点信息；第二步，从麦克斯韦方程方程组出发，对等式采用伽辽金法测试，强加电场边界条件，得到电场波动方程，进而求解得到各节点TN时刻的电场及电流分布；第三步，由电场及电流分布得出各节点电磁损耗，包括介质损耗及导体损耗；第四步，建立微波无源电路的热传导方程，将电磁损耗作为热源项代入该方程中，求解得到各节点温度分布；第五步，判断TN时刻是否为预定加热时刻，若是预定加热时间则停止计算，若不是预定加热时间，则更新电场波动方程中的介电常数和电导率，计算TN1的电。
4、场分布及电流分布。2根据权利要求1所述的基于时域谱元法的微波无源电路电磁热一体化分析方法，其特征在于第一步中，建立微波无源电路求解模型时忽略金属贴片厚度，第二步数值分析计算时，在原本模型金属贴片位置处强加金属边界条件。3根据权利要求1所述的基于时域谱元法的微波无源电路电磁热一体化分析的数值方法，其特征在于第三步中，计算微波电路的介质损耗和金属带上的导体损耗由计算得到各节点介质损耗，其中2F，为复介电常数的虚部，F为信号源中心频率，为电导率，由得到各节点的导体损耗，RS为导体的表面电阻，为表面电流。权利要求书CN104050307A1/6页3基于时域谱元法的微波无源电路电磁热一体化分析方法技术领。
5、域0001本发明属于微波无源器件的电磁热一体化分析，特别是针对SIW滤波器设计的数值分析方法。背景技术0002上世纪70年代以来，由于计算机技术的飞速发展，计算电磁学获得了迅猛发展。相对于解析方法，数值方法能够解决很大一类计算巨大，结构复杂而解析方法难以或无法得到精确结果的问题。针对于微波无源电路，如滤波器、微带线的热效应作用过程，可以单独考虑为电磁求解过程和热模型求解过程。即分别对应于麦克斯韦电场波动方程和热传导方程的求解。研究微波热模型就是选择合适的数值分析方法求解分析这两大过程。对于电磁模型数值求解，目前主流方法按时频域分为时域方法如时域有限差分法和频域方法如有限元，矩量法等等。在应用上。
6、述方法分析微波热模型中的电磁问题时，有限元，矩量法，时域有限差分法FDTD等当前热门数值算法都可以使用，但考虑到微波热模型中电参数为时变函数和热模型时间延续性的特点，采用时域方法更为合适。一般FDTD，FEM等方法更为普遍。然而由于FDTD的YEE网格特性决定了其无法模拟结构复杂的模型。FEM应用到时域时每个时间步都涉及到对线性方程组的求解，计算量非常庞大，很浪费时间。发明内容0003本发明的目的在于提供一种微波无源电路的电磁热一体化分析，同时考虑介质损耗和导体损耗（忽略辐射损耗）从而实现快速得到器件内部温度分布的方法。0004实现本发明目的的技术解决方案为一种微波无源器件电热一体化分析方法，。
7、步骤表述如下0005第一步，建立求解模型及网格剖分。并采用曲六面体对模型进行整体剖分，得到模型的结构信息，包括每个曲六面体单元的结点编号和坐标等。剖分网格的尺寸大于满足精度所需的剖分尺寸。0006第二步，从麦克斯韦方程方程组出发，对等式采用伽辽金法测试，强加边界条件，得到电场波动方程，进而求解得到各节点TN时刻的电场及电流分布。0007第三步，由以上步骤得到的电场及电流分布得出各点的导体损耗及介质损耗（忽略辐射损耗）。在求解介质损耗时用等效电导率代替，即2F。其中，为复介电常数的虚部，F为工作频率。然后通过公式计算得到单位体积内介质损耗所转化的热能。求解导体损耗的时候，首先通过麦克斯韦方程由磁。
8、场分布计算得到电流密度分布。由于良导体中电流分布主要集中在导体表面。假设微波电路所敷金属贴片厚度即为电磁波进入金属的趋肤深度，计算得到表面电阻。由得到导体损耗。0008第四步，建立微波无源电路的热传导方程，将电磁损耗作为热源项代入该方程中，说明书CN104050307A2/6页4求解得到各节点温度分布；0009第五步，判断TN时刻是否为预定加热时刻，若是预定加热时间则停止计算，若不是是预定加热时间则更新电场波动方程中的介电常数和电导率，计算TN1的电场分布及电流密度分布。重复步骤二、三、四，五步，如此反复循环，直到完成预定加热时间；0010本发明与现有技术相比，其显著优点（1）SETD采用曲六。
9、面体剖分，建模灵活，剖分方便，使用用特定的正交多项式作为基函数，随着多项式阶数的提高，计算误差将呈指数下降。（2）对于电磁场场计算空间，得到的质量矩阵为块对角阵；可用块对角矩阵的求逆方法事先求出质量矩阵的逆，使整个方程的求解变为显式，降低计算量（3）微波器件的电特性与热特性是一体化分析的，没有将其割裂开来。而且下一步可以将电热之间的相互影响联系起来。（4）热传导方程形成的矩阵方程性态较好，直接求逆相当方便，使用直接解法可快速求得器件内部的温度分布。附图说明0011图1是SIW带通滤波器的模型图。0012图2是本发明计算得到带通滤波器在通过高斯调制脉冲时温度分布图。具体实施方式0013下面结合附。
10、图对本发明作进一步详细描述。0014一、传导方程简单推导0015在热媒质内X,Y,Z点，T时刻热量传输满足导热偏微分方程，即热传导方程HTE，方程的形式为00160017式中，T为物体的瞬态温度；TS为过程进行的时间；KTW/M为材料的导热系数；MKG/M3为材料的密度；CMJ/KG为材料的定压比热；VSW/M3为冷却流的热容流积；TA为冷却流的温度；PDW/M3为内部热源的功率密度。0018为简化问题，假定M、CM和KT均不是时间、空间和温度的函数，且不考虑冷却流，在直角坐标系下，式1改写为00190020上式中00210022DT定义为热扩散率。需要说明的是，导热偏微分方程是描述导热过程共。
11、性的通用表达式，适用于任何导热过程。其中热源项PD我们仅考虑介质损耗和导体损耗，忽略辐射损耗，将在第三节详细说明。0023传导方程的差分形式存在数值稳定性条件。为了防止随着迭代次数的增加而导致场值结果的溢出，在求解热模型时，时间步长和物体单元的网格大小必须要满足一定的准说明书CN104050307A3/6页5则。具体表达式为00240025需要指出的是，当同时求解麦克斯韦方程和热传导方程时，必须选择最小时间步的值。0026在计算热方程时，需要边界上的温度场值，同样可以采用类似电磁场中MUR吸收边界的方法，利用边界附近网格的温度场去构建边界面上的温度场值。为了得到固体导热偏微分方程的唯一解，必须。
12、附加边界条件和初始条件，统称为定解条件，这里介绍初始条件和三类边界条件。0027初始条件已知是过程开始时物体整个区域中所具有的温度为已知值，用公式表示为0028T|T0T0150029式中，T0为一已知常数，表示物体初温是均匀的。0030第一类边界条件是指物体边界上的温度函数为已知，用公式表示为0031T|TW160032式中，为物体边界，其方向是逆时针方向；TW为已知壁面的温度（常数）。0033第二类边界条件是指物体边界上的热流密度为已知，用公式表示为00340035式中，Q0W/M2为已知热流密度（常数）；当Q00时，称为绝热边界条件。绝热边界条件本身不能唯一确定温度场，需要别的边界条件或。
13、初始条件配合下才能确定温度场。0036第三类边界条件是指与物体相接触的环境温度TF和热对流系数H为已知，用公式表示为00370038式中，H与TF可以是常数，也可以是某种随时间和位置而变化的函数，如果H和TF不是常数，则在数值计算中经常分段取其平均值作为常数。0039二、时域谱元法的简单介绍0040谱元法结合了有限元法和谱方法，求解过程与有限元法类似，其主要步骤是1首先把计算的区域分成许多子域（单元），每个子域由若干节点组成；2在每个子域上把近似解表示成截断的切比雪夫正交多项式或勒让德多项式展开；3用伽辽金法对波动方程进行测试从而形成矩阵方程；4求解方程。它的基本思想是在每一个单元上GAUSS。
14、LOBATTOLEGENDREGLL基函数离散。在一维标准参考单元1,1中，N阶GLL基函数为00410042其中，J0,1,N，LN是N阶勒让德多项式，LN是其导数。GLL积分点说明书CN104050307A4/6页6的选取为方程式1J2LNJ0的根，它们均定义在式1,1内。因而温度T可以表示成0043利用变分原理和GALERKIN法，将热传导方程在离散空间中转换成其残量形式0044RINIRD0230045其中NI为加权函数，残量R为00460047将22式和24式代入23式并利用矢量恒等式及后得到00480049在不同媒质交界面上，温度场满足一类特定的关系。以第三类边界条件为例，温度场满。
15、足以下关系00500051其中，H为对流传热系数，一般取H10W/M2；N为边界面的单位法向矢量，TAIR为媒质外部温度即环境温度，TSURF为媒质的表面温度。将26式代入25式得到，00520053其紧凑格式为00540055其中00560057TIJNINJDV210说明书CN104050307A5/6页70058005900600061采用前向差分格式得0062TTN1TTSRTNTRQTF2140063其中，T为质量矩阵，S为刚度矩阵，R边界积分矩阵，F为源矢量矩阵，RQ由边界条件所引起的边界源矢量矩阵。由于单元基函数的正交特性，我们所得到的质量矩阵为对角阵或块对角阵，这样我们可利用块。
16、对角矩阵求逆的方法事先求出质量矩阵的逆，使得被求解方程为显式方程，降低了计算量，提高了计算效率。这是谱元法的一大重要特性。0064三、热源项分析0065微带线的辐射损耗相对于介质损耗和导体损耗比较小决定了辐射损耗可以忽略不计。在SETD公式推导中，当只记入介质损耗时，此时用等效电导率代替0066310067其中，为电导率，为复介电常数的虚部；、分别对应待求电场值处的复介电常数实部和等效电导率。在求解电磁问题时，假定介质中的场分布经过激励源加载后的短暂瞬时过程后，达到一个相对稳定的分布。一旦得到电场值分布，据此就可以计算出媒质损耗电磁能以及转化到热模型中的热能。每个单元处的单位体积损耗能转化为热。
17、能功率的计算式如下00680069上式中，是电场E的均方根值，R是复相对介电系数的虚部。0070在SETD公式推导中，只记入导体损耗时，由麦克斯韦第二方程求解得到介质基片表面磁场分布，再由麦克斯韦第一方程求解得到电流然而电磁波在良导体中的衰减很快，在传播很短的一段距离之后就几乎衰减完了。因此，良导体中的电磁波局限于导体表面附近的区域，即导体的趋肤效应。导体的趋肤深度非常小，以至于可以认为电流仅分布于导体表面很薄的一层内，也即介质基片表面一层。这与恒定电流或低频电流均匀分布于导体横截面上的情况不同，在高频时，导体的实际载流截面减小了，因而导体的高频电阻大于直流或低频电阻，它与电导率和趋肤深度有关。
18、，可表示成即厚度为的导体每平方米的电阻，称为导体的表面电阻率，简称表面电阻。最后由公式得出导体损耗。0071四、电磁热耦合的基本理论0072数值方法模拟微波加热问题的总体思想是首先，在TN时刻根据初始的电特性条件，输入一定的微波功率，计算出各点的微波耗散功率，在模拟仿真时，假设电磁损耗能全说明书CN104050307A6/6页8部转化为热能。然后，将得到的热能代入到热传导方程中，得到各点的温度。根据物体各点处不同的温度分布，利用介电参数随温度变化的关系式得到下一时刻TN1的介质材料的电特性参数并再次重复计算电磁场方程，计算耗散功率，如此反复循环，直到完成预定的加热时间。在热场和电磁场的求解过程。
19、中采用相同的剖分网格单元。0073以图1真分析的带通滤波器为例带通滤波器的通带范围为12GHZ125GHZ。介质基板采用ROGERSRT5880，板厚0508MM，具体尺寸由于涉密不方便给出。0074在进行热模拟仿真时需要具体知道微带电路的各部分性能参数电导率，相对介电常数等。一般情况下，这些性能参数都会受到随周围的环境温度，压强的影响，这些参数都比较难以精确的获得，这也是数值模拟过程和实际升温过程存在差异的重要原因。0075表1，表2列出了一般情况下ROGERSRT5880介质基板的性能参数，及覆铜片的的性能参数00760077表1介质基板的性能参数00780079表2铜的性能参数0080将周围环境温度设为20，由于微波电路本身接触空气，故将周围表面设为第三类边界条件即为热对流形式。在微波电路的输入端注入幅度为10MV/M的调制高斯脉冲，中心频率设为125GHZ，脉冲相关量设为9。图2为脉冲通过微波电路003US时的瞬态热分布图，可以清楚的看到，能量聚集在输入端口，这是由于调制高斯脉冲的频谱范围覆盖比较广，只有处在通带范围内的一部分频谱分量通过微波电路，到达输出端口。说明书CN104050307A1/1页9图1图2说明书附图CN104050307A。