Sources/Lexikon/Hermitesche-Polynome

Aus numpedia
Version vom 21. April 2021, 06:41 Uhr von Mechaniker (Diskussion | Beiträge)
(Unterschied) ← Nächstältere Version | Aktuelle Version (Unterschied) | Nächstjüngere Version → (Unterschied)
Zur Navigation springen Zur Suche springen

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.