Simple finite difference solver to solve the Euler Bernoulli beam equation. More...
#include <BeamSolver.h>
Simple finite difference solver to solve the Euler Bernoulli beam equation.
2 4 d w d w rho*h --- + E I --- = -p(x,t) 2 4 dt dx
where I is the moment of inertia divided by the beam width and p(x,t) the pressure loading. Typical boundary conditions are supported. The equation is transformed into a first-order system and temporally discretized with Crank-Nicholson. A straightforward LR decomposition is used for direct solution of the set of linear equations.