Main Page | Report this Page
 
   
Science Forum Index  »  Math - Numerical Analysis Forum  »  Looking for an implementation of the Newmark scheme for nume
Page 1 of 1    
Author Message
mr.d.poon@gmail.com
Posted: Wed Jan 24, 2007 8:18 am
Guest
Hi everyone

This is my first post ever to this group. I don't lurk here, but after
trawling through the archives, I think Im safe in saying that this
question has not been asked before.

Ive got a highly nonlinear system, and I understand that the Newmark
scheme the traditional choice for solving these. Searching Netlib and
the like for an existing implementation of this scheme does not yield
any results. Does anybody have an implementation in C or Fortran that
they are willing to share?

Alternatively, are there any suggestions for any other schemes that
will cope nicely with a highly nonlinear system and have an
implemenation out there for ready download (free or otherwise)?

Cheers

Daniel Poon
www.romaxtech.com
Evgenii Rudnyi
Posted: Thu Jan 25, 2007 5:15 pm
Guest
Quote:
Ive got a highly nonlinear system, and I understand that the Newmark
scheme the traditional choice for solving these. Searching Netlib and
the like for an existing implementation of this scheme does not yield
any results. Does anybody have an implementation in C or Fortran that
they are willing to share?

I have an implementation in Mathematica at

http://modelreduction.com/soft/Post4MOR/Post4MOR.m

Search for AnsysTransientSolution that ends with /;
SecondOrderSystemQ[sys]. As you will see, the function is not that
long, so it should not be too difficult to program it. You need just to
have at hand basic linear algebra subroutines.

Best wishes,

Evgenii Rudnyi

erudnyi at cadfem point de

http://modelreduction.com/
mr.d.poon@gmail.com
Posted: Fri Jan 26, 2007 8:31 am
Guest
Quote:
http://modelreduction.com/soft/Post4MOR/Post4MOR.m

Search for AnsysTransientSolution that ends with /;
SecondOrderSystemQ[sys]. As you will see, the function is not that
long, so it should not be too difficult to program it. You need just to
have at hand basic linear algebra subroutines.

Best wishes,

Evgenii Rudnyi

Thanks, I'll look into it.

btw have you considered making it available using the LGPL, rather than
GPL?

Daniel
 
Page 1 of 1       All times are GMT - 5 Hours
The time now is Thu Jan 08, 2009 6:27 pm