Pressure Vessel Thickness Calculator

def cylindrical_thickness(P, D, S, E, CA=0, min_t=0):
    """
    Calculate minimum wall thickness for a cylindrical pressure vessel.
    
    P: Internal pressure (MPa)
    D: Inside diameter (mm)
    S: Allowable stress (MPa)
    E: Joint efficiency (0-1)
    CA: Corrosion allowance (mm)
    min_t: Minimum required thickness (mm)
    """
    if P <= 0 or D <= 0 or S <= 0 or E <= 0 or E > 1:
        raise ValueError("Invalid input parameters")
    if CA < 0 or min_t < 0:
        raise ValueError("Corrosion allowance and min thickness cannot be negative")
    
    numerator = P * D
    denominator = 2 * S * E - 1.2 * P
    if denominator <= 0:
        raise ValueError("Denominator non-positive. Check inputs.")
    t = numerator / denominator + CA
    return max(t, min_t)

# Example
P = 1.5
D = 1000
S = 138
E = 0.85
CA = 3
min_t = 0

t_cyl = cylindrical_thickness(P, D, S, E, CA, min_t)
print(f"Required thickness: {t_cyl:.2f} mm")

function t = cylindrical_thickness(P, D, S, E, CA, min_t)
    % Calculate minimum wall thickness for a cylindrical pressure vessel.
    if nargin < 5, CA = 0; end
    if nargin < 6, min_t = 0; end
    if P <= 0 || D <= 0 || S <= 0 || E <= 0 || E > 1
        error('Invalid input parameters');
    end
    if CA < 0 || min_t < 0
        error('Corrosion allowance and min thickness cannot be negative');
    end
    
    numerator = P * D;
    denominator = 2 * S * E - 1.2 * P;
    if denominator <= 0
        error('Denominator non-positive. Check inputs.');
    end
    t = numerator / denominator + CA;
    t = max(t, min_t);
end

% Example
P = 1.5;
D = 1000;
S = 138;
E = 0.85;
CA = 3;
min_t = 0;

t_cyl = cylindrical_thickness(P, D, S, E, CA, min_t);
fprintf('Required thickness: %.2f mm\n', t_cyl);