In finite element, I can calculate the Lagrangian interpolation shape functions for each degree of freedom in an element, from the the number of nodal degrees of freedom and the number of nodes required in the element, if I know the data points ($x$,$y$ coordinates) required to generate linear algebraic conditions, which is in number of total degrees of freedom in the element.

In my problem, for a 2D element I have 6 degree of freedom in each node, namely $\phi, \phi _x, \phi _y, \phi _{xx}, \phi _{yy}$ and $\phi _{xy}$. As the name indicate, the five degrees of freedom are derivative of the $\phi$ degree of freedom.

If I want to calculate, shape functions for degree of freedom in a quadrilateral of four nodes, I would select a incomplete polynomial of degree five and $24$ terms. But the complexity will increase, if the number of nodes are more.

I would like to know that is there any programming library available which can be used for generation of these shape functions for higher nodes?

  • $\begingroup$ do you know how to compute the dual basis of $P_1(T)$ where $T =$ simplex of $R^d$, as in this problem: scicomp.stackexchange.com/questions/25719/… $\endgroup$
    – user177196
    Commented Dec 5, 2016 at 15:24
  • $\begingroup$ @user177196 Sorry, I'm very beginner in computational science $\endgroup$
    – user294664
    Commented Dec 5, 2016 at 15:51
  • $\begingroup$ For what kind of problem do you want to do this? $\endgroup$
    – P. G.
    Commented Dec 5, 2016 at 20:14
  • $\begingroup$ @P.G., trying to solve non-linear analysis problem by optimising implicit constitutive relation involving Airy's stress potential. These are the stress degree of freedom given in the problem, derived from Airy's stress potential $\endgroup$
    – user294664
    Commented Dec 6, 2016 at 5:57
  • $\begingroup$ @user294664 ok, not what I expected. But anyway, I looked up Airy's stress function. To me it looks a bit like the equations for plate bending. If you have four nodes, you can't have Lagrangian shape functions with 24 terms. You should look into hermitian shape fuctions for this problem. $\endgroup$
    – P. G.
    Commented Dec 6, 2016 at 11:45


Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge you have read our privacy policy.