CN113191061A - 一种基于曲面特征识别的有限元网格变换方法 - Google Patents
一种基于曲面特征识别的有限元网格变换方法 Download PDFInfo
- Publication number
- CN113191061A CN113191061A CN202110707275.9A CN202110707275A CN113191061A CN 113191061 A CN113191061 A CN 113191061A CN 202110707275 A CN202110707275 A CN 202110707275A CN 113191061 A CN113191061 A CN 113191061A
- Authority
- CN
- China
- Prior art keywords
- geometric
- subjected
- transformed
- feature
- geometric body
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
- G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/15—Vehicle, aircraft or watercraft design
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Automation & Control Theory (AREA)
- Aviation & Aerospace Engineering (AREA)
- Image Analysis (AREA)
Abstract
本发明提出了一种基于曲面特征识别的有限元网格变换方法,通过提出待进行变换的结合体以及待进行特征识别的几何体的外部曲面,通过曲面特征进行相似度的判断,对于判断为一致的待进行变换的结合体以及待进行特征识别的几何体,通过设置相对应的几何特征点,并通过几何特征点完成有限元网格的变换,实现自动化的识别和变化。本发明方法能够快速为一件多用零件自动划分网格形式完全一致的有限元模型,相比传统的手动建模方法能够大幅缩减建模过程,有效避免了零件遗漏、网格不匹配等常见的人因错误。
Description
技术领域
本发明属于有限元结构建模技术领域,具体地说,一种基于曲面特征识别的有限元网格变换方法。
背景技术
在航空、航天、车辆、船舶、机械等专业领域,同一结构零件在整个系统部件中时常被多次使用,这种设计能够降低系统的零部件数量和CAE(Computer Aided Engineering,计算机辅助工程)结构数模的内存占用,并且有利于工艺加工和装配,这类零件通常被称为一件多用零件。
结构强度设计规范要求一件多用零件的有限元网格必须完全一致。目前对于一件多用零件的有限元建模方法通常为:对初始零件划分有限元网格,人工查找与该零件几何特征相一致的零件,分别在前后两个零件上按照一定的规则选取几何特征点,以几何特征点为参考点手动将之前划分好的网格备份并移动到新的位置。这一过程操作繁多且易于出错,一件多用零件查找不全或者几何特征点选取错误等问题很难避免,特别是在结构设计发生迭代更改时,一件多用零件有限元网格的更新过程将变得异常繁琐,时间和人力成本巨大。
发明内容
本发明针对现有技术的上述缺陷,提出了一种基于曲面特征识别的有限元网格变换方法,通过提出待进行变换的结合体以及待进行特征识别的几何体的外部曲面,通过曲面特征进行相似度的判断,对于判断为一致的待进行变换的结合体以及待进行特征识别的几何体,通过设置相对应的几何特征点,并通过几何特征点完成有限元网格的变换,实现自动化的识别和变化。
本发明具体实现内容如下:
本发明提出了一种基于曲面特征识别的有限元网格变换方法,用于将有限元网格中待进行变换的几何体变换到待进行特征识别的几何体,所述有限元网格变换方法包括以下步骤:
步骤S1:在有限元网格中对待进行变换的几何体获取并进行分组标记;
步骤S2:对待进行变换的几何体进行外部曲面的提取,并在待进行特征识别的几何体上进行几何特征点的标记;
步骤S3:提取待进行变换的几何体的所有外部曲面;
步骤S4:汇总所有待进行特征识别的几何体,建立循环列表,依照循环列表依次筛选出一个待进行特征识别的几何体,对筛选出的待进行特征识别的几何体的每一个外部曲面进行提取;
步骤S5:计算待进行变换的几何体的外部曲面和筛选出的待进行特征识别的几何体的外部曲面之间的相似度系数;
步骤S6:根据步骤S5的计算结果汇总得到待进行变换的几何体和筛选出的待进行特征识别的几何体的之间的相似度矩阵;并根据相似度矩阵判断待进行变换的几何体和筛选出的待进行特征识别的几何体是否完全一致;
步骤S7:对于判断为完全一致的待进行变换的几何体和筛选出的待进行特征识别的几何体,提取出在待进行变换的几何体的外部曲面上标记的几何特征点,并在筛选出的待进行特征识别的几何体上设置对应的几何特征点,然后执行步骤S8的操作;对于判断为不完全一致的待进行变换的几何体和筛选出的待进行特征识别的几何体,回到步骤S4,筛选循环列表中的下一个待进行特征识别的几何体进行步骤S4及之后的操作,直至处理完整个循环列表;
步骤S8:以待进行变换的几何体上的几何特征点为初始参考点,以待进行特征识别的几何体上的几何特征点为终止参考点,将待进行变换的几何体在有限元网格上进行变换。
为了更好地实现本发明,进一步地,所述步骤S1具体包括:
步骤S1:获取待进行变换的有限元网格F_0对应的待进行变换的几何体Part_0所在分组并记为Comp_0。
为了更好地实现本发明,进一步地,所述步骤S2具体包括:
为了更好地实现本发明,进一步地,所述步骤S3具体包括:
步骤S3:将待进行特征识别的几何体记为Part_x (x=1,2,⋯,k),将待进行特征识别的几何体Part_x所在的分组记为Comp_x (x=1,2,⋯,k),其中k为待进行特征识别的几何体的总数。
为了更好地实现本发明,进一步地,所述步骤S4具体包括:
为了更好地实现本发明,进一步地,所述步骤S5具体包括:
步骤S5:将待进行特征识别的几何体Part_x的每一个曲面分别和待进行变换的Part_0的每一个曲面进行对比并将对比结果记录为相似度系数;如果曲面和曲面几何特征完全一致,则相似度系数,否则相似度系数。
为了更好地实现本发明,进一步地,所述步骤S6具体包括:
为了更好地实现本发明,进一步地,所述步骤S7具体包括:
步骤S7:对判断为和待进行变换的几何体Part_0完全一致的待进行特征识别的几何体Part_x,设置和几何特征点对应的几何特征点,然后执行步骤S8;对判断为不一致的几何体Part_x,令x=x+1,然后执行步骤S4。
为了更好地实现本发明,进一步地,所述步骤S8具体包括:
本发明与现有技术相比具有以下优点及有益效果:
1. 本发明提供了一种基于曲面特征识别的有限元网格变换方法,该方法能够快速为一件多用零件自动划分网格形式完全一致的有限元模型,相比传统的手动建模方法能够大幅缩减建模过程,有效避免了零件遗漏、网格不匹配等常见的人因错误;
2. 本发明提供的一种基于曲面特征识别的有限元网格变换方法是一种适用于一件多用零件有限元建模的自动化快速建模方法,对于重复零部件使用广泛以及设计迭代频繁的结构设计领域,能够大幅缩短设计周期,保证设计质量。
附图说明
图1是本发明的方法流程框图;
图2是典型的使用一件多用零件的零组件示意图;
图3显示的是图2所示零件Part_0的面网格有限元模型;
图4显示的是按特定规则标记的图2所示零件Part_0的3个几何特征点;
图5显示的是采用本发明方法变换得到的有限元网格模型。
具体实施方式
为了更清楚地说明本发明实施例的技术方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,应当理解,所描述的实施例仅仅是本发明的一部分实施例,而不是全部的实施例,因此不应被看作是对保护范围的限定。基于本发明中的实施例,本领域普通技术工作人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
在本发明的描述中,需要说明的是,除非另有明确的规定和限定,术语“设置”、“相连”、“连接”应做广义理解,例如,可以是固定连接,也可以是可拆卸连接,或一体地连接;可以是机械连接,也可以是电连接;也可以是直接相连,也可以是通过中间媒介间接相连,可以是两个元件内部的连通。对于本领域的普通技术人员而言,可以具体情况理解上述术语在本发明中的具体含义。
实施例1:
本实施例提出了一种基于曲面特征识别的有限元网格变换方法,用于将有限元网格中待进行变换的几何体变换到待进行特征识别的几何体,如图1所示,所述有限元网格变换方法包括以下步骤:
步骤S1:在有限元网格中对待进行变换的几何体获取并进行分组标记;
步骤S2:对待进行变换的几何体进行外部曲面的提取,并在待进行特征识别的几何体上进行几何特征点的标记;
步骤S3:提取待进行变换的几何体的所有外部曲面;
步骤S4:汇总所有待进行特征识别的几何体,建立循环列表,依照循环列表依次筛选出一个待进行特征识别的几何体,对筛选出的待进行特征识别的几何体的每一个外部曲面进行提取;
步骤S5:计算待进行变换的几何体的外部曲面和筛选出的待进行特征识别的几何体的外部曲面之间的相似度系数;
步骤S6:根据步骤S5的计算结果汇总得到待进行变换的几何体和筛选出的待进行特征识别的几何体的之间的相似度矩阵;并根据相似度矩阵判断待进行变换的几何体和筛选出的待进行特征识别的几何体是否完全一致;
步骤S7:对于判断为完全一致的待进行变换的几何体和筛选出的待进行特征识别的几何体,提取出在待进行变换的几何体的外部曲面上标记的几何特征点,并在筛选出的待进行特征识别的几何体上设置对应的几何特征点,然后执行步骤S8的操作;对于判断为不完全一致的待进行变换的几何体和筛选出的待进行特征识别的几何体,回到步骤S4,筛选循环列表中的下一个待进行特征识别的几何体进行步骤S4及之后的操作,直至处理完整个循环列表;
步骤S8:以待进行变换的几何体上的几何特征点为初始参考点,以待进行特征识别的几何体上的几何特征点为终止参考点,将待进行变换的几何体在有限元网格上进行变换。
实施例2:
本实施例在上述实施例1的基础上,为了更好地实现本发明,进一步地,如图1所示,所述方法具体包括:
步骤S1:获取待进行变换的有限元网格F_0对应的待进行变换的几何体Part_0所在分组并记为Comp_0;
步骤S3:将待进行特征识别的几何体记为Part_x (x=1,2,⋯,k),将待进行特征识别的几何体Part_x所在的分组记为Comp_x (x=1,2,⋯,k),其中k为待进行特征识别的几何体的总数。
步骤S5:将待进行特征识别的几何体Part_x的每一个曲面分别和待进行变换的Part_0的每一个曲面进行对比并将对比结果记录为相似度系数;如果曲面和曲面几何特征完全一致,则相似度系数,否则相似度系数;
步骤S7:对判断为和待进行变换的几何体Part_0完全一致的待进行特征识别的几何体Part_x,设置和几何特征点对应的几何特征点,然后执行步骤S8;对判断为不一致的几何体Part_x,令x=x+1,然后执行步骤S4;
本实施例的其他部分与上述实施例1相同,故不再赘述。
实施例3:
本实施例在上述实施例1-2任一项的基础上,如图1、图2、图3、图4、图5所示,本方法采用Tcl语言进行程序开发,在HyperMesh软件平台下运行。下面结合附图,对本发明的具体实施方式作进一步的详细描述。
根据本发明的一个方面,提出了一种基于曲面特征识别的有限元网格变换方法,所述方法以零件几何体的所有外部曲面为对象,通过识别曲面特征参数来计算零件的相似度矩阵,并以此判断待识别的零件几何体与原零件几何体特征是否一致,对于曲面特征完全一致的零件,通过识别相应的几何特征点,以此为参考完成零件几何体有限元网格的跟随变换。具体方法流程如图1所示。
图2显示的是一种典型的使用一件多用零件的零组件示意图,该零组件由6个零件构成,分别记为Part_x (x=0,1,2,⋯,5),其中Part_x (x=0,1,2,3) 4个几何体为一件多用零件,几何特征完全一致,Part_4和 Part_5与其余4个零件几何特征不同。初始有限元模型是为Part_0而划分的面网格单元,具体如图3所示。下面结合此实例来说明本发明的具体实施方式:
步骤S1:获取待进行变换的有限元网格F_0对应的待进行变换的几何体Part_0所在分组并记为Comp_0 。
步骤S2:提取待进行变换的几何体Part_0的所有外部曲面,待进行变换的几何体Part_0共有19个外部曲面;按照一定的规则标记3个待进行变换的几何体Part_0中的几何特征点,标记点记录为,如图4所示。
步骤S3:选择需要进行特征识别的待进行特征识别的几何体Part_x (x=1,2,⋯,5),将待进行特征识别的几何体Part_x所在的分组记为Comp_x (x=1,2,⋯,5)。
判定待进行特征识别的几何体Part_x和待进行变换的几何体Part_0几何特征完全一致;否则,判定待进行特征识别的几何体Part_x和待进行变换的几何体Part_0几何特征不一致。
步骤S8:以为初始参考点,为终止参考点,将有限元网格F_0从Comp_0变换至Comp_x,变换得到的有限元网格模型如图5所示,可以看到,与Part_0几何特征完全一致的Part_x (x=1,2,3)被识别了出来且自动生成了有限元网格,而Part_x (x=4,5)与Part_0几何特征不一致,则没有生成有限元网格。
本发明方法能够快速为一件多用零件自动划分网格形式完全一致的有限元模型,相比传统的手动建模方法能够大幅缩减建模过程,有效避免了零件遗漏、网格不匹配等常见的人因错误。为了使本发明的具体实施方法说明更为简洁明了,本实例选取的零组件零件数目相对较少,结构形式也较为简单,实际工程中应用的一件多用零件数目更大,复杂程度也更高,本发明方法带来的效率提升和设计质量提升效果将更为显著。
本实施例的其他部分与上述实施例1-2任一项相同,故不再赘述。
以上所述,仅是本发明的较佳实施例,并非对本发明做任何形式上的限制,凡是依据本发明的技术实质对以上实施例所作的任何简单修改、等同变化,均落入本发明的保护范围之内。
Claims (9)
1.一种基于曲面特征识别的有限元网格变换方法,用于将有限元网格中待进行变换的几何体变换到待进行特征识别的几何体,其特征在于,包括以下步骤:
步骤S1:在有限元网格中对待进行变换的几何体获取并进行分组标记;
步骤S2:对待进行变换的几何体进行外部曲面的提取,并在待进行特征识别的几何体上进行几何特征点的标记;
步骤S3:提取待进行变换的几何体的所有外部曲面;
步骤S4:汇总所有待进行特征识别的几何体,建立循环列表,依照循环列表依次筛选出一个待进行特征识别的几何体,对筛选出的待进行特征识别的几何体的每一个外部曲面进行提取;
步骤S5:计算待进行变换的几何体的外部曲面和筛选出的待进行特征识别的几何体的外部曲面之间的相似度系数;
步骤S6:根据步骤S5的计算结果汇总得到待进行变换的几何体和筛选出的待进行特征识别的几何体的之间的相似度矩阵;并根据相似度矩阵判断待进行变换的几何体和筛选出的待进行特征识别的几何体是否完全一致;
步骤S7:对于判断为完全一致的待进行变换的几何体和筛选出的待进行特征识别的几何体,提取出在待进行变换的几何体的外部曲面上标记的几何特征点,并在筛选出的待进行特征识别的几何体上设置对应的几何特征点,然后执行步骤S8的操作;对于判断为不完全一致的待进行变换的几何体和筛选出的待进行特征识别的几何体,回到步骤S4,筛选循环列表中的下一个待进行特征识别的几何体进行步骤S4及之后的操作,直至处理完整个循环列表;
步骤S8:以待进行变换的几何体上的几何特征点为初始参考点,以待进行特征识别的几何体上的几何特征点为终止参考点,将待进行变换的几何体在有限元网格上进行变换。
2.如权利要求1所述的一种基于曲面特征识别的有限元网格变换方法,其特征在于,所述步骤S1具体包括:
步骤S1.1:获取待进行变换的有限元网格F_0对应的待进行变换的几何体Part_0所在分组并记为Comp_0。
4.如权利要求3所述的一种基于曲面特征识别的有限元网格变换方法,其特征在于,所述步骤S3具体包括:
步骤S3.1:将待进行特征识别的几何体记为Part_x (x=1,2,⋯,k),将待进行特征识别的几何体Part_x所在的分组记为Comp_x (x=1,2,⋯,k),其中k为待进行特征识别的几何体的总数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110707275.9A CN113191061B (zh) | 2021-06-25 | 2021-06-25 | 一种基于曲面特征识别的有限元网格变换方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110707275.9A CN113191061B (zh) | 2021-06-25 | 2021-06-25 | 一种基于曲面特征识别的有限元网格变换方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113191061A true CN113191061A (zh) | 2021-07-30 |
CN113191061B CN113191061B (zh) | 2021-10-08 |
Family
ID=76976836
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110707275.9A Active CN113191061B (zh) | 2021-06-25 | 2021-06-25 | 一种基于曲面特征识别的有限元网格变换方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113191061B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101702185A (zh) * | 2009-11-05 | 2010-05-05 | 西北工业大学 | 一种飞机框肋类钣金零件制造工件模型设计方法 |
US20140324400A1 (en) * | 2013-04-30 | 2014-10-30 | Marquette University | Gesture-Based Visualization System for Biomedical Imaging and Scientific Datasets |
CN107529641A (zh) * | 2016-06-24 | 2018-01-02 | 波音公司 | 飞行器机翼的前缘翼肋的建模与分析 |
CN108053480A (zh) * | 2017-12-08 | 2018-05-18 | 东华大学 | 基于逆向工程技术的三维全尺度着装人体网格构建方法 |
-
2021
- 2021-06-25 CN CN202110707275.9A patent/CN113191061B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101702185A (zh) * | 2009-11-05 | 2010-05-05 | 西北工业大学 | 一种飞机框肋类钣金零件制造工件模型设计方法 |
US20140324400A1 (en) * | 2013-04-30 | 2014-10-30 | Marquette University | Gesture-Based Visualization System for Biomedical Imaging and Scientific Datasets |
CN107529641A (zh) * | 2016-06-24 | 2018-01-02 | 波音公司 | 飞行器机翼的前缘翼肋的建模与分析 |
CN108053480A (zh) * | 2017-12-08 | 2018-05-18 | 东华大学 | 基于逆向工程技术的三维全尺度着装人体网格构建方法 |
Non-Patent Citations (2)
Title |
---|
CHRISTOPH LEHRENFELD: "High order unfitted finite element methods on level set domains using isoparametric mappings", 《COMPUTER METHODS IN APPLIED MECHANICS AND ENGINEERING》 * |
楼航迪 等: "面部软组织三维有限元建模和变换", 《系统仿真学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN113191061B (zh) | 2021-10-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111832102B (zh) | 一种高维随机场条件下的新型复合材料结构优化设计方法 | |
CN110084221B (zh) | 一种基于深度学习的带中继监督的序列化人脸关键点检测方法 | |
CN106960088B (zh) | 面向工序更改的三维工序模型自动更新方法 | |
CN111914480B (zh) | 一种基于点云语义分割的加工特征智能识别方法 | |
US7991595B2 (en) | Adaptive refinement tools for tetrahedral unstructured grids | |
CN112257197B (zh) | 一种大型铸锻件微缺陷工作应力评估方法 | |
CN101751499A (zh) | 用于三维船舶建模零件的自动装配方法 | |
CN113191061B (zh) | 一种基于曲面特征识别的有限元网格变换方法 | |
CN113642129B (zh) | 一种工件校正载荷快速施加及变形有限元获取方法 | |
CN115274006A (zh) | 基于蛋白质-配体复合物结构信息的亲和力预测方法 | |
CN112765876B (zh) | 机械设计方案推荐方法、装置及电子设备 | |
Song et al. | Automatic recognition and suppression of holes on mold bases for finite element applications | |
CN115587480A (zh) | 数字化仿真方法以及数字化仿真装置 | |
Ando et al. | Aerodynamic performance evaluation system at the early concept stage of automotive styling development based on CFD | |
CN108089136B (zh) | 一种燃料电池电堆极化曲线测试数据自动切片方法 | |
CN113361630B (zh) | 一种基于曲面特征对比的零件异同自动识别方法 | |
KR100993297B1 (ko) | Catia를 이용한 패널 코드 전처리 방법 | |
CN117252993B (zh) | 特征点提取算法的验证方法、装置、电子设备及存储介质 | |
Lai et al. | Recognition and simplification of holes in CAD models of an injection mold for mold flow analysis | |
CN111062167A (zh) | 一种用于叶片精锻的叶片-工装系统设计的分析方法 | |
CN112597727B (zh) | 一种新型快速高效的滤波器小样本建模及优化方法 | |
CN113705034B (zh) | 一种仿真结果处理方法和装置 | |
EP1083499B1 (en) | Method for thermal fluid analysis and thermal stress analysis of fluid machine | |
CN118194428A (zh) | 一种舰船目标易损性模型的快速建模方法 | |
CN117854652B (zh) | 焊缝微观组织晶体塑性有限元模型的建模方法及系统 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |