Télécharger g_cas.procedur

Retour à la liste

* G_CAS     PROCEDUR  JB251061  20/11/04    21:17:46     10763          
* =============================================================================
*               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 ;
                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
OBJPOSS = MOTS 'J' 'C*' 'C*H' 'DJ/DA' 'J_DYNA' 'DECOUPLAGE' ;
SI (NON (EXIS OBJPOSS GOBJ)) ;
        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 ;
                        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 ;
 
 
* =============================================================================
 
* 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