CN108804821A - 一种基于空间网格的三维点阵结构设计方法 - Google Patents
一种基于空间网格的三维点阵结构设计方法 Download PDFInfo
- Publication number
- CN108804821A CN108804821A CN201810595037.1A CN201810595037A CN108804821A CN 108804821 A CN108804821 A CN 108804821A CN 201810595037 A CN201810595037 A CN 201810595037A CN 108804821 A CN108804821 A CN 108804821A
- Authority
- CN
- China
- Prior art keywords
- lattice
- design
- lattice structure
- memeber
- dimensional
- 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
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
- 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/17—Mechanical parametric or variational design
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
Abstract
本发明一种基于空间网格的三维点阵结构设计方法:1、导入结构初始三维实体模型;2、分解结构并选择设计区域;3、选择功能表面,删除设计区域实体与非功能表面;4、确定所有子区域点阵结构的胞元构型及胞元密度;5、在子区域生成点阵结构拓扑信息;6、构成三维点阵结构信息并计算杆件数量;7、输入增材制造工艺要求参数;8、对所有杆件进行有效性判断;9、计算失效杆件数量;10、计算失效杆件比;11、判断全局有效性;12、构建点阵结构几何模型;13、对非设计区域、功能表面与点阵结构进行合并运算,得到包含点阵结构的设计模型。本发明设计结果与增材制造工艺要求及数据接口匹配良好;采用批处理,避免较高技巧和繁琐工作。
Description
技术领域:
本发明涉及一种基于空间网格的三维点阵结构设计方法,属于机械结构设计的技术领域。
背景技术:
作为航空、航天、汽车、医疗等多个领域的新型结构类型,类桁架式点阵结构因其轻质、高比刚度、高比强度的机械特性,减震、吸能等功能特征而受到广泛关注与应用。三维点阵结构通常采用设计区域内直接填充的方法,该方法设计自由度小,设计师参与程度低,设计人员无法对点阵结构胞元构型、胞元密度、分布形式以及结构外观进行掌控,另一方面,直接填充的方法对设计区域的边缘位置处理不足,通常会产生不完整的胞元结构或者独立的杆件,形成结构上的缺陷。工程上希望设计人员可以对设计区域内的点阵结构构型、密度等多个方面进行控制,并希望设计方法可以保证设计区域内设计产生规则且完整的三维点阵结构。由于点阵结构几何特征复杂,工程上常采用增材制造技术生产加工点阵结构,因此工程上需考虑到增材制造的工艺特点,保证结构可实现性。另一方面,工程上点阵结构通常不会独立作为结构出现,需要配合面板、蒙皮、实体等特征才能形成完整的结构,因此在结构外观基本确定的情况下,设计方法应提供点阵结构设计区域及功能表面的选择,结构的设计结果中可以同时包括点阵结构与曲面、实体等特征。
发明内容:
为了实现三维点阵结构的交互式设计,本发明的目的是提供一种基于空间网格的三维点阵结构设计方法,该方法利用有限元分析软件平台,首先导入结构初始三维实体模型,然后通过人机交互的方式对结构进行划分,之后由设计人员选择点阵结构设计区域与功能表面。然后由设计人员制定胞元构型以及胞元密度,随后利用空间网格划分的算法,对设计区域划分网格并提取网格信息。根据网格信息建立点阵结构并根据增材制造工艺要求对点阵结构进行有效性验证。之后将非设计区域与功能表面与点阵结构进行几何合并布尔运算,形成完整的结构模型。本方法人机交互操作少,设计人员设定好初始需求后,网格划分与点阵结构模型建立过程可以实现完全自动处理,通用性强,适用范围广。
本发明为一种基于空间网格的三维点阵结构设计方法,包括以下步骤:
步骤1:在有限元软件平台内导入结构初始三维实体模型;
步骤2:分解结构并选择设计区域;通过人机交互方式选择若干设计子区域。
步骤3:通过人机交互的方式选择功能表面,删除设计区域实体与非功能表面;
步骤4:通过人机交互的方式输入所有子区域点阵结构的胞元构型以及胞元密度;
步骤5:在子区域内生成空间网格并产生点阵结构拓扑信息;
步骤6:基于空间网格信息构成三维点阵结构信息。通过遍历的方式,计算三维点阵结构信息中所包含的所有点阵杆件数量N。
步骤7:通过人机交互的方式输入增材制造工艺要求参数;
步骤8:对所有杆件进行有效性判断;计算各杆件特征并根据工艺参数,判定杆件是否有效。
步骤9:计算失效杆件数量Nf:遍历所有杆件,若杆件被判定为失效,则失效杆件数量计数加1。
步骤10:计算失效杆件比R;失效杆件比R为失效杆件数量Nf与总体杆件数量N的比值,计算公式为
步骤11:若R<Rmax,执行步骤12,否则返回步骤4。
步骤12:根据点阵结构信息,构建点阵结构几何模型;
步骤13:对非设计区域、功能表面与点阵结构进行合并运算,得到包含点阵结构的设计模型。
本发明的有益效果是:提供了一种基于空间网格的三维点阵结构设计方法,通过该方法可以在原有三维实体的基础上,设计并生成边缘完整,可设计性高的三维点阵结构模型。设计过程中考虑增材制造工艺要求,对点阵结构进行有效性验证,同时设计结果为实体与曲面、曲线相结合的抽象模型,该模型与增材制造设备匹配良好,从而保证结构可以利用增材制造技术生产加工。设计方法中采用算法批处理实现,减少了大量重复性人工操作。本发明在航空、航天、汽车、医疗等领域的机械结构设计中具有很强的工程实用价值。
附图说明
图1.本发明所述方法流程图。
图2.实施例初始模型。
图3.实施例结构分区与编号(主视方向)。
图4.实施例保留的实体与表面。
图5.杆件示意图。
图6.制造方向示意图。
图7.实施例区域1处的点阵结构(主视方向)。
图8.实施例设计结果。
图中符号说明如下:
N为杆件数量;
Nf为失效杆件数量;
R为失效比;
Rmax为失效可行值;
P1,P2为杆件起止点;
为杆件方向向量;
l为杆件长度;
为增材制造方向;
θ为单位方向与杆件向量的夹角;
具体实施方式
下面结合附图和实施例对本发明作详细说明。
本发明一种基于空间网格的三维点阵结构设计方法,见图1所示,它包括以下步骤:
步骤1:本专利所提出方法是在有限元软件平台内实现的,但不拘束于某一个特定的有限元软件。通过软件的模型导入功能即可以三维实体模型的导入,实施例初始模型如图2所示;
步骤2:通过几何分解功能对结构进行分解,通过人机交互的方式选择设计区域并统计设计子区域数目Nd。如图3所示,实施例中选择区域1,2,3为设计区域,区域4,5为非设计区域。设计区域或非设计区域可能分别包含多个子区域,统计设计子区域数量为Nd,实施例中设计子区域数量Nd=3;
步骤3:通过人机交互的方式选择功能表面。通过几何编辑功能删除设计区域实体与非功能表面,并保留非设计区域实体与功能表面。实施例中保留的实体与表面如图4所示;
步骤4:通过人机交互的方式输入所有子区域点阵结构的胞元构型以及胞元密度,构型及密度由设计人员确定。;
步骤5:在子区域内生成空间网格并产生点阵结构拓扑信息;
步骤6:基于点阵结构拓扑信息构成三维点阵结构杆件信息。通过遍历的方式,统计三维点阵结构信息中所包含的所有点阵杆件数量N。
步骤7:通过人机交互的方式输入增材制造工艺要求参数;
步骤8:计算杆件长度和夹角,根据向量模长计算方法,计算杆件长度公式为如图5所示。根据向量的数量积公式计算制造单位方向与杆件向量的夹角θ,由于则此处规定夹角θ满足:0≤θ≤90°,则有如图6所示。并对所有杆件进行有效性判断;
步骤9:遍历所有杆件,若杆件被判定为失效,则失效杆件数量计数加1;
步骤10:失效杆件比R为失效杆件数量Nf与总体杆件数量N的比值,计算公式为
步骤11:根据加工工艺需求和质量要求,确定Rmax。若R<Rmax,执行步骤12,否则返回步骤4。
步骤12:根据点阵结构信息的每一根杆件的起止点坐标,在有限元软件或三维建模软件中建立杆件特征,从而完成点阵结构的建立,实施例区域1处的点阵结构如图7所示;
步骤13:将步骤2选择的非设计区域及功能表面与点阵结构合并,所得几何模型即为设计结果,如图8所示。
通过以上各步骤,可以在原有三维实体的基础上,设计并生成边缘完整,可设计性高的三维点阵结构模型。设计过程中考虑增材制造工艺要求,对点阵结构进行有效性验证,同时设计结果为实体与曲面、曲线相结合的抽象模型,该模型与增材制造设备匹配良好,从而保证结构可以利用增材制造技术生产加工。设计方法中采用算法批处理实现,减少了大量重复性人工操作。
Claims (1)
1.一种基于空间网格的三维点阵结构设计方法,其特征在于:该方法包括以下步骤:
步骤1:在有限元软件平台内导入结构初始三维实体模型;
步骤2:分解结构并选择设计区域;通过人机交互方式选择若干设计子区域;
步骤3:通过人机交互的方式选择功能表面,删除设计区域实体与非功能表面;
步骤4:通过人机交互的方式输入所有子区域点阵结构的胞元构型以及胞元密度;
步骤5:在子区域内生成空间网格并产生点阵结构拓扑信息;
步骤6:基于空间网格信息构成三维点阵结构信息;通过遍历的方式,计算三维点阵结构信息中所包含的所有点阵杆件数量N;
步骤7:通过人机交互的方式输入增材制造工艺要求参数;
步骤8:对所有杆件进行有效性判断;计算各杆件特征并根据工艺参数,判定杆件是否有效;
步骤9:计算失效杆件数量Nf:遍历所有杆件,若杆件被判定为失效,则失效杆件数量计数加1;
步骤10:计算失效杆件比R;失效杆件比R为失效杆件数量Nf与总体杆件数量N的比值,计算公式为
步骤11:若R<Rmax,执行步骤12,否则返回步骤4;
步骤12:根据点阵结构信息,构建点阵结构几何模型;
步骤13:对非设计区域、功能表面与点阵结构进行合并运算,得到包含点阵结构的设计模型。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810595037.1A CN108804821A (zh) | 2018-06-11 | 2018-06-11 | 一种基于空间网格的三维点阵结构设计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810595037.1A CN108804821A (zh) | 2018-06-11 | 2018-06-11 | 一种基于空间网格的三维点阵结构设计方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108804821A true CN108804821A (zh) | 2018-11-13 |
Family
ID=64089037
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810595037.1A Pending CN108804821A (zh) | 2018-06-11 | 2018-06-11 | 一种基于空间网格的三维点阵结构设计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108804821A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110570526A (zh) * | 2019-09-06 | 2019-12-13 | 南京航空航天大学 | 拓扑信息引导的三维点阵结构件内部几何数据获取方法 |
CN114329770A (zh) * | 2021-12-14 | 2022-04-12 | 北京理工大学 | 基于等参变换的点阵结构随形填充式结构及其设计方法 |
CN114309658A (zh) * | 2021-11-15 | 2022-04-12 | 上海工程技术大学 | 一种基于非均匀点阵结构的增材制造方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104732028A (zh) * | 2015-03-30 | 2015-06-24 | 沈阳飞机工业(集团)有限公司 | 一种基于空间点阵的飞机结构件铺粉烧结支撑结构生成方法 |
CN105205213A (zh) * | 2015-08-24 | 2015-12-30 | 哈尔滨工业大学 | 一种点阵材料等效力学性能分析系统 |
CN107138726A (zh) * | 2017-05-12 | 2017-09-08 | 中国航发北京航空材料研究院 | 一种具有点阵冷却结构的导向叶片制备方法 |
CN107599445A (zh) * | 2017-09-08 | 2018-01-19 | 哈尔滨工程大学 | 一种复合材料点阵结构及制备方法 |
WO2018031904A1 (en) * | 2016-08-12 | 2018-02-15 | Siemens Product Lifecycle Management Software Inc. | Computer aided design with high resolution lattice structures using graphics processing units (gpu) |
CN107742014A (zh) * | 2017-09-29 | 2018-02-27 | 北京空间飞行器总体设计部 | 基于增材制造的相变储能装置点阵夹层结构的设计方法 |
CN107742047A (zh) * | 2017-11-10 | 2018-02-27 | 北京航空航天大学 | 一种变相对密度octet点阵结构的设计方法 |
-
2018
- 2018-06-11 CN CN201810595037.1A patent/CN108804821A/zh active Pending
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104732028A (zh) * | 2015-03-30 | 2015-06-24 | 沈阳飞机工业(集团)有限公司 | 一种基于空间点阵的飞机结构件铺粉烧结支撑结构生成方法 |
CN105205213A (zh) * | 2015-08-24 | 2015-12-30 | 哈尔滨工业大学 | 一种点阵材料等效力学性能分析系统 |
WO2018031904A1 (en) * | 2016-08-12 | 2018-02-15 | Siemens Product Lifecycle Management Software Inc. | Computer aided design with high resolution lattice structures using graphics processing units (gpu) |
CN107138726A (zh) * | 2017-05-12 | 2017-09-08 | 中国航发北京航空材料研究院 | 一种具有点阵冷却结构的导向叶片制备方法 |
CN107599445A (zh) * | 2017-09-08 | 2018-01-19 | 哈尔滨工程大学 | 一种复合材料点阵结构及制备方法 |
CN107742014A (zh) * | 2017-09-29 | 2018-02-27 | 北京空间飞行器总体设计部 | 基于增材制造的相变储能装置点阵夹层结构的设计方法 |
CN107742047A (zh) * | 2017-11-10 | 2018-02-27 | 北京航空航天大学 | 一种变相对密度octet点阵结构的设计方法 |
Non-Patent Citations (4)
Title |
---|
ZENG F , DENG H , ZHANG X , ET AL.: "Research on the Design Method of Lightweight Lattice Sandwich Structure for Additive Manufacturing", 《INTERNATIONAL CONFERENCE ON MECHANICAL DESIGN(ICMD 2017):ADVANCES IN MECHANICAL DESIGN》 * |
付志方; 王春洁: "节点位置不确定下桁架结构稳健拓扑优化", 《机械工程学报》 * |
冀宾; 顾铖璋; 韩涵; 宋林郁; 吴春雷: "几种变截面点阵承受面外压缩载荷的力学行为", 《固体力学学报》 * |
段尊义; 阎军; 张瑞; 赵国忠: "复合材料框架结构刚度几何双尺度优化设计", 《计算力学学报》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110570526A (zh) * | 2019-09-06 | 2019-12-13 | 南京航空航天大学 | 拓扑信息引导的三维点阵结构件内部几何数据获取方法 |
CN110570526B (zh) * | 2019-09-06 | 2021-02-12 | 南京航空航天大学 | 拓扑信息引导的三维点阵结构件内部几何数据获取方法 |
CN114309658A (zh) * | 2021-11-15 | 2022-04-12 | 上海工程技术大学 | 一种基于非均匀点阵结构的增材制造方法 |
CN114309658B (zh) * | 2021-11-15 | 2023-06-02 | 上海工程技术大学 | 一种基于非均匀点阵结构的增材制造方法 |
CN114329770A (zh) * | 2021-12-14 | 2022-04-12 | 北京理工大学 | 基于等参变换的点阵结构随形填充式结构及其设计方法 |
CN114329770B (zh) * | 2021-12-14 | 2024-07-02 | 北京理工大学 | 基于等参变换的点阵结构随形填充式结构及其设计方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108804821A (zh) | 一种基于空间网格的三维点阵结构设计方法 | |
CN112836411B (zh) | 加筋板壳结构的优化方法、装置、计算机设备和存储介质 | |
CN107748824A (zh) | 一种基于bim技术的混凝土框构桥设计方法 | |
CN104851126B (zh) | 基于广义圆柱体的三维模型分割方法及装置 | |
CN106994483A (zh) | 一种汽车覆盖件拉延模具精确型面加工的方法 | |
CN112818470B (zh) | 一种基结构的优化方法、装置、计算机设备和存储介质 | |
CN103150416B (zh) | 设计由深度元素表示的建模体 | |
CN110334450A (zh) | 一种多块结构网格生成中物面投影错误的修复方法 | |
CN109063405A (zh) | 一种基于形貌优化的空调顶盖板加强筋设计方法 | |
CN109101671A (zh) | 一种变密度与变构型三维点阵结构建模方法 | |
CN102930589B (zh) | 一种非结构直角网格相贯线修正方法 | |
Cui et al. | Shape design of curved surface of membrane structure using developable surface | |
CN111740424A (zh) | 一种基于运算树gpu并行加速模型的电力系统潮流计算方法 | |
CN116522660A (zh) | 一种用于爆炸问题的数值模拟方法及系统 | |
Hou et al. | Multilevel optimisation method for vehicle body in conceptual design | |
Ferlemann et al. | Parametric geometry, structured grid generation, and initial design study for REST-Class hypersonic inlets | |
JP7257472B2 (ja) | 設計支援システム、設計支援方法及び設計支援プログラム | |
CN113011014B (zh) | 一种复合材料铺层优化方法及系统 | |
CN114372426A (zh) | 一种基于ffd方法的隔离段的优化设计方法及系统 | |
CN104156502B (zh) | 一种基于位置的衣物褶皱几何生成方法 | |
Leidinger et al. | Explicit isogeometric b-rep analysis on trimmed nurbs-based multi-patch cad models in ls-dyna | |
Liu et al. | Heterogeneous CPU-GPU Accelerated Parallel Subgridding FDTD Algorithm | |
CN110502819A (zh) | 混凝土桥梁条形构件外封闭箍筋的计算机辅助设计方法 | |
Xie et al. | A CAD/CAE integration technology and its application in hydraulic engineering | |
Patel et al. | Automated Synthesis of Sheet Metal Parts by Optimizing a Fabrication Based Graph Topology |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20181113 |