5. Exemples

5.1. Degrés de liberté imposés en mécanique

chcine =    AFFE_CHAR_CINE    (MODELE = mo,

MECA_IMPO=( _F (TOUT = “OUI” , DRZ = 0.), _F (GROUP_NO = “bord1”, DX = 0., DY = 0., DZ = 0., DRX = 0., DRY = 0.,)))

Pour ce problème de plaque dans le plan \(\mathrm{XY}\) , on bloque tous les degrés de liberté de rotation autour de \(Z\) et on encastre la plaque sur son bord \(\mathrm{bord1}\) .

5.2. Utilisation comparée des charges cinématiques et « ordinaires »

5.2.1. Commandes globales

ch1        =    AFFE_CHAR_THER    ( ... )

ch2        =    AFFE_CHAR_CINE_F    ( TEMP_IMPO =    _F ( ... ))

evoth    =    THER_LINEAIRE        (    ...

EXCIT = ( _F (CHARGE = ch1), _F (CHARGE = ch2),) …)

Il n’y a pas de différence.

5.2.2. Calcul « pas à pas »

Charges ordinaires

ch1 = AFFE_CHAR_MECA ( … ) mel = CALC_MATR_ELEM ( … OPTION = “RIGI_MECA” , CHARGE = ch1) matas = ASSE_MATRICE ( MATR_ELEM = mel …) matas = FACTORISER ( reuse = matas, MATR_ASSE = matas ) U = RESOUDRE ( MATR = matas , CHAM_NO = F )

Charges cinématiques

ch1 = AFFE_CHAR_CINE ( … ) mel = CALC_MATR_ELEM ( … OPTION = “RIGI_MECA” ) matas = ASSE_MATRICE ( MATR_ELEM = mel, …, CHAR_CINE = ch1 ) matas = FACTORISER ( reuse = matas, MATR_ASSE = matas, ) vcine = CALC_CHAR_CINE ( …, CHAR_CINE = ch2, ) U = RESOUDRE ( MATR = matas , CHAM_NO = F, CHAM_CINE = vcine)

Les termes induits par les charges cinématiques sont reportés au second membre ce qui nécessite le calcul d’un champ aux nœuds supplémentaire vcinepar la commande CALC_CHAR_CINE[U4.61.03].