6. Example of use#
# CONSTRUCTION DES MATRICES DU MODELE MECANIQUE
ASSEMBLAGE (MODELE = MODNUM,
CARA_ELEM = CARE,
CHARGE = CHARGE_L,
NUME_DDL =CO ("NAKED"),
SOLVEUR =_F (METHODE = 'MUMPS'),
MATR_ASSE =( _F (MATRICE = CO ("K"), OPTION = 'RIGI_MECA'),
_F (MATRICE = CO ("M"), OPTION = 'MASS_MECA'),),
);
# CONSTRUCTION OF THE MATRICE REDUITE OF GUYAN SERVANT OF MATRICE_NORME
# CL POUR THE CONSTRUCTION OF THE MATRICE REDUITE OF GUYAN
CHAR_G = AFFE_CHAR_MECA (MODELE = MODNUM, DDL_IMPO =( _F (_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 =( 'X',),),),);
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 (for example) LA MATRICE NORME COMME UNE COMBINAISON DES
MATRICES REDUITES FROM MASSE AND 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 OF THE MATRICE OF PROJECTION ENTRE MODELE NUMERIQUE AND 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 = ('X',),),);
### RESOLUTION FROM PROBLEME FROM ERC
CHOIX DES FREQUENCES SUR LESQUELLES WE'RE GOING RESOUDRE THE 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 “,),);