本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顯示