FDEs
Moderator: Mike Everman
SAGEacious
minor rant, deleted  life's too short (and too long) to spent it working with anything but finest tools.

 Posts: 3716
 Joined: Tue Dec 07, 2004 6:51 pm
 Antipspambot question: 0
 Location: 41d 1' N 80d 22' W
what is a derivative? (re: FDEs)
This was at the bottom of the last page and I feel it was prematurely
advanced to the 'back' of this thread due to last post.
This is important to know because of what's to come next.
advanced to the 'back' of this thread due to last post.
This is important to know because of what's to come next.
Last edited by WebPilot on Tue Feb 12, 2008 8:30 pm, edited 1 time in total.

 Posts: 3716
 Joined: Tue Dec 07, 2004 6:51 pm
 Antipspambot question: 0
 Location: 41d 1' N 80d 22' W
central difference formulae (re: FDEs)
So, here we have formulas for f'(xj) and f''(xj) whose values at
a point xj can be determined by evaluating the function at
the points xj1, xj and xj+1.
The error decreases as one makes smaller the interval, h. The
'run time' will increase though.
Last edited by WebPilot on Tue Feb 12, 2008 8:32 pm, edited 1 time in total.

 Posts: 4930
 Joined: Fri Oct 31, 2003 7:25 am
 Antipspambot question: 0
 Location: santa barbara, CA
 Contact:
Is it just me, or are the picture links dead?
Mike
__________________________
Follow my technical science blog at: http://mikeeverman.com/
Get alerts for the above on twitter at: http://twitter.com/mikeeverman
__________________________
Follow my technical science blog at: http://mikeeverman.com/
Get alerts for the above on twitter at: http://twitter.com/mikeeverman

 Posts: 3716
 Joined: Tue Dec 07, 2004 6:51 pm
 Antipspambot question: 0
 Location: 41d 1' N 80d 22' W
central difference formulae (re: FDEs)
Appears to be, Mike.
Here is a repeat of the last two:
Here is a repeat of the last two:
When working with structural components calculating the deflection and such, may give an idea about how the structure will behave. Depending upon the Factor of Safety "allowed" may tell more as to why a structure appears very rigged for certain applications.
Sometimes you find yourself in weird areas, such as narrow load bearing columns where standard equations miss a buckling effect. Euler has some interesting formulas.
http://en.wikipedia.org/wiki/Buckling
Sometimes you find yourself in weird areas, such as narrow load bearing columns where standard equations miss a buckling effect. Euler has some interesting formulas.
http://en.wikipedia.org/wiki/Buckling

 Posts: 3716
 Joined: Tue Dec 07, 2004 6:51 pm
 Antipspambot question: 0
 Location: 41d 1' N 80d 22' W
discretization of the deflection (re: FDEs)
Pics are back online ...
~~~~~~~~~~~~~~~~~~~~~~~~~~~~
In those 'weird areas' where you cannot find a 'canned' formula
that applies, a numerical method is your only resort.
Euler was a prolific man; I see his name in numerical methods,
engineering mechanics and fluid mechanics books all the time.
~~~~~~~~~~~~~~~~~~~~~~~~~~~~
This is the beginning of the numerical approach, for which some of
you have patiently been waiting. I shall give it to you, the reader,
in 'bite size' chunks. If I have 'done a good job' in my preliminary
posts, I think some of you will 'pick it up'.
NOTE: M(x) means M is a function of x, not M times x ... I have
tried to be consistent and where two terms are multiplied together, I
have used a 'centered dot'.
NOTE #2: the notation, yj+1, means y(xj+1), etc.
NOTE #3: the f(x) term in the middle on the right hand side of the
central difference formula should have been written, f(xj) or more
compactly, fj.
NOTE #4: j and i are 'dummy variables' and can take on the values
0, 1, 2, 3, 4 ... n.
~~~~~~~~~~~~~~~~~~~~~~~~~~~~
Thanks PyroJoe,PyroJoe wrote:
When working with structural components calculating the deflection and such, may give an idea about how the structure will behave. Depending upon the Factor of Safety "allowed" may tell more as to why a structure appears very rigged for certain applications.
Sometimes you find yourself in weird areas, such as narrow load bearing columns where standard equations miss a buckling effect. Euler has some interesting formulas.
http://en.wikipedia.org/wiki/Buckling
In those 'weird areas' where you cannot find a 'canned' formula
that applies, a numerical method is your only resort.
Euler was a prolific man; I see his name in numerical methods,
engineering mechanics and fluid mechanics books all the time.
~~~~~~~~~~~~~~~~~~~~~~~~~~~~
This is the beginning of the numerical approach, for which some of
you have patiently been waiting. I shall give it to you, the reader,
in 'bite size' chunks. If I have 'done a good job' in my preliminary
posts, I think some of you will 'pick it up'.
NOTE: M(x) means M is a function of x, not M times x ... I have
tried to be consistent and where two terms are multiplied together, I
have used a 'centered dot'.
NOTE #2: the notation, yj+1, means y(xj+1), etc.
NOTE #3: the f(x) term in the middle on the right hand side of the
central difference formula should have been written, f(xj) or more
compactly, fj.
NOTE #4: j and i are 'dummy variables' and can take on the values
0, 1, 2, 3, 4 ... n.

 Posts: 3716
 Joined: Tue Dec 07, 2004 6:51 pm
 Antipspambot question: 0
 Location: 41d 1' N 80d 22' W
Matrix inversion in 2 simple steps (re: FDEs)
Excel has some powerful tools.
Let's invert the matrix, A.
Define Matrix, A
1. select cells C14:M24 and click Insert  Name  Define  A
Select where you wish the inverse to be located
2. select cells C27:M37 and type in window
 minverse(A)
 and then press CtrlShiftEnter (3 keys simultaneously)
(c.t.)
Voilà! there it is! Rather easy, n'estce pas ?
Let's invert the matrix, A.
Define Matrix, A
1. select cells C14:M24 and click Insert  Name  Define  A
Select where you wish the inverse to be located
2. select cells C27:M37 and type in window
 minverse(A)
 and then press CtrlShiftEnter (3 keys simultaneously)
(c.t.)
Voilà! there it is! Rather easy, n'estce pas ?

 Posts: 3716
 Joined: Tue Dec 07, 2004 6:51 pm
 Antipspambot question: 0
 Location: 41d 1' N 80d 22' W
compare 'fdm' and 'exact' plots {re:FDEs}
For plotting purposes, we need to compute xi/L (the station numbers).
Please see previous pic.
6. type xi/L in cell F39
7. type =B14/L in cell F40
8. select and copy this into cells F41:F50
 this is done by left mouse clicking cell F:40 and 'dragging'
its lower right corner down to cell F50
We are ready to plot the numerical solution against the
analytical one that I (we) did previously.
(c.t.)
The analytical (exact) maximum deflection is 2.396 and
the numerical (fdm) value is 2.415 . The fdm method differs
from the exact at this location by a mere +0.793% (< 1%).
If you wish to experiment, change the value of h.
Please see previous pic.
6. type xi/L in cell F39
7. type =B14/L in cell F40
8. select and copy this into cells F41:F50
 this is done by left mouse clicking cell F:40 and 'dragging'
its lower right corner down to cell F50
We are ready to plot the numerical solution against the
analytical one that I (we) did previously.
(c.t.)
The analytical (exact) maximum deflection is 2.396 and
the numerical (fdm) value is 2.415 . The fdm method differs
from the exact at this location by a mere +0.793% (< 1%).
If you wish to experiment, change the value of h.