Scribd
Upload
201 views

Solving Delay Differential Equations (Maple)

The document discusses several examples of delay differential equations modeled in Maple. It introduces simple harmonic motion with lag, compares it to a growth model, and discusses implementation details like interpolation methods and handling of small delays. It then presents a predator-prey model with delay, an example showing chaotic behavior, and models with multiple delays and history-dependent solutions. Further examples include a suitcase rocking model with human response delay, showing the system becomes unstable with the introduction of delay.

Uploaded by

peanut pain
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
201 views

Solving Delay Differential Equations (Maple)

The document discusses several examples of delay differential equations modeled in Maple. It introduces simple harmonic motion with lag, compares it to a growth model, and discusses implementation details like interpolation methods and handling of small delays. It then presents a predator-prey model with delay, an example showing chaotic behavior, and models with multiple delays and history-dependent solutions. Further examples include a suitcase rocking model with human response delay, showing the system becomes unstable with the introduction of delay.

Uploaded by

peanut pain
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 12

Delay Differential Equations in Maple

Allan Wittkopf
Maplesoft

Delay Example
Modeling simple harmonic motion with lag:

>

(1.1)

>
(1.2)
>
Compare to growth:

>

(1.3)

>
(1.4)
>
Implementation Details
Mechanism: Natural C[1] interpolant: rkf45, rosenbrock, ck45
Initial Conditions: Use constant assumption
Delay storage: fixed ('delaypts' option)
Variable delay: supported (use 'delaymax' option)
Small delays: extrapolation
Derivatives of delay: 1st order only (C[1] interpolant)

Preadator-Prey Model
Adapted from Hutchinson model: delay accomodates differences in resource consumption between
young and adult members of population:

>
>

(3.1)

>

>
>
>
Wille and Baker Example
Demonstrates chaotic behavior for simple first order ODE with delay
>

(4.1)
>
Multiple Delays and History
Solve:
>

(5.1)

With initial condition for


Constant assumption: use piecewise
>

(5.2)
(5.2)

>

Suitcase Model
This is the model from the paper:

S. Suherman, R.H. Plaut, L.T. Watson, S. Thompson,


"Human delayed response time in correcting the side-to-side motion of a two wheeled suitcase.
J. Sound Vibration 207 (1997).

http://www.researchgate.
net/publication/243364611_EFFECT_OF_HUMAN_RESPONSE_TIME_ON_ROCKING_INSTAB
ILITY_OF_A_TWO-WHEELED_SUITCASE

and describes correction in the side-to-side motion of a 2-wheeled suitcase with a human delay in
the response time.
The dde is as follows:

>

(6.1)

Where:

- effective moment of inertia of suitcase rocking about either wheel


- product of weight and the eff. width of the suitcase between
wheels
- product of weight and height of suitcase
- coefficient of the restoring moment
- amplitude, freq, phase of excitation moment
In addition, when the angle passes through 0, there is a loss of energy when one of the wheels
impacts the ground, and this is described by a decrease in the velocity based on a coefficient of
restitution, which we choose to have the value ,
We choose the following parameter values and initial conditions:
>

(6.2)
>
(6.3)

where the delay has been left unspecified.


The energy loss of the wheel striking the ground is handled through the following event that states
that when passes through 0, the velocity is reduced by :
>

(6.4)

Now consider the behavior of the system if there is no delay in the response time:
>
so we can see that the angle varies between approximately -0.92 and 1.16.

If, however, we introduce a 0.1 sec. delay in the response time, the situation is quite different:
>
which shows that with the delay, the system is visibly unstable.

You might also like