[发明专利]一种关于大长径比圆柱壳结构拓扑优化的过滤变量设计方法

摘要

本发明提供一种关于大长径比圆柱壳结构拓扑优化的过滤变量设计方法。本发明用于解决大长径比圆柱壳结构的拓扑优化中考虑稳定性条件的技术问题，通过引入等效长度概念设计过滤函数，过滤出影响结构稳定性的重要网格单元，并将过滤后的单元密度作为设计变量进行优化设计。解决了增加稳定性约束条件优化结果不佳问题，解决了在拓扑优化中引入大量参数问题，并将稳定性条件引入到设计变量中，从而实现满足稳定性条件的拓扑优化分析。

主权项

1.一种关于大长径比圆柱壳结构拓扑优化的过滤变量设计方法，其特征在于，包括如下步骤：S1、有限元建模：建立大长径比的圆柱壳结构初始设计模型，定义载荷和边界条件；S2、过滤设计变量：对于载荷方向上处于同一列的单元，其伪密度在载荷反方向的顺序依次为ρ(x₁)、ρ(x₂)、…、ρ(x_n)，通过设计过滤函数S将部分单元过滤出去，其中过滤函数为：S(δ,m,i)＝δm+i                 (1)式中，δ为等效长度，m＝0，1，…，M，M为过滤单元的段数，i＝0，1，…，I，I为每段的过滤单元列数，以上参数的确定方法如下：原始等效长度δ₀：式中，H为有限元模型载荷方向的长度，F_cr为体积约束下柔度最小化优化模型的临界力，F为有限元模型所受载荷，β为单元长度，过滤单元的段数M：式中，为取整函数，等效长度δ：每段的过滤单元列数I，依次赋值为1、2、3…，令i＝1，i＝i+1，通过设计停止准则取值I＝i，停止准则为：if ψ_i+1(ρ(x_n))≥ψ_i(ρ(x_n))            (5)通过式(1)～式(5)可得S的一系列数值，用来过滤有相同下标的x_s，其过滤的单元集合为：x_s＝{x_δm+i|m＝0，1，…，M，i＝0，1，…，I}过滤出的x_s为有限元模型满足稳定性的重要单元，则令ρ(x_s)＝1；S3、建立拓扑优化模型：在S2中筛选出x_s后，x_i＝x_n‑x_s为即将参与拓扑优化设计的单元变量，从而建立拓扑优化模型，式中，ρ(x_i)为设计变量的密度,x_n为全部单元集合，x_s为过滤的单元，ψ(ρ(x_n))为目标函数，K为有限元模型总体刚度矩阵，U为节点整体位移向量，F为节点等效载荷向量，H_j(x_n)为第j个约束函数，为第j个约束函数的上限，j为约束的数量，ε＝10^‑3以避免刚度奇异值；通过上述优化模型对大长径比圆柱壳结构进行稳定性过滤和优化设计得到最优结果。