Download visufour1.dgibi

Back to the list

* fichier : visufour1.dgibi
************************************************************************
************************************************************************
*                                                                      *
*             Exemple d'utilisation de four2tri                        *
*         Visualisation 3D de résultats de calcul Fourier              *
*               avec possibilité de recombinaison des harmoniques      *
*                                                                      *
*                      D. Combescure - Janvier 2007                    *
*                      Laboratoire DYN - CEA Saclay                    *
************************************************************************
************************************************************************
*
*FLAGVISU = VRAI;
*FLAGLONG = VRAI;
*
FLAGVISU = FAUX;
FLAGLONG = FAUX;
*
nf = 1;
opti dime 2 elem tri3 mode four nf;
*
p1 = 10 0.01;
p2 = 11 0.01;
p3 = 20 0.01;
p4 = 30 0.01;
p5 = 30 1;
v1 = 0 1;
*
lig1 = d 3 p1 p2;
mesh2d = tran 3 lig1 v1;
opti elem seg2;
lig2 = d 3 p2 p3;
mod1 = modl mesh2d mecanique elastique;
mod2 = modl lig2 mecanique elastique coq2;
mod3 = modl (p5 et p4) mecanique elastique cerc;
* CAS 1 modèle
*modt = mod1;
*modt = mod2;
* Cas 2 modèles
modt = mod1 et mod2;
*
opti dime 2 elem tri3 mode four nf;
U2D1 = manu chpo (extr (mod1) mail) 
     3 UR 1. UT -1. UZ 1. ;
F2D1 = manu chpo (extr (mod1) mail) 
     3 FR 0.001 FT 0.001 FZ 1. ;
U2D2 = manu chpo (extr (mod2) mail) 
     3 UR 1. UT -1. UZ 1. ;
F2D2 = manu chpo (extr (mod2) mail) 
     3 FR 0.001 FT 0.001 FZ 1.;
U2D = manu chpo (extr (mod1 et mod2) mail) 
     3 UR 1. UT -1. UZ 1. ;
F2D = manu chpo (extr (mod1 et mod2) mail) 
     3 FR 0.001 FT 0.001 FZ 1. ;
opti mode four 0;
PRE1 = manu chpo (extr (mod1 et mod2) mail) 
     2 P 0.001 EPAI 0.001 ;
*
*opti trac open;
*
* Ensemble
*
*
tab1 = table;
mesh3d1 = four2tri (mod1 et mod2)
        tab1 (prog 0. pas 10. 360.);
SI FLAGVISU;
 trac mesh3D1;
FINSI;
*
tab1.'DEPLACEMENTS' = TABLE;
tab1.'DEPLACEMENTS'. 1 = U2D;
mesh3d = four2tri (mod1 et mod2) tab1 nf ;
dep3D = (tab1.deplacements_3d. 1);
SI FLAGVISU;
 trac (defo mesh3D dep3D) cach;
FINSI;
*
oubli tab1 'DEPLACEMENTS';
tab1.'CHAMPS_POINTS_SYMETRIQUES' = TABLE;
tab1.'CHAMPS_POINTS_SYMETRIQUES'. 1 = PRE1;
mesh3d = four2tri (mod1 et mod2) tab1 0 ;
*
SI FLAGVISU;
 trac mesh3d (exco (tab1.CHAMPS_POINTS_SYME_3D. 1) P);
FINSI;
*
opti mode four 1;
FPONC = manu chpo p3 1 FR 1. nature discrete;
opti mode four 2;
FPONC = FPONC et (manu chpo p3 1 FR 1. nature discrete);
opti mode four 3;
FPONC = FPONC et (manu chpo p3 1 FR 1. nature discrete);
opti mode four 4;
FPONC = FPONC et (manu chpo p3 1 FR 1. nature discrete);
tab1 = table;
tab1.'EFFORTS' = TABLE;
tab1.'EFFORTS'. 1 = FPONC;
mesh3d = four2tri (mod1 et mod2) tab1 3
       (prog 0. pas 10. 360.);
FVECT = vecteur (tab1. 'EFFORTS_3D'. 1) FX FY FZ;
SI FLAGVISU;
trac mesh3d FVECT;
FINSI;
*
FIN;
*

 

© Cast3M 2003 - All rights reserved.
Disclaimer