一种板式卫星结构共节点网格快速生成方法技术

一种板式卫星结构共节点网格快速生成方法，它包括以下步骤：步骤1：对导入有限元建模软件的卫星结构几何模型提取中面；步骤2：获取需要处理的中面集合及其元素数目N；步骤3：对所有的面对(A,B)进行有关操作；步骤4：对所有的面进行缝合；步骤5：对结构的几何模型划分有限元网格。通过该方法能对卫星结构三维实体模型抽取中面得到的几何进行自动处理，消除面与面之间的间隙并且删除多余的部分，最终生成的几何模型可以划分高质量的共节点有限元网格。避免了大量的人机交互操作，需要的时间少，通用性强，适用范围广泛，具有很强的工程实用价值。

【技术实现步骤摘要】

：本专利技术涉及一种板式卫星结构共节点网格快速生成方法，属于航天器结构分析的

技术介绍
：目前卫星结构多采用蜂窝夹层板承受和传递力学载荷，并为仪器设备提供安装支撑。为了验证卫星结构的设计是否满足要求，工程中一般采用有限元方法对其进行力学性能分析。卫星结构的几何模型一般采用三维实体进行建模，但在有限元建模时一般采用板壳结构单元进行力学模拟。如果直接由三维实体抽取中面后划分有限元网格，则由于相邻面之间的间隙会导致有限元网格的不连续；如果采用刚性连接单元进行单元连接，则一方面要求相邻面上的网格比较匹配以利于创建刚性连接单元，另一方面刚性连接单元的创建会对结构引入附加的刚度；除此之外，刚性连接单元的创建还会增加有限元分析的计算量，特别是对于卫星结构分析中常见的模态分析和频率响应分析等动力学分析情形。因此工程中希望能生成共节点的有限元网格。卫星结构共节点网格快速生成的主要难点是：由于蜂窝夹层板厚度的存在，从三维实体模型提取中面后得到的几何模型中，相邻面之间存在间隙，无法直接生成有限元网格。为了实现共节点网格的划分，需要首先对面进行延伸或裁剪，从而消除面与面之间的间隙或多余的部分。由于卫星结构中一般有几十块结构板，因此只依靠商业有限元软件自身的几何处理功能并通过人机交互方式实现整个卫星结构共节点有限元网格的生成将需要较高的技巧、大量的时间和烦琐的工作。
技术实现思路
：1、目的：为了实现卫星结构有限元模型的高效生成，本专利技术的目的是提供一种板式卫星结构共节点网格快速生成方法。该方法首先利用有限元软件的中面提取功能得到卫星结构板的中面几何模型，然后通过遍历获取需要进行处理的面对，其次对每组面对依次进行自动处理，从而获得整个卫星结构的几何模型，再次利用有限元软件的几何缝合功能实现将卫星结构的几何模型缝合为一个整体。在此基础上进行网格划分，得到基于板壳单元进行模拟的高质量共节点有限元网格。由于本方法实现自动处理，因此避免了大量的人机交互操作，需要的时间少，通用性强，适用范围广泛。2、技术方案：本专利技术解决其技术问题所采用的技术方案是：本专利技术一种板式卫星结构共节点网格快速生成方法，其特征在于它包括以下步骤：步骤1：对导入有限元建模软件的卫星结构几何模型提取中面；步骤2：获取需要处理的中面集合及其元素数目N：步骤3：对所有的面对(A,B)进行以下操作：步骤3.1：计算A和B的法向量之间的夹角θAB；步骤3.2：若10°＜θAB＜170°，则执行步骤3.3至3.10；步骤3.3：找出A中所有需要向B延伸的边及其数目NA；步骤3.4：若NA＞0，则执行步骤3.5至3.6；否则执行步骤3.7至3.10；步骤3.5：计算A中的边需要向B延伸的长度δA；步骤3.6：将A中所有需要向B延伸的边延伸δA；步骤3.7：找出B中所有需要向A延伸的边及其数目NB；步骤3.8：若NB＞0，则执行步骤3.7；否则执行步骤3.9；步骤3.9：计算B中的边需要向A延伸的长度δB；步骤3.10：将B中所有需要向A延伸的边延伸δB；步骤3.11：利用B分割A；步骤3.12：删除A中面积较小的面；步骤3.13：利用A分割B；步骤3.14：删除B中面积较小的面；步骤4：对所有的面进行缝合；步骤5：对结构的几何模型划分有限元网格；其中，在步骤1中所述的“对导入有限元建模软件的卫星结构几何模型提取中面”的方法为：现有的有限元软件一般都支持薄壁零件的中面提取功能，并且卫星蜂窝夹层板结构一般都是厚度均匀的，因此直接通过有限元软件的用户界面选择卫星模型中蜂窝夹层结构板的几何模型进行中面提取即可得到相应的中面几何模型；其中，在步骤2中所述的“获取需要处理的中面集合及其元素数目N”的方法为：现有的有限元软件一般都提供二次开发接口使得用户可以开发新的功能，因此可以利用相关的函数接口通过二次开发得到模型中所有中面所组成的集合以及其元素数目。获取中面集合所需要的二次开发接口函数可以通过查询有限元软件的相关帮助文档得到；其中，在步骤3.1中所述的“计算A和B的法向量之间的夹角θAB”的方法为：假设A和B的单位法向量分别为nA＝(nAx,nAy,nAz)和nB＝(nBx,nBy,nBz)，则它们之间夹角的余弦值为cosθAB＝nAxnBx+nAynBy+nAznBz，从而利用初等代数中的反三角函数可以得到两个向量之间夹角为θAB＝arccos(nAxnBx+nAynBy+nAznBz)；其中，在步骤3.3中所述的“找出A中所有需要向B延伸的边及其数目NA”的方法为：现有的有限元软件一般都提供二次开发接口使得用户可以开发新的功能，因此可以利用相关的函数接口通过二次开发并按照以下方法实现：1)获取A所有的边EA1,…,EAn和B的一个顶点VB；2)对每条边EAj(j＝1,…,n)，假设其自然坐标表示为Pj(s)(0≤s≤1)，均匀地获取Nj+1个点Pjk＝Pj(k/Nj)(k＝0,1,…,Nj)，其中Nj≥1000；3)取dmax＝-999999,d＝0，然后按照初等解析几何中的公式djk＝nBx(xjk-xB)+nBy(yjk-yB)+nBz(zjk-zB)依次计算每个点Pjk到B的有向距离，若|djk|＞dmax，则dmax＝|djk|,d＝djk；其中(xjk,yjk,zjk)为Pjk在直角坐标系中的坐标，(xB,yB,zB)为VB的坐标，(nBx,nBy,nBz)为B的单位法向量；4)取abd＝999999，对每个点Pjk，若d·djk＞0且|djk|＜abd，则abd＝|djk|；5)若abd＜0.04，则δmax＝abd+0.01；否则δmax＝0.05；6)NA＝0，对每条边EAj：若对k＝0,1,…,Nj均有|djk|≤δmax，则EAj标记为1，NA＝NA+1；否则EAj标记为0；在本步骤中获取A所有的边和B的顶点均可以通过软件的接口实现，并且在有限元软件中边或曲线一般会采用自然坐标进行表示，给定其自然坐标后可以通过接口得到其在直角坐标系下的坐标表示。本步骤中所需要的二次开发接口函数可以通过查询有限元软件的相关帮助文档得到；其中，在步骤3.5中所述的“计算A中的边需要向B延伸的长度δA”的方法为：δA＝δmax/sinθAB，其中δmax的值已经在步骤3.3中得到，θAB的值已经在步骤3.1中得到；其中，在步骤3.6中所述的“将A中所有需要向B延伸的边延伸δA”的方法为：现有的有限元软件一般都提供二次开发接口使得用户可以开发新的功能，因此本文档来自技高网...

【技术保护点】
一种板式卫星结构共节点网格快速生成方法，其特征在于：它包括以下步骤：步骤1：对导入有限元建模软件的卫星结构几何模型提取中面；步骤2：获取需要处理的中面集合及其元素数目N：步骤3：对所有的面对(A,B)进行以下操作：步骤3.1：计算A和B的法向量之间的夹角θAB；步骤3.2：若10°＜θAB＜170°，则执行步骤3.3至3.10；步骤3.3：找出A中所有需要向B延伸的边及其数目NA；步骤3.4：若NA＞0，则执行步骤3.5至3.6；否则执行步骤3.7至3.10；步骤3.5：计算A中的边需要向B延伸的长度δA；步骤3.6：将A中所有需要向B延伸的边延伸δA；步骤3.7：找出B中所有需要向A延伸的边及其数目NB；步骤3.8：若NB＞0，则执行步骤3.7；否则执行步骤3.9；步骤3.9：计算B中的边需要向A延伸的长度δB；步骤3.10：将B中所有需要向A延伸的边延伸δB；步骤3.11：利用B分割A；步骤3.12：删除A中面积较小的面；步骤3.13：利用A分割B；步骤3.14：删除B中面积较小的面；步骤4：对所有的面进行缝合；步骤5：对结构的几何模型划分有限元网格。

【技术特征摘要】
1.一种板式卫星结构共节点网格快速生成方法，其特征在于：它包括以下步骤：
步骤1：对导入有限元建模软件的卫星结构几何模型提取中面；
步骤2：获取需要处理的中面集合及其元素数目N：
步骤3：对所有的面对(A,B)进行以下操作：
步骤3.1：计算A和B的法向量之间的夹角θAB；
步骤3.2：若10°＜θAB＜170°，则执行步骤3.3至3.10；
步骤3.3：找出A中所有需要向B延伸的边及其数目NA；
步骤3.4：若NA＞0，则执行步骤3.5至3.6；否则执行步骤3.7至3.10；
步骤3.5：计算A中的边需要向B延伸的长度δA；
步骤3.6：将A中所有需要向B延伸的边延伸δA；
步骤3.7：找出B中所有需要向A延伸的边及其数目NB；
步骤3.8：若NB＞0，则执行步骤3.7；否则执行步骤3.9；
步骤3.9：计算B中的边需要向A延伸的长度δB；
步骤3.10：将B中所有需要向A延伸的边延伸δB；
步骤3.11：利用B分割A；
步骤3.12：删除A中面积较小的面；
步骤3.13：利用A分割B；
步骤3.14：删除B中面积较小的面；
步骤4：对所有的面进行缝合；
步骤5：对结构的几何模型划分有限元网格。
2.根据权利要求1所述的一种板式卫星结构共节点网格快速生成方法，其特征在于：
在步骤3.3中所述的“找出A中所有需要向B延伸的边及其数目NA”的方法为：现
有的有限元软件一般都提供二次开发接口使得用户能开发新的功能，因此利用相关的
函数接口通过二次开发并按照以下方法实现：1)获取A所有的边EA1,…,EAn和B的

\t一个顶点VB；2)对每条边EAj(j＝1,…,n)，假设其自然坐标表示为Pj(s)(0≤s≤1)，
均匀地获取Nj+1个点Pjk＝Pj(k/Nj)(k＝0,1,…,Nj)，其中Nj≥1000；3)
取dmax＝-999999,d＝0，然后按照初等解析几何中的公式
djk＝nBx(xjk-xB)+nBy(yjk-yB)+nBz(zjk-zB)依次计算每个点Pjk到B的有向距离，若
|djk|＞dmax，则dmax＝|djk|,d＝djk；其中(xjk,yjk,zjk)为Pjk在直角坐标系中的坐标，
(xB,yB,zB)为VB的坐标，(nBx,nBy,nBz)为B的单位法向量；4)取abd＝999999，对每
个点Pjk，若d·djk＞0且|djk|＜abd，则abd＝|djk|；5)若abd＜0.04，
则δmax＝abd+0.01；否则δmax＝0.05；6)NA＝0，对每条边EAj：若对k＝0,1,…,Nj均有|djk|≤δmax，则EAj标记为1，NA＝NA+1；否则EAj标记为0；在本步骤中获取A
所有的边和B的顶点均通过软件的接口实现，并且在有限元软件中边及曲线一般会采
用自然坐标进行表示，给定其自然坐标后通过接口得到其在直角坐标系下的坐标表
示；本步骤中所需要的二次开发接口函数通过查询有限元软件的相关帮助文档得到。
3.根据权利要求1所述的一种板式卫星结构共节点网格快速生成方法，其特征在于：
在步骤3.5中所述的“计算A中的边需要向B延伸的长度δA”的方法为：
δA＝δmax/sinθAB，其中δmax的值已经在步骤3.3中得到，θAB的值已经在步骤3.1中
得到。
4.根据权利要求1所述的一种板式卫星结构共节点网格快速生成方法，其特征在于：
在步骤3.6中所述的“将A中所有需要向B延伸的边延伸δA”的方法为：现有的有
限元软件一般都提供二次开发接口使得用户能开发新的功能，因此利用相关的函数接
口通过二次开发将A中所有需要向B延伸的边延伸δA；对一个面中的各条边进行延
伸所需要的二次开发接口函数通过查询有限元软件的相关帮助文档得到。
5.根据权利要求1所述的一种板式卫...

【专利技术属性】
技术研发人员：赵军鹏，刘腾达，付志方，王春洁，
申请(专利权)人：北京航空航天大学，
类型：发明
国别省市：北京;11