[发明专利]一种RDPRQCG的大规模结构拓扑频率优化方法

说明书

本申请公开了一种RDPRQCG的大规模结构拓扑频率优化方法，涉及结构优化技术领域。该方法包括：采用GPU与提出的变缩比和对角预处理瑞利商共轭梯度算法进行并行计算，获得大规模结构少量固有频率和振型；基于计算得到的固有频率和振型，以结构柔顺度最小为目标函数，同时考虑频率和体积约束，给出大规模结构拓扑优化模型、求解方法和典型案例。本申请相比于传统对角预处理瑞利商共轭梯度的结构拓扑频率优化方法可以提升计算速度，大大节省计算时间，节省计算内存，有效地避免规模较大时设备内存不足的问题。

技术领域

本申请属于结构优化技术领域，具体涉及一种RDPRQCG(Reduction ratio anddiagonal preconditioning Rayleigh quotient conjugate gradient，变缩比和对角预处理的瑞利商共轭梯度)的大规模结构拓扑频率优化方法。

背景技术

在考虑频率要求的大规模结构拓扑优化研究中，不可避免地涉及到大规模结构少量固有频率和振型计算(即广义特征值问题进行求解)。大规模结构总自由度数N很大,大规模结构拓扑优化过程中,广义特征值问题的系数矩阵常出现病态现象。求解此类广义特征值问题的部分特征对的难度是巨大的,且传统对角预处理器可以看作是提取K的对角元素，忽略该未知元素(用这个对角元素表示)对其他未知元素可能存在的耦合作用。基于传统对角预处理器的结构拓扑优化设计的效率不高，且难于适应规模结构拓扑优化过程中含有病态系数矩阵的结构固有频率和振型计算。

另外，当采用传统的串行拓扑优化算法对结构进行优化设计时，组装整体刚度矩阵的效率直接受到单元数量的影响。随着单元数量的增加，基于整体刚度矩阵和质量矩阵求解广义特征值问题的部分特征对的时间急剧上升，每个优化迭代步需要几十分钟甚至更久。因此，借助最新发展的具有多线程并行计算的GPU，研发既能减少迭代次数，且又能适应含有病态系数矩阵的结构固有频率和振型并行计算的方法是很有必要的。基于GPU的大规模结构拓扑频率优化方法会极大地提高计算效率，满足大规模工程结构设计要求。

发明内容

本申请实施的目的是提供一种RDPRQCG的大规模结构拓扑频率优化方法。首先，提出了一种新的变缩比和对角预处理器；并将该预处理技术与GPU并行计算方式和瑞利商共轭梯度法结合，提出了一种基于GPU的RDPRQCG大规模结构拓扑频率优化方法和其并行计算技术，该方法和其技术显著地提高了结构拓扑频率优化效率，进而可以解决背景技术中涉及的大规模结构拓扑频率优化问题。

为了解决上述技术问题，本申请是这样实现的：

本申请实施例提供了一种RDPRQCG的大规模结构拓扑频率优化方法，包括：

步骤S1、采用GPU与提出的变缩比和对角预处理的正交加速瑞利商共轭梯度算法进行并行计算，获得大规模结构少量固有频率和振型；

步骤S2、基于并行计算获得的大规模结构少量固有频率和振型，以结构柔顺度最小为目标函数，同时考虑频率和体积约束，给出大规模结构拓扑优化模型、并行求解方法和典型案例。

可选的，所述采用GPU与提出的RDPRQCG算法进行并行计算，获得大规模结构少量固有频率和振型，包括：

步骤S11、依据结构总刚度矩阵K，采用式(1)构建新的对角线预处理矩阵的对角元素

其中，是总体坐标系中扩维的结构单元刚度矩阵K^(e)的第i行第j列元素，m为结构单元个数,N为结构自由度数；

步骤S12、采用式(1)构建由下列对角正定矩阵定义的对角线预处理矩阵

其中，

其中，δ是一个小量，取δ＝10^-6；