On Thu, Mar 19, 2009 at 11:07 AM, Vasilis Vavourakis <[email protected]> wrote: > i don't know if this might sound a silly question BUT i want to > perform the inversion of a "DenseMatrix" of a system... > > to be more specific, i want to implement the "arc-length" procedure > for non-linear structural problems. > for this topic see: Crisfield M.A. - Non-Linear Finite Element > Analysis of Solids and Structures Vol. 1 (see pages: 266-286) > one of the solutions proposed in the book is to invert the "tangent > stiffness matrix" (the Jacobian matrix in other words)
Yes, this is a standard technique but the Jacobian is not dense if you are using a standard finite element discretization. Usually it is "bordered" which means there is one dense row and column augmenting the original system. -- John ------------------------------------------------------------------------------ Apps built with the Adobe(R) Flex(R) framework and Flex Builder(TM) are powering Web 2.0 with engaging, cross-platform capabilities. Quickly and easily build your RIAs with Flex Builder, the Eclipse(TM)based development software that enables intelligent coding and step-through debugging. Download the free 60 day trial. http://p.sf.net/sfu/www-adobe-com _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
