混合基函数的天线模型载荷位移测量方法.pdf

《混合基函数的天线模型载荷位移测量方法.pdf》由会员分享，可在线阅读，更多相关《混合基函数的天线模型载荷位移测量方法.pdf（13页完成版）》请在专利查询网上搜索。

1、(19)中华人民共和国国家知识产权局 (12)发明专利申请 (10)申请公布号 (43)申请公布日 (21)申请号 202010014864.4 (22)申请日 2020.01.07 (71)申请人 浙江大学 地址 310058 浙江省杭州市西湖区余杭塘 路866号 (72)发明人 刘振宇袁沛杨明龙裘辿 谭建荣 (74)专利代理机构 杭州求是专利事务所有限公 司 33200 代理人 林超 (51)Int.Cl. G06F 30/20(2020.01) G06F 111/10(2020.01) (54)发明名称 一种混合基函数的天线模型载荷位移测量 方法 (57)摘要 本发明公开了一种混合基函数。

2、的天线模型 载荷位移测量方法。 针对梁利用NURBS基函数构 造梁单元， 计算基本单元的贝塞尔提取算子， 耦 合 “梁-梁” 单元与 “梁-壳” 单元； 并建立和整理约 束方程， 然后根据不同基本单元之间的耦合关系 进行判断设位移系数矩阵设置， 采用罚函数方法 将位移约束方程融入到整个等几何模型的系统 刚度矩阵中， 建立考虑带罚函数的修正后的势能 泛函， 求解上述最终系统平衡方程获得天线控制 顶点的位移。 本发明采用不同的样条基函数离散 阵元天线复杂产品的不同结构类型构件， 提高了 产品模型建模的灵活性。 权利要求书2页 说明书5页 附图5页 CN 111241672 A 2020.06.05。

3、 CN 111241672 A 1.一种混合基函数的天线模型载荷位移测量方法， 其特征在于， 针对天线的等几何模 型， 线分为支撑框架和反射面， 支撑框架是由多根梁连接构成， 反射面是由多个壳连接构 成； 方法包括以下步骤： (1)针对梁利用NURBS基函数构造梁单元， 针对壳利用非结构化T样条基函数构造壳单 元， 梁单元和壳单元均作为基本单元； 通过计算基本单元的贝塞尔提取算子， 将NURBS基函 数和非结构化T样条基函数转换为定义在-1,1空间的伯恩斯坦多项式的线性组合， 这样 避免了高斯求积时的参数映射； (2)耦合 “梁-梁” 单元与 “梁-壳” 单元，“梁-梁” 单元是指两个梁单元相。

4、连接构成的单 元，“梁-壳” 单元是指梁单元和壳单元相连接构成的单元； 对于 “梁-梁” 单元耦合建立以下位移约束方程： u1,1un,2 其中， u1,1和un,2分别表示其中一个梁单元的控制顶点P11和另一个梁单元的控制顶点Pn2 的位移向量； 对于 “梁-壳” 单元耦合， 分相切与相交分别建立以下位移约束方程： uIuJ 其中， uI和uJ分别表示梁单元与壳单元上控制顶点的位移矢量， 下标I和J分别表示控制 顶点的编号， p表示NURBS基函数阶次， ue,j表示壳单元e对应的控制顶点的位移矢量， j表示 控制顶点在壳单元e中的局部编号，表示壳单元与梁单元的交点对应的基函数， 是 壳单元。

5、与梁单元的交点的参数坐标； (3)整理约束方程信息 将位移约束方程整理成如下形式， 然后根据不同基本单元之间的耦合关系进行判断设 位移系数矩阵设置C1和C2： (C1-C2)u0 其中， u表示整个等几何模型所有N个控制顶点的位移矢量， C1、 C2分别表示基本单元的 控制顶点的位移系数矩阵； (4)采用罚函数方法将位移约束方程融入到整个等几何模型的系统刚度矩阵K中； 建立 考虑带罚函数的修正后的势能泛函*： 其中， S为罚函数方法的惩罚参数， u表示控制顶点的位移， F表示施加到天线的等几何 模型上的载荷； 对势能泛函求取极小值， 令： 整理得以下最终系统平衡方程： K+S(C1-C2)T(。

6、C1-C2)uF 权利要求书 1/2 页 2 CN 111241672 A 2 K*uF K*K+S(C1-C2)T(C1-C2) 其中， K*表示融入位移约束方程后的等几何模型的系统刚度矩阵K， T表示矩阵转置； 求解上述最终系统平衡方程获得天线控制顶点的位移u(K*)-1F。 2.根据权利要求1所述的一种混合基函数的天线模型载荷位移测量方法， 其特征在于： 所述的步骤(3)中， 根据不同基本单元之间的耦合关系进行判断设位移系数矩阵设置， 具体 为： 对于 “梁-梁” 耦合情况， C1、 C2分别表示两根梁单元的端点对应的控制顶点的位移系数 矩阵： 其中， C1、 C2均为33N阶矩阵， E。

7、为三阶单位矩阵， I、 J分别表示两根梁单元控制顶点的 编号； 对于 “梁-壳” 相切耦合情况， 即梁单元和壳单元相切， C1， C2分别表示梁单元和壳单元的 端点对应的控制顶点的位移系数矩阵： 其中， C1、 C2均为33N阶矩阵， E为三阶单位矩阵， I、 J分别表示梁单元和壳单元控制顶 点的编号； 对于 “梁-壳” 相交耦合情况， 即梁单元和壳单元相交， C1， C2分别表示梁单元和壳单元相 交的交点所在梁单元与壳单元处端点对应的控制顶点的位移系数矩阵： 其中， J1,J(p+1)2表示梁单元和壳单元相交的交点所在壳单元的各个控制顶点的编 号；表示壳单元与梁单元的交点对应的基函数。 3.。

8、根据权利要求1所述的一种混合基函数的天线模型载荷位移测量方法， 其特征在于： 所述步骤(4)中， 求解上述最终系统平衡方程获得控制顶点的位移u(K*)-1F。 4.根据权利要求1所述的一种混合基函数的天线模型载荷位移测量方法， 其特征在于： 所述的天线具体为有源相控阵天线。 权利要求书 2/2 页 3 CN 111241672 A 3 一种混合基函数的天线模型载荷位移测量方法 技术领域 0001 本发明涉及了一种复杂产品模型载荷位移测量方法， 尤其是涉及了一种混合基函 数的天线模型载荷位移测量方法。 背景技术 0002 天线属于多种结构类型组成的复杂工业产品， 尺寸庞大， 包含壳、 梁、 杆等。

9、结构。 不 同结构的几何复杂度不同， 其等几何模型的构建难度也不同。 如梁杆结构往往可以通过简 化， 由其中性层曲线来表示整体的物理性质， 因此， NURBS就足以构建梁杆的等几何模型。 而 壳结构通常无法用单个NURBS或Hermite曲面片描述， T样条更适合做壳结构的等几何基函 数。 0003 对于由不同结构类型构件组成的天线产品的等几何模型， 各类结构件基本单元间 的耦合拼接(如梁结构的线单元、 壳结构的面单元等)需要额外的处理， 否则各类结构件的 网格单元无法构成统一的整体， 无法进行天线模型整体位移的精确分析。 替代方法则将天 线产品的不同结构单独进行分析， 将支架结构与反射面结构。

10、的连接点作为发射面结构的固 定支撑点， 限制了固定支撑点的自由度， 忽略了支架结构在实际情况中的变形， 使得天线发 射面的位移计算结果与真实值出现偏差， 不够精确。 发明内容 0004 为了解决背景技术中存在的问题， 本发明提出了一种混合基函数的天线等几何模 型载荷位移测量方法。 本发明采用不同的样条基函数离散阵元天线复杂产品的不同结构类 型构件， 提高了产品模型建模的灵活性和测量准确性。 0005 本发明所采用的技术方案是： 0006 本发明针对天线的等几何模型， 等几何模型具有系统刚度矩阵K， 天线分为支撑框 架和反射面， 支撑框架是由多根梁连接构成， 反射面是由多个壳连接构成； 方法包括。

11、以下步 骤： 0007 (1)针对梁利用NURBS基函数构造梁单元， 针对壳利用非结构化T样条基函数构造 壳单元， 梁单元和壳单元均作为基本单元； 通过计算基本单元的贝塞尔(Bzier)提取算子， 将NURBS基函数和非结构化T样条基函数转换为定义在-1,1空间的伯恩斯坦(Bernstein) 多项式的线性组合， 这样避免了高斯求积时的参数映射； 0008 (2)耦合 “梁-梁” 单元与 “梁-壳” 单元，“梁-梁” 单元是指两个梁单元相连接构成的 单元，“梁-壳” 单元是指梁单元和壳单元相连接构成的单元； 0009 对于 “梁-梁” 单元耦合建立以下位移约束方程： 0010 u1,1un,2。

12、 0011 其中， u1,1和un,2分别表示其中一个梁单元的控制顶点P11和另一个梁单元的控制顶 点Pn2的位移向量； 0012 对于 “梁-壳” 单元耦合， 分相切与相交分别建立以下位移约束方程： 说明书 1/5 页 4 CN 111241672 A 4 0013 uIuJ(相切) 0014(相交) 0015 其中， uI和uJ分别表示梁单元与壳单元上控制顶点的位移矢量， 下标I和J分别表示 控制顶点的编号， p表示NURBS基函数阶次， ue,j表示壳单元e对应的控制顶点的位移矢量， j 表示控制顶点在壳单元e中的局部编号，表示壳单元与梁单元的交点对应的基函数， 是壳单元与梁单元的交点的。

13、参数坐标； 0016 (3)整理约束方程信息 0017 将位移约束方程整理成如下形式， 然后根据不同基本单元之间的耦合关系进行判 断设位移系数矩阵设置C1和C2： 0018 (C1-C2)u0 0019 其中， u表示整个等几何模型所有N个控制顶点的位移矢量， C1、 C2分别表示基本单 元的控制顶点的位移系数矩阵； 0020 (4)采用罚函数方法将位移约束方程融入到整个等几何模型的系统刚度矩阵K中， 即融合约束方程信息到系统刚度矩阵； 建立考虑带罚函数的修正后的势能泛函*： 0021 0022 其中， S为罚函数方法的惩罚参数， u表示控制顶点的位移， F表示施加到天线的等 几何模型上的载荷。

14、； 0023 对势能泛函求取极小值， 令： 0024 0025 整理得以下最终系统平衡方程： 0026 K+S(C1-C2)T(C1-C2)uF 0027 K*uF 0028 K*K+S(C1-C2)T(C1-C2) 0029 其中， K*表示融入位移约束方程后的等几何模型的系统刚度矩阵K， T表示矩阵转 置； 0030 梁与梁、 梁与壳的耦合仅改变了系统的刚度矩阵， 上述变化后的刚度矩阵K*仍然 是对称的， 其维数也与原刚度矩阵相同， 给最终系统平衡模型求解带来极大的好处和便利 性。 0031 求解上述最终系统平衡方程获得天线控制顶点的位移u(K*)-1F。 0032 所述的步骤(3)中， 。

15、根据不同基本单元之间的耦合关系进行判断设位移系数矩阵 设置， 具体为： 0033 对于 “梁-梁” 耦合情况， C1、 C2分别表示两根梁单元的端点对应的控制顶点的位移 系数矩阵： 0034 说明书 2/5 页 5 CN 111241672 A 5 0035 0036 其中， C1、 C2均为33N阶矩阵， E为三阶单位矩阵， I、 J分别表示两根梁单元控制顶 点的编号； 0037 对于 “梁-壳” 相切耦合情况， 即梁单元和壳单元相切， C1， C2分别表示梁单元和壳单 元的端点对应的控制顶点的位移系数矩阵： 0038 0039 0040 其中， C1、 C2均为33N阶矩阵， E为三阶单位。

16、矩阵， I、 J分别表示梁单元和壳单元控 制顶点的编号； 0041 对于 “梁-壳” 相交耦合情况， 即梁单元和壳单元相交， C1， C2分别表示梁单元和壳单 元相交的交点所在梁单元与壳单元处端点对应的控制顶点的位移系数矩阵： 0042 0043 0044 其中， J1,J(p+1)2表示梁单元和壳单元相交的交点所在壳单元的各个控制顶点 的编号；表示壳单元与梁单元的交点对应的基函数。 0045 所述步骤(4)中， 求解上述最终系统平衡方程获得控制顶点的位移u(K*)-1F。 0046 所述的天线具体为有源相控阵天线。 0047 本发明根据步骤(2)、 (3)将基于不同基函数的等几何梁壳单元融合。

17、到同一整体。 0048 本发明根据步骤(4)对载荷作用区域进行适当的数值积分处理， 实现等几何分析 的任意区域载荷施加。 0049 本发明具有的有益效果是： 0050 本发明将不同基函数的基本单元耦合为一个整体， 实现了天线模型不同结构的整 体分析， 实现了天线产品整体位移的精确测量。 0051 本发明提出了将由不同基函数表述的等几何模型精确构建成整体的一种混合基 函数的天线等几何模型载荷位移测量方法， 构建了基于NURBS基函数的三维空间欧拉-伯努 利(Euler-Bernoulli)梁单元之间的强耦合， 考虑了梁单元与基于非结构化T样条的克希霍 夫-拉夫(Kirchhoff-Love)壳单。

18、元之间的两种耦合情况， 将不同基函数的基本单元的耦合 信息进行融合， 实现混合基函数等几何模型的刚度矩阵构造， 实现了天线模型不同结构的 整体分析， 最终实现天线产品整体位移的精确测量， 提高了精确度。 附图说明 0052 图1是本发明的方法流程图。 说明书 3/5 页 6 CN 111241672 A 6 0053 图2是实施例的天线的结构原理图。 0054 图3是45度工作角时的天线及其外载示意图。 0055 图4是天线支架梁截面的3种类型示意图。 0056 图5是混合基函数等几何模型和传统有限元分析模型对比示意图。 0057 图5(a)是天线的混合基函数等几何模型图。 0058 图5(b。

19、)是天线的有限元分析模型图。 0059 图6是混合基函数等几何分析得到的位移云图。 0060 图7是用于电磁分析的天线反射面等几何模型取点示意图。 0061 图8是phi0度极化平面上的2D增益方向图。 具体实施方式 0062 下面结合天线模型对本发明作进一步说明。 0063 本发明的实施例如下： 0064 以图2所示的天线为例。 天线结构在自重、 风载等外力作用下产生的精度误差会改 变天线孔径的电磁振幅与相位的分布， 图中红色虚线表示变形后的反射面， 形状误差将 最终影响天线的远场电性能。 0065 (1)分别用基于NURBS的Euler-Bernoulli梁单元和基于非结构化T样条的 Ki。

20、rchhoff-Love壳单元构建天线结构的桁架和反射面， 如图3， 黄色部分为主反射面， 蓝色 部分为天线支架。 分析在45度工作角时， 受重力、 20m/s风力载荷作用的反射面天线形面误 差。 0066 其中 “梁-梁” 单元耦合方法用于天线支架梁结构间的耦合连接，“梁-壳” 单元耦合 方法用于天线支架梁结构与天线反射面间的耦合连接。 0067 (2)天线支架梁结构共包含3种梁截面形式， 如图4所示为各梁截面形式。 0068 (3)分别采用混合基函数等几何分析与传统有限元分析方法计算反射面天线在重 力和风力载荷下的最大位移和反射面形面误差， 如图5分别为混合基函数等几何模型和传 统有限元分。

21、析模型。 0069 (4)取不同网格粗细程度， 对比混合基函数等几何分析与传统有限元分析的计算 结果随节点数目变化的收敛情况， 如表1， 其中N表示等几何分析的控制顶点数和传统有限 元分析的单元节点数。 0070 表1最大位移与形面误差的收敛分析 0071 0072 为了达到相同的计算精度， 混合基函数等几何分析仅用了不到五分之一的计算规 模。 如图6为混合基函数等几何分析得到的位移云图。 0073 (5)在等几何模型上取足够细密的点， 使其最小距离满足电磁性能分析所需的小 于0.3倍波长的要求， 如图7， 以这些点生成包含位移信息的三角网格模型， 分析反射面形变 说明书 4/5 页 7 CN。

22、 111241672 A 7 误差对增益的影响。 0074 (6)计算得到变形反射面在phi0度极化平面上的2D增益方向图， 如图8， 变形后， 天线增益也发生变化， 且由混合基函数等几何分析模型计算结果得到增益图与最细网格的 传统有限元方法得到的增益图相近， 极大地避免了因几何离散误差导致的误差放大问题。 说明书 5/5 页 8 CN 111241672 A 8 图1 说明书附图 1/5 页 9 CN 111241672 A 9 图2 图3 说明书附图 2/5 页 10 CN 111241672 A 10 图4 图5 说明书附图 3/5 页 11 CN 111241672 A 11 图6 图7 说明书附图 4/5 页 12 CN 111241672 A 12 图8 说明书附图 5/5 页 13 CN 111241672 A 13 。