Accueil | Cast3M

Télécharger g_defaut_circonferentiel_1.dgibi
* fichier : g_defaut_circonferentiel_1.dgibi
****************************************************
****************************************************
*                                                  *
*      VERIFICATION DE LA PROCEDURE G_THETA        *
*    POUR UN DEFAUT CIRCONFERENTIEL DEBOUCHANT     *
* DANS UN TUYAU (FRONT FERME, 2 LEVRES MODELISEES) *
*                        -                         *
*  COMPARAISON ENTRE LE CALCUL 3D ET UN CALCUL     *
*                2D AXISYMETRIQUE                  *
*                                                  *
****************************************************
 
* I - CALCUL 3D COMPLET
* ---------------------
 
OPTI 'DIME' 3 'ELEM' 'CU20' 'DENS' 1. ;
BTRAC = FAUX ;
 
* MAILLAGE DE LA TRANCHE
L1 = 10. ;
L2 = 20. ;
H1 = 5. ;
SURF1 = (DROI (L1 0. 0.) (L2 0. 0.)) TRAN (0. H1 0.) ;
ANG1 = 360. ;
VOLU1 = SURF1 VOLU 10 'ROTA' ANG1 (0. 1. 0.) (0. 0. 0.) ;
ELIM VOLU1 1.E-10 ;
VOLU2 = VOLU1 SYME 'PLAN' (0. 0. 0.) (1. 0. 0.) (0. 0. 1.) ;
OMEGA1 = VOLU1 ET VOLU2 ;
 
* DEFINITION DE LA FISSURE ET FUSION DES NOEUDS COINCIDENTS
* PARTIE SUP
ENVE1 = ENVE VOLU1 ;
POIN1 = (COOR 2 ENVE1) POIN 'MINI' ;
ELEM1 = ENVE1 ELEM 'APPUYE' POIN1 ;
R1 = (((COOR 1 ELEM1)**2) + ((COOR 3 ELEM1)**2))**0.5 ;
POIN1 = R1 POIN 'EGINFE' ((L1 + L2)/2) ;
LVSUP = ELEM1 ELEM 'APPUYE' POIN1 ;
LIGSUP = ELEM1 DIFF LVSUP ;
SI BTRAC ;
        TRAC 'CACH' (VOLU1 ET (LVSUP COUL 'ROUG') ET (LIGSUP COUL 'VERT')) ;
FINSI ;
* PARTIE INF
ENVE2 = ENVE VOLU2 ;
POIN2 = (COOR 2 ENVE2) POIN 'MAXI' ;
ELEM2 = ENVE2 ELEM 'APPUYE' POIN2 ;
R2 = (((COOR 1 ELEM2)**2) + ((COOR 3 ELEM2)**2))**0.5 ;
POIN2 = R2 POIN 'EGINFE' ((L1 + L2)/2) ;
LVINF = ELEM2 ELEM 'APPUYE' POIN2 ;
LIGINF = ELEM2 DIFF LVINF ;
SI BTRAC ;
        TRAC 'CACH' (VOLU2 ET (LVINF COUL 'ROUG') ET (LIGINF COUL 'VERT')) ;
FINSI ;
* ELIMINATION
ELIM LIGSUP LIGINF 1.E-10 ;
SI BTRAC ;
        TRAC 'CACH' OMEGA1 ;
FINSI ;
* DEFINITION DU FRONT DE FISSURE
FRONFISS = (CONT LIGSUP) INTE (CONT LVSUP) ;
 
* MODELE ET MATERIAU
MOD1 = MODE OMEGA1 'MECANIQUE' 'ELASTIQUE' ;
MAT1 = MATE MOD1 'YOUN' 2.E5 'NU' 0.3 ;
RIG1 = RIGI MOD1 MAT1 ;
 
* CL
PSUP = (COOR 2 OMEGA1) POIN 'MAXI' ;
BLOQ1 = BLOQ (MOTS 'UX' 'UZ') PSUP ;
BLOQ2 = BLOQ 'UY' PSUP ;
PINF = (COOR 2 OMEGA1) POIN 'MINI' ;
BLOQ3 = BLOQ (MOTS 'UX' 'UZ') PINF ;
BLOQ4 = BLOQ 'UY' PINF ;
BLOQ0 = BLOQ1 ET BLOQ2 ET BLOQ3 ET BLOQ4 ;
DEP1 = DEPI BLOQ2 1. ;
DEP2 = DEPI BLOQ4 (0.-1.) ;
DEP0 = DEP1 ET DEP2 ;
 
* RESO
U1 = RESO (RIG1 ET BLOQ0) DEP0 ;
 
* G_THETA
SUPTAB = TABL ;
SUPTAB.'MODELE' = MOD1 ;
SUPTAB.'CARACTERISTIQUES' = MAT1 ;
SUPTAB.'BLOCAGES_MECANIQUES' = BLOQ0 ;
SUPTAB.'CHARGEMENTS_MECANIQUES' = DEP0 ;
SUPTAB.'OBJECTIF' = 'J' ;
SUPTAB.'LEVRE_SUPERIEURE' = LVSUP ;
SUPTAB.'LEVRE_INFERIEURE' = LVINF ;
SUPTAB.'COUCHE' = 4 ;
SUPTAB.'FRONT_FISSURE' = FRONFISS ;
SUPTAB.'SOLUTION_RESO' = U1 ;
 
G_THETA SUPTAB ;
G3D = PROG ;
IND1 = INDE SUPTAB.'RESULTATS' ;
NFRON = DIME IND1 ;
REPE IFRON NFRON ;
        PFRON = IND1.&IFRON ;
        G3D = G3D ET (PROG SUPTAB.'RESULTATS'.PFRON) ;
FIN IFRON ;
 
 
* II - CALCUL 2D AXISYMETRIQUE
* ----------------------------
 
OPTI 'DIME' 2 'MODE' 'AXIS' 'ELEM' 'QUA8' ;
 
SURF2 = SURF1 SYME 'DROI' (0. 0.) (1. 0.) ;
OMEGA2 = SURF1 ET SURF2 ;
 
* DEFINITION DE LA FISSURE ET FUSION DES NOEUDS COINCIDENTS
C1 = CONT SURF1 ;
POIN1 = (COOR 2 C1) POIN 'MINI' ;
ELEM1 = C1 ELEM 'APPUYE' POIN1 ;
POIN1 = (COOR 1 POIN1) POIN 'EGINFE' ((L1 + L2)/2) ;
LVSUP = C1 ELEM 'APPUYE' POIN1 ;
LIGSUP = ELEM1 DIFF LVSUP ;
C2 = CONT SURF2 ;
POIN2 = (COOR 2 C2) POIN 'MAXI' ;
ELEM2 = C2 ELEM 'APPUYE' POIN2 ;
POIN2 = (COOR 1 POIN2) POIN 'EGINFE' ((L1 + L2)/2) ;
LVINF = C2 ELEM 'APPUYE' POIN2 ;
LIGINF = ELEM2 DIFF LVINF ;
ELIM LIGSUP LIGINF 1.E-10 ;
SI BTRAC ;
        TRAC (OMEGA2 ET (LVSUP COUL 'ROUG')) ;
FINSI ;
 
* DEFINITION DU FRONT DE FISSURE
FRONFISS = ((COOR 1 LVSUP) POIN 'MAXI') POIN 1 ;
 
* MODELE ET MATERIAU
MOD2 = MODE OMEGA2 'MECANIQUE' 'ELASTIQUE' ;
MAT2 = MATE MOD2 'YOUN' 2.E5 'NU' 0.3 ;
RIG2 = RIGI MOD2 MAT2 ;
 
* CL
PSUP = (COOR 2 OMEGA2) POIN 'MAXI' ;
BLOQ1 = BLOQ 'UR' PSUP ;
BLOQ2 = BLOQ 'UZ' PSUP ;
PINF = (COOR 2 OMEGA2) POIN 'MINI' ;
BLOQ3 = BLOQ 'UR' PINF ;
BLOQ4 = BLOQ 'UZ' PINF ;
BLOQ0 = BLOQ1 ET BLOQ2 ET BLOQ3 ET BLOQ4 ;
DEP1 = DEPI BLOQ2 1. ;
DEP2 = DEPI BLOQ4 (0.-1.) ;
DEP0 = DEP1 ET DEP2 ;
 
* RESO
U2 = RESO (RIG2 ET BLOQ0) DEP0 ;
 
* G_THETA
SUPTAB = TABL ;
SUPTAB.'MODELE' = MOD2 ;
SUPTAB.'CARACTERISTIQUES' = MAT2 ;
SUPTAB.'BLOCAGES_MECANIQUES' = BLOQ0 ;
SUPTAB.'CHARGEMENTS_MECANIQUES' = DEP0 ;
SUPTAB.'OBJECTIF' = 'J' ;
SUPTAB.'LEVRE_SUPERIEURE' = LVSUP ;
SUPTAB.'LEVRE_INFERIEURE' = LVINF ;
SUPTAB.'COUCHE' = 4 ;
SUPTAB.'FRONT_FISSURE' = FRONFISS ;
SUPTAB.'SOLUTION_RESO' = U2 ;
 
G_THETA SUPTAB ;
G2D = PROG NFRON * SUPTAB.'RESULTATS' ;
 
 
* III - COMPARAISON
* -----------------
 
CRIT1 = MAXI (G3D - G2D) 'ABS' ;
CMAX = MAXI G2D 'ABS' ;
SI (CRIT1 > (CMAX*5.E-3)) ;
        MESS 'LA SOLUTION 3D EST TROP ELOIGNEE DE LA SOLUTION 2D AXIS' ;
        ERRE ;
FINSI ;
 
FIN ;