基于载荷不确定性的材料结构多尺度拓扑优化设计方法与流程

[0001]
本发明涉及材料结构多尺度优化设计相关技术领域，特别是涉及基于载荷不确定性的材料结构多尺度拓扑优化设计方法。


背景技术：

[0002]
为了满足现代工业产品对“结构超轻量化”、“功能特殊化”与“性能集成化”的需求，同时考虑材料微观结构优化和宏观结构优化进行材料结构多尺度拓扑优化设计，获取最优的材料微观结构和宏观结构多尺度设计，获得周期性多空复合材料结构；它具有高比刚度、高比强度、吸能等优异的力学性能，广泛应用于航空航天、汽车工业、建筑行业等相关领域。
[0003]
现有的材料结构多尺度拓扑优化设计方法是考虑确定性载荷工况，即载荷的大小与方向均为确定的，但在实际工程中不可避免存在载荷不确定性因素。因此，目前急需一种基于载荷不确定的材料结构多尺度拓扑优化设计方法，以解决上述问题。


技术实现要素：

[0004]
本发明的目的是提供基于载荷不确定的材料结构一体化拓扑优化方法，采用一阶摄动法对不确定性载荷进行量化，有效地减低统计性指标计算的复杂程度，提高计算效率，实现了宏观结构和材料微观结构的多尺度化稳健性拓扑优化设计，获得的周期性多空复合材料结构具有更好的稳健性。
[0005]
本发明提供基于载荷不确定性的材料结构多尺度拓扑优化设计方法，包括如下步骤：
[0006]
s1、定义设计域，给定初始测试应变条件，设定宏观结构与材料微观结构的尺寸、边界、材料参数和不确定性载荷参数，设置设计变量的过滤半径，设置迭代准则，初始化宏观结构与材料微观结构的设计变量；
[0007]
s2、对微观结构边界添加周期性边界条件，在给定的初始测试应变条件下进行微观结构有限元分析，计算材料微观结构的宏观等效弹性张量；
[0008]
s3、采用摄动方法对不确定性载荷的参数进行量化，采用随机有限元法得到宏观结构的单元节点位移矢量，计算宏观结构柔顺度的期望和方差；
[0009]
s4、将宏观结构柔顺度的期望值与方差进行加权得到加权和，以加权和最小化作为优化目标，以宏观结构和材料微观结构的体积份数作为约束，采用变密度法与能量均匀化方法构建考虑载荷不确定性的材料结构多尺度拓扑优化模型；
[0010]
s5、根据所述材料结构多尺度拓扑优化模型，计算宏观结构与材料微观结构优化设计的目标函数和约束灵敏度，采用heaviside映射过滤方法对目标函数和约束灵敏度进行修正，获得宏观结构和材料微观结构拓扑构型；
[0011]
s6、采用优化准则法oc更新宏观结构与材料微观结构设计变量，判断是否满足收敛准则；若不满足，转至步骤s2，若满足，输出最优的宏观结构与材料微观结构拓扑型构；
[0012]
优选地，所述s2包括以下步骤：
[0013]
s2-1、在微观结构边界添加周期性边界条件，将宏观位移场和微观材料内的周期性波动场相加得到材料微观结构的位移场；
[0014]
s2-2、在给定的初始测试应变条件下，将位移场分为四个部分，建立材料微观结构的有限元平衡方程进行有限元分析，并采用能量均匀化方法求得材料微观结构的宏观等效弹性张量；
[0015]
优选地，所述s3包括以下步骤：
[0016]
s3-1、采用摄动法，以一阶泰勒展开式表示具有不确定性大小向量和方向向量的作用载荷向量，对不确定性载荷的参数进行量化，得到不确定性量化参数；
[0017]
s3-2、采用固体各向同性材料惩罚模型simp表示宏观结构的材料分布，由材料微观结构的宏观等效弹性张量建立宏观结构有限元平衡方程，采用随机有限元法，通过一阶泰勒展开式求得宏观结构的单元节点位移矢量；
[0018]
s3-3、将不确定性量化参数和高阶函数的平均值代入结构柔顺度的不确定性表达式中，求得柔顺度的期望e(c)和方差var(c)；
[0019]
优选地，所述s4构建基于载荷不确定的材料结构多尺度拓扑优化模型为如下方法：
[0020]
find{x y}
[0021][0022]
s.t.k(x,y)u
i
＝f
i
(i＝0,1,2)
[0023][0024][0025][0026]
0≤y
j
≤1 j＝1,...,n
[0027]
为宏观结构的柔顺度，x和y分别为宏观设计变量和微观设计变量，统一表示为单元密度，u
i
为宏观结构的位移场，f0为载荷列阵的均值矩阵，f1为不确定载荷大小对随机变量的偏导数，f2为不确定载荷方向对随机变量的偏导数；α为加权系数，k(x,y)为宏观结构的总刚度矩阵，ω0和y0分别为宏观和微观结构的初始体积，ω
e
和y
e
分布为宏观和微观结构的单元体积，和分别为宏观和微观结构允许的材料体积比，n和n分别为宏观和微观结构的有限元的数目；
[0028]
优选地，所述s5包括以下步骤：
[0029]
s5-1、采用拉格朗日乘子法计算宏观结构柔顺度相对宏观设计变量的灵敏度，计算柔顺度对微观结构设计变量的灵敏度，并且求解宏观结构和材料微观结构体积约束灵敏度；
[0030]
s5-2、采用基于heaviside函数的映射过滤方法对目标函数和宏观结构材料微观结构约束灵敏度进行修正，首先采用密度过滤方法对单元密度进行过滤，通过相邻单元密
度的加权平均方法，获得过滤后的单元密度；
[0031]
s5-3、采用heaviside映射过滤法对过滤后的单元密度进行修正，使单元密度值向0和1两端集中，获得宏观结构和材料微观结构拓扑构型；
[0032]
优选地，所述s6采用优化准则法oc算法更新宏观结构与材料微观结构设计变量，求解拓扑优化问题，判断是否满足收敛准则，若不满足，转至步骤s2，若满足，则材料结构多尺度拓扑优化过程结束，输出最优的宏观结构和材料微观结构拓扑构型。
[0033]
本发明公开了以下技术效果：
[0034]
(1)本发明基于载荷不确定性的材料结构多尺度拓扑优化设计，采用一阶摄动法对不确定性载荷进行量化；采用变密度法与能量均匀化方法构建基于载荷不确定性的材料结构多尺度拓扑优化模型；采用伴随矩阵技术进行目标函数和约束灵敏度分析；获得最优的宏观结构拓扑和材料微观结构拓扑；
[0035]
(2)采用一阶摄动法对不确定性载荷进行量化，计算考虑载荷大小和方向不确定条件下的柔顺度期望值和方差表达式，有效地减低统计性指标计算的复杂程度，提高计算效率；
[0036]
(3)考虑载荷的方向和大小的不确定，进行宏观结构和材料微观结构的多尺度化稳健性拓扑优化设计，获得的周期性多空复合材料结构，具有更好的稳健性。
附图说明
[0037]
图1为本发明基于载荷不确定性的材料结构多尺度拓扑优化设计方法流程示意图；
[0038]
图2为载荷确定性条件下的材料结构多尺度拓扑优化设计获得的二维立柱结构设计域图；
[0039]
图3为载荷确定性条件下的材料结构多尺度拓扑优化设计获得的立柱的宏观结构拓扑构型图；
[0040]
图4为载荷确定性条件下的材料结构多尺度拓扑优化设计获得的立柱的材料微观结构拓扑图；
[0041]
图5为载荷不确定性的材料结构多尺度拓扑优化设计获得的立柱的3
×
3微观结构图；
[0042]
图6为载荷不确定性的材料结构多尺度拓扑优化设计获得的立柱的宏观结构拓扑构型；
[0043]
图7为载荷不确定性的材料结构多尺度拓扑优化设计获得的材料微观结构拓扑；
[0044]
图8为载荷不确定性的材料结构多尺度拓扑优化设计获得的3
×
3微观结构图。
具体实施方式
[0045]
下面将结合本实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有付出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
[0046]
参照图1所示,本发明提供基于载荷不确定性的材料结构多尺度拓扑优化设计方
法，包括如下步骤：
[0047]
s1、定义设计域，给定初始测试应变条件，设定宏观结构与材料微观结构的尺寸、边界条件、材料参数和不确定载荷的大小和方向参数，设置设计变量的过滤半径，设置迭代准则，初始化宏观结构与材料微观结构的设计变量；
[0048]
s2、对微观结构边界添加周期性边界，给定的初始测试应变条件下进行微观结构有限元分析，计算材料微观结构的宏观等效弹性张量；
[0049]
s2-1、在微结构边界施加周期性边界，将宏观位移场和微观材料内的周期性波动场相加得到微观结构的位移场：
[0050][0051]
u
i
为微观结构的位移场，i,j表示索引向量值，为给定的初始测试应变，为微观材料内的周期性波动位移场，y
j
为微观结构的局部坐标，
[0052]
求得材料微观结构所有法线方向上边界的位移场为：
[0053][0054]
k+与k-表示微观结构内一对相反的两个边界的法线方向,由上式相减消除材料微观结构内部的周期性波动位移场，可得：
[0055][0056]
s2-2、将位移场全局位移分为四个部分，建立材料微观结构的有限元平衡方程，进行有限元分析：
[0057][0058]
k
kl
(k,l＝1,2,3,4)表示材料微观结构刚度矩阵分块矩阵，表示给定的位移，χ2表示微观结构的内部节点的未知位移，χ3和χ4表示位于相反边界的节点的未知位移，为材料微观结构相反边界的位移差值，且满足其中f1为给定位移节点出的反力，f2＝0，由周期假设有f3+f4＝0,且k
ij
＝k
ji
，则上式得：
[0059][0060]
采用能量均匀化方法求得材料微观结构的宏观等效弹性张量：
[0061][0062]
为材料微观结构的宏观等效弹性张量，y0表示材料微观结构的面积，y
j
表示微观材料结构有限元单元的面积，dy表示材料结构中的微元体积，t表示矩阵转置的符号，
χ
0(ij)
表示由施加在微观结构上的初始测试应变引起的位移，χ
ij
表示为未知的位移，由初始测试应变引起的微观结构内部产生的位移；b为应变位移矩阵，c(y)为基材料的弹性矩阵，n为材料微观结构的有限元的数目；
[0063]
s3、考虑载荷大小和方向不确定性，采用摄动方法对不确定性载荷的参数进行量化，采用随机有限元法求解宏观结构的单元节点位移矢量，计算宏观结构柔顺度的期望和方差；
[0064]
s3-1、采用摄动法，将一阶泰勒展开式表示具有不确定性大小向量μ和方向向量θ的作用载荷向量f：
[0065]
f(μ,θ)≈f0+f1·
δμ+f2·
δθ
[0066]
δμ和δθ分别为载荷大小和方向随机变化量称为载荷不确定参数，f0为载荷列阵的均值矩阵，f1为不确定载荷大小对随机变量的偏导数，f2为不确定载荷方向对随机变量的偏导数；
[0067]
s3-2、采用simp(solid isotropic material with penalization，固体各向同性材料惩罚模型)表示宏观结构的材料分布，并建立宏观结构有限元平衡方程：
[0068]
ku(μ,θ)＝f(μ,θ)
[0069]
k为宏观结构的总装刚度矩阵，u为宏观结构的单元节点位移矢量，宏观结构的节点位移矢量通过一阶泰勒展开式求得为：
[0070]
u(μ,θ)＝u0+u1·
δμ+u2·
δθ
[0071]
u0为位移的均值矢量，u1和u2分别位移矢量对载荷大小和方向随机变量的偏导数，由下式求得：
[0072]
u
i
＝k-1
f
i
(i＝0,1,2)；
[0073]
s3-3、将载荷不确定参数δμ和δθ高阶函数的平均值代入结构柔顺度c的不确定性表达式中，求得柔顺度的期望e(c)和方差var(c)：
[0074][0075]
e(δμ)2为载荷大小的高阶函数的均值，e(δθ)2为载荷方向的高阶函数的均值，u2为位移矢量对随机变量的偏导数：
[0076][0077]
e(c2)为柔顺度平方的均值，e(c)2柔顺度均值的平方，var(δμ)为不确定载荷大小的方差，var(δθ)为不确定载荷方向的方差，var(δμδθ)为不确定载荷大小与方向乘积的方差；
[0078]
s4、以结构柔顺度的期望值与标准差进行加权得到加权和，以加权和最小化作为优化目标，以宏观结构和材料微观结构的体积份数作为约束，采用变密度法与能量均匀化方法构建基于载荷不确定的材料结构多尺度拓扑优化模型，构建基于载荷不确定的材料结构多尺度拓扑优化模型方法如下：
[0079]
find{x y}
[0080][0081]
s.t.k(x,y)u
i
＝f
i
(i＝0,1,2)
[0082][0083][0084][0085]
0≤y
j
≤1 j＝1,...,n，
[0086]
为宏观结构的柔顺度，x和y分别为宏观设计变量和微观设计变量，u
i
为宏观结构的位移场，α为加权系数，k(x,y)为宏观结构的总刚度矩阵，ω0和y0分别为宏观和微观结构的初始体积，ω
e
和y
e
分别为宏观和微观结构的单元体积，和分别为宏观和微观结构允许的材料体积比，n和n分别为宏观和微观结构的有限元的数目；
[0087]
s5、根据所述材料结构多尺度拓扑优化模型，计算宏观结构与材料微观结构优化设计的目标函数和约束灵敏度，采用heaviside映射过滤方法对目标函数和约束灵敏度进行修正，获得清晰的宏观结构和材料微观结构拓扑构型；
[0088]
s5-1、引入拉格朗日乘子计算宏观结构柔顺度约束灵敏度：
[0089][0090]
α为加权系数，λ
i
(i＝0,1,2)表示为拉格朗日乘子，结构柔顺度对宏观设计变量的约束灵敏度求得为：
[0091][0092]
其中a0＝f
0t
+f
t1
·
e(δμ)+f
2t
·
e(δθ)，
[0093]
a1＝f
0t
·
e(δμ)+f
t1
·
e(δμ)2+f
2t
·
e(δθ)e(δμ)，
[0094]
a2＝f
0t
·
e(δθ)+f
t1
·
e(δμ)e(δθ)+f
2t
·
e(δθ)2，
[0095]
b0＝f1·
var(δμ)
·
f
1t
+f2·
var(δθ)
·
f
2t
，
[0096]
b1＝f0·
var(δμ)
·
f
0t
+f1·
var(δμ)2·
f
1t
+f2·
var(δμδθ)
·
f
2t
，
[0097]
b2＝f0·
var(δθ)
·
f
0t
+f1·
var(δμδθ)
·
f
1t
+f2·
var(δθ)2·
f
2t
，
[0098]
为了消除项，令下列等式成立：
[0099][0100]
结构柔顺度的约束灵敏度公式为：
[0101][0102]
将拉格朗日乘子λ
i
代入结构柔顺度的约束灵敏度公式得：
[0103][0104]
其中
[0105][0106]
p是惩罚因子，e0是材料的弹性模量，e
min
空洞材料的弹性模量，为充满材料的单元刚度矩阵与弹性模量的比值；
[0107]
宏观体积约束对宏观设计变量的约束灵敏度为：
[0108][0109]
结构柔顺度对微观结构设计变量的约束灵敏度为：
[0110][0111]
其中
[0112][0113]
微观体积约束对微观设计变量的约束灵敏度为：
[0114][0115]
5-2、采用基于heaviside函数的映射过滤方法对目标函数和宏观结构、材料微观结构的约束灵敏度进行修正，首先采用密度过滤方法对单元密度进行过滤，通过相邻单元密度x
j
的加权平均，获得过滤后的单元密度x％
e
：
[0116][0117]
ω
e
表示与单元e质心距离小于过滤半径r
min
的所有单元的设计变量集合，w
j
为权重系数，其表达式为
[0118]
w(x
j
)＝r-|x
j-x
i
|
[0119]
r为最小过滤半径，x
j
和x
i
分别表示单元j和单元i的坐标；
[0120]
5-3、采用密度过滤方法的拓扑优化结果中仍存在灰度单元，采用heaviside映射
过滤法对过滤后的单元密度进行修正，使单元的密度值向0和1两端集中，得到修正后的单元密度，获取清晰的宏观结构和材料微观结构拓扑构型；
[0121][0122]
为单元物质密度；β为控制光滑映射的参数；η为阈值，取值为0.5；
[0123]
s6、采用优化准则法oc算法更新宏观结构与材料微观结构设计变量，求解拓扑优化问题，迭代至满足收敛为止，获得最优的宏观结构拓扑和材料微观结构拓扑，获取宏观结构和材料微观结构的多尺度化稳健性拓扑优化设计，获得的周期性多空复合材料结构具有更好的稳健性，判断是否满足收敛准则，若不满足，转至步骤s2，若满足，则材料结构多尺度拓扑优化过程结束，输出最优的宏观结构和材料微观结构拓扑构型。
[0124]
为了进一步验证本发明基于载荷不确定性的材料结构多尺度拓扑优化设计方法的有效性，本实施例以二维立柱结构为例对本发明进行解释；
[0125]
立柱结构的设计域、边界条件如图2所示，宏观结构的设计域尺寸长
×
高为500mm
×
250mm，设计域下端为固定边界，载荷作用f在宏观结构上边界的中心处，所用材料的弹性模量e为1gpa，泊松比v为0.3，宏观设计域离散为50
×
50个平面四边形单元，材料微观结构的设计域长
×
高为0.01mm
×
0.01mm，材料微观结构的设计域50
×
50个平面四边形单元，载荷f的大小为1kn，作用方向与水平方向的夹角为-π/2(垂直向下)，载荷大小标准差为0，载荷方向标准差为0.3。宏观结构和材料微观结构的允许的体积份数为25％，最小过滤半径r
min
为2倍单元尺寸，惩罚系数p为3。
[0126]
图3、图4和图5分别为载荷确定性条件下的材料结构多尺度拓扑优化设计获得的立柱的宏观结构拓扑构型、材料微观结构拓扑和3
×
3微观结构图，立柱宏观结构的材料主要分布在竖直载荷线附近，形成一种柱状结构，以承受宏观结构的竖直方向的承载，材料微观结构的材料主要分布在两根竖直柱状结构。图6、图7和图8分别为基于载荷不确定性的材料结构多尺度拓扑优化设计获得的立柱的宏观结构拓扑构型、材料微观结构拓扑和3
×
3微观结构图。与确定性载荷条件下的材料结构多尺度拓扑优化结果相比，基于载荷不确定性的材料结构多尺度拓扑优化获得的立柱宏观结构底部出现的一定角度的开口，以此来提高立柱结构的水平方向的承载能力；基于载荷不确定性的材料结构多尺度拓扑优化获得的材料微观结构拓扑构型形成更多网格结构，材料分布会更加均匀，并非集中在竖直的两根柱状结构附近，这种材料微观结构使得立柱宏观结构具有更好水平方向的承载能力。
[0127]
以上所述的实施例仅是对本发明的优选方式进行描述，并非对本发明的范围进行限定，在不脱离本发明设计精神的前提下，本领域普通技术人员对本发明的技术方案做出的各种变形和改进，均应落入本发明权利要求书确定的保护范围内。