Accueil | Cast3M

* fichier : cubp-cubt.dgibi
************************************************************************
************************************************************************
************************************************************************
* NOM         : CUBP-CUBT
* DESCRIPTION : Utilisation des operateurs de construction de cubiques
*
*
*
* LANGAGE     : GIBIANE-CAST3M
* AUTEUR      : Stephane GOUNAND (CEA/DES/ISAS/DM2S/SEMT/LTA)
*               mel : stephane.gounand@cea.fr
**********************************************************************
* VERSION    : v1, 14/10/2024, version initiale
* HISTORIQUE : v1, 14/10/2024, création
* HISTORIQUE :
* HISTORIQUE :
************************************************************************
*
*
interact= FAUX ;
graph = faux ;
*
'SI' ('NON' interact) ;
  'OPTION' 'TRAC' 'PSC' 'ECHO' 0 ;
'SINON' ;
  'OPTION' 'TRAC' 'X' ;
'FINSI' ;
 
'OPTI' 'DENS' 0.1 ;
*
************************************************************************
*                     E N   D I M E N S I O N   2                      *
************************************************************************
 
OPTI 'DIME' 2 ;
p1 = 0. 0. ;
p2 = 1. 0. ;
p3 = 1. 1. ;
p4 = 1. 0.5 ;
p5 = 2. 0.5 ;
nl = 10 ;
 
* Éléments de maillage SEG2
OPTI 'ELEM' 'SEG2' ;
l1 = 'DROI' p1 p2 ;
l2 = 'DROI' p2 p3 ;
l3 = 'DROI' p4 p5 ;
*
ll1   = l1 'ET' l2 ;
*'OPTI' impi 1 ;
v1 = 'MOIN' p2 p1 ; v2 = 'MOIN' p3 p2 ;
lct1  = 'CUBT' nl p1 v1 v2 p3 'COUL' 'OLIV' ;
lct1u  = 'CUBT' nl p1 v1 v2 p3 'UNIF' 'COUL' 'ROUG' ;
*
'SI' graph ;
   TRAC 'QUAL' 'NOEU' (ll1 'ET' lct1 'ET' lct1u) 'TITR'
   'Lignes et cubiques en angle SEG2' ;
'FINS' ;
ll2 = l1 'ET' l3 ;
lcp1  = 'CUBP' nl p1 p2 p4 p5 'COUL' 'OLIV' ;
lcp1u  = 'CUBP' nl p1 p2 p4 p5 'UNIF' 'COUL' 'ROUG' ;
'SI' graph ;
   TRAC 'QUAL' 'NOEU' (ll2 'ET' lcp1 'ET' lcp1u) 'TITR'
   'Lignes et cubiques en parallele SEG2' ;
'FINS' ;
 
* Éléments de maillage SEG3
OPTI 'ELEM' 'SEG3' ;
lct2  = 'CUBT' nl p1 v1 v2 p3 'COUL' 'OLIV' ;
lct2u  = 'CUBT' nl p1 v1 v2 p3 'UNIF' 'COUL' 'ROUG' ;
*
'SI' graph ;
   TRAC 'QUAL' 'NOEU' (ll1 'ET' lct2 'ET' lct2u) 'TITR'
   'Lignes et cubiques en angle SEG3' ;
'FINS' ;
lcp2  = 'CUBP' nl p1 p2 p4 p5 'COUL' 'OLIV' ;
lcp2u  = 'CUBP' nl p1 p2 p4 p5 'UNIF' 'COUL' 'ROUG' ;
'SI' graph ;
   TRAC 'QUAL' 'NOEU' (ll2 'ET' lcp2 'ET' lcp2u) 'TITR'
   'Lignes et cubiques en parallele SEG3' ;
'FINS' ;
*
mlct1 = 'MESU' lct1 ;
mlct1u = 'MESU' lct1u ;
mlct2 = 'MESU' lct2 ;
mlct2u = 'MESU' lct2u ;
lmlt = 'PROG' mlct1 mlct1u mlct2 mlct2u ;
milmlt = 'MINI' lmlt ;
malmlt = 'MAXI' lmlt ;
*
* Test sur la longueur de la courbe a defaut de mieux
*
lok = vrai ;
'SI' (milmlt '<' 1.48) ; lok = lok 'ET' faux ; 'FINS' ;
'SI' (malmlt '>' 1.49) ; lok = lok 'ET' faux ; 'FINS' ;
*
mlcp1 = 'MESU' lcp1 ;
mlcp1u = 'MESU' lcp1u ;
mlcp2 = 'MESU' lcp2 ;
mlcp2u = 'MESU' lcp2u ;
lmlp = 'PROG' mlcp1 mlcp1u mlcp2 mlcp2u ;
milmlp = 'MINI' lmlp ;
malmlp = 'MAXI' lmlp ;
'SI' (milmlp '<' 2.57) ; lok = lok 'ET' faux ; 'FINS' ;
'SI' (malmlp '>' 2.61) ; lok = lok 'ET' faux ; 'FINS' ;
 
************************************************************************
*                     E N   D I M E N S I O N   3                      *
************************************************************************
 
OPTI 'DIME' 3 ;
p1 = 0. 0. 0. ;
p2 = 0. 0. 1. ;
p3 = 0. 1. 1. ;
p4 = 0.5 0. 1. ;
p5 = 0.5 0. 2. ;
mm = @REPERE ('PROG' 0.25 0.25 0.25) 'TURQ' ;
 
* Éléments de maillage SEG2
OPTI 'ELEM' 'SEG2' ;
l1 = 'DROI' p1 p2 ;
l2 = 'DROI' p2 p3 ;
l3 = 'DROI' p4 p5 ;
*
ll1   = l1 'ET' l2 ;
*'OPTI' impi 1 ;
v1 = 'MOIN' p2 p1 ; v2 = 'MOIN' p3 p2 ;
lct1  = 'CUBT' nl p1 v1 v2 p3 'COUL' 'OLIV' ;
lct1u  = 'CUBT' nl p1 v1 v2 p3 'UNIF' 'COUL' 'ROUG' ;
*
'SI' graph ;
   TRAC 'QUAL'  (ll1 'ET' lct1 'ET' lct1u 'ET' mm) 'TITR'
   'Lignes et cubiques en angle SEG2' ;
'FINS' ;
ll2 = l1 'ET' l3 ;
lcp1  = 'CUBP' nl p1 p2 p4 p5 'COUL' 'OLIV' ;
lcp1u  = 'CUBP' nl p1 p2 p4 p5 'UNIF' 'COUL' 'ROUG' ;
'SI' graph ;
   TRAC 'QUAL'  (ll2 'ET' lcp1 'ET' lcp1u 'ET' mm) 'TITR'
   'Lignes et cubiques en parallele SEG2' ;
'FINS' ;
 
* Éléments de maillage SEG3
OPTI 'ELEM' 'SEG3' ;
lct2  = 'CUBT' nl p1 v1 v2 p3 'COUL' 'OLIV' ;
lct2u  = 'CUBT' nl p1 v1 v2 p3 'UNIF' 'COUL' 'ROUG' ;
*
'SI' graph ;
   TRAC 'QUAL'  (ll1 'ET' lct2 'ET' lct2u 'ET' mm) 'TITR'
   'Lignes et cubiques en angle SEG3' ;
'FINS' ;
lcp2  = 'CUBP' nl p1 p2 p4 p5 'COUL' 'OLIV' ;
lcp2u  = 'CUBP' nl p1 p2 p4 p5 'UNIF' 'COUL' 'ROUG' ;
'SI' graph ;
   TRAC 'QUAL'  (ll2 'ET' lcp2 'ET' lcp2u 'ET' mm) 'TITR'
   'Lignes et cubiques en parallele SEG3' ;
'FINS' ;
*
mlct1 = 'MESU' lct1 ;
mlct1u = 'MESU' lct1u ;
mlct2 = 'MESU' lct2 ;
mlct2u = 'MESU' lct2u ;
lmlt = 'PROG' mlct1 mlct1u mlct2 mlct2u ;
milmlt = 'MINI' lmlt ;
malmlt = 'MAXI' lmlt ;
'SI' (milmlt '<' 1.48) ; lok = lok 'ET' faux ; 'FINS' ;
'SI' (malmlt '>' 1.49) ; lok = lok 'ET' faux ; 'FINS' ;
*
mlcp1 = 'MESU' lcp1 ;
mlcp1u = 'MESU' lcp1u ;
mlcp2 = 'MESU' lcp2 ;
mlcp2u = 'MESU' lcp2u ;
lmlp = 'PROG' mlcp1 mlcp1u mlcp2 mlcp2u ;
milmlp = 'MINI' lmlp ;
malmlp = 'MAXI' lmlp ;
'SI' (milmlp '<' 2.57) ; lok = lok 'ET' faux ; 'FINS' ;
'SI' (malmlp '>' 2.61) ; lok = lok 'ET' faux ; 'FINS' ;
*
'SI' ('NON' lok) ;
   'MESS' 'Il y a eu des erreurs' ;
'SINO' ;
   'MESS' 'Tout sest bien passe' ;
'FINS' ;
'SI' interact ;
    'OPTI' 'DONN'  5 ;
'FINS' ;
'SI' ('NON' lok) ;
   'ERRE' 5 ;
'FINS' ;
*
* End of dgibi file CUBP-CUBT
*
'FIN' ;