Тут для нужд систем магнитной фокусировки взялся я рассчитывать магнитное поле однородно намагниченных цилиндрических магнитов. Причем экспериментаторы мои намерены использовать магниты, намагниченные перпендикулярно оси, но общий случай не сложнее, так что изложу его.
Используя верную вне магнита формулу
и разложение фундаментального решения по функциям Бесселя
(в правой части --- цилиндрические координаты), дошел вот до какого выражения
Здесь радиус цилиндра обозначен R, высота h, ,
--- сводится к набору экспонент.
Видно, что возникает три эталонных интеграла, первый из них
В Прудникове---Брычкове---Маричеве (том 2, раздел 2.12.38) приводится такой ответ
где
верхний знак (и верхний член в фигурных скобках) для , нижний --- для . Определение функции я еле нашел на стр. 662
причем на стр. 657 и 659 расшифрованы стандартные эллиптические интегралы
А вот дальше проблема: при прямом подсчете интеграла для (a,b,c)=(1,1,2) получаем 0.05185272, а при расчете по формуле из ПБМ получаем -0.266674632. При этом другая формула, для интеграла с двумя , правильная. В чем тут может быть дело?
P. S. Тему можно в некотором смысле считать продолжением этой.
Магнитное поле однородно намагниченного цилиндра, как проверить формулу из ПБМ2?
-
- Сообщений: 620
- Зарегистрирован: 29 дек 2015, 13:17
-
- Сообщений: 620
- Зарегистрирован: 29 дек 2015, 13:17
Магнитное поле однородно намагниченного цилиндра, как проверить формулу из ПБМ2?
Я, тем временем, обнаружил, что Maple, похоже, неверно преобразует эллиптические интегралы. Потому что на вопрос посчитать , где
он отвечает значением , тогда как правильный ответ .
Взялся я считать сам. Прежде всего представляем в виде интерала
Тогда
Теперь пользуемся лапласовским образом
в результате чего и приходим к эллиптическому интегралу
Теперь задача в том, чтобы привести это к стандартным эллиптическим интегралам. Подынтегральное выражение имеет точки ветвления , , лежащие, в силу симметрии, на одной окружности. Обозначая нижнюю точку пересечения этой окружности с мнимой осью через , сконструируем дробно-линейное преобразование , которое оставляло бы неподвижными, а переводило бы в бесконечность. Тогда указанная окружность на плоскости переменной будет вещественной прямой, а интеграл приведется к стандартному виду. Указанное преобразование выписывается элементарно
Величину проще всего найти, подставив преобразование в интеграл и потребовав, чтобы под корнем были только четные степени
Сам же интеграл принимает вид
и выражается через стандартные интегралы
здесь
И вот это выражение уже дает правильное значение на Maple.
Только я теперь ума не приложу, как Прудникову и Ко удалось выразить интеграл третьего рода через первый и второй... Или есть возможность как-то по другому провести преобразование интеграла?
он отвечает значением , тогда как правильный ответ .
Взялся я считать сам. Прежде всего представляем в виде интерала
Тогда
Теперь пользуемся лапласовским образом
в результате чего и приходим к эллиптическому интегралу
Теперь задача в том, чтобы привести это к стандартным эллиптическим интегралам. Подынтегральное выражение имеет точки ветвления , , лежащие, в силу симметрии, на одной окружности. Обозначая нижнюю точку пересечения этой окружности с мнимой осью через , сконструируем дробно-линейное преобразование , которое оставляло бы неподвижными, а переводило бы в бесконечность. Тогда указанная окружность на плоскости переменной будет вещественной прямой, а интеграл приведется к стандартному виду. Указанное преобразование выписывается элементарно
Величину проще всего найти, подставив преобразование в интеграл и потребовав, чтобы под корнем были только четные степени
Сам же интеграл принимает вид
и выражается через стандартные интегралы
здесь
И вот это выражение уже дает правильное значение на Maple.
Только я теперь ума не приложу, как Прудникову и Ко удалось выразить интеграл третьего рода через первый и второй... Или есть возможность как-то по другому провести преобразование интеграла?
-
- Сообщений: 620
- Зарегистрирован: 29 дек 2015, 13:17
Магнитное поле однородно намагниченного цилиндра, как проверить формулу из ПБМ2?
Ничего не понимаю... Читаю русскую Википедию, и вижу формулы, очень похожие на Прудникова и Ко. Но в английской версии страницы этих формул нет!
С другой стороны, на английской странице есть ссылка на симметричную форму Карлсона, где утверждается, что все эллиптические интегралы (и неполные тоже) выражаются через два стандартных (а не три)...
P. S. Обнаружил, откуда все списано: Абрамовиц---Стиган, раздел 17.7. Но по прежнему не понял, выражается ли интеграл третьего рода через первый и второй... Может быть, он линейно не выражается?
С другой стороны, на английской странице есть ссылка на симметричную форму Карлсона, где утверждается, что все эллиптические интегралы (и неполные тоже) выражаются через два стандартных (а не три)...
P. S. Обнаружил, откуда все списано: Абрамовиц---Стиган, раздел 17.7. Но по прежнему не понял, выражается ли интеграл третьего рода через первый и второй... Может быть, он линейно не выражается?
-
- Сообщений: 620
- Зарегистрирован: 29 дек 2015, 13:17
Магнитное поле однородно намагниченного цилиндра, как проверить формулу из ПБМ2?
Так, я проверил формулы из Абрамовиц---Стиган. Оказывается, они работают, а интегал 3-го рода выражается через интегралы 1-го и 2-го родов! Я по-прежнему не понимаю, как это получается, но Maple подтверждает это численно.
Более-менее стандартное определение интеграла 3-го рода
Заодно определение интеграла 1-го рода
Аналогично для интеграла 2-го рода.
Справедливы такие формулы.
Для
Для
Для вне интервала существуют формулы приведения к этому интервалу. Например, для интересного нам случая
Более-менее стандартное определение интеграла 3-го рода
Заодно определение интеграла 1-го рода
Аналогично для интеграла 2-го рода.
Справедливы такие формулы.
Для
Для
Для вне интервала существуют формулы приведения к этому интервалу. Например, для интересного нам случая
-
- Сообщений: 620
- Зарегистрирован: 29 дек 2015, 13:17
Магнитное поле однородно намагниченного цилиндра, как проверить формулу из ПБМ2?
Свой правильный ответ у меня есть, но что же делать с ПБМ? Путь к их формуле был долог и труден.
Прежде всего нужно сделать другую подстановку: подгадать так, чтобы одна из точек ветвления, для определенности , отображалась в бесконечность, а еще одна, для определенности --- в нуль. Тогда после преобразования под корнем останется многочлен третьей степени, причем один из корней у него будет . Последующей заменой получим стандартный эллиптический интеграл. В общем, замена теперь имеет вид
После такой замены получается интеграл
Тут уже вырисовываются модуль и параметр
сам же интеграл выражается через стандартные
Однако остаются комплексные числа в формуле, притом что сам результат вещественный. Вот изгнание этих мнимостей и заняло большую часть времени.
Прежде всего я обнаружил, что для нашего работает приведенная выше формула
причем множитель перед скобками равен , а аргумент ---
Интеграл принимает вид
Чтобы выделить отсюда вещественную часть явно, нужно сосчитать , что сводится к суммированию неполных интегралов 1-го и 2-го родов с разными амплитудами. Для таких сумм есть формулы
Однако в нашем случае они напрямую неприменимы, потому что оказывается и получается неопределенность 0/0. Лопиталя ее по , находим формулу для случая
В нашем случае результирующая амплитуда
а . Окончательно
То есть, собственно говоря, в ПБМ лишние квадратные скобки, а так все правильно. Последний член () описывает постоянную добавку к полю внутри магнита, а знак во втором члене описывает ветвление при (загиб контура интегрирования за край разреза).
Прежде всего нужно сделать другую подстановку: подгадать так, чтобы одна из точек ветвления, для определенности , отображалась в бесконечность, а еще одна, для определенности --- в нуль. Тогда после преобразования под корнем останется многочлен третьей степени, причем один из корней у него будет . Последующей заменой получим стандартный эллиптический интеграл. В общем, замена теперь имеет вид
После такой замены получается интеграл
Тут уже вырисовываются модуль и параметр
сам же интеграл выражается через стандартные
Однако остаются комплексные числа в формуле, притом что сам результат вещественный. Вот изгнание этих мнимостей и заняло большую часть времени.
Прежде всего я обнаружил, что для нашего работает приведенная выше формула
причем множитель перед скобками равен , а аргумент ---
Интеграл принимает вид
Чтобы выделить отсюда вещественную часть явно, нужно сосчитать , что сводится к суммированию неполных интегралов 1-го и 2-го родов с разными амплитудами. Для таких сумм есть формулы
Однако в нашем случае они напрямую неприменимы, потому что оказывается и получается неопределенность 0/0. Лопиталя ее по , находим формулу для случая
В нашем случае результирующая амплитуда
а . Окончательно
То есть, собственно говоря, в ПБМ лишние квадратные скобки, а так все правильно. Последний член () описывает постоянную добавку к полю внутри магнита, а знак во втором члене описывает ветвление при (загиб контура интегрирования за край разреза).
Кто сейчас на форуме
Количество пользователей, которые сейчас просматривают этот форум: нет зарегистрированных пользователей и 4 гостей