电磁波导的半解析辛分析lowast;.pdfVIP

电磁波导的半解析辛分析∗ 钟万勰 (大连理工大学 工程力学研究所, 辽宁大连 116023) 提 要 根据电磁波导的 Hamilton 体系 辛几何可用于任意各向异性材料 而且便于处理不同区段的界面 条件 横向的电场和磁场构成了对偶向量 基于 Hamilton 变分原理用半解析横向离散应当保持体系的辛结 构 对于半解析有限元之间插值函数的变分犯规 引入了变分原理的补偿项 离散后可以运用应用力学的 有效算法 求解了其辛本征值问题 每段波导可以引入两端 Riccati 矩阵 精细积分法可以用于求解其方程 组 关键词: 波导; 哈密顿体系; 辛 半解析法 1 频域 Maxwell 方程的变分原理 各向同性材料的波导 其分析总是由横向电波(TE)或横向磁波(TM)组成 然而一般各向 异性材料的波导并不能简单地用 TE 波及 TM 波来区分 见[1~4]等 通过变分原理寻求一般 条件下的解法 将体系导入 Hamilton 对偶方程的辛几何形式 是有效的途径 但除个别情 况外要寻找波导的解析解很困难 因此采用在横截面离散的半解析解数值方法是自然的 离 散后的辛本征解或精细积分可运用结构力学与控制中的方法 本文将做介绍 坐标 其横向是 的直角坐标 介电(permittivity)矩阵 是 将波导的长度方向取为 z x ,y ε 3 ×3 对称正定阵 可以是常矩阵也可以是横向坐标x ,y 的函数 给出了横向的函数 还可 以有长度方向的变化 例如是周期性的变化 在光纤调制中的光栅就是折射率对长度坐标z 周期变化的 即周期变化介电常数的光纤材料 Maxwell 方程是在时域列式的 频域列式可以表示为 ω − ω −i t i t H he E iee (1.1) 其中e(x ,y ,z ,ω) h(x ,y ,z ,ω) 待求 频域 Maxwell 方程为 ωµ h R ⋅e , ωεe R ⋅h (1.2a,b) 0 其中光学材料的介磁(permeability)矩阵可看成µ µ I µ 是真空中的介电常数 而 0 3 0  0 −∂ ∂z ∂ ∂y  R  ∂ ∂z 0 −∂ ∂x  (1.3) −∂ ∂y ∂ ∂x 0  是算子矩阵 设有限域V (谐振腔)其边界为S 理想导体的边界条件是 n ×e 0 , n i l +i m +i n 在边界 S (1.4) x y z 其中i ,i ,i 是沿坐标轴的单位向量 如果是无限区域 则要求向无穷远衰减足够快