5. Examples#

5.1. Degrees of freedom imposed in mechanics#

dog = AFFE_CHAR_CINE (MODELE = me,

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

For this plate problem in plane \(\mathrm{XY}\), we block all the degrees of freedom of rotation around \(Z\) and we embed the plate on its edge \(\mathrm{bord1}\).

5.2. Comparative use of kinematic and « ordinary » loads#

5.2.1. Global orders#

ch1 = AFFE_CHAR_THER (...)

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

evoth = THER_LINEAIRE (...

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

There is no difference.

5.2.2. « Step by step » calculation#

Ordinary expenses

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

Kinematic loads

honey = 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 = field, CHAM_NO = F, CHAM_CINE = vcine)

The terms induced by the kinematic loads are carried over to the second member, which requires the calculation of a field at the additional nodes vcine by the command CALC_CHAR_CINE [U4.61.03].