一种RDPRQCG的大规模结构拓扑频率优化方法技术

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

【技术实现步骤摘要】
一种RDPRQCG的大规模结构拓扑频率优化方法


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

技术介绍

[0002]在考虑频率要求的大规模结构拓扑优化研究中,不可避免地涉及到大规模结构少量固有频率和振型计算(即广义特征值问题进行求解)。大规模结构总自由度数N很大,大规模结构拓扑优化过程中,广义特征值问题的系数矩阵常出现病态现象。求解此类广义特征值问题的部分特征对的难度是巨大的,且传统对角预处理器可以看作是提取K的对角元素,忽略该未知元素(用这个对角元素表示)对其他未知元素可能存在的耦合作用。基于传统对角预处理器的结构拓扑优化设计的效率不高,且难于适应规模结构拓扑优化过程中含有病态系数矩阵的结构固有频率和振型计算。
[0003]另外,当采用传统的串行拓扑优化算法对结构进行优化设计时,组装整体刚度矩阵的效率直接受到单元数量的影响。随着单元数量的增加,基于整体刚度矩阵和质量矩阵求解广义特征值问题的部分特征对的时间急剧上升,每个优化迭代步需要几十分钟甚至更久。因此,借助最新发展的具有多线程并行计算的GPU,研发既能减少迭代次数,且又能适应含有病态系数矩阵的结构固有频率和振型并行计算的方法是很有必要的。基于GPU的大规模结构拓扑频率优化方法会极大地提高计算效率,满足大规模工程结构设计要求。

技术实现思路

[0004]本申请实施的目的是提供一种RDPRQCG的大规模结构拓扑频率优化方法。首先,提出了一种新的变缩比和对角预处理器;并将该预处理技术与GPU并行计算方式和瑞利商共轭梯度法结合,提出了一种基于GPU的RDPRQCG大规模结构拓扑频率优化方法和其并行计算技术,该方法和其技术显著地提高了结构拓扑频率优化效率,进而可以解决
技术介绍
中涉及的大规模结构拓扑频率优化问题。
[0005]为了解决上述技术问题,本申请是这样实现的:
[0006]本申请实施例提供了一种RDPRQCG的大规模结构拓扑频率优化方法,包括:
[0007]步骤S1、采用GPU与提出的变缩比和对角预处理的正交加速瑞利商共轭梯度算法进行并行计算,获得大规模结构少量固有频率和振型;
[0008]步骤S2、基于并行计算获得的大规模结构少量固有频率和振型,以结构柔顺度最小为目标函数,同时考虑频率和体积约束,给出大规模结构拓扑优化模型、并行求解方法和典型案例。
[0009]可选的,所述采用GPU与提出的RDPRQCG算法进行并行计算,获得大规模结构少量固有频率和振型,包括:
[0010]步骤S11、依据结构总刚度矩阵K,采用式(1)构建新的对角线预处理矩阵的对角
元素
[0011][0012]其中,是总体坐标系中扩维的结构单元刚度矩阵K
(e)
的第i行第j列元素,m为结构单元个数,N为结构自由度数;
[0013]步骤S12、采用式(1)构建由下列对角正定矩阵定义的对角线预处理矩阵
[0014][0015]其中,
[0016][0017]其中,δ是一个小量,取δ=10
‑6;
[0018]步骤S13、在第k次迭代步,引入一个缩放因子γ
k
,构建下列变缩比和对角预处理矩阵W
k

[0019][0020]式中,为变缩比,为变缩比,
[0021]步骤S14、设定计划求广义特征值问题的前p个最小特征值及其对应的特征向量,该最小特征值即为固有频率,该特征向量即为固有振型,包括:
[0022]选取N
×
p阶初始迭代向量且满足其中I
p
是一个p阶的单位矩阵;设定容许误差TOLL,最大迭代次数NITMAX,重新启动步数NREST和平均残余量然后由式(5)计算初始残量矩阵:
[0023][0024]其中,M为质量矩阵,G
(0)
表示初始残差矩阵,p为p个初始残差向量,
[0025]令迭代步k=0;若k比NITMAX小,且比TOLL大,则执行迭代步骤S15

S17,否则退出迭代,执行步骤S16;
[0026]步骤S15、根据当前迭代步数判断k是否进行Ritz投影步,Ritz投影步相当于小规模广义特征值问题,即:
[0027][0028]式中,y是一个p维列向量,且X
(k)
满足
[0029]求解小规模广义特征值问题,得到p个特征向量y
i
(i=1,

,p),以及其对应的特征值μ
i
(i=1,

,p);
[0030]将各个y
i
(i=1,

,p)单位化,得:Q={q1,

,q
p
},则矩阵Q满足:
[0031][0032]令则是一个对角阵,且满足由式(8)和(9)求得P
(k)
,即:
[0033][0034][0035]步骤S16、新的矩阵由式(10)计算求得:
[0036][0037]对于j=1,

,p,执行以下步骤:
[0038]将向量和关于矩阵M正交化得到具体方法如下:
[0039]若j=1,则
[0040]若j≠1,则由式(11)求出γ
i

[0041][0042]通过γ
i
和式(12)求出即:
[0043][0044]通过最小化瑞利商求得α
j
的值;
[0045]其中,α
j
为标量,其值由方程(bc

ad)t2+(nd

mb)t+ma

nc=0确定,其中,
[0046][0047]其解为:
[0048][0049]式中,δ=(nd

mb)2‑
4(bc

ad)(ma

nc);
[0050]由式(14)求得为:
[0051][0052]将向量和关于矩阵M正交规范化得到向量具体方法如下:
[0053]若则直接进行正交化处理;否则,先令直接进行正交化处理;
[0054]步骤S17、将X
(k
+1)代入式(8)计算余量矩阵G
(k+1)

[0055]将代入式(15)

式(17)计算平均残余量具体方法如下:
[0056][0057][0058][0059]令k=k+1,进行迭代;
[0060]收敛判断,若比TOLL小,停止迭代;则则是符合要求的前p个特征值的本文档来自技高网
...

【技术保护点】

【技术特征摘要】
1.一种RDPRQCG的大规模结构拓扑频率优化方法,其特征在于,包括:步骤S1、采用GPU与提出的变缩比和对角预处理的正交加速瑞利商共轭梯度算法进行并行计算,获得大规模结构少量固有频率和振型;步骤S2、基于并行计算获得的大规模结构少量固有频率和振型,以结构柔顺度最小为目标函数,同时考虑频率和体积约束,给出大规模结构拓扑优化模型、并行求解方法和典型案例。2.根据权利要求1所述的方法,其特征在于,所述采用GPU与提出的变缩比和对角预处理的正交加速瑞利商共轭梯度算法进行并行计算,获得大规模结构少量固有频率和振型,包括:步骤S11、依据结构总刚度矩阵K,采用式(1)构建新的对角线预处理矩阵的对角元素的对角元素其中,是总体坐标系中扩维的结构单元刚度矩阵K
(e)
的第i行第j列元素,m为结构单元个数,N为结构自由度数;步骤S12、采用式(1)构建由下列对角正定矩阵定义的对角线预处理矩阵步骤S12、采用式(1)构建由下列对角正定矩阵定义的对角线预处理矩阵其中,其中,δ是一个小量,取δ=10
‑6;步骤S13、在第k次迭代步,引入一个缩放因子γ
k
,构建下列变缩比和对角预处理矩阵W
k
:式中,为变缩比,为变缩比,步骤S14、设定计划求广义特征值问题的前p个最小特征值及其对应的特征向量,该最小特征值即为固有频率,该特征向量即为固有振型,包括:选取N
×
p阶初始迭代向量且满足其中I
p
是一个p阶的单位矩阵;设定容许误差TOLL,最大迭代次数NITMAX,重新启动步数NREST和平均残余量然后由式(5)计算初始残量矩阵:其中,M为质量矩阵,G
(0)
表示初始残差矩阵,为p个初始残差向量,
令迭代步k=0;若k比NITMAX小,且比TOLL大,则执行迭代步骤S15

S17,否则退出迭代,执行步骤S16;步骤S15、根据当前迭代步数判断k是否进行Ritz投影步,Ritz投影步相当于小规模广义特征值问题,即:式中,y是一个p维列向量,且X
(k)
满足求解小规模广义特征值问题,得到p个特征向量y
i
(i=1,

,p),以及其对应的特征值μ
i
(i=1,

,p);将各个y
i
(i=1,

,p)单位化,得:Q={q1,

,q
p
},则矩阵Q满足:令则是一个对角阵,且满足由式(8)和(9)求得P
(k)
,即:,即:步骤S16、新的矩阵由式(10)计算求得:对于j=1,

,p,执行以下步骤:将向量和关于矩阵M正交化得到具体方法如下:若j=1,则若j≠1,则由式(11)求出γ
i
;通过γ
i
和式(12)求出即:通过最小化瑞利商求得α
j
的值;其中,α
j
为标量,其值由方程(bc

ad)t2...

【专利技术属性】
技术研发人员:周泉钟宇陈仙阁李水生宋英杰荣见华蔡金虎赵磊冯俊杰
申请(专利权)人:中国建筑第五工程局有限公司
类型:发明
国别省市:

网友询问留言 已有0条评论
  • 还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。

1