CN109101671A - 一种变密度与变构型三维点阵结构建模方法 - Google Patents

一种变密度与变构型三维点阵结构建模方法 Download PDF

Info

Publication number
CN109101671A
CN109101671A CN201810593903.3A CN201810593903A CN109101671A CN 109101671 A CN109101671 A CN 109101671A CN 201810593903 A CN201810593903 A CN 201810593903A CN 109101671 A CN109101671 A CN 109101671A
Authority
CN
China
Prior art keywords
unit
face
establish
lattice structure
rod piece
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
Application number
CN201810593903.3A
Other languages
English (en)
Other versions
CN109101671B (zh
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.)
Beihang University
Original Assignee
Beihang 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 Beihang University filed Critical Beihang University
Priority to CN201810593903.3A priority Critical patent/CN109101671B/zh
Publication of CN109101671A publication Critical patent/CN109101671A/zh
Application granted granted Critical
Publication of CN109101671B publication Critical patent/CN109101671B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • G06F30/23Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/17Mechanical parametric or variational design

Abstract

一种变密度与变构型三维点阵结构建模方法,它包括以下步骤:步骤1:导入结构实体模型;步骤2:分解结构;步骤3:统计结构棱边数E和面数F;步骤4:对所有子结构计算点阵结构布局;步骤5:建立面内单元;步骤6:对所有杆件判定是否建立并建立杆上单元;步骤7:单元联结;步骤8:模型检查;通过以上各步骤,建立了共节点的杆、梁与壳单元,避免了刚性连接和多余约束,实现了变密度与变构型三维点阵结构有限元建模。该方法摆脱了几何模型的约束,避免了几何模型精度带来的建模失败的问题。同时减少了人为操作,提高了建模效率,适用性广,具有较高的工程实用性。

Description

一种变密度与变构型三维点阵结构建模方法
技术领域:
本发明涉及一种变密度与变构型三维点阵结构建模方法,属于机械结构设计与分析的技术领域。
背景技术:
作为航空、航天、汽车、医疗等多个领域的新型结构类型,类桁架式点阵结构因其轻质、高比刚度、高比强度的机械特性,减震、吸能等功能特征而受到广泛关注与应用。三维点阵结构通常采用先建立三维几何模型,再由几何模型生成有限元模型,但由于几何模型精度问题,经常出现有限元模型存在缺陷甚至无法建立的问题。另一方面,工程实际中对点阵结构的设计自由度需求不断提高,经常同一个结构中存在有不同密度与构型的点阵结构,变密度与变构型点阵结构连接成为有限元模型建立中必须考虑的问题,现有方法中采用刚性联结的方式连接非共节点单元,但该方法额外引入了刚性连接单点,增加了人工操作量与模型计算量,工程中希望实现点阵结构共节点有限元模型的自动建立。
发明内容:
为了实现变密度与变构型三维点阵结构建模方法,本发明的目的是提供一种变密度与变构型三维点阵结构建模方法,通过该方法可以解决几何模型精度不足导致三维点阵结构有限元建模失效甚至失败的问题,同时避免引入额外的刚性连接,实现变密度与变构型的三维点阵结构有限元单元共节点连接。采用批处理方式,减少人工操作量并提高建模效率。
本发明一种变密度与变构型三维点阵结构建模方法,解决其技术问题所采用的技术方案是:
它包括以下步骤:
步骤1:导入结构实体模型;
步骤2:分解结构;通过几何分解方法将结构分为多个子结构。
步骤3:统计结构棱边数E和面数F;
步骤4:对所有子结构Si计算点阵结构布局;通过遍历方法统计子结构数量Nsub,通过人机交互方法确定各子结构点阵结构胞元构型与密度,计算点阵结构杆件分布。
步骤5:建立面内单元;遍历结构的所有面,根据面上的胞元分布,建立面上壳单元。若面上的包括多组胞元分布数据,则只建立单元数量多的一组壳单元。
步骤6:对所有杆件Bj判定是否建立并建立杆上单元;若杆件存在于面内或者杆件之间已经建立单元,则Bj不需要建立单元。对需要建立单元的杆件进行杆上单元建立。
步骤7:单元联结;在变密度与变构型子结构之间联结单元。
步骤8:模型检查;对有限元模型进行检查,检查内容包括重合节点、畸形等。根据检查结果,调整或删除错误节点及单元。
通过以上各步骤,实现了变密度与变构型三维点阵结构共节点有限元模型的快速建模,消除了直接划分网格时几何模型处理精度不足对有限元模型建模的影响。生成的网格为连续的且杆、梁与壳共享节点,从而无需额外引入刚性连接单元。另外,由于采用批处理实现,因此避免了大量人工操作、提高了建模精确度与建模效率。
本发明的有益效果是:提供了变密度与变构型三维点阵结构建模方法,消除了直接划分网格时几何模型处理精度不足对有限元模型建模的影响,避免建模失效或者失败的问题。解决了变密度与变构型三维点阵结构有限元模型之间联结困难的问题,通过生成杆、梁与壳共享节点的有限元单元,避免了额外刚性连接单元的引用,提高了有限元模型的计算精度和效率。本方法实现自动批处理,因此避免了大量人工操作、提高了建模效率。另外,由于三维点阵结构应用越来越广泛,工程中对其设计与分析需求不断提高,本发明在航空、航天、汽车、医疗等领域的机械结构设计中具有很强的工程实用价值。
附图说明
图1.本发明方法流程图。
图2.实施例结构分解为子结构。
图3.杆件重合示意图。
图4.杆件在面内示意图。
图5.实施例一维单元建立结果。
图6.实施例建模结果。
图7.本发明实施例步骤4中基于坐标的三维点阵构型设计方法流程图。
图中符号说明如下:
Bj为杆件;
Pm与Pn为杆件起止点;
(xm,ym,zm)、(xn,yn,zn)为点阵结构杆件Bj三维坐标值;
LAT1、LAT2为点阵胞元结构;
fl为结构面
具体实施方式
下面结合附图和实施例对本发明作详细说明:
本发明一种变密度与变构型三维点阵结构建模方法,见图1所示,它包括以下步骤:
步骤1:现有的有限元软件一般都支持三维模型导入功能,通过该功能在软件中导入结构的三维实体模型。
步骤2:现有的有限元软件一般都支持几何模型编辑功能,通过该功能将结构分为若干子结构。如图2所示,实施例的原有结构分解为4个子结构。
步骤3:根据通过遍历方法计算结构几何体的基本构成元素:棱边以及面的数量,并分别记为棱边数E和面数F。实施例中,棱边数为41,面数为22。
步骤4:设实施例中所有子结构为Si。通过遍历方法统计子结构数量Nsub,通过人机交互方法确定点阵结构胞元构型与密度,计算点阵结构杆件分布。子结构之间的点阵结构胞元构型与密度可以不同,从而形成变密度与变构型的点阵结构。
步骤5:遍历结构的所有面,根据面上的胞元分布,建立面上壳单元。若面上包括多组胞元分布数据,则只建立单元数量多的一组壳单元。
步骤6:判断杆件Bj是否需要建立单元,由于同一位置的杆件可能被不同的点阵胞元共有,则同一位置已经建立的杆件单元不需要重新建立,如图3所示,B1和B2在同一位置,若B1的单元已经建立,则B2不需要重新建立。若杆件起止点均存在于同一个面中,则杆件Bj不需要建立,如图4所示,Bj在面fl内。对需要建立单元的杆件进行杆上单元建立。
步骤7:对相关单元进行单元联结,使得有限元模型具有共节点连续性。其中,所述的相关单元为几何关系上具有共用结点或共用棱边等相关性的单元。
步骤8:对实施例中有限元模型进行检查,检查内容包括重合节点、畸形等。根据检查结果,调整或删除错误节点及单元。实施例建模结果如图5及图6所示。
其中,步骤4所述的特定方法为一种基于坐标的三维点阵构型设计方法,具体如下(流程如图7所示):
步骤40:建立三维点阵构型框架;
步骤41:建立三维点阵构型表达式;
步骤42:获取初始点坐标值Gi(xgi,ygi,zgi)与拓扑属性TGi(i=1,2,3...);
步骤43:计算特征点坐标值Pj(xpj,ypj,zpj)与特征点拓扑属性TPj
步骤44:获取杆件起止点坐标值Bk1(xk1,yk1,zk1),Bk2(xk2,yk2,zk2)和杆拓扑属性TBk
步骤45:输出点阵结构信息
步骤46:若需继续产生点阵结构信息,跳至步骤42,否则结束。

Claims (1)

1.一种变密度与变构型三维点阵结构建模方法,其特征在于:该方法包括以下步骤:
步骤1、导入结构实体模型;
步骤2、分解结构:通过几何分解方法将结构分为多个子结构;
步骤3、统计结构棱边数E和面数F;
步骤4、对所有子结构Si计算点阵结构布局:通过遍历方法统计子结构数量Nsub,通过人机交互方法确定各子结构点阵结构胞元构型与密度,计算点阵结构杆件分布;
步骤5、建立面内单元:遍历结构的所有面,根据面上的胞元分布,建立面上壳单元;若面上的包括多组胞元分布数据,则只建立单元数量多的一组壳单元;
步骤6、对所有杆件Bj判定是否建立并建立杆上单元:若杆件存在于面内或者杆件之间已经建立单元,则Bj不需要建立单元;对需要建立单元的杆件进行杆上单元建立;
步骤7、单元联结:在变密度与变构型子结构之间联结单元;
步骤8、模型检查:对有限元模型进行检查,检查内容包括重合节点、畸形;根据检查结果,调整或删除错误节点及单元。
CN201810593903.3A 2018-06-11 2018-06-11 一种变密度与变构型三维点阵结构建模方法 Active CN109101671B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810593903.3A CN109101671B (zh) 2018-06-11 2018-06-11 一种变密度与变构型三维点阵结构建模方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810593903.3A CN109101671B (zh) 2018-06-11 2018-06-11 一种变密度与变构型三维点阵结构建模方法

Publications (2)

Publication Number Publication Date
CN109101671A true CN109101671A (zh) 2018-12-28
CN109101671B CN109101671B (zh) 2023-01-17

Family

ID=64796800

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810593903.3A Active CN109101671B (zh) 2018-06-11 2018-06-11 一种变密度与变构型三维点阵结构建模方法

Country Status (1)

Country Link
CN (1) CN109101671B (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110570526A (zh) * 2019-09-06 2019-12-13 南京航空航天大学 拓扑信息引导的三维点阵结构件内部几何数据获取方法
CN111310296A (zh) * 2019-10-12 2020-06-19 南京航空航天大学 一种应用于复杂曲面的变密度金属板材点阵芯体设计方法
CN113051805A (zh) * 2021-06-02 2021-06-29 湖南工业职业技术学院 一种基于变密度点阵结构的航空齿轮及其设计方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150190971A1 (en) * 2014-01-09 2015-07-09 Siemens Product Lifecycle Management Software Inc. Method for structure preserving topology optimization of lattice structures for additive manufacturing
CN106557638A (zh) * 2016-12-07 2017-04-05 重庆科技学院 焊接混合双向过渡单元网格模型的建立方法
CN107742014A (zh) * 2017-09-29 2018-02-27 北京空间飞行器总体设计部 基于增材制造的相变储能装置点阵夹层结构的设计方法
CN107742047A (zh) * 2017-11-10 2018-02-27 北京航空航天大学 一种变相对密度octet点阵结构的设计方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150190971A1 (en) * 2014-01-09 2015-07-09 Siemens Product Lifecycle Management Software Inc. Method for structure preserving topology optimization of lattice structures for additive manufacturing
CN106557638A (zh) * 2016-12-07 2017-04-05 重庆科技学院 焊接混合双向过渡单元网格模型的建立方法
CN107742014A (zh) * 2017-09-29 2018-02-27 北京空间飞行器总体设计部 基于增材制造的相变储能装置点阵夹层结构的设计方法
CN107742047A (zh) * 2017-11-10 2018-02-27 北京航空航天大学 一种变相对密度octet点阵结构的设计方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
FUMING ZENG等: ""Research on the Design Method of Lightweight Lattice Sandwich Structure for Additive Manufacturing"", 《INTERNATIONAL CONFERENCE ON MECHANICAL DESIGN, ICMD 2017: ADVANCES IN MECHANICAL DESIGN》, 15 November 2017 (2017-11-15), pages 301 - 316 *
付志方 等: ""节点位置不确定下桁架结构稳健拓扑优化"", 《机械工程学报》, vol. 53, no. 2, 31 January 2017 (2017-01-31), pages 135 - 142 *
阎军 等: ""基于拓展多尺度有限元的点阵材料结构最小柔顺性设计"", 《固体力学学报》, vol. 34, no. 1, 28 February 2013 (2013-02-28), pages 47 - 54 *

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110570526A (zh) * 2019-09-06 2019-12-13 南京航空航天大学 拓扑信息引导的三维点阵结构件内部几何数据获取方法
CN110570526B (zh) * 2019-09-06 2021-02-12 南京航空航天大学 拓扑信息引导的三维点阵结构件内部几何数据获取方法
CN111310296A (zh) * 2019-10-12 2020-06-19 南京航空航天大学 一种应用于复杂曲面的变密度金属板材点阵芯体设计方法
CN111310296B (zh) * 2019-10-12 2024-03-29 南京航空航天大学 一种应用于复杂曲面的变密度金属板材点阵芯体设计方法
CN113051805A (zh) * 2021-06-02 2021-06-29 湖南工业职业技术学院 一种基于变密度点阵结构的航空齿轮及其设计方法
CN113051805B (zh) * 2021-06-02 2021-08-10 湖南工业职业技术学院 一种基于变密度点阵结构的航空齿轮及其设计方法

Also Published As

Publication number Publication date
CN109101671B (zh) 2023-01-17

Similar Documents

Publication Publication Date Title
US11010504B2 (en) Lattice structures for additive manufacturing
Park et al. Unstructured grid adaptation: status, potential impacts, and recommended investments towards CFD 2030
US8831913B2 (en) Method of design optimisation
CN106373184B (zh) 一种三维打印模型摆放所需支撑量快速估算方法
CN109101671A (zh) 一种变密度与变构型三维点阵结构建模方法
CN108647370B (zh) 基于双环迭代的无人直升机气动外形优化设计方法
EP2813958A1 (en) A method and apparatus for preparation of a computational fluid dynamics model
Park Anisotropic output-based adaptation with tetrahedral cut cells for compressible flows
CN104268322B (zh) Weno差分方法的一种边界处理技术
CN107229794B (zh) 一种基于cae和vr的模型构建系统及其管理方法
CN107958103B (zh) 基于妥协决策的零件结构拓扑优化设计方法
CN106971043B (zh) 一种嵌套网格的建立方法及装置
CN112380650A (zh) 工作装置结构件设计方法
Gaitonde et al. A three-dimensional moving mesh method for the calculation of unsteady transonic flows
CN103065015B (zh) 一种基于内力路径几何形态的承载结构低碳节材设计方法
CN110334450A (zh) 一种多块结构网格生成中物面投影错误的修复方法
CN105138750A (zh) 一种可扩展多线交叉机身线框模型构建方法
CN113673186A (zh) 一种基于stl文件的笛卡尔网格快速生成方法
Coirier et al. A Cartesian, cell-based approach for adaptively-refined solutions of the Euler and Navier-Stokes equations
CN110489907A (zh) 一种轨道交通车辆车体数字样机的优化设计方法
JP7257472B2 (ja) 設計支援システム、設計支援方法及び設計支援プログラム
CN115630542A (zh) 一种薄壁加筋结构的加筋布局优化方法
Kumar et al. A multi-GPU based accurate algorithm for simulations of gas-liquid flows
CN101625708A (zh) 可视化结构分析系统
CN104657557A (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