[发明专利]一种面向各向异性材料的连续体结构双材料拓扑优化方法

权利要求书

1.一种面向各向异性材料的连续体结构双材料拓扑优化方法，其特征在于：实现步骤如下：

步骤一：基于经典的带罚函数的材料插值模型，构建适用于各向异性材料的双材料插值模型：

其中为插值后第i个单元的单元刚度矩阵，为第i个单元对应于材料1的单元刚度矩阵，为第i个单元对应于材料2的单元刚度矩阵，x_1,i和x_2,i分别为第i个单元的设计变量1和设计变量2，p为惩罚因子，p＞1；

步骤二：利用弹性常数的插值模型去间接构建两种材料刚度矩阵相融合的插值模型：

其中，E_i,x1、E_i,x2、E_i,y1、E_i,y2、E_i,z1、E_i,z2、μ_i,xy1、μ_i,xy2、μ_i,yz1、μ_i,yz2、μ_i,xz1、μ_i,xz2、G_i,xy1、G_i,xy2、G_i,yz1、G_i,yz2、G_i,xz1、G_i,xz2分别是对应于材料1和材料2的弹性系数，依次为三个坐标方向的弹性模量E、泊松比μ和剪切模量G，其中，下标i表示第i个单元，下标1对应于材料1，下标2对应于材料2；

步骤三：基于拓扑优化数学模型，以最小化结构柔度作为目标，以结构质量为约束，建立各向同性与各向异性一体的连续体结构双材料拓扑优化的数学模型：

其中，C是结构柔度，V_i为第i个单元的体积，u_j为结构的第j个位移，f_j为结构的第j个载荷，m为载荷的个数，n为设计域划分的单元总数，V_t1和V_t2分别为材料1和材料2的目标体积，x_1,i为第i个单元的设计变量1，x_2,i为第i个单元的设计变量2；x₁和x₂分别为设计变量1和设计变量2的下界，和分别为设计变量1和设计变量2的上界；

步骤四：基于伴随向量法求解目标函数结构柔度对设计变量的灵敏度：

其中，u_e,i为第i个单元的单元位移向量；

步骤五：使用移动渐近线优化算法(Method of Moving Asymptotes)，以最小化结构柔度为目标，以结构相对质量分数为约束，利用柔度和相对质量分数对涉及变量的灵敏度进行迭代求解，在迭代过程中，如果当前设计不满足相对质量约束，或前后两次迭代的设计变量变化绝对值之和大于预设值ε时，则返回步骤二进行新一轮的迭代优化，否则，进行步骤六；

步骤六：如果当前设计不满足相对质量约束，或前后两次迭代的设计变量变化绝对值之和大于预设值ε时，则迭代结束，得到双材料连续体结构非概率可靠性拓扑优化的最优构型。

2.根据权利要求1所述的一种面向各向异性材料的连续体结构双材料拓扑优化方法，其特征在于：所述步骤一中的材料插值模型采用了适用于各向异性材料的双材料插值模型。

3.根据权利要求1所述的一种面向各向异性材料的连续体结构双材料拓扑优化方法，其特征在于：所述步骤二中利用两种材料的弹性常数插值模型去间接构建两种材料刚度矩阵相融合的材料插值模型。

4.根据权利要求1所述的一种面向各向异性材料的连续体结构双材料拓扑优化方法，其特征在于：所述步骤三中构建的各向同性与各向异性一体的连续体结构双材料拓扑优化的数学模型以最小化结构柔度作为目标，以结构质量为约束。

5.根据权利要求1所述的一种面向各向异性材料的连续体结构双材料拓扑优化方法，其特征在于：所述步骤四基于伴随向量法求解目标函数对设计变量的灵敏度。