Improve Yourself

Just another WordPress.com weblog

  • Top Clicks

    • None

Analytical Dynamics: Virtual Work Principle

Posted by yunusbillah on December 7, 2008

For static case;

  1. Virtual displacement is an imaginary displacement, there is no actual displacement occur in the system (the system is still in equilibrium)
  2. If a system forces is in equilibrium, the work done by externally applied forces through virtual displacements are ZERO
  3. The virtual displacements do not violate the system constraints.

Example

Consider mechanical system bellow

conservative dynamic system

conservative dynamic system

Posted in Analytical Dynamics | Leave a Comment »

Analytical Dynamics: Introduction

Posted by yunusbillah on December 7, 2008

Analytical mechanics or analytical dynamic principle;

  1. Consider a system as a whole, not as individual components
  2. Excludes the reaction and constraint forces
  3. Formulates dynamic problems in two scalar function; kinetic and potential energy
  4. Formulates equation of motion in generalized coordinate and generalized forces which are not necessary physical coordinates and forces.

Definitions:

Degree of Freedom

Set of independent coordinates that required to completely describe the system’s configuration.

Generalized Coordinate

A set of minimum number of coordinates necessary and sufficent to uniquely define the system’s configuration.

Constraints

Any mathematical expression relates the cordinate and limits the system freedom.

Example

Describe the DoF number, generalized coordinates, and constraint functions of the system bellow!

conservative dynamic system

conservative dynamic system

The motion of each mass can be described mathematically;

- pulley (rotating) position: \theta

- disk (rolling):  x=\theta \cdot{r}_{p}

- block (translation): y=\theta \cdot 2{r}_{p}

Since the motion of each mass in the system is function of \theta, thus the system is only 1 DoF dynamic system. The mathematical expression of disk and block position are holonomic constraints. In order to simplified the equation of motion, x (translation coordinate) is chosen as the generalized coordinate.

Posted in Analytical Dynamics, Fundamentals of Vibration | Leave a Comment »

Solve The Eigenvalue Problem Using Matlab, It’s so Easy!

Posted by yunusbillah on December 4, 2008

When you have any eigenvalue problem bellow;

\left(\left[A\right]-\lambda\left[I\right]\right)=0

you can solve it easily using matlab program. This software is really powerfull!

Example:

Given a matrix \left[A\right]=\begin{bmatrix} 2 & 1 & -1\\ 1 & 4 & 3\\-1 & 3 & 4\end{bmatrix}

Calculate all its eigenvalues and the corresponding eigenvectors using matlab!

Step 1: Describe matrix A

Start your matlab and type matrix A as follow:

>> A=[2 1 -1;1 4 3;-1 3 4]    (enter)

A =

2     1    -1
1     4     3
-1     3     4

Step 2: Calculate the eigenvalues and the corresponding eigenvectors

>> [V D]=eig(A)    (enter)

V =

-0.5774    0.8165    0.0000
0.5774    0.4082    0.7071
-0.5774   -0.4082    0.7071

D =

0.0000         0         0
0    3.0000         0
0         0    7.0000

You will get new matrices D and V. Matrix D represent the eigenvalues of matrix A with matrix V as its corresponding eigenvectors (normalized eigenvectors).

You can practice by yourself, because it is so easy. Try to solve several problem using matlab!

Posted in Matlab, Software | Leave a Comment »

Solve The Eigenvalue Problem, It is Easy!

Posted by yunusbillah on December 4, 2008

There are so many problems can be modeled as the eigenvalue problems. You can search this on search engine and you will find million articles about this or you can find many books related to this topic.

When it said the eigenvalue problem?

Any equations that can be written in matrix form bellow it is said as the eigenvalue problem:

\left[A\right]\vec{x}=\lambda \vec{x}

where \left[A\right] is n x n constant matrix and \vec{x}=\left[{x}_{1},{x}_{2},...,{x}_{n} \right] is called as eigenvector, and \lambda is an any scalar that we called eigenvalue. The solution of eigenvector \vec{x} depend on eigenvalue (\lambda).

How to Solve It?

It is easy!. This problem has “Trivial Solution” \vec{x}=\vec{0}, but of course we are not interest in that! Than, how could we find the “Non-trivial Solution”?

If we multiplied the right hand side of the eigenvalue problem by identity matrix \left[I\right], we can modified the equation without change it value into the standard eigenvalue problem bellow

\left(\left[A\right]-\lambda\left[I\right]\right)\vec{x}=0

This standard form has a unique characteristic. If  \left|\left[A\right]-\lambda\left[I\right]\right| = 0 (determinant = 0) then it will have non-trivial solution.

Example 1:

\left[A\right] = \begin{bmatrix}1 & 1\\ 1 & 1\end{bmatrix}

\left|\begin{bmatrix}1 & 1\\ 1 & 1\end{bmatrix}-\lambda\begin{bmatrix}1 & 0\\ 0 & 1\end{bmatrix}\right|=0

\left|\begin{bmatrix}\left(1-\lambda\right) & 1\\ 1 & \left(1-\lambda\right)\end{bmatrix}\right|=0

\left(\lambda-2\right)\lambda=0

Solve for the \lambda, then we get {\lambda}_{1}=2 and {\lambda}_{2}=0

The eigenvector for the corresponding eigenvalue is

for {\lambda}_{1}=2, then \begin{bmatrix}-1 & 1\\ 1 & -1\end{bmatrix}\begin{bmatrix}{x}_{1} \\ {x}_{2} \end{bmatrix}=\begin{bmatrix}0 \\ 0 \end{bmatrix}

Use Gauss-Jordan elimination to reduce the equation in the matrix (become single line or triangle)

\begin{bmatrix}-1 & 1\\ 0 & 0\end{bmatrix}\begin{bmatrix}{x}_{1} \\ {x}_{2} \end{bmatrix}=\begin{bmatrix} 0 \\ 0 \end{bmatrix}

let {x}_{1}=\alpha, from the matrix we know that {x}_{1}={x}_{2}=\alpha

hence, the corresponding eigenvector is {\vec{x}}_{\lambda1}=\alpha\begin{bmatrix}1 \\ 1 \end{bmatrix}

Do the same for {\lambda}_{2}=0 then we get

{\vec{x}}_{\lambda2}=\beta\begin{bmatrix}1 \\ -1 \end{bmatrix}

That’s so easy isn’t it?! Try with bigger matrix size and different number you will enjoy while your ability improve! You can use mathematics software like matlab, maple, or mathematica to check your calculation.

You can also learn and practice to use those software in the articles bellow:

- Solve The Eigenvalue Problem Using Matlab (click!)

- Solve The Eigenvalue Problem Using Mathematica

- Solve The Eigenvalue Problem Using Maple

After you mastered this, then we can try to use the eigenvalue problem to solve so many practical problem in various field.

Posted in Academic, Basic Mathematics, Mathematics | Leave a Comment »

Hello world!

Posted by yunusbillah on December 4, 2008

Welcome to WordPress.com. This is your first post. Edit or delete it and start blogging!

Posted in Academic | 1 Comment »