电源地板系统中的同步开关噪声仿真.docVIP

电子科技大学物理电子学院 标 准 实 验 报 告 （实验）课程名称 电磁仿真综合实践 电子科技大学教务处制表 电 子 科 技 大 学 实 验 报 告 学生姓名：胡仕虹 学 号：2804203021 指导教师： 王秉中、王晓华 实验地点： 实验时间：2011年5月—6月 一、实验室名称：物理电子学院仿真实验室 二、实验项目名称：电磁仿真综合实践I级 三、实验学时：36学时 四、实验原理：电磁仿真的时域有限差分法 五、实验目的：加强学生的计算机综合应用能力、尤其是运用计算机分析和解决专业问题的能力培养，使学生对独立进行科学研究有初步实践；初步掌握一种纯数值电磁仿真方法——时域有限差分法；初步学会综合应用一种程序设计语言进行科学与工程计算；增强科技报告写作能力，学会相关软件使用。 六、实验内容： 电源/地板系统中的同步开关噪声仿真 七、实验器材（设备、元器件）： 电子计算机 八、实验步骤： 1．Mur二阶吸收边界条件 Engquist和Majda导出的二阶近似解析吸收边界条件为： 对于矢量麦克斯韦方程的FDTD模拟，吸收边界条件式～中的代表位于网格边界的或的各个切向分量。 由于教材上是以边界为例来推导Mur二阶吸收边界条件的，因此此处推导边界的Mur二阶吸收边界条件。如图1所示，令代表位于网格边界的或的某一切向分量，然后在距离网格边界半个 步长的辅助网格点，对分量，将式中的偏微分用中点差分来代替。 图1.MUR吸收边界条件示意图 首先，将式中关于的混合偏导用中点差分写成 然后，式中对的二阶偏导数可以写成相邻两点和处时间偏导的平均值。 当然，对的二阶偏导数可以写成相邻两点和处对偏导的平均值。 将有限差分表达式～代入式，可以解出，就得到在网格边界处的二阶吸收边界条件。 对于另外三个边界处的吸收条件，就不详细写推导过程了，直接给出结果： 边界： 边界： 边界： 以上就是根据Engquist和Majda的二阶近似解析吸收边界推导出来的Mur吸收边界条件。 2．电源/地板系统中的同步开关噪声仿真。 电源地板系统示意图如下： 图2.电源地板系统示意图 电源/地板系统简化模型(平行板模型，电磁场不随z变化) 图中的 设位于处，且 其中， 由于场不随变化，因此有；再加上考虑到导体表面切向电场 为0，即因此平行板之间的场分量只有三个分量，网格的划分如下图所示： 图3.YEE网格的划分 相应的有源FDTD方程差分格式为： 对于互阻抗的提取，应根据互阻抗的定义式： 来计算，其中表示观测点电压的频谱值，由电压观测值进行傅氏变换得到；表示源点的电流的频谱值，由源点电流时域值进行傅氏变换得到，源点电流时域值可由安培环路定理计算得到： 如果要计算点到 点的互阻抗，那么需要在点加一定激励源，然后利用安培环路定理计算点处的电流观测值，以及点处的电压观测值，然后进行后续计算得到互阻抗。 图4.安培环路定理计算电流示意图 3．测试题 我们这一组的测试题是第12题，电源地板的结构图如下，其中尺寸参数是自行设置的（中间的小黑点是设置的激励源位置）。 图5.测试题12电源地板系统 由于开路边界是当作磁壁处理的，所以有必要推导一下磁壁的处理方法。处理磁壁的边界，可以采用镜像法来处理，对于一边是磁壁的情况，示意图如下： 图6.磁壁的网格划分示意图 可以在磁壁以外假设一个磁场分量，它和它关于磁壁对称的那个磁场分量是大小相等符号相反的，这样，电场分量周围就有了四个磁场分量，就可以计算了。 对于处在两条正交磁场边界交点上的电场值，可以假设出两个镜像的磁场分量的点，然后进行计算，由于此处测试题里不会遇到这种情况，就不附示意图了。 为了计算方便以及结果的对称性，选择把激励源加在坐中心处，坐标为(30mm,20mm)，同时把中间所有的开路边界都取在整格上，这样就没有磁场切向分量会因为处在磁壁上而为0，就可以统一进行计算。 测试题提取互阻抗的过程与前边的例题是一样的，就不累述了，只在后边给出结果。 九、实验数据及结果分析： 1. 电源/地板系统中的同步开关噪声仿真（四周匹配） 在实验中，取了三个观测点，坐标分别为：(40mm,25mm), (30mm,20mm),(15mm,15mm),激励源电流加在点(10mm,10mm)处，下边是仿真结果。 图7.电源地板三个观测点电压时域波形 这是三个观测点电压以及源点电流的时域波形图，由三个观测点的电压及源点的电流时域波形图可以看出，即使四周加的是吸收边界，但是由计算导致的虚假反射还是很明显的，而且离观测点越远（或者说离边界越近）