In der Methode der Finiten Elemente verwenden wir Trial-Functions, die eindeutig den kinematischen Koordinaten an den Element-Knoten zugeordnet werden. Diese besonderen Trial-Functions heißen Hermite-Polynome.
So können wir als Formfunktion für ein Finites Element ein Polynom dritten Grades
wählen, die gesuchten generalisierten Koordinaten des Systems sind dann
- .
Äquivalent - aber viel praktischer - ist die Darstellung durch die neuen Koordinaten
Die resultierende Formfunktion ist wiederum ein Polynom dritten Grades:
Maxima Source-Code
Für kubische Polynome finden Sie hier die Herleitung.
/* Maxima */
C : [[a[0], a [1],a[2], a [3]],
[W[0],Phi[0],W[1],Phi[1]]];
ansatz: [w[1] = sum(a[i]*x^(i),i,0,3)];
bc : [subst([x=0], subst(ansatz,w[1]) )= W [0],
subst([x=0],diff(subst(ansatz,w[1]),x))=Phi[0],
subst([x=l], subst(ansatz,w[1]) )= W [1],
subst([x=l],diff(subst(ansatz,w[1]),x))=Phi[1]];
sol[1] : solve(bc,C[1]);
ansatz: append(ansatz, [w[2] = sum(coeff(expand(subst(sol[1],subst(ansatz,w[1]))),C[2][i])*C[2][i],i,1,4)]);
Die Funktionen
heißen Hermite-Polynome.