Matematické Fórum

Nevíte-li si rady s jakýmkoliv matematickým problémem, toto místo je pro vás jako dělané.

Nástěnka
22. 8. 2021 (L) Přecházíme zpět na doménu forum.matweb.cz!
04.11.2016 (Jel.) Čtete, prosím, před vložení dotazu, děkuji!
23.10.2013 (Jel.) Zkuste před zadáním dotazu použít některý z online-nástrojů, konzultovat použití můžete v sekci CAS.

Nejste přihlášen(a). Přihlásit

#1 24. 01. 2015 18:26

pajam
Zelenáč
Příspěvky: 14
Reputace:   
 

Metoda konečných prvků

Hezký večer,
potřebuji pomocí MKP spočítat dif. rovnici $-u'' - u =- x^2    , x\in (0,1)$
s okrajovými podmínkami : $u(0) = 0 , u'(1)=1.$


Zasekla jsem se u slabé formulace, protože přenásobením funkcí $v\in \{H^{1}(0,1), v(0)=0\}$ a integrací, dostanu rovnici  $-\int_{0}^{1} u''v - \int_{0}^{1}uv= \int_{0}^{1}-x^2v$ a při úpravě prvního integrálu pomocí per partes zůstanu u tvaru $-v(1)+\int_{0}^{1} u'v' $.

Asi by to mělo být bez té zavislosti na hodnotě funkce v.. pouze integrál... Můžu tedy do definice funkce $v$ přidat podmínku, aby byla i v pravém krajním bodě nulová, když tuto podmínku z okrajových podmínek nevím???

Offline

  • (téma jako vyřešené označil(a) pajam)

#2 25. 01. 2015 10:21

jardofpr
Příspěvky: 1241
Reputace:   88 
 

Re: Metoda konečných prvků

Zdravím,

povedal by som že to $v(1)$ pôjde k funkcionálu na pravej strane,

po diskretizácii a vyrobení bázy priestoru v ktorom sa hľadá riešenie by mali byť k dispozícii
konkrétne bázové funkcie $v_1,v_2,\dots,v_N$, z ktorých jediná, ktorá nebude mať
nulovú hodnotu v bode $x=1$ bude práve tá, ktorá je s ním asociovaná, t.j. hodnota $v(1)$
nakoniec ovplyvní len príslušnú jednu hodnotu vektora na pravej strane

Offline

 

Zápatí

Powered by PunBB
© Copyright 2002–2005 Rickard Andersson