一种隐式曲面多孔结构的优化方法技术

本发明专利技术公开了一种隐式曲面多孔结构的优化方法，包括输入优化隐式曲面函数表达式，隐式曲面多孔结构相对密度，粒子群规模，最大迭代次数；初始化粒子，每个粒子位置坐标包含隐式曲面多孔结构的曲率参数和壁厚参数；计算迭代粒子的约束适应度；更新粒子群中的局部最优位置和全局最优位置；线性加权计算每个粒子的更新位置；满足突变要求时粒子位置中的曲率参数取相反数；当超过最大迭代次数时停止迭代，输出最优粒子位置，生成最优参数对应的隐式曲面多孔结构。本方法以线性加权方式进行粒子迭代，同时添加突变概率，保证粒子的随机行为，方法迭代速度较快，稳定可靠，迭代得到曲率参数和壁厚参数可以生成满足多功能需求的隐式曲面多孔结构。

【技术实现步骤摘要】
一种隐式曲面多孔结构的优化方法
本专利技术涉及计算机辅助设计与多功能优化领域，尤其是涉及一种隐式曲面多孔结构的优化方法。
技术介绍
当前在计算机辅助设计领域，实体结构的建模与分析优化已经研究较为成熟，大量的专业建模软件和分析软件被开发出来用于结构设计与优化。事实上，在自然界中，大部分天然结构中存在大量的多孔结构，例如生物的骨骼就是典型的多孔结构。多孔结构一般拥有错综复杂的内部孔洞，这些孔洞的结构大多是自然选择的结果，随着外界条件的变化，多孔结构在不断生成变化中成为具有最理想性能的复杂结构。为了精确描述多孔结构的形态特征，数学家开发了很多工具来实现参数化的设计建模。其中，隐式曲面是一种优秀的数学建模工具，运用简洁的数学函数表达式，就可以准确描绘多孔结构的复杂特征。近年来，研究人员提出了多种优秀的隐式曲面用于建模自然界常见的多孔结构，例如Gyroid曲面、Diamond曲面以及Schwarz曲面等。通过改变曲面的曲率参数以及形成多孔实体的偏置壁厚，可以方便地控制隐式曲面多孔结构的关键性能，满足不同的工程需求。针对隐式曲面多孔结构优良的孔洞特征，国内外学者研究了不同的应用案例。Abueidda等人研究了隐式曲面多孔结构的基本力学性能(参见Abueidda,DiabW.,etal."Mechanicalpropertiesof3Dprintedpolymericcellularmaterialswithtriplyperiodicminimalsurfacearchitectures."Materials&Design122(2017):255-267.)，论证了将多孔结构用于承受外部载荷的应用前景。Abueidda等人研究了隐式曲面多孔结构的声学带隙(参见AbueiddaDW,JasiukI,SobhNA.AcousticbandgapsandelasticstiffnessofPMMAcellularsolidsbasedontriplyperiodicminimalsurfaces[J].Materials&Design,2018,145:20-27.)，可以将此类结构用作吸音装置。Melchels等人将隐式曲面多孔结构用作组织工程支架(参见Melchels,FerryPW,etal."Effectsofthearchitectureoftissueengineeringscaffoldsoncellseedingandculturing."Actabiomaterialia6.11(2010):4208-4217.)，内部相互贯通的孔洞以及较大的表面积可以为细胞增殖生长提供理想的环境。一般来说，评价隐式曲面多孔结构主要有两个指标，首先是结构的力学性能，根据Gibson的经典多孔力学理论，多孔结构的力学性能和相对密度成正相关关系。在一般的应用环境中，需要结构有一定的力学强度，即在设计过程中，要保证多孔结构的相对密度达到一定程度。其次是多孔的结构性能，如常见的连通率、渗透率以及比表面积，其中比表面积应用最为广泛，是散热、吸声、组织培养应用中最为关心的性能指标。随着应用环境的日益复杂，工程实际对多孔结构的性能提出了越来越高的要求，如何让多孔结构满足多功能化的需求成为了研究的难点。根据文献分析可知，当前对隐式曲面多孔结构的性能研究主要在单一应用性能的探索，多功能一体化的设计与优化目前还是研究的空白。然而，目前多孔结构不仅应具有一定的力学性能，同时还应具有理想的结构特性，如何优化隐式曲面多孔结构的参数以获得多功能化的特性目前还没有理想的解决方案。此外，未发现任何关于约束粒子群的隐式曲面多孔结构多功能优化的文献。
技术实现思路
为了解决现有多孔结构多功能优化困难的问题，本专利技术提供了一种隐式曲面多孔结构的优化方法。本专利技术的技术方案为：一种隐式曲面多孔结构的优化方法，包括以下步骤：(1)输入待优化隐式曲面函数表达式φ(x，y，z)＝c，c为曲率参数；(2)初始化L个粒子，第i个粒子初始位置坐标为其中，为隐式曲面的曲率参数，为隐式曲面多孔结构的壁厚参数，i∈[1，L]；(3)根据粒子当前位置，计算获得第k次迭代中第i个粒子的约束适应度k∈[0，K]，K为最大迭代次数；(4)根据粒子的约束适应度更新粒子群中的局部最优位置和全局最优位置(5)根据粒子当前位置局部最优位置和全局最优位置线性加权计算粒子更新位置(6)在步骤(5)之后，保持粒子位置中壁厚参数项不变，根据变异条件对粒子更新位置进行选择变异；(7)迭代执行步骤(3)～(6)，直到迭代次数达到最大迭代次数K为止，输出最后一次迭代中获得的全局最优位置，该全局最优位置即为最优曲率与壁厚参数；(8)根据最优曲率参数对应的最优隐式曲面函数和多孔实体壁厚生成相对密度一定、比表面积最大的隐式曲面多孔结构。本专利技术提供的隐式曲面多孔结构的优化方法实现方便，稳定可靠，可以以较快的收敛速度找到隐式曲面多孔结构的合理参数组合，进而能够获得指定相对密度、比表面积最大化的隐式曲面多孔结构。优选地，步骤(3)中，根据粒子当前位置利用移动立方体算法生成隐式曲面对应的网格曲面，并偏置成壁厚为的多孔实体结构，再根据该多孔实体结构计算第k次迭代中第i个粒子的约束适应度移动立方体算法基本思想是在体数据的每一个立方体单元中根据其八个顶点的数据值与给定数据值的关系在单元的12条边上寻找等值点，然后用三角形将等值点连成等值面。移动立方体算法对感兴趣的等值面可以产生清晰的图像，因此能够快速稳定获得三角网格曲面。具体地，所述根据该多孔实体结构计算第k次迭代中第i个粒子的约束适应度包括：首先，计算多孔实体结构的相对密度和比表面积然后，计算相对密度误差值和比表面积数值数量级差值n，其中，ρr为隐式曲面多孔结构的预定相对密度，为相对密度误差值，[·]表示取整；最后，计算第k次迭代中第i个粒子的约束适应度为：本专利技术中，根据多孔实体结构的相对密度和比表面积计算约束适应度，即通过给粒子适应度函数增加约束项保证在达到指定相对密度的同时使结构的比表面积最大化。迭代寻优过程是使粒子的约束适应度达到最大值，即当为零时，且达到最大值时，适应度为最大结果，也就是当相对密度等于预设相对密度ρr时，寻找多孔实体结构的比表面积达到最大值的对应参数。粒子适应度函数约束项中的10n能够加快迭代收敛的速度，使迭代优化结果更稳定。优选地，步骤(4)包括：如果粒子当前位置的适应度大于局部最优位置的适应度，则将局部最优位置更新为粒子当前位置即如果当前粒子位置的适应度大于全局最优位置的适应度，则将全局最优位置更新为粒子当前位置即步骤(5)中，根据以下公式计算得到粒子更新位置其中，随机因子α∈[0，1]，随机因子β∈[0，1]。在本专利技术中，在获得粒子更新位置时，采用三项加权的方法，即局部最优位置全局最优位置以及粒子当前位置这三项前分别乘以权重α、β以及1-α-β，再以三项加权的和作为粒子更新位置，这样保证了粒子更新位置的稳定性和相关性。具体地，步骤(6)包括：保持粒子位置中壁厚参数项不变；设置变异概率g∈[0，1]，随机变异因子δ∈[0，1]；当δ≥g时，粒子更新位置发生变异，即粒子更新位置变异为当δ＜g时，粒子更新位置不发生变异。本专利技术本文档来自技高网...

【技术保护点】
1.一种隐式曲面多孔结构的优化方法，包括以下步骤：(1)输入待优化隐式曲面函数表达式φ(x，y，z)＝c，c为曲率参数；(2)初始化L个粒子，第i个粒子初始位置坐标为

【技术特征摘要】
1.一种隐式曲面多孔结构的优化方法，包括以下步骤：(1)输入待优化隐式曲面函数表达式φ(x，y，z)＝c，c为曲率参数；(2)初始化L个粒子，第i个粒子初始位置坐标为其中，为隐式曲面的曲率参数，为隐式曲面多孔结构的壁厚参数，i∈[1，L]；(3)根据粒子当前位置，计算获得第k次迭代中第i个粒子的约束适应度k∈[0，K]，K为最大迭代次数；(4)根据粒子的约束适应度更新粒子群中的局部最优位置和全局最优位置(5)根据粒子当前位置局部最优位置和全局最优位置线性加权计算粒子更新位置(6)在步骤(5)之后，保持粒子位置中壁厚参数项不变，根据变异条件对粒子更新位置进行选择变异；(7)迭代执行步骤(3)～(6)，直到迭代次数达到最大迭代次数K为止，输出最后一次迭代中获得的全局最优位置，该全局最优位置即为最优曲率与壁厚参数；(8)根据最优曲率参数对应的最优隐式曲面函数和多孔实体壁厚生成相对密度一定、比表面积最大的隐式曲面多孔结构。2.如权利要求1所述的隐式曲面多孔结构的优化方法，其特征在于，步骤(3)中，根据粒子当前位置利用移动立方体算法生成隐式曲面对应的网格曲面，并偏置成壁厚为的多孔实体结构，再根据该多孔实体结构计算第k次迭代中第i个粒...

【专利技术属性】
技术研发人员：冯嘉炜，傅建中，林志伟，商策，
申请(专利权)人：浙江大学，
类型：发明
国别省市：浙江,33