采用基于NURBS曲面建模的物理光学方法计算散射近场.pdf

2005’全国微波毫米波会议论文集 采用基于NURBS曲面建模的物理光学方法 计算散射近场 陈铭张玉梁晶洪 西安电子科技大学天线与微波技术国家重点实验室710071 摘要：本文采用基于非均匀有理B样条(NURBS)曲面建模的物理光学方法(PO)计算电 大尺寸理想导体的散射近场。文章推导了在平面波入射的情况下基于NURBS建模的散射近 场计算公式。应用驻相法计算了有理贝齐尔曲面上的物理光学积分。通过与三角形平面片建 模的物理光学方法得到的计算结果对比，说明了公式和方法的有效性。 关键字：NuRBS；物理光学；驻相法；近场 1引言 非均匀有理B样条(NURBs)建模技术采用NURBS曲面描述散射体，具有建模精度高、 面片数量少的优点，在1991年已经成为工业产品几何定义的国际标准。将先进的NURBS建 模技术引入电磁计算不仅可以提高电磁建模的精度而且可以提高计算的速度。因此，采用基 于NURBS建模技术的高频电磁计算已经得到广泛的研究[1。3】。 到目前为止，大部分的研究工作关注的是雷达散射截面的计算【1’31，而雷达散射截面的 计算关注的是远区散射场。据作者所知，尚没有关于NURBS建模的近场计算的相关文献发 表。由于近场问题对于电磁兼容分析十分关键，因此本文主要研究采用基于NURBS曲面建 模的物理光学方法计算电大尺寸目标的散射近场。 2平面波照射下的近场公式 图I NURBS曲面被平面波照射 图l为NURBS曲面在平面波入射下的散射情况。图中尹为场点坐标，J。是曲面表面 的感应电流，R=I R I是曲面上一点到场点的距离。 由电磁理论可以推导出散射场的磁矢位： 确=卷肛竿幽’ ㈣ 由散射场的磁矢位并应用严格的自由空间格林函数可以得到散射近场豆=志j[半e-J从R x(R啦)+半P一脚(_2J,杠2， 2005’全国微波毫米波会议论文集 假设平面波入射的情况下得到入射电场以及入射磁场： 雹j：琶Oe-Jk’。 豆i：!云×豆i：三云×豆。e一∥．， 刁 刁 根据物理光学近似可以得到理想导体表面的感应电流歹。： 了s：2h x旌i：兰Ax琵x雹oe前1：ge前1 77 最终可以得到散射近场的表达式为： 豆：—L．c[G，辰×眩×霞)+G2(_2霞灌一肛c胄托砷幽·j a)s4；w p 1 、 7 ‘、 71 其中a1=坐等巡，G：=半 (3) (4) 3采用驻相法计算物理光学积分 如果对目标采用NURBS曲面建模那就需要在NURBS曲面上计算积分(6)式。NURBS 曲面的数学表示式的基函数是一组递推函数，对其进行求导等运算时比较复杂且不稳定。因 此目前通用的做法是将NURBS曲面分解为有理贝齐尔曲面，有理贝齐尔曲面采用伯恩斯坦 多项式作为基函数，对其进行数学运算比较简单且稳定。问题转化为在有理贝齐尔曲面上计 算积分(6)式。有理贝齐尔曲面的数学表示式为： ∑∑嵋，，巨，JBTO,W；(O 一， 、 i=O，=0 r I“。1，J 2——』—————————————一 ∑EW。筇(”)B；(V) i=O，=O 其中，6f，，为曲的控制点，wf，J为对应于各个控制点的权，m和n为曲面的次数，Em(f)and 彰n(，)是伯恩斯坦多项式。对于参数曲丽，积分面元可以写成ds’=阮×元|枞，其中 艺：姿，瓦：譬。那么，散射场(6)可以简单写成如下形式 E(尹)2衰石， (8)10’sq冗 其中 7=仨￡季(州)eJkf(u,v)dudv (9) 季(材，y)=[G。辰×伍×霞)+G：(_2霞)]J尹·@)×产’(v)I (10) 1278 2005’全国微波毫米波会议论文集 f(u，V)=一(R+k·尹’) (11) 可见，散射场的求解关键就是计算积分式(9)。驻相法可以被应用于计算如(9)式这种类 型的积分。驻相法认为类似(9)式的积分在k1时，积分的主要贡献来自于某些关键点。 关键点可以分为三类： 1曲面内部驻相点，这类关键点满足的条件为 ∞浆0。 m，【工(“o，1，o)= 、’ 2边界上的关键点，这类关键点满足下列四个条件之一： 工(0，yo)=0；fv(1，Vo)=0；无(Uo,0)=0；无(Uo，1)=0； (13) 3曲面角点 这三类驻相点对于积分的贡献可以分别由下面三式计算 7=季(甜o，vo)P彤‘‰％’ 2no尼佩 (14) H_1)“丽jg(uo而,Vo)eJ矿(uo,vo)√蒜焉 ∞，、7巩(甜o，1，o) 1『％(甜o，Vo) 、’ ，=一(一-)Uo+Vo eJ妒(uo,vo)乏i7_石善‰ c-6， 劫的特征值，矩阵元素分别为口=厶、∥=厶俨厶 盯1O，盯20 盯lO，盯20 (17) 口l∥20 式(15)中的口和∥根据边界驻相点处于不同的边界分别代表U或者，。将所有关键点对 积分的贡献相加就可以计算得到