CN110210085A - 内凹六边形负泊松比点阵结构参数化有限元建模方法 - Google Patents
内凹六边形负泊松比点阵结构参数化有限元建模方法 Download PDFInfo
- Publication number
- CN110210085A CN110210085A CN201910417817.1A CN201910417817A CN110210085A CN 110210085 A CN110210085 A CN 110210085A CN 201910417817 A CN201910417817 A CN 201910417817A CN 110210085 A CN110210085 A CN 110210085A
- Authority
- CN
- China
- Prior art keywords
- node
- finite element
- matrix
- cell element
- lattice structure
- 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
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]
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
- Y02T90/00—Enabling technologies or technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Complex Calculations (AREA)
Abstract
本发明公开一种内凹六边形负泊松比点阵结构参数化有限元建模方法,包括如下步骤:(10)独立设计参数获取:获取点阵结构的独立设计参数;(20)有限元模型节点编号矩阵建立:制定参数化建模节点编号规则并建立有限元模型节点编号矩阵;(30)节点及其坐标对应关系矩阵建立:定义坐标原点,建立点阵结构有限元模型的节点‑坐标关系矩阵;(40)单元节点关系矩阵建立:基于梁单元制定参数化建模单元编号规则并建立有限元模型单元及其从属节点间的对应关系矩阵;(50)计算文件编写:将上述步骤编制成脚本程序,并输出Radioss计算文件,实现参数化有限元建模。本发明的有限元建模方法,建模效率高、仿真精度高。
Description
技术领域
本发明属于周期性点阵材料/结构有限元仿真分析领域,特别是一种建模效率高、仿 真精度高的内凹六边形负泊松比点阵结构参数化有限元建模方法。
背景技术
具有负泊松比效应的点阵结构是一种典型的机械超结构,在受到单轴拉伸时发生侧 向膨胀(单轴压缩时侧向收缩)。伴随着独特的拉胀(压缩)现象,负泊松比结构呈现 出多孔轻质、高比强度、抗冲击和隔热等优良性能,不仅能够满足超大型结构的承载、 降噪、冷却等特殊功能性要求,更能在此基础上满足超轻结构的设计理念,即具有高的 比刚度和比强度以及良好的韧性,对航空航天、船舶、汽车以及建筑材料等领域的结构 轻量化和多功能化发展具有重大意义。
负泊松比结构的核心是高效的微观结构设计和尺寸控制,因此其宏观力学性能也强 烈依赖于结构参数的优化设计。内凹六边形负泊松比点阵结构作为一种典型负泊松比结 构具有更优异的承载性能,但其微观结构复杂,设计参数较多,进行优化设计必须依靠有限元手段。
现阶段,在点阵结构优化设计的迭代过程中,需不断调整各独立设计参数值,并对调整后的结构进行性能评估,依靠人工进行逐次调整将使得优化效率低下。
因此在内凹六边形负泊松比点阵结构的优化设计中,需要提出一种高效的参数化有 限元建模方法,以提高建模效率与仿真精度。
发明内容
本发明的目的在于提供一种内凹六边形负泊松比点阵结构参数化有限元建模方法, 建模效率高、仿真精度高。
实现本发明目的的技术方案为:一种内凹六边形负泊松比点阵结构参数化有限元建模方 法,包括如下步骤:
(10)独立设计参数获取:根据周期性胞元的几何特征,获取内凹六边形负泊松比点阵结构的独立设计参数;
(20)有限元模型节点编号矩阵建立:结合点阵结构周期性和对称性特征,制定参数化建模节点编号规则并建立有限元模型节点编号矩阵;
(30)节点及其坐标对应关系矩阵建立:根据微观胞元结构参数,定义坐标原点,按照起始胞元→同一行胞元→同一层胞元→不同层胞元的顺序计算节点的空间坐标,建立点阵结构有限元模型的节点-坐标关系矩阵;
(40)单元节点关系矩阵建立:基于梁单元制定参数化建模单元编号规则并建立有限元模型单元及其从属节点间的对应关系矩阵;
(50)计算文件编写:将上述步骤编制成脚本程序,并输出Radioss计算文件,实 现内凹六边形负泊松比点阵结构的参数化有限元建模。
与现有技术相比,本发明的显著优点为:建模效率高、仿真精度高。
本发明结合点阵结构周期性和对称性条件,采用参数化方法建立结构的单元-节点- 坐标信息矩阵,解决了具有不同微观胞元结构参数和宏观周期性参数的内凹六边形负泊 松比点阵结构的参数化建模问题,具有快速、准确和适应性好的优点,提高了基于有限元分析的结构优化设计效率。
附图说明
图1为本发明内凹六边形负泊松比点阵结构参数化有限元建模方法的主流程图。
图2为内凹六边形负泊松比点阵结构示意图。
图3为图1中有限元模型节点编号矩阵建立步骤的流程图。
图4为点阵结构节点编号规则示意图。
其中,图4a为层间胞元编号规则,图4b为基础层胞元分类编码,图4c为起始胞 元节点编号规则示意图。
图5为图1中单元节点关系矩阵建立步骤的流程图。
图6为点阵结构中起始胞元单元编号规则示意图。
具体实施方式
如图1所示,本发明内凹六边形负泊松比点阵结构参数化有限元建模方法,包括如下步骤:
(10)独立设计参数获取:根据周期性胞元的几何特征,获取内凹六边形负泊松比点阵结构的独立设计参数;
所述(10)独立设计参数获取步骤中,所述内凹六边形负泊松比点阵结构是由多个内凹六边形胞元在三维方向阵列排布连接而成;如图2所示。
所述点阵结构的独立设计参数主要包括:伞状胞臂长度L、垂直胞壁长度H、胞臂截面系数α和胞元夹角θ,胞元在正交三个方向上的数量分别是n_X×n_Y×n_Z;所 述独立设计参数填写在参数化建模脚本程序的开头部分。
(20)有限元模型节点编号矩阵建立:结合点阵结构周期性和对称性特征,制定参数化建模节点编号规则并建立有限元模型节点编号矩阵;
如图3所示,所述(20)有限元模型节点编号矩阵建立步骤包括:
(21)点阵结构拆分:将点阵结构沿Y向拆分为n_Y个胞元层,每层胞元节点编 号通过基础层对应节点累加而来,所述基础层胞元根据X、Z向周期性特点,区分为垂 直胞臂和伞状胞臂两大类;
(22)大类编码:先对垂直胞臂和伞状胞臂依次进行大类编码;
(23)节点编号:选取起始垂直胞臂和伞状胞臂进行节点编号;
(24)编号递推:依据大类编码进行胞元之间的编号递推,形成有限元模型节点编号矩阵。
如图4所示,在对节点进行编号时将点阵结构沿Y向拆分为由n_Y个胞元层组成,每层胞元节点编号通过基础层对应节点累加而来;基础层胞元根据X、Z向周期性特点 区分为垂直胞臂和伞状胞臂两大类;先对垂直胞臂和伞状胞臂依次进行大类编码,再选 取起始垂直胞臂和伞状胞臂进行节点编号,而后依据大类编码进行胞元之间的编号递 推,形成有限元模型节点编号矩阵;
(30)节点及其坐标对应关系矩阵建立:根据微观胞元结构参数,定义坐标原点,按照起始胞元→同一行胞元→同一层胞元→不同层胞元的顺序计算节点的空间坐标,建立点阵结构有限元模型的节点-坐标关系矩阵;
此处以起始胞元中#1垂直胞臂为例描述节点-坐标信息矩阵建立方式,节点在图示 笛卡尔坐标系下的坐标可由下列方法求得:图4-c中#1垂直胞臂(基础层初始胞臂)的节点拥有相同的x、z坐标,垂直胞臂的单元数固定为6个,将节点1定为坐标原点, 则#1胞臂的1~7号节点的坐标可用矩阵N#1表示:
其中,行号代表节点的节点编号,矩阵每一行三个元素分别代表x、y、z坐标。图 4-b中与#1胞臂处于同一行的#2、#3、#4胞臂可视为#1胞臂沿Z轴阵列而来,阵列总 数即为内芯在Z轴方向的胞元数量nZ,因此第一行垂直胞臂的节点坐标矩阵Nrow#1可表 示为:
其中,R称为Z向间距矩阵,其第三列元素均为2L sinθ,其余元素全为0,矩阵维 度与N#1一致。其余垂直胞臂行均可由第一行沿X轴阵列而来,采用相同方法计算基础 层胞元中伞状胞臂的节点坐标,此处不具体展开。
(40)单元节点关系矩阵建立:基于梁单元制定参数化建模单元编号规则并建立有限元模型单元及其从属节点间的对应关系矩阵;
如图5所示,所述(40)单元节点关系矩阵建立步骤包括:
(41)构建梁单元:在生成所有节点基础上,在相邻两节点之间构建Euler-Bernoulli 梁单元;
(42)单元编号:按照起始胞元→同一行胞元→同一层胞元→不同层胞元的顺序,对单元进行编号;
(43)编号匹配:将构成单元的从属节点编号与对应编号的单元匹配,建立点阵结构有限元模型的单元-节点关系矩阵。
基于梁单元制定参数化建模单元编号规则并建立有限元模型单元及其从属节点间 的对应关系矩阵;生成所有节点基础上,在相邻两节点之间构建Euler-Bernoulli梁单元; 根据图6所示起始胞元单元编号规则,按照起始胞元→同一行胞元→同一层胞元→不同 层胞元的顺序对单元进行编号,将构成单元的从属节点编号与单元匹配,建立点阵结构有限元模型的单元-节点矩阵;
此处以起始单元为例描述单元-节点信息矩阵建立方式,如图4所示,#1垂直胞臂中由相邻节点构成的1~6号梁单元节点矩阵E#1可表示为:
其中,行号代表单元的单元编号,矩阵中每一行的两个元素表示构成该单元的节点 的编号。其余垂直胞臂的单元节点矩阵均可由E#1转化而来,因此基础层胞元所有垂直胞臂的单元节点矩阵Elayer#1可表示为:
其中,I是维度与E#1一致的矩阵且矩阵元素全为1。
采用相同方法计算基础层胞元中伞状胞臂的单元节点信息,需要注意的是,在对节 点进行编号时曾将垂直胞臂和伞状胞臂分为两个独立的节点群,而在构建梁单元时会出 现伞状胞臂中的部分单元需要同时用两个节点群的节点进行表征,如#1伞状胞臂中的4号、8号、12号等单元,因此在构建单元节点矩阵时将不同节点群的节点编号进行合并 描述,此处不具体展开。
(50)计算文件编写:将上述步骤编制成脚本程序,并输出Radioss计算文件,实 现内凹六边形负泊松比点阵结构的参数化有限元建模。
所述(50)计算文件编写步骤具体为:
采用Tcl/Tk语言二次开发,模型文件中按照Radioss文件格式写入单元节点编号信 息、单元-节点信息、节点-坐标信息,进一步写入材料、属性信息以及约束、载荷信息,选择求解控制卡片和输出结果类型,完成计算文件编写后可直接提交Radioss求解器进 行数值计算。
本发明解决了具有不同微观胞元结构参数和宏观周期性参数的内凹六边形负泊松比点 阵结构的参数化建模问题,具有快速、准确和适应性好的优点,提高了基于有限元分析 的结构优化设计效率。
Claims (5)
1.一种内凹六边形负泊松比点阵结构参数化有限元建模方法,其特征在于,包括如下步骤:
(10)独立设计参数获取:根据周期性胞元的几何特征,获取内凹六边形负泊松比点阵结构的独立设计参数;
(20)有限元模型节点编号矩阵建立:结合点阵结构周期性和对称性特征,制定参数化建模节点编号规则并建立有限元模型节点编号矩阵;
(30)节点及其坐标对应关系矩阵建立:根据微观胞元结构参数,定义坐标原点,按照起始胞元→同一行胞元→同一层胞元→不同层胞元的顺序计算节点的空间坐标,建立点阵结构有限元模型的节点-坐标关系矩阵;
(40)单元节点关系矩阵建立:基于梁单元制定参数化建模单元编号规则并建立有限元模型单元及其从属节点间的对应关系矩阵;
(50)计算文件编写:将上述步骤编制成脚本程序,并输出Radioss计算文件,实现内凹六边形负泊松比点阵结构的参数化有限元建模。
2.根据权利要求1所述的有限元建模方法,其特征在于:
所述(10)独立设计参数获取步骤中,所述内凹六边形负泊松比点阵结构是由多个内凹六边形胞元在三维方向阵列排布连接而成;所述点阵结构的独立设计参数主要包括:伞状胞臂长度L、垂直胞壁长度H、胞臂截面系数α和胞元夹角θ,胞元在正交三个方向上的数量分别是nx×nY×nz;所述独立设计参数填写在参数化建模脚本程序的开头部分。
3.根据权利要求1所述的有限元建模方法,其特征在于,所述(20)有限元模型节点编号矩阵建立步骤包括:
(21)点阵结构拆分:将点阵结构沿Y向拆分为nY个胞元层,每层胞元节点编号通过基础层对应节点累加而来,所述基础层胞元根据X、Z向周期性特点,区分为垂直胞臂和伞状胞臂两大类;
(22)大类编码:先对垂直胞臂和伞状胞臂依次进行大类编码;
(23)节点编号:选取起始垂直胞臂和伞状胞臂进行节点编号;
(24)编号递推:依据大类编码进行胞元之间的编号递推,形成有限元模型节点编号矩阵。
4.根据权利要求1所述的有限元建模方法,其特征在于,所述(40)单元节点关系矩阵建立步骤包括:
(41)构建梁单元:在生成所有节点基础上,在相邻两节点之间构建Euler-Bernoulli梁单元;
(42)单元编号:按照起始胞元→同一行胞元→同一层胞元→不同层胞元的顺序,对单元进行编号;
(43)编号匹配:将构成单元的从属节点编号与对应编号的单元匹配,建立点阵结构有限元模型的单元-节点关系矩阵。
5.根据权利要求1所述的有限元建模方法,其特征在于,所述(50)计算文件编写步骤具体为:
采用Tcl/Tk语言二次开发,模型文件中按照Radioss文件格式写入单元节点编号信息、单元-节点信息、节点-坐标信息,进一步写入材料、属性信息以及约束、载荷信息,选择求解控制卡片和输出结果类型,完成计算文件编写后可直接提交Radioss求解器进行数值计算。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910417817.1A CN110210085B (zh) | 2019-05-20 | 2019-05-20 | 内凹六边形负泊松比点阵结构参数化有限元建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910417817.1A CN110210085B (zh) | 2019-05-20 | 2019-05-20 | 内凹六边形负泊松比点阵结构参数化有限元建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110210085A true CN110210085A (zh) | 2019-09-06 |
CN110210085B CN110210085B (zh) | 2022-09-06 |
Family
ID=67787797
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910417817.1A Active CN110210085B (zh) | 2019-05-20 | 2019-05-20 | 内凹六边形负泊松比点阵结构参数化有限元建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110210085B (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110598319A (zh) * | 2019-09-11 | 2019-12-20 | 杭州群核信息技术有限公司 | 一种基于可变参数的空间划分的规则语法 |
CN112380741A (zh) * | 2020-10-31 | 2021-02-19 | 华南理工大学 | 一种面向随机纤维网格的有限元模型建模方法 |
CN113190905A (zh) * | 2021-05-13 | 2021-07-30 | 中国电子工程设计院有限公司 | 一种建筑模型的分析方法、装置及存储介质 |
CN115427172A (zh) * | 2020-04-15 | 2022-12-02 | 西门子能源全球有限两合公司 | 用在增材制造应用中的拉胀三维结构 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107577875A (zh) * | 2017-09-07 | 2018-01-12 | 东南大学 | 一种缝合式夹芯板有限元参数化建模方法 |
CN109766656A (zh) * | 2019-01-25 | 2019-05-17 | 北京航空航天大学 | 一种基于拓扑优化的梯度点阵结构设计方法 |
-
2019
- 2019-05-20 CN CN201910417817.1A patent/CN110210085B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107577875A (zh) * | 2017-09-07 | 2018-01-12 | 东南大学 | 一种缝合式夹芯板有限元参数化建模方法 |
CN109766656A (zh) * | 2019-01-25 | 2019-05-17 | 北京航空航天大学 | 一种基于拓扑优化的梯度点阵结构设计方法 |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110598319A (zh) * | 2019-09-11 | 2019-12-20 | 杭州群核信息技术有限公司 | 一种基于可变参数的空间划分的规则语法 |
CN115427172A (zh) * | 2020-04-15 | 2022-12-02 | 西门子能源全球有限两合公司 | 用在增材制造应用中的拉胀三维结构 |
CN112380741A (zh) * | 2020-10-31 | 2021-02-19 | 华南理工大学 | 一种面向随机纤维网格的有限元模型建模方法 |
CN113190905A (zh) * | 2021-05-13 | 2021-07-30 | 中国电子工程设计院有限公司 | 一种建筑模型的分析方法、装置及存储介质 |
CN113190905B (zh) * | 2021-05-13 | 2024-05-10 | 中国电子工程设计院有限公司 | 一种建筑模型的分析方法、装置及存储介质 |
Also Published As
Publication number | Publication date |
---|---|
CN110210085B (zh) | 2022-09-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN110210085A (zh) | 内凹六边形负泊松比点阵结构参数化有限元建模方法 | |
CN107844676B (zh) | 一种基于多性能约束的结构拓扑优化设计方法 | |
Lebaal et al. | Optimised lattice structure configuration for additive manufacturing | |
CN103034766B (zh) | 一种确定层合板复合材料的铺层角方向与厚度的方法 | |
CN109145393B (zh) | 一种结构仿生轻量化设计方法 | |
CN107301287A (zh) | 自动提取车身接附点动刚度的方法 | |
Hadi et al. | Design configurations and creation of lattice structures for metallic additive manufacturing | |
CN108763827B (zh) | 一种输电塔有限元模型建立方法及装置 | |
CN106202841B (zh) | 一种机床床身三点支撑位置的优化设计方法 | |
CN110852000B (zh) | 一种车身结构优化方法 | |
CN109408939A (zh) | 一种兼顾应力和位移约束的薄板结构加强筋分布优化的改进方法 | |
CN110795835A (zh) | 一种基于自动同步建模的三维工序模型逆向生成方法 | |
CN108170947B (zh) | 一种基于萤火虫算法获取新型点阵结构的方法 | |
CN103065015B (zh) | 一种基于内力路径几何形态的承载结构低碳节材设计方法 | |
Rosso et al. | Optimization approaches in design for additive manufacturing | |
CN111539064B (zh) | 一种中线及厚度提取方法及有限元网格自动生成方法 | |
CN110517350B (zh) | 一种蒙特卡洛粒子输运计算程序的预处理方法及处理器 | |
CN116118196A (zh) | 一种基于力流管载荷路径的连续纤维3d打印路径设计方法 | |
CN114840926B (zh) | 一种轮胎复杂花纹三维有限元网格的生成方法 | |
CN113326582B (zh) | 一种基于应力分布的变密度点阵结构及其设计方法 | |
CN113111612B (zh) | 一种基于自适应空间剖分的离散点云重复点快速查找方法 | |
CN113066523B (zh) | 基于空间三角函数的鳞翅目微纳结构统一表征方法和系统 | |
CN115221622A (zh) | 一种大尺寸复合材料机身壁板装配定位夹持布局优化方法 | |
CN113343512A (zh) | 一种基于Mobile-U-Net的多尺度拓扑优化设计方法 | |
Jin et al. | Optimal design and modeling of 3D variable-density lattice structures |
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 |