Télécharger la plaquette Cast3M


fr en

Télécharger deadjaco.procedur

Retour à la liste

Numérotation des lignes : 

  1. * DEADJACO  PROCEDUR  GOUNAND   26/01/12    21:15:02     12448          
  2. ************************************************************************
  3. * NOM         : DEADJACO
  4. * DESCRIPTION : Le déterminant de la matrice jacobienne
  5. *               de la transformation géométrique (par élément)
  6. *
  7. *
  8. *
  9. * LANGAGE     : GIBIANE-CAST3M
  10. * AUTEUR      : Stéphane GOUNAND (CEA/DEN/DM2S/SFME/LTMF)
  11. *               mél : gounand@semt2.smts.cea.fr
  12. **********************************************************************
  13. * VERSION    : v1, 05/04/2006, version initiale
  14. * HISTORIQUE : v1, 05/04/2006, création
  15. * HISTORIQUE : 2025/12/29 : simplification donnees NLIN
  16. * HISTORIQUE :
  17. * HISTORIQUE :
  18. ************************************************************************
  19. * Prière de PRENDRE LE TEMPS de compléter les commentaires
  20. * en cas de modification de ce sous-programme afin de faciliter
  21. * la maintenance !
  22. ************************************************************************
  23. *
  24. *
  25. 'DEBPROC' DEADJACO ;
  26. 'ARGUMENT' _mt*'MAILLAGE' ;
  27. 'ARGUMENT' moterrj/'MOT' ;
  28.  
  29. lerrj = 'EXIS' moterrj ;
  30. 'SI' lerrj ;
  31.    'SI' ('NEG' moterrj 'ERRJ') ;
  32. *1052 2
  33. *Mot-cle incorrect "%M1:4". Voici la liste des valeurs admises : %M5:40
  34.       'ERRE' 1052 'AVEC' moterrj 'ERRJ' ;
  35.    'FINS' ;
  36. 'FINS' ;
  37.  
  38. gdisc = DEADUTIL 'TYPM' _mt ;
  39. 'SI' ('EGA' gdisc 'LINE') ;
  40.    methgau = 'GAM1' ;
  41. 'SINON' ;
  42.    methgau = 'GAM2' ;
  43. 'FINSI' ;
  44. *methgau = 'GAU7' ;
  45. *
  46. * Calcul du résidu
  47. *
  48. vdim  = DEADUTIL 'DIMM' _mt ;
  49. numop  = 1 ;
  50. numvar = 1 ;
  51. numder = vdim ;
  52. numdat = 0 ;
  53. numcof = 1 ;
  54. *
  55.    A = ININLIN numop numvar numdat numcof numder ;
  56.    A . 'VAR' . 1 . 'VALEUR' = 1.D0 ;
  57.    A . 'COF' . 1 . 'COMPOR' = 'CHAINE' 'VOLORI' ;
  58. *
  59.    A . 1 . 1 . 0  = 1 ;
  60. *
  61.    numcof = 0 ;
  62.    B = ININLIN numop numvar numdat numcof numder ;
  63. * Le Jacobien est nul ou presque sur les noeuds sommet du TRI6
  64. * donc ce serait contraignant
  65. *   B . 'VAR' . 1 . 'NOMDDL' = 'SCAL' ;
  66. *   B . 'VAR' . 1 . 'DISC'   = gdisc ;
  67.    B . 'VAR' . 1 . 'VALEUR' = 1.D0 ;
  68.    B . 1 . 1 . 0  = 0 ;
  69. *
  70. 'SI' lerrj ;
  71.    detj = 'NLIN' gdisc _mt A B 'EREF' 'ERRJ' 'CHAM' methgau ;
  72. 'SINO' ;
  73.    detj = 'NLIN' gdisc _mt A B 'EREF' 'CHAM' methgau ;
  74. 'FINS' ;
  75. *
  76. 'RESPRO' detj ;
  77. *
  78. * End of procedure file DEADJACO
  79. *
  80. 'FINPROC' ;
  81.  
  82.  

© Cast3M 2003 - Tous droits réservés.
Mentions légales