Solving system of 2nd order coupled ODEs


by guest1234
Tags: coupled, odes, order, solving
guest1234
guest1234 is offline
#1
Nov25-12, 10:00 AM
P: 24
I have to derive equations of motion from Lagrangian and stumbled upon the following system of equations (constants are simplified, that information is unneeded)

[itex]
\begin{cases}
\ddot{x}-A\dot{y}+Bx=0 \\
\ddot{y}+A\dot{x}+Dy=0
\end{cases}
[/itex]

This is an extension of a simpler problem where B=D. There I just multiplied the second equation by i, added equations together, substituted z=x+iy and solved for z.

But doing the same with current system doesn't help much (can't substitute z=x+iy, b/c [itex]B\ne{D}[/itex]).

When choosing E=(B+D)/2 and F=E-B then

[itex]
\begin{cases}
\ddot{x}-A\dot{y}+Ex=Fx\\
\ddot{y}+A\dot{x}+Ey=-Fy
\end{cases}
[/itex]

isn't much help either.

I know there are some techniques for solving coupled differential equations like writing a reciprocal matrix for the system but it seems that it applies only to 1st order ODEs.

Is there any analytic solutions to this?
When punching it into Maple, it throws a huge block of square roots, although simpler problem gave me a combination of exponentials (for z). I know this problem must give similar answer with little modifications, but don't know how to tackle it.

Help and pointers much appreciated
Phys.Org News Partner Science news on Phys.org
SensaBubble: It's a bubble, but not as we know it (w/ video)
The hemihelix: Scientists discover a new shape using rubber bands (w/ video)
Microbes provide insights into evolution of human language
arildno
arildno is offline
#2
Nov25-12, 10:11 AM
Sci Advisor
HW Helper
PF Gold
P: 12,016
Differentiate the first equation, and get an equation for doubledotted y in terms of the x-terms.
Insert this into your second equation, to get y expressed in terms of x-terms.
Differentiate THIS for an expression for singledotted y, and insert this into your first equation.
That is a fourth-order differential equation for x(t)
guest1234
guest1234 is offline
#3
Nov25-12, 02:16 PM
P: 24
Thanks, worked.
But it's strange that equations for x and y have both the same solution (although it shouldn't). I dropped fourth order derivative and solved second order ODE.
Case closed for now.

arildno
arildno is offline
#4
Nov25-12, 02:19 PM
Sci Advisor
HW Helper
PF Gold
P: 12,016

Solving system of 2nd order coupled ODEs


You can't drop the fourth order derivative.
You have no reason to assume it is sufficiently small, relative to the other terms
guest1234
guest1234 is offline
#5
Nov25-12, 02:26 PM
P: 24
Well, depends on the problem (I am using small angle approximation as my problem inherently allows me to do that)
arildno
arildno is offline
#6
Nov25-12, 04:30 PM
Sci Advisor
HW Helper
PF Gold
P: 12,016
Quote Quote by guest1234 View Post
Well, depends on the problem (I am using small angle approximation as my problem inherently allows me to do that)
Not really. The angle might be small, even though the higher rates of changes aren't.

If you have positive B, for example, you will have two vibrational modes, not a single harmonic motion.

Only if (1+A^2+B)^2-4B is negligible, relative to (1+A^2+B) can you make that simplification.
guest1234
guest1234 is offline
#7
Nov25-12, 06:23 PM
P: 24
Context: problem I need to solve is about pendulum in a rotating frame
Well, I already dropped all components that are Big O(angle ^3) and beyond in Lagrangian, for the sake of finding analytic solution (otherwise I'd have system of three coupled 2nd order diff equation with variable coefficient, no good). I believe that including further analytic accuracy won't weigh up real life effects, although interesting effects may occour.
I already included centrifugal, Coriolis and Eötvös effects (the reason I was surprised that equations of motion are the same for x and y; will redo my math tomorrow). Two vibrational modes are kind of things that are beyond me. Is it something like oscillating around its trajectory or..? Eager to know.
arildno
arildno is offline
#8
Nov25-12, 06:37 PM
Sci Advisor
HW Helper
PF Gold
P: 12,016
IGNORE THIS POST!!!
It means you have two distinct angular frequencies.
(Your four solutions is therefore a linear combination of the sines and cosines of those two frequencies)

Assuming B>0, those frequencies will be:
[tex]\omega_{1}=\sqrt{\frac{(1+A^{2}+B)+\sqrt{(1+A^{2}+B)^{2}-4B}}{2}}[/tex]
[tex]\omega_{2}=\sqrt{\frac{(1+A^{2}+B)-\sqrt{(1+A^{2}+B)^{2}-4B}}{2}}[/tex]


OOPS!
I just see that in my mental calulation here, I managed to treat A=D
That makes my solution all awry, but my point about the two frequencies is still valid.
Now, it is night time here in Oslo, so I'll get back to this tomorrow.
arildno
arildno is offline
#9
Nov25-12, 07:03 PM
Sci Advisor
HW Helper
PF Gold
P: 12,016
Couldn't sleep yet, so here we go:
The first equation yields us:
[tex]A\ddot{y}=\dddot{x}+B\dot{x}[/tex]
Thus, the second equation yields:
[tex]ADy=-\dddot{x}-(B+A^{2})\dot{x}[/tex]
so that:
[tex]-A\dot{y}=\frac{1}{D}(\ddddot{x}+(B+A^{2})\ddot{x})[/tex]
Thus, the first equation can be rewritten as:
[tex]\ddddot{x}+(B+D+A^{2})\ddot{x}+BDx=0[/tex]

An exponential trial solution [itex]x(t)=Ke^{rt}[/itex] yields the characteristic equation:
[tex]r^{4}+(B+D+A^{2})r^{2}+BD=0[/tex]
Or, solving for r squared:
[tex]r^{2}=\frac{-(B+D+A^{2})\pm\sqrt{(B+D+A^{2})^{2}-4BD}}{2}[/tex]
(my error was D=1, not D=A).
Now, note that for non-negative B and D, r^2 is necessarily negative (and the innermost radicand can be shown to be positive).
the plus/minus solutions yield the two angular frequencies, when finding imaginary solutions r:
[tex]\omega_{1}=\sqrt{\frac{(B+D+A^{2})+\sqrt{(B+D+A^{2})^{2}-4BD}}{2}}[/tex]
[tex]\omega_{2}=\sqrt{\frac{(B+D+A^{2})-\sqrt{(B+D+A^{2})^{2}-4BD}}{2}}[/tex]

You may also rewrite these two frequencies as:
[tex]\omega_{1}=\frac{\hat{\omega}}{\sqrt{2}}\sqrt{1+{\sqrt{1-\epsilon}}}, \hat{\omega}=\sqrt{(B+D+A^{2})}, \epsilon=\frac{4BD}{(B+D+A^{2})^{2}}[/tex]
[tex]\omega_{2}=\frac{\hat{\omega}}{\sqrt{2}}\sqrt{1-\sqrt{1-\epsilon}}[/tex]

----------------------
Your general solution for x(t) is, therefore, up to 4 constants determined by initial conditions:
[tex]x(t)=C_{1}\cos(\omega_{1}t)+C_{2}\sin(\omega_{1}t)+C_{3}\cos(\omega_{2} t)+C_{4}\sin(\omega_{2}t)[/tex]
----------------------------------
Note that in the limit, [itex]\epsilon\to{1}[/itex], we get a single frequency,
[tex]\omega=\frac{\hat{\omega}}{\sqrt{2}}[/tex], whereas
in the limit [itex]\epsilon\to{0}[/itex], we get TWO frequencies,
[tex]\omega_{1}=\hat{\omega}[/tex]
[tex]\omega_{2}=\sqrt{\frac{BD}{B+D+A^{2}}}[/tex]
By omitting the fourth order derivative, you lose the high frequency, [itex]\omega_{1}[/itex]
-----
That your problem should have TWO fundamental frequencies, rather than just one, ought to be obvious, because you have two independent sources of "periodic" motion: The pendulum's own frequency of oscillation, and the rotating frame's rotation rate.
arildno
arildno is offline
#10
Nov26-12, 06:23 AM
Sci Advisor
HW Helper
PF Gold
P: 12,016
Having nothing better to do today, I'll derive the full solution for x(t)
1. Setting:
[tex]x(0)=x^{(0)}[/tex]
[tex]\dot{x}(0)=x^{(1)}[/tex]
[tex]\ddot{x}(0)=x^{(2)}[/tex]
[tex]\dddot{x}(0)=x^{(3)}[/tex]
2. We gain, by differentiating thrice the general solution, evaluating at t=0:
[tex]x^{(0)}=C_{1}+C_{3}[/tex]
[tex]x^{(2)}=-\omega_{1}^{2}C_{1}-\omega_{2}^{2}C_{3}[/tex]
[tex]x^{(1)}=\omega_{1}C_{2}+\omega_{2}C_{4}[/tex]
[tex]x^{(3)}=-\omega_{1}^{3}C_{2}-\omega_{2}^{3}C_{4}[/tex]
3. Therefore:
[tex]C_{1}=\frac{\omega_{2}^{2}x^{(0)}+x^{(2)}}{\omega_{2}^{2}-\omega_{1}^{2}}[/tex]
[tex]C_{3}=\frac{\omega_{1}^{2}x^{(0)}+x^{(2)}}{\omega_{1}^{2}-\omega_{2}^{2}}[/tex]
[tex]C_{2}=\frac{\omega_{2}^{2}x^{(1)}+x^{(3)}}{\omega_{1}(\omega_{2}^{2}-\omega_{1}^{2})}[/tex]
[tex]C_{4}=\frac{\omega_{1}^{2}x^{(1)}+x^{(3)}}{\omega_{2}(\omega_{1}^{2}-\omega_{2}^{2})}[/tex]
4. Simplifications:
We have:
a) [tex]\omega_{1}^{2}-\omega_{2}^{2}=\hat{\omega}^{2}\sqrt{1-\epsilon}[/tex]
b)
[tex]x^{(2)}=Ay^{(1)}-Bx^{(0)}[/tex]
c)
[tex]x^{(3)}=-(A^{2}+B)x^{(1)}-ADy^{(0)}[/tex]
5. That is:
[tex]C_{1}=-\frac{(\omega_{2}^{2}-B)x^{(0)}+Ay^{(1)}}{\hat{\omega}^{2}\sqrt{1-\epsilon}}[/tex]
[tex]C_{3}=\frac{(\omega_{1}^{2}-B)x^{(0)}+Ay^{(1)}}{\hat{\omega}^{2}\sqrt{1-\epsilon}}[/tex]
[tex]C_{2}=\frac{(\omega_{2}^{2}-A^{2}-B)x^{(1)}+ADy^{(0)}}{\omega_{1}\hat{\omega}^{2}{\sqrt{1-\epsilon}}}[/tex]
[tex]C_{4}=-\frac{(\omega_{1}^{2}-A^{2}-B)x^{(1)}+ADy^{(0)}}{\omega_{2}\hat{\omega}^{2}{\sqrt{1-\epsilon}}}[/tex]
Collected together, this is the solution for x(t), relative to the 7 arbitrary constants the problem includes.
guest1234
guest1234 is offline
#11
Dec1-12, 09:45 AM
P: 24
Thank you!

I did my math by myself and verify your solution, although I searched solution in the form [itex]x(t)=x_0e^{i\omega{t}}[/itex], so it's just a matter of adding an imaginary unit here and there. It's now clear to me, why one cannot drop fourth derivative - when dealing with rotating reference frames, the solution must include two different frequencies.

For now, I must find conditions for coefficients in equation of motion that satisfy small angle approximation, but that's a whole different topic. Some questions have arisen tho, but I will formulate them after deadline in another topic.
arildno
arildno is offline
#12
Dec1-12, 10:36 AM
Sci Advisor
HW Helper
PF Gold
P: 12,016
Quote Quote by guest1234 View Post
Thank you!

I did my math by myself and verify your solution, although I searched solution in the form [itex]x(t)=x_0e^{i\omega{t}}[/itex], so it's just a matter of adding an imaginary unit here and there. It's now clear to me, why one cannot drop fourth derivative - when dealing with rotating reference frames, the solution must include two different frequencies.
You're welcome!


Register to reply

Related Discussions
Solving a system of first order ODEs Differential Equations 0
Solving a system of two nonlinear second order ODEs (Mechanical vibrations) Differential Equations 6
Coupled Second order odes Differential Equations 0
Solving Coupled System of ODEs in MATLAB Engineering, Comp Sci, & Technology Homework 2
Solving coupled ODEs with Matlab Math & Science Software 0