C1=40.0 # Overall alloy composition
Cb=98.0
Ca=10.0
Wa=(Cb-C1)/(Cb-Ca)
Wb=(C1-Ca)/(Cb-Ca)
d_Sn=7.24 # in g/cm**3 density of tin
d_Pb=11.23 # in g/cm**3 density of lead
Ca_Sn=10.0
Ca_Pb=90.0
Cb_Sn=98.0
Cb_Pb=2.0
d_a=100/((Ca_Sn/d_Sn)+(Ca_Pb/d_Pb))
d_b=100/((Cb_Sn/d_Sn)+(Cb_Pb/d_Pb))
Va=Wa/(d_a*((Wa/d_a)+(Wb/d_b)))
Vb=Wb/(d_b*((Wa/d_a)+(Wb/d_b)))
print"Mass fractions for alpha and beta phases are respectively",round(Wa,2),"and",round(Wb,2)
print"Density of alpha phase is ",round(d_a,2),"g/cm**3"
print"Density of beta phase is ",round(d_b,2),"g/cm**3"
print"Volume fraction of alpha phase ",round(Va,2),"g/cm**3"
print"Volume fraction of beta phase ",round(Vb,2)
C0=0.35
Ca=0.022
C_Fe3C=6.7
Wa=(C_Fe3C-C0)/(C_Fe3C-Ca)
W_Fe3C=(C0-Ca)/(C_Fe3C-Ca)
C_p=0.76
Wp=(C0-Ca)/(C_p-Ca)
W_a=(C_p-C0)/(C_p-Ca)
Wp=(C0-Ca)/(C_p-Ca)
W_a=(C_p-C0)/(C_p-Ca)
Wa=(C_Fe3C-C0)/(C_Fe3C-Ca)
Wae=Wa-W_a
print"Mass fraction of total ferritic phase is",round(Wa,2)
print"Mass fraction of Fe3C is",round(W_Fe3C,2)
print"Mass fraction of Pearlite is",round(Wp,2)
print"ass fraction of proeutectoid ferrite is",round(W_a,2)
print"Mass fraction of eutectoid ferrite is",round(Wae,2)