Something More for Research

Explorer of Research #HEMBAD

Pendulum in viscous media

Posted by Hemprasad Y. Badgujar on October 12, 2014

Pendulum in viscous media

Category. Microcontroller firmware development
Note. The sample is a part of a production project

Download MatLab pendulum demo (zip, 2 Kb)


The article was written as result of R&D for a project related to development of ultrasound probe positioning mechanism. The problem considered in this article can be outlined as follows: sonic transmitter moves along arc in viscous media with some predefined velocity and the task is to switch the motion direction in minimum time.

The demo package includes 2 files — MatLab script and description:

  • pendulumdemo.m — MatLab script,
  • readme.txt — description.

Some technical notes on demo:

  • To run the sample MatLab should be installed on your computer.
  • Unpack the sample and copy pendulumdemo.m file into MatLab work directory.
  • Start upMatLab and type in:

    >> pendulumdemo

Problem consideration

Mathematically the problem can be posed as follows.

The equation of transmitter motion can be written down in the form

Equation (1) (1)

where Ω is angular displacement and U is external torque. The equation is equivalent to the equation of pendulum in viscous liquid. External torque is limited by some value:

Condition (2) (2)

Our boundary conditions are:

Boundary conditions (3) (3)

Here time T is unknown and minimum. These boundary conditions say we should return in the same point with opposite velocity. So our task is to find function U(t) which returns the transmitter in same point with opposite velocity in minimum time T. In other words we have problem of optimal control.

To solve the problem let us write down equation (1) in phase variables:

Phase variables

that gives

Equation (4) (4)

Now we can write Bellman equation:

Equation (5) (5)

from which taking into account restriction (2) we can conclude that:

Conditions (6) (6)

From system (4) we obtain

Equation (7) (7)

This equation has solution:

Solution (8) (8)

where C is an arbitrary constant. Now we accept conditions (8) and depending on sign of U we have two families of phase curves — fig. 1.

Phase planeFig. 1. Phase plane.

As one can see to solve the problem it is necessary to apply maximum negative driving torque and then switch it to maximum positive. Now the only point we should make clear is time instance tsw of switching driving moment from –Umax to Umax. To find time instance tsw let us write down solution of equation (1) for U that does not depend on time:

Solution (9) (9)

where C0 and C1 are arbitrary constants. The angular velocity respectively is

Solution (10) (10)

Accepting boundary conditions (3) we find time instance tsw:

Solution (11) (11)

Now knowing switch instance tsw we can build solution of the problem in phase plane — fig. 2. Here switching process is highlighted in green.

Velocity switchingFig. 2. Velocity switching.

Dependence of angular displacement, velocity and acceleration on time is depicted in figures 3—5 respectively.

Angular displacementFig. 3. Angular displacement.Angular velocityFig. 4. Angular velocity.Angular accelerationFig. 6. Angular acceleration.


Leave a Reply

Fill in your details below or click an icon to log in: Logo

You are commenting using your account. Log Out /  Change )

Google+ photo

You are commenting using your Google+ account. Log Out /  Change )

Twitter picture

You are commenting using your Twitter account. Log Out /  Change )

Facebook photo

You are commenting using your Facebook account. Log Out /  Change )

Connecting to %s

Extracts from a Personal Diary

dedicated to the life of a silent girl who eventually learnt to open up

Num3ri v 2.0

I miei numeri - seconda versione


Just another site

Algunos Intereses de Abraham Zamudio Chauca

Matematica, Linux , Programacion Serial , Programacion Paralela (CPU - GPU) , Cluster de Computadores , Software Cientifico




A great site

Travel tips

Travel tips

Experience the real life.....!!!

Shurwaat achi honi chahiye ...

Ronzii's Blog

Just your average geek's blog

Karan Jitendra Thakkar

Everything I think. Everything I do. Right here.


News About Tech, Money and Innovation

Chetan Solanki

Helpful to u, if u need it.....


Explorer of Research #HEMBAD


Explorer of Research #HEMBAD


A great site


This is My Space so Dont Mess With IT !!

%d bloggers like this: