Numerické modelování elektromagnetických polí se s rozvojem výpočetní techniky a neustále rostoucí výkonností počítačů stalo spolu s optimalizačními technikami nepostradatelnou složkou návrhu konstrukcí nových elektrotechnických a elektronickýchzařízení i zařízení z ostatních oblastí technické praxe. Numerické modelování je také bezesporu nedílnou součástí komplexních analýz chování časoprostorových polí, které jsou důležité pro posouzení nových požadavků na kvalitu zařízení jako je elektromagnetická kompatibilita. Složité problémy řešené v současné technické praxi nelze zvládnout ve většině případů jinými prostředky než pomocí vhodných numerických metod za použití výkonných počítačů.
Analogicky
pro parabolický trojúhelník určený body třeba použít aproximaci polynomem 6
konstantami
fyexdxycybxaN
e
+++++= 22)(
1
Obvykle pro všechny typy prvků používají vhodně normované souřadnice.
Čtveřice souřadnic čtyřúhelníkového prvku vyžaduje tvarovou funkci čtyřmi konstantami.7a). 4. 4. Nulová mimo prvek Postupem stejným
jako lineárního trojúhelníku odvodíme, konstanty jsou řešením soustavy rovnic
. Lze něm zavést lineární tvarovou funkci čtyřmi
konstantami
dczybxaN
e
j +++=
)(
Tvarová funkce N1
(e)
vrcholu prvku rovna tomto vrcholu nulová třech
zbývajících uzlech celé protější základně.
Aproximace potenciálu úloze
Základní prvky čtyřstěn, pětistěn šestistěn jsou Obr.Modelování elektromagnetických polí 41
Ze všech tvarových funkcí uzlu lze sestrojit analogicky úlohou aproximační funkci
tohoto uzlu jako součet všech tvarových funkcí uzlu Princip patrný Obr.2b).7b), dostaneme globální aproximaci potenciálu celé oblasti
∑
=
=
NU
j
jj yxN
1
a ),(φφ
Trojice souřadnic vrcholů trojúhelníku umožnila definovat tři konstanty určující rovinu. Nejjednodušší lineární
čtyřstěn určen čtyřmi vrcholy. 4.
Sestavíme-li takovouto aproximační funkci pro každý uzel sítě NU, kde počet
uzlů viz Obr. to
jehlan jednotkové výšky vrcholem uzlu hranami spojujícími uzly sousedící uzlem j.
Použitá aproximace tvar
),(),(),(),(
)(
44
)(
33
)(
22
)(
11
)(
yxNyxNyxNyxN
eeeee
φφφφφ +++=
kde např.
∑=
jP
e
jj )(
kde počet prvků společným uzlem j. bilineární funkce
dcybxaxyN
e
+++=
)(
1
Výrazy pro koeficienty lze odvodit postupem stejným jako pro trojúhelník. funkce N1
(e)
už není lineární jako předchozí, ale tzv
