From: Yusuke Sakamoto <ysakamoto@ic...>  20111202 02:32:29

Hi, I am working on the problem with large deformation, where the mesh is updated (deformed) at each time step quasistatically. It means that I solve for velocities at each time step and use that to compute displacements. Deforming the mesh is easy. I only had to map every node to new configuration according to its displacement. The connectivity of the nodes is unchanged. 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? Thanks, Yusuke 