* fichier :  test_proi.dgibi
******************************************************************
* Projection de champs                                           *
*    Syntaxe :                                                   *
*    1) Champ par point aux noeuds du maillage GEO2              *
*    2) Champ par element aux noeuds ou pt d'integration de MOD2 *
*    3) Projection au sens des moindres carrés                   *
*    4) Projection d'un champ par point et calcul de derivees    *
*   Note : ce test verifie la fonctionnalite, pas le resultat    *
*                                                                *
* G. Folzan, Octobre 2013                                        *
*                                                                *
******************************************************************

'OPTI' 'DIME' 2 'ELEM' 'QUA4' ;
*                                                                    *

* Creation des points
PT_A       = 0   0 ;
PT_B       = 0   1 ;
PT_G       = 0   0 ;
PT_H       = 0   1 ;
PT_C       = 1   1 ;
PT_D       = 1   0 ;
PT_E       = 0.5 0 ;
PT_F       = 0.5 1 ;

                    
*Creation des cotes du maillage                                       *
SEG_AB       = 'DROI' 3 PT_A PT_B ;
SEG_GH       = 'DROI' 3 PT_G PT_H ;
SEG_BF       = 'DROI' 3 PT_B PT_F ;
SEG_FC       = 'DROI' 3 PT_F PT_C ;
SEG_CD       = 'DROI' 3 PT_C PT_D ;
SEG_DE       = 'DROI' 3 PT_D PT_E ;
SEG_EA       = 'DROI' 3 PT_E PT_A ;
SEG_FE       = 'DROI' 3 PT_F PT_E ;


*Creation des contours
CONT_1       = SEG_AB et SEG_BF et SEG_FE et SEG_EA;
CONT_2       = SEG_FC et SEG_CD et SEG_DE et SEG_FE;


*Creation  de maillages differents
*     1) 10 couches d'elements
*     2) 13 couches d'elements
*     3) non structure: quadrangle et triangle

MAIL_1       = SEG_AB 'TRAN' 10 PT_D ;

MAIL_2       = (SEG_GH 'TRAN' 13 PT_D) 'COUL' 'VERT' ;

mail3a       = surf cont_1;
'OPTI' 'DIME' 2 'ELEM' 'TRI3' ;
mail3b       = surf cont_2;
MAIL_3       = (MAIL3A et MAIL3B) coul 'ROUGE';

*                                                                      *
*'TRAC' (MAIL_1 'ET' MAIL_2 'ET' MAIL_3) ;
*                                                                      *

*Creation de modeles associes a chaque maillage

MOD1     = 'MODE' MAIL_1 'THERMIQUE' ;
MOD2     = 'MODE' MAIL_2 'MECANIQUE' ;
MOD3     = 'MODE' MAIL_3 'MECANIQUE' ;


*Creation du champ de depart: champ par element defini au noeud
CHPO1    = 'MANU' 'CHPO' MAIL_1 1 'T' 100. ;
CHAM1    = 'CHAN' 'CHAM' CHPO1 MOD1 ;
*                                                              

*Test syntaxe 1
CHPO2 = 'PROI' CHAM1 MAIL_2;
CHPO2 = 'PROI' MAIL_2 CHAM1;
CHPO3 = 'PROI' CHAM1 MAIL_3;
'MENAGE' 'OBLI';

*Test syntaxe 2
CHAM2    = 'PROI' MOD2 CHAM1;
CHAM2    = 'PROI' CHAM1 MOD2;
CHAM2    = 'PROI' 'GRAVITE' MOD2 CHAM1 ;
CHAM2    = 'PROI' MOD2 'RIGIDITE' CHAM1 ;
CHAM2    = 'PROI' CHAM1 MOD2 'MASSES' ;
CHAM2    = 'PROI' CHAM1 'STRESSES' MOD2;
CHAM3    = 'PROI' MOD3 CHAM1;
CHAM3    = 'PROI' MOD3 CHAM1 'GRAVITE' ;
CHAM3    = 'PROI' MOD3 CHAM1 'RIGIDITE' ;
CHAM3    = 'PROI' MOD3 CHAM1 'MASSES' ;
CHAM3    = 'PROI' MOD3 CHAM1 'STRESSES' ;
'MENAGE' 'OBLI';

*Test syntaxe 3
CHAM2    = 'PROI' 'MINI' MOD2 CHAM1;
CHAM2    = 'PROI' 'MINI' CHAM1 MOD2;
CHAM2    = 'PROI' 'MINI' MOD2 CHAM1 5;
CHAM3    = 'PROI' MOD3 CHAM1 'MINI' ;
CHAM3    = 'PROI' 5 MOD3 'MINI' CHAM1 ;
'MENAGE' 'OBLI';

*Test syntaxe 4
CHPO2 = PROI 'POLY' MAIL_1 MAIL_2 CHPO1 1 'PLAN';
CHPO2 = PROI 'POLY' CHPO1 MAIL_1 MAIL_2 1 'PLAN';
CHPO2 = PROI 'POLY' MAIL_1 MAIL_2 CHPO1 1 'PLAN' 'POID' 5 9;
CHPO2 = PROI 'POLY' MAIL_1 CHPO1 MAIL_2 1 'AXIS';
'MENAGE' 'OBLI';

*  TEST DE PROJECTION SUR DES ELEMENTS COQUES                          *
OPTI DIME 3 elem cub8;

*Creation du maillage et champ de depart (element massif)
vec1 = 0. 0. 0.1;
vol1 = mail_1 VOLU 1 'TRAN' vec1;
vec1 = 0. 0. -0.05;
depl vol1 'PLUS' vec1;
mod1b =  'MODE' vol1 'THERMIQUE' ;
CHPO1    = 'MANU' 'CHPO' vol1 1 'T' 100. ;
CHAM1    = 'CHAN' 'CHAM' CHPO1 MOD1b ;

*Creation du maillage de reception (element coque)
MOD2b    = 'MODE' MAIL_2 'MECANIQUE' 'ELASTIQUE' 'ISOTROPE' 'COQ4';
CAR2 = CARA MOD2b EPAI 0.1 ;

*Test syntaxe 2
CHAM2    = 'PROI' CAR2 MOD2b CHAM1 ;
CHAM2    = 'PROI' MOD2b CAR2 CHAM1 ;
CHAM2    = 'PROI' CHAM1 MOD2b CAR2 ;
CHAM2    = 'PROI' CAR2 MOD2b CHAM1 'GRAVITE' ;
CHAM2    = 'PROI' CAR2 MOD2b CHAM1 'RIGIDITE' ;
CHAM2    = 'PROI' CAR2 MOD2b CHAM1 'STRESSES' ;
CHAM2    = 'PROI' CAR2 MOD2b CHAM1 'STRESSES' ;
'MENAGE' 'OBLI';

*Test syntaxe 3
CHAM2    = 'PROI' 'MINI' CAR2 MOD2b CHAM1 ;
CHAM2    = 'PROI' 'MINI' MOD2b CHAM1 CAR2 ;
CHAM2    = 'PROI' 'MINI' CHAM1 CAR2 MOD2b ;
CHAM2    = 'PROI' 'MINI' CAR2 MOD2b CHAM1 5; 
'MENAGE' 'OBLI';

fin;