CN108170947A - 一种基于萤火虫算法获取新型点阵结构的方法 - Google Patents

一种基于萤火虫算法获取新型点阵结构的方法 Download PDF

Info

Publication number
CN108170947A
CN108170947A CN201711439036.XA CN201711439036A CN108170947A CN 108170947 A CN108170947 A CN 108170947A CN 201711439036 A CN201711439036 A CN 201711439036A CN 108170947 A CN108170947 A CN 108170947A
Authority
CN
China
Prior art keywords
lattice structure
swarm algorithm
glowworm swarm
rod piece
obtained based
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
CN201711439036.XA
Other languages
English (en)
Other versions
CN108170947B (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.)
Chongqing University
Original Assignee
Chongqing 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 Chongqing University filed Critical Chongqing University
Priority to CN201711439036.XA priority Critical patent/CN108170947B/zh
Publication of CN108170947A publication Critical patent/CN108170947A/zh
Application granted granted Critical
Publication of CN108170947B publication Critical patent/CN108170947B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N3/00—Computing arrangements based on biological models
    • G06N3/004—Artificial life, i.e. computing arrangements simulating life
    • 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/06—Power analysis or power optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Theoretical Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • General Health & Medical Sciences (AREA)
  • Artificial Intelligence (AREA)
  • Data Mining & Analysis (AREA)
  • Biophysics (AREA)
  • Biomedical Technology (AREA)
  • Molecular Biology (AREA)
  • Computing Systems (AREA)
  • Computational Linguistics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mathematical Physics (AREA)
  • Software Systems (AREA)
  • Health & Medical Sciences (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明公开了一种基于萤火虫算法获取新型点阵结构的方法，包括以下步骤：s1:建立点阵结构的基结构；s2:建立以所述基结构中杆单元的截面面积为设计变量、以杆件的轴向力作为约束函数、以最小体积作为评价函数的数学模型，采用萤火虫算法获得在一定约束条件下的最优点阵结构；本发明不采用模拟分子构型制造多孔材料，而是采用基结构法作为杆系结构的拓扑生成方法与萤火虫算法结合完成结构拓扑优化，获得轻质高强的点阵材料。

Description

一种基于萤火虫算法获取新型点阵结构的方法
技术领域
本发明涉及新型点阵结构的方法，具体地，涉及一种基于萤火虫算法获取新型点阵结构的方法。
背景技术
点阵结构的构型设计是研究点阵结构最基础的工作，它将直接影响后面选用怎样的制备方法和具有怎样的力学性能和特殊功能。目前点阵结构的单胞结构是模拟分子点阵构型制造出的，但由于不同的分子结构具有不同的优缺点，力学性能不能达到最好。而且在母本材料与加工参数相同的情况下，对于不同的微结构构型的点阵材料具有不同的力学性能。
发明内容
有鉴于此，本发明提供一种不采用模拟分子构型制造多孔材料，而是采用基结构法作为杆系结构的拓扑生成方法与萤火虫算法结合完成结构拓扑优化，获得轻质高强的点阵材料。
本发明的基于萤火虫算法获取新型点阵结构的方法，包括以下步骤：s1:建立点阵结构的基结构；s2:建立以所述基结构中杆单元的截面面积为设计变量、以杆件的轴向力作为约束函数、以最小体积作为评价函数的数学模型，采用萤火虫算法获得在一定约束条件下的最优点阵结构；
进一步，步骤s1中，将所有可能节点进行两两连接形成初始基础结构，并在所述初始结构的基础上，优化算法选择每根杆件的最小截面或删除其中不必要的杆件及节点为基结构；
进一步，步骤s2中,通过所述数学模型获得最优点阵结构的步骤包括：
s21.对所述模型进行参数设置；
s22.赋予杆件界面的初始值，生成初始变量；
s23.计算目标函数；
s24.判断总体刚度矩阵是否正定以及杆件轴向力约束与否，确定出惩罚因子；
S25.进行个体评价；
s26.更新各杆的截面积；
S27.重复步骤s23-s25直至达到终止条件；
S28.输出各杆的最终截面积；
进一步，步骤s23中，所述目标函数为：
式中：li为第i个杆单元的杆长，xi为第i个杆单元的截面面积；
进一步，步骤s24中，当所述总体刚度矩阵正定、且轴向力约束时，惩罚因子为1，否则，惩罚因子为1000。
本发明的有益效果：
本发明采用基结构法作为杆系结构的拓扑生成方法与萤火虫算法结合完成结构拓扑优化，获得轻质高强的点阵材料。
附图说明
下面结合附图和实施例对本发明作进一步描述。
图1为本实施例中建立的点阵结构的基结构的模型；
图2为本实施例中通过数学模型获得最优点阵结构的流程图；
图3为本实施例中优化得到的最优点阵结构的单胞模型；
图4为本实施例中优化出的新型结构与BCC仿真中加载的位移与承载力的变化曲线；
图5为本实施例中优化出的新型结构与BCC结构的吸能对比图。
具体实施方式
本实施例的基于萤火虫算法获取新型点阵结构的方法，包括以下步骤，
s1:如图1所示，建立一个由正方体8个顶点Di(i＝1,2……,8)，六个的面心点Si(i＝1,2……,6)与1个体心B点组成的节点集合，任意两个节点之间用离散杆单元相联的基结构。由于在多孔泡沫结构中，当含有贯穿整个模型的直杆时，模型的弹性模量与相对密度成线性关系，使纵向的拉伸与压缩成为唯一的变形形式，故，大多数多孔泡沫不包含直杆，令梁的弯曲变形发挥作用。本文选取的基结构中不考虑横杆与竖杆。在基结构的基础上，优化算法选择每根杆件的最小截面或删除其中不必要的杆件及节点为基结构法；
s2:通过基结构法作为杆系结构拓扑生成方法与萤火虫算法结合。建立以基结构中杆单元的截面面积为设计变量、以杆件的轴向力作为约束函数、以最小体积作为评价函数的数学模型，采用萤火虫算法获得在一定约束条件下的最优点阵结构。进行基结构优化的流程图如图2所示。其步骤包括：
s21.对所述模型进行参数设置；参数包括潜在解的个数、最大迭代次数M、步长因子α、最大吸引力β、光吸收系数γ；
s22.赋予杆件界面的初始值，生成初始变量。其中设计变量为：
式中：xi为第i个杆单元的截面面积。
s23.计算目标函数：
式中：li为第i个杆单元的杆长；
s24.判断总体刚度矩阵是否正定以及杆件轴向力是否满足约束条件，确定出惩罚因子；
S25.进行个体评价：
|FN|≤xi[σi](i＝1,2,…,32)
式中：|FN|为杆单元的轴向力；xi为第i个杆单元的截面面积；[σi]为第i根杆的许用应力。
s26.更新各杆的截面积；
S27.重复步骤s23-s25直至达到终止条件；
S28.输出各杆的最终截面积。
在步骤s24中，当所述总体刚度矩阵正定、且满足轴向力约束条件时，惩罚因子为1，否则，惩罚因子为1000。
优化后得到的新型模型如图3所示；
使用ABAQUS仿真软件对新型结构与移除模型即经典模型BCC结构进行准静态压缩实验仿真。仿真结果如图4与图5所示，
以移除模型即经典模型BCC点阵结构为参考分析对象，两种模型在相同约束条件与相同的相对密度情况下，新型结构所能承受的最大载荷大5000N，弹性模量增加了3111MPa，在吸能上，新型结构也优于BCC点阵结构。在相同相对密度条件下，新型结构显示出力学性能与吸能特性上的明显优势，为点阵材料构型的设计提供了理论参考。
最后说明的是，以上实施例仅用以说明本发明的技术方案而非限制，尽管参照较佳实施例对本发明进行了详细说明，本领域的普通技术人员应当理解，可以对本发明的技术方案进行修改或者等同替换，而不脱离本发明技术方案的宗旨和范围，其均应涵盖在本发明的权利要求范围当中。

Claims (5)

1.一种基于萤火虫算法获取新型点阵结构的方法，其特征在于，包括以下步骤：
s1:建立点阵结构的基结构；
s2:建立以所述基结构中杆单元的截面面积为设计变量、以杆件的轴向力作为约束函数、以最小体积作为评价函数的数学模型，采用萤火虫算法获得在一定约束条件下的最优点阵结构。
2.根据权利要求1所述的基于萤火虫算法获取新型点阵结构的方法，其特征在于：步骤s1中，将所有可能节点进行两两连接形成初始基础结构，并在所述初始结构的基础上，优化算法选择每根杆件的最小截面或删除其中不必要的杆件及节点为基结构。
3.根据权利要求2所述的基于萤火虫算法获取新型点阵结构的方法，其特征在于：步骤s2中,通过所述数学模型获得最优点阵结构的步骤包括：
s21.对所述模型进行参数设置；
s22.赋予杆件界面的初始值，生成初始变量；
s23.计算目标函数；
s24.判断总体刚度矩阵是否正定以及杆件轴向力是否满足约束条件，确定出惩罚因子；
S25.进行个体评价；
s26.更新各杆的截面积；
S27.重复步骤s23-s25直至达到终止条件；
S28.输出各杆的最终截面积。
4.根据权利要求3所述的基于萤火虫算法获取新型点阵结构的方法，其特征在于：步骤s23中，所述目标函数为：
式中：li为第i个杆单元的杆长，xi为第i个杆单元的截面面积。
5.根据权利要求4所述的基于萤火虫算法获取新型点阵结构的方法，其特征在于：步骤s24中，当所述总体刚度矩阵正定、且轴向力满足约束条件时，惩罚因子为1，否则，惩罚因子为1000。
CN201711439036.XA 2017-12-27 2017-12-27 一种基于萤火虫算法获取新型点阵结构的方法 Active CN108170947B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711439036.XA CN108170947B (zh) 2017-12-27 2017-12-27 一种基于萤火虫算法获取新型点阵结构的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711439036.XA CN108170947B (zh) 2017-12-27 2017-12-27 一种基于萤火虫算法获取新型点阵结构的方法

Publications (2)

Publication Number Publication Date
CN108170947A true CN108170947A (zh) 2018-06-15
CN108170947B CN108170947B (zh) 2021-03-26

Family

ID=62521851

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711439036.XA Active CN108170947B (zh) 2017-12-27 2017-12-27 一种基于萤火虫算法获取新型点阵结构的方法

Country Status (1)

Country Link
CN (1) CN108170947B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109829261A (zh) * 2019-04-03 2019-05-31 大连理工大学 一种由变截面圆管组成的轻质点阵结构及其优化设计方法
CN111639399A (zh) * 2020-05-08 2020-09-08 燕山大学 考虑吸能刚度的点阵筛选和相对强度点阵填充方法
CN112800553A (zh) * 2021-01-04 2021-05-14 中国科学院力学研究所 多级可控渐进吸能点阵结构
CN112861252A (zh) * 2020-12-24 2021-05-28 中国航空工业集团公司成都飞机设计研究所 一种自定义点阵标准单元及点阵结构

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101195197A (zh) * 2007-12-18 2008-06-11 西安交通大学 一种蜂窝环形点阵的成型方法
ES2345521B1 (es) * 2008-03-13 2011-09-16 Grande Grupo Andaluz De Estudios, S.L. Camara de combustion.
CN106940148A (zh) * 2016-11-26 2017-07-11 西南电子技术研究所（中国电子科技集团公司第十研究所） 变梯度分形点阵夹芯强化相变热沉
CN107016196A (zh) * 2017-04-10 2017-08-04 合肥工业大学 一种基于萤火虫算法的窄基角钢输电塔结构优化设计方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101195197A (zh) * 2007-12-18 2008-06-11 西安交通大学 一种蜂窝环形点阵的成型方法
ES2345521B1 (es) * 2008-03-13 2011-09-16 Grande Grupo Andaluz De Estudios, S.L. Camara de combustion.
CN106940148A (zh) * 2016-11-26 2017-07-11 西南电子技术研究所（中国电子科技集团公司第十研究所） 变梯度分形点阵夹芯强化相变热沉
CN107016196A (zh) * 2017-04-10 2017-08-04 合肥工业大学 一种基于萤火虫算法的窄基角钢输电塔结构优化设计方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
仲梁维等: "轻质点阵机构的参数化建模及力学性能研究", 《中国机械工程》 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109829261A (zh) * 2019-04-03 2019-05-31 大连理工大学 一种由变截面圆管组成的轻质点阵结构及其优化设计方法
CN111639399A (zh) * 2020-05-08 2020-09-08 燕山大学 考虑吸能刚度的点阵筛选和相对强度点阵填充方法
CN111639399B (zh) * 2020-05-08 2022-03-22 燕山大学 考虑吸能刚度的点阵筛选和相对强度点阵填充方法
CN112861252A (zh) * 2020-12-24 2021-05-28 中国航空工业集团公司成都飞机设计研究所 一种自定义点阵标准单元及点阵结构
CN112861252B (zh) * 2020-12-24 2022-07-12 中国航空工业集团公司成都飞机设计研究所 用于飞行器结构设计和建模的点阵标准单元及点阵结构
CN112800553A (zh) * 2021-01-04 2021-05-14 中国科学院力学研究所 多级可控渐进吸能点阵结构
CN112800553B (zh) * 2021-01-04 2024-01-30 中国科学院力学研究所 多级可控渐进吸能点阵结构

Also Published As

Publication number Publication date
CN108170947B (zh) 2021-03-26

Similar Documents

Publication Publication Date Title
CN108170947A (zh) 一种基于萤火虫算法获取新型点阵结构的方法
CN107103138B (zh) 一种激光喷丸变刚度轻量化方法
CN110147599B (zh) 一种基于遗传算法的索杆张力结构质量优化方法和系统
Hasançebi et al. Layout optimisation of trusses using simulated annealing
CN108038318A (zh) 一种变截面金属点阵结构初始刚度及塑性破坏强度计算算法
CN102867101B (zh) 一种确定桁架结构参数的方法
Çarbaş et al. Optimum topology design of various geometrically nonlinear latticed domes using improved harmony search method
CN105243460A (zh) 一种基于改进人工鱼群算法的输电塔塔腿辅材拓扑结构优化方法
CN107016196A (zh) 一种基于萤火虫算法的窄基角钢输电塔结构优化设计方法
CN111353246B (zh) 一种混凝土构件设计的静动力多目标拓扑演化方法
CN107742047A (zh) 一种变相对密度octet点阵结构的设计方法
Talatahari et al. Improved bat algorithm for optimum design of large-scale truss structures
CN114999591A (zh) 一种多构型点阵结构的拓扑优化方法
CN110188471A (zh) 一种考虑风荷载的铝合金屋架优化方法
CN103065015A (zh) 一种基于内力路径几何形态的承载结构低碳节材设计方法
CN103150444B (zh) 一种基于两层面承载安全性的工程结构设计与优化方法
CN112883472A (zh) 一种基于数值逆吊法的树状结构拓扑优化设计方法
CN109255142A (zh) 基于小生境遗传算法的环形张拉整体结构拓扑优化方法
Hunt et al. A design tool for spatial tree structures
CN111274624A (zh) 一种基于rbf代理模型的多工况异形节点拓扑优化设计方法
CN111539138B (zh) 基于阶跃函数的结构动力学峰值时域响应灵敏度求解方法
Shallan et al. Optimization of plane and space trusses using genetic algorithms
CN111523162A (zh) 基于两种极限状态的深梁分离式配筋多目标拓扑设计方法
CN109543226A (zh) 一种用于空间结构的混合智能优化方法
Baradaran et al. Application of an improved genetic algorithm for optimal design of planar steel frames

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