d5.03.03 Introduire un nouveau degré de liberté et les conditions aux limites associées#
Résumé :
Que faut-il faire lorsqu’on ajoute un nouveau degré de liberté dans l’une des grandeurs DEPL_R, TEMP_R ou PRES_C?
En particulier, que doit-on faire pour que l’on puisse «contraindre» ce degré de liberté par des relations linéaires dualisées (AFFE_CHAR_MECA par exemple) ou éliminées (AFFE_CHAR_CINE) ?
Exemple réunissant les trois étapes pour introduire le degré de liberté, GRX (gauchissement des poutres).#
Catalogue des grandeurs Commons/physical_quantities.py#
Il suffit de rajouter la nouvelle composante dans le tuple:
list_cmp_depl=(
“DX”, “DY”, “DZ”, “DRX”, “DRY”, “DRZ”,
“GRX” , “GLIS”, “PRES”, “PRE[2]”, “TEMP”, “PHI”,
“DH”, “H1X”, “H1Y”, “H1Z”, “H1PRE1”, “H2X”,
Modification du catalogue d’élément : Elements/d_depl_r.py#
Il suffit de rajouter la nouvelle composante dans le tuple:
for cmp in (
“DH”, “DRX”, “DRY”, “DRZ”, “DX”, “DY”, “DZ”, “E1X”, “E1Y”, “E1Z”, “E2X”, “E2Y”, “E2Z”,
“E3X”, “E3Y”, “E3Z”, “E4X”, “E4Y”, “E4Z”, “GLIS”, “GONF”, “GRX” , “H1PRE1”, …
Modification des catalogues de commande (AFFE_CHAR_XXX/_C/_F et AFFE_CHAR_CINE/_F)#
Prenons l’exemple de la commande AFFE_CHAR_MECA. Si on veut permettre de bloquer le degré de liberté GRX , il faut rajouter ce nom ( GRX ) à différents endroits du catalogue, comme on peut le voir ci-dessous.
#& MODIF COMMANDE DATE 20/07/2010
…
AFFE_CHAR_MECA=OPER(nom= »AFFE_CHAR_MECA »,op= 7,sd_prod=char_meca,
fr= »Affectation de charges et conditions aux limites …
reentrant=”n”,
…
DDL_IMPO =FACT(statut=”f”,max=”**”,
…
regles=(AU_MOINS_UN(“TOUT”,”GROUP_MA”,”MAILLE”,”GROUP_NO”,”NOEUD”),
AU_MOINS_UN(“DX”,”DY”,”DZ”,”DRX”,”DRY”,”DRZ”,” GRX “,”PRES”,”PHI”,
“TEMP”,”PRE1”,”PRE2”,”UI2”,”UI3”,”VI2”,”VI3”,”WI2”,
…
DRX =SIMP(statut=”f”,typ=”R” ),
DRZ =SIMP(statut=”f”,typ=”R” ),
GRX =SIMP(statut=”f”,typ=”R” ),
…
FACE_IMPO =FACT(statut=”f”,max=”**”,
…
regles=(UN_PARMI(“GROUP_MA”,”MAILLE”,),
AU_MOINS_UN(“DX”,”DY”,”DZ”,”DRX”,”DRY”,”DRZ”,” GRX “,”PRES”,”PHI”,
“TEMP”,”PRE1”,”PRE2”,”DNOR”,”DTAN”),
…
DRX =SIMP(statut=”f”,typ=”R” ),
DRY =SIMP(statut=”f”,typ=”R” ),
DRZ =SIMP(statut=”f”,typ=”R” ),
DNOR =SIMP(statut=”f”,typ=”R” ),
DTAN =SIMP(statut=”f”,typ=”R” ),
GRX =SIMP(statut=”f”,typ=”R” ),
…