Jon, > I'm working on an axi-symmetric formulation of elasticity in cylindrical > coordinates, and I'm a bit stuck when it comes to computing the gradients. > The problem is that I need to compute gradients of the shape functions that > contain values of shape functions as well. > > For example, for a vector-valued shape function v^i, with components > [v^i_r, v^i_\phi, v^i_z],
I'm a little confused. You say that you are using an axi-symmetric formulation, but then you typically only have the r- and z-components of the solution, because the solution is constant in phi direction, no? In that case, the actual 3d gradient in x-y-z-coordinates at a position (r,z) of the x-z-plane would be [\partial_r v^i, 0, \partial_z v^i] Or are you instead trying to work in cylindrical coordinates without the assumption of symmetry? W. ------------------------------------------------------------------------- Wolfgang Bangerth email: [email protected] www: http://www.math.tamu.edu/~bangerth/ _______________________________________________ dealii mailing list http://poisson.dealii.org/mailman/listinfo/dealii
