1.一种考虑增材制造倾角约束的结构拓扑优化方法,其包括以下步骤:步骤一:定义初始条件即结构设计域、载荷条件、约束条件及材料的相关属性,将材料均匀的离散到每一个单元,得到初始单元密度分布;
步骤二:按照增材制造倾角约束的单元密度更新迭代规则得到满足倾角约束的可打印密度分布;
步骤三:对结构进行有限元分析,获取结构的整体位移场;
步骤四:构建考虑增材制造倾角约束的拓扑优化数学模型,计算得到目标函数值和关于设计变量的敏度;
步骤五:采用OC算法(Modified Optimization Criteria methods)更新设计变量,得到最新的密度分布;
步骤六:根据收敛条件判断目标函数是否收敛,如果满足收敛条件,则输出拓扑构型,结束迭代计算,否则,转入步骤二继续进行迭代计算;
所述考虑增材制造倾角约束的单元密度更新迭代规则,数学模型为:其中,i,j分别为单元在竖直方向和水平方向的位置,ξ(i,j)为每次迭代后(i,j)处单元的可打印密度值;x(i,j)为每次迭代前(i,j)处单元的初始密度值;Γ(i,j)为(i,j)处单元支撑区域内各单元可打印密度的最大值;Ψ(i,j)为(i,j)处单元支撑区域内各层最大可打印密度中的最小值,n为支撑区域的层数。
2.根据权利要求1所述的一种考虑增材制造倾角约束的结构拓扑优化方法,其特征在于:支撑区域选择一层得到的倾角约束为45度,其迭代模型为:ξ(i,j)=min(x(i,j),Γ(i,j))(2)Γ(i,j)=max(ξ(i‑1,j‑1),ξ(i‑1,j),ξ(i‑1,j+1))(3)。
3.根据权利要求2所述的一种考虑增材制造倾角约束的结构拓扑优化方法,其特征在于:支撑区域选择两层得到的倾角约束为63度,其迭代模型为:ξ(i,j)=min(x(i,j),Ψ(i,j))(4)
1 2
Ψ(i,j)=min(Γ(i,j),Γ(i,j)) (5)1
Γ(i,j)=max(ξ(i‑1,j‑1),ξ(i‑1,j),ξ(i‑1,j+1)) (6)2
Γ(i,j)=max(ξ(i‑2,j‑1),ξ(i‑2,j),ξ(i‑2,j+1)) (7)n
其中,Γ(i,j)为(i,j)处单元支撑区域内i‑n层中各单元可打印密度的最大值。
4.根据权利要求3所述的一种考虑增材制造倾角约束的结构拓扑优化方法,其特征在于:选取支撑区域为两层且每一层的单元数为三个,得到63度倾角约束的迭代模型为:
1 2
其中,x是单元初始密度,Ψ是单元支撑区域内各层最大可打印密度的较小值,Γi、Γi是每层中各单元可打印密度的最大值,ξ1、ξ2、ξ3、 是可打印结构密度,ns是单元支撑区域每一层的单元数,ns=3;ε和P是调节光滑近似模型的参数,ξ0取0.5。
5.根据权利要求4所述的一种考虑增材制造倾角约束的结构拓扑优化方法,其特征在于:考虑增材制造倾角约束的模型为:ξ(i,j)=smin(x(i,j),Ψ(i,j)) (12)
1 2
Ψ(i,j)=smin(Γ(i,j),Γ(i‑1,j)) (13)1
Γ(i,j)=smax(ξ(i‑1,j‑1),ξ(i‑1,j),ξ(i‑1,j+1)) (14)2
Γ(i,j)=smax(ξ(i‑2,j‑1),ξ(i‑2,j),ξ(i‑2,j+1)) (15)其中,smin和smax分别采用式(8)和(10)的近似形式,需要指出的是第i层的可打印密度与第1层到第i‑1层的可打印密度均密切相关,表示smin符号的简写形式,单一下标表示层数的索引符号。
6.根据权利要求5所述的一种考虑增材制造倾角约束的结构拓扑优化方法,其特征在于:考虑增材制造倾角约束的拓扑优化数学模型为:其中,c为目标函数,即结构的整体柔度值;K为整体刚度矩阵;U为整体位移;ue为单元位移矩阵;ke为单元刚度矩阵;F为外载荷矢量;V(x)为结构优化后的体积;V0为结构的初始体积;f为给定体积分数;ξi为结构的倾角约束等式即增材制造的倾角约束,表示smin符号的简写形式,单一下标表示层数的索引符号。
7.根据权利要求6所述的一种考虑增材制造倾角约束的结构拓扑优化方法,其特征在于:目标函数关于设计变量的敏度,可由以下公式进行推导:其中,λi为拉格朗日乘子矢量,对于第一层和第二层,定义 且 则I为单位矩阵。
8.根据权利要求7所述的一种考虑增材制造倾角约束的结构拓扑优化方法,其特征在于:采用OC算法更新设计变量:new
其中,m为一个正的移动界限, 为尺度因子,η为数值阻尼系数,xmin为取值下限;xe为经过优化准则算子迭代优化后的更新解。