Задания

Задание 4.1

Найти и нарисовать распределение скалярного потенциала электростатического поля между двумя концентрическими цилиндрами (в предположении бесконечно протяженной структуры). Граничные условия: .

а. a = 0.01 м, b = 0.04 м, V0 = 6 В.

б. a = 0.02 м, b = 0.04 м, V0 = 6 В.

в. a = 0.02 м, b = 0.08 м, V0 = 6 В.

г. a = 0.01 м, b = 0.04 м, V0 = 9 В.

д. a = 0.02 м, b = 0.04 м, V0 = 9 В.

е. a = 0.02 м, b = 0.08 м, V0 = 9 В.

ж. a = 0.01 м, b = 0.04 м, V0 = 12 В.

з. a = 0.02 м, b = 0.04 м, V0 = 12 В.

и. a = 0.02 м, b = 0.08 м, V0 = 12 В.

Решение (стр.197, пример 4.7)

Краевая задача формулируется как задача Дирихле (4.9). В силу симметрии решаем уравнение Лапласа в цилиндрической системе координат

. (1)

Показать, что , где .

В силу симметрии и . В результате (1) принимает вид:

. (2)

Решая последовательно это уравнение имеем:

, (3)

где константы определяются из граничных условий

(4)

Откуда . Соответственно из (3) потенциал электростатического поля между двумя концентрическими цилиндрами равен

, где . (5)

Используя MATLAB строим распределение скалярного потенциала электростатического поля между двумя концентрическими цилиндрами, полагая, что .

clear;

clf

a = 0.01;

b = 0.04;

V0 = 6;

radius = a: 0.001: b;

voltage = V0 * log(radius / b) ./ log(a / b);

plot (radius, voltage)

xlabel ('\rho')

ylabel ('V')

grid on

Задание 4.2

Двухмерное распределение потенциала аппроксимируется выражением . Показать, что эта функция удовлетворяет уравнению Пуассона и построить поверхности распределения потенциала и электрического заряда , Кл/м³. Поверхности строятся в интералах , с шагом 0.2 м.

а.

б.

в.

г.

д.

е.

ж.

Решение (стр.190, пример 4.4)

Так как , то заданное распределение потенциала удовлетворяет двухмерному уравнению Пуассона .

Используя MATLAB строим распределение скалярного потенциала электростатического поля u и электрического заряда .

clear;

clf

xmin = -5;

xmax = 5;

ymin = -5;

ymax = 5;

step = 0.2;

[x, y] = meshgrid(xmin: step: xmax, xmin: step: xmax);

voltage = (x .^ 2 + y .^ 2) / 4;

surf(x, y, voltage)

hold on

charge = del2(voltage);

mesh (x, y, charge)

view(-37.5, 20)

xlabel ('x', 'fontsize', 12)

ylabel ('y', 'fontsize', 12)

zlabel ('\rho_v, V(x, y)', 'fontsize', 12)

Задание 4.3

Двухмерное распределение потенциала в районе и описывается выражением . Продифференцировать аналитически, а также, используя MATLAB, реализовать численное дифференцирование уравнения Пуассона при , и построить графики плотности электрических зарядов , Кл/м² , а также распределение потенциала. Графики строятся с шагом 0.1 по осям.

а. =1, V₀ = 1 В.

б. =2, V₀ = 1 В.

в. =3, V₀ = 1 В.

г. =4, V₀ = 1 В.

д. =4, V₀ = 0.1 В.

е. =4, V₀ = 0.5 В.

ж. =4, V₀ = 0.2 В.

з =4, V₀ = 0.3 В.

Решение (стр.224, пример 4.16)

Уравнение Пуассона имеет вид . Использовать в MATLAB функцию del2.

clear;

clf;

V0 = 1;

eps = 1;

xmin = -2;

xmax = 2;

ymin = -2;

ymax = 2;

step = 0.1;

eps_0 = 8.854187817e-12;

eps_a = eps * eps_0;

[x, y] = meshgrid (xmin: step : xmax, ymin: step: ymax);

z = V0 * exp(-x .^ 2 - y .^ 2);

% Умножаем на 4, т.к. del2 делит оператор Лапласа на 2 * N, где N -

% размерность матрицы

ro = -del2(z, step) * 4 * eps_a;

subplot (2, 1, 1)

surf(x, y, z)

view(-45, 10)

xlabel('x')

ylabel('y')

zlabel('u')

subplot (2, 1, 2)

mesh(x, y, ro)

view(-45, 10)

xlabel('x')

ylabel('y')

zlabel('\rho_s')

% Аналитическая проверка

ro_anal = 4 * V0 * exp (-x .^ 2 - y .^ 2) .* (1 - x .^ 2 - y .^ 2) * eps_a;

% Расчет разности между аналитикой и тем, что рассчитывалось с помощью del2

delta = abs (ro - ro_anal);

delta_max = max (max(delta))