CN115408796A - 一种多孔结构建模方法、装置、设备及应用 - Google Patents

一种多孔结构建模方法、装置、设备及应用 Download PDF

Info

Publication number
CN115408796A
CN115408796A CN202211084299.4A CN202211084299A CN115408796A CN 115408796 A CN115408796 A CN 115408796A CN 202211084299 A CN202211084299 A CN 202211084299A CN 115408796 A CN115408796 A CN 115408796A
Authority
CN
China
Prior art keywords
axis
target part
along
coordinate
scaling factor
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
CN202211084299.4A
Other languages
English (en)
Inventor
李文利
吕先林
邢占文
刘卫卫
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.)
Suzhou University
Original Assignee
Suzhou University
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 Suzhou University filed Critical Suzhou University
Priority to CN202211084299.4A priority Critical patent/CN115408796A/zh
Publication of CN115408796A publication Critical patent/CN115408796A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • 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
    • 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
    • G06T17/205—Re-meshing
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T2200/00—Indexing scheme for image data processing or generation, in general
    • G06T2200/04—Indexing scheme for image data processing or generation, in general involving 3D image data

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)
  • Computer Graphics (AREA)
  • Software Systems (AREA)
  • Prostheses (AREA)

Abstract

本发明公开了一种多孔结构建模方法、装置、设备及应用，涉及计算机辅助设计和增材制造技术领域，包括根据目标零件结构构建长方体区域，分别沿x轴、y轴改变等值面顶点之间的距离，完成等值面沿各个方向轮廓的修改，基于修改后的等值面构建目标零件，保证了多孔支架边界处仍保持极小曲面单胞的完整性，实现了利用极小曲面建立边界完整复杂多孔模型，为具有复杂轮廓的多孔支架设计提供了新方向。

Description

一种多孔结构建模方法、装置、设备及应用
技术领域
本发明涉及计算机辅助设计和增材制造技术领域，特别是涉及一种多孔结构建模方法、装置、设备及应用。
背景技术
为了获得理想生物医用支架所需的性能，多孔结构脱颖而出，其具有比表面积大，密度较小等特点，可以避免应力遮挡等问题。同时，还可以对多孔结构的表面进行修饰或进行细胞、生长因子等活性物质的填充。随着3D打印的发展，使得突破基于减材制造的传统设计方法，实现个性化多孔结构设计、制造、应用成为可能。借助点阵结构的优异机械力学性能，研究者开展了大量设计、制造尝试及性能评价。但是在应力分布分析及增材制造成型过程中，传统的点阵结构展现出承受外加载荷下过早应力集中，以及由于3D打印技术原理约束引发早期失效等问题。
极小曲面(Triply Periodic Minimal Surface，TPMS)结构因为其平均曲率为零的特点，能够降低应用集中，且内部贯通率为100％，能够在不牺牲孔隙率的前提下实现更高的减重目标值，是一种性能优越的多孔结构模型。同时，该结构在增材制造过程中能够保持稳定的自支撑能力以及较为一致的横截面积，具有更佳可打印性。
对于X,Y,Z三个方向上的周期T相同的情况下，极小曲面结构为正方体域内单胞。然而对于实际应用，需要建立具有复杂轮廓外形的多孔支架，现有的方法为：一是支架与极小曲面结构进行布尔求交得到多孔支架；二是对支架进行六面体网格划分并且利用等参单元映射。然而上述方法具有较大的局限性。对于第一种方法，布尔运算必然会导致支架外轮廓边界处不完整单胞的产生，而这种边界单胞的不完整性，会对支架承载产生不良影响。对于第二种方法，其一为对于复杂模型进行合适单胞大小的全六面体网格划分较为困难；其二是并非所有类型的TPMS都能通过这种方法得到复杂模型。
如图1为周期T为2的P，D，G三种单胞，由于P，D单元六个端面的轮廓和尺寸一致，如图2和图3以及图4分别为P，D，G的中心切片,|x|＝1截面处切片,|y|＝1截面处切片,|z|＝1截面处切片。由其截面切片可以看到，P，D截面轮廓形状一致，其中任意一个截面均可以和其它五个截面连接，形成连续的结构。而对于G型极小曲面，其各个截面的轮廓不一致，例如4中x＝1截面和x＝-1截面轮廓一致，而和其它四个截面不一致。这种类型的单胞和周围的单胞连接形成一个整体的时候，由于每个单胞摆放方式会受周围单胞的限制，所以每个单胞只有一种摆放方式才能和其它单胞连接起来，而对于P,D型极小曲面,任意摆放方式都能和周围单胞进行连接。在形函数映射的过程中，有些六面体单元和周围单元所连接的截面会发生变化，例如可能y＝1截面需要和z＝1或者x＝1截面连接，这样就会导致G型单胞和周围单胞断开，无法连接成为整体，因此无法采用六面体网格结合形状函数方式映射获得G型多孔支架。对于六个截面的形状轮廓和方向完全相同的的P，D单胞，即使连接截面发生变化，也能和周围单胞进行连接。采用圆形截面网格，分别利用G，P进行映射。如图5可以看到，对于P单胞，即使邻近单元连接截面发生变化之后仍连接良好。
由图5展示的G单胞而言，在有些位置处某一单胞和周围单胞的连接截面会发生变化，而这种变化导致其不能和周围单胞正确连接起来。因此，对于利用六面体网格划分结合形状函数进行单元变换方式进行多孔结构建模存在两个限制：(1)复杂模型难以进行高质量的六面体网格的划分；(2)利用形状函数变换对于G单胞并不适用。因此，该方法对TPMS单胞类型的选择也具有局限性。
综上所述可以看出，如何找到适合所有类型极小曲面结构单胞，并且边界具有完整单胞的多孔支架是目前有待解决的问题。
发明内容
本发明的目的是提供一种多孔结构建模方法、装置及设备，以解决现有技术中没有适合所有类型的极小曲面结构，且针对不同类型结构无法保证边界单胞完整性的问题。
为解决上述技术问题，本发明提供一种多孔结构建模方法，包括：
S1.基于目标零件结构，确定尺寸大小与所述目标零件结构最小包围范围的长方体区域，将尺寸最大的方向作为z轴，另外两个方向分别作为x轴与y轴，建立所述长方体区域内第一等值面；
S2.沿所述y轴将所述第一等值面点的y坐标乘以第一缩放系数修改所述第一等值面中顶点间距离，得到第二等值面；
S3.沿所述x轴将所述第二等值面点的x坐标乘以第二缩放系数修改所述第二等值面中顶点间距离，得到第三等值面；
S4.判断沿z轴方向轮廓是否变换，若沿所述z轴方向轮廓截面一致，则基于所述第三等值面生成目标零件，若沿所述z轴方向轮廓截面均匀变换，根据不同z坐标对截面进行缩放，生成目标零件，若沿所述z轴方向轮廓截面大小形状均变换，则根据不同z坐标重复步骤S2、S3、S4，直到得到截面大小和形状均发生变换的第四等值面，利用所述第四等值面生成目标零件。
优选地，所述建立所述长方体区域内第一等值面后包括：
基于所述目标零件结构，确定其内部长方形区域，确保长方形的两个点，或者四个点落在轮廓上，由长方形边长结合z向高度，确定极小曲面单胞尺寸。
优选地，所述沿所述y轴将所述第一等值面点的y坐标乘以第一缩放系数修改所述第一等值面中顶点间距离，得到第二等值面包括：
所述第一缩放系数由通过所述第一等值面并与y轴平行的直线在所述长方形边界内侧线段的距离除以长方形的边长确定。
优选地，所述沿所述x轴将所述第二等值面点的x坐标乘以第二缩放系数修改所述第二等值面中顶点间距离，得到第三等值面包括：
所述第二缩放系数由通过所述第二等值面并与x轴平行的直线在所述长方形边界内侧线段的距离除以长方形的边长确定。
优选地，所述第一缩放系数与所述第二缩放系数计算公式为：
Figure BDA0003834803630000031
其中，K1为第一缩放系数，D1为目标零件截面轮廓曲线与直线X＝x的割线位于轮廓内部长度，b为第一等值面所在长方体区域沿y轴方向长度；
Figure BDA0003834803630000041
其中，K2为第二缩放系数，D2为目标零件截面轮廓曲线与直线Y＝y的割线位于轮廓内部长度，c为第二等值面所在长方体区域沿x轴方向长度，Y1＝K1*Y，Y1为原始y坐标变换后的坐标，Y为原始等值面的y坐标。
优选地，所述第一等值面为所述长方体内的点组成的极小曲面。
优选地，所述第一等值面建模公式为：
F＝sin(ax)*cos(ay)+sin(ay)*cos(az)+sin(az)*cos(ax)
其中，F为第一等值面建模表达式，a为周期，x，y，z为函数变量。
本发明还提供一种多孔结构建模装置，包括：
模型建立模块，用于基于目标零件结构，确定尺寸大小与所述目标零件结构匹配的长方体区域，将尺寸最大的方向作为z轴，另外两个方向分别作为x轴与y轴，建立所述长方体区域内第一等值面；
y轴轮廓修改模块，用于沿所述y轴将所述第一等值面点的y坐标乘以第一缩放系数修改所述第一等值面中顶点间距离，得到第二等值面；
x轴轮廓修改模块，用于沿所述x轴将所述第二等值面点的x坐标乘以第二缩放系数修改所述第二等值面中顶点间距离，得到第三等值面；
目标零件生成模块，用于判断沿z轴方向轮廓是否变换，若沿所述z轴方向轮廓截面一致，则基于所述第三等值面生成目标零件，若沿所述z轴方向轮廓截面均匀变换，根据不同z坐标对截面进行缩放，生成目标零件，若沿所述z轴方向轮廓截面大小形状均变换，则根据不同z坐标重复步骤S2、S3、S4，直到得到截面大小和形状均发生变换的第四等值面，利用所述第四等值面生成目标零件。
本发明还提供一种多孔结构建模设备，包括：
存储器，用于存储计算机程序；
处理器，用于执行所述计算机程序时实现上述所述一种多孔结构建模方法的步骤。
本发明还提供一种如上述所述的多孔结构建模方法在生物医用支架的应用。
本发明所提供的一种多孔结构建模方法，根据目标零件结构构建长方体区域，分别沿x轴、y轴改变等值面顶点之间的距离，完成等值面沿各个方向轮廓的修改，基于修改后的等值面构建目标零件，保证了多孔支架边界处仍保持极小曲面单胞的完整性，实现了利用极小曲面建立边界完整复杂多孔模型，为具有复杂轮廓的多孔支架设计提供了新方向。
附图说明
为了更清楚的说明本发明实施例或现有技术的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单的介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为P，D，G三种单胞结构图；
图2为P单胞截面切片图；
图3为D单胞截面切片图；
图4为G单胞截面切片图；
图5为六面体网格及G,P单胞映射图；
图6为本发明所提供的一种多孔结构建模方法的第一种具体实施例的流程图；
图7为一种多孔结构建模方法的流程图；
图8为长方体区域内等值面图；
图9为两条水平边变为圆弧的等值面图；
图10为两条竖直边变为圆弧后的等值面图；
图11为对等值面光顺偏移得到具有厚度的圆柱体STL模型图；
图12为调整两条水平边间距离图；
图13为调整长方体区域为正六边形图；
图14为调整等值面沿z向截面大小图；
图15为对等值面光顺偏移得到具有厚度的角锥STL模型图；
图16为本发明实施例提供的一种多孔结构建模装置的结构框图。
具体实施方式
本发明的核心是提供一种多孔结构建模方法、装置、设备及应用，实现了利用极小曲面建立边界完整复杂多孔模型。
为了使本技术领域的人员更好地理解本发明方案，下面结合附图和具体实施方式对本发明作进一步的详细说明。显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
请参考图6、图7，图6为本发明所提供的一种多孔结构建模方法的第一种具体实施例的流程图；具体操作步骤如下：
S1.基于目标零件结构，确定尺寸大小与所述目标零件结构最小包围范围的长方体区域，将尺寸最大的方向作为z轴，另外两个方向分别作为x轴与y轴，建立所述长方体区域内第一等值面；
确定长方体区域，首先调整目标零件的位置，使尺寸最大的方向作为z轴，另外两个方向分别作为x轴与y轴，根据零件x，y轮廓确定位于其内部的长方形区域。该步骤需要根据实际截面确定长方形位置和边长。例如圆截面，则为其内接正方形；对于正六面体，则选择其边长的正方形区域；对于更加复杂的轮廓，确保长方形的两个点，或者四个点落在轮廓上即可。由长方向边长结合z向高度，确定一个长方体区域。由确定的长方体区域的边长确定极小曲面单胞尺寸，利用MATLAB软件中的meshgrid与isosurface函数建立长方体区域内等值面；
所述长方形区域由目标零件的最小包围区域，或者规则模型的内接长方形，并且可以根据极小曲面合适的单胞尺寸以及周期个数对长方形的长宽做出合适的调整，仍能保证使用该方法得到目标模型；
F＝sin(ax)*cos(ay)+sin(ay)*cos(az)+sin(az)*cos(ax)
其中，F为第一等值面建模表达式，a为周期，x，y，z为函数变量。
S2.沿所述y轴将所述第一等值面点的y坐标乘以第一缩放系数修改所述第一等值面中顶点间距离，得到第二等值面；
根据支架截面边界轮廓，沿y方向通过对等值面上所有点的y坐标乘以缩放系数K1修改等值面中顶点之间的距离，其缩放系数K1由通过该点并且与y平行的直线在长方形边界内侧线段的距离除以正方形的边长确定。通过MATLAB软件给予不同位置的顶点不同的缩放系数K1，实现对等值面的沿y方向轮廓的修改；
Figure BDA0003834803630000071
其中，K1为第一缩放系数，D1为目标零件截面轮廓曲线与直线X＝x的割线位于轮廓内部长度，b为第一等值面所在长方体区域沿y轴方向长度；
S3.沿所述x轴将所述第二等值面点的x坐标乘以第二缩放系数修改所述第二等值面中顶点间距离，得到第三等值面；
根据支架截面边界轮廓，沿x方向通过对等值面上部分顶点的x坐标乘以缩放系数K2修改等值面中顶点之间的距离，操作对象并不是所有顶点，而是第二步之后，仍保留在正方形区域内的所有顶点，这些顶点的缩放系数K2由通过该点并且与x平行的直线在长方形边界内侧线段的距离除以正方形的边长确定。通过MATLAB软件给予不同位置的顶点不同的缩放系数K2，实现对等值面的沿x轮廓的修改；
Figure BDA0003834803630000072
其中，K2为第二缩放系数，D2为目标零件截面轮廓曲线与直线Y＝y的割线位于轮廓内部长度，c为第二等值面所在长方体区域沿x轴方向长度，Y1＝K1*Y，Y1为原始y坐标变换后的坐标，Y为原始等值面的y坐标。
S4.判断沿z轴方向轮廓是否变换，若沿所述z轴方向轮廓截面一致，则基于所述第三等值面生成目标零件，若沿所述z轴方向轮廓截面均匀变换，根据不同z坐标对截面进行缩放，生成目标零件，若沿所述z轴方向轮廓截面大小形状均变换，则根据不同z坐标重复步骤S2、S3、S4，直到得到截面大小和形状均发生变换的第四等值面，利用所述第四等值面生成目标零件；
对于该种z截面一致的模型，不进行操作，跳过此步骤。
对于z向各个截面大小发生变化而轮廓不变，即z截面为相似图形的截面：
对于该种模型，只需要根据线性插值来确定每个z截面的缩放系数，来对第三等值面的点的x，y坐标进行缩放操作，得到第四等值面，通过最大截面面积S1,最小截面面积S2,以及当前截面面积S,通过插值方法确定缩放系数K3；
Figure BDA0003834803630000081
其中，z为截面的z坐标，X2，Y2为等到的第四等值面x，y坐标，X2＝K3*X1，Y2＝K3*Y1。
对于这种复杂模型，截面既不相同也不相似，无法采用缩放方式来调整截面轮廓。故需要对模型沿z进行切片处理，选择层厚t，对模型切片，再对每一个切片重复第二步第三步操作，最后把每个截面的点再组合起来得到新的第四等值面。
本实施例提供一种多孔结构建模方法，针对极小曲面建立的复杂轮廓的多孔结构，由长方形区域等值面，通过改变顶点之间的距离，实现对复杂轮廓建模的新方法，保证了多孔支架边界处仍保持极小曲面单胞的完整性，实现了利用极小曲面建立边界完整复杂多孔模型。
基于上述实施例，本实施例利用具体实验对本方法进行阐述，具体操作步骤如下：
①直径为10mm，高为15mm的圆柱体多孔支架；
第一步：如图8所示，根据直径10mm,确定单胞大小为2.5mm,并且长方体区域内x,y方向为三个周期单胞，z方向为六个周期单胞。长方体区域内尺寸为7.5mm*7.5mm*15mm，通过MATLAB采用meshgrid函数划分三维网格，并生成三维体数据，利用isosurface函数提取该区域内Gyroid类型单胞等值面；
第二步：如图9所示，以第一步中生成的等值面中所有点的x坐标为判断条件，以通过x处的圆的弦长比上正方形边长作为缩放系数K1，对所有点的y坐标乘以缩放系数K1，以对顶点之间的距离进行调整，将其两条水平边变形为两段圆弧。
第三步：如图10所示，以第二步中生成的等值面中所有坐标仍在1中长方体区域内的顶点的y坐标为判断条件，对所有点的x坐标乘以缩放系数K2，将其两条垂直边变形为两段圆弧。
第四步：将圆直径
Figure BDA0003834803630000091
缩放至10，缩放系数K3为直径10与
Figure BDA0003834803630000092
的比值，对所有x,y乘以缩放系数K3,得到最终等值面。
第五步：如图11所示，以STL形式导出第四步中等值面，进行后续光顺偏移处理用于3D打印。
②上底直径为8mm,下底直径为10mm,高为10mm，棱边数为6的角锥；
第一步：根据角锥截面确定长方体域内G型等值面分别为x范围为-2.5～2.5，y向为-3.75～3.75，z向为-5～5，即为5mm*7.5mm*10mm，通过MATLAB采用meshgrid函数划分三维网格，并生成三维体数据，利用isosurface函数提取该区域内G型单胞等值面。
第二步：如图12所示，以第一步中生成的等值面中所有点的x坐标为判断条件，对所有点的y坐标乘以缩放系数K1进行调整，将其两条水平边距离由7.5变为
Figure BDA0003834803630000093
第三步：如图13所示，以第二步中生成的等值面中所有点的y标为判断条件，对所有点的x坐标乘以缩放系数K2进行调整，改变其两条水平边距离，将截面变为正六边形。
第四步：如图14所示，截面大小沿z方向由下底至上底，直径由10mm变为8mm，对截面进行缩放，方法为所有点x坐标，y坐标分别乘以相应截面的缩放系数
Figure BDA0003834803630000101
采用该缩放系数对x,y坐标进行缩小处理，得到第三等值面，最终得到截面沿z变化的模型。
第五步：如图15所示，以STL形式导出第四步中等值面，进行后续光顺偏移处理用于3D打印。
本实施例提供一种多孔结构建模方法，利用具体数据对本方法进行实验操作，根据目标零件结构构建长方体区域，分别沿x轴、y轴改变等值面顶点之间的距离，完成等值面沿各个方向轮廓的修改，基于修改后的等值面构建目标零件，保证了多孔支架边界处仍保持极小曲面单胞的完整性，实现了利用极小曲面建立边界完整复杂多孔模型，为具有复杂轮廓的多孔支架设计提供了新方向。
请参考图16，图16为本发明实施例提供的一种多孔结构建模装置的结构框图；具体装置可以包括：
模型建立模块100，用于基于目标零件结构，确定尺寸大小与所述目标零件结构匹配的长方体区域，将尺寸最大的方向作为z轴，另外两个方向分别作为x轴与y轴，建立所述长方体区域内第一等值面；
y轴轮廓修改模块200，用于沿所述y轴将所述第一等值面点的y坐标乘以第一缩放系数修改所述第一等值面中顶点间距离，得到第二等值面；
x轴轮廓修改模块300，用于沿所述x轴将所述第二等值面点的x坐标乘以第二缩放系数修改所述第二等值面中顶点间距离，得到第三等值面；
目标零件生成模块400，用于判断沿z轴方向轮廓是否变换，若沿所述z轴方向轮廓截面一致，则基于所述第三等值面生成目标零件，若沿所述z轴方向轮廓截面均匀变换，根据不同z坐标对截面进行缩放，生成目标零件，若沿所述z轴方向轮廓截面大小形状均变换，则根据不同z坐标重复步骤S2、S3、S4，直到得到截面大小和形状均发生变换的第四等值面，利用所述第四等值面生成目标零件。
本实施例的一种多孔结构建模装置用于实现前述的一种多孔结构建模方法，因此一种多孔结构建模装置中的具体实施方式可见前文中的一种多孔结构建模方法的实施例部分，例如，模型建立模块100，y轴轮廓修改模块200，x轴轮廓修改模块300，目标零件生成模块400，分别用于实现上述一种多孔结构建模方法中步骤S1，S2，S3和S4，所以，其具体实施方式可以参照相应的各个部分实施例的描述，在此不再赘述。
本发明具体实施例还提供了一种多孔结构建模设备，包括：存储器，用于存储计算机程序；处理器，用于执行所述计算机程序时实现上述一种多孔结构建模方法的步骤。
本发明具体实施例还提供了一种多孔结构建模方法在生物医用支架的应用。
本说明书中各个实施例采用递进的方式描述，每个实施例重点说明的都是与其它实施例的不同之处，各个实施例之间相同或相似部分互相参见即可。对于实施例公开的装置而言，由于其与实施例公开的方法相对应，所以描述的比较简单，相关之处参见方法部分说明即可。
专业人员还可以进一步意识到，结合本文中所公开的实施例描述的各示例的单元及算法步骤，能够以电子硬件、计算机软件或者二者的结合来实现，为了清楚地说明硬件和软件的可互换性，在上述说明中已经按照功能一般性地描述了各示例的组成及步骤。这些功能究竟以硬件还是软件方式来执行，取决于技术方案的特定应用和设计约束条件。专业技术人员可以对每个特定的应用来使用不同方法来实现所描述的功能，但是这种实现不应认为超出本发明的范围。
结合本文中所公开的实施例描述的方法或算法的步骤可以直接用硬件、处理器执行的软件模块，或者二者的结合来实施。软件模块可以置于随机存储器(RAM)、内存、只读存储器(ROM)、电可编程ROM、电可擦除可编程ROM、寄存器、硬盘、可移动磁盘、CD-ROM、或技术领域内所公知的任意其它形式的存储介质中。
以上对本发明所提供的一种多孔结构建模方法、装置、设备及应用进行了详细介绍。本文中应用了具体个例对本发明的原理及实施方式进行了阐述，以上实施例的说明只是用于帮助理解本发明的方法及其核心思想。应当指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下，还可以对本发明进行若干改进和修饰，这些改进和修饰也落入本发明权利要求的保护范围内。

Claims (10)

1.一种多孔结构建模方法，其特征在于，包括：
S1.基于目标零件结构，确定尺寸大小与所述目标零件结构最小包围范围的长方体区域，将尺寸最大的方向作为z轴，另外两个方向分别作为x轴与y轴，建立所述长方体区域内第一等值面；
S2.沿所述y轴将所述第一等值面点的y坐标乘以第一缩放系数修改所述第一等值面中顶点间距离，得到第二等值面；
S3.沿所述x轴将所述第二等值面点的x坐标乘以第二缩放系数修改所述第二等值面中顶点间距离，得到第三等值面；
S4.判断沿z轴方向轮廓是否变换，若沿所述z轴方向轮廓截面一致，则基于所述第三等值面生成目标零件，若沿所述z轴方向轮廓截面均匀变换，根据不同z坐标对截面进行缩放，生成目标零件，若沿所述z轴方向轮廓截面大小形状均变换，则根据不同z坐标重复步骤S2、S3、S4，直到得到截面大小和形状均发生变换的第四等值面，利用所述第四等值面生成目标零件。
2.如权利要求1所述的多孔建模方法，其特征在于，所述建立所述长方体区域内第一等值面后包括：
基于所述目标零件结构，确定其内部长方形区域，确保长方形的两个点，或者四个点落在轮廓上，由长方形边长结合z向高度，确定极小曲面单胞尺寸。
3.如权利要求2所述的多孔建模方法，其特征在于，所述沿所述y轴将所述第一等值面点的y坐标乘以第一缩放系数修改所述第一等值面中顶点间距离，得到第二等值面包括：
所述第一缩放系数由通过所述第一等值面并与y轴平行的直线在所述长方形边界内侧线段的距离除以长方形的边长确定。
4.如权利要求2所述的多孔建模方法，其特征在于，所述沿所述x轴将所述第二等值面点的x坐标乘以第二缩放系数修改所述第二等值面中顶点间距离，得到第三等值面包括：
所述第二缩放系数由通过所述第二等值面并与x轴平行的直线在所述长方形边界内侧线段的距离除以长方形的边长确定。
5.如权利要求1所述的多孔建模方法，其特征在于，所述第一缩放系数与所述第二缩放系数计算公式为：
Figure FDA0003834803620000021
其中，K1为第一缩放系数，D1为目标零件截面轮廓曲线与直线X＝x的割线位于轮廓内部长度，b为第一等值面所在长方体区域沿y轴方向长度；
Figure FDA0003834803620000022
其中，K2为第二缩放系数，D2为目标零件截面轮廓曲线与直线Y＝y的割线位于轮廓内部长度，c为第二等值面所在长方体区域沿x轴方向长度，Y1＝K1*Y，Y1为原始y坐标变换后的坐标，Y为原始等值面的y坐标。
6.如权利要求1所述的多孔建模方法，其特征在于，所述第一等值面为所述长方体内的点组成的极小曲面。
7.如权利要求6所述的多孔建模方法，其特征在于，所述第一等值面建模公式为：
F＝sin(ax)*cos(ay)+sin(ay)*cos(az)+sin(az)*cos(ax)
其中，F为第一等值面建模表达式，a为周期，x，y，z为函数变量。
8.一种多孔结构建模装置，其特征在于，包括：
模型建立模块，用于基于目标零件结构，确定尺寸大小与所述目标零件结构匹配的长方体区域，将尺寸最大的方向作为z轴，另外两个方向分别作为x轴与y轴，建立所述长方体区域内第一等值面；
y轴轮廓修改模块，用于沿所述y轴将所述第一等值面点的y坐标乘以第一缩放系数修改所述第一等值面中顶点间距离，得到第二等值面；
x轴轮廓修改模块，用于沿所述x轴将所述第二等值面点的x坐标乘以第二缩放系数修改所述第二等值面中顶点间距离，得到第三等值面；
目标零件生成模块，用于判断沿z轴方向轮廓是否变换，若沿所述z轴方向轮廓截面一致，则基于所述第三等值面生成目标零件，若沿所述z轴方向轮廓截面均匀变换，根据不同z坐标对截面进行缩放，生成目标零件，若沿所述z轴方向轮廓截面大小形状均变换，则根据不同z坐标重复步骤S2、S3、S4，直到得到截面大小和形状均发生变换的第四等值面，利用所述第四等值面生成目标零件。
9.一种多孔结构建模设备，其特征在于，包括：
存储器，用于存储计算机程序；
处理器，用于执行所述计算机程序时实现如权利要求1至7任一项所述一种多孔结构建模方法的步骤。
10.一种如权利要求1-7任一项所述的多孔结构建模方法在生物医用支架的应用。
CN202211084299.4A 2022-09-06 2022-09-06 一种多孔结构建模方法、装置、设备及应用 Pending CN115408796A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211084299.4A CN115408796A (zh) 2022-09-06 2022-09-06 一种多孔结构建模方法、装置、设备及应用

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211084299.4A CN115408796A (zh) 2022-09-06 2022-09-06 一种多孔结构建模方法、装置、设备及应用

Publications (1)

Publication Number Publication Date
CN115408796A true CN115408796A (zh) 2022-11-29

Family

ID=84163811

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211084299.4A Pending CN115408796A (zh) 2022-09-06 2022-09-06 一种多孔结构建模方法、装置、设备及应用

Country Status (1)

Country Link
CN (1) CN115408796A (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116386788A (zh) * 2023-04-10 2023-07-04 精创石溪科技(成都)有限公司 一种基于多目标优化的变密度多孔网格结构参数化建模方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116386788A (zh) * 2023-04-10 2023-07-04 精创石溪科技(成都)有限公司 一种基于多目标优化的变密度多孔网格结构参数化建模方法
CN116386788B (zh) * 2023-04-10 2023-11-21 精创石溪科技(成都)有限公司 一种基于多目标优化的变密度多孔网格结构参数化建模方法

Similar Documents

Publication Publication Date Title
Shimada Physically-based mesh generation: automated triangulation of surfaces and volumes via bubble packing
CN111797555B (zh) 一种基于有限元模型的几何重构方法
AU2018203479B2 (en) High fatigue strength porous structure
US10121279B1 (en) Systems and methods for generating a mesh
CN107886569B (zh) 一种基于离散李导数的测度可控的曲面参数化方法及系统
CN106650085B (zh) 基于有限元网格的实体模型生成镂空模型的方法
CN111581776B (zh) 一种基于几何重建模型的等几何分析方法
Zhang et al. Adaptive hexahedral mesh generation based on local domain curvature and thickness using a modified grid-based method
CN115408796A (zh) 一种多孔结构建模方法、装置、设备及应用
CN110188395B (zh) 一种基于线面体的维度增加式计算流体网格生成方法
CN108389263B (zh) 面向边界元方法求解计算的iges面网格快速生成方法
CN111079326B (zh) 二维各向异性网格单元度量张量场光滑化方法
Choi et al. Support-free hollowing with spheroids and efficient 3D printing utilizing circular printing motions based on Voronoi diagrams
CN115994410A (zh) 基于八叉树细化四面体网格的飞行器仿真驱动设计方法
Zhang et al. Adaptive generation of hexahedral element mesh using an improved grid-based method
CN111079353B (zh) 应用于复杂流体分析的快速均匀网格划分的方法及装置
CN110796735A (zh) Nurbs曲面有限元板壳网格划分方法及计算机实现系统
CN115170734A (zh) 一种基于截面切片的三维回转结构重构方法及系统
CN111047687B (zh) 一种基于三维t样条的异质材料实体建模方法
CN115130340A (zh) 基于分数布朗运动的管道建模方法
Lin et al. Trivariate B-spline solid construction by pillow operation and geometric iterative fitting
CN113806951A (zh) 一种基于半边数据结构的自然邻近点搜索的弹性仿真方法
CN110765506A (zh) 实体模型的多分辨率等几何拓扑优化方法
CN116992513B (zh) 一种适应复杂几何的模拟仿真软件贴体粒子生成方法
Zhao et al. Geometry-adaptive generation algorithm and boundary match method for initial hexahedral element mesh

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