Télécharger g_cas.procedur

Retour à la liste

* G_CAS     PROCEDUR  CB215821  21/06/10    21:15:25     11029          
* =============================================================================
*       PROCEDURE D'ANALYSE DE LA TABLE D'ENTREE FOURNIE A G_THETA
*       ----------------------------------------------------------
*
* DESCRIPTION : DETERMINE SI LE CALCUL DEMANDE EST PREVU PAR LA VERSION
*               ACTUELLE DE G_THETA, ET PREVIENT L'UTILISATEUR SI LE CAS
*               ENVISAGE N'EST PAS TESTE OU VALIDE PAR LA BASE DE CAS TESTS
*               ACTUELLE.
*               ON EN PROFITE EGALEMENT POUR CREER CERTAINS OBJETS QUI SERONT
*               UTILISES PAR G_THETA.
* =============================================================================
DEBP G_CAS SUPTAB*'TABLE' BOOL*'TABLE' ;
 
* I - DETERMINATION DU CAS DEMANDE ET VERIFICATION DES DONNEES OBLIGATOIRES
* -------------------------------------------------------------------------
 
SAUT 'LIGNE' ;
MESS 'IDENTIFICATION DU CAS DEMANDE PAR L''UTILISATEUR :' ;
SAUT 'LIGNE' ;
NDEC = 20 ;
 
* I.1 - LECTURE DES DONNEES FOURNIES A SUPTAB
* *******************************************
 
* I.1.A - OBJECTIF
SI (NON (EXIS SUPTAB 'OBJECTIF')) ;
    MESS 'ERREUR : IL FAUT SPECIFIER L INTEGRALE' ;
    MESS '      A CALCULER PAR UN MOT' ;
    ERRE 21 ; QUIT G_CAS ;
SINON ;
    SI (NEG (TYPE SUPTAB.'OBJECTIF') 'MOT') ;
        MESS 'ERREUR : OBJECTIF INCORRECT. ON ATTEND UN ''MOT''.' ;
        ERRE 21 ; QUIT G_CAS ;
    FINSI ;
    GOBJ = MOT SUPTAB.'OBJECTIF' ;
FINSI ;
* MESSAGE
MESSOBJ = CHAI 'OBJECTIF :' (MOT GOBJ)/NDEC ;
MESS MESSOBJ ;
 
* I.1.B - DIMENSION ET MODE DE CALCUL
GDIME = VALE 'DIME' ;
GMODE = VALE 'MODE' ;
* MESSAGE
MESSDIME = CHAI 'DIMENSION :' GDIME/NDEC ;
MESS MESSDIME ;
MESSMODE = CHAI 'MODE :' GMODE/NDEC ;
MESS MESSMODE ;
 
* I.1.C - TYPE DE SOLUTION DONNEE
SI (NON ((EXIS SUPTAB 'SOLUTION_RESO') OU (EXIS SUPTAB 'SOLUTION_PASAPAS'))) ;
    MESS 'ERREUR : IL FAUT UNE SOLUTION PROVENANT DE PASAPAS' ;
    MESS '       OU DE RESO POUR DETERMINER L''INTEGRALE' ;
    ERRE 21 ; QUIT G_CAS ;
SINON ;
    SI (EXIS SUPTAB 'SOLUTION_RESO') ;
        SI (NEG (TYPE SUPTAB.'SOLUTION_RESO') 'CHPOINT') ;
            MESS 'ERREUR : SOLUTION_RESO INCORRECTE.'
            MESS '       ON ATTEND UN ''CHPOINT''.' ;
            ERRE 21 ; QUIT G_CAS ;
        FINSI ;
        GSOLU = 'LINEAIRE' ;
    SINON ;
        GSOLU = MOT 'PASAPAS' ;
    FINSI ;
FINSI ;
 
* I.1.D - MODELE MECANIQUE, MATERIAU ET HYPOTHESE DE CALCUL
SI (EGA GSOLU 'LINEAIRE') ;
    SI (NON (EXIS SUPTAB 'MODELE')) ;
        MESS 'ERREUR : UNE SOLUTION ISSUE DE RESO NECESSITE UN MODELE' ;
        ERRE 21 ; QUIT G_CAS ;
    FINSI ;
    SI (NON (EXIS SUPTAB 'CARACTERISTIQUES')) ;
        MESS 'ERREUR : IL FAUT DONNER LE CHAMP CARACTERISTIQUE' ;
        ERRE 21 ; QUIT G_CAS ;
    FINSI ;
    MOD_TOT = SUPTAB.'MODELE' ;
    MOD_MEC = EXTR MOD_TOT 'FORM' 'MECANIQUE' ;
    MAT_MEC = REDU (SUPTAB.'CARACTERISTIQUES') MOD_MEC ;
    GHYPO = 'HPP' ;
SINON ;
    SI (EXIS SUPTAB.'SOLUTION_PASAPAS' 'WTABLE') ;
        MOD_TOT = SUPTAB.'SOLUTION_PASAPAS'.'MODELE' ;
        WTAB = SUPTAB.'SOLUTION_PASAPAS'.'WTABLE' ;
        MOD_MEC = WTAB.'MOD_MEC' ;
        MAT_MEC = WTAB.'MAT_MEC' ;
    SINON ;
        MESS 'Absence de WTABLE !  l execution continue ...' ;
        MOD_TOT = SUPTAB.'SOLUTION_PASAPAS'.'MODELE' ;
        WTAB = SUPTAB.'SOLUTION_PASAPAS' ;
        MOD_MEC = EXTR MOD_TOT 'FORM' 'MECANIQUE' ;
        MAT_MEC = REDU (WTAB.'CARACTERISTIQUES') MOD_MEC ;
    FINSI ;
    SUPTAB.'MODELE' = MOD_MEC ;
    SI WTAB.'GRANDS_DEPLACEMENTS' ;
        GHYPO = 'GRANDS_DEPLACEMENTS' ;
    SINON ;
        GHYPO = 'HPP' ;
    FINSI ;
    SI WTAB.'DYNAMIQUE' ;
        GSOLU = CHAI GSOLU ' ' 'DYNAMIQUE' ;
    SINON ;
        GSOLU = CHAI GSOLU ' ' 'STATIQUE' ;
    FINSI ;
FINSI ;
GMTYPE = TABL ;
* COMPORTEMENT ELASTIQUE
SI (EXIS MOD_MEC 'MATE' 'ISOTROPE') ;
    GMTYPE.'ELASTIQUE' = MOT 'ELASTIQUE ISOTROPE' ;
SINON ;
    GMTYPE.'ELASTIQUE' = MOT 'ELASTIQUE ANISOTROPE' ;
FINSI ;
* COMPORTEMENT PLASTIQUE
MPLAS = TABL ;
MPLAS.(1) MPLAS.(2) MPLAS.(3) = 'ISOTROPE' 'PARFAIT' 'CINEMATIQUE' ;
SI (EXIS MOD_MEC 'MATE' 'PLASTIQUE') ;
    GMTYPE.'PLASTIQUE' = MOT 'PLASTIQUE' ;
    REPE I (DIME MPLAS) ;
        SI (EXIS MOD_MEC 'MATE' MPLAS.&I) ;
            GMTYPE.'PLASTIQUE' = CHAI GMTYPE.'PLASTIQUE' ' ' MPLAS.&I ;
        FINSI ;
    FIN I ;
FINSI ;
* COMPORTEMENT VISCOPLASTIQUE
MVISC = TABL ;
MVISC.(1) MVISC.(2) MVISC.(3) MVISC.(4) MVISC.(5) MVISC.(6) =
    'NORTON' 'POLYNOMIAL' 'BLACKBURN' 'RCCMR_316' 'RCCMR_304' 'LEMAITRE' ;
SI (EXIS MOD_MEC 'MATE' 'FLUAGE') ;
    GMTYPE.'FLUAGE' = MOT 'FLUAGE' ;
    REPE I (DIME MVISC) ;
        SI (EXIS MOD_MEC 'MATE' MVISC.&I) ;
            GMTYPE.'FLUAGE' = CHAI GMTYPE.'FLUAGE' ' ' MVISC.&I ;
        FINSI ;
    FIN I ;
FINSI ;
 
* PROPRIETES QUI DEPENDENT D'UN PARAMETRE
CPARA = EXTR MAT_MEC 'COVA' ;
PPARA = MOTS ;
REPE I (DIME CPARA) ;
    PPARA = PPARA ET (EXTR CPARA &I) ;
FIN I ;
* PROPRIETES QUI DEPENDENT DE LA POSITION
MCOMP = EXTR MAT_MEC 'COMP' ;
PPOSI = MOTS ;
REPE I (DIME MCOMP) ;
    ICOMP = EXTR MCOMP &I ;
    SI (NON (EXIS PPARA ICOMP)) ;
        ICHAM = EXCO MAT_MEC ICOMP ;
 
*       On teste le Type car 'COVA' a des exceptions (comme 'ECRO' qui est un EVOLUTION mais consideree constante)
*       Amelioration a prevoir : boucle sur les SOUS-ZONES et s'assurer qu'il y a bien au moins 1 élément...
*       On aimerait bien avoir un operateur qui repond par VRAI ou FAUX si un objet est constant !!!
        ValTyp='TYPE' ('EXTR' ICHAM ICOMP 1 1 1);
        SI ('NEG' ValTyp ('MOT' 'FLOTTANT')); 'ITER' I ; 'FINS';
        SI (NEG (MINI ICHAM) (MAXI ICHAM)) ;
            PPOSI = PPOSI ET ICOMP ;
        FINSI ;
    FINSI ;
FIN I ;
* PROPRIETES A GRADIENT
PGRAD = PPARA ET PPOSI ;
* DISCONTINUITES MATERIELLES
SI (EXIS SUPTAB 'MODELES_COMPOSITES') ;
    DISMAT = VRAI ;
    TABMOD = SUPTAB.'MODELES_COMPOSITES' ;
SINON ;
    DISMAT = FAUX ;
FINSI ;
* MESSAGE TYPE DE SOLUTION
MESSSOLU = CHAI 'CALCUL :' GSOLU/NDEC ;
MESS MESSSOLU ;
* MESSAGE HYPOTHESE DE CALCUL
MESSHYPO = CHAI 'HYPOTHESE :' GHYPO/NDEC ;
MESS MESSHYPO ;
* MESSAGE MATERIAU
MESSMAT = CHAI 'COMPORTEMENT :' GMTYPE.'ELASTIQUE'/NDEC ;
SI (EXIS GMTYPE 'PLASTIQUE') ;
    MESSMAT = CHAI MESSMAT ' ' GMTYPE.'PLASTIQUE' ;
FINSI ;
SI (EXIS GMTYPE 'FLUAGE') ;
    MESSMAT = CHAI MESSMAT ' ' GMTYPE.'FLUAGE' ;
FINSI ;
MESS MESSMAT ;
* MESSAGE VARIATION DES PROPRIETES
SI DISMAT ;
    MESSVAR = CHAI 'PROPRIETES :' 'HOMOGENES PAR MORCEAUX'/NDEC ;
SINON ;
    SI ((DIME PGRAD) > 0) ;
        MESSVAR = CHAI 'PROPRIETES :' 'HETEROGENES A GRADIENT'/NDEC ;
    SINON ;
        MESSVAR = CHAI 'PROPRIETES :' 'HOMOGENES'/NDEC ;
    FINSI ;
FINSI ;
MESS MESSVAR ;
 
* I.1.E - MAILLAGE ET TYPE D'ELEMENTS
MAILLAGE = EXTR MOD_MEC 'MAIL' ;
SUPTAB.'FISSURE' = VIDE 'MAILLAGE' ;
SI (EXIS SUPTAB 'LEVRE_SUPERIEURE') ;
    SUPTAB.'FISSURE' = SUPTAB.'FISSURE' ET SUPTAB.'LEVRE_SUPERIEURE' ;
FINSI ;
SI (EXIS SUPTAB 'LEVRE_INFERIEURE') ;
    SUPTAB.'FISSURE' = SUPTAB.'FISSURE' ET SUPTAB.'LEVRE_INFERIEURE' ;
FINSI ;
* ELEMENTS LINEAIRES OU NONLINEAIRES
ELTYPE = MAILLAGE ELEM 'TYPE' ;
ELLIN = MOTS 'TRI3' 'QUA4' 'CUB8' 'PRI6' 'PYR5' 'TET4' ;
ELQUA = MOTS 'TRI6' 'QUA8' 'CU20' 'PR15' 'PY13' 'TE10' ;
IENLE = LECT ;
REPE I (DIME ELTYPE) ;
    SI (NON (EXIS (ELLIN ET ELQUA) (EXTR ELTYPE &I))) ;
        IENLE = IENLE ET (LECT &I) ;
    FINSI ;
FIN I ;
SI ((DIME IENLE) NEG 0) ;
    MESS 'ATTENTION : IL Y A DES ELEMENTS ATYPIQUES DANS LE MAILLAGE...' ;
    MESS '          LE CALCUL SE POURSUIT QUAND MEME.' ;
    ELTYPE = ENLE ELTYPE IENLE ;
FINSI ;
SI (NON ((EXIS ELLIN ELTYPE 'ET') OU (EXIS ELQUA ELTYPE 'ET'))) ;
    MESS 'ATTENTION : IL Y A A LA FOIS DES ELEMENTS LINEAIRES ET DES' ;
    MESS '          ELEMENTS QUADRATIQUES DANS LE MAILLAGE' ;
    GELDEG = 1.5 ;
SINON ;
    SI (EXIS ELLIN ELTYPE 'ET') ;
        GELDEG = 1 ;
    SINON ;
        GELDEG = 2 ;
    FINSI ;
FINSI ;
* ELEMENTS XFEM OU STANDARD (MASSIF OU COQUE)
SI (EXIS MOD_MEC 'ELEM' 'XQ4R' 'XC8R') ;
    GELTYP = MOT 'XFEM' ;
SINON ;
    EL2D = MOTS 'TRI3' 'TRI6' 'QUA4' 'QUA8' ;
    EL3D = MOTS 'CUB8' 'CU20' 'PRI6' 'PR15' 'PYR5' 'PY13' 'TET4' 'TE10' ;
    SI (NON ((EXIS EL2D ELTYPE 'ET') OU (EXIS EL3D ELTYPE 'ET'))) ;
        MESS 'ERREUR : LE MAILLAGE SUPPORT DU MODELE EST CONSTITUE' ;
        MESS '       D''ELEMENTS 2D ET 3D.' ;
        ERRE 21 ; QUIT G_CAS ;
    SINON ;
        SI ((EGA GDIME 3) ET (EXIS EL2D ELTYPE 'ET')) ;
            GELTYP = 'STD COQUE' ;
        SINON ;
            GELTYP = MOT 'STD MASSIF' ;
        FINSI ;
    FINSI ;
FINSI ;
* MESSAGE
MESSELEM = CHAI 'ELEMENTS :' GELTYP/NDEC ;
MESSELEM = CHAI MESSELEM ' DEGRE' ' ' GELDEG ;
MESS MESSELEM ;
 
* I.1.F - CHARGEMENT
GCHAR = TABL ;
SI (EXIS MOD_TOT 'MATE' 'PRESSION') ;
    SI (EGA GELTYP 'XFEM') ;
        GCHAR.'PRESSION' = 'HORS LEVRES' ;
    SINON ;
        MODCHA = EXTR MOD_TOT 'FORM' 'CHARGEMENT' ;
        MODPRE = EXTR MODCHA 'COMP' 'PRESSION' ;
        MAILPTOT = EXTR MODPRE 'MAIL' ;
        MAILPF = INTE MAILPTOT (SUPTAB.'FISSURE') ;
        SI ((NBEL MAILPF) > 0) ;
            GCHAR.'PRESSION' = 'SUR LEVRES' ;
        SINON ;
            GCHAR.'PRESSION' = 'HORS LEVRES' ;
        FINSI ;
    FINSI ;
SINON ;
    GCHAR.'PRESSION' = 'AUCUN' ;
FINSI ;
SI (EGA GSOLU 'LINEAIRE') ;
    GCHAR.'MECANIQUE' = EXIS SUPTAB 'CHARGEMENTS_MECANIQUES' ;
    SI (EXIS SUPTAB 'TEMPERATURES') ;
        GCHAR.'THERMIQUE' = MOT 'DECOUPLE' ;
    SINON ;
        GCHAR.'THERMIQUE' = 'AUCUN' ;
    FINSI ;
    GCHAR.'DEFORMATIONS IMPOSEES' = EXIS SUPTAB 'DEFORMATIONS_IMPOSEES' ;
SINON ;
    CHAR1 = SUPTAB.'SOLUTION_PASAPAS'.'CHARGEMENT' ;
    GCHAR.'MECANIQUE' = (EXIS CHAR1 'DIMP') OU (EXIS CHAR1 'MECA') ;
    SI (EXIS CHAR1 'T') ;
        GCHAR.'THERMIQUE' = MOT 'DECOUPLE' ;
    SINON ;
        SI ((EXIS CHAR1 'TIMP') OU (EXIS CHAR1 'Q') OU
            (EXIS CHAR1 'TECO') OU (EXIS CHAR1 'TERA')) ;
            GCHAR.'THERMIQUE' = MOT 'COUPLE' ;
        SINON ;
            GCHAR.'THERMIQUE' = 'AUCUN' ;
        FINSI ;
    FINSI ;
    GCHAR.'DEFORMATIONS IMPOSEES' = EXIS CHAR1 'DEFI' ;
FINSI ;
NMESS = 0 ;
MESSCHAR = MOT 'CHARGEMENT :' ;
ICHAR = INDE GCHAR ;
REPE I (DIME ICHAR) ;
    MCHAR = MOT ICHAR.&I ;
    SI (NEG NMESS 0) ;
        MESSCHAR = ' ' ;
    FINSI ;
    SI (EGA (TYPE GCHAR.(MOT ICHAR.&I)) 'LOGIQUE') ;
        SI GCHAR.MCHAR ;
            MESSCHAR = CHAI MESSCHAR MCHAR/NDEC ;
            MESS MESSCHAR ;
            NMESS = NMESS + 1 ;
        FINSI ;
    FINSI ;
    SI (EGA (TYPE GCHAR.MCHAR) 'MOT') ;
        SI (NEG GCHAR.MCHAR 'AUCUN') ;
            MESSCHAR = CHAI MESSCHAR MCHAR/NDEC ;
            MESSCHAR = CHAI MESSCHAR ' ' GCHAR.MCHAR ;
            MESS MESSCHAR ;
            NMESS = NMESS + 1 ;
        FINSI ;
    FINSI ;
FIN I ;
 
* I.1.G - CONTACT ENTRE LEVRES DE FISSURE
BCONT = (EXIS SUPTAB 'MODELE_FISSURE') ;
MESSCONT = CHAI 'CONTACT LEVRES :' BCONT/NDEC ;
MESS MESSCONT ;
 
 
* I.2 - VERIFICATION QUE LE CAS EST TRAITABLE
* *******************************************
 
* I.2.A - VERIFICATION QUE L'ON CONNAIT L'OBJECTIF
SI ((NEG GOBJ 'J') ET (NEG GOBJ 'C*') ET (NEG GOBJ 'C*H')
        ET (NEG GOBJ 'DJ/DA') ET (NEG GOBJ 'J_DYNA')
            ET (NEG GOBJ 'DECOUPLAGE')) ;
    MESS 'ERREUR : ON NE CONNAIT PAS L INTEGRALE SPECIFIEE' ;
    MESS '       A CALCULER' ;
    ERRE 21 ; QUIT G_CAS ;
FINSI ;
 
* I.2.B - VERIFICATION QUE LE MATERIAU EST ISOTROPE
SI (EGA GMTYPE.'ELASTIQUE' (MOT 'ELASTIQUE ANISOTROPE')) ;
    MESS 'ERREUR : ON NE PEUT PAS TRAITER LE CAS DE MATERIAUX ANISOTROPES.' ;
    ERRE 21 ; QUIT G_CAS ;
FINSI ;
 
* I.2.C - VERIFICATION QUE L'ON A PAS A LA FOIS DES INTERFACES ET DES GRADIENTS
*        DE PROPRIETES
SI (DISMAT ET ((DIME PGRAD) > 0)) ;
    MESS 'ERREUR : ON NE PEUT PAS TRAITER LE CAS DE DISCONTINUITES' ;
    MESS '       MATERIELLES DANS UN SOLIDE A GRADIENT DE PROPRIETES' ;
    ERRE 21 ; QUIT G_CAS ;
FINSI ;
 
* I.2.D - INTERFACES (2D SEULEMENT, ET ON DOIT AVOIR AU MOINS DEUX MODELES)
SI DISMAT ;
    SI (EGA GDIME 3) ;
        MESS 'ERREUR : ON NE PEUT PAS ENCORE TRAITER LES DISCONTINUITES DE' ;
        MESS '       PROPRIETES EN 3D' ;
        ERRE 21 ; QUIT G_CAS ;
    FINSI ;
    SI ((DIME SUPTAB.'MODELES_COMPOSITES') < 2) ;
        MESS 'ERREUR : IL FAUT AU MOINS DEUX MODELES POUR' ;
        MESS '       DETERMINER LA LIGNE COMMUNE D''INTERFACE' ;
        ERRE 21 ; QUIT G_CAS ;
    FINSI ;
FINSI ;
 
* I.2.E - OBJECTIF : J
SI (EGA GOBJ 'J') ;
    SI (EXIS GMTYPE 'FLUAGE') ;
        MESS 'ERREUR : ON ATTEND UN COMPORTEMENT ELASTOPLASTIQUE POUR' ;
        MESS '       CALCULER L''INTEGRALE J' ;
        ERRE 21 ; QUIT G_CAS ;
    FINSI ;
FINSI ;
 
* I.2.F - OBJECTIF : C* OU C*H
SI ((EGA GOBJ 'C*') OU (EGA GOBJ 'C*H')) ;
    SI (NON (EXIS GMTYPE 'FLUAGE')) ;
        MESS 'ERREUR : ON ATTEND UN COMPORTEMENT VISCOPLASTIQUE POUR' ;
        MESS '       CALCULER L''INTEGRALE' ' ' GOBJ ;
        ERRE 21 ; QUIT G_CAS ;
    FINSI ;
    SI (NEG GCHAR.'THERMIQUE' 'AUCUN') ;
        MESS 'ERREUR : ON NE PEUT PAS ENCORE CALCULER L''INTEGRALE' ' ' GOBJ ;
        MESS '       DANS LE CAS D''UN CHARGEMENT THERMIQUE' ;
        ERRE 21 ; QUIT G_CAS ;
    FINSI ;
    SI (NEG GHYPO 'HPP') ;
        MESS 'ERREUR : ON NE PEUT PAS ENCORE CALCULER L''INTEGRALE' ' ' GOBJ ;
        MESS '       EN GRANDS DEPLACEMENTS' ;
        ERRE 21 ; QUIT G_CAS ;
    FINSI ;
    SI (EGA GOBJ 'C*H') ;
        REPE I (DIME MVISC) ;
            SI ((NEG MVISC.&I 'NORTON') ET (EXIS MOD_MEC 'MATE' MVISC.&I)) ;
                MESS 'ERREUR : IL FAUT UN MODELE DE FLUAGE NORTON' ;
                MESS '       SEUL POUR CALCULER L INTEGRALE' ' ' GOBJ ;
                                ERRE 21 ; QUIT G_CAS ;
            FINSI ;
        FIN I ;
    FINSI ;
FINSI ;
 
* I.2.G - OBJECTIF : DJ/DA
SI (EGA GOBJ 'DJ/DA') ;
    SI (EXIS GMTYPE 'FLUAGE') ;
        MESS 'ERREUR : ON ATTEND UN COMPORTEMENT ELASTOPLASTIQUE POUR' ;
        MESS '       CALCULER L''INTEGRALE DJ/DA' ;
        ERRE 21 ; QUIT G_CAS ;
    FINSI ;
    SI (DISMAT OU ((DIME PGRAD) > 0)) ;
        MESS 'ERREUR : ON NE PEUT CALCULER L''INTEGRALE DJ/DA QUE' ;
        MESS '       DANS LE CAS DE PROPRIETES HOMOGENES' ;
        ERRE 21 ; QUIT G_CAS ;
    FINSI ;
    SI (NEG GELTYP 'STD MASSIF') ;
        MESS 'ERREUR : ON NE PEUT CALCULER L''INTEGRALE DJ/DA QU''AVEC' ;
        MESS '       DES ELEMENTS STANDARDS MASSIFS' ;
        ERRE 21 ; QUIT G_CAS ;
    FINSI ;
    SI (NEG GHYPO 'HPP') ;
        MESS 'ERREUR : ON NE PEUT CALCULER L''INTEGRALE DJ/DA QUE' ;
        MESS '       SOUS L''HYPOTHESE DES PETITES PERTURBATIONS' ;
        ERRE 21 ; QUIT G_CAS ;
    FINSI ;
    SI (NON (EXIS SUPTAB 'BLOCAGES_MECANIQUES')) ;
        MESS 'ERREUR : IL FAUT DONNER LES BLOCAGES MECANIQUES' ;
        MESS '       POUR CALCULER L''INTEGRALE DJ/DA'
        ERRE 21 ; QUIT G_CAS ;
    FINSI ;
FINSI ;
 
* I.2.H - OBJECTIF : J_DYNA
SI (EGA GOBJ 'J_DYNA') ;
    SI ((EXIS GMTYPE 'PLASTIQUE') OU (EXIS GMTYPE 'FLUAGE')) ;
        MESS 'ERREUR : ON ATTEND UN COMPORTEMENT ELASTIQUE POUR CALCULER' ;
        MESS '       L''INTEGRALE J_DYNA' ;
        ERRE 21 ; QUIT G_CAS ;
    FINSI ;
    SI (NEG GSOLU (MOT 'PASAPAS DYNAMIQUE')) ;
        MESS 'ERREUR : ON ATTEND UNE SOLUTION DE PASAPAS DYNAMIQUE POUR' ;
        MESS '       CALCULER L''INTEGRALE J_DYNA' ;
        ERRE 21 ; QUIT G_CAS ;
    FINSI ;
FINSI ;
 
* I.2.I - OBJECTIF : DECOUPLAGE
SI (EGA GOBJ (MOT 'DECOUPLAGE')) ;
    SI ((EXIS GMTYPE 'PLASTIQUE') OU (EXIS GMTYPE 'FLUAGE')) ;
        MESS 'ERREUR : ON ATTEND UN COMPORTEMENT ELASTIQUE POUR UTILISER' ;
        MESS '       L''OPTION DECOUPLAGE' ;
        ERRE 21 ; QUIT G_CAS ;
    FINSI ;
    SI ((DIME PGRAD) > 0) ;
        MESS 'ERREUR : ON NE PEUT PAS UTILISER L''OPTION DECOUPLAGE POUR UN' ;
        MESS '       SOLIDE A GRADIENT DE PROPRIETES' ;
        ERRE 21 ; QUIT G_CAS ;
    FINSI ;
    SI (EGA GELTYP 'STD COQUE') ;
        MESS 'ERREUR : ON NE PEUT PAS UTILISER L''OPTION DECOUPLAGE DANS LE' ;
        MESS '       CAS DE COQUES' ;
        ERRE 21 ; QUIT G_CAS ;
    FINSI ;
*   SI IL Y A DES INTERFACES IL FAUT VERIFIER QU'ELLES SONT PARALLELES A LA
*   FISSURE, MAIS LA DIRECTION DE CELLE-CI EST OBTENUE APRES AVOIR CALCULE
*   LE CHAMP THETA, DONC CETTE VERIFICATION EST FAITE DANS G_THETA
FINSI ;
 
* I.2.J - CONTACT ENTRE LES LEVRES : UNIQUEMENT AVEC RESO ET XFEM
SI BCONT ;
    SI (NEG GSOLU 'LINEAIRE') ;
        MESS 'ERREUR : CONTACT ENTRE LES LEVRES NON PREVU AVEC UNE SOLUTION' ;
        MESS '       DE PASAPAS' ;
        ERRE 21 ; QUIT G_CAS ;
    FINSI ;
    SI (NEG GELTYP (MOT 'XFEM')) ;
        MESS 'ERREUR : CONTACT ENTRE LES LEVRES NON PREVU AVEC DES ELEMENTS' ;
        MESS '       STANDARDS' ;
        ERRE 21 ; QUIT G_CAS ;
    FINSI ;
FINSI ;
 
* I.2.K - ON VERIFIE QU'IL N'Y A PAS DE GRADIENT DE COEFFICIENT DE POISSON
SI (EXIS PGRAD 'NU') ;
    MESS 'ERREUR : ON NE SAIT PAS ENCORE TRAITER LA PRESENCE D''UN GRADIENT' ;
    MESS '       DE COEFFICIENT DE POISSON' ;
    ERRE 21 ; QUIT G_CAS ;
FINSI ;
 
 
* I.3 - VERIFICATION DES DONNEES NECESSAIRES
* ******************************************
 
* I.3.A - FRONT DE FISSURE
SI (NON (EXIS SUPTAB 'FRONT_FISSURE')) ;
    MESS 'ERREUR : ON VEUT LE FRONT DE LA FISSURE' ;
    ERRE 21 ; QUIT G_CAS ;
FINSI ;
M_FRONT = SUPTAB.'FRONT_FISSURE' ;
SI ((NEG (TYPE M_FRONT) 'MAILLAGE') ET (NEG (TYPE M_FRONT) 'POINT')) ;
    MESS 'ERREUR : LE FRONT DE FISSURE DOIT ETRE DE TYPE MAILLAGE OU POINT' ;
    ERRE 21 ; QUIT G_CAS ;
FINSI ;
SI (EGA (TYPE M_FRONT) 'POINT') ;
    M_FRONT = MANU 'POI1' M_FRONT ;
FINSI ;
M_FRONT = M_FRONT COUL 'OLIV' ;
 
* I.3.B - FISSURE_2 ET SON FRONT POUR DJ/DA
SI (EGA GOBJ 'DJ/DA') ;
    SI ((EXIS SUPTAB 'FISSURE_2') OU (EXIS SUPTAB 'FRONT_FISSURE_2')) ;
        SI (NON (EXIS SUPTAB 'FISSURE_2')) ;
            MESS 'ERREUR : ON VEUT AUSSI LA FISSURE 2 POUR CALCULER' ;
            MESS '       LES TERMES CROISES DE LA MATRICE' ;
            ERRE 21 ; QUIT G_CAS ;
        FINSI ;
        SI (NON (EXIS SUPTAB 'FRONT_FISSURE_2')) ;
            MESS 'ERREUR : ON VEUT AUSSI LE FROND DE LA FISSURE 2 POUR' ;
            MESS '       CALCULER LES TERMES CROISES DE LA MATRICE' ;
            ERRE 21 ; QUIT G_CAS ;
        FINSI ;
    SINON ;
        SI (EGA SUPTAB.'COUCHE' 0) ;
            MESS 'ERREUR : LE NOMBRE DE COUCHES DOIT ETRE SUPERIEUR A' ;
            MESS '       0 POUR LE CALCUL DU TERME PRINCIPAL DJi/DAi' ;
            ERRE 21 ; QUIT G_CAS ;
        FINSI ;
    FINSI ;
FINSI ;
 
* I.3.C - PRESSION OBLIGATOIRE SI CALCUL LINEAIRE ET MODELE DE PRESSION
SI ((NEG GCHAR.'PRESSION' 'AUCUN') ET (EGA GSOLU 'LINEAIRE')) ;
    SI (NON (EXIS SUPTAB 'PRESSION')) ;
        MESS 'ERREUR : IL MANQUE LA DONNEE DU CHAMP DE PRESSION' ;
        ERRE 21 ; QUIT G_CAS ;
    FINSI ;
FINSI ;
 
* I.3.D - VERIFICATION QU'ON A LE CHARGEMENT SI SOLUTION LINEAIRE
SI (EGA GSOLU 'LINEAIRE') ;
    SI ((EGA GCHAR.'THERMIQUE' 'AUCUN')
            ET (EGA GCHAR.'PRESSION' 'AUCUN')
            ET (NON GCHAR.'MECANIQUE')
            ET (NON GCHAR.'DEFORMATIONS IMPOSEES')) ;
        MESS 'ERREUR : IL FAUT LES CHARGEMENTS APPLIQUES :' ;
        MESS '       MECANIQUES, THERMIQUES OU LES DEUX' ;
        ERRE 21 ; QUIT G_CAS ;
    FINSI ;
FINSI ;
 
* I.3.E - DEFINITION DE LA FISSURE NECESSAIRE
SI (EGA GELTYP (MOT 'XFEM')) ;
    SI (NON ((EXIS SUPTAB 'PSI') ET (EXIS SUPTAB 'PHI'))) ;
        MESS 'ERREUR : ON VEUT PSI et PHI LEVELSET DE LA FISSURE' ;
        ERRE 21 ; QUIT G_CAS ;
    FINSI ;
SINON ;
    SI (EGA GT_OBJ 'DECOUPLAGE') ;
        SI (NON ((EXIS SUPTAB 'LEVRE_SUPERIEURE') ET
                    (EXIS SUPTAB 'LEVRE_INFERIEURE'))) ;
            MESS 'ERREUR : L''OBJECTIF DECOUPLAGE NECESSITE LES DEUX' ;
            MESS '       LEVRES DE LA FISSURE' ;
            ERRE 21 ; QUIT G_CAS ;
        FINSI ;
    SINON ;
        SI (NON ((EXIS SUPTAB 'LEVRE_SUPERIEURE') OU
                    (EXIS SUPTAB 'LEVRE_INFERIEURE'))) ;
            MESS 'ERREUR : IL FAUT DONNER LA FISSURE' ;
            MESS '(LEVRE_SUPERIEURE OU LEVRE_INFERIEURE OU LES 2)' ;
            ERRE 21 ; QUIT G_CAS ;
        FINSI ;
    FINSI ;
FINSI ;
 
* I.3.F - ON VERIFIE QU'ON A OU BIEN LE NOMBRE DE COUCHES OU BIEN CHAM_THETA
SI (NON ((EXIS SUPTAB 'COUCHE') OU (EXIS SUPTAB 'CHAMP_THETA'))) ;
    MESS 'ERREUR : ON VEUT LE NOMBRE DE COUCHES D ELEMENTS' ;
    MESS '       AUTOUR DE LA FISSURE QUI SE DEPLACE' ;
    MESS '       OU LE CHAMP_THETA' ;
    MESS '       POUR SIMULER LA PROPAGATION DE LA FISSURE' ;
    ERRE 21 ; QUIT G_CAS ;
FINSI ;
 
* I.3.G - ON NE CALCULE PAS CORRECTEMENT LE CAS THERMIQUE AVEC DES COQUES
*        DONC POUR L'INSTANT ON DEBRANCHE L'OPTION
SI ((EGA GELTYP 'STD COQUE') ET (NEG GCHAR.'THERMIQUE' 'AUCUN')) ;
    MESS 'ERREUR : UNE CORRECTION EST NECESSAIRE POUR TRAITER UN CHARGEMENT' ;
    MESS '       THERMIQUE AVEC LES ELEMENTS COQUE' ;
    ERRE 21 ; QUIT G_CAS ;
FINSI ;
 
 
* =============================================================================
 
* II - STOCKAGE DES INFORMATIONS POUR TRANSFERT A G_THETA
* -------------------------------------------------------
 
* II.1 - OBJETS ET INFOS UTILES DANS OBJUTI
* *****************************************
 
OBJUTI = TABL ;
OBJUTI.'DIMENSION' = GDIME ;
OBJUTI.'MODE' = GMODE ;
OBJUTI.'MOD_TOT' = MOD_TOT ;
OBJUTI.'MOD_MEC' = MOD_MEC ;
OBJUTI.'MAT_MEC' = MAT_MEC ;
OBJUTI.'MAILLAGE' = MAILLAGE ;
OBJUTI.'FRONT' = M_FRONT ;
SI (NEG GSOLU 'LINEAIRE') ;
    OBJUTI.'WTAB' = WTAB ;
FINSI ;
 
* II.2 - BOOLEENS DANS BOOL
* *************************
 
* OBJECTIF
BOOL.'J' = EGA SUPTAB.'OBJECTIF' 'J' ;
BOOL.'C*' = EGA SUPTAB.'OBJECTIF' 'C*' ;
BOOL.'C*H' = EGA SUPTAB.'OBJECTIF' 'C*H' ;
BOOL.'DJ/DA' = EGA SUPTAB.'OBJECTIF' 'DJ/DA' ;
BOOL.'J_DYNA' = EGA SUPTAB.'OBJECTIF' 'J_DYNA' ;
BOOL.'DECOUPLAGE' = EGA SUPTAB.'OBJECTIF' (MOT 'DECOUPLAGE') ;
* PASAPAS
BOOL.'PASAPAS' = NEG GSOLU 'LINEAIRE' ;
* GRANDS DEPLACEMENTS
BOOL.'GRANDS_DEP' = EGA GHYPO 'GRANDS_DEPLACEMENTS' ;
* MODELE DE PRESSION
BOOL.'MODE_PRES' = NEG GCHAR.'PRESSION' 'AUCUN' ;
* PRESSION SUR FISSURE
BOOL.'PRES_FISS' = EGA GCHAR.'PRESSION' 'SUR LEVRES' ;
* CHARGEMENTS THERMIQUES
BOOL.'THER_DECO' = EGA GCHAR.'THERMIQUE' (MOT 'DECOUPLE') ;
BOOL.'THER' = NEG GCHAR.'THERMIQUE' 'AUCUN' ;
* DEFORMATIONS IMPOSEES
BOOL.'DEF_IMP' = GCHAR.'DEFORMATIONS IMPOSEES' ;
* ORDRE DES ELEMENTS
BOOL.'EL_LIN' = EGA GELDEG 1 ;
BOOL.'EL_QUA' = EGA GELDEG 2 ;
* ELEMENTS XFEM
BOOL.'XFEM' = EXIS MOD_MEC 'ELEM' 'XQ4R' 'XC8R' ;
* ELEMENTS COQUE
BOOL.'COQ' = EGA GELTYP 'STD COQUE' ;
* ROTATION DE CORPS RIGIDE
BOOL.'ROT_RIG' = EXIS SUPTAB 'ROTATION_RIGIDIFIANTE' ;
* PERSO1
SI (EXIS SUPTAB 'PERSO1') ;
    BOOL.'PERSO1' = SUPTAB.'PERSO1' ;
SINON ;
    BOOL.'PERSO1' = FAUX ;
FINSI ;
* FROTTEMENT ENTRE LEVRES DE FISSURE
BOOL.'FROT' = EXIS SUPTAB 'MODELE_FISSURE' ;
 
 
* SAUT DE LIGNE AVANT DE RETOURNER DANS G_THETA
SAUT 'LIGNE' ;
 
FINP OBJUTI ;
 
 
 
 

© Cast3M 2003 - Tous droits réservés.
Mentions légales