On Thu, 1 Dec 2011, Yusuke Sakamoto wrote: > Now, I would like to compute the deformation gradient in the deformed > domain when I assemble the stiffness matrix/load vector. In order to do > so, I need the shape functions generated from the current mesh, and the > shape functions generated from the initial mesh. It seems that > EquationSystems class can only hold one mesh. Do I need to create two > EquationSystems objects? If so, how can they communicate with each other?
If you're using Lagrange finite elements, then the easiest way to do this is probably to use a single Mesh, kept in the undeformed state, then store your deformed mesh coordinates as variables, then each time you're doing computations on an element you reinit one FE on the undeformed element, deform the Elem temporarily, reinit another FE on the deformed element, then deform the Elem back when you're done. --- Roy ------------------------------------------------------------------------------ All the data continuously generated in your IT infrastructure contains a definitive record of customers, application performance, security threats, fraudulent activity, and more. Splunk takes this data and makes sense of it. IT sense. And common sense. http://p.sf.net/sfu/splunk-novd2d _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
