flansch.het revision 624e6701e0deb7ac6c03c0cba0190fbc5033cf93
library flansch
logic Reduce
spec Flansch1 =
. d_G1:=500.0
. d_G2:=550.0
. d_4:=730.0
. d_0:=488.0
. d_3:=660.0
. d_2:=534.0
. d_1:=498.0
. d_Be:=29.72
. d_5:=36
. A_F:=5094
. e_G:=2.0
. E_F0:=199964.0
. e_1:=10.0
. e_2:=35.0
. e_P:=42.10
. dBs:=31.0
. n_B:=20
. f_B0:=20
. F_A0:=0
. M_A0:=0
. l_H:=61.0
. Q_maxy:=160.0
. phi_S:=0
then
. d_3e:=d_3 * (1-2/n_B^2)
. A_B:=min(d_Be,dBs)^2 * n_B * Pi/4
. F_R0:=F_A0+(4/d_3e)*M_A0
. F_G[0]:=A_B+f_B[0]/3 -F_R[0]
. b_Gt:=(d_G2 - d_G1)/2
. d_F:=(d_4 + d_0)/2
. beta:=e_2/e_1
. p_B:=Pi * d_3/n_B
. d_5e:=d_5 * sqrt(d_5/p_B)
. b_F:=(d_4 - d_0)/2 - d_5e
. e_E:=e_1 * (1 + (beta-1)*l_H/((beta/3) * sqrt(d_1*e_1) + l_H))
. d_E:=(min(d_1-e_1+e_E,d_2+e_2-e_E) + max(d_1+e_1-e_E,d_2-e_2+e_E))/2
. gamma:=e_E * d_F/(b_F * d_E * cos(phi_S))
%% . e_F:=2 * A_F/(d_4 - d_0)
. A_F:=e_F * (d_4 - d_0) / 2
. llambda:=1 - e_P/e_F
. theta:=0.55 * cos(phi_S) * sqrt(d_E * e_E)/e_F
. E_Gm:=530.0
. c_F:=(1 + gamma * theta) / (1 + gamma * theta * (4 * (1 - 3 * llambda + 3 * llambda^2) + 6 * (1 - 2 * llambda) * theta + 6 * theta^2) + 3 * gamma^2 * theta^4)
. Z_F:=3*d_F * c_f/(Pi * b_F * e_F^3)
. b_Gi:=b_Gt
. repeat
b_Ge:=min(b_Gi,b_Gt);
d_Ge:=d_G2 - b_Ge;
A_Ge:=Pi * d_Ge * b_Ge;
h_G0:=(d_3e-d_Ge)/2;
b_Gi:=sqrt(e_G/(Pi * d_Ge * E_Gm)/(h_G0 * Z_F/E_F0 + h_G0 * Z_F/E_F0) + (F_G0 / (Pi * d_Ge * Q_maxy))^2)
until convergence(0.001,b_Ge) %% && F_G0req <= F_G0
end