Hi, I'm currently working on a linear elastic problem using the finite element method; in short, all I do is assemble the system matrix 'K' (consisting of the single element stiffness matrices for flat shell elements), the right hand side 'F' and solve the system 'K * u = F'. I'm using a LinearImplicitSolver in libmesh to do this (I oriented myself on example 4 of the equations systems).
Now I have to expand the model to be time dependent and structural dynamic which alters the equation to this: 'M * a'' + C * a' + K * a = F(t)' with M being the mass matrix, C the damping matrix, K still the stiffness matrix, and a'' = d^2a/dt^2, a'=da/dt, Time t. Can you please give me a hint how I can solve this equation with libmesh (what solver I should use and if there is perhaps an example/documentation of how this solver works and what I need to provide). If there are too less details in order to help me, I can explain further; but since the structure of the equation is the problem, I think it should suffice to comment on it. Thanks. Stephan ------------------------------------------------------------------------------ _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
