u4.55.03 Opérateur ELIM_LAGR#

Syntaxe#

Détail de la syntaxe
matr_asse_elim_r = ELIM_LAGR(
    ◆ MATR_RIGI = matr_asse_depl_r,
    ◇ MATR_ASSE = matr_asse_depl_r,
    ◇ INFO = / 1,
             / 2,
)


◆ : obligatoire
◇ : optionnel
⟐ : présent par défaut
& : ensemble
/ : un seul parmi
| : plusieurs choix possibles

Opérandes#

Opérande MATR_RIGI#

♦    MATR_RIGI    =    K1,

Nom de la matrice assemblée de rigidité (avec des relations linéaires à éliminer)

Si le mot clé MATR_ASSE n’est pas utilisé, c’est la matrice K1 que l’on «réduit» pour créer la matrice résultat (M2).

Opérande MATR_ASSE#

◊    MATR_ASSE    =    M1,

Nom de la matrice assemblée masse, amortissement, … (pas de rigidité) que l’on veut réduire.

Si ce mot clé est utilisé, c’est la matrice M1 que l’on «réduit» pour créer la matrice résultat (M2).

Opérande TITRE#

◊    TITRE    =    titr ,

Titre que l’on veut donner au résultat produit [U4.03.01].

Opérande INFO#

◊    INFO    =

1: pas d’impression.

2: impressions

Calcul sur plusieurs processeurs#

Lorsqu’on utilise plusieurs processeurs, lesdonnées du modèle doivent être répliquées sur tous les processeurs participant au calcul en utilisant le mot-clé DISTRIBUTION de la commande AFFE_MODELE

MODEL=AFFE_MODELE(…

DISTRIBUTION=_F(METHODE='CENTRALISE'))

Le calcul peut alors s’effectuer sur plusieurs processeurs. Attention, tous les processeurs effectuant les mêmes opérations, il n’y a pas de gain de calcul pour l’étape d’élimination des multiplicateurs de Lagrange. L’intérêt est d’effectuer ensuite un calcul modal sur plusieurs processeurs.

Exemple#

Calcul modal sur des matrices avec ou sans équations de Lagrange#

K1 = ASSE_MATRICE(NUME_DDL=NU, MATR_ELEM=KEL,)

M1 = ASSE_MATRICE(NUME_DDL=NU, MATR_ELEM=MEL,)


# 1. calcul avec les matrices complètes :

#-------------------------------------------------------


mode1 = CALC_MODES( MATR_RIGI=K1,

MATR_MASS=M1,

OPTION='BANDE',

CALC_FREQ=_F( FREQ=(-2, 30),) )


# 2. calcul avec les matrices réduites par ELIM_LAGR :

#-------------------------------------------------------


K2=ELIM_LAGR(MATR_RIGI=K1, )

M2=ELIM_LAGR(MATR_RIGI=K1, MATR_ASSE=M1)


mode2 = CALC_MODES( MATR_RIGI=K2,

MATR_MASS=M2,

OPTION='BANDE',

CALC_FREQ=_F( FREQ=(-2, 30),) )

mode3 = REST_GENE_PHYS(RESU_GENE=mode2)