应力约束下三维拓扑优化MATLAB程序

本MATLAB程序为考虑应力约束的3D拓扑优化程序，可为开展进一步研究提供基础和平台。

输入参数如下：

nelx=200;   %X单元数

nely=60;    %Y单元数

nelz=1;     %Z单元数

volfrac=0.3;%体分比

penal=3;    %密度惩罚因子

q=0.5;      %应力惩罚因子

p=10;       %凝聚函数参数

rmin=2.5;   %过滤半径

优化目标：凝聚应力最小；

灵敏度分析方法：链式法则

优化算法：MMA。

MMA算法可直接调用子程序mmasub和subsolv ，因版权原因，请向Krister Svanberg发邮件获取（krille@math.kth.se），人很nice，都会给的。网上也有相关资源。可私信我免费分享，作它用请向Krister Svanberg发邮件获取许可。

输出：最大应力变化曲线、结构3D显示、Von  Mises应力显示。

最大应力变化曲线

结构3D显示