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].