Применение метода A.B. Драгилева к расчету положений и анимации механизмов - Ber7 - Сообщения
#21 Опубликовано: 21.04.2015 16:38:26
Если в octave взять н.к. в [0 sqrt(2) 0] он тоже зависает. Относительно легко избавился от пологового участка с условием log10(abs(x(:,1)))>-2), этих точек 28%. Сделал 3-ый вариант, с++ без зависания, координаты [1 1 1 2.7321 2.7321 2.7321], тут условие более сложное log10(abs(del_x))>log10(h), del_x=xi+1-xi, h- шаг, доля тоже 28%. На рисунке координаты по x c условием.
1 пользователям понравился этот пост
Andrey Ivashov 02.11.2015 18:27:00
#22 Опубликовано: 09.10.2015 12:59:57
Умер жесткий диск, поэтому перерегистрировался. Небольшое объявление. Предлагаю сотрудничество на github. Метод реализован на с++, вывод в файл и на экран. Буду переходить к 3д и сферическим механизмам.
example.txt
example.txt
2 пользователям понравился этот пост
#23 Опубликовано: 20.10.2015 14:32:38
Попробовал перейти к цилиндрическим координатам, обратный переход делал при помощи функции octave pol2cart().
example.txt
example.txt
1 пользователям понравился этот пост
Andrey Ivashov 02.11.2015 18:27:00
#24 Опубликовано: 23.10.2015 16:51:57
Параллельный робот "Дельта".Кинематическая схема и размеры взяты из книги
Williams II, “The Delta Parallel Robot: Kinematics Solutions”, Internet Publication, April
2015.
http://www.ohio.edu/people/williar4/html/pdf/DeltaKin.pdf
Delta2.smz
Williams II, “The Delta Parallel Robot: Kinematics Solutions”, Internet Publication, April
2015.
http://www.ohio.edu/people/williar4/html/pdf/DeltaKin.pdf
Delta2.smz
Метод расчета рычажных механизмов :
http://www.exponenta.ru/educat/systemat/selitskiy-ivanov/index.asp
3 пользователям понравился этот пост
#25 Опубликовано: 26.10.2015 14:55:58
#26 Опубликовано: 27.10.2015 07:29:33
Вариант с конической спиралью
https://smath.com/wiki/GetFile.aspx?File=Delta3.smz
https://smath.com/wiki/GetFile.aspx?File=Delta3.smz
Метод расчета рычажных механизмов :
http://www.exponenta.ru/educat/systemat/selitskiy-ivanov/index.asp
1 пользователям понравился этот пост
Davide Carpi 22.11.2015 06:32:00
#27 Опубликовано: 27.10.2015 12:12:16
Спираль это хорошо. Сделал демонстарционный пример для декартовой и полярной ск. Странно, что файлы нельзя прикладывать. Тогда только по ссылке.
Файлы для публикации можно хранить в каталоге Wiki,как я сделал выше (Требуется дополнительная регистрация).
Метод расчета рычажных механизмов :
http://www.exponenta.ru/educat/systemat/selitskiy-ivanov/index.asp
#28 Опубликовано: 16.11.2015 14:24:51
Сделал пример с аттрактором Лоренца. 3d система, описываемая д.у. связана с линией или с окружностью. Символьные вычисления делал в wxmaxima, замена автоматически не получалась.
1 пользователям понравился этот пост
Fridel Selitsky 16.11.2015 17:08:00
#29 Опубликовано: 21.11.2015 15:32:47
Заставил wxmaxima делать замену. lisp_3.txt Получил решение для системы из 2-ух точек и стрежнем между ними. Скорость точки должны зависить только от ее координаты или быть постоянной. Можно получать системы, зная только их скорости. D1,D2,D3 - скорости 1-ой точки по оси x,y,z; аналогично и для D3,D4,D5. Пример два аттрактора Лоренца.
#30 Опубликовано: 19.12.2015 13:29:09
Применим метод к математическому маятнику. Составим две эквивалентные системы в Декартовой с.к.
1-ая система
x[0]-x
x[1]-y
x[2]-V_x
x[3]-V_y
x[0]^2+x[1]^2 (1)
9.81*2*x[1]+x[2]^2+x[3]^2 (2) з-н сохр. энергии
x[0]*x[2]+x[1]*x[3] (3) получен из (1) дифференцированием по t
2-ая система
x[0]-x
x[1]-y
x[2]-V
x[0]^2+x[1]^2
9.81*2*x[1]+x[2]^2 з-н сохр. энергии
Сравним скорости, полученные через взятие производной с введенным скоростями. Качествено похожи, отличаются на dt.
Правая часть 1-ой системы:
f[0] = -( 4.0*x[3]*x[0]*x[1]+-4.0*x[2]*(x[1]*x[1]));
f[1] = -( -4.0*x[3]*(x[0]*x[0])+4.0*x[2]*x[0]*x[1]);
f[2] = -( -3.9240000000000002e+01*x[0]*x[1]+-4.0*x[2]*x[3]*x[1]+4.0*(x[3]*x[3])*x[0]);
f[3] = -( 3.9240000000000002e+01*(x[0]*x[0])+4.0*(x[2]*x[2])*x[1]+-4.0*x[2]*x[3]*x[0]);
Видно, что f[0] не равно x[2], f[1] не равно x[3].
mayatnik.mkv.zip
1-ая система
x[0]-x
x[1]-y
x[2]-V_x
x[3]-V_y
x[0]^2+x[1]^2 (1)
9.81*2*x[1]+x[2]^2+x[3]^2 (2) з-н сохр. энергии
x[0]*x[2]+x[1]*x[3] (3) получен из (1) дифференцированием по t
2-ая система
x[0]-x
x[1]-y
x[2]-V
x[0]^2+x[1]^2
9.81*2*x[1]+x[2]^2 з-н сохр. энергии
Сравним скорости, полученные через взятие производной с введенным скоростями. Качествено похожи, отличаются на dt.
Правая часть 1-ой системы:
f[0] = -( 4.0*x[3]*x[0]*x[1]+-4.0*x[2]*(x[1]*x[1]));
f[1] = -( -4.0*x[3]*(x[0]*x[0])+4.0*x[2]*x[0]*x[1]);
f[2] = -( -3.9240000000000002e+01*x[0]*x[1]+-4.0*x[2]*x[3]*x[1]+4.0*(x[3]*x[3])*x[0]);
f[3] = -( 3.9240000000000002e+01*(x[0]*x[0])+4.0*(x[2]*x[2])*x[1]+-4.0*x[2]*x[3]*x[0]);
Видно, что f[0] не равно x[2], f[1] не равно x[3].
mayatnik.mkv.zip
#31 Опубликовано: 08.02.2016 14:41:41
Ожерелье из кубиков.
1 пользователям понравился этот пост
Fridel Selitsky 08.02.2016 16:25:00
#32 Опубликовано: 09.02.2016 13:49:09
Точка из примеры выше, подчиняется зсэ. Возможно ли получение общего интеграла из x''+f(x)=0 и использование в методе?
#33 Опубликовано: 10.02.2016 13:06:36
Классический пример с зсэ. Можно наверное сделать и аттрактор с зсэ, но это уже по желанию.
#34 Опубликовано: 12.02.2016 14:29:06
Не забываем о преобразовании скорости при переходе в цск
Уравнения без учета зсэ
x[1]^2+x[2]^2+cos(2*x[0])
x[0]+10*x[2]
Уравнения без учета зсэ
x[1]^2+x[2]^2+cos(2*x[0])
x[0]+10*x[2]
#35 Опубликовано: 20.02.2016 11:05:22
При помощи метода получается правая часть ду. Она тоже может изучаться методом. Правая часть аттрактора Лоренца небольших степеней, поэтому получены точки бифуркации, собственные числа "точно" в квадратурах в комплексном виде. При степени алгебраического многочлена больше четвертого, тем более для трансцендентного уравнения получение собственных числе в комплексной форме реально только методом. Расширение идеи Алексея Борисовича о равномерном движении позволило решать с сну диф. ур. сохр. эн. При этом делитель(скорость) необходимо разделить на "введенную скорость".
-
Новые сообщения
-
Нет новых сообщений