Mathematica 9 is now available
Services & Resources / Wolfram Forums
-----
 /
MathGroup Archive
2004
*January
*February
*March
*April
*May
*June
*July
*August
*September
*October
*November
*December
*Archive Index
*Ask about this page
*Print this page
*Give us feedback
*Sign up for the Wolfram Insider

MathGroup Archive 2004

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: NDSolve Repeated convergence test failure

  • To: mathgroup at smc.vnet.net
  • Subject: [mg47122] Re: NDSolve Repeated convergence test failure
  • From: Paul Abbott <paul at physics.uwa.edu.au>
  • Date: Fri, 26 Mar 2004 03:56:26 -0500 (EST)
  • Organization: The University of Western Australia
  • References: <c3uem8$9ub$1@smc.vnet.net>
  • Sender: owner-wri-mathgroup at wolfram.com

In article <c3uem8$9ub$1 at smc.vnet.net>,
 Steven Rey Ortiz <sro7516 at cs.tamu.edu> wrote:

> Howdy!  I at trying to use NDSolve to evaluate a set of three ODEs (two
> second order and one first order).  However, by construction the third
> function value will eventually eliminate the second order terms of the
> other two equations, and when it does this, Mathematica returns a
> "Repeated convergence test failure" and stops evaluating.  Is there a way
> to continue the integration process beyond this point?  The differential
> equations are still meaningful.  Following is a link to the notebook with
> the specific equations.  The last line shows the error.
> 
> http://students.cs.tamu.edu/sro7516/double-pendulum.nb

Having a look at your Notebook, one sees that theta2 stablilizes before 
t=2. That is its first and second derivatives are effectively zero (the 
first and second derivatives are less than 10^(-6) in absolute value by 
t = 1.92). 

From this time onwards, the equation for theta2 is satisfied and the 
equation for theta1 effectively reduces to an exactly solvable trivial 
first order equation:

  3 (Pi (Cos[t] + Sin[t]) + 8 theta1'[t]) == 0

Substituting the solution to this equation into the one for mtilde2 one 
finds that it is also satisfied. Hence you have the solution, in closed 
form, after the singularity.

> Problem Background: I am experimenting with an adaptive controller for a
> double pendulum system.  The notebook first defines the equations of
> motion for the system.  The controller assumes an initial value for the
> mass of the second link, but adjusts this value based on feedback from the
> system.  When it has the correct value for the mass, the control law will
> cancel the secend order terms from the equations of motion (and correctly
> track the desired trajectory).

The question, then, is how to get NDSolve to do this automatically. 
Perhaps Rob Knapp from WRI will be able to offer some help here.

Cheers,
Paul

-- 
Paul Abbott                                   Phone: +61 8 9380 2734
School of Physics, M013                         Fax: +61 8 9380 1014
The University of Western Australia      (CRICOS Provider No 00126G)         
35 Stirling Highway
Crawley WA 6009                      mailto:paul at physics.uwa.edu.au 
AUSTRALIA                            http://physics.uwa.edu.au/~paul


  • Prev by Date: Re: Yet another Mathematica docs whopper
  • Next by Date: Re: Yet another Mathematica docs whopper
  • Previous by thread: NDSolve Repeated convergence test failure
  • Next by thread: Re: NDSolve Repeated convergence test failure