[发明专利]一种基于凸模型不确定性的连续体结构非概率可靠性拓扑优化方法

权利要求书

1.一种基于凸模型不确定性的连续体结构非概率可靠性拓扑优化方法，其特征在于：实现步骤如下：

第一步：以结构的体积最小即重量最小作为优化目标，以结构关心的重要位置的位移作为约束条件，采用带惩罚因子的固体各向同性材料插值模型(Solid IsotropicMaterial with Penalization，简称SIMP)，建立拓扑优化模型如下：

优化目标:

约束条件:Ku＝F

u_i≤u_c,i＝1,2,…,NC

0＜ρ_min≤ρ_e≤1

其中，Vol是优化区域的体积，ρ＝[ρ₁,ρ₂,...,ρ_NE]^T是相对密度向量，ρ_e和V_e分别为第e个单元的相对密度和体积，NE为优化区域划分的单元总数，ρ_min为单元相对密度的下限；K为单元的总体刚度矩阵，u为单元的总体位移列向量，F为总体载荷列向量；u_i是第i个位移约束点的实际位移值，u_c是第i个位移约束的容许位移值，NC为位移约束的个数，对于SIMP模型,单元的弹性模量是材料相对密度的函数：

其中P＞1是惩罚因子，用于实现对中间密度单元的惩罚，按照经验，一般取P＝3，E₀是完全实心材料的弹性模量；

第二步：以凸模型描述优化模型的不确定性参数，通过静力平衡方程K(b)u(b)＝F(b)以及凸集合理论方法可以得到设计域节点位移响应，进而建立优化模型的非概率可靠性指标，可建立如下非概率拓扑优化模型：

优化目标:

约束条件:K(b)u(b)＝F(b)

R_i≥R_c,i＝1,2,…,NC

0＜ρ_min≤ρ_e≤1

其中，R_i凸模型可靠度，

第三步：在原非概率可靠性指标的基础上建立可靠性位移差指标，从而改善原非概率可靠度指标对设计变量求灵敏度的缺陷，计算相应的可靠性位移差，利用可靠性位移差可以将原非概率可靠性优化模型修改如下：

优化目标:

约束条件:K(b)u(b)＝F(b)

0＜ρ_min≤ρ_e≤1

b∈Ω_C＝{b|(b-b^c)^TW(b^r,ρ_b)(b-b^c)≤θ²}

其中，为可靠性位移差，b为m维不确定向量，b^c为椭球中心值向量，b^r椭球半长短轴向量，W(b^r,ρ_b)∈R^m×m是一个对称正定矩阵，由向量b^r和相关系数矩阵ρ_b＝(ρ(b_j1,b_j2))决定；参数θ量化了向量b的不确定程度，d_i表示从坐标原点到失效超平面的距离，d_i(R_c)是从原点到给定的非概率可靠度R_c确定的极限状态超平面的距离；

第四步：根据第三步得到的可靠性位移差的表达式，通过复合函数求导法则对设计变量进行求导，并运用伴随向量法和有限差分法得到可靠性位移差对设计变量的灵敏度；

第五步：利用第三步中得到的可靠性位移差以及第四步中得到的可靠性位移差对设计变量的灵敏度代入MMA算法中对非概率可靠性优化模型进行求解，得到新的设计变量；

第六步：判定新的设计变量是否满足收敛性条件，若不满足收敛性条件，则将已经完成的迭代次数增加1，并返回第二步；否则，迭代过程结束。

2.根据权利要求1所述的一种基于凸模型不确定性的连续体结构非概率可靠性拓扑优化方法，其特征在于：以凸模型来描述参数的不确定性，以结构的体积最小即重量最小作为优化目标，以结构关心的重要位置的位移作为约束条件，采用带惩罚因子的固体各向同性材料插值模型建立拓扑优化模型。