Расчет электростатического поля несоосных параллельных цилиндров (Еr = 1; l = 7 м; U = 2000, D = 0,42 м), страница 2

Напряжённость в определённой точке пространства в нашем случае определяется, как суперпозиция двух этих напряжённостей в данной точке:
Image(1)
α+ и α- - углы между радиус-векторами и вектором оси для положительно и отрицательно заряженной оси, соответственно.

Если выбираемые точки лежат на одной прямой, перпендикулярной отрезку, соединяющему нити, и делящей его пополам, то: 

Image
и конечное выражение для вектора напряжённости можно записать в форме:
Image

а так как направление вектора в данной точке детерминировано, то можно перейти к скаляру
Image
Image- y - координата точки, отсчитываемая от одной из нитей (*).
E(y) - функция напряжённости электрического поля в точках данной поверхности.
Выражения потока примет следующий вид:
Image
y для каждой силовой линии в 1 и 2-й четверти чертежа будет координатой пересечения соответствующей окружности и оси Y в положительной области.

y для каждой силовой линии в 3 и 4-й четверти чертежа будет координатой пересечения соответствующей окружности и оси Y в отрицательной области. 

Phi = 43319.26060*arctan(6.232618586*y0)/Pi 

Номер линии

1

2

3

4

5

6

Ф, В*м 

18049.69191 

14439.75353 

10829.81515 

7219.876764 

3609.938377

 0.

Θ, рад

1/6*Pi 

1/3*Pi 

1/2*Pi 

2/3*Pi 

5/6*Pi 

Pi 

Y, м

0.2779009792

0.9263365973e-1

0

-0.9263365973e-1 

-0.2779009792 

R, м 

0.3208924102

0.1852673194

0.1604462051

0.1852673194

0.3208924102 

Выражение потенциала в произвольной точке пространства записывается, как:
Image, где
r- - расстояние от отрицательной электрической оси до точки;
r+ - расстояние от положительной электрической оси до точки;

Подставляя в эту формулу значения координат, получим:
Image, где
x- - координата отрицательно заряженной оси;
x+ - координата положительно заряженной оси;
x,y - координаты точки.
 

phi(x, y) = 1547.116450*ln(((.1508033479-x)^2+y^2)/((-.1700890623-x)^2+y^2))/Pi 

Построение эквипотенциальных поверхностей можно произвести на основе семейства окружностей. Радиусы и центры этих окружностей находятся из условий:
Image, где
Δф - шаг следования эквипотенциальных поверхностей;
ф- - значение потенциала на поверхности отрицательно заряженного цилиндра;
ф+ - значение потенциала на поверхности положительно заряженного цилиндра.

Учитывая следствия из теоремы Гаусса, построим чертёж, с нанесёнными силовыми линиями и эквипотенциальными поверхностями:
 

Х, м 

R, м 

ф, В 

-0.2003571431

-0.2193836084
-0.2497248338
-0.3014858424
-0.4008446235
-0.6433766393
-1.880357181
1.921802819
0.6478943580
0.4024105440
0.3022410546
0.2501490641
0.2196428571 

0.1200000004

0.1496201287

0.1913622426

0.2552464074
0.3673328565
0.6230493685
1.873499438
1.915093494
0.6277134015 

0.3690410021
0.2561379908
0.1919155271
0.1499999998 

1084.246054
917.5793870 

750.9127202
584.2460535
417.5793867
250.9127201
84.24605342
-82.42061323
-249.0872799
-415.7539464
-582.4206131
-749.0872799
-915.7539466 



3. Вычислим и представим в виде графиков изменение потенциала и напряжённости элктростатического поля на линии, соединяющей наиболее близкие точки цилиндров, и распределение плотности заряда на поверхности цилиндра меньшего радиуса.

График изменения потенциала на линии A-B (функция Image):
 

phi(x, 0) = 1547.116450*ln((.1508033479-x)^2/(-.1700890623-x)^2)/Pi 

Plot 

Аналитическую функцию для напряжённости поля найдём из следующих соображений (на линии A-B):
Image

График напряжённости поля на отрезке A-B:
 

E(x, 0) = 3094.232900*(1/(.1700890623+x)-1/(-.1508033479+x))/Pi 

Plot 

Построим зависимость плотности заряда на положительном цилиндре от угла поворота радиус-вектора, соединяющего точку на цилиндре и его геометрическую ось.
Вектор напряжённости поля, создаваемого двумя заряженными нитями: 

Image

Координаты x,y определяются параметрическими выражениями:
Image, где

α - вышеуказанный угол;

R+ - радиус положительно заряженного цилиндра;

Xc+ - координата по оси X центра цилиндра.
 (Typesetting:-mprintslash)([E(x, y) = Vector[row]([3094.232900*((.1700890623+x)/((.1700890623+x)^2+y^2)-(-.1508033479+x)/((-.1508033479+x)^2+y^2))/Pi, 3094.232900*(y/((.1700890623+x)^2+y^2)-y/((-.1508...
(Typesetting:-mprintslash)([E(x, y) = Vector[row]([3094.232900*((.1700890623+x)/((.1700890623+x)^2+y^2)-(-.1508033479+x)/((-.1508033479+x)^2+y^2))/Pi, 3094.232900*(y/((.1700890623+x)^2+y^2)-y/((-.1508...
(Typesetting:-mprintslash)([E(x, y) = Vector[row]([3094.232900*((.1700890623+x)/((.1700890623+x)^2+y^2)-(-.1508033479+x)/((-.1508033479+x)^2+y^2))/Pi, 3094.232900*(y/((.1700890623+x)^2+y^2)-y/((-.1508...
(Typesetting:-mprintslash)([E(x, y) = Vector[row]([3094.232900*((.1700890623+x)/((.1700890623+x)^2+y^2)-(-.1508033479+x)/((-.1508033479+x)^2+y^2))/Pi, 3094.232900*(y/((.1700890623+x)^2+y^2)-y/((-.1508...
(Typesetting:-mprintslash)([E(x, y) = Vector[row]([3094.232900*((.1700890623+x)/((.1700890623+x)^2+y^2)-(-.1508033479+x)/((-.1508033479+x)^2+y^2))/Pi, 3094.232900*(y/((.1700890623+x)^2+y^2)-y/((-.1508...
(Typesetting:-mprintslash)([E(x, y) = Vector[row]([3094.232900*((.1700890623+x)/((.1700890623+x)^2+y^2)-(-.1508033479+x)/((-.1508033479+x)^2+y^2))/Pi, 3094.232900*(y/((.1700890623+x)^2+y^2)-y/((-.1508... 

x = -.2100000000+.12*cos(alpha) 

y = .12*sin(alpha) 

(Typesetting:-mprintslash)([E(x, y) = Vector[row]([3094.232900*((-0.399109377e-1+.12*cos(alpha))/((-0.399109377e-1+.12*cos(alpha))^2+0.144e-1*sin(alpha)^2)-(-.3608033479+.12*cos(alpha))/((-.3608033479...
(Typesetting:-mprintslash)([E(x, y) = Vector[row]([3094.232900*((-0.399109377e-1+.12*cos(alpha))/((-0.399109377e-1+.12*cos(alpha))^2+0.144e-1*sin(alpha)^2)-(-.3608033479+.12*cos(alpha))/((-.3608033479...
(Typesetting:-mprintslash)([E(x, y) = Vector[row]([3094.232900*((-0.399109377e-1+.12*cos(alpha))/((-0.399109377e-1+.12*cos(alpha))^2+0.144e-1*sin(alpha)^2)-(-.3608033479+.12*cos(alpha))/((-.3608033479...
(Typesetting:-mprintslash)([E(x, y) = Vector[row]([3094.232900*((-0.399109377e-1+.12*cos(alpha))/((-0.399109377e-1+.12*cos(alpha))^2+0.144e-1*sin(alpha)^2)-(-.3608033479+.12*cos(alpha))/((-.3608033479...
(Typesetting:-mprintslash)([E(x, y) = Vector[row]([3094.232900*((-0.399109377e-1+.12*cos(alpha))/((-0.399109377e-1+.12*cos(alpha))^2+0.144e-1*sin(alpha)^2)-(-.3608033479+.12*cos(alpha))/((-.3608033479...
(Typesetting:-mprintslash)([E(x, y) = Vector[row]([3094.232900*((-0.399109377e-1+.12*cos(alpha))/((-0.399109377e-1+.12*cos(alpha))^2+0.144e-1*sin(alpha)^2)-(-.3608033479+.12*cos(alpha))/((-.3608033479...
(Typesetting:-mprintslash)([E(x, y) = Vector[row]([3094.232900*((-0.399109377e-1+.12*cos(alpha))/((-0.399109377e-1+.12*cos(alpha))^2+0.144e-1*sin(alpha)^2)-(-.3608033479+.12*cos(alpha))/((-.3608033479...
(Typesetting:-mprintslash)([E(x, y) = Vector[row]([3094.232900*((-0.399109377e-1+.12*cos(alpha))/((-0.399109377e-1+.12*cos(alpha))^2+0.144e-1*sin(alpha)^2)-(-.3608033479+.12*cos(alpha))/((-.3608033479... 

abs(E(x, y)) = (9574277.239*((-0.399109377e-1+.12*cos(alpha))/((-0.399109377e-1+.12*cos(alpha))^2+0.144e-1*sin(alpha)^2)-(-.3608033479+.12*cos(alpha))/((-.3608033479+.12*cos(alpha))^2+0.144e-1*sin(alp...
abs(E(x, y)) = (9574277.239*((-0.399109377e-1+.12*cos(alpha))/((-0.399109377e-1+.12*cos(alpha))^2+0.144e-1*sin(alpha)^2)-(-.3608033479+.12*cos(alpha))/((-.3608033479+.12*cos(alpha))^2+0.144e-1*sin(alp...
abs(E(x, y)) = (9574277.239*((-0.399109377e-1+.12*cos(alpha))/((-0.399109377e-1+.12*cos(alpha))^2+0.144e-1*sin(alpha)^2)-(-.3608033479+.12*cos(alpha))/((-.3608033479+.12*cos(alpha))^2+0.144e-1*sin(alp...
abs(E(x, y)) = (9574277.239*((-0.399109377e-1+.12*cos(alpha))/((-0.399109377e-1+.12*cos(alpha))^2+0.144e-1*sin(alpha)^2)-(-.3608033479+.12*cos(alpha))/((-.3608033479+.12*cos(alpha))^2+0.144e-1*sin(alp...
abs(E(x, y)) = (9574277.239*((-0.399109377e-1+.12*cos(alpha))/((-0.399109377e-1+.12*cos(alpha))^2+0.144e-1*sin(alpha)^2)-(-.3608033479+.12*cos(alpha))/((-.3608033479+.12*cos(alpha))^2+0.144e-1*sin(alp...
abs(E(x, y)) = (9574277.239*((-0.399109377e-1+.12*cos(alpha))/((-0.399109377e-1+.12*cos(alpha))^2+0.144e-1*sin(alpha)^2)-(-.3608033479+.12*cos(alpha))/((-.3608033479+.12*cos(alpha))^2+0.144e-1*sin(alp...
abs(E(x, y)) = (9574277.239*((-0.399109377e-1+.12*cos(alpha))/((-0.399109377e-1+.12*cos(alpha))^2+0.144e-1*sin(alpha)^2)-(-.3608033479+.12*cos(alpha))/((-.3608033479+.12*cos(alpha))^2+0.144e-1*sin(alp... 

sigma = epsilon[r]*epsilon[0]*E 

sigma = 0.8850000000e-11*(9574277.239*((-0.399109377e-1+.12*cos(alpha))/((-0.399109377e-1+.12*cos(alpha))^2+0.144e-1*sin(alpha)^2)-(-.3608033479+.12*cos(alpha))/((-.3608033479+.12*cos(alpha))^2+0.144e...
sigma = 0.8850000000e-11*(9574277.239*((-0.399109377e-1+.12*cos(alpha))/((-0.399109377e-1+.12*cos(alpha))^2+0.144e-1*sin(alpha)^2)-(-.3608033479+.12*cos(alpha))/((-.3608033479+.12*cos(alpha))^2+0.144e...
sigma = 0.8850000000e-11*(9574277.239*((-0.399109377e-1+.12*cos(alpha))/((-0.399109377e-1+.12*cos(alpha))^2+0.144e-1*sin(alpha)^2)-(-.3608033479+.12*cos(alpha))/((-.3608033479+.12*cos(alpha))^2+0.144e...
sigma = 0.8850000000e-11*(9574277.239*((-0.399109377e-1+.12*cos(alpha))/((-0.399109377e-1+.12*cos(alpha))^2+0.144e-1*sin(alpha)^2)-(-.3608033479+.12*cos(alpha))/((-.3608033479+.12*cos(alpha))^2+0.144e...
sigma = 0.8850000000e-11*(9574277.239*((-0.399109377e-1+.12*cos(alpha))/((-0.399109377e-1+.12*cos(alpha))^2+0.144e-1*sin(alpha)^2)-(-.3608033479+.12*cos(alpha))/((-.3608033479+.12*cos(alpha))^2+0.144e...
sigma = 0.8850000000e-11*(9574277.239*((-0.399109377e-1+.12*cos(alpha))/((-0.399109377e-1+.12*cos(alpha))^2+0.144e-1*sin(alpha)^2)-(-.3608033479+.12*cos(alpha))/((-.3608033479+.12*cos(alpha))^2+0.144e... 



Построим зависимость в декартовых координатах:
 

Plot 



Построим зависимость в полярных координатах:
 

Plot