Hi Guys

I am currently developing my own Finite element toolbox on Matlab.
The idea is to incorporate truss, beams and shell elements in the toolbox and eventually develop something that can interface with gmsh and/or develop solutions to standardized components via txt files developed from an excel based system or something similar.

The one aspect of the solvers consists of the transient (direct integration) solutions.

I have previously developed my own Runge Kutta for a specified problem, in addition I have also used ode45 to develop solutions to specified problems. Has anyone managed to develop an algorithm using ode45 that can be be generalized to any "n" degrees of freedom?

I am running into some problems to generalize the algorithm...

Has anyone worked with the Newmark method? It seems that most software uses this method to develop solutions as I think it is far easier to generalize as it avoids the state space representation that ode45 would require.

Kind regards
Mishal Mohanlal