Добрый день.
Вот вообщем суть проблемы: Веду расчет поляризованного электромагнитного привода в Femm, задача осисемметричная. Как известно электромагнитная сила представляет собой производную энергии магнитного поля по зазору (F = dW/dx), в описании к Femm указано, что электромагнитная сила вычисляется как разность коэнергии в текущем положении и коэнергии на предыдущем положении деленная разность текущем положением и предыдущем положением якоря (вообщем величина на которую якорь переместился).

Путем использования блока интеграции я определил электромагнитную силу действующую на якорь в текущем положении, а также с помощью того же блока определил коэнергию на текущем и предыдущем положении якоря и согласно вышеизложенного алгоритма вычислил электромагнитную силу. Оказалось что они не соотвествую друг другу, как же так??? Пытался варьировать величиной на который перемещается якорь, результат тот же - силы разные. Может быть кто нибудь сталкивался с подобной проблемой.......
Заранее благодарен.

Может, программа правильно делает, использую предельный переход? А Вы получаете среднюю силу.
Кто-то из вас ошибается...

Добрый день.
Немного не понимаю ваш вопрос , не могли бы Вы подробней пояснить.


На счет средней силы я задумывался, ну тогда бы на мой взгляд уменьшение величины на которую перемещается якорь вело бы к уменьшению расхождения между получаемыми значениями электромагнитной силы, разве не так???

Это так. А какое у Вас расхождение? В разы? Или проценты?

В первоначальном положении сила составляет 644 Н, коэнергия 0,946552 Дж, перемещаю якорь на 0,01мм соотвественно сила равна 630 Н, коэнергия 0,942165 Дж.
Вычисляю силу через коэнергию (0,946552 - 0,942165)/0,00001 = 438,7 Н, относительная погрешность равна 32%