DIPLOMOVÁ PRÁCA - Žilinská univerzita

ccc1=(((F/(Ic_B/Ic))^(2*Ic_B/Ic))*2*Ic_B/Ic)*0.5*(1+sign(Ic_B/Ic-
F));
ccc2=(((F/(Ic_B/Ic))^(-2*Ic_B/Ic))*2*Ic_B/Ic)*0.5*(1+sign(F-
Ic_B/Ic));
Qmag=0.4*(1+20*alfa^2.1)*(ccc1+ccc2);
DPmag(k,m)=(Qmag*4e-7*pi*(I1_Ia_max^2)*f)/(2*pi);
DP(k,m)=DPmag(k,m)+DPr(k,m);
DP_SC=DP_SC+DP(k,m)*dlzka_zav(k);
end;
zavit=k
end;
for k=1:(N1s/2)
for m=1:21
poloha(1,1)=vrstvy1(k);
poloha(2,1)=LV_y(m);
Br_LV(k,m)=postinterp(fem0,'Bx_qa',poloha);
F=I_LV_max/Ic;
alfa=(abs(Br_LV(k,m))/Bp)/F;
n=12*exp(-F*alfa*5e-4);
Ic_B=Ic*(1.1*((1+((alfa+0.1)/2)^2)^-0.5)-0.1*exp(-alfa/0.7));
DPr(k,m)=((1e-4/(Ic_B^n))*I_LV_max^n)*I_LV_max;
ccc1=(((F/(Ic_B/Ic))^(2*Ic_B/Ic))*2*Ic_B/Ic)*0.5*(1+sign(Ic_B/Ic-
F));
ccc2=(((F/(Ic_B/Ic))^(-2*Ic_B/Ic))*2*Ic_B/Ic)*0.5*(1+sign(F-
Ic_B/Ic));
Qmag=0.4*(1+20*alfa^2.1)*(ccc1+ccc2);
DPmag(k,m)=(Qmag*4e-7*pi*(I_LV_max^2)*f)/(2*pi);
DP(k,m)=DPmag(k,m)+DPr(k,m);
DP_SC=DP_SC+DP(k,m)*dlzka_zav(k);
end;
zavit=k
end;
for k=1:(Nks/2)
for m=1:8
poloha(1,1)=vrstvy1(k);
poloha(2,1)=Kur_y(m);
Br_Kur(k,m)=postinterp(fem0,'Bx_qa',poloha);
F=I_kur_max/Ick;
alfa=(abs(Br_Kur(k,m))/Bp)/F;
n=12*exp(-F*alfa*5e-4);
Ic_B=Ick*(1.1*((1+((alfa+0.1)/2)^2)^-0.5)-0.1*exp(-alfa/0.7));
DPr(k,m)=((1e-4/(Ic_B^n))*I_kur_max^n)*I_kur_max;
ccc1=(((F/(Ic_B/Ick))^(2*Ic_B/Ick))*2*Ic_B/Ick)*0.5*(1+sign(Ic_B/Ick-
F));
ccc2=(((F/(Ic_B/Ick))^(-2*Ic_B/Ick))*2*Ic_B/Ick)*0.5*(1+sign(F-
Ic_B/Ick));
Qmag=0.4*(1+20*alfa^2.1)*(ccc1+ccc2);
DPmag(k,m)=(Qmag*4e-7*pi*(I_kur_max^2)*f)/(2*pi);
DP(k,m)=DPmag(k,m)+DPr(k,m);
DP_SC=DP_SC+DP(k,m)*dlzka_zav(k);
end;
zavit=k
end;