On Thu, Jan 24, 2013 at 4:01 PM, Robert <libm...@robertweidlich.de> wrote:
>
> As I still need quantities on the reference configuration (mapping
> jacobian), on the current configuration (shape functions,
> their derivations, qp positions, normals) and the mapping between
> undeformed and deformed configuration I currently deform each
> element manually and reinit the fe myself.
>
The way I would do this (and have done in other non-libMesh codes) is use
the deformation gradient to be able to compute in reference configuration.
For traditional nonlinear elasticity (i.e. no ALE), you should never need
to store the current configuration mesh. If you want to visualize it, then
sure, but your finite element formulation, and any postprocessing of
reference and current configuration quantities, can be computed in the
reference configuration using the deformation gradient.
HTH,
Paul
------------------------------------------------------------------------------
Master Visual Studio, SharePoint, SQL, ASP.NET, C# 2012, HTML5, CSS,
MVC, Windows 8 Apps, JavaScript and much more. Keep your skills current
with LearnDevNow - 3,200 step-by-step video tutorials by Microsoft
MVPs and experts. ON SALE this month only -- learn more at:
http://p.sf.net/sfu/learnnow-d2d
_______________________________________________
Libmesh-devel mailing list
Libmesh-devel@lists.sourceforge.net
https://lists.sourceforge.net/lists/listinfo/libmesh-devel