d1 = {μFe -> 
            500*4*Pi*(10^(-7))*(1 - 0.005*T[r, t])*(1 - 
                  0.09*B[r, t]), γFe -> 
            (10^6)*(1 - 0.02*T[r, t]), rmin -> 0.001, 
          rmax -> 2*(10^(-2)), ki -> 1, 
          u[t] -> 200  Sin[314  t], 
          tmax ->  .1};
d2 = {λFe -> 5 - 0.01*T[r, t], 
          T2 -> 20, α -> 5, ϵσ -> 
            4*(10^(-8)), 
          Tv -> 20, ρcFe -> 1*(10^5)*(1 - 0.03*T[r, t]), 
          T0 -> 20};
d3 = {γCu -> 
            57*(10^6)*(1 - 0.02*Tind[r, t]), ϵσext -> 
            4*(10^(-8)), 
          kiext -> 0.1, ρcCu -> 
            5*(10^5) (1 - 0.03*Tind[r, t]), kR -> 4*(10^12), 
          kS -> 0.1, αext -> 5, T0ext -> 20, 
          Text -> 20, λCu -> 
            10000};
d4 = {R -> 1, L0 -> 1, 
          k -> 1};

d1234 = 
        Union[d1, d2, d3, 
          d4];

rceI = 
        R*i[r, t] + L0*Derivative[0, 1][i][r, t] + 
            (10^20)*2*Pi*
              (rmax/γ)Fe*Derivative[1, 0][B[r, t]/μFe][
                rmax, t] == 
          u[t] + 10000*D[i[r, t], {r, 2}];
okr1I = Derivative[1, 0][i][rmin, t] == 
          0;
okr2I = Derivative[1, 0][i][rmax, t] == 
          0;
pocI = i[r, 0] == 0;
argI = {rceI, okr1I, okr2I, pocI} /. 
          d1234;

rceB = 
        Simplify[
          D[B[r, t], t] - 
              (1/r)*D[(r/γ)Fe*D[B[r, t]/μFe, r], r] == 
            0];
okr1B = Derivative[1, 0][B][rmin, t] == 
          0;
okr2B = 
        B[rmax, t] - ki*i[rmax, t] == 0;
pocB = B[r, 0] == 0;
oborer = {r, rmin, rmax};
oborte = {t, 0, tmax};
argB = {rceB, okr1B, okr2B, pocB} /. 
          d1234;


rceT = (1/r) D[r*λFe*D[T[r, t], r], r] - ρcFe*
              D[T[r, t], t] + 
            (1/(γFe*μ(Fe^2)))*D[B[r, t], r]^2 == 
          0;
okr1T = Derivative[1, 0][T][rmin, t] == 
          0;
okr2T = Derivative[1, 0][T][rmax, 
            t] == α*(T[rmax, t] - 
                  Tv) + ϵσ*((((T[rmax, t] + 273))^4) - (((T2 + 273))^4));
pocT = T[r, 0] == T0;
argT = {rceT, okr1T, okr2T, pocT} /. 
          d1234;


rceTind = ρcCu*D[Tind[r, t], t] == 
          i[r, t]^2*(kR/γ)Cu - 
            kS*αext*(Tind[r, t] - Text) - 
            kS*ϵσext*((((Tind[r, t]))^4) - (((Text + 273))^4)) + 
            kS*ϵσ*((((T[r, t] + 273))^4) - (Tind[r, t] + 
                        273)^4) + 
            kS*α*(T[r, t] - Tind[r, t]) + λCu*
              D[Tind[r, t], {r, 2}];
okr1Tind = Derivative[1, 0][Tind][rmin, t] == 
          0;
okr2Tind = Derivative[1, 0][Tind][rmax, t] == 
          0;
pocTind = Tind[r, 0] == T0;
argTind = {rceTind, okr1Tind, okr2Tind, pocTind} /. 
          d1234;

arg = 
        Union[argI, argB, argT, 
          argTind];

resBTTind = NDSolve[arg, {i[r, t], T[r, t], B[r, t], Tind[r, t]}, 
            Evaluate[oborer //. d1, oborte //. d1], 
            MaxSteps -> 
              (10^5)][[1]];

proud[r_, t_] = 
        Evaluate[i[r, t] /. resBTTind];
Plot3D[proud[r, t], Evaluate[oborer //. d1, oborte //. d1], 
      PlotPoints -> 50, 
      PlotRange -> All]

indukce[r_, t_] = 
        Evaluate[B[r, t] /. resBTTind];
Plot3D[indukce[r, t], Evaluate[oborer //. d1, oborte //. d1], 
      PlotPoints -> 50, 
      PlotRange -> All]

teplota[r_, t_] = 
        Evaluate[T[r, t] /. resBTTind];
Plot3D[teplota[r, t], Evaluate[oborer //. d1, oborte //. d1], 
      PlotPoints -> 50, 
      PlotRange -> All]
teplotaind[r_, t_] = 
        Evaluate[Tind[r, t] /. resBTTind];
Plot3D[teplotaind[r, t], Evaluate[oborer //. d1, oborte //. d1], 
      PlotPoints -> 50, 
      PlotRange -> 
        All]