Программный комплекс расчета и моделирования катушки управления электромагнита постоянного тока (стр. 2 из 4)
, (13)где мощность тепловых потерь ОУ
.Сопротивление ОУ в горячем состоянии
(14)где ρ0 - удельное сопротивление материала обмоточного провода ОУ при температуре окружающей среды Q0, a - температурный коэффициент сопротивления материала обмоточного провода при температуре окружающей среды Q0. Коэффициент теплоотдачи kт в выражении (13) также является функцией от температуры. С учетом этого на основании (13) получается нелинейное уравнение вида
, (15)Алгоритмы и программа расчета обмотки управления
 | ||||
 | ||||
 |  | |
 | ||
Программа расчета обмотки управления.
program kr;
uses graph , crt;
const m0=4*3.1415926535897932385*0.00000010;pi=3.1415926535897932385;n1=33;
label 2,3,4,5,6,7,8,9,10,12,14,15,16,25,26,20,29;
var F,Unom,Umax,Umin,dc,e1,e2,h,Sn,km,y1,y2,dh,dx,pgor,Tdop,T0,T,T1,Ta,Tb,
kz,Q,Soxl,lcp,km1,l,n,d0,m,avn,p0,a,hnp,Svn,ha,hb,z,dp,R,W,dci,e,P,Rxc,Rgc,
t1dop,Pp,Ik:real;
D1,D,K:array[1..n1] of real; i,j:integer;
begin
clrscr;
{ishodnie dannie}
F:=4800;Unom:=110;dc:=24*0.001;Tdop:=105;T0:=40;
n:=1;Umax:=1.15*Unom;Umin:=0.9*Unom;kz:=0.5;m:=5;
avn:=1.7; p0:=1.62*0.00000001;a:=4.3*0.001; d0:=dc+0.0001;
pgor:=p0*(1+a*Tdop);
t1dop:=Tdop-T0;
{Raschet minimalnih razmerov OU metodom polovinnogo delenia}
{write('vvedite ha:= ');readln(ha);
write('vvedite hb:= ');readln(hb);}
ha:=0.001;hb:=0.100;
e1:=0.000001;e2:=0.000001; T:=Tdop;
5:
h:=ha;
l:=m*h;
Sn:=pi*(d0+2*h)*l;
Svn:=pi*d0*l;
Q:=h*l;
Soxl:= Sn+avn*Svn;
lcp:=pi*(d0+h);
if (Soxl>=0.0001) and (Soxl<=0.01) then km:=2.1*(1+0.005*(T-T0))/exp(1/3*ln(Soxl));
if (Soxl>0.01) and (Soxl<=0.5) then km:=3.6*(1+0.005*(T-T0))/exp(1/5*ln(Soxl));
if Soxl>0.5 then goto 2;
if Soxl<0.0001 then goto 3;
y1:=sqr(Umax/Umin)*sqr(F)*pgor/n/km/t1dop/kz-Q*Soxl/lcp;
h:=(ha+hb)/2;
l:=m*h;
Sn:=pi*(d0+2*h)*l;
Svn:=pi*d0*l;
Q:=h*l;
Soxl:= Sn+avn*Svn;
lcp:=pi*(d0+h);
if (Soxl>=0.0001) and (Soxl<=0.01) then km:=2.1*(1+0.005*(T-T0))/exp(1/3*ln(Soxl));
if (Soxl>0.01)and(Soxl<=0.5) then km:=3.6*(1+0.005*(T-T0))/exp(1/5*ln(Soxl));
if Soxl>0.5 then goto 2;
if Soxl<0.0001 then goto 3;
y2:=sqr(Umax/Umin)*sqr(F)*pgor/n/km/t1dop/kz-Q*Soxl/lcp;
if abs(y2)<e1 then goto 4 else begin
if y1*y2<0 then hb:=h else ha:=h;
if abs(hb-ha)<e2 then goto 4 else goto 5 end;
{Raschet i vibor diametra provoda iz standartnogo rada,raschet parametrov
OU v holodnom sostoanii}
4:
D[1]:=0.050*0.001; K[1]:=0.250; D1[1]:=0.080*0.001;
D[2]:=0.063*0.001; K[2]:=0.290; D1[2]:=0.090*0.001;
D[3]:=0.071*0.001; K[3]:=0.340; D1[3]:=0.100*0.001;
D[4]:=0.080*0.001; K[4]:=0.380; D1[4]:=0.110*0.001;
D[5]:=0.090*0.001; K[5]:=0.410; D1[5]:=0.120*0.001;
D[6]:=0.100*0.001; K[6]:=0.440; D1[6]:=0.130*0.001;
D[7]:=0.112*0.001; K[7]:=0.460; D1[7]:=0.140*0.001;
D[8]:=0.125*0.001; K[8]:=0.480; D1[8]:=0.155*0.001;
D[9]:=0.140*0.001; K[9]:=0.495; D1[9]:=0.170*0.001;
D[10]:=0.160*0.001;K[10]:=0.510;D1[10]:=0.200*0.001;
D[11]:=0.180*0.001;K[11]:=0.527;D1[11]:=0.220*0.001;
D[12]:=0.200*0.001;K[12]:=0.538;D1[12]:=0.240*0.001;
D[13]:=0.224*0.001;K[13]:=0.550;D1[13]:=0.270*0.001;
D[14]:=0.250*0.001;K[14]:=0.560;D1[14]:=0.300*0.001;
D[15]:=0.280*0.001;K[15]:=0.572;D1[15]:=0.330*0.001;
D[16]:=0.315*0.001;K[16]:=0.580;D1[16]:=0.365*0.001;