Am 16.10.2015 um 20:17 schrieb Roy Stogner: > > 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 Hey, super! I think that's exactly what I need. Thanks for your quick help. It's absolutely amazing what libMesh has on offer.
Bye, Stephan ------------------------------------------------------------------------------ _______________________________________________ Libmesh-users mailing list Libmesh-users@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-users
