> On Oct 16, 2015, at 11:39 AM, Stephan Herb <[email protected]> > wrote: > > 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.
I think you can use the Newmark System for this...but I'm not too familiar with that class. ------------------------------------------------------------------------------ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
