本文属于类型a：单一原创研究（original research）的学术论文。以下是根据该文档内容撰写的学术报告。
作者与发表信息
本文由Olli Luukkonen、Constantin Simovski、Gérard Granet、George Goussetis、Dmitri Lioubtchenko、Antti V. Räisänen和Sergei A. Tretyakov共同完成。主要研究机构包括芬兰赫尔辛基理工大学（TKK Helsinki University of Technology）无线电科学与工程系、俄罗斯圣彼得堡精密机械与光学研究所、法国Université Blaise Pascal的LASMEA实验室以及英国赫瑞瓦特大学。该研究发表于IEEE Transactions on Antennas and Propagation期刊，2008年6月，第56卷第6期。
研究背景与目的
本研究属于电磁场与天线传播领域，具体聚焦于平面网格（planar grids）和高阻抗表面（high-impedance surfaces, HIS）的解析建模。高阻抗表面通常由金属条带（metal strips）或方形贴片（square patches）阵列置于接地介质基板上构成，在微波和毫米波器件中具有重要应用，如人工磁导体（artificial magnetic conductors）、反射阵列天线（reflect-array antennas）、可调移相器（tunable phase shifters）等。此前的研究中，对于金属条带阵列在垂直入射条件下的等效电路模型已有报道，但对于方形贴片阵列以及斜入射（oblique incidence）条件下的精确且简单解析模型尚缺乏。尤其是当相邻金属单元之间的间隙非常窄、电容耦合强时，现有模型如Holloway等人的模型基于偶极子相互作用假设，精度显著下降。本文的目标是从已知的条带网格解析模型出发，结合近似巴比涅原理（approximate Babinet principle），推导出电密排布的方形贴片阵列的网格阻抗（grid impedance）解析表达式，并进一步获得基于此类阵列的高阻抗表面的表面阻抗（surface impedance）和反射系数的解析公式，实现对大角度斜入射的高精度预测。
研究流程与方法
本研究主要包括三个递进的解析建模与数值验证阶段。
第一阶段的工作对象是均匀介质中的方形贴片阵列。作者首先从Kontorovich提出的平均边界条件（averaged boundary conditions, ABC）出发，给出了由沿x轴和y轴的金属条带构成的网格（mesh of strips）在TE和TM极化下的网格阻抗表达式。对于位于介质界面处的网格，采用有效介电常数近似，即取基板相对介电常数与自由空间介电常数的平均值，从而将均匀介质模型推广到介质界面情况。然后，利用近似巴比涅原理，将网格阻抗经公式 (Z_g^{TM} Zg^{TE} = \eta{eff}^2 /4) 变换，得到互补结构——贴片阵列的网格阻抗。对于TM极化，贴片阵列的网格阻抗与平行电容条带网格相同；对于TE极化，二者相差一个与波数有关的因子，该因子取决于沿网格平面内波矢分量与自由空间波矢的关系。由此获得了贴片阵列在自由空间中的TE和TM网格阻抗表达式。
第二阶段研究电容条带和贴片阵列的反射与透射系数。利用传输线模型将自由空间波阻抗与网格阻抗并联，计算输入阻抗进而得到反射系数。其中自由空间TE波阻抗为 (\eta_0/\cos\theta)，TM波阻抗为 (\eta_0\cos\theta)。为验证模型精度，作者采用基于Floquet展开的傅里叶模态法（Fourier modal method）进行全波数值模拟。该方法在金属边缘处提高空间分辨率，能有效处理场奇异性，且算法稳定性良好。作者将本文模型与Holloway等人的模型进行了对比。Holloway模型给出的TE极化网格阻抗具有不同的角度依赖形式，其网格参数采用不同近似，但本文模型与之在角度依赖性上有明显差异。
第三阶段将模型推广至高阻抗表面。高阻抗表面由电容网格或贴片阵列置于厚度为h、相对介电常数为(\varepsilon_r)的接地介质基板上构成。其输入阻抗是网格阻抗与接地介质层表面阻抗的并联。接地介质层的表面阻抗对TE和TM极化分别通过公式 (Z_s^{TE} = j\omega\mu_0 h)（近似）和 (Z_s^{TM} = j\omega\mu_0 h \cos^2\theta) 或其他相关表达式给出。对于贴片阵列，TM极化输入阻抗与电容条带相同，TE极化则因贴片阵列网格阻抗中的额外因子而不同。最终得到高阻抗表面在斜入射下的反射相位随频率变化的理论预测。数值验证采用了Ansoft HFSS 10.1.1全波仿真和傅里叶模态法，模拟对象为无限周期结构。验证用的高阻抗表面参数为：周期d=2 mm，金属宽度w=0.2 mm，基板厚度h=1 mm，相对介电常数(\varepsilon_r=10.2)。入射角从0°到近掠入射角均进行了对比。
主要结果
在第一阶段中，通过解析推导得到了贴片阵列的网格阻抗表达式。关键结果是，对于TE极化，贴片阵列的网格阻抗与平行电容条带阵列的网格阻抗不同，前者分母中多了一个与波数相关的因子。该差异来源于贴片阵列在二维平面内的周期性，而平行条带阵列只在一个方向上具有周期性。这一角度依赖性是本文模型与Holloway模型的主要区别之一。
在反射系数的数值验证中，图4给出了TE和TM极化反射系数随入射角变化的曲线。对于周期为波长的十分之一、金属宽度为周期的十分之一的贴片阵列，本文模型在整个入射角范围内与傅里叶模态法结果吻合很好，而Holloway模型（图中以(r_H)标记）在TE极化下偏差明显较大。这表明本文模型在狭窄间隙、强电容耦合条件下具有更高的精度。
在高阻抗表面反射相位验证中，图6和图7分别给出了电容条带高阻抗表面在TE和TM极化下的反射相位图。图8和图9则对应方形贴片高阻抗表面。所有图中，解析结果与HFSS及傅里叶模态法结果均高度一致，特别是在主谐振频率附近。对于TM极化，谐振带宽随入射角增大而增大；对于TE极化，谐振带宽随入射角增大而减小。作者指出，这一现象并非由表面阻抗本身随角度变化直接引起，而是由于自由空间波阻抗随角度变化，改变了其与高阻抗表面表面阻抗的比值，进而影响反射系数的频率响应。
第四阶段研究了沿高阻抗表面传播的表面波。利用横向谐振条件 (Z_s + Z_0 = 0) 求解TE和TM表面波的色散关系。在无损耗情况下，表面阻抗为纯虚数，自由空间阻抗也需为纯虚数，因此入射角需取虚数，对应于沿表面传播、垂直表面方向衰减的表面波。图10给出了TE和TM表面波的色散曲线，解析结果与HFSS仿真结果吻合良好。这进一步证明了本文模型不仅适用于平面波激励，也适用于表面波传播特性的预测。
结论与价值
本文推导了电密排布的方形贴片阵列的网格阻抗解析表达式，并基于该结果建立了高阻抗表面在斜入射下的解析模型。通过与傅里叶模态法和HFSS仿真结果的系统对比，证明该模型在主谐振频率范围内对所有入射角（直至近掠入射）和两种极化均具有很高的精度。与文献[27]中的模型相比，本文模型在贴片间距狭窄、电容耦合强的条件下精度更高。此外，模型还能准确预测高阻抗表面的表面波传播特性。
在应用价值方面，该模型可用于人工阻抗表面、人工磁导体、反射阵列天线、可调移相器以及吸收体等器件的设计。由于解析表达式简单、计算量小，适合嵌入到复杂系统的优化设计流程中。此外，该模型还可用于推导更高阶的阻抗边界条件（higher-order impedance boundary conditions），为人工高阻抗表面的进一步理论研究提供基础。
研究亮点
本研究的核心贡献在于：第一，通过严格应用近似巴比涅原理，首次给出了方形贴片阵列在斜入射下的网格阻抗解析表达式，并揭示其与平行电容条带阵列在TE极化下的角度依赖性差异；第二，系统验证了有效介电常数近似在介质界面处电容网格建模中的适用性；第三，证明了简单的传输线模型结合巴比涅原理能够以高精度预测高阻抗表面的平面波反射和表面波传播特性；第四，与已有解析模型相比，在强电容耦合条件下具有显著更高的精度，克服了偶极子近似模型的局限性。