一种连续体结构拓扑与形状尺寸综合优化方法技术

本发明专利技术涉及一种连续体结构拓扑与形状尺寸综合优化方法。利用有限元思想对设计域进行网格划分并提取单元的节点信息；创建优化目标函数、添加约束和荷载，通过计算目标函数灵敏度并完成灵敏度过滤，来迭代更新设计域单元相对密度值；判断达到收敛条件时，输出具有灰度单元的拓扑构型；然后将灰度单元数值置为1，同时将结构内部细小孔洞填满，形成离散0‑1矩阵形式的初始拓扑构型；引入边界函数对初始拓扑构型进行形状尺寸优化，最终获得满足约束条件的最优设计结构。本发明专利技术解决了传统拓扑优化中灰度单元对结构可制造性的影响，同时引入边界函数对初始拓扑构型进行形状尺寸优化以消除嵌入结构内部的细小孔洞，减小了连续体结构在制造过程中的复杂度。

A comprehensive optimization method of topology and shape size of continuum structure

【技术实现步骤摘要】
一种连续体结构拓扑与形状尺寸综合优化方法
本专利技术涉及结构优化设计的相关
，具体涉及一种连续体结构拓扑与形状尺寸综合优化方法。
技术介绍
在对连续体结构进行拓扑优化设计中，常用的方法有基于有限元思想的密度法和基于边界传播思想的水平集法。依据密度法进行拓扑优化，是将连续体结构设计域划分成有限个单元，并对每个单元赋以相对密度值作为优化设计的变量。如果将单元相对密度值离散成0-1的阶跃函数描述，并按一定进化比来删除或生成单元，得到的优化结果将总是依赖于划分网格的疏密程度，并在很大程度上无法获得最优解。如果将单元相对密度值用[0,1]范围内的连续函数描述，得到的优化结果将出现大量灰度单元。而实际工程结构在制造中往往都是使用单一材料进行加工制作，因此灰度单元的存在将影响连续体结构的真实形状尺寸。依据水平集法进行拓扑优化，将连续体结构设计域的边界轮廓用高维标量函数描述，并对其求取时间梯度以获得边界演化的速度方程来反映结构的变化情况。但水平集法是对连续体结构的边界轮廓进行演化，无法有利描述结构内部的演化过程，需要对初始设计域创建许多内部孔洞来协助完成优化工作，因此连续体结构设计域的初始构型将在很大程度上影响最终的优化结果。
技术实现思路
本专利技术针对现连续体结构拓扑优化设计中的上述缺点与不足，提供了一种连续体结构拓扑与形状尺寸综合优化方法，通过密度法获得拓扑构型，然后对其结果中的所有灰度单元数值均置为1，同时填补结构内部的细小孔洞，得到单元相对密度值为离散0-1形式且无细小孔洞的拓扑构型，再通过构造边界函数进行二次拓扑，以完成连续体结构设计域的形状尺寸优化，最终得到具有可制造性的最优结构。本专利技术的上述目的是通过以下技术方案得以实施的：一种连续体结构拓扑与形状尺寸综合优化方法，包括：S1、为连续体结构创建设计域，进行结构离散化和有限元网格划分，得到节点与单元的编号信息；S2、通过有限元分析，得到单元的节点位移与节点力信息；S3、创建需要优化的目标函数，并添加实际工况下的约束条件和荷载信息；S4、计算单元目标函数的灵敏度，并完成灵敏度过滤，根据新的灵敏度数值来更新设计域的单元相对密度；S5、迭代收敛判断，当优化的目标达到收敛标准时，优化迭代结束，执行步骤S6，否则重复步骤S2至步骤S5；S6、提取设计域单元相对密度值，将所有灰度单元的相对密度值和嵌入结构内部的细小孔洞单元相对密度值均置为1，以形成离散0-1矩阵形式；S7、根据获得的拓扑构型创建边界函数；S8、通过有限元分析，得到单元的节点位移与节点力信息；S9、计算当前结构的形状灵敏度和拓扑灵敏度；S10、迭代收敛判断，当优化的目标达到收敛标准时，优化迭代结束，执行步骤S12，否则执行步骤S11；S11、根据新的灵敏度数值来更新设计域的单元相对密度，并重复执行步骤S7至步骤S10；S12、输出优化结果，得到具有可制造性的最优连续体结构。作为进一步优选地，所述步骤S1中，在对设计域进行网格划分时，应留意增加网格带来的经济性，根据实际应用比较前后两次网格划分的计算结果，如果两次计算结果相差较大，应当继续增加网格并重新进行拓扑优化的计算，相反则满足要求停止计算。作为进一步优选地，所述步骤S2和步骤S8中，获得的节点位移信息包括但不限于节点编号及坐标、单元的节点位移分量大小及方向信息。作为进一步优选地，所述步骤S2和步骤S8中，获得的节点力信息包括但不限于单元边界上的表面力、单元上的体积力和集中力信息通过等效方式移植到节点上去，利用等效节点力来替代所有作用在单元上的力。作为进一步优选地，所述步骤S3中，目标函数的优化目标包括但不限于结构刚度特性优化、结构固有频率性能优化、结构动态特性优化。作为进一步优选地，所述步骤S3中，实际工况的约束条件包括但不限于结构位移约束条件、结构材料属性约束条件、结构制造加工精度约束条件。作为进一步优选地，所述步骤S4中，在单元相对密度值中引入惩罚因子p和最小过滤半径rmin，采用下式计算目标函数灵敏度：同时，通过下式完成目标函数灵敏度过滤：式中，F为目标函数，x为单元相对密度，[U]为单元节点位移矩阵，[K]为单元刚度矩阵，dist(e,f)为两节点之间的最短距离，e、f为单元节点编号，N为连续体结构进行有限元划分所得到的总单元数目，Hf为卷积算子。作为进一步优选地，所述步骤S4和步骤S11中，对结构优化的数学模型构造成下式的拉格朗日函数形式，采用标准KKT(Karush-Kuhn-Tucker)最优化条件完成迭代优化，来更新设计域的单元相对密度数值：式中，λi为常系数，hi(x)为结构优化中添加的约束条件。作为进一步优选地，所述步骤S5和步骤S10中，通过下式对前后两次优化中的设计域体积改变量进行计算，当改变量小于设定标准值时达到收敛标准，从而跳出迭代循环：|Vnew-Vreq|≤ε式中，Vnew为当前结构优化得到的拓扑构型体积量，Vreq为上一次结构优化得到的拓扑构型体积量，ε为设定的前后两次优化结果体积改变量的标准值。作为进一步优选地，所述步骤S6中，借助但不限于编程的方式，实现将所有非空洞区域的单元相对密度值和嵌入结构内部的细小孔洞单元相对密度值均置为1，将大孔洞区域的单元相对密度数值均置为0，形成0-1离散形式的单元相对密度矩阵。作为进一步优选地，所述步骤S7中，通过下式提取初始拓扑构型中的边界轮廓，完成边界函数的创建：式中，为高维标量函数ψ对时间t进行求导以反映结构的边界演化情况，v为单元的相应敏感度。作为进一步优选地，所述步骤S9中，通过下式计算形状灵敏度和拓扑灵敏度：式中，Ω为结构设计域，λ和μ为固体材料的拉梅常数。与现有技术相比，本专利技术的优点是解决了密度法带来的灰度单元影响和水平集方法创建新孔困难的缺点，通过密度法获得初始拓扑构型，并将灰度单元转换成0-1离散数值矩阵，同时填补结构内部的细小孔洞，通过边界函数完成对初始拓扑构型的形状尺寸优化，从而确保得到的优化结果具有工程可制造性。附图说明图1为本专利技术一种连续体结构拓扑与形状尺寸综合优化方法的流程图；图2为本专利技术实施例1，二维矩形初始设计域的示意图；图3为本专利技术实施例1的迭代过程中，设计域具有灰度单元和细小孔洞的拓扑构型图；图4为本专利技术实施例1的迭代过程中，设计域转为0-1离散形式的拓扑构型图；图5为本专利技术实施例1的迭代终止时，完成结构拓扑与形状尺寸优化的最优结果图；图6为本专利技术实施例2，二维矩形初始设计域的示意图；图7为本专利技术实施例2的迭代过程中，设计域具有灰度单元的拓扑构型图；图8为本专利技术实施例2的迭代过程中，设计域为0-1离散形式的拓扑构型图；图9为本专利技术实施例2的迭代终止时，完成结构拓扑与本文档来自技高网...

【技术保护点】
1.一种连续体结构拓扑与形状尺寸综合优化方法，其特征在于，包括如下步骤：/nS1、为连续体结构创建设计域，进行结构离散化和有限元网格划分，得到节点与单元的编号信息；/nS2、通过有限元分析，得到单元的节点位移与节点力信息；/nS3、创建需要优化的目标函数，并添加实际工况下的约束条件和荷载信息；/nS4、计算单元目标函数的灵敏度，并完成灵敏度过滤，根据新的灵敏度数值来更新设计域的单元相对密度；/nS5、迭代收敛判断，当优化的目标达到收敛标准时，优化迭代结束，执行步骤S6，否则重复步骤S2至步骤S5；/nS6、提取设计域单元相对密度值，将所有灰度单元的相对密度值和嵌入结构内部的细小孔洞单元相对密度值均置为1，以形成离散0-1矩阵形式；/nS7、根据获得的拓扑构型创建边界函数；/nS8、通过有限元分析，得到单元的节点位移与节点力信息；/nS9、计算当前结构的形状灵敏度和拓扑灵敏度；/nS10、迭代收敛判断，当优化的目标达到收敛标准时，优化迭代结束，执行步骤S12，否则执行步骤S11；/nS11、根据新的灵敏度数值来更新设计域的单元相对密度，并重复执行步骤S7至步骤S10；/nS12、输出优化结果，得到具有可制造性的最优连续体结构。/n...

【技术特征摘要】
1.一种连续体结构拓扑与形状尺寸综合优化方法，其特征在于，包括如下步骤：
S1、为连续体结构创建设计域，进行结构离散化和有限元网格划分，得到节点与单元的编号信息；
S2、通过有限元分析，得到单元的节点位移与节点力信息；
S3、创建需要优化的目标函数，并添加实际工况下的约束条件和荷载信息；
S4、计算单元目标函数的灵敏度，并完成灵敏度过滤，根据新的灵敏度数值来更新设计域的单元相对密度；
S5、迭代收敛判断，当优化的目标达到收敛标准时，优化迭代结束，执行步骤S6，否则重复步骤S2至步骤S5；
S6、提取设计域单元相对密度值，将所有灰度单元的相对密度值和嵌入结构内部的细小孔洞单元相对密度值均置为1，以形成离散0-1矩阵形式；
S7、根据获得的拓扑构型创建边界函数；
S8、通过有限元分析，得到单元的节点位移与节点力信息；
S9、计算当前结构的形状灵敏度和拓扑灵敏度；
S10、迭代收敛判断，当优化的目标达到收敛标准时，优化迭代结束，执行步骤S12，否则执行步骤S11；
S11、根据新的灵敏度数值来更新设计域的单元相对密度，并重复执行步骤S7至步骤S10；
S12、输出优化结果，得到具有可制造性的最优连续体结构。


2.如权利要求1所述的连续体结构拓扑与形状尺寸综合优化方法，其特征在于，在步骤S1中，在对设计域进行网格划分时，应留意增加网格带来的经济性，根据实际应用比较前后两次网格划分的计算结果，如果两次计算结果相差较大，应当继续增加网格并重新进行拓扑优化的计算，相反则满足要求停止计算。


3.如权利要求1所述的连续体结构拓扑与形状尺寸综合优化方法，其特征在于，在步骤S2和步骤S8中，获得的节点位移信息包括但不限于节点编号及坐标、单元的节点位移分量大小及方向信息。


4.如权利要求1所述的连续体结构拓扑与形状尺寸综合优化方法，其特征在于，在步骤S2和步骤S8中，获得的节点力信息包括但不限于单元边界上的表面力、单元上的体积力和集中力信息通过等效方式移植到节点上去，利用等效节点力来替代所有作用在单元上的力。


5.如权利要求1所述的连续体结构拓扑与形状尺寸综合优化方法，其特征在于，在步骤S3中，目标函数的优化目标包括但不限于结构刚度特性优化、结构固有频率性能优化、结构动态特性优化。


6.如权利要求1所述的连续体结构拓扑与形状尺...

【专利技术属性】
技术研发人员：张军峰，刘恩海，廖靖宇，朱自发，杨金利，雷华舟，
申请(专利权)人：中国科学院光电技术研究所，
类型：发明
国别省市：四川;51