SCUBA DIVER OPTIMIZATION ALGORITHM: Pressure Vessel Design

SCUBA DIVER OPTIMIZATION ALGORITHM: Pressure Vessel Design A NOVEL METAHEURISTIC ALGORITHM FOR CONSTRAINED ENGINEERING DESIGN

% SCUBA DIVER OPTIMIZATION ALGORITHM: Pressure Vessel Design

%

% A NOVEL METAHEURISTIC ALGORITHM FOR CONSTRAINED ENGINEERING DESIGN

% ----------------------------------------------------------------------- %

% Author: Assistant Professor Saman M. Almufti

% Novel Metaheuristic Algorithm (PhD Thesis Contribution)

% ----------------------------------------------------------------------- %

% Application: Pressure Vessel Design Optimization

% This implementation adapts SDOA to minimize the total material, forming,

% and welding cost of a cylindrical pressure vessel while satisfying the

% shell-thickness, head-thickness, volume, and length constraints.

%

% Design variables:

% x(1) = Ts : shell thickness (discrete, multiple of 0.0625 in)

% x(2) = Th : head thickness (discrete, multiple of 0.0625 in)

% x(3) = R : inner radius (continuous, in)

% x(4) = L : cylinder length (continuous, in)

% ----------------------------------------------------------------------- %

% E-mail: [e-mail]

% Google Scholar: https://scholar.google.com/citations?user=bYCPWOgAAAAJ&hl=en

% ORCID: https://orcid.org/0000-0002-1843-745X

% LinkedIn: https://www.linkedin.com/in/saman-m-almufti/

% ----------------------------------------------------------------------- %

% Standard benchmark formulation:

%

% Minimize

% f(x) = 0.6224*Ts*R*L + 1.7781*Th*R^2

% + 3.1661*Ts^2*L + 19.84*Ts^2*R

%

% Subject to

% g1(x) = -Ts + 0.0193*R <= 0

% g2(x) = -Th + 0.00954*R <= 0

% g3(x) = -pi*R^2*L - (4/3)*pi*R^3 + 1296000 <= 0

% g4(x) = L - 240 <= 0

%

% Bounds

% 0.0625 <= Ts, Th <= 6.1875

% 10 <= R, L <= 200

%

% Reference global solution for the mixed-discrete benchmark:

% [Ts, Th, R, L] = [0.8125, 0.4375, 42.0984455958549,

% 176.6365958424394]

% f* = 6059.714335048436

% ----------------------------------------------------------------------- %

% File: PV_SDOA.m

% Run: >> PV_SDOA

% ======================================================================= %