
STORX — Shape and Topology Optimization for Research and Experimentation — ni mfumo huria wa MATLAB unaokusudia kuunganisha uboreshaji wa umbo kwa vigezo, uboreshaji wa umbo wa level-set, na familia kadhaa za uboreshaji wa topolojia ndani ya usanifu mmoja wa programu unaolenga vitu. Mchango wake mkuu si algoriti moja mpya, bali ni usanifu unaotenganisha jiometri, uchanganuzi wa elementi kikomo, mlinganyo wa hali, functional za lengo na vizuizi, uchanganuzi wa sensitivity, sasisho la design na vizuizi vya utengenezaji kupitia interface za msingi za abstract.
Ndani ya msingi huo huo, mfumo unaonyesha parametric shape optimization, Hamilton–Jacobi level-set, density optimization kwa SIMP na RAMP, ESO, BESO na PareTO tracing. Pia unaonyesha multiple load scenarios, self-weight, steady-state heat conduction, stress minimization, local volume fraction na fluid topology optimization inayotumia steady laminar incompressible Navier–Stokes pamoja na Brinkman penalization.
Matokeo hayasemi kuwa STORX ndiyo code ya MATLAB yenye kasi zaidi. Ikilinganishwa na top88, overhead ya usanifu wake wa jumla hufanya runtime kwa iteration kuwa karibu 65–68% ndefu katika resolutions za kati; kwenye 320×160 tofauti huwa karibu 60%. Kwa takriban elementi 50,000, iteration moja bila visualization huchukua karibu sekunde 0.4. Kwa hiyo hoja kuu ya utafiti ni modularity, transparency na extensibility, si ushindi wa kasi.
Tatizo linalolengwa
Katika shape optimization, design inaweza kudhibitiwa na radius ya tundu, nafasi ya notch, chamfer au boundary motion. Katika topology optimization, pseudo-density ya elementi, level-set field au topological sensitivity huamua wapi material ibaki na kama connectivity ibadilike. Ingawa representation hizi ni tofauti, zote hutumia mnyororo wa msingi: geometry → discretization → state solve → objective/constraint evaluation → sensitivity → regularization/filter → design update.
Usanifu wa object-oriented
STORX hutenganisha brep2d, gridMesher, fea2d, simulation2d, functional, mfgConstraints, parameterOpt2d, density2d, levelset2d, evolutionary2d na pareto2d. Functional mpya huhitaji hasa evaluate na gradient; manufacturing constraint mpya hutumia filterDesign na filterSensitivity.
B-Rep na FEA
Boundary Representation ya 2D hutumia vertices pamoja na line/arc segments kuunda domains zisizolazimika kuwa rectangular. Hii inaruhusu loads, supports na retained regions kufungamanishwa na boundaries za maana za jiometri.
Kwa mfano wa gripper: karibu elementi 4,000, \(E=2\) GPa, \(\nu=0.35\), nguvu 10 N. Triangular FEA inatoa \(1.00\times10^{-6}\) m na 6.87 MPa; grid FEA inatoa \(1.05\times10^{-6}\) m na 6.31 MPa.
Parametric shape optimization
Semi-analytic FD: \(p_0=(0.20,0.15,1.20,0.10)\) hadi \(p_{FD}=(0.363,0.224,1.500,0.076)\); iterations 118, FEA runs 1,070; area 1.906→1.7995 m²; compliance 3.91→3.95 N·m; maximum deflection karibu 0.04 mm.
GlobalSearch: \(p_{GS}=(0.310,0.226,1.248,0.141)\), area 1.79 m², compliance 4.022 N·m, FEA 305. MultiStart: \(p_{MS}=(0.327,0.138,1.466,0.155)\), area 1.798 m², compliance 4.001 N·m, FEA 2,545. Random Search: \(p_{RS}=(0.359,0.150,1.337,0.150)\), area 1.779 m², compliance 4.021 N·m, FEA 22.
Level-set shape optimization
Reinitialization hurejesha field karibu na signed-distance function ili \(|\nabla\psi|\approx1\) bila kubadilisha zero-level-set.
Density-based topology optimization
Framework inaruhusu SIMP au RAMP pamoja na OC, MMA au GCMMA. Continuation ya penalty inaweza kutumiwa kusukuma intermediate densities kuelekea solid/void.
Level-set TO na topological sensitivity
Sehemu ya modified HJE inaandika forcing term \(g=-sign(\psi)T\), lakini sentensi hiyo pia inataja positive weight factor \(w\). \(w\) haipo kwenye formula iliyoonyeshwa, hivyo tofauti hiyo ya chanzo haijasahihishwa kimya kimya.
ESO, BESO na PareTO
ESO huondoa material bila kuirudisha. BESO inaruhusu removal na re-introduction. PareTO inalenga local Pareto optimality kwa volume levels zinazofuatana na hutumia fixed-point iteration katika kila outer step.
Minimum feature size na physical-density projection
Filter ya minimum feature size hupunguza checkerboard na thin members. Heaviside projection huongeza ukali wa mabadiliko ya material; katika chanzo \(\eta=0.5\) na mifano ya \(\beta=1,8,64\) inaonyeshwa.
Retained regions
Gripper yenye elementi 10,000, PETG \(E=2\) GPa, \(\nu=0.35\), force 10 N na volume fraction 0.65: bila retain \(C/C_0=0.32\); kwa circular regions retained \(C/C_0=0.33\).
Multiple loads, self-weight na heat conduction
Cantilever yenye elementi 3,200, \(E=100\) GPa, \(\nu=0.3\), volume fraction 0.5; loads ni \(F_1=+200\) kN na \(F_2=-100\) kN. Average compliance: LSSO 2.11, standard HJE 1.84, modified HJE 1.68, SIMP 1.89, PareTO 1.70 N·m.
Bridge self-weight: elementi 3,000, \(E=100\) GPa, \(\nu=0.3\), density ya material 1000 kg/m³, \(g=[0,-10]\) m/s².
Mifano ya joto hutumia elementi 10,000 na volume fraction 0.5. Boundary heat flux: \(\kappa=1\), \(q_n=1\) W/m², \(T=0\) K, \(f=0\), SIMP p=3, OC. Internal heat generation: \(f=0.01\) W/m³, RAMP q=5, OC.
Stress minimization, local volume na fluid
Stress example: elementi 6,000, SIMP p=3, \(E=100\) GPa, \(\nu=0.3\), \(p_{vm}=6\); peak von Mises inapungua kutoka 60.6 MPa hadi 23.0 MPa.
Local-volume beam ina elementi 80,000, local fraction 0.5 na radius 6. Prose ya chanzo inasema p-norm p=6 lakini code ina localPNorm=16. Gripper local-volume ina elementi 80,000 active, local fraction 0.65, radius 10 na prose p=6.
Fluid wrapper: steady laminar incompressible Navier–Stokes + Brinkman. Pipe bend/double pipe hutumia elementi 2,000, SIMP p=3, MMA, \(V_f=0.3\), \(U_{in}=1\) m/s, \(\rho=1\) kg/m³, \(\mu=1\) kg/(m·s). Wind tunnel hutumia elementi 40,000, Re=10, active design rectangle yenye center [1.35,0.5] m, width 1.5 m na height 0.5 m. Objectives ni energy dissipation, drag minimization na lift maximization.
Umuhimu unaowezekana katika Afrika Mashariki
Utafiti haujafanywa kwa data ya Afrika Mashariki na hauonyeshi kuwa designs hizi zimehakikiwa moja kwa moja katika viwanda au vyuo vya eneo hilo. Hata hivyo framework huria ya MATLAB inaweza kuwa muhimu kwa elimu ya computational mechanics, mechanical/aerospace engineering, additive manufacturing na thermal-fluid design katika eneo hilo. Matumizi ya ndani yatahitaji kuthibitisha material properties, manufacturing tolerances, upatikanaji wa MATLAB, uwezo wa computation na experimental validation kulingana na mazingira ya taasisi au sekta husika.
Utafiti unasema nini na hauseni nini?
Unaonyesha kwamba familia nyingi za SO/TO zinaweza kuunganishwa katika architecture moja ya programu na kwamba new functionals/constraints zinaweza kuongezwa bila kuandika upya core optimizer. Hauseti kuwa STORX ni fastest solver, kwamba PareTO ni bora kwa kila tatizo, kwamba designs ni global optima, au kwamba benchmark za 2D ni uthibitisho wa moja kwa moja wa industrial-scale performance.
Mbinu na Matokeo ya Utafiti
| Mbinu / mfano | C (N·m) | δmax (m) | σvm,max (MPa) |
|---|---|---|---|
| LSSO Cantilever bottom | 7.85 | 9.37e-05 | 10.5 |
| LSSO Cantilever middle | 12.5 | 9.99e-05 | 4.93 |
| LSSO L-bracket top | 25.8 | 2.76e-04 | 37.8 |
| LSSO L-bracket mid | 25.7 | 2.89e-04 | 54.3 |
| LSSO MBB | 22.7 | 2.48e-04 | 6.18 |
| SIMP-OC Cantilever bottom | 6.98 | 7.46e-05 | 6.94 |
| SIMP-OC Cantilever middle | 10.6 | 8.50e-05 | 8.83 |
| SIMP-OC L-bracket top | 22.9 | 2.46e-04 | 21.0 |
| SIMP-OC L-bracket mid | 21.7 | 2.42e-04 | 23.0 |
| SIMP-OC MBB | 15.0 | 1.64e-04 | 6.98 |
| Standard HJE Cantilever bottom | 6.38 | 6.85e-05 | 2.34 |
| Standard HJE Cantilever middle | 9.68 | 7.84e-05 | 2.83 |
| Standard HJE L-bracket top | 21.2 | 2.27e-04 | 25.2 |
| Standard HJE L-bracket mid | 19.7 | 2.21e-04 | 36.1 |
| Standard HJE MBB | 14.0 | 1.53e-04 | 3.98 |
| Modified HJE Cantilever bottom | 6.31 | 6.70e-05 | 2.59 |
| Modified HJE Cantilever middle | 9.65 | 7.76e-05 | 3.06 |
| Modified HJE L-bracket top | 20.1 | 2.20e-04 | 9.49 |
| Modified HJE L-bracket mid | 19.3 | 2.19e-04 | 29.5 |
| Modified HJE MBB | 13.7 | 1.50e-04 | 4.74 |
3D-printed gripper compliance: LSSO \(1.04\times10^{-5}\), standard HJE \(8.49\times10^{-6}\), modified HJE \(8.33\times10^{-6}\), SIMP \(8.66\times10^{-6}\), PareTO \(8.29\times10^{-6}\) N·m. Waandishi wanaonya kuwa hizi si ranking ya jumla ya methods.
Runtime benchmark ilifanywa kwenye Intel Core Ultra 9 285K 3.70 GHz na RAM 128 GB. STORX ina overhead ya karibu 65–68% kwa resolutions za kati dhidi ya top88, karibu 60% kwenye 320×160, lakini iteration ya takriban elementi 50,000 ni karibu 0.4 s bila plotting.
Maelezo ya Chanzo na Mbinu
Kichwa asilia: STORX: An Open-Source Object-Oriented Framework for Shape and Topology Optimization in MATLAB.
Waandishi: Amir M. Mirzendehdel; Krishnan Suresh.
Aina ya chanzo: Preprint ya computational engineering na open-source software framework.
arXiv: 2606.17291v2 [cs.CE].
Tarehe ya PDF: 25 Julai 2026.
DOI: 10.48550/arXiv.2606.17291.
Leseni: CC BY 4.0.
Peer review: Utafiti huu ni preprint ambayo haijapitia peer review; matokeo yake yanapaswa kutathminiwa kwa kuzingatia hatua hiyo ya uchapishaji.
Code: https://github.com/DEL-KU/storx.
Conflict of interest: Waandishi wanasema hawana conflict of interest.
Generative AI: Ilitumika kuboresha clarity na readability ya baadhi ya sehemu za manuscript; technical content, results na interpretations zilitengenezwa na kuthibitishwa na waandishi.
Limitations: Core framework inalenga hasa 2D; benchmark results si universal ranking ya methods; object-oriented abstraction ina measurable runtime overhead; thermal na fluid cases ni demonstrations za extensibility na si industrial validation kamili.
Internal inconsistencies: Local-volume prose p=6 dhidi ya localPNorm=16, modified-HJE w factor na malformed active-area expression kwenye wind-tunnel zimehifadhiwa bila kuzisahihisha kimya kimya.

Acha maoni
Anwani yako ya barua pepe haitachapishwa. Sehemu za lazima zimewekewa alama ya *