Double pendulum model and animation

#1 Posted: 2025/10/29 09:15:19

Valery Ochkov

64 likes in 670 posts.

Group: User

Hi friends!
Now I am on the forum again!
Anime2.sm (55.59 KiB) downloaded 257 time(s).

But the model is not correct. See please Alvaro model!

Edited 2025/10/30 19:22:14

#2 Posted: 2025/10/29 10:45:45

Вячеслав Мезенцев

1,532 likes in 1,838 posts.

Group: Moderator

Using X-Y Plot

Anime2-xyplot.sm (57.84 KiB) downloaded 266 time(s).

Russia ☭ forever, Viacheslav N. Mezentsev

2 users liked this post

sergio 2025/10/29 14:50:17, Ricardo_Herrera 2026/1/3 14:15:17

#3 Posted: 2025/10/29 11:14:52

Valery Ochkov

64 likes in 670 posts.

Group: User

Thanks, Слава!

Edited 2025/10/29 13:53:54

1 users liked this post

Alvaro 2025/10/30 07:04:33

#4 Posted: 2025/10/29 13:56:31

Valery Ochkov

64 likes in 670 posts.

Group: User

Anime-x-y-plot.sm (18.24 KiB) downloaded 248 time(s).
And where here is a bug?

#5 Posted: 2025/10/29 14:01:44

Вячеслав Мезенцев

1,532 likes in 1,838 posts.

Group: Moderator

Vector t must be defined below before the component.
The animation feature is still being developed when I have time for it.

Edited 2025/10/29 14:05:45

Russia ☭ forever, Viacheslav N. Mezentsev

1 users liked this post

Alvaro 2025/10/30 07:04:42

#6 Posted: 2025/10/29 17:02:10

Valery Ochkov

64 likes in 670 posts.

Group: User

Слава, а когда этот новая версия плагина появилась?

#7 Posted: 2025/10/29 18:19:29

Вячеслав Мезенцев

1,532 likes in 1,838 posts.

Group: Moderator

Quote
Слава, а когда этот новая версия плагина появилась?

В конце лета.

Edited 2025/10/29 18:21:02

Russia ☭ forever, Viacheslav N. Mezentsev

#8 Posted: 2025/10/30 07:07:54

Alvaro

1,078 likes in 1,749 posts.

Group: User

Hi. My approach to solving the double pendulum problem.

double_pendulum.sm (75.92 KiB) downloaded 304 time(s).

Best regards.
Alvaro.

3 users liked this post

Вячеслав Мезенцев 2025/10/30 09:13:40, sergio 2025/10/30 12:18:26, Valery Ochkov 2025/10/30 07:45:32

#9 Posted: 2025/10/30 13:32:00

Valery Ochkov

64 likes in 670 posts.

Group: User

Newton vs Lagrange
Pendulum-2.xmcd (119.12 KiB) downloaded 237 time(s).

#10 Posted: 2025/10/30 18:31:11

Valery Ochkov

64 likes in 670 posts.

Group: User

My little edit

1 users liked this post

Alvaro 2025/10/31 04:29:45

#11 Posted: 2025/10/30 19:33:11

Valery Ochkov

64 likes in 670 posts.

Group: User

Triple pendulum - see errors at the second haft of time interval.
Pendulum-3-for-SMath.xmcd (220.62 KiB) downloaded 267 time(s).

Edited 2025/10/30 19:39:16

1 users liked this post

Alvaro 2025/10/31 05:05:16

#12 Posted: 2025/10/31 05:03:31

Alvaro

1,078 likes in 1,749 posts.

Group: User

Hi. Using pure SMath code to deduce the equations of motion from the Lagrangian, with the functions in this post: https://smath.com/en-US/forum/topic/AeeMSM/Using-differentials.

double_pendulum_lagrangian.sm (174.64 KiB) downloaded 245 time(s).

double_pendulum_lagrangian.pdf (503.4 KiB) downloaded 151 time(s).

Best regards.
Alvaro.

2 users liked this post

Valery Ochkov 2025/10/31 07:37:56, sergio 2025/10/31 08:48:37

#13 Posted: 2025/10/31 05:10:58

Alvaro

1,078 likes in 1,749 posts.

Group: User

Wrote
Triple pendulum - see errors at the second haft of time interval.

I no longer use Mathcad, except on rare occasions, when I rely on Mathcad version 11, which I consider the last decent version. So I can't see the triple pendulum file. And I suppose this is the situation for many SMath users who no longer have Mathcad readily available. But more than worrying about the conservation of energy, which seems to be a problem with the numerical method used, I'm more surprised to see how the future affects the past, since the animation shouldn't alter the trajectories already described by the particles as time passes.

This is my approach to the triple pendulum.

triple_pendulum_lagrangian.sm (135.01 KiB) downloaded 242 time(s).

triple_pendulum_lagrangian.pdf (248.09 KiB) downloaded 178 time(s).

Best regards.
Alvaro.

Edited 2025/10/31 06:13:23

2 users liked this post

Valery Ochkov 2025/10/31 07:38:10, Вячеслав Мезенцев 2025/10/31 07:18:43

#14 Posted: 2025/10/31 07:35:30

Valery Ochkov

64 likes in 670 posts.

Group: User

Wrote

I no longer use Mathcad,

See pls pdf
Pendulum-3-for-SMath.pdf (589.25 KiB) downloaded 151 time(s).

Edited 2025/10/31 07:38:52

#15 Posted: 2025/10/31 08:37:45

Alvaro

1,078 likes in 1,749 posts.

Group: User

Wrote
See pls pdf

With those values:

triple_pendulum_lagrangian.sm (159.26 KiB) downloaded 249 time(s).

Best regards.
Alvaro.

1 users liked this post

Valery Ochkov 2025/10/31 10:39:04

#16 Posted: 2025/10/31 10:38:47

Valery Ochkov

64 likes in 670 posts.

Group: User

Wrote

Wrote

I'm more surprised to see how the future affects the past, since the animation shouldn't alter the trajectories already described by the particles as time passes.
Best regards.
Alvaro.

For me too!

#17 Posted: 2025/11/1 01:52:48

Valery Ochkov

64 likes in 670 posts.

Group: User

Wrote

Wrote
See pls pdf

With those values:

Best regards.
Alvaro.

One my plot:
Change of three types of energy of a swinging pendulum on an elastic suspension - SE - the spring energy

Edited 2025/11/1 01:55:02

#18 Posted: 2025/11/5 07:01:02

Valery Ochkov

64 likes in 670 posts.

Group: User

Alvaro, hi!
Is it posibble to solve the single pendulun problem with Lagrangian and rkfixed but without the al_nleqsolve function?

Edited 2025/11/5 08:29:15

#19 Posted: 2025/11/6 04:03:01

Alvaro

1,078 likes in 1,749 posts.

Group: User

Hi Valery. For the single pendulum you can get the lagrangian from the olf post about differentials ( https://smath.com/en-US/forum/topic/AeeMSM/Using-differentials ) It's very simple and doesn't require much processing to use rkfixed.

For the double pendulum, you can use Newton-Raphson to solve the system within the function D. Here, a very simple version of the method is used.

double_pendulum_lagrangian-NR.sm (178.15 KiB) downloaded 229 time(s).

Something similar applies to the pendulum with a spring. Note that the Lagrange equations become somewhat complicated because I use the x,y coordinates to calculate the kinetic energy, and SMath cannot simplify trigonometric functions. However, if you use r, θ' (in green), the equations of motion are greatly simplified. But that breaks the principle that the program should do the calculations, not the user.

Best regards.
Alvaro.

1 users liked this post

Valery Ochkov 2025/11/6 05:16:18

#20 Posted: 2025/11/6 11:47:20

Valery Ochkov

64 likes in 670 posts.

Group: User

Wrote
Hi Valery. For the single pendulum you can get the lagrangian from the olf post about differentials ( https://smath.com/en-US/forum/topic/AeeMSM/Using-differentials ) It's very simple and doesn't require much processing to use rkfixed.
Alvaro.

Hi Alvaro!

Can you create this simple sm-file with animation? Single pendulum!
Thanks!

SMath Forum

Double pendulum model and animation - Is the model correct? - Messages