Компьютерные технологии в физике. Часть 1. Компьютерное моделирование физических процессов. Красов В.И - 71 стр.

UptoLike

ВУЗ: 

ИГУ | Иркутск

Составители: 

Рубрика: 

74
поля
z
B обратится в 0 в силу аксиальной симметрии. Для оставшихся
компонент
x
B
и
y
B получим:
∫∫
ππ
αα
−
α
=
2
0
3
2
0
3
2
cos
r
d
yRI
r
d
RIB
x
(4.16)
∫
π
αα
=
2
0
3
cos
r
d
xRIB
x
(4.17)
Уравнения (4.15), (4.16), (4.17) определяют магнитную индукцию
кольца с током в произвольной точке плоскости
(
)
YX , в безразмерном ви-
де.
Для вычисления интегралов можно использовать метод трапеций.
Соответствующая подпрограмма может выглядеть так:
{В данной процедуре используются глобальные константы: x0,y0 - коор-
динаты центра кольца; R0 - радиус кольца}
Procedure Bfield(x,y: real);
const
N = 100; {число шагов интегрирования}
var
j : integer;
r,alpha,delta,s1,s2 : real:
begin
alpha := 0;
s1:=0;
s2:=0;
delta:=2*Pi/N;
for j:=1 to N do
begin
r:=sqrt(sqr(x-x0)+sqr(y-y0)+sqr(R0)-2*(y-y0)*R0*cos(alpha);
s1:=s1+delta/(2*r*r*r);
s2:=s2+cos(alpha)*delta/(2*r*r*r);
(a3)
alpha:=alpha+delta;
r:=sqrt(sqr(x-x0)+sqr(y-y0)+sqr(R0)-2*(y-y0)*R0*cos(alpha);
s1:=s1+delta/(2*r*r*r);
s2:=s2+cos(alpha)*delta/(2*r*r*r);
end;
Bx:=I*sqr(R0)*s1-I*(y-y0)*R0*s2;
By:=I*(x-x0)*R0*s2;
PDF created with FinePrint pdfFactory Pro trial version http://www.fineprint.com
              поля B z обратится в 0 в силу аксиальной симметрии. Для оставшихся
              компонент B x и B y получим:

                                        2π              π
                                       dα       2
                                                   cos α dα
                           B x = I R ∫ 3 − I yR ∫
                                    2
                                                                                   (4.16)
                                     0 r         0    r3

                                        2π
                                             cos α dα
                           B x = I xR   ∫      r3
                                                                                   (4.17)
                                        0


                   Уравнения (4.15), (4.16), (4.17) определяют магнитную индукцию
              кольца с током в произвольной точке плоскости ( X , Y ) в безразмерном ви-
              де.
                   Для вычисления интегралов можно использовать метод трапеций.
              Соответствующая подпрограмма может выглядеть так:

               {В данной процедуре используются глобальные константы: x0,y0 - коор-
                        динаты центра кольца; R0 - радиус кольца}
               Procedure Bfield(x,y: real);
               const
                   N = 100;       {число шагов интегрирования}
               var
                   j                   : integer;
                   r,alpha,delta,s1,s2 : real:
               begin
                 alpha := 0;
                 s1:=0;
                 s2:=0;
                 delta:=2*Pi/N;
                 for j:=1 to N do
                 begin
                    r:=sqrt(sqr(x-x0)+sqr(y-y0)+sqr(R0)-2*(y-y0)*R0*cos(alpha);
                    s1:=s1+delta/(2*r*r*r);
                    s2:=s2+cos(alpha)*delta/(2*r*r*r);                      (a3)
                    alpha:=alpha+delta;
                    r:=sqrt(sqr(x-x0)+sqr(y-y0)+sqr(R0)-2*(y-y0)*R0*cos(alpha);
                    s1:=s1+delta/(2*r*r*r);
                    s2:=s2+cos(alpha)*delta/(2*r*r*r);
                end;
                Bx:=I*sqr(R0)*s1-I*(y-y0)*R0*s2;
                By:=I*(x-x0)*R0*s2;

                                                            74


PDF created with FinePrint pdfFactory Pro trial version http://www.fineprint.com

Страницы