λ = 0.82;
ρ = 1810;
cp = 1000;
To = -10;
d = 0.4;
αext = 20;
αint = 5;
Texter = 
        To + 10*Sin[((2  Pi)/(3600*24))*t];
Texint = -10 + 30*(t/(t + 1));
rce = ρ*cp*D[T[x, t], t] == λ*
            D[T[x, t], {x, 2}];
poc = T[x, 0] == To;
okrleft = -λ*
              D[T[x, t], x] == αext*(Texter - T[x, t]) /. 
          x -> 0;
okrright = -λ*
              D[T[x, t], x] == αint*(-Texint + T[x, t]) /. 
          x -> d;
tmax = 3600*24*3;
res = NDSolve[{rce, poc, okrleft, okrright}, 
            T[x, t], {x, 0, d}, {t, 0, 
              tmax}][[1]];
teplota = T[x, t] /. res;
Plot3D[teplota, {x, 0, d}, {t, 0, tmax}, 
      AxesLabel -> {"x[m]", "t[s]", "T[°C]"}, 
      PlotRange -> All]

q = αint*(-Texint + T[x, t]) /. res /. 
          x -> d;
Tinzed = T[x, t] /. res /. 
          x -> d;
\(Plot[q, {t, 0, tmax}, 
      AxesLabel -> {"t[s]", \*""q[W.\!\(m\^\(-2\)\)]""}, 
      PlotRange -> All]\)
Plot[Tinzed, {t, 0, tmax}, 
      AxesLabel -> {"t[s]", "T[°C]"}, 
      PlotRange -> All]