CN109101671B - 一种变密度与变构型三维点阵结构建模方法 - Google Patents
一种变密度与变构型三维点阵结构建模方法 Download PDFInfo
- Publication number
- CN109101671B CN109101671B CN201810593903.3A CN201810593903A CN109101671B CN 109101671 B CN109101671 B CN 109101671B CN 201810593903 A CN201810593903 A CN 201810593903A CN 109101671 B CN109101671 B CN 109101671B
- Authority
- CN
- China
- Prior art keywords
- lattice structure
- unit
- variable
- model
- configuration
- 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.)
- Active
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/17—Mechanical parametric or variational design
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- General Physics & Mathematics (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
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、模型检查:对有限元模型进行检查,检查内容包括重合节点、畸形;根据检查结果,调整或删除错误节点及单元。
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 CN109101671A (zh) | 2018-12-28 |
CN109101671B true 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) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110570526B (zh) * | 2019-09-06 | 2021-02-12 | 南京航空航天大学 | 拓扑信息引导的三维点阵结构件内部几何数据获取方法 |
CN111310296B (zh) * | 2019-10-12 | 2024-03-29 | 南京航空航天大学 | 一种应用于复杂曲面的变密度金属板材点阵芯体设计方法 |
CN113051805B (zh) * | 2021-06-02 | 2021-08-10 | 湖南工业职业技术学院 | 一种基于变密度点阵结构的航空齿轮及其设计方法 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106557638A (zh) * | 2016-12-07 | 2017-04-05 | 重庆科技学院 | 焊接混合双向过渡单元网格模型的建立方法 |
CN107742047A (zh) * | 2017-11-10 | 2018-02-27 | 北京航空航天大学 | 一种变相对密度octet点阵结构的设计方法 |
CN107742014A (zh) * | 2017-09-29 | 2018-02-27 | 北京空间飞行器总体设计部 | 基于增材制造的相变储能装置点阵夹层结构的设计方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9902114B2 (en) * | 2014-01-09 | 2018-02-27 | Siemens Product Lifecycle Management Software Inc. | Method for creating three dimensional lattice structures in computer-aided design models for additive manufacturing |
-
2018
- 2018-06-11 CN CN201810593903.3A patent/CN109101671B/zh active Active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
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)
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》.2017,第301-316页. * |
付志方 等."节点位置不确定下桁架结构稳健拓扑优化".《机械工程学报》.2017,第53卷(第2期),第135-142页. * |
阎军 等."基于拓展多尺度有限元的点阵材料结构最小柔顺性设计".《固体力学学报》.2013,第34卷(第1期),第47-54页. * |
Also Published As
Publication number | Publication date |
---|---|
CN109101671A (zh) | 2018-12-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11429765B2 (en) | Meshless method for solid mechanics simulation, electronic device, and storage medium | |
CN109101671B (zh) | 一种变密度与变构型三维点阵结构建模方法 | |
Park et al. | Unstructured grid adaptation: status, potential impacts, and recommended investments towards CFD 2030 | |
Secco et al. | Component-based geometry manipulation for aerodynamic shape optimization with overset meshes | |
CN102063555B (zh) | 基于网格结构的有限元数值模型查错方法 | |
CN108351906B (zh) | 用于具有晶格结构的部件的建模的系统和方法 | |
Arora et al. | Volumetric michell trusses for parametric design & fabrication | |
US20190325098A1 (en) | System, method, and computer program for part model generation and analysis and part production and validation | |
CN116629079B (zh) | 混合有限元空间构造及求解线弹性力学问题的方法及装置 | |
CN117150858B (zh) | 含裂纹连续纤维增强树脂基复合材料三维有限元建模方法 | |
Suchde | Conservation and accuracy in meshfree generalized finite difference methods | |
CN112231863B (zh) | 太阳翼电池阵基板建模方法、装置、设备及存储介质 | |
CN114022641A (zh) | 一种含凹坑缺陷蜂窝夹芯板的有限元网格生成方法 | |
Bournival et al. | A mesh-geometry based method for coupling 1D and 3D elements | |
Wu et al. | A global approach to multi-axis swept mesh generation | |
CN105138750A (zh) | 一种可扩展多线交叉机身线框模型构建方法 | |
Rendall et al. | Conservative unsteady aerodynamic simulation of arbitrary boundary motion using structured and unstructured meshes in time | |
Liu et al. | Small polyhedron reconnection for mesh improvement and its implementation based on advancing front technique | |
Brezillon et al. | Development and application of a flexible and efficient environment for aerodynamic shape optimisation | |
Jain et al. | PostBL: Post-mesh boundary layer mesh generation tool | |
Mei et al. | IGA based bi-layer fiber angle optimization method for variable stiffness composites | |
CN113806951A (zh) | 一种基于半边数据结构的自然邻近点搜索的弹性仿真方法 | |
Montoya-Zapata et al. | Sensitivity Analysis in Shape Optimization using Voxel Density Penalization. | |
Georgiou | Interactive Structural Analysis and Form-Finding | |
Wang et al. | Meshfree simulation of flow around airfoil using different turbulent models |
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 |