本发明专利技术公开了一种有限元模型中粗细网格单元的匹配方法及装置。匹配方法包括：获取有限元模型中的未匹配单元；在未匹配单元中选取任一单元作为待匹配单元，搜索与待匹配单元相邻的未匹配单元和已匹配单元，将待匹配单元以及与待匹配单元相邻的未匹配单元组成连通组；根据连通组的相邻已匹配单元确定连通组的类型，连通组的类型包括：延伸组和延展组；按照连通组的类型为待匹配单元匹配相应的粗网格单元。本发明专利技术可实现细网格模型中细网格单元与粗网格单元的全面匹配，解决了由于细网格模型与粗网格模型形状不同，而导致细网格单元不能全部与粗网格单元匹配的问题。部与粗网格单元匹配的问题。部与粗网格单元匹配的问题。
【技术实现步骤摘要】
一种有限元模型中粗细网格单元的匹配方法及装置
[0001]本专利技术涉及有限元模型分析
，特别涉及一种有限元模型中粗细网格单元的匹配方法及装置。
技术介绍
[0002]在有限元模型中，模型分析的精度与所用的有限元网格直接相关。有限元网格将CAD模型分割为很多单元，在这些单元上进行分析。随着网格的不断细化，单元变得越来越小，从而使模型分析的求解结果越来越接近真实解。
[0003]在模型分析中，基于粗糙网格完成求解后，开始进行网格细化。网格细化简而言之是一个用不断细化的网格来解析模型，并对不同网格下得到的结果进行比较的过程。通过对模型中一个或多个点上的物理场进行分析，或者在某些域或边界上对某个物理场进行积分计算，可以实现结果比较。
[0004]在网格细化的过程中需要对粗细网格进行映射，然而当前粗细网格映射仅能匹配表达相同的几何形状的粗细网格，而对于几何形状存在差异的模型，则不能完全匹配。例如，存在压弯型折角的细网格，由于细化后的网格描述了实际的几何形状，而粗网格是简化网格，对于几何形状是简化描述，与实际几何形状有较大差别，因此无法实现压弯型折角的粗细网格映射。
技术实现思路
[0005]为了解决上述现有技术中存在的至少一个技术问题，本专利技术实施例提供了一种有限元模型中粗细网格单元的匹配方法及装置。所述技术方案如下：
[0006]第一方面，提供了一种一种有限元模型中粗细网格单元的匹配方法，包括：
[0007]获取有限元模型中的未匹配单元，未匹配单元为所述有限元模型中无法与粗网格单元建立映射关系的细网格单元；
[0008]在未匹配单元中选取任一单元作为待匹配单元，搜索与所述待匹配单元相邻的未匹配单元和已匹配单元，将所述待匹配单元以及与所述待匹配单元相邻的未匹配单元组成连通组，与所述待匹配单元相邻的已匹配单元为所述连通组的相邻已匹配单元；
[0009]根据所述连通组的相邻已匹配单元确定所述连通组的类型，所述连通组的类型包括：延伸组和延展组；
[0010]按照所述连通组的类型为所述待匹配单元匹配相应的粗网格单元。
[0011]进一步地，所述搜索所述待匹配单元相邻的未匹配单元和已匹配单元，包括：
[0012]按照所述有限元模型中各未匹配单元的法向量对未匹配单元分组；
[0013]选取一单元组内任一未匹配单元作为所述待匹配单元，搜索所述待匹配单元相邻的未匹配单元和已匹配单元；
[0014]将搜索到的与所述待匹配单元相邻的未匹配单元添加至所述连通组中并将其从所述单元组中删除。
[0015]进一步地，所述按照所述有限元模型中各未匹配单元的法向量对未匹配单元分组，包括：
[0016]计算未匹配单元的法向量；
[0017]在相同坐标系中若未匹配单元的法向量的角度差符合同组阈值，则确定未匹配单元为同一单元组。
[0018]进一步地，所述根据所述连通组的相邻已匹配单元确定所述连通组的类型，包括：
[0019]计算所述连通组内各所述相邻已匹配单元的法向量；
[0020]若存在所述相邻已匹配单元的法向量与所述连通组内任一未匹配单元的法向量的角度差符合延伸组角度阈值，则所述连通组为所述延伸组；
[0021]若存在所述相邻已匹配单元的法向量与所述连通组内任一未匹配单元的法向量的夹角符合延展组角度阈值，则所述连通组为所述延展组。
[0022]进一步地，所述延展组角度阈值的确定，包括：
[0023]计算所述未匹配单元的形心与相邻的已匹配单元的形心的距离；
[0024]所述延展组角度阈值α为：其中，S为未匹配单元的形心与相邻的已匹配单元的形心的距离，R为最小曲率半径。
[0025]进一步地，所述按照所述连通组的类型为所述待匹配单元匹配所述粗网格单元，包括：
[0026]确定与所述延伸组中所述待匹配单元具有公共边的已匹配边界单元；
[0027]计算所述延伸组中所述待匹配单元与所述已匹配边界单元之间的距离，确定距离最小的所述已匹配边界单元对应的粗网格单元为所述待匹配单元匹配的粗网格单元。
[0028]进一步地，所述按照所述连通组的类型为所述待匹配单元匹配所述粗网格单元，包括：
[0029]计算经过所述延展组中待匹配单元形心的法向量，在法向量与所述待匹配单元的法向量相交的粗网格单元中，交点与所述待匹配单元的形心距离最小的粗网格单元为所述延展组中所述待匹配单元匹配的粗网格单元。
[0030]第二方面，提供了一种有限元模型中粗细网格单元的匹配装置，包括：
[0031]单元获取模块，用于获取有限元模型中的未匹配单元，未匹配单元为所述有限元模型中未与粗网格单元建立映射关系的细网格单元；
[0032]单元搜索模块，用于在未匹配单元中选取任一单元作为待匹配单元，搜索与所述待匹配单元相邻的未匹配单元和已匹配单元，将所述待匹配单元以及与所述待匹配单元相邻的未匹配单元组成连通组，与所述待匹配单元相邻的以匹配单元未所述连通组的相邻已匹配单元；
[0033]类型确定模块，用于根据所述连通组的相邻已匹配单元确定所述连通组的类型，所述连通组的类型包括：延伸组和延展组；
[0034]匹配模块，用于按照所述连通组的类型为所述待匹配单元匹配相应的粗网格单元。
[0035]进一步地，所述单元搜索模块，包括：
[0036]分组模块，用于按照所述有限元模型中各未匹配单元的法向量对未匹配单元分
组；
[0037]搜索模块，用于选取一单元组内任一未匹配单元作为所述待匹配单元，搜索所述待匹配单元相邻的未匹配单元和已匹配单元；
[0038]以及，将搜索到的与所述待匹配单元相邻的未匹配单元添加至所述连通组中并将其从所述单元组中删除。
[0039]进一步地，所述分组模块，具体用于：
[0040]用于计算未匹配单元的法向量；
[0041]在相同坐标系中若未匹配单元的法向量的角度差符合同组阈值，则确定未匹配单元为同一单元组。
[0042]进一步地，所述类型确定模块，包括：
[0043]计算模块，用于计算所述连通组内各所述相邻已匹配单元的法向量；
[0044]延伸组确定模块，用于若存在所述相邻已匹配单元的法向量与所述连通组内任一未匹配单元的法向量的角度差符合延伸组角度阈值，则所述连通组为所述延伸组；
[0045]延展组确定模块，用于若存在所述相邻已匹配单元的法向量与所述连通组内任一未匹配单元的法向量的夹角符合延展组角度阈值，则所述连通组为所述延展组。
[0046]进一步地，所述类型确定模块，还包括：
[0047]延展组角度阈值确定模块，具体用于：
[0048]计算所述未匹配单元的形心与相邻的已匹配单元的形心的距离；
[0049]所述延展组角度阈值α为：其中，S为未匹配单元的形心与相邻的已匹配单元的形心的距离，R为最小曲率半径。
[0050]进一步地，所述匹配本文档来自技高网...
【技术保护点】
【技术特征摘要】
1.一种有限元模型中粗细网格单元的匹配方法，其特征在于，包括：获取有限元模型中的未匹配单元，未匹配单元为所述有限元模型中无法与粗网格单元建立映射关系的细网格单元；在未匹配单元中选取任一单元作为待匹配单元，搜索与所述待匹配单元相邻的未匹配单元和已匹配单元，将所述待匹配单元以及与所述待匹配单元相邻的未匹配单元组成连通组，与所述待匹配单元相邻的已匹配单元为所述连通组的相邻已匹配单元；根据所述连通组的相邻已匹配单元确定所述连通组的类型，所述连通组的类型包括：延伸组和延展组；按照所述连通组的类型为所述待匹配单元匹配相应的粗网格单元。2.如权利要求1所述的方法，其特征在于，所述搜索所述待匹配单元相邻的未匹配单元和已匹配单元，包括：按照所述有限元模型中各未匹配单元的法向量对未匹配单元分组；选取一单元组内任一未匹配单元作为所述待匹配单元，搜索所述待匹配单元相邻的未匹配单元和已匹配单元；将搜索到的与所述待匹配单元相邻的未匹配单元添加至所述连通组中并将其从所述单元组中删除。3.如权利要求2所述的方法，其特征在于，所述按照所述有限元模型中各未匹配单元的法向量对未匹配单元分组，包括：计算未匹配单元的法向量；在相同坐标系中若未匹配单元的法向量的角度差符合同组阈值，则确定未匹配单元为同一单元组。4.如权利要求1所述的方法，其特征在于，所述根据所述连通组的相邻已匹配单元确定所述连通组的类型，包括：计算所述连通组内各所述相邻已匹配单元的法向量；若存在所述相邻已匹配单元的法向量与所述连通组内任一未匹配单元的法向量的角度差符合延伸组角度阈值，则所述连通组为所述延伸组；若存在所述相邻已匹配单元的法向量与所述连通组内任一未匹配单元的法向量的夹角符合延展组角度阈值，则所述连通组为所述延展组。5.如权利要求4所述的方法，其特征在于，所述延展组角度阈值的确定，包括：计算所述未匹配单元的形心与相邻的已匹配单元的形心的距离；所述延展组角度阈值α为：其中，S为未匹配单元的形心...
【专利技术属性】
技术研发人员：胡丰梁，孙安林，胡欣，刘玉川，张道坤，吴薇，王梦泽，陈永念，孟凡冲，李明，王丽荣，孙磊，王思源，
申请(专利权)人：中国船级社，
类型：发明
国别省市：
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。