CN110059332A - 包含晶格结构的3d部件的3d有限元网格的设计 - Google Patents

包含晶格结构的3d部件的3d有限元网格的设计 Download PDF

Info

Publication number
CN110059332A
CN110059332A CN201811580554.8A CN201811580554A CN110059332A CN 110059332 A CN110059332 A CN 110059332A CN 201811580554 A CN201811580554 A CN 201811580554A CN 110059332 A CN110059332 A CN 110059332A
Authority
CN
China
Prior art keywords
unit
group
component
entity
finite element
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.)
Pending
Application number
CN201811580554.8A
Other languages
English (en)
Inventor
M-P·施密特
J-M·拉绍姆
D·L·邦纳
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Dassault Systemes SE
Dassault Systemes of America Corp
Original Assignee
Dassault Systemes of America Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Dassault Systemes of America Corp filed Critical Dassault Systemes of America Corp
Publication of CN110059332A publication Critical patent/CN110059332A/zh
Pending legal-status Critical Current

Links

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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
    • G06T17/20—Finite element generation, e.g. wire-frame surface description, tesselation
    • B—PERFORMING OPERATIONS; TRANSPORTING
    • B33—ADDITIVE MANUFACTURING TECHNOLOGY
    • B33Y—ADDITIVE MANUFACTURING, i.e. MANUFACTURING OF THREE-DIMENSIONAL [3-D] OBJECTS BY ADDITIVE DEPOSITION, ADDITIVE AGGLOMERATION OR ADDITIVE LAYERING, e.g. BY 3-D PRINTING, STEREOLITHOGRAPHY OR SELECTIVE LASER SINTERING
    • B33Y50/00—Data acquisition or data processing for additive manufacturing
    • G—PHYSICS
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B19/00—Programme-control systems
    • G05B19/02—Programme-control systems electric
    • G05B19/18—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form
    • G05B19/4097—Numerical control [NC], i.e. automatically operating machines, in particular machine tools, e.g. in a manufacturing environment, so as to execute positioning, movement or co-ordinated operations by means of programme data in numerical form characterised by using design data to control NC machines, e.g. CAD/CAM
    • G05B19/4099—Surface or curve machining, making 3D objects, e.g. desktop manufacturing
    • 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
    • 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
    • G05—CONTROLLING; REGULATING
    • G05B—CONTROL OR REGULATING SYSTEMS IN GENERAL; FUNCTIONAL ELEMENTS OF SUCH SYSTEMS; MONITORING OR TESTING ARRANGEMENTS FOR SUCH SYSTEMS OR ELEMENTS
    • G05B2219/00—Program-control systems
    • G05B2219/30—Nc systems
    • G05B2219/49—Nc machine tool, till multiple
    • G05B2219/49023—3-D printing, layer of powder, add drops of binder in layer, new powder

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Manufacturing & Machinery (AREA)
  • General Engineering & Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Software Systems (AREA)
  • Computer Graphics (AREA)
  • Human Computer Interaction (AREA)
  • Automation & Control Theory (AREA)
  • Chemical & Material Sciences (AREA)
  • Materials Engineering (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)

Abstract

本发明尤其涉及一种用于设计包含晶格结构的三维(3D)部件的3D有限元网格的计算机实现的方法。该方法包括提供3D部件的实体表示；提供单元的规则拼接，每个单元包括有限元网格；将单元的规则拼接与3D部件的实体表示相叠加；将单元划分(S40)成两组，其中，第一组单元中的每一个单元与3D部件的实体表示相接触，以及第二组单元中没有一个单元与3D部件的实体表示相接触。该方法还包括：计算第一组单元与3D部件的实体表示的布尔并集，所述布尔并集形成体积；对所计算的布尔并集的体积进行有限元网格划分，同时保留第一组单元的与第二组单元共享的面的集合；以及对第二组单元中的单元的有限元网格和所计算的布尔并集的经网格划分的体积进行合并。

Description

包含晶格结构的3D部件的3D有限元网格的设计
技术领域
本发明涉及计算机程序和***领域，并更具体地，涉及用于设计包含晶格结构(lattice structure)的3D部件的3D有限元网格(mesh)的方法、***和程序。
背景技术
市场上存在用于对象的设计、工程规划和制造的多种***和程序。CAD是“计算机辅助设计”的缩略词，例如，其涉及用于设计对象的软件解决方案。CAE是“计算机辅助工程”的缩略词，例如，其涉及用于仿真未来产品的物理行为的软件解决方案。CAM是“计算机辅助制造”的缩略词，例如，其涉及用于规定制造过程和操作的软件解决方案。在这些计算机辅助设计***中，图形用户界面针对技术的高效性起到了重要作用。这些技术可以被嵌入到产品生命周期管理(PLM)***中。PLM指的是在扩展型企业概念上的一种商业策略，其用来帮助各个公司针对从产品概念提出到产品生命周期结束的产品研发，来共享产品数据、应用公共过程、利用团体知识。由达索***(Dassault Systèmes)(使用商标CATIA、ENOVIA和DELMIA)提供的PLM解决方案提供了工程中心(其对产品工程知识进行组织)、制造中心(其管理制造工程知识)以及企业中心(其使得企业集成和连接到工程中心和制造中心)。所有这些使得该***产生了将产品、过程、资源链接在一起的开放式对象模型，实现了动态的基于知识的产品创建和决策支持，其推动了优化的产品定义、制造准备、生产以及服务。
CAD/CAE***用于实现物理***的仿真，该物理***可以属于但不限于航空、汽车、海事以及城市工程。尤其是，这些***用于使用有限元方法(FEM)来执行对该物理***上的物理现象的仿真。这种仿真允许预测使用中的物理***的性能(热学性能、机器性能或任何其他物理性能)并对设计进行优化以改善所述性能。其还允许预测部件的可制造性。在3D打印的情况下，可制造性可能是成败的关键，例如，对具有诸如在直接金属打印期间部件中的热分布之类的现象的打印过程进行仿真。事实上，诸如残余应力之类的现象在3D打印中可能存在严重问题，并且可能使得物理***没有资格被用于诸如航空等关键应用。因此，能够仿真在打印过程中出现的现象是非常重要的。
晶格结构经常被用于设计所建模的物理***的一个或多个部件。晶格结构表示在建模物理***的部件中重复交替的实体和空心的图案。这种晶格图案的目的可以是调整部件的重量、机械阻力、传热或其他属性。因此，晶格结构改善了部件的整体性能；尤其是，由于晶格结构，机械应力能够在整个部件上更好地分布。由于对3D打印的兴趣的提高，晶格实体(lattice solid)的设计和分析变得非常重要。3D打印机对形状的复杂性几乎没有任何限制。因此，设计人员可以自由地创建以前无法制造的非常复杂的形状。特别是，3D打印机可以容易地打印晶格实体。
为了用有限元分析(FEA)对建模的物理***进行分析，必须将其分解成多个有限元。将建模的物理***的一部分分解成元素的过程称为网格划分。这种分解越精细，有限元分析就越精确。当前，看到单个机械部件分解成几百万个元素并不罕见。在有限元分析开始时，网格的布局是手工创建的，但这是一个痛苦的过程，是实现物理仿真真正效益的主要瓶颈，因为每个新的或改善的设计都涉及重新创建这个有限元网格的任务。近20年来，自动化3D网格划分逐渐取代手工网格划分，使得有限元分析的整体成本大大降低。在有限元网格划分中，两种主要元素形状是占大多数的：具有4条边的四面体(Tet)以及具有6条边的六面体(Hex)。对于大多数类型的问题，有限元分析使用六面体元素将更有效的。因此，当用手工进行网格划分时，工程师将引入大比例的这些六面体元素。然而，3D六面体网格划分已经被证明是很难实现自动化的。因此，目前使用的大多数自动化3D网格划分是基于四面体的。由于自动化3D网格划分带来的整体生产力的巨大提高，计算效率的牺牲对于用户而言是可接受的。
通常，自动化3D网格划分试图用有限元填充一个体积，而不对体积的性质作出太多的假设。由于要填充的体积的所有潜在形状变化，该算法必然是复杂且执行时间长的。这些算法的执行时间中最重要的因素是要填充的形状的边界复杂度。
在晶格结构的情况下，自动化3D网格划分会遇到问题。第一个问题是晶格实体像海绵一样是微孔的。这意味着它们具有特别高的边界复杂性，并且非常不适用于自动化3D网格划分，3D网格划分对于边界的情况复杂性会增加。对于具有数万或数十万个单位单元(unit cells)的晶格，边界变得过于复杂，使得使用自动化3D网格划分变得不切实际。另一个问题是，自动化3D网格划分不会在每个晶格单元(lattice cell)中产生相同的网格，从而产生不均匀的结果和局部伪影，这在分析中造成了噪声和不精确性。为了避免这种伪影，用户可能会选择生成更精细的网格，但是这构成了恶性循环，因为其将加剧上面提到的第一个问题的复杂度增加问题。此外，为了选择有利于最有效计算的元素类型，3D网格划分将不会利用晶格实体的重复性质的优势。特别是，自动化3D网格划分将通常用四面体元素填充实体，而有限元分析技术(在大多数情况下)在六面体占大多数的网格(也就是其中大多数元素是六面体的网格)上表现最好。
在这种情况下，仍然存在对于包含晶格结构的3D部件的3D有限元网格的3D结构的改善设计的需求。
发明内容
因此，提供了一种用于设计包含晶格结构的三维(3D)部件的3D有限元网格的计算机实现的方法。该方法包括：提供3D部件的实体表示(solid representation)；提供单元的规则拼接(tiling)，每个单元包括有限元网格。该方法还包括：将所述单元的规则拼接与3D部件的实体表示相叠加；将所述单元划分成两组：(i)第一组单元，其中每一个单元与3D部件的实体表示相接触，以及(ii)第二组单元，其中没有一个单元与3D部件的实体表示相接触。该方法还包括：计算第一组单元与3D部件的实体表示的布尔并集，该布尔并集形成体积；对所计算的布尔并集的所述体积进行有限元网格划分，同时保留第一组单元的与第二组单元共享的面的集合；以及对第二组单元中的单元的有限元网格和所计算的布尔并集的经网格划分的体积进行合并。
该方法可以包括以下一项或几项：
-所述规则拼接的单元在拓扑上相同；
-对所计算的第一组单元与3D部件的实体表示的布尔并集进行修剪；
-保留第一组单元的与第二组单元共享的面的集合进一步包括：存储第一组单元和第二组单元的交界面，所述交界面包括第一组单元的与第二组单元共享的面的集合；
-对所计算的布尔并集的所述体积进行有限元网格划分的步骤进一步包括：将所存储的交界面合并到所计算的布尔并集上；
-将高级前沿(advanced front)网格生成器应用于所计算的布尔并集的经网格划分的体积；
-将所合并的交界面的一个或多个自由边用作被应用于所合并的交界面的高级前沿网格生成器的初始前沿；
-从被应用于所合并的交界面的高级前沿网格生成器保留所合并的交界面的一个或多个自由边；
-在划分之后，移除第一组单元中的单元的内部晶格结构；
-在划分之后，实例化第二组单元中每个单元的有限元网格。
还提供了一种用于增材制造的计算机实现的方法，该方法包括：提供根据上述方法设计的包含晶格结构的三维(3D)部件的3D有限元网格；以及对该3D部件进行增材制造。
该用于增材制造的方法可以包括以下一项或多项：
-提供步骤包括：对包含晶格结构的三维(3D)部件的3D有限元网格进行仿真；
-对包含晶格结构的三维(3D)部件的3D有限元网格进行仿真的步骤包括：仿真在增材制造过程中发生的至少一个物理现象；
-所仿真的所述至少一个物理现象是从如下中选择的：在直接金属打印过程中3D部件中的热分布、3D部件中的残余应力。
还提供了一种包括用于执行所述方法的指令的计算机程序。
还提供了一种其上记录有所述计算机程序的计算机可读存储介质。
还提供了一种***，其包括耦合到存储器的处理器，所述存储器中记录有所述计算机程序。
还提供了一种通过上述用于增材制造的方法可以获得的所制作的3D部件。
附图说明
现在将借助于非限定性实施例以及参考附图来描述本发明的实施例，其中：
图1示出了方法的示例的流程图；
图2示出了CAD/CAE***的图形用户界面的示例；
图3示出了CAD/CAE***的示例；
图4-12示出了图1的流程图的步骤；
图13示出了两个晶格单元；
图14-16是图1的流程图的步骤的另一图示说明；
图17-19示出了边界表面的示例。
具体实施方式
参照图1的流程图，提出一种用于设计包含晶格结构的三维(3D)部件的3D有限元网格的计算机实现的方法。该方法包括提供3D部件的实体表示以及提供单元的规则拼接。该规则拼接的每个单元包括有限元网格。该方法还包括将单元的规则拼接与3D部件的实体表示相叠加。然后，将所述单元分成两个组，即第一组和第二组。第一组单元包括与3D部件的实体表示接触的单元，并且第二组单元包括与3D部件的实体表示不接触的单元。接下来，计算第一组单元和3D部件的实体表示的布尔并集。该布尔并集形成体积。然后，该方法包括对所计算的布尔并集的体积进行有限元网格划分，同时保留第一组单元的与第二组单元共享的面的集合。然后，将第二组单元中的单元的有限元网格与所计算的布尔并集的经网格划分的体积进行合并。
这种方法改善了3D有限元网格的创建，并且允许用有限元方法来分析至少部分地由晶格组成的实体的结构性能。本发明通过首先创建与每个晶格单元相关联的有限元网格来实现该目的。当将各个晶格单元放置在一起时，这些晶格单元的相关联的有限元被添加到网格中，使得生成整个网格的过程非常有效。所有有限元在所述单元的公共面上节点对节点地完全匹配。然后，采用以适应性方式创建网格的填充物网格划分，将晶格区域的由晶格生成的有限元附接到一个(或多个)非晶格区域。这是可行的，因为第一组单元的与第二组单元共享的面的集合被保留。因此，本发明允许将由晶格生成的有限元附接到非晶格区域中使用的“填充物(filler)”元素上。生成了一个连贯的整体网格。
因此，本发明提出了一种混合网格划分方法，该方法针对全部晶格区域使用预先计算的有限元网格的重复，并且使用填充物网格划分来填充非晶格区域。填充物网格划分相当于执行以自适应方式创建网格。本发明使预先网格划分的晶格与填充物网格划分的非晶格区域正确匹配。有趣的是，本发明允许处理重叠的单元。最终结果是一种用于对晶格实体进行网格划分的方法，其性能、可靠性和精度高于现有技术能够达到的性能、可靠性和精度。此外，由于本发明依赖于每个晶格单元中的有限元的预先计算的布局，因此有可能对该单元使用六面体占大多数的布局。对于许多应用来说，六面体占大多数的网格优于现有技术中的自动化3D网格划分所生成的四面体占大多数的网格。
该方法是计算机实现的。这意味着该方法的步骤(或基本上所有步骤)由至少一个计算机或任何类似***执行。因此，该方法的步骤由计算机执行，可能完全自动地或半自动地执行。在示例中，可以通过用户-计算机交互来执行对该方法的至少一些步骤的触发。所需的用户-计算机交互的级别可能取决于所预见的自动化级别，并与实现用户愿望的需要保持平衡。在示例中，该级别可以是用户定义的和/或预定义的。
例如，可以在用户操作时触发用于提供的步骤S10、S20。例如，在图1的步骤S10，用户选择实体表示，诸如3D部件，并且作为该选择的结果执行步骤S10。
方法的计算机实现的典型示例是使用适用于此目的的***执行该方法。该***包括耦合到存储器的处理器。可选地，该***可以包括：用于显示图形用户界面(GUI)的显示器，存储器具有在其上记录的计算机程序，计算机程序包括用于执行该方法的指令。存储器还可以存储数据库。存储器是适用于这种存储的任何硬件，可能包括几个物理上不同的部分(例如，一个用于程序，并可能一个用于数据库)。
所谓“数据库”，是指为搜索和检索而组织的任何数据(即，信息)集合(例如，关系数据库，例如，基于预定结构化语言，例如，SQL)。当存储在存储器中时，数据库允许计算机进行快速搜索和检索。数据库的结构确实是为了便于与各种数据处理操作相结合地存储、检索、修改和删除数据。数据库可以由文件或文件集合组成，文件或文件集合可以被分解为记录，每个记录由一个或多个字段组成。字段是数据存储的基本单元。用户可以主要通过查询来检索数据。通过使用关键字和排序命令，用户可以根据所使用的数据库管理***的规则快速地搜索、重新排列、分组和选择许多记录中的字段，以检索或创建关于特定数据集合的报告。
该方法通常处理建模对象。建模对象是由存储在例如数据库中的数据定义的任何对象。通过扩展，表述“建模对象(modeled object)”表示该数据本身。根据***的类型，建模对象可以由不同类型的数据定义。该***实际上可以是CAD***、CAE***、CAM***、PDM***和/或PLM***的任何组合。在这些不同的***中，建模对象由相应的数据定义。本领域技术人员可以相应地谈论CAD对象、PLM对象、PDM对象、CAE对象、CAM对象、CAD数据、PLM数据、PDM数据、CAM数据、CAE数据。然而，这些***并非互斥的，因为建模对象可以由对应于这些***的任何组合的数据定义。因此，***可以是CAD和PLM***，从下面提供的这些***的定义中可以明显看出。
“CAD***”还表示任何至少适用于基于建模对象的图形表示来设计该建模对象的***，诸如CATIA。在这种情况下，定义建模对象的数据包括允许建模对象的表示的数据。在某些具有面(face)或表面(surface)的情况下，CAD***可以例如使用边或线来提供CAD建模对象的表示。线、边或表面可以用各种方式来表示，例如非均匀有理B样条(NURBS)。具体而言，CAD文件包含规格(specification)，可以从规格中生成几何图形，这又允许生成表示。建模对象的规格可以存储在单个CAD文件中，也可以存储在多个CAD文件中。在CAD***中表示建模对象的文件的典型大小是在每部件1兆字节的范围内。并且，建模对象通常可以是数千个部件的组装件。
在CAD的上下文中，建模对象通常可以是3D建模对象，例如，表示产品，诸如部件或部件的组装件，或者可能表示多个产品的组装件。所谓“3D建模对象”，是指任何由允许其3D表示的数据建模的对象。3D表示允许从各个角度观察部件。例如，当对3D建模对象进行3D表示时，可以围绕3D建模对象的任何轴对其进行处理和旋转，或者围绕显示该表示的屏幕中的任何轴旋转。这尤其排除了2D图标，它们不是3D建模的。3D表示的显示便于设计(即，提高设计人员在统计上完成其任务的速度)。这加快了工业中的制造过程，因为产品的设计是制造过程的一部分。
3D建模对象可以表示在完成其虚拟设计之后在真实世界中要制造的产品的几何形状，该虚拟设计使用例如CAD软件解决方案或CAD***，诸如(例如，机械)部件或部件的组装件(或等同于部件的组装件，因为从该方法的角度来看，部件的组装件可以被视为一个部件本身，或者该方法可以独立地应用于组装件的每个部件)，或者更一般而言，应用于任何刚性体组装件(例如，移动机构)。CAD软件解决方案允许在各种的并且无限制的工业领域中设计产品，包括：航空航天、建筑、建造、消费品、高科技设备、工业设备、运输、海事和/或海上油/气生产或运输。由此，由该方法设计的3D建模对象可以表示工业产品，其可以是任何机械部件，诸如地面车辆(包括例如汽车和轻型卡车设备、赛车、摩托车、卡车和摩托设备、卡车和公共汽车、火车)的部件、飞行器(包括例如机身设备、航空航天设备、推进设备、国防产品、航空设备、空间设备)的部件、海洋船只(包括例如海军设备、商船、离岸设备、游艇和工作艇、海洋设备)的部件、一般机械部件(包括例如工业制造机械、重型移动机械或设备、安装的设备、工业设备产品、制造的金属产品、轮胎制造产品)、机电或电子部件(包括例如消费电子、安全和/或控制和/或仪表产品、计算和通信设备、半导体、医疗设备和装备)、消费商品(包括例如家具、家庭和花园产品、休闲用品、时尚产品、耐用品零售商的产品、非耐用品零售商的产品)、包装(包括例如食品和饮料以及香烟、美容和个人护理、家用产品包装)。
PLM***还意味着任何适用于管理表示物理制造的产品(或将要制造的产品)的建模对象的***。因此，在PLM***中，建模对象由适于制造物理对象的数据定义。这些数据通常可以是尺寸值和/或公差值。对于一个对象的正确制造，拥有这样的值确实更好。
CAM解决方案还意味着任何适用于管理产品的制造数据的解决方案、硬件软件。制造数据通常包括与要制造的产品、制造过程和所需资源有关的数据。CAM解决方案用于规划和优化产品的整个制造过程。例如，其可以向CAM用户提供关于制造过程的可行性、持续时间或可在制造过程的特定步骤中使用的资源数量(诸如特定机器人)的信息；并且从而允许就管理或所需投资作出决定。CAM是在CAD过程和潜在的CAE过程之后的后续过程。这种CAM解决方案是由达索***使用商标来提供的。
CAE解决方案还意味着适合于分析建模对象的物理行为的任何解决方案、硬件软件。一种众所周知且广泛使用的CAE技术是有限元方法(FEM)，该方法通常涉及将建模对象划分成可以通过方程计算和仿真其物理行为的元素。这种CAE解决方案由达索***使用商标来提供。另一种正在发展的CAE技术涉及对由来自不同物理领域的多个组件组成的复杂***的建模和分析，而没有CAD几何数据。CAE解决方案允许对要制造的产品进行仿真，并且从而允许进行优化、改善和验证。这种CAE解决方案是由达索***公司使用商标提供的。
PDM代表产品数据管理。PDM解决方案意味着适合于管理与特定产品相关的所有类型的数据的任何解决方案、硬件软件。PDM解决方案可以被在产品生命周期中所有参与者使用：主要是工程师，但也包括项目经理、财务人员、销售人员和买家。PDM解决方案通常基于面向产品的数据库。其允许参与者共享关于其产品的一致数据，并且从而防止参与者使用有差异的数据。这种PDM解决方案是由达索***使用商标提供的。
图2示出了***的GUI的示例，其中该***是CAD***。
GUI 2100可以是典型的类CAD接口，具有标准菜单栏2110、2120以及底部和侧部工具栏2140、2150。这些菜单栏和工具栏包含一组用户可选择的图标，每个图标与一个或多个操作或功能相关联，如本领域已知的。这些图标中的一些图标与软件工具相关联，所述软件工具适用于对在GUI2100中显示的3D建模对象2000进行编辑和/或处理。软件工具可以分组到工作台(workbench)中。每个工作台包括软件工具的子集。特别地，其中一个工作台是编辑工作台，适用于对被建模产品2000的几何特征进行编辑。在操作中，设计人员可以例如预先选择对象2000的一部分，然后通过选择适当的图标来发起操作(例如，改变尺寸、颜色等)或编辑几何约束。例如，典型的CAD操作是在屏幕上显示的3D建模对象的冲压或折叠的建模。GUI可以例如显示与所显示的产品2000相关的数据2500。在该图的示例中，被显示为“特征树”的数据2500及其3D表示2000涉及包括制动卡钳和制动盘的制动组装件。该GUI可以进一步显示各种类型的图形工具2130、2070、2080，例如用于促进对象的3D定向、用于触发对所编辑的产品的操作的仿真或呈现所显示的产品2000的各种属性。光标2060可以由触觉设备控制以允许用户与图形工具交互。
图3示出了***的示例，其中该***是客户端计算机***，例如用户的工作站。
该示例的客户端计算机包括连接到内部通信总线1000的中央处理单元(CPU)1010、同样连接到总线的随机存取存储器(RAM)1070。客户端计算机还设有图形处理单元(GPU)1110，图形处理单元(GPU)1110与连接到总线的视频随机存取存储器1100相关联。在本领域中，视频RAM1100也被称为帧缓冲器。大容量存储设备控制器1020管理对大容量存储设备(诸如硬盘驱动器1030)的访问。适用于有形地体现计算机程序指令和数据的大容量存储设备包括所有形式的非易失性存储器，包括：例如，半导体存储设备，诸如EPROM、EEPROM和闪存设备；磁盘，诸如内部硬盘和可移动磁盘；磁光盘；以及CD-ROM盘1040。上述任一项可以附带有特别设计的ASIC(专用集成电路)，或者合并到特别设计的ASIC(专用集成电路)中。网络适配器1050管理对网络1060的访问。客户端计算机还可以包括触觉设备1090，诸如光标控制设备、键盘等。在客户端计算机中使用光标控制设备，以允许用户选择性地将光标定位在显示器1080上的任何期望位置处。此外，光标控制设备允许用户选择各种命令，并输入控制信号。光标控制设备包括用于向***输入控制信号的许多信号产生设备。通常，光标控制设备可以是鼠标，鼠标的按钮用于生成信号。可选地或附加地，客户端计算机***可以包括触摸板和/或触摸屏。
计算机程序可以包括可由计算机执行的指令，指令包括用于使上述***执行该方法的单元。程序可以在任何数据存储介质上记录，包括***的存储器。程序可以例如以数字电子电路实现，或者以计算机硬件、固件、软件实现，或者以它们的组合实现。程序可被实现为装置，例如有形地体现在机器可读存储设备中的产品，用于由可编程处理器执行。方法步骤可由可编程处理器执行，可编程处理器执行指令的程序以通过对输入数据进行操作并生成输出来执行方法的功能。处理器因此可以是可编程的并且被耦合以从数据存储***、至少一个输入设备和至少一个输出设备接收数据和指令以及向其发送数据和指令。应用程序可以用高级过程性或面向对象的编程语言实现，或者如果需要，可以用汇编语言或机器语言实现。在任何情况下，语言可以是被编译或解释的语言。程序可以是完整的安装程序或更新程序。在***上应用该程序在任何情况下都会产生执行该方法的指令。
“设计3D建模对象”指的是作为详细描述3D建模对象的过程的至少一部分的任何动作或动作系列。因此，该方法可以包括从零开始创建3D建模对象。可替换地，该方法可以包括提供先前创建的3D建模对象，然后修改该3D建模对象。
该方法可以包括在制造过程中，所述制造过程可以包括：在执行该方法之后，生产与建模对象相对应的物理产品。在任何情况下，由该方法设计的建模对象可以表示制造对象。因此，建模对象可以是被建模的实体(即，表示实体的建模对象)。制造对象可以是产品，诸如部件，或者部件的组装件。由于该方法改善了建模对象的设计，所以该方法也改善了产品的制造，并且从而提高了制造过程的生产率。
尤其是，包括该方法的制造过程可以包括通过FEM仿真物理现象。例如，所述3D表示对飞机机翼进行建模，并且FEM被用来仿真机翼在飞行过程中的行为。这种FEM仿真允许预测使用中的物理***的性能(热性能、机械性能或任何其它物理性能)，以及优化设计以提高所述性能。制造过程可以包括预测部件的可制造性，并最终制造部件。在3D打印的情况下，可制造性可能是成败的关键，例如，对具有诸如在直接金属打印期间部件中的热分布之类的现象的打印过程进行仿真。事实上，诸如残余应力之类的现象在3D打印中可能存在严重问题，并且可能使得物理***没有资格被用于诸如航空之类的关键应用。“3D部件的3D打印”也称为“3D部件的增材制造”或“3D建模对象的增材制造”。增材制造(additivemanufacturing，AM)指的是通过将能量和/或材料顺序输送到空间中的指定点以生产实体3D部件来制造该3D部件的过程。所述材料可以是但不限于：逐层添加的材料，不论该材料是塑料、金属、混凝土…。生产与建模对象(所提供的3D部件)相对应的物理产品。增材制造的对象可以是产品，诸如部件、或者部件的组装件。该方法改善了所提供的3D部件的3D表示的可制造性，并且因此该方法也改善了产品的制造，并且从而提高了制造过程的生产率。
现在参考图1的流程图，讨论3D部件的3D有限元网格的设计方法的示例。该3D部件包括晶格结构。
该方法可以从提供(S10)3D部件的实体表示开始。3D部件通常对物理***进行建模，例如要制造的真实世界的产品。如前所述，3D部件可以是3D建模对象。该实体表示包括至少一个要用晶格填充的空腔。空腔是可以是封闭或不封闭的体积。空腔是在所提供的3D部件中的中空；因此空腔没有任何材料。更一般地说，空腔可以从实体表示的边界的任何连接的子集创建。此后，连接的子集将被称为空腔表面。在空腔表面形成一个封闭的体积的情况中，这就明确地定义了要由晶格填充的空间。如果不是，即如果空腔不形成封闭的体积，则空腔表面可以包括一个开口，例如它类似于碗；空腔甚至可以具有两个或更多个开口，例如具有两个开口的管道。在实践中，为了从空腔表面构建封闭的体积，设计人员将指定一个或多个盖表面(cap surface)。盖表面的存在可以堵住一个开口，就像罐子的盖子。盖表面可以通过合并将会位于开口位置的内边界和外边界来获得。一个或多个盖表面的集合与空腔表面一起限定了一个封闭的体积。
加盖(capping)过程开始于定义一个或多个边缘(rim)，这些边缘是在实体表示的边界上形成一个闭环的一维(1D)曲线或线的集合。该一个或多个边缘类似于罐盖附着的罐子的边缘。然后，将盖表面创建为将该边缘作为外边界的表面。这种操作可以通过“填充表面(fill surface)”来实现，“填充表面”是现代CAD/CAM***中常见的功能。对开放空腔加盖的过程是自愿的；设计人员可以在本发明的意义上，自由地决定将形状中的中空或管道(如果有的话)或封闭体积(如果有的话)中的哪一者作为用于由晶格填充的空腔来对待。例如，这在图18和19中表示。图18示出了图19的加盖实体190。加盖实体190包括外边界196、由边缘192限定的开口，该开口由盖表面194填充。如图18所示，加盖实体190包括空腔180，空腔180可以由晶格填充。因此，图18示出了包括一个开口的空腔的3D部件的实体表示，其中空腔已用盖表面处理并填充有晶格。图17还示出了包括已填充有晶格的一个封闭空腔的3D部件的实体表示。
晶格是多孔的—也可以说是肺泡结构，其使得降低材料的密度成为可能。术语“晶格”表示3D实体部件中交替出现的实体和空心的图案；换句话说，实体和空心的图案可能是重复的。晶格的目的是改善部件的机械阻力、稳定性、传热或其他性能。在非限制性示例中，由晶格结构引入的材料的孔隙率等于或大于50％：由材料定义的体积的50％或更多是空的(即，没有材料)。应当理解，本发明可以用具有任意百分比孔隙率的任何类型的晶格结构来实现。
仍然参照步骤S10，可以显示或不显示所提供的3D部件的实体表示。当显示时，实体表示显示在3D场景中，如本领域已知的。
三维部件的实体表示可以是边界表示(称为B-REP)，如本领域已知的。B-REP使用拓扑(顶点、边、表面)和几何(点、曲线、表面)信息来计算所述表示。
3D部件的实体表示可以是几何体，其是点、段、多边形面或更复杂的元素(例如，基于NURBS的表面)的集合。更一般地说，曲面表示是网格。不需要关于这些元素的方向或连通性的信息。例如，如果用NURBS对3D部件建模，则所述表示将被细分(例如，所述表示被细分为三角形)。
如图18所示，实体表示190包括外边界196和一个或多个空腔表面182，其形成了实体表示的表皮或包络。一个或多个空腔表面182与盖表面194(如果有的话)组合限定了实体表示的一个或多个空腔180；这等于说，空腔表面是实体表示的表皮的与空腔接触的部分。外边界196是实体表示的表皮的与空腔不接触的部分。在只具有封闭空腔的实体的情况中，如图17所示，空腔表面不与外边界连接。在具有一个或多个开口空腔的实体的情况下，如图18所示，外边界连接到该开口空腔的空腔表面。空腔表面与外边界之间的边界由一个或多个边缘192形成。
提供实体表示意味着用于显示实体表示或至少用于计算实体表示所需的信息可用于执行该方法的计算机化***。例如，这些信息存储在计算机化***的存储器中。实体表示是在用户操作时提供的，例如，用户在本领域公知的图形用户界面上选择文件。可以自动提供实体表示；在这种情况下，不需要用户的干预。
图4的左边部分示出了在步骤S10中提供的实体表示的示例。在此示例中，3D部件40的表示包括空腔42，并且该空腔由大致具有矩形形状的材料限定。材料具有厚度48。因此，所提供的3D部件的实体表示具有空腔表面46和外边界44，所述空腔表面46和外边界44形成充满材料48的实体表示的部分。在这个例子中，没有盖表面，并且因此空腔表面完全限定了空腔42：空腔表面包括面的集合，并且该面的集合中的每个面都有助于限定空腔。
在图1的步骤S20，提供单元的规则拼接。单元的拼接使体积离散化，并且每个单元表示离散化体积的子体积；因此，该离散化体积内的任何点位于(或封闭于)一个单元中。因此，单元即是3D对象。
拼接是规则拼接，也就是说，存在一个由诸如(等边)三角形、正方形、六边形等等的全等的规则的多边形构成的边对边的拼接。实际上，拼接包括是四面体或六面体的单元，并且全等的四面体或六面体共享至少一个匹配在一起的公共面(面重合)：公共面具有相同的形状和大小，并且公共面在拼接中彼此重叠。例如，立方体的规则拼接形成体素的规则网格。
每个单元包括一个有限元网格：在方法的这个步骤中，有限元网格可以被实例化或不被实例化。如果网格没有被实例化，那么其将在稍后的阶段(S52)。FEM是基于从问题域的精化或“离散化”中获得的数值模型，并结合了物理问题的完整定义所需的附加信息。有限元模型是有限元网格与仿真属性的集合的关联。如本领域所公知，执行每个单元的网格生成过程，可以使用生成网格的任何过程而不会脱离本发明，包括特定布局的手动交互创建。在佛罗里达州博卡拉顿的CRC出版社2013年的S.W.Cheng，T.Dey和J.Shewchuk著作的“Delaunay Mesh Generation”精装本(ISBN978-1-58488-730-0)中，给出了可用于本发明的自动算法的示例，要理解的是，可以使用任何其他网格生成过程。在示例中，为了在每个单元内生成网格，设计人员可能希望手动地对元素进行排列。虽然这个过程缓慢，但是将给出最佳类型的排列，并且由于相同的单元可能被本发明重复使用数千次甚至数百万次，因此额外的努力可能是值得的。因此，每个单元包括更小的有限元。每个单元的FEM网格通常是四面体或六面体，但其最好是六面体网格，因为FEM分析在六面体占优势的网格化单元上表现最好。在本发明中，包括有限元网格意味着单元被分解成有限元，或者单元可分解成有限元并且在根据本发明的方法的后期阶段将被分解成有限元网格。
规则拼接的单元在拓扑上可以是相同的。这使得对单元的重新使用成为可能；最终可能需要一个或多个单元的变形。单元可以是拓扑和几何上相同的，并且两个单元的迭合面(congruent faces)必须彼此完全匹配；通过这种方式，拼接保持规则。尤其是，如果遵守正确的间距，则在拼接的单元之间没有可能的相交。
拼接的单元是晶格单元。晶格单元对表示晶格结构的图案进行建模；换言之，晶格单元是包括多孔材料的单元。因此，晶格单元不充满材料，而是部分地没有材料。例如，晶格单元可以是正六面体(立方体)，其仅包括沿着该立方体的对角线的材料，例如连接该立方体的相对边的材料条。作为另一个示例，晶格单元可以是立方体，该立方体仅包括沿着连接该立方体的两个相对面的中心的轴线的材料。
选择晶格单元的图案以使两个单元的迭合面的交界面处的材料匹配。换句话说，该图案按照相邻面的对进行重复；因此相邻面是迭合面。例如，晶格立方体最多有三个面图案：针对该立方体的每对相对面，有一个面图案。
规则拼接的所有单元可以呈现相同的图案，这涉及所有单元具有相同的有限元网格；因此，单元的规则拼接中的单元的网格生成更快，因为为一个单元生成的网格将应用于其他单元。
图4的右边部分显示了形成拼接的单元集合的示例。为了便于理解该图，该图用2D表示拼接；本例中的拼接的单元(标记为a到l)是具有相同拓扑结构的立方体。每个立方体都是具有基本晶格结构的晶格单元，单元中心附近的体积是没有任何材料的。例如，这在图5中图示，其示出了图4中标记为a到l的单元之一的示例。在图5中，“空(void)”表示没有材料，“满(full)”表示充满材料。
图6示出了图5的单元，其中已经生成或计算了有限元网格。这个2D示例中的FEM网格是三角形的，并且类似于具有四面体网格的3D例子。对于图5的“空”体积没有计算FEM网格。
图13显示了图1中标记为a和b的晶格单元的细节。单元a和单元b具有相同的拓扑结构和相同的晶格图案。单元a和b各具有匹配的迭合面，分别是a1和b1。也就是说，顶点1与顶点5、顶点2与顶点6、顶点4与顶点7、顶点3与顶点8分别相符。这里，晶格图案是凹槽(recess)，其包括三个通孔，每个通孔都居中于其穿过的两个相对面上。例如，单元a的面a1上的孔a2与单元b的面b1上的孔2b相同，使得当单元a和b属于拼接时，面a1和1b的FEM网格的三角形的顶点和边相符。
返回参考图1，在步骤S22，将单元的规则拼接与3D部件的实体表示相叠加。规则拼接的尺寸可以等于或大于要由在步骤S10处提供的实体表示的晶格填充的至少一个空腔。将规则拼接与实体表示相叠加意味着规则拼接与实体表示相符；更具体地说，规则拼接与实体表示的至少一个空腔相符。作为步骤S22的结果，规则拼接与实体表示的至少一个空腔相交；用规则拼接填充空腔。如果拼接比空腔大，那么规则拼接的一部分将与实体表示相交(或冲突)。
正如所讨论的，规则拼接可能大于空腔。但是，规则拼接甚至可以大于加盖的实体表示。由于这个原因，作为步骤S22的叠加的结果，拼接的单元可能位于加盖的实体表示的外部。这些在加盖的实体表示的外部的单元与加盖的实体表示的至少一个空腔的外部一样是无用的。在步骤S30，删除规则拼接中的完全位于加盖的实体表示外部的单元。因此，保持规则拼接中的、与加盖的实体表示接触(例如通过冲突)和/或包含在加盖的实体表示的至少一个空腔中的单元。
在步骤S30中未删除的单元被划分成两组单元(步骤S40)。分成两组就意味着这些单元被分为第一组单元和第二组单元。
第一组单元被创建(S42)，并且包括与3D部件的实体表示接触的单元。“接触”意味着：(i)单元和实体表示之间存在冲突；单元和实体表示之间至少有一个公共点。
第二组单元被创建(S44)，并且包括与3D部件的实体表示没有接触的单元。应当理解，第二组可以是空的，尽管这种情况并不完全利用本发明的优点。
步骤S42和S44可以同时执行，也可以一个接一个地执行。
图7示出了单元的规则拼接与图4所示的实体表示的叠加。在这个示例中，标注为a、d、i和l的单元与3D部件的实体表示相接触70—它们冲突。因此，这四个单元被添加到第一组单元中。拼接的其他单元被添加到第二组单元中，这些单元与实体表示没有接触(或没有冲突)。图8示出作为步骤S44的结果而获得的没有冲突的晶格单元的集合的表示。
当执行步骤S40(包括子步骤S42、S44)时，或者在已经执行步骤S40之后，计算第一组单元中的单元和第二组单元中的单元之间的交界面，并最终进行存储。交界面可以例如存储在执行本发明的方法的***的存储器中。
所述交界面包括第一组单元和第二组单元共有的、有限元网格的面的集合。应该理解，第一组单元和第二组单元共有的、有限元网格的这些面依赖于第一组单元和第二组单元中的单元的匹配的迭合面。因此，交界面包括作为有限元网格划分的结果而获得的二维(2D)表面。因此，这些2D表面是有限元的面，典型地是三角形(tri)或四边形(quad)，并且根据其排列，其可以形成一个或多个表面。我们将这些表面称为2D有限元表面，或2D有限元网格。
交界面的边界可以包括2D有限元网格的自由边。自由边是表面的与孔接触的边。在这种情况下，交界面的每一组被连接的表面将产生一条闭合折线(poly-line)。这个闭合折线(也可以称为环)限定了交界面的外边界。对于所产生的每个外边界，定义描述交界面的内边界的零条或更多条闭合折线。
回到图13，标注为a和b的两个单元是拓扑上相同的晶格单元，并且因此其具有相同的FEM网格。因此，两个面a1、b1的三角形相符并形成交界面。
在图1的步骤S50，去除第一组单元中的单元的内部网格结构。换言之，第一组单元中的单元的体积有限元被忘记。忘记第一组单元的内部网格结构涉及简单地删除元素。如果在提供步骤S20处没有实例化单元的有限元网格，则不执行步骤S50。
在步骤S52，实例化第二组单元中的每个单元的有限元网格。在具有拼接的步骤S20中提供的单元的有限元网格没有被实例化(即没有被计算)的情况下，执行此步骤。按照本领域所公知的，来执行实例化。
图10例示了第二组单元中的单元及其实例化的有限元网格。
通常刚好在执行了步骤S42(S44)之后执行步骤S50(相应的S52)。可替换地，可以同时执行步骤S50、S52。
然后，在步骤S60，计算第一组单元和3D部件的实体表示的布尔并集。布尔并集运算是按照本领域公知的来执行的，并且可以使用任何算法来计算部分之间的布尔并集。作为布尔并集的结果，得到一个新的部分，并且这个新的部分形成一个体积。这等于说，布尔并集形成了一个体积。尤其是，这个体积仍然包括原始提供的表示的空腔，除非第二组单元是空的—也就是说，没有单元属于第二组，并且拼接的所有单元都属于第一组单元—这涉及所提供的拼接的单元的大小与要用晶格结构填充的空腔的体积相比太粗糙—拼接的每个单元总是与实体表示至少有一个接触。如果出现这种情况，根据本发明的方法可以用包括更小单元(体积更小的单元)的新拼接重新启动。
第一组单元和3D部件的实体表示的经修剪的布尔并集包括至少两个边界。第一个边界是经修剪的并集的外边界，其通常由一个或多个两流形(two-manifold)外壳组成，每个外壳具有零个或更多个两流形内壳。第二个边界是由在步骤S46计算并存储的交界面最初形成的边界；该交界面边界包括一个或多个单流形外环，每个外环具有零个或更多个单流形内环。
一旦已经执行步骤S60，就可以在该布尔并集内区分出仅来自第二组单元中的单元的材料、仅来自3D部件的实体表示的材料以及来自实体表示和第一组单元中的单元的并集的材料。对于来自第一组单元中的单元的材料，可以区分出位于3D部件外边界的实体表示内部的单元部分和3D外边界的实体表示外部的单元部分。位于3D外边界的实体表示外部的单元部分是不希望的，因为它们造成了所提供的实体表示的原始形状的变形。
为了保持实体表示的原始形状—即保持实体表示的外边界—可以执行对所计算的第一组单元和3D部件的实体表示的布尔并集进行修剪(S70)的步骤。经修剪的并集保持与加盖的实体表示相同的外边界。换言之，所提供的实体表示的原始外边界未被修改，并且只有实体表示的空腔表面(即，与所述表示的与空腔接触的表皮或包络的子集)可作为步骤S60的并集的结果被修改。
图9示出了图4的实体表示(左部分)与第一组单元a、d、i和l的布尔并集。在这个示例中，只要第一组单元中的单元不会影响实体的外包络，就不需要步骤S70的修剪。因此，修剪步骤是可选的，并且可以在已经检测到或识别出对实体表示的原始外边界的修改的情况下执行。
接下来，在步骤S80，以完成后的表面网格包含在步骤S46处所存储的交界面的方式，计算经修剪的并集的边界的有限元表面网格。所得到的表面网格是用于计算所计算的布尔并集的体积的有限元网格的输入。因此，精确地形成交界面的2D面由此被保存在所提供的实体表示的空腔表面上，该空腔表面已经通过与第一组单元的布尔并集而被改变。换句话说，保留了第一组单元的与第二组单元共享的面的集合。可以使用任何有限元网格生成器；特别地，对于表面网格划分，可以使用任何三角形或四边形网格生成器，并且对于体积网格划分，可以使用任何四面体或六面体网格生成器。
注意，在封闭空腔(图17)的情况下，空腔表面的网格划分和布尔并集的外边界的网格划分将是分开的操作，因为两者(空腔表面和外边界)没有连接。相比之下，对于图18中所描述的开放空腔，空腔表面的网格划分也需要在相同操作中的外边界的网格划分，因为二者是连接的。在这两种情况下，必须首先对所有的表面、空腔和外边界进行网格划分，以便允许在布尔并集上创建3D有限元。
图11示出了经过有限元网格划分的、图9的布尔并集的体积的示例。虚线100表示已经保留的步骤S46的交界面。有趣的是，该交界面100同样呈现在图10上，因为图11和10上的两个交界面是相同的：布尔并集(经修剪或未修改)的内边界与第二组单元中的单元集合形成的外边界相同。
在布尔体积的计算中，可能出现对布尔并集的空腔表面的修改。尤其是，这些修改可以涉及将一个单元与实体表示的空腔表面连接的面(例如，三角形)：该单元的一个或多个面可以被修改，以创建与在S10处提供的体积的空腔表面的一个或多个面的连接。换言之，体积的交界面边界可以与步骤S46处存储的交界面稍有不同。结果，布尔并集的内边界的有限元表面网格可能不遵循步骤S46处所存储的交界面，例如，有限元表面网格的点可能位于修改后的三角形上。
为了确保布尔并集(S60、S70)的空腔表面的有限元表面网格遵循第一组单元和第二组单元的面的集合之间的交界面，可以在执行布尔并集的体积的有限元网格划分之前，将所存储的交界面合并到所计算的布尔并集的有限元表面网格上。所存储的交界面的合并可能需要执行下面定义的两个策略中的至少一个策略。这些策略教导了在计算(经修剪的)布尔并集的有限元网格时，如何能够遵循原始交界面(S46)。
在示例中，这两种策略可用于计算所计算的布尔并集的体积的有限元网格。这两种策略都依赖于使用高级前沿网格生成器，该生成器将被应用于所计算的布尔并集的经网格划分的体积。可以使用任何高级前沿网格生成器。这两种策略的主要区别在于，对所计算的布尔并集的经网格划分的体积的约束将如何修改前沿的传播。
第一种策略使用所合并的交界面的一个或多个自由边作为应用于布尔体积的边界的高级前沿网格生成器的初始前沿。换言之，经修剪的体积的边界的重建是从描述交界面的1D边界的闭合折线(如果有的话)开始的。一旦已经重建2D表面边界，就可以用新计算的有限元表面网格进行布尔并集的有限元体积网格划分。因此，该第一种策略可以包括以下步骤。首先，将在步骤S46处存储的交界面边界合并到经修剪的并集的边界上；例如，将交界面叠加在经修剪的并集的内边界上，因为在构造上交界面相符。接下来，对经修剪的布尔并集的边界应用前沿生成器，其中描述内边界的折线(也称为环)用作前沿生成器的初始前沿。其结果是将这些环合并为嵌入的有限元的边的有限元表面网格。换句话说，经修剪的布尔并集的孔(自由边)已经被填充。然后，由于这些描述内边界的环的嵌入，将被计算的有限元网格(S80)将包含其边界与交界面边界(S46)完全相符的连接的有限元子集。这个有限元子集可以被删除。最后，将交界面的2D有限元与前沿网格生成器生成的其余有限元进行组装(或连接)。其结果将是准确地合并了交界面的2D有限元的一个2D有限元表面网格。
第二种策略是类似的，但并非是创建覆盖边界的网格并随后删除一些元素，而是选择性地创建网格，以使得经修剪的体积的边界的每个孔(由自由边的集合定义)不会被填充(例如不会由前沿网格生成器计算的三角形填充)。该第二种策略可以包括以下步骤。首先，将在步骤S46处存储的交界面边界合并到经修剪的并集的边界上；可以按照第一种策略的方式，来执行该合并。接下来，在经修剪的布尔并集的边界上应用前沿生成器，其中，将描述交界面的外边界的折线(也称为环)用作前沿生成器的初始前沿，反之亦然。只在该前沿的一侧应用前沿发生器，意思是(自动或经由用户输入)使前沿远离孔或环的内部。结果是其中经修剪的布尔并集的孔(自由边)没有被填充的一个有限元表面网格。然后，由于在S46处获得的交界面的嵌入，将被计算(S80)的有限元网格将包括其边界与交界面边界完全相符的连接的有限元子集。换言之，其结果将是准确地合并了交界面的2D有限元的一个2D有限元表面网格。
在步骤S90，对第二组单元中的单元的有限元网格和所计算的布尔并集的经网格划分的体积进行合并。结果，得到了一个有限元网格。该合并可以在不需要对两个网格进行任何进一步修改的情况下进行：实际上，所计算的布尔并集的经网格划分的体积的边界和第二组单元中的单元的外边界呈现相同的有限元表面网格，从而使得该合并相当于识别两个有限元表面网格上的公共点对。步骤S90的合并通常是按照本领域所公知地来执行。
虽然已经参照本发明的示例性实施例具体示出和描述了本发明，但是本领域技术人员将理解，在不脱离所附权利要求所涵盖的本发明范围的情况下，可以在形式和细节上作出各种改变。例如，图14至16示出了与参考附图4到13讨论的示例类似的本发明的示例。在此示例中，第一组单元与3D部件的实体表示的布尔并集包括位于实体表示的外边界之外的材料；在这种情况下，为了保持所提供的实体表示的外边界，需要修剪(S70)(即，不是可选的)。
本发明还可以包括对包含晶格结构的3D部件进行增材制造。该3D部件包括根据前面讨论的方法设计的三维(3D)有限元网格。然后将3D部件提供给3D打印机，3D打印机执行增材制造或3D打印。3D打印机可以与提供3D部件的***可通信地耦合。3D打印机通常包括用于存储3D部件的存储器。可以在增材制造之前对所提供的3D部件进行仿真；该仿真通常是用有限元方法仿真至少一个物理现象。所述至少一个物理现象可以是在增材制造期间发生的物理现象。例如，可能重要的是，仿真在金属直接打印过程中3D部件中的热分布或打印过程中的残余应力。因此，3D部件受益于用于预测在使用中的3D部件的性能和/或预测对3D部件的增材制造的可制造性的仿真。因此，所制造的3D部件提供了改善的物理和/或技术特性。

Claims (15)

1.一种用于设计包含晶格结构的三维(3D)部件的3D有限元网格的计算机实现的方法，该方法包括：
-提供(S10)所述3D部件的实体表示；
-提供(S20)单元的规则拼接，每个单元包括有限元网格；
-将所述单元的规则拼接与所述3D部件的实体表示相叠加(S22)；
-将所述单元划分(S40)成两组：
-第一组单元(S42)，其中每一个单元与所述3D部件的实体表示相接触，以及
-第二组单元(S44)，其中没有一个单元与所述3D部件的实体表示相接触；
-计算(S60)所述第一组单元与所述3D部件的实体表示的布尔并集，所述布尔并集形成体积；
-对所计算的布尔并集的所述体积进行有限元网格划分(S80)，同时保留所述第一组单元的与所述第二组单元共享的面的集合；以及
-对所述第二组单元中的单元的有限元网格和所计算的布尔并集的经网格划分的体积进行合并(S90)。
2.如权利要求1所述的计算机实现的方法，其中，所述规则拼接的单元在拓扑结构上相同。
3.如权利要求1或2所述的计算机实现的方法，还包括：
-对所计算的所述第一组单元与所述3D部件的实体表示的布尔并集进行修剪(S70)。
4.如权利要求1-3中任一项所述的计算机实现的方法，其中，保留所述第一组单元的与所述第二组单元共享的面的集合还包括：
-存储(S46)所述第一组单元和所述第二组单元之间的交界面，所述交界面包括所述第一组单元的与所述第二组单元共享的面的集合。
5.如权利要求4所述的计算机实现的方法，其中，对所计算的布尔并集的所述体积进行有限元网格划分的步骤还包括：
-将所存储的交界面合并(S82)到所计算的布尔并集上。
6.如权利要求1-5中任一项所述的计算机实现的方法，还包括：
-将高级前沿网格生成器应用于所计算的布尔并集的经网格划分的体积。
7.如权利要求5和6的组合所述的计算机实现的方法，其中，将所合并的交界面的一个或多个自由边用作被应用于所合并的交界面的所述高级前沿网格生成器的初始前沿。
8.如权利要求5和6的组合所述的计算机实现的方法，其中，所合并的交界面的一个或多个自由边是从被应用于所合并的交界面的所述高级前沿网格生成器保留的。
9.如权利要求1-8中任一项所述的计算机实现的方法，还包括：在所述划分之后，
-移除(S50)所述第一组单元中的单元的内部网格结构。
10.如权利要求1-9中任一项所述的计算机实现的方法，还包括：在所述划分之后，
-实例化(S52)所述第二组单元中的每个单元的所述有限元网格。
11.一种增材制造的方法，该方法包括：
-由计算机对根据如权利要求1-10中任一项所设计的包含晶格结构的三维(3D)部件的3D有限元网格进行仿真；
-对所述3D部件进行增材制造。
12.一种包括用于执行如权利要求1-10中的任一项所述的方法的指令的计算机程序。
13.一种其上记录有如权利要求12所述的计算机程序的计算机可读存储介质。
14.一种***，其包括耦合到存储器的处理器和图形用户接口，所述存储器上记录有如权利要求12所述的计算机程序。
15.一种通过如权利要求11所述的增材制造的方法可获得的所制作的3D部件。
CN201811580554.8A 2017-12-24 2018-12-24 包含晶格结构的3d部件的3d有限元网格的设计 Pending CN110059332A (zh)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
EP17306931.1A EP3503040B1 (en) 2017-12-24 2017-12-24 Design of a 3d finite element mesh of a 3d part that comprises a lattice structure
EP17306931.1 2017-12-24

Publications (1)

Publication Number Publication Date
CN110059332A true CN110059332A (zh) 2019-07-26

Family

ID=60957110

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811580554.8A Pending CN110059332A (zh) 2017-12-24 2018-12-24 包含晶格结构的3d部件的3d有限元网格的设计

Country Status (4)

Country Link
US (1) US10796039B2 (zh)
EP (1) EP3503040B1 (zh)
JP (1) JP7248420B2 (zh)
CN (1) CN110059332A (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112182798A (zh) * 2020-09-11 2021-01-05 厦门大学 一种嵌入板结构的连接构件的网格构建方法和结构
CN113239473A (zh) * 2021-07-12 2021-08-10 中国空气动力研究与发展中心计算空气动力研究所 用于复合材料性能预测的升力体标模设计方法及飞行器
CN113600832A (zh) * 2021-06-30 2021-11-05 西安航天发动机有限公司 一种用于激光选区熔化成形块状支撑的省粉结构设计方法

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11610366B2 (en) * 2020-05-21 2023-03-21 General Lattice, Inc. Systems and methods for generating non-approximated three-dimensional representations
CN111625905B (zh) * 2020-05-29 2023-08-04 中国航空工业集团公司西安飞机设计研究所 一种螺旋桨飞机舱内噪声高频仿真方法
US20230152779A1 (en) * 2020-06-17 2023-05-18 Hewlett-Packard Development Company, L.P. Perforations in a membrane for a lattice structure
WO2022093223A1 (en) * 2020-10-29 2022-05-05 Hewlett-Packard Development Company, L.P. Surface lattice structure formation
CN114248519A (zh) * 2021-12-24 2022-03-29 中国舰船研究设计中心 一种基于3d打印的层板结构隐身材料及其制备方法

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1978487A2 (en) * 2007-04-06 2008-10-08 Hitachi, Ltd. Numerical analysis mesh generating method and apparatus
CN101814103A (zh) * 2010-04-01 2010-08-25 西北工业大学 基于超单元的多组件布局建模与结构优化设计方法
CN102306396A (zh) * 2011-09-15 2012-01-04 山东大学 一种三维实体模型表面有限元网格自动生成方法
US20120232848A1 (en) * 2011-03-07 2012-09-13 Masato Naito Method for creating finite element model of rubber composite
WO2013063433A1 (en) * 2011-10-26 2013-05-02 Engine Simulation Partners Method and apparatus for modeling interactions of the fluid with system boundaries in fluid dynamic systems
EP2680232A2 (en) * 2012-06-25 2014-01-01 Fujitsu Limited Decoupled parallel meshing in computer aided design
EP2985708A1 (en) * 2014-08-15 2016-02-17 Dassault Systèmes Navigation through the hierarchical structure of a three-dimensional modeled assembly of parts
US20160096318A1 (en) * 2014-10-03 2016-04-07 Disney Enterprises, Inc. Three dimensional (3d) printer system and method for printing 3d objects with user-defined material parameters
CN105528494A (zh) * 2015-12-29 2016-04-27 上海大学 基于三维元胞自动机的轻量化模型生成及优化方法
CN105874510A (zh) * 2014-01-09 2016-08-17 西门子产品生命周期管理软件公司 用于增材制造的晶格结构的结构保留拓扑优化的方法
CN106650124A (zh) * 2016-12-28 2017-05-10 谢亿民工程科技（常州）有限公司 一种基于有限元结果的连续体表面轻量化方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US8259104B2 (en) * 2008-08-28 2012-09-04 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Domain decomposition by the advancing-partition method for parallel unstructured grid generation
GB201003065D0 (en) * 2010-02-23 2010-04-07 Simpleware Ltd Image processing method and method of three-dimensional printing incorporating the same
CN108351906B (zh) * 2015-11-25 2022-05-17 西门子工业软件有限公司 用于具有晶格结构的部件的建模的***和方法
US10503149B2 (en) * 2016-09-30 2019-12-10 Honeywell Federal Manufacturing & Technologies, Llc System, method, and computer program for creating united cellular lattice structure

Patent Citations (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP1978487A2 (en) * 2007-04-06 2008-10-08 Hitachi, Ltd. Numerical analysis mesh generating method and apparatus
CN101814103A (zh) * 2010-04-01 2010-08-25 西北工业大学 基于超单元的多组件布局建模与结构优化设计方法
US20120232848A1 (en) * 2011-03-07 2012-09-13 Masato Naito Method for creating finite element model of rubber composite
CN102306396A (zh) * 2011-09-15 2012-01-04 山东大学 一种三维实体模型表面有限元网格自动生成方法
WO2013063433A1 (en) * 2011-10-26 2013-05-02 Engine Simulation Partners Method and apparatus for modeling interactions of the fluid with system boundaries in fluid dynamic systems
EP2680232A2 (en) * 2012-06-25 2014-01-01 Fujitsu Limited Decoupled parallel meshing in computer aided design
CN105874510A (zh) * 2014-01-09 2016-08-17 西门子产品生命周期管理软件公司 用于增材制造的晶格结构的结构保留拓扑优化的方法
CN105900145A (zh) * 2014-01-09 2016-08-24 西门子产品生命周期管理软件公司 用于在计算机辅助设计模型中创建三维网格结构以进行增材制造的方法
EP2985708A1 (en) * 2014-08-15 2016-02-17 Dassault Systèmes Navigation through the hierarchical structure of a three-dimensional modeled assembly of parts
US20160096318A1 (en) * 2014-10-03 2016-04-07 Disney Enterprises, Inc. Three dimensional (3d) printer system and method for printing 3d objects with user-defined material parameters
CN105528494A (zh) * 2015-12-29 2016-04-27 上海大学 基于三维元胞自动机的轻量化模型生成及优化方法
CN106650124A (zh) * 2016-12-28 2017-05-10 谢亿民工程科技（常州）有限公司 一种基于有限元结果的连续体表面轻量化方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
杨明;陈涛;廖敦明;曹流;唐玉龙;: "基于尖锐边界等值面填充算法的型壳有限元网格自动划分", 铸造, no. 11, 10 November 2016 (2016-11-10), pages 41 - 44 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112182798A (zh) * 2020-09-11 2021-01-05 厦门大学 一种嵌入板结构的连接构件的网格构建方法和结构
CN112182798B (zh) * 2020-09-11 2022-06-21 厦门大学 一种嵌入板结构的连接构件的网格构建方法和结构
CN113600832A (zh) * 2021-06-30 2021-11-05 西安航天发动机有限公司 一种用于激光选区熔化成形块状支撑的省粉结构设计方法
CN113239473A (zh) * 2021-07-12 2021-08-10 中国空气动力研究与发展中心计算空气动力研究所 用于复合材料性能预测的升力体标模设计方法及飞行器

Also Published As

Publication number Publication date
US20190197205A1 (en) 2019-06-27
EP3503040A1 (en) 2019-06-26
JP2019114260A (ja) 2019-07-11
JP7248420B2 (ja) 2023-03-29
US10796039B2 (en) 2020-10-06
EP3503040B1 (en) 2024-06-05

Similar Documents

Publication Publication Date Title
CN110059332A (zh) 包含晶格结构的3d部件的3d有限元网格的设计
JP7235478B2 (ja) ３ｄ部品の付加製造
US11256832B2 (en) Replica selection
CN110008495A (zh) 包含晶格结构的3d部件的3d有限元网格的设计
CN105761303B (zh) 在3d建模组件上创建包围盒
CN105760570B (zh) 选择对象集的视点
CN110134982A (zh) 可通过铣削操作制造的部件的设计
EP3321817A1 (en) Querying a database based on a parametric view function
US20170371948A1 (en) Querying A Database With Morphology Criterion
US11604902B2 (en) Modeling using a weak type definition
US20230177229A1 (en) Partitioning 3d cad model
US20200211296A1 (en) Flexible modeling using a weak type definition
JP2021131852A (ja) 機械部品の構造シミュレーション
JP2021111376A (ja) 製品の物理的なプロトタイプの３ｄモデル化オブジェクト

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