shell.th 3[mm] "Thickness, shell"
shell.Do 75[mm] "Outer diameter, shell"
shell.Ro shell.Do/2 "Outer radius, shell"
shell.Ri shell.Ro-shell.th "Inner radius, shell"
shell.hsf 3.5*shell.th "Straight flange height, shell"
shell.Rk 0.2*shell.Do "Internal knuckle radius, shell"
shell.Rc 1.4*shell.Ri "Internal crown radius, shell"
shell.hi shell.Rc-sqrt((shell.Rc-shell.Ri)*(shell.Rc+shell.Ri-2*shell.Rk)) "Internal head height, shell"
shell.alpha_k atan((shell.Ri-shell.Rk)/(shell.Rc-shell.hi)) "Angle at intersection crown-knuckle, shell"
shell.alpha_bp 20[deg] "Bottom plate intersection angle"
dia.th shell.th/2 "Thickness, diaphragm"
dia.thb dia.th*1.8 "Bead thickness, diaphragm"
dia.Zb dia.th/2 "Bead position, diaphragm"
dia.Ro shell.Ri-shell.th*2 "Outer radius, diaphragm"
dia.Ri dia.Ro-dia.th "Inner radius, diaphragm"
dia.alpha 12[deg] "Angle parameter, diaphragm"