Wrote

... drawbacks. In linux smath does not work very well. I have all calculations going on the right side of de system. I calculate the coordinates of the tangent point and respectively the rotation speed of the oval. At the tangent point the velocity is zero and a tick icon is observed. ...

Hi grelkin2. I don't have a linux rigth now, and can tell which issue could be that. About the numeric precission, there are two big sorces of numerical errors. First, the arc length, which was "calculated" by eye looking at the graph when the curve closes, given the value of 7.3. This value can be improved numerically by solving an equation that sets the square of the norm between the final and initial points equal to zero. This problem is not difficult to solve.

As comment, this is why I notate in other posts under this topic that one need to divide with the norme of the Dragile system to parametrize about the arc lenght.

The second major source of error is considering the minimum of the entire curve only by taking into account a few points: only 31 for each rotated figure, from which the minimum is selected at each iteration. It could be improved by considering many more points, but it would be very slow and impractical, so it would be necessary to interpolate assuming a smooth curve, but it would require many steps to achieve this.
Finally, rkfixed is used, which is the least accurate of the wide variety of ode solvers available.

Best regards.
Alvaro.

Wrote

Hi grelkin2. I don't have a linux rigth now, and can tell which issue could be that. About the numeric precission, there are two big sorces of numerical errors. First, the arc length, which was "calculated" by eye looking at the graph when the curve closes, given the value of 7.3. This value can be improved numerically by solving an equation that sets the square of the norm between the final and initial points equal to zero. This problem is not difficult to solve.

As comment, this is why I notate in other posts under this topic that one need to divide with the norme of the Dragile system to parametrize about the arc lenght.

The second major source of error is considering the minimum of the entire curve only by taking into account a few points: only 31 for each rotated figure, from which the minimum is selected at each iteration. It could be improved by considering many more points, but it would be very slow and impractical, so it would be necessary to interpolate assuming a smooth curve, but it would require many steps to achieve this.
Finally, rkfixed is used, which is the least accurate of the wide variety of ode solvers available.

Best regards.
Alvaro.

Arc length is about 7.335. I changed N is your sm file to 121 points for smooth curves. rkfixed is a reliable basic method, the other methods are needed for optimization or experimentation or have specific uses.
Below is a plot of the translational velocity of the oval at which the velocity at the tangent point is zero.

Hi grelkin2

Wrote

Arc length is about 7.335. ...

This is how Dragilev's method can be used to calculate the arc length of a closed curve numerically. My improvement on Dragilev's method consists in dividing the system of differential equations by its norm to obtain the parameterization of the curves along their arc length. I have not seen this observation anywhere else.



Arc length of a planar closed curve - Dragilev Method.sm (29 KiB) downloaded 51 time(s).

Wrote

... I changed N is your sm file to 121 points for smooth curves. rkfixed is a reliable basic method, the other methods are needed for optimization or experimentation or have specific uses. ...

If you say ...

Wrote

Below is a plot of the translational velocity of the oval at which the velocity at the tangent point is zero.

Yes, like my approach with 31 points, except that your zero is more zero than mine, which is only an approximation, but a pretty good one, because to detect it you have to zoom in on the graph. There is no numerical method that can withstand a sufficiently large zoom on its graph. It looks like you use more than 3000 points, something SMath couldn't handle in this case. One question: what's the interpration of the values of the x-axis?

Best regards.
Alvaro.

You can configure the method using a rkfixed variable. rkfixed := "RK1", …, rkfixed := “RK6”.
In Mathcad, it seems to me this method with a variable step. When I figure it out as this is done, I will make an update.

Wrote

You can configure the method using a rkfixed variable. rkfixed := "RK1", …, rkfixed := “RK6”.
In Mathcad, it seems to me this method with a variable step. When I figure it out as this is done, I will make an update.

Будет замечательно если ты подготовишь модель и рассмотришь зависимость времени вычисления от порядка при заданной точности.

It will be great if you prepare a model and consider the dependence of the calculation time on the order at a given accuracy.

Wrote

Yes, like my approach with 31 points, except that your zero is more zero than mine, which is only an approximation, but a pretty good one, because to detect it you have to zoom in on the graph. There is no numerical method that can withstand a sufficiently large zoom on its graph. It looks like you use more than 3000 points, something SMath couldn't handle in this case. One question: what's the interpration of the values of the x-axis?

Best regards.
Alvaro.

You could have sufficiently large zoom in octave.
olv2.mov (161 KiB) downloaded 53 time(s).
Value of the x-axis is the number of points at a given accuracy. I changed the accuracy and, accordingly, the number and combined everything in one animation.
olv3.mov (365 KiB) downloaded 53 time(s).

Wrote

You could have sufficiently large zoom in octave.
olv2.mov (161 KiB) downloaded 53 time(s).

If 1000x is enough zoom for you, and make this point zero, it is ok for me too.



With and without sliding.





Rolling Curve - Dragilev Method.sm (41 KiB) downloaded 61 time(s).
Rolling Curve - Dragilev Method.pdf (262 KiB) downloaded 62 time(s).

Best regards.
Alvaro.

Wrote

Best regards.
Alvaro.

[RUS]В этих замечательных анимациях можно наблюдать движение замкнутых кривых трех видов в зависимости от соотношения скорости движения центра и скорости вращения в точке касания. Движение с пробуксовкой, движение без проскальзывания и движение с проскальзыванием получаются при скорости движения центра меньше скорости вращения в точке касания, равно и больше соотвественно. Владея методом Драгилева можно задавать любой закон движения.
Ниже анимация с ускорением скорости движения, при которой можно наблюдать плавный переход между тремя видами движения, описываемых раннее. Также в эту анимацию включил движение без проскальзывания с циклоидоподобными кривыми.[/RUS]
[ENG]In these wonderful animations, you can observe the movement of closed curves of three types, depending on the ratio of the speed of movement of the center and the speed of rotation at the point of contact. Movement with stalling, movement without slippage and movement with slippage are obtained when the speed of movement of the center is less than the speed of rotation at the point of contact, equally and more respectively. Using the Draghilev method, you can set any law of motion.
Below is an animation with an acceleration of the speed of movement, at which you can observe a smooth transition between the three types of movement described earlier. Also included in this animation is a non-slip motion with cycloid-like curves.[/ENG]

olvm.mov (1 MiB) downloaded 52 time(s).

Rolling forcing constant speeds.





Rolling Curve on Curve - Dragilev Method.sm (27 KiB) downloaded 68 time(s).

Rolling Curve on Curve - Dragilev Method.pdf (181 KiB) downloaded 76 time(s).

Best regards.
Alvaro.

Yesterday I asked the author of one publication, which mentioned the use of the Method when constructing graphs in Mathcad, to write and publish a scientific paper in a peer-reviewed journal.

1
1`

2
2`

File not found.File not found.
[RUS]Движение эллипса без проскальзывания с поворотом относительно эллипса. Анимация грубая, погрешности большие. Решение получил не точное из системы ду, а комбинированное(сумма нескольких движений)[/RUS]
[ENG]IMovement of the ellipse without slippage with rotation relative to the ellipse. The animation is rough, the errors are large. I got an inaccurate solution from the system of de's, but a combined one (the sum of several movements)[/ENG]
elnel.mov (1 MiB) downloaded 57 time(s).

[RUS]Движение эллипса без проскальзывания с поворотом относительно заданной кривой. Улучшенная способ расчета движения точек эллипса только при помощи системы ду. Видно, что центр эллипса движется по некой кривой, зависимой от эллипса и эквидистанты к неподвижной кривой. При окружности центр бугет двигаться по эквидистанте.[/RUS]
[ENG]Motion of an ellipse without slip with rotation with respect to a given curve. An improved way of calculating the motion of the points of an ellipse using the de system only. It can be seen that the center of the ellipse moves along some curve dependent on the ellipse and equidistant to the fixed curve. For a circle, the center of the ellipse moves along the equidistance.[/ENG]
elnkr.mov (811 KiB) downloaded 55 time(s).

Rolling of a paraboloid clone in its own image.

Directly the picture.

Picture with equations.