基于近似与灵敏度分析的反射面天线机电集成优化设计方法
【技术领域】
[0001] 本发明属于雷达技术领域，具体涉及雷达天线领域中的一种基于近似与灵敏度分 析的反射面天线机电集成优化设计方法。
【背景技术】
[0002] 反射面天线广泛应用在通信、雷达、射电天文学、微波通信、卫星通信和跟踪以及 遥感等各个领域。反射面天线结构是典型的机电一体化结构，其机械结构性能与电性能相 互影响、相互制约。为了设计出高性能的反射面天线，需要从学科交叉、机电集成的角度出 发，对反射面天线进行机电集成设计。
[0003] 段宝岩等人在中国专利"基于拟合变形反射面的天线电性能预测方法"中，公开了 一种基于拟合变形反射面的天线电性能预测方法。该方法利用实际变形面构造一相近的 拟合变形反射面，并在拟合变形反射面上进行电性能计算。但该方法存在的不足是，采用拟 合变形反射面代替实际反射面进行电性能计算，其计算量比较大；同时，由于该方法无法提 供较为准确的优化迭代搜索方向，使得优化迭代费时，优化时间较长。郑飞等人在中国专利 "基于误差因素的反射面天线机电综合分析方法"中，公开了一种基于误差因素的反射面天 线机电综合分析方法。该方法通过将结构网格转化为电磁分析网格，进行电性能计算。但 该方法同样存在计算时间长，无法提供较为准确的优化迭代搜索方向，优化迭代耗时的问 题。因此有必要针对上述方法计算时间长、迭代搜索方向不准确的问题，进行基于近似与灵 敏度分析的反射面天线机电集成优化设计方法研究。
【发明内容】
[0004] 本发明的目的是克服上述现有技术的不足，提供一种基于近似与灵敏度分析的反 射面天线机电集成优化设计方法。通过采用近似计算方法获得变形反射面天线远区电场， 通过灵敏度信息构建优化迭代搜索方向，进而实现反射面天线机电集成优化设计。
[0005] 本发明的技术方案是：基于近似与灵敏度分析的反射面天线机电集成优化设计方 法，包括如下步骤：
[0006] (1)输入反射面天线结构参数和电参数
[0007] 输入用户提供的反射面天线的结构参数和电参数信息，其中结构参数包含口径、 焦距、反射面板尺寸与厚度参数、背架辐射梁尺寸与截面参数、中心体尺寸参数和载荷参 数，电参数包含工作波长、馈源总辐射功率参数和要求的电性能参数；
[0008] (2)建立机电集成优化模型
[0009] 从反射面天线结构参数信息中提取反射面背架辐射梁尺寸与截面参数，将背架辐 射梁尺寸与截面参数作为优化模型的设计变量，以反射面天线要求的电性能参数为目标函 数，依此建立反射面天线机电集成优化设计模型：
[0010] Find X = U1, X2,…，xM}T
[0011] Min D(X)
[0012] s.t. G(x)<0
[0013] 其中，Find表示迭代运算，x表示结构设计参数列向量，x2、. . .、xM依次表示编 号为1、2、. ..、M的结构设计参数，M表示结构设计参数总数，上标T表示向量转置运算；Min 表示最小化运算，D(x)表示反射面天线要求的电性能参数，s. t.表示约束运算，G(x)表示 根据设计要求添加的约束函数；
[0014] ⑶建立结构有限元模型
[0015] 根据用户提供的结构参数，计算节点坐标，并根据天线反射面板、背架辐射梁、中 心体的结构参数选择梁单元、壳单元，获得梁单元和壳单元的尺寸、截面与厚度，利用有限 元软件建立结构有限元模型；
[0016] (4)求解结构有限元模型
[0017] 针对已建立的结构有限元模型，添加结构位移、自由度约束或者边界条件；根据反 射面天线的载荷参数，在结构有限元模型上施加工作载荷；在此基础上，利用有限元软件生 成结构刚度矩阵，对结构有限元模型进行求解，获得节点位移、单元应力；
[0018] (5)采用近似方法计算远区电场
[0019] 在步骤（3)和步骤⑷的基础上，采用近似方法计算远区电场；
[0020] (6)判断电性能指标是否满足要求
[0021] 判断远区电场是否满足用户在步骤（1)中指定的电性能要求，如果满足要求，则 转至步骤（11)，否则转至步骤（7);
[0022] (7)结构灵敏度分析
[0023] 在步骤（4)的基础上，计算节点位移对结构设计参数的灵敏度，通过下式获得：
[0024]
.觸:
[0025] 其中，y表示节点位移对结构设计参数的灵敏度，U表示步骤⑷生成的节点位 OX： 移列向量，X表示结构设计参数列向量，.9表示求偏导数运算，K表示步骤（4)生成的结构刚 度矩阵，上标-1表示矩阵求逆运算，P表示步骤（4)中的工作载荷列向量；
[0026] (8)电性能灵敏度分析
[0027] 在步骤（4)的基础上，采用直接微分法计算电性能对反射面节点位移的灵敏度；
[0028] (9)电性能对结构设计参数的灵敏度
[0029] 在步骤（7)和步骤（8)的基础上，通过下式获得电性能参数对结构设计参数的灵 敏度：
[0030]
[0031] 其中，^表示电性能对结构设计参数的灵敏度，D表示电性能参数，X表示结构设 .? 计参数列向量，@表示求偏导数运算：
表示由电性能参数对节点轴向位移的灵敏度信 息组成的矩阵，Z表示从反射面天线结构节点位移信息中提取的节点轴向位移量，表示 节点位移对结构设计参数的灵敏度，U表示节点位移列向量；
[0032] (10)更新结构设计参数
[0033] 采用基于灵敏度信息优化迭代方法，得到下次迭代的结构设计参数，更新结构设 计参数，转至步骤（3);
[0034] (11)输出天线结构设计方案。
[0035] 上述步骤（5)所述的近似方法是一种基于近似技术的反射面天线电性能计算方 法，计算过程如下：
[0036] 5a)计算理想反射面天线的远区电场，计算公式如下：
[0037]
[0038]
[0039] 其中，1&表示理想反射面天线远区电场，肩表示远场观察点位置矢量，j表示虚数 单位，k表示自由空间波数，η表示自由空间波阻抗，exp表示自然对数的指数运算，R表示 远场观察点位置矢量幅度，π表示圆周率，I r表示单位并矢，表示单位矢量左的并矢， Σ表示反射曲面，為F)表示反射面上位置矢量^;处的面电流密度，^表示反射面位置矢量， 表示远场观察点的单位矢量，σ表示投影口面，I.表示单位法向矢量，/7(/=)表示反射面 位置矢量？处的入射磁场。
[0040] 5b)以建立的结构有限元模型为基础，提取有限元模型中处于电磁波照射下的反 射面部分的节点、单元和形函数信息；
[0041] 5c)通过下式计算单元一次系数矩阵：
[0042]
[0043]
[0044] 其中，Af表示单元e的一次系数矩阵，上标e表示结构有限元模型中某一单元，下 标i表示位于单元e上的节点编号，表示单元e的一次系数矩阵的第i个分量，符号e 表示从属关系，NUM表示单元e上的节点总数，表示单元e的法向矢量，P(F)表示反射面 位置矢量P处的入射磁场，？:::表示反射面位置矢量，exp表示自然对数的指数运算，j表示虚 数单位，k表示自由空间波数，身表示远场观察点的单位矢量，Qi表示相对于第i个节点的 形函数，0s表示位置矢量Γ在馈源坐标系下的俯仰角，下标s表示馈源坐标系，Θ表示远 场观察点俯仰角，σ 3表示单元e在口径面内的投影面积；
[0045] 5d)通过下式计算单元二次系数矩阵：
[0046]
[0047]
[0048] 其中，表示单元e的二次系数矩阵，上标e表示结构有限元模型中某一单元，u 和V分别表示位于单元e上的节点编号，表示由节点u和V构成的单元e的二次系数 矩阵分量，符号e表示从属关系，NUM表示单元e上的节点总数，k表示自由空间波数，#表 示单元e的法向矢量，#(?|表示反射面位置矢量^处的入射磁场，表示反射面位置矢量， exp表示自然对数的指数运算，j表示虚数单位，及表示远场观察点的单位矢量，Qu表示相 对于第u个节点的形函数，Qv表示相对于第V个节点的形函数，Θ s表示位置矢量P在馈源 坐标系下的俯仰角，下标s表示馈源坐标系，Θ表示远场观察点俯仰角，〇；3表示单元e在 口径面内的投影面积；
[0049] 5e)通过下式组集总体一次系数矩阵：
[0050]
[0051] 其中，H1表示总体一次系数矩阵，/<表示单元e的一次系数矩阵，上标e表示结构 有限元模型中某一单元，m表示单元总数，A表示有限元组集运算；
[0052] 5f)通过下式组集总体二次系数矩阵：
[0053]
[0054] 其