Exemple d'utilisation
=====================


.. code-block:: text

    # 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',),);