%% Thick Wall Assumption Hoop Stress
syms ro
pi = 2000;
po = 14.7;
ri = 4/2;
yield = 35000;
sf = 2;
hoop_str_thick = sf*((pi*ri^2-po*ro^2)-(ro^2*(po-pi)))/(ro^2-ri^2) == yield;
sol = double(solve(hoop_str_thick, ro))

%% Thick Wall Assumption Axial Stress
syms ro
pi = 2000;
po = 14.7;
ri = 4/2;
yield = 35000;
sf = 2;
hoop_str_thick = 2*(pi*ri^2-po*ro^2)/(ro^2-ri^2) == yield;
sol = double(solve(hoop_str_thick, ro))
%% Thin Wall Assumption Hoop Stress
syms t
p = 2000;
yield = 35000;
r = 2; %OD
sf = 2;
hoop_str_thin = sf*p == yield*t/r
sol = double(solve(hoop_str_thin, t))

%% Thin Wall Assumption Axial Stress
syms t
p = 2000;
yield = 35000;
r = 2; %OD
sf = 2;
hoop_str_thin = sf*p == 2*yield*t/r
sol = double(solve(hoop_str_thin, t))