6. Exemple d’utilisation#
# CONSTRUCTION DES MATRICES DU MODELE MECANIQUE
ASSEMBLAGE(MODELE=MODNUM,
CARA_ELEM=CARE,
CHARGE=CHARGE_L,
NUME_DDL=CO("NU"),
SOLVEUR=_F( METHODE = 'MUMPS'),
MATR_ASSE=(_F( MATRICE = CO("K"), OPTION = 'RIGI_MECA'),
_F( MATRICE = CO("M"), OPTION = 'MASS_MECA'),),
);
# CONSTRUCTION DE LA MATRICE REDUITE DE GUYAN SERVANT DE MATRICE_NORME
# CL POUR LA CONSTRUCTION DE LA MATRICE REDUITE DE GUYAN
CHAR_G=AFFE_CHAR_MECA(MODELE=MODNUM,DDL_IMPO=(_F( GROUP_NO = 'OBSPOINT', DX=0. ),),);
ASSEMBLAGE(MODELE=MODNUM,
CARA_ELEM=CARE,
CHARGE=(CHAR_G,CHARGE_L),
NUME_DDL=CO("NU_G"),
SOLVEUR=_F( METHODE = 'MUMPS'),
MATR_ASSE=(_F( MATRICE = CO("K1"), OPTION = 'RIGI_MECA'),
_F( MATRICE = CO("M1"), OPTION = 'MASS_MECA'),),
);
MOD_STA=MODE_STATIQUE(MATR_RIGI=K1,
MATR_MASS=M1,
MODE_STAT=(_F(GROUP_NO='OBSPOINT',
AVEC_CMP=('DX',),),),);
NUME_RED=NUME_DDL_GENE(BASE=MOD_STA,STOCKAGE='PLEIN');
MA_RE=PROJ_MATR_BASE(BASE=MOD_STA,NUME_DDL_GENE=NUME_RED,MATR_ASSE=M);
K_RE=PROJ_MATR_BASE(BASE=MOD_STA,NUME_DDL_GENE=NUME_RED,MATR_ASSE=K);
# ON CHOISIT (par exemple)LA MATRICE NORME COMME UNE COMBINAISON DES
# MATRICES REDUITES DE MASSE ET RAIDEUR
G_GUY=COMB_MATR_ASSE(COMB_R=(_F(MATR_ASSE = MA_RE,COEF_R = 1.),
_F(MATR_ASSE = K_RE, COEF_R = 1.),),);
## CALCUL DE LA MATRICE DE PROJECTION ENTRE MODELE NUMERIQUE ET EXPERIMENTAL ## (concept corresp_2_mailla)
## METHODE COLLOCATION OBLIGATOIRE POUR CALC_ERC_DYN
MATPROJ = PROJ_CHAMP(METHODE='COLLOCATION',
PROJECTION='NON',
MODELE_1=MODNUM,
MODELE_2=MODEXP);
### CONDITIONNEMENET DES MESURES (concept MODES) SUR LES DDL QUI NOUS
### INTERESSENT
OBS = OBSERVATION( RESULTAT = MODES,
MODELE_1 = MODNUM,
MODELE_2 = MODEXP,
PROJECTION = 'OUI',
TOUT_ORDRE = 'OUI',
MATR_RIGI = KASEXP,
MATR_MASS = MASEXP,
NOM_CHAM = 'DEPL',
FILTRE = _F( GROUP_MA = 'MESURE',
NOM_CHAM = ('DEPL'),
DDL_ACTIF = ('DX',),),);
### RESOLUTION DU PROBLEME DE L'ERC
# CHOIX DES FREQUENCES SUR LESQUELLES ON VA RESOUDRE L'ERC:
LFREQ=DEFI_LIST_REEL(VALE=[0.12181191980055407,1.25*0.22507907903927654],);
# CHOIX DES FREQUENCES SUR LESQUELLES
ERC=CALC_ERC_DYN( EVAL_FONCTIONNELLE=”OUI”,
MATR_PROJECTION=MATPROJ,
MESURE=OBS,
CHAMP_MESURE=”DEPL”,
MATR_MASS=M,
MATR_RIGI=K,
MATR_NORME=G_GUY,
LIST_FREQ=LFREQ,
GAMMA=0.5,
ALPHA=0.5,
INFO=1,
SOLVEUR=_F(METHODE=”MUMPS”,),);