Применение метода A.B. Драгилева к расчету положений и анимации механизмов - Ber7 - Messages
#21 Posted: 4/21/2015 4:38:26 PM
Если в 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 users liked this post
Andrey Ivashov 11/2/2015 6:27:00 PM
#22 Posted: 10/9/2015 12:59:57 PM
Умер жесткий диск, поэтому перерегистрировался. Небольшое объявление. Предлагаю сотрудничество на github. Метод реализован на с++, вывод в файл и на экран. Буду переходить к 3д и сферическим механизмам.
example.txt
example.txt
2 users liked this post
#23 Posted: 10/20/2015 2:32:38 PM
Попробовал перейти к цилиндрическим координатам, обратный переход делал при помощи функции octave pol2cart().
example.txt
example.txt
1 users liked this post
Andrey Ivashov 11/2/2015 6:27:00 PM
#24 Posted: 10/23/2015 4:51:57 PM
Параллельный робот "Дельта".Кинематическая схема и размеры взяты из книги
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 users liked this post
Davide Carpi 11/22/2015 6:32:00 AM, grelkin2 10/26/2015 2:54:00 PM, Andrey Ivashov 10/26/2015 8:50:00 AM
#25 Posted: 10/26/2015 2:55:58 PM
#26 Posted: 10/27/2015 7:29:33 AM
Вариант с конической спиралью
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 users liked this post
Davide Carpi 11/22/2015 6:32:00 AM
#27 Posted: 10/27/2015 12:12:16 PM
Спираль это хорошо. Сделал демонстарционный пример для декартовой и полярной ск. Странно, что файлы нельзя прикладывать. Тогда только по ссылке.
Файлы для публикации можно хранить в каталоге Wiki,как я сделал выше (Требуется дополнительная регистрация).
Метод расчета рычажных механизмов :
http://www.exponenta.ru/educat/systemat/selitskiy-ivanov/index.asp
#28 Posted: 11/16/2015 2:24:51 PM
Сделал пример с аттрактором Лоренца. 3d система, описываемая д.у. связана с линией или с окружностью. Символьные вычисления делал в wxmaxima, замена автоматически не получалась.
1 users liked this post
Fridel Selitsky 11/16/2015 5:08:00 PM
#29 Posted: 11/21/2015 3:32:47 PM
Заставил wxmaxima делать замену. lisp_3.txt Получил решение для системы из 2-ух точек и стрежнем между ними. Скорость точки должны зависить только от ее координаты или быть постоянной. Можно получать системы, зная только их скорости. D1,D2,D3 - скорости 1-ой точки по оси x,y,z; аналогично и для D3,D4,D5. Пример два аттрактора Лоренца.
#30 Posted: 12/19/2015 1:29:09 PM
Применим метод к математическому маятнику. Составим две эквивалентные системы в Декартовой с.к.
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 Posted: 2/8/2016 2:41:41 PM
Ожерелье из кубиков.
1 users liked this post
Fridel Selitsky 2/8/2016 4:25:00 PM
#32 Posted: 2/9/2016 1:49:09 PM
Точка из примеры выше, подчиняется зсэ. Возможно ли получение общего интеграла из x''+f(x)=0 и использование в методе?
#33 Posted: 2/10/2016 1:06:36 PM
Классический пример с зсэ. Можно наверное сделать и аттрактор с зсэ, но это уже по желанию.
#34 Posted: 2/12/2016 2:29:06 PM
Не забываем о преобразовании скорости при переходе в цск
Уравнения без учета зсэ
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 Posted: 2/20/2016 11:05:22 AM
При помощи метода получается правая часть ду. Она тоже может изучаться методом. Правая часть аттрактора Лоренца небольших степеней, поэтому получены точки бифуркации, собственные числа "точно" в квадратурах в комплексном виде. При степени алгебраического многочлена больше четвертого, тем более для трансцендентного уравнения получение собственных числе в комплексной форме реально только методом. Расширение идеи Алексея Борисовича о равномерном движении позволило решать с сну диф. ур. сохр. эн. При этом делитель(скорость) необходимо разделить на "введенную скорость".
-
New Posts
-
No New Posts