TP4 - Aides : 3. 2 - Les termes constants /

TP4 - Aides : Etape 3 : schéma en temps : 3. 2 - Les termes constants / $t$

3. 2 - Les termes constants / $t$

construire les termes indépendants du temps :

Ici encore, les matrices $M$ , $K$ , $A=M +\theta c^2 \Delta t^2 K$ , la source et les CL sont indépendantes du temps, donc

$\color{mycolor}{\triangleright }$ on construit les matrices et $Fc$ ,
$\color{mycolor}{\triangleright }$ on fait la décomposition LU de $A$ ,

avant la boucle en temps (seules les expressions de $A$ et $B$ diffèrent du TP3).

Initialiser les 1ère et 2ème colonnes du tableau $\color{mycolor}{U}$ :
la 1ère colonne = $\color{mycolor}{Uo}$ et la 2ième = $\color{mycolor}{Uo+\Delta t \, V}$
Sous scilab : U( : ,1:2) = [Uo , Uo+ $d t$ *V ]