* 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;
© Cast3M 2003 - All rights reserved.
Disclaimer