MathCAD - Ошибка "Эту функцию здесь использовать нельзя" (rfixed), Нахождение токов в электрической цепи по диф. уравнениям Опции

[Andbiz]

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

Может ли кто-нибудь найти ошибку в выложенном мною файле?

 Calculation3.rar ( 65.73 килобайт ) Кол-во скачиваний: 266

[Andbiz]

Сейчас дошел до следующей проблемы. Уравнения упростил. Необходимо добиться, чтобы маткад посчитал диф.уравнения. Пробовал применять кроме функции rkfixed функцию Bulstoer, но считает очень долго и требует уменьшить шаг расчета. Уменьшаю до максимума, но все равно не получаю результата.
Может кто-нибудь взглянуть - в чем ошибка и как их можно посчитать?
 __________._______________.rar ( 17.9 килобайт ) Кол-во скачиваний: 83

[fider]

Может ли кто-нибудь найти ошибку в выложенном мною файле?

Рискну предположить, что проблема в исходных данных - их сочетание вызывает расчетные трудности.

1. Судя по е(t)=37sin(11t) период колебаний, навязываемых источником - 0,57 секунды.
2. С другой стороны, величины L и C Такие, что резонансная частота должна быть 0,3 МГц при периоде 3*10е-6 секунды.
3. Однако сопротивление R1 = 6,25 кОм - процесс для собственных колебаний явно апериодический.

То есть для расчета хотя бы трех колебаний на частоте е(t) надо брать интервал времени 1,5...2 секунды.
Но из-за переходного процесса в LC цепи вроде бы надо брать интервал 10е-7 секунды.

Формулы посмотрел на бегу - похоже на правду. И второй вариант - не с контурными токами, а с общим током и напряжением на C.

Прикрепляю свои версии файлов с результатими расчета для предложенных данных и еще дополнительно со своими L, C, R1.

 Calculation3_modified.zip ( 77.17 килобайт ) Кол-во скачиваний: 98

 Calculation3_modified_2.zip ( 61.97 килобайт ) Кол-во скачиваний: 78

 RTF.zip ( 79.83 килобайт ) Кол-во скачиваний: 61


Сообщение отредактировал fider - Nov 25 2014, 03:43

[Andbiz]

Рискну предположить, что проблема в исходных данных - их сочетание вызывает расчетные трудности.

1. Судя по е(t)=37sin(11t) период колебаний, навязываемых источником - 0,57 секунды.
2. С другой стороны, величины L и C Такие, что резонансная частота должна быть 0,3 МГц при периоде 3*10е-6 секунды.
3. Однако сопротивление R1 = 6,25 кОм - процесс для собственных колебаний явно апериодический.

То есть для расчета хотя бы трех колебаний на частоте е(t) надо брать интервал времени 1,5...2 секунды.
Но из-за переходного процесса в LC цепи вроде бы надо брать интервал 10е-7 секунды.

Формулы посмотрел на бегу - похоже на правду. И второй вариант - не с контурными токами, а с общим током и напряжением на C.

Прикрепляю свои версии файлов с результатими расчета для предложенных данных и еще дополнительно со своими L, C, R1.

 Calculation3_modified.zip ( 77.17 килобайт ) Кол-во скачиваний: 98

 Calculation3_modified_2.zip ( 61.97 килобайт ) Кол-во скачиваний: 78

 RTF.zip ( 79.83 килобайт ) Кол-во скачиваний: 61

Спасибо, посмотрел. Проблему решил иначе. Прочитал про функцию BDF. Она смогла посчитать значения емкости и индуктивности с номиналами нано и микро.