Code_Aster ®
Version
8.2
Titre :
Opérateur RESO_LDLT
Date :
31/01/06
Auteur(s) :
J. PELLET Clé
:
U4.55.02-I1 Page
: 1/4
Organisme(s) : EDF-R&D/AMA
Manuel d'Utilisation
Fascicule U4.5- : Méthodes de résolution
Document : U4.55.02
Opérateur RESO_LDLT
1 But
Résoudre un système d'équations linéaires par une méthode "directe".
L'opérateur permet des résolutions réelles ou complexes. La matrice doit avoir été préalablement
factorisée par la commande FACT_LDLT [U4.55.01].
Malgré son nom, cette commande correspond aux trois méthodes de résolution "directe" implantées
dans Aster : la résolution LDLT ,la méthode multifrontale et la méthode MUMPS. Le choix effectif de la
méthode se fait au travers de la commande NUME_DDL [U4.61.11].
La résolution par gradient conjugué se fait par la commande RESO_GRAD [U4.55.04].
Produit une structure de données de type cham_no_*.
Manuel d'Utilisation
Fascicule U4.5- : Méthodes de résolution
HT-62/06/004/A
Code_Aster ®
Version
8.2
Titre :
Opérateur RESO_LDLT
Date :
31/01/06
Auteur(s) :
J. PELLET Clé
:
U4.55.02-I1 Page
: 2/4
2 Syntaxe
U
[cham_no_*]
=
RESO_LDLT
( reuse = U,
MATR_FACT
=
A,
/
[matr_asse_DEPL_R]
/
[matr_asse_DEPL_C]
/
[matr_asse_TEMP_R]
/
[matr_asse_TEMP_C]
/
[matr_asse_PRES_R]
/
[matr_asse_PRES_C]
CHAM_NO
=
B,
/
[cham_no_DEPL_R]
/
[cham_no_DEPL_C]
/
[cham_no_TEMP_R]
/
[cham_no_TEMP_C]
/
[cham_no_PRES_R]
/
[cham_no_PRES_C]
CHAM_CINE
= vcine , / [cham_no_DEPL_R]
/
[cham_no_TEMP_R]
/
[cham_no_PRES_C]
# si solveur MUMPS :
RESI_RELA
= / 1.e-6 , [DEFAUT]
/
eps
, [R]
TITRE
= titr ,
[l_K80]
INFO
=
/ 1 , [DEFAUT]
/
2
,
)
si CHAM_NO :
[cham_no_DEPL_R] alors [*]
->
DEPL_R
[cham_no_DEPL_C]
DEPL_C
[cham_no_TEMP_R]
TEMP_R
[cham_no_TEMP_C]
TEMP_C
[cham_no_PRES_R]
PRES_R
[cham_no_PRES_C]
PRES_C
Manuel d'Utilisation
Fascicule U4.5- : Méthodes de résolution
HT-62/06/004/A
Code_Aster ®
Version
8.2
Titre :
Opérateur RESO_LDLT
Date :
31/01/06
Auteur(s) :
J. PELLET Clé
:
U4.55.02-I1 Page
: 3/4
3 Généralités
Cette commande permet de résoudre les systèmes linéaires AX = B où A est une matrice
préalablement "factorisée" par la commande FACT_LDLT [U4.51.01].
La résolution est possible pour des conditions aux limites de DIRICHLET (conditions aux limites
cinématiques) dualisées ou éliminées [U2.01.02]. Dans ce dernier cas, si le chargement X = X0
sur le "bord" 0 a été traduit par une charge cinématique (opérateur AFFE_CHAR_CINE [U4.44.03]
prise en compte dans la matrice assemblée (opérateur ASSE_MATRICE [U4.61.22], la "valeur" de ce
chargement (X0) , calculée par l'opérateur CALC_CHAR_CINE [U4.61.03] doit être fournie par le mot
clé CHAM_CINE.
4 Opérandes
4.1 Opérande
MATR_FACT
MATR_FACT = A ,
Nom de la matrice factorisée, concept produit par l'opérateur FACT_LDLT. Cette matrice peut être
réelle ou complexe, symétrique ou non.
4.2 Opérande
CHAM_NO
CHAM_NO = B
,
Nom du vecteur second membre (en général obtenu par la commande ASSE_VECTEUR).
4.3 Opérande
CHAM_CINE
CHAM_CINE = vcine ,
Nom du vecteur représentant la "valeur" des conditions aux limites de type "DIRICHLET" traduites
sous forme de chargement cinématique (c'est à dire par utilisation d'une des commandes
AFFE_CHAR_CINE ou AFFE_CHAR_CINE_F).
Ce cham_no provient de l'exécution de l'opérateur CALC_CHAR_CINE sur la liste des char_cine
(chargements cinématiques) associée à la matrice assemblée A [U2.01.02].
4.4 Opérande
RESI_RELA
RESI_RELA = 1.e-6 (défaut) / eps
Ce mot clé est décrit dans [U4.50.01]
4.5 Opérande
TITRE
TITRE =
titr ,
Titre que l'on veut donner au résultat produit [U4.03.01].
4.6 Opérande
INFO
INFO =
1 : pas d'impression.
Manuel d'Utilisation
Fascicule U4.5- : Méthodes de résolution
HT-62/06/004/A
Code_Aster ®
Version
8.2
Titre :
Opérateur RESO_LDLT
Date :
31/01/06
Auteur(s) :
J. PELLET Clé
:
U4.55.02-I1 Page
: 4/4
5 Exemples
·
Constitution des matrices assemblées :
On a calculé auparavant les termes élémentaires Kel , Fel.
nu = NUME_DDL ( MATR_RIGI = Kel )
Kass =
ASSE_MATRICE
(MATR_ELEM = Kel,
NUME_DDL = nu, )
Fass =
ASSE_VECTEUR
(MATR_ELEM = Fel,
NUME_DDL = nu, )
·
Factorisation en place :
Kass
= FACT_LDLT
(reuse = Kass,
MATR_ASSE = Kass,
)
·
Résolution hors-place :
Uass
= RESO_LDLT
(MATR_FACT = Kass, CHAM_NO = Fass, )
·
Résolution en place :
Fass =
RESO_LDLT
(reuse=Fass,
MATR_FACT=Kass,CHAM_NO=Fass,)
·
pour l'utilisation des charges cinématiques (avec élimination des degrés de liberté imposés), voir
l'exemple donné dans la commande AFFE_CHAR_CINE [U4.44.03].
Manuel d'Utilisation
Fascicule U4.5- : Méthodes de résolution
HT-62/06/004/A
Document Outline