Correct me if I am wrong but wouldn't the Newmark system assume that damping force is of the form : [F]=Factor*d[u]/dt, I would like to have [F] that is not of this form but a more generalized expression that I can create.
Regards Amal ________________________________________ From: John Peterson [[email protected]] Sent: Friday, June 20, 2014 10:09 AM To: Sahai, Amal (ARC-TN)[Universities Space Research Association (USRA)] Cc: [email protected] Subject: Re: [Libmesh-users] 2nd Order PDE in time On Fri, Jun 20, 2014 at 10:53 AM, Sahai, Amal (ARC-TN)[Universities Space Research Association (USRA)] <[email protected]> wrote: > Hi what would be the best way of solving a second order linear partial > differential equation in time of the type in LibMesh: > > [M] d^2u/dt^2 = [K] u + [F] > > I could obviously break it down into sets of equations - one solving for the > first derivative and other for u itself but I'd like to know if there is an > easier way to implement this. Also, how does one account for [M]? I think this is what the NewmarkSystem is good for? It allows you to specify mass, damping, and stiffness matrices... see transient_ex2 for an example. -- John ------------------------------------------------------------------------------ HPCC Systems Open Source Big Data Platform from LexisNexis Risk Solutions Find What Matters Most in Your Big Data with HPCC Systems Open Source. Fast. Scalable. Simple. Ideal for Dirty Data. Leverages Graph Analysis for Fast Processing & Easy Data Exploration http://p.sf.net/sfu/hpccsystems _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
