利索能及
我要发布
收藏
专利号: 2021104386576
申请人: 湖南科技大学
专利类型:发明专利
专利状态:已下证
更新日期:2026-06-16
缴费截止日期: 暂无
联系人

摘要:

权利要求书:

1.一种材料多尺度双向演化结构优化算法,其特征在于,具体包括如下步骤:S1:使用有限元网格离散整个设计域;

S2:定义初始参数,包括目标体积率Vobj、各级材料的平均值调整系数αj、升级率er与降级率rr、材料尺度数n和各材料的弹性模量物理参数;

S3:进行有限元分析;

S4:计算出所有单元的灵敏度;

S5:进行灵敏度的过滤与更新:S6:计算每种材料的升级指标Pj,1和降级指标Pj,2;

S7:进行相邻材料等级间的升降级操作;

S8:判断是否达到目标体积率Vobj,是则优化结束,否则回到S3。

2.根据权利要求1所述的材料多尺度双向演化结构优化算法,其特征在于,S2中,定义初始参数后,使每种材料通过各异的弹性模量赋值来体现等级差别,目标函数为最小化柔顺度并施加体积约束,数学模型描述为:式中C为结构柔顺度指标;K、u、P分别为结构总体刚度矩阵、位移向量及荷载向量;Ei、ui、E0、ki分别为第i个单元的弹性模量、第i个单元的位移向量、基准弹性模量,以基准弹性模量为基础确定其他材料的弹性模量、单元刚度矩阵;V、Vobj分别为结构当前体积率及目标体积率。

3.根据权利要求2所述的材料多尺度双向演化结构优化算法,其特征在于,S4中,单元灵敏度的计算方法为:

结构的刚度通过平均应变能间接表示,其公式为:不同材料的弹性模量为E0的线性内值,其表达式为:Ei=E0xi                                  (3)式中xi表示第i个单元的单元状态,1表示最高等级材料,其余材料根据材料尺度在0‑1之间线性内插;

根据式(3),当某结构中某一个单元状态xi发生改变时,结构刚度的改变为:ΔK=K′‑K=‑ΔxiE0ki     (4)式中K′为单元状态改变后的刚度矩阵;

在优化中,结构的荷载向量是与单元状态无关的常量,因此单元状态变化时:KΔu+ΔKu+ΔKΔu=0         (5)忽略相对高阶的微量ΔKΔu,单元状态改变时结构位移的改变为:Δu=‑KΔKu       (6)根据式(2)和(4)可以得到,结构平均应变能的改变为:因此基于平均柔顺度的单元灵敏度可以被定义为:

4.根据权利要求3所述的材料多尺度双向演化结构优化算法,其特征在于,S6中,升降级指标Pj,1和Pj,2的计算方法为式中Pj,1和Pj,2分别指第j种材料的升级指标和降级指标,αj为第j种材料的平均值调整系数, 为第j种材料中第i个单元的灵敏度, 为第j种材料所有单元的灵敏度平均值。

5.根据权利要求4所述的材料多尺度双向演化结构优化算法,其特征在于,所述αj>1,或αj<1。

6.根据权利要求1所述的材料多尺度双向演化结构优化算法,其特征在于,S2中:定义初始参数还包括定义折算体积:

式中Vi为第i个单元的体积,Ei为该单元当前所属材料的弹性模量,E0为基准弹性模量。

7.根据权利要求4所述的材料多尺度双向演化结构优化算法,其特征在于,S7‑S8的具体步骤为:

S7.1:定义j=n;

S7.2:判断是否j>1,如果判断结果为否,则进行S7.3,如果判断结果为是,则进行S7.4;

S7.3:判断是否V=Vobj,如果判断结果为否,则返回至S3,如果判断结果为是,则结束算法;

S7.4:判断是否Pj,2>Pj‑1,1,如果判断结果为否,则对第j‑1种材料开展升级操作,令j=j‑1,返回S7.2;如果判断结果为是,则对第j种材料开展降级操作,令j=j‑1,返回S7.2。