FINISH RESUME,ed_gray-exper,db *IF,ARG1,EQ,1,THEN,*GO,:POST !--------------------------------------------------------------------------------------------------- epr = 3.000E+06 ! [-]. . . . . . . . . . . . . . . . cores relative permitivity mur = 3.000E+03 ! [-]. . . . . . . . . . . . . . . . cores relative permeability res = 1.715E+01 ! [Ohm.m]. . . . . . . . . . . . . . cores electric resistivity Rx = 8.880E+02 ! [Ohm]. . . . . . . . . . . . . . . winding resistivity Lx = 2.820E+00 ! [Henry]. . . . . . . . . . . . . . coil inductance Cx = 4.700E-05 ! [Farad]. . . . . . . . . . . . . . source capacity Nt = 3200 ! [-]. . . . . . . . . . . . . . . . number of turns Sc = 4.750E-04 ! [m2] . . . . . . . . . . . . . . . cross section area U = 4.000E+03 ! [Volt] . . . . . . . . . . . . . . initial voltage l = 3.000E-02 ! [second] . . . . . . . . . . . . . half of the pulses period q = 6.000E-03 ! [second] . . . . . . . . . . . . . pulses length Start = 52000 ! [Hz] . . . . . . . . . . . . . . . start frequency Stop = 98000 ! [Hz] . . . . . . . . . . . . . . . stop frequency !--------------------------------------------------------------------------------------------------- /PREP7 MP,MURX,2,mur MP,PERX,2,epr MP,RSVX,2,res FINISH *afu,rad PI = ACOS(-1) p1 = -Rx/2/Lx+SQRT(Rx**2/4/(Lx**2)-1/Lx/Cx) $p2 = -Rx/2/Lx-SQRT(Rx**2/4/(Lx**2)-1/Lx/Cx) AK= $BK= $FK= *DIM,AK,,700 *DIM,BK,,700 *DIM,FK,,700 *DO,I,1,700 k = I+Start*l*2 x1 = PI**2*k**2+l**2*p1**2 $x2 = PI**2*k**2+l**2*p2**2 t = q akq = U*(l*(exp(p1*t)*p1*x2-exp(p2*t)*p2*x1)*COS(k*PI*t/l)-k*PI*(exp(p2*t)*x1-exp(p1*t)*x2)*SIN(k*PI*t/l)) bkq = U*(l*(exp(p1*t)*p1*x2-exp(p2*t)*p2*x1)*SIN(k*PI*t/l)+k*PI*(exp(p2*t)*x1-exp(p1*t)*x2)*COS(k*PI*t/l)) t = 0 ak0 = U*(l*(exp(p1*t)*p1*x2-exp(p2*t)*p2*x1)*COS(k*PI*t/l)-k*PI*(exp(p2*t)*x1-exp(p1*t)*x2)*SIN(k*PI*t/l)) bk0 = U*(l*(exp(p1*t)*p1*x2-exp(p2*t)*p2*x1)*SIN(k*PI*t/l)+k*PI*(exp(p2*t)*x1-exp(p1*t)*x2)*COS(k*PI*t/l)) jm = Lx*(p1-p2)*x1*x2 ak(I) = (akq-ak0)/jm $bk(I) = (bkq-bk0)/jm fk(I) = k/l/2 *IF,fk(I),LT,Stop,THEN /SOLUTION ANTYPE,HARMIC HROUT,ON HARFRQ,fk(I),fk(I) $NSUBST,1 $KBC,1 BFV,8,JS,0,ak(I)*Nt/Sc,0,-90 $BFV,23,JS,0,ak(I)*Nt/Sc,0,-90 SOLVE BFV,8,JS,0,bk(I)*Nt/Sc,0,0 $BFV,23,JS,0,bk(I)*Nt/Sc,0,0 SOLVE /COPY,file,rmg,,ed_gray-exper%I%,rmg FINISH *ENDIF *ENDDO /DELETE,ed_gray-response,rmg !------------------------------------------------------------------ :POST /eof /POST1 *DO,I,1,10 tx = I/100 *DO,J,1,5 FILE,ed_gray-exper%J%,rmg SET,1,,1,0 *GET,fx,ACTIVE,0,SET,FREQ LCDEF,1,1,1,0 LCDEF,2,1,1,1 LCDEF,3,1,1,0 LCDEF,4,1,1,1 LCFACT,1,SIN(2*PI*fx*tx) LCFACT,2,COS(2*PI*fx*tx) LCFACT,3,SIN(2*PI*fx*tx) LCFACT,4,COS(2*PI*fx*tx) LCASE,1 LCOPER,ADD,2 LCASE,3 LCOPER,ADD,4 LCASE,1 LCOPER,ADD,3 *ENDDO FILE,ed_gray-response,rmg RAPPND,1,tx *ENDDO