结构拓扑优化程序

结构拓扑优化程序如果开发者对于本文件有需要的可以参考。

while change > 0.001: #在循环体中，首先执行FEA以获得目标函数和元素灵敏度ae;然后ae被过滤(如果存在有效的rmin> 0)
iter += 1 # 计数器用于计算当前迭代次数。
# Run FEA
oh.append(FEA(iter, mddb, xe, ae)) # oh.append 函数增加字符串
# Process sensitivities
if rmin > 0: fltAe(ae, fm)
if iter > 0: ae = dict([(k, (ae[k] + oae[k]) / 2.0) for k in ae.keys()])
oae = ae.copy()
# BESO optimization
vh.append(sum(xe.values()) / len(xe))
nv = max(vf, vh[-1] * (1.0 - ert)) #由公式10计算
BESO(nv, xe, ae, part, elmts)
if iter > 10: change = math.fabs((sum(oh[iter - 4:iter + 1]) - sum(oh[iter - 9:iter - 4])) / sum(oh[iter - 9:iter - 4]))
# 结果保存，历史信息列表oh和vh存储在customData中，customData保存在名为“Final_design.cae”的最终设计模型数据库，可以通过ABAQUS打开。