Пособие по практике "Теория электромагнитного поля"

ПРИЛОЖЕНИЕ П12 Программа расчета разветвленной магнитной системы с интерполяцией кривой намагничивания PROGRAMMAGNECIR; USES GRAPH; VAR ВВ,В 11 ,В 12,В 13 ,В21 ,В22,В23 ,ВЗ 1 ,BDE, ВМАХ,С1 ,C2,DB,DEB,DEL,DF1 ,DF2,DF3, FF1 ,FF2 ,FF3 ,FL1 ,FL2 ,FL3 ,G,H,HH,L 11 ,L 12, L13 ,L21 ,L22,L23,L31 ,LM,UO,MUO,MX,MY, Q,UM1,UM2,UM3,S,S11,S12,S13,S21,S22,S23, S31 ,SDE,X,X0,X1 ,X2,X3,X12,Y,Y0,Y1 ,Y2,YM: REAL; I,J,N,NN,GD,GM,X0,XG,XS,Y0,YG,YS: INTEGER; LOG: BOOLEAN; MA: ARRAY[-L.20] OF REAL; MAX,MAY: ARRAY[-20..20,L.3] OF INTEGER; MACO: ARRAY[1..3] OF INTEGER; FIMA: TEXT; FUNCTION FHH(BB: REAL): REAL; BEGIN LOG:=FALSE; IF BB<0 THEN BEGIN BB:=-BB; LOG:=TRUE END; IF BB<=1.8 THEN BEGIN S:=BB/DB; J:=TRUNC(S); H:=S-J; YM:=MA[J-1]; YO:=MA[J]; Y1:=MA[J+1]; Y2:=MA[J+2]; FHH:=YO+H*0.5*(Y1-YM)+ H*H*(YM-2.5*YO+2*Y1-0.5*Y2)+ H*H*H*0.5*(Y2-3*Y1+3*YO-YM) END ELSE FHH:=-49450+32000*BB; IF LOG THEN FHH:=-FHH(BB) END; {$I D:\BP\BIN\KOORSET.PAS'} BEGIN ASSIGN(FIMA;D:\BP\BIN\FIMA.TXT'); REWRITE(FIMA); MA[-l]:=-80; MA[0]:=0; MA[1]:=80; MA[2]:=160; MA[3]:=240; MA[4]:=320; MA[5]:=400; MA[6]:=488; MA[7]:=584; MA[8]:=682; MA[9]:=798; MA[10]:=920; MA[11]:=1090; MA[12]:=1298; MA[13]:=1600; MA[14]:=2100; MA[15]:=2900; MA[16]:=4100; MA[17]:=5800; MA[18]:=8150; MA[19]:=11350; MAC0[1]:=12; MAC0[2]:=11; MACO[3]:=14; DB:=0.1; N:=200; G:=(20*DB)/N; GD:=DETECT; 35