Sources/Lexikon/Hermitesche-Polynome: Unterschied zwischen den Versionen

Aus numpedia
Zur Navigation springen Zur Suche springen
(Die Seite wurde neu angelegt: „nmn“)
 
Keine Bearbeitungszusammenfassung
 
Zeile 1: Zeile 1:
nmn
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
 
::<math>\displaystyle w(x) = \sum_{i=0}^3 a_{i}\cdot x^{i}</math>
 
wählen, die gesuchten generalisierten Koordinaten des Systems sind dann
 
::<math>a_0, a_1, a_2, a_3</math>.
 
Äquivalent - aber viel praktischer - ist die Darstellung durch die neuen Koordinaten
 
::<math>\begin{array}{ll}W_0 &=w(0)\\ \Phi_0 &=w'(0)\\W_1 &=w(\ell)\\ \Phi_1 &=w'(\ell) \end{array}</math>
 
Die resultierende Formfunktion ist wiederum ein Polynom dritten Grades:
 
::<math>\displaystyle  w(x)=  \underbrace{\left( 1-\frac{3\cdot {{x}^{2}}}{{{\ell}^{2}}}+\frac{2\cdot {{x}^{3}}}{{{\ell}^{3}}}\right)}_{\displaystyle =: \phi_1(\frac{x}{\ell})} \cdot {{W}_{0}}      + \underbrace{\left( x-\frac{2\cdot {{x}^{2}}}{\ell}+\frac{{{x}^{3}}}{{{\ell}^{2}}}\right)}_{\displaystyle =: \phi_2(\frac{x}{\ell})} \cdot {{\Phi}_{0}}      + \underbrace{\left( \frac{3\cdot {{x}^{2}}}{{{\ell}^{2}}}-\frac{2\cdot {{x}^{3}}}{{{\ell}^{3}}}\right)}_{\displaystyle =: \phi_3(\frac{x}{\ell})} \cdot  {{W}_{1}}      + \underbrace{\left( \frac{{{x}^{3}}}{{{\ell}^{2}}}-\frac{{{x}^{2}}}{\ell}\right)}_{\displaystyle =: \phi_4(\frac{x}{\ell})} \cdot {{\Phi}_{1}}</math>
 
{{MyCodeBlock|title=Maxima Source-Code
|text=Für kubische Polynome finden Sie hier die Herleitung.
|code=
/* 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)]);
<syntaxhighlight lang="lisp" line start=1>
 
</syntaxhighlight>
}}
Die Funktionen
 
::<math>\begin{array}{ccc}    \phi_1(\xi) &=& 1-3\,\xi^2+2\,\xi^3\\  \phi_2(\xi) &=& \ell \left(\xi-2\,\xi^2+\xi^3\right)\\  \phi_3(\xi) &=& 3\,\xi^2-2\,\xi^3\\  \phi_4(\xi) &=& \ell \left(\xi^3-\ell\,\xi^2\right)\\ \end{array}</math>
 
heißen Hermite-Polynome.

Aktuelle Version vom 21. April 2021, 06:41 Uhr

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.