CN109145408B - 一种基于改进粒子群算法的隐式曲面多孔结构优化方法 - Google Patents
一种基于改进粒子群算法的隐式曲面多孔结构优化方法 Download PDFInfo
- Publication number
- CN109145408B CN109145408B CN201810865521.1A CN201810865521A CN109145408B CN 109145408 B CN109145408 B CN 109145408B CN 201810865521 A CN201810865521 A CN 201810865521A CN 109145408 B CN109145408 B CN 109145408B
- Authority
- CN
- China
- Prior art keywords
- particle
- curved surface
- porous structure
- optimal
- particle swarm
- 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/10—Geometric CAD
- G06F30/17—Mechanical parametric or variational design
-
- 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
- G06N3/006—Artificial life, i.e. computing arrangements simulating life based on simulated virtual individual or collective life forms, e.g. social simulations or particle swarm optimisation [PSO]
Abstract
本发明公开了一种基于改进粒子群算法的隐式曲面多孔结构优化方法,包括输入待优化隐式曲面函数表达式,多孔实体壁厚,粒子群规模,最大迭代次数;初始化粒子,每个粒子包含一个位置项;计算每个粒子在迭代中的适应度;更新粒子群中的局部最优位置和全局最优位置;根据粒子当前位置、局部最优位置和全局最优位置,计算每个粒子的更新位置;当满足变异条件时,将粒子更新位置进行变异;不断迭代直到达到迭代最大次数,输出最优位置生成力学性能最优隐式曲面多孔结构。本方法通过修改粒子位置移动规则,避免了标准粒子群算法在定义域限制下的缺点,提高了局部最优和全局最优搜索能力。本方法稳定可靠,特别适合隐式曲面多孔结构的参数优化。
Description
技术领域
本发明涉及计算机辅助设计与优化领域,尤其是涉及一种基于改进粒子群算法的隐式曲面多孔结构优化方法。
背景技术
隐式曲面是一种拥有确定函数表达式的数学曲面,又被称为等值面。不同于工业领域广泛应用的参数曲面或者离散网格曲面,隐式曲面可以表示任意复杂拓扑结构的三维曲面,此外复杂的布尔运算可以在隐式曲面函数式上进行快速的计算,隐式曲面在经过几何操作后仍然可以方便地表示为隐函数形式。由于这些重要的优点,越来越多的复杂结构采用隐式曲面来构造。
多孔结构是一种具有复杂拓扑的几何结构,根据不同的分类标准,可以分为开孔结构、闭孔结构、周期性多孔结构、非周期性多孔结构等。这类复杂多孔结构在工程领域有着广泛的应用,错综复杂的内部孔洞可以为热交换、分子吸附提供充足的空间,同时复杂的孔壁可以在电磁波吸收领域发挥重要作用。不同于传统工业中的实体零部件,多孔结构可以在达到力学要求的同时大大降低材料的消耗,实现轻量化的目的,这对于降低制造成本以及后续应用中的资源消耗具有重要的意义。
随着研究人员的深入探索,越来越多形态各异的隐式曲面被提出用于多孔结构的建模与应用。为了由隐式曲面生成具有一定壁厚的多孔实体结构,一般需要将曲面在三维空间进行加厚偏置,Kapfer等人研究了不同实体构造策略导致的力学性能差异(参见Kapfer S C,Hyde S T,Mecke K,et al.Minimal surface scaffold designs for tissueengineering[J].Biomaterials,2011,32(29):6875-6882.),结果证明偏置形成的薄片型隐式曲面多孔结构具有更加理想的力学性能。通过改变隐式曲面函数式的基本参数以及构建多孔实体结构的壁厚,可以方便地控制多孔结构的性能。Almeida和Bártolo离散分析了不同参数下隐式曲面多孔结构的基本力学性能(参见Almeida H A,Bártolo P J.Designof tissue engineering scaffolds based on hyperbolic surfaces:Structuralnumerical evaluation[J].Medical engineering&physics,2014,36(8):1033-1040.),为参数的选择提供了基本的参考。Yan等人金属打印制造了具有不同参数的Gyroid隐式曲面多孔结构(参见Yan C,Hao L,Hussein A,et al.Evaluations of cellular latticestructures manufactured using selective laser melting[J].InternationalJournal of Machine Tools and Manufacture,2012,62:32-38.),验证了孔径大小与压缩性能的关系。
粒子群算法是一种从自然界中得到启发的智能优化算法,具有实现简单、收敛较快的优点,在工程优化问题中得到了广泛的应用。但是传统的标准粒子群迭代算法由于依赖于粒子的位置项与速度项,当粒子位置存在定义域时经常会出现粒子超出可行域的情况,导致粒子在迭代多次后仍处于混沌状态,降低了算法的性能。在一些实际问题中需要对标准粒子群算法进行改进处理。
根据文献分析可知,当前针对隐式曲面多孔结构的研究主要聚焦在几何设计与实验应用方面,一些对基本参数性能影响的研究大多数是基于离散采样的结果,如何在给定条件下寻找最优化的隐式曲面多孔结构参数目前还是研究的空白。此外,未发现任何关于隐式曲面多孔结构改进粒子群算法优化的文献。
发明内容
为了解决在给定条件下隐式曲面多孔结构优化的问题,以及改善标准粒子群算法在存在定义域约束时迭代性能较差的问题,本发明提供了一种基于改进粒子群算法的隐式曲面多孔结构优化方法。本方法稳定可靠,通过不断迭代,可以高效准确地找到使隐式曲面多孔结构相对密度最大化的曲率参数,进而得到最佳的力学性能。
本发明提供的技术方案为:
一种基于改进粒子群算法的隐式曲面多孔结构优化方法,包括以下步骤:
(1)输入待优化隐式曲面函数表达式f(x,y,z)=c,c为曲率参数;
(2)对标准粒子群算法进行改进,具体为:抛弃标准粒子群算法中的速度项,改进粒子位置更新公式,并引入变异概率对更新粒子位置进行变异,获得改进粒子群算法;
(3)采用改进粒子群算法对曲率参数c进行优化,以获得最优曲率参数;
(4)根据最优曲率参数对应的最优隐式曲面函数和多孔实体壁厚生成力学性能最优隐式曲面多孔结构。
改进的粒子群算法,抛弃了速度项,改进了粒子位置更新公式,这样在更新粒子位置时,能够限制粒子更新距离,进而能够保证粒子始终在定义域范围内迭代;同时,为了避免粒子群算法中局部最优位置收敛太快,陷入局部最优,定义一定概率的粒子变异,以此来保证粒子的全局和局部搜索能力。
其中,采用改进的粒子群算法对曲率参数c进行优化的具体步骤为:
(a)初始化N个粒子;
(b)根据隐式曲面函数,计算粒子在第m次迭代的适应度;
(c)根据粒子当前位置的适应度更新粒子群中的局部最优位置和全局最优位置;
(d)根据粒子当前位置、局部最优位置以及全局最优位置对粒子当前位置进行更新,获得粒子更新位置;
(e)在步骤(d)之后,根据变异条件对粒子更新位置进行选择变异;
(f)迭代执行步骤(b)~(e),直到迭代次数达到预设阈值M为止,输出最后一次迭代中获得全局最优位置,该全局最优位置即为最优曲率参数。
具体地,步骤(b)具体为:
然后,再对所述网络曲面进行偏置生成壁厚为t的多孔结构实体;
本发明中将粒子当前位置带入到入隐式曲面函数表达式中,巧妙地将粒子当前位置与曲率参数结合起来,通过该粒子当前位置对应的多孔结构实体体积和包络体体积来确定粒子当前位置的适应度,为更新粒子位置、局部最优位置以及全局最优位置提供基础。
步骤(c)具体为:
步骤(d)具体包括:
步骤(d)中,在更新粒子位置时,抛弃了标准粒子群算法中的速度项,利用当前位置、全局最优位置、局部最优位置以及给定的随机因子来更新粒子位置,由于限定了随机因子取值在0~1之间,即限定了更新距离,确保了粒子位置在定义域范围内,即保证优化参数的可行域。
步骤(e)具体为:
设置变异概率P,变异因子r5∈[0,1];
本发明提供的基于改进粒子群算法的隐式曲面多孔结构优化方法的有益效果表现为:
通过取消粒子速度项,限制粒子在定义域范围内迭代,改进粒子位置更新公式,避免了标准粒子群算法可能出现的混沌情况,同时规避了迭代速度选择的难题,保证了理想的迭代收敛速度。添加粒子位置变异,有效避免算法陷入局部最优,保证了算法整体理想的全局搜索与局部搜索能力。通过不断迭代,可以高效寻找隐式曲面多孔结构最佳的曲率参数,更好地满足实际应用的需求。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图做简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动前提下,还可以根据这些附图获得其他附图。
图1为本发明提供的一种基于改进粒子群算法的隐式曲面多孔结构优化方法流程图;
图2为改进粒子群算法原理示意图;
图3为实施例改进粒子群算法与标准粒子群算法迭代比较结果;
图4为实施例力学压缩应力-应变曲线;
图5为实施例力学性能指标比较结果。
具体实施方式
为使本发明的目的、技术方案及优点更加清楚明白,以下结合附图及实施例对本发明进行进一步的详细说明。应当理解,此处所描述的具体实施方式仅仅用以解释本发明,并不限定本发明的保护范围。
本发明一种基于改进粒子群算法的隐式曲面多孔结构优化方法流程图如图1所示,具体实施步骤如下:
步骤101:输入待优化隐式曲面函数表达式f(x,y,z)=c,x∈[xmin,xmax],y∈[ymin,ymax],z∈[zmin,zmax],其中曲率参数c为待优化的参数,且c∈[cmin,cmax],多孔实体壁厚t,粒子群规模N,最大迭代次数M;
步骤103:计算每个粒子在第m次迭代的适应度m∈[1,M],将第i个粒子在第m次迭代的位置带入隐式曲面函数表达式,利用移动立方体算法生成对应的网格曲面后,再将曲面进行偏置生成壁厚为t的多孔结构实体,计算多孔实体的体积和对应包络体的体积则适应度为多孔结构的相对密度
步骤107:如果迭代次数大于M,输出全局最优位置,该全局最优位置即为最优曲率参数,否则跳至步骤3。
步骤108:根据最优曲率参数对应的最优隐式曲面函数和多孔实体壁厚生成力学性能最优隐式曲面多孔结构。
本发明的典型实施实例如下:
输入隐式曲面I-WP表达式f(x,y,z)=2[cos(0.25πx)cos(0.25πy)+cos(0.25πy)cos(0.25πz)+cos(0.25πz)cos(0.25πx)]-[cos(0.5πx)+cos(0.5πy)+cos(0.5πz)]=c,c∈[-0.6,0.6],x∈[0,8],y∈[0,8],z∈[0,8],多孔实体壁厚0.1mm,粒子群规模30,最大迭代次数30。在英特尔至强3.4GHz处理器上进行粒子群算法迭代,本发明提出的改进粒子群算法和不同速度的标准粒子群算法迭代结果如图3所示,改进粒子群算法在30次迭代内得到了更加理想的结果,最优位置copt=0.00236185。不同速度的设置会导致标准粒子群算法差异明显的搜索性能,本发明提出的改进粒子群算法避免了速度选择的难题,具有较强的局部最优和全局最优搜索能力。最优曲率参数和其他离散参数的力学压缩应力-应变曲线如图4所示,具体曲率参数下的力学性能指标如图5所示,利用最优参数copt生成的隐式曲面多孔结构显然具有最大的弹性模量和平台应力,验证了本发明方法的有效性。
以上所述的具体实施方式对本发明的技术方案和有益效果进行了详细说明,应理解的是以上所述仅为本发明的最优选实施例,并不用于限制本发明,凡在本发明的原则范围内所做的任何修改、补充和等同替换等,均应包含在本发明的保护范围之内。
Claims (6)
1.一种基于改进粒子群算法的隐式曲面多孔结构优化方法,该方法应用在工程领域,包括以下步骤:
(1)输入待优化隐式曲面函数表达式f(x,y,z)=c,c为曲率参数;
(2)对标准粒子群算法进行改进,具体为:抛弃标准粒子群算法中的速度项,改进粒子位置更新公式,并引入变异概率对更新粒子位置进行变异,获得改进粒子群算法;
(3)采用改进粒子群算法对曲率参数c进行优化,以获得最优曲率参数;
(4)根据最优曲率参数对应的最优隐式曲面函数和多孔实体壁厚生成力学性能最优隐式曲面多孔结构,在达到力学要求的同时降低材料的消耗,实现多孔结构的轻量化。
2.如权利要求1所述的基于改进粒子群算法的隐式曲面多孔结构优化方法,其特征在于,采用改进的粒子群算法对曲率参数c进行优化的具体步骤为:
(a)初始化N个粒子;
(b)根据隐式曲面函数,计算粒子在第m次迭代的适应度;
(c)根据粒子当前位置的适应度更新粒子群中的局部最优位置和全局最优位置;
(d)根据粒子当前位置、局部最优位置以及全局最优位置对粒子当前位置进行更新,获得粒子更新位置;
(e)在步骤(d)之后,根据变异条件对粒子更新位置进行选择变异;
(f)迭代执行步骤(b)~(e),直到迭代次数达到预设阈值M为止,输出最后一次迭代中获得全局最优位置,该全局最优位置即为最优曲率参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810865521.1A CN109145408B (zh) | 2018-08-01 | 2018-08-01 | 一种基于改进粒子群算法的隐式曲面多孔结构优化方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810865521.1A CN109145408B (zh) | 2018-08-01 | 2018-08-01 | 一种基于改进粒子群算法的隐式曲面多孔结构优化方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109145408A CN109145408A (zh) | 2019-01-04 |
CN109145408B true CN109145408B (zh) | 2020-12-25 |
Family
ID=64799466
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810865521.1A Active CN109145408B (zh) | 2018-08-01 | 2018-08-01 | 一种基于改进粒子群算法的隐式曲面多孔结构优化方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109145408B (zh) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106229964A (zh) * | 2016-07-22 | 2016-12-14 | 南京工程学院 | 一种基于改进二进制粒子群算法的配电网故障定位方法 |
CN107530026A (zh) * | 2015-05-12 | 2018-01-02 | 海珀菲纳研究股份有限公司 | 射频线圈方法和装置 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8874473B2 (en) * | 2008-10-31 | 2014-10-28 | At&T Intellectual Property I, Lp | System and method for managing e-commerce transaction |
-
2018
- 2018-08-01 CN CN201810865521.1A patent/CN109145408B/zh active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107530026A (zh) * | 2015-05-12 | 2018-01-02 | 海珀菲纳研究股份有限公司 | 射频线圈方法和装置 |
CN106229964A (zh) * | 2016-07-22 | 2016-12-14 | 南京工程学院 | 一种基于改进二进制粒子群算法的配电网故障定位方法 |
Non-Patent Citations (2)
Title |
---|
《Acoustic band gaps and elastic stiffness of PMMA cellular solids based on triply periodic minimal surfaces》;Abueidda etc;《MATERIALS & DESIGN》;20180505;第145卷;全文 * |
《基于粒子群算法的页岩孔隙结构反演及横波速度预测》;刘财 等;《地球物理学进展》;20171231;第32卷(第2期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN109145408A (zh) | 2019-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Pehlivanoglu | A new particle swarm optimization method enhanced with a periodic mutation strategy and neural networks | |
Domingues et al. | Space–time adaptive multiresolution methods for hyperbolic conservation laws: Applications to compressible Euler equations | |
Yang et al. | Polynomial chaos expansion for probabilistic uncertainty propagation | |
CN109145408B (zh) | 一种基于改进粒子群算法的隐式曲面多孔结构优化方法 | |
Zhou et al. | A comparative study of genetic and firefly algorithms for sensor placement in structural health monitoring | |
CN109492816B (zh) | 一种基于混合智能的煤与瓦斯突出动态预测方法 | |
Chu et al. | Reliability based optimization with metaheuristic algorithms and Latin hypercube sampling based surrogate models | |
Horie et al. | Physics-embedded neural networks: E (n)-equivariant graph neural pde solvers | |
CN105490836B (zh) | 一种复杂网络可靠度的蒙特卡罗评估方法 | |
CN116522542A (zh) | 罐体结构的确定方法、装置、电子设备和存储介质 | |
CN115453867A (zh) | 一种鲁棒的自适应大规模气力输送控制方法 | |
CN109145409A (zh) | 一种隐式曲面多孔结构的优化方法 | |
Friedrich | A new method for generating inner points of triangulations in two dimensions | |
CN117010260A (zh) | 一种裂缝性油藏自动历史拟合模型预测方法、系统及设备 | |
Fang et al. | Network Traffic Prediction Model Based on Catfish-PSO-SVM. | |
Wang et al. | Application of BA-BP neural network in surface reconstruction | |
Zhong et al. | An information geometry algorithm for distribution control | |
Song et al. | A New Hybrid Update Scheme for an Evolutionary Search Strategy Using Genetic Algorithm and Kriging | |
Chen et al. | A quadrilateral mesh generating method using mean value interpolation | |
US20230351065A1 (en) | System and method for heat exchanger shape optimization | |
Nie et al. | Multi-swarm hybrid optimization algorithm with prediction strategy for dynamic optimization problems | |
ZHanG et al. | Water resources big data classification based on multi-objective optimization for mining area. | |
Li et al. | EHPSO: An Enhanced Hybrid Particle Swarm Optimization Algorithm for Internet of Things. | |
Kobayashi et al. | Shell topology optimization based on level set method | |
CN117556557A (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 |