//============================================================================
// nom : coef_dilatation_acier.sce
// auteur : Christophe Dang Ngoc Chan
// date de création : 2013-07-10
// dates de modification : ...
//----------------------------------------------------------------------------
// version de Scilab : 5.4.1
// module Atoms requis : aucun
//----------------------------------------------------------------------------
// Objectif : trace des courbes alpha = f(T)
// Entrées : aucune (codé en dur)
// Sorties : fenêtre graphique
//============================================================================

// valeurs en 10^-6 K^-1 (source : OTUA)

carbone = [ 100  200  300  400  500  600 
        11.1 12.1 12.9 13.5 13.9 14.1]; // C35E 1.1181
inox_ferr = [   100  200  300  400  500
                10.4 10.8 11.2 11.6 11.9]; // X2CrNi12 1.4003
inox_aust = [   100 200 300 400 500
                16  17  17  18  18]; // X2CrNiMo17-12-2 1.4404
inox_aust_ferr = [  100 200  300
                    13  13.5 14]; // X2CrNiMoN22-5-3 1.4462
inox_martens = [    100  200 300  400
                    10.5 11  11.5 12]; // X20Cr13 1.4021

// tracé

cmap = coolcolormap(5);
xset("colormap", cmap)
plot2d(inox_ferr(1,:), inox_ferr(2,:), style=1);
plot2d(inox_martens(1,:), inox_martens(2,:), style = 2);
plot2d(carbone(1,:), carbone(2,:), style = 3);
plot2d(inox_aust_ferr(1,:), inox_aust_ferr(2,:), style = 4);
plot2d(inox_aust(1,:), inox_aust(2,:), style = 5);
xtitle("Coefficients de dilatation linéiques",...
    "$\mathrm{T}\ (°\mathrm{C})$", "$\alpha\ (10^{-6}\mathrm{K}^{-1})$")

legend("X2CrNi12 (1.4003, 403)", ...
    "X20Cr13 (1.4021, 420)",...
    "C35E (1.1181, 1035)",... 
    "X2CrNiMoN22-5-3 (1.4462, 2205)",...
    "X2CrNiMo17-12-2 (1.4404, 316L)");