On Fri, 16 Oct 2015, John Peterson wrote: >> On Oct 16, 2015, at 11:39 AM, Stephan Herb <inf74...@stud.uni-stuttgart.de> >> wrote: > >> 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.
That's what it was designed for, IIRC, but I've likewise never used it. At least we've got an example with it: transient_ex2. Paul Bauman recently added a NewmarkSolver option to the FEMSystem framework, too; there's an example of that in fem_system_ex3. --- Roy ------------------------------------------------------------------------------ _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
