DEBUT() MAYA=LIRE_MAILLAGE() MAYA=DEFI_GROUP( reuse=MAYA, MAILLAGE=MAYA, CREA_GROUP_MA=_F( NOM = 'TOUT', TOUT = 'OUI')) BARRE1=AFFE_MODELE( MAILLAGE=MAYA, AFFE=_F( GROUP_MA='SS1', MODELISATION = 'POU_D_E', PHENOMENE = 'MECANIQUE')) # MATERIO1=DEFI_MATERIAU( ELAS=_F( RHO = 1.E08, NU = 0.3, E = 1.E10)) MATERIO2=DEFI_MATERIAU( ELAS=_F( RHO = 0., NU = 0., E = 1.E15)) # CHMAT=AFFE_MATERIAU( MAILLAGE=MAYA,AFFE=( #_F( TOUT = 'OUI', MATER = MATERIO1), _F( GROUP_MA = 'SS1', MATER = MATERIO1), _F( GROUP_MA = 'SS2', MATER = MATERIO1), _F( GROUP_MA = 'SS3', MATER = MATERIO1), _F( GROUP_MA = 'RIGIDE', MATER = MATERIO2), )) # commentaire test CARA1=AFFE_CARA_ELEM( MODELE=BARRE1, POUTRE=_F( GROUP_MA = 'TOUT', SECTION = 'CERCLE', CARA = ('R', 'EP',), VALE = (0.1, 0.01,))) GUIDAGE1=AFFE_CHAR_MECA( MODELE=BARRE1,DDL_IMPO=( _F( GROUP_MA='SS1', DX = 0., DZ = 0., DRX = 0., DRY = 0.), _F( GROUP_NO = 'L1', DY = 0., DRZ = 0.), )) K_ELEM1=CALC_MATR_ELEM( MODELE=BARRE1, CARA_ELEM=CARA1, CHAM_MATER=CHMAT, OPTION='RIGI_MECA', CHARGE=GUIDAGE1) NUM1=NUME_DDL( MATR_RIGI=K_ELEM1) # K_ASSE1=ASSE_MATRICE( MATR_ELEM=K_ELEM1, NUME_DDL=NUM1) MODESTA1=MODE_STATIQUE( MATR_RIGI=K_ASSE1, MODE_STAT=_F( GROUP_NO=('L12'), AVEC_CMP = ('DY','DRZ'), )) # parse: -affectation DS1=[None]*5 DS2=[None]*5 DS3=[None]*5 DS4=[None]*5 CHS1=[None]*5 CHS2=[None]*5 # parse: +affectation #for k in range(1,5): # DS1[k] = CREA_CHAMP( OPERATION='EXTR', TYPE_CHAM='NOEU_DEPL_R', # RESULTAT= MODESTA1, NUME_ORDRE=k, NOM_CHAM = 'DEPL'); FIN()