Télécharger la plaquette Cast3M


fr en

Télécharger relmod.eso

Retour à la liste

Numérotation des lignes : 

relmod

  1. C RELMOD    SOURCE    CB215821  24/04/12    21:17:06     11897          
  2.       SUBROUTINE RELMOD(IMODEL) 
  3. C***********************************************************************
  4. C   cet operateur impose des relation linéaire entre les inconnes 
  5. C   d'un modele pour le cas particulier des modelel SURE et XFEM
  6. C
  7. C   syntaxe :
  8. C    rig1 = rela   mod1
  9. C
  10. C    mod1 modele pouvant eventuellement contenit des sous 
  11. C           modeles avec des relations linéaires à imposer : 
  12. C
  13. C           MODELE de SURE : 
  14. C                           elements géométrique  itypel= 48
  15. C                           elements du modeles   nefmod = 259
  16. C                           formulation           mfr=1
  17. C
  18. C           MODELE de XFEM (XQ4R):
  19. C                           elements géométrique  itypel= 8
  20. C                           elements du modeles   nefmod = 263
  21. C                           formulation           mfr=63   
  22. C           MODELE de XFEM (XC8R):
  23. C                           elements géométrique  itypel= 14
  24. C                           elements du modeles   nefmod = 264
  25. C                           formulation           mfr=63
  26. C
  27. C                                                                     
  28. C   ENTREES :                                                         
  29. C   ________                                                          
  30. C                                                                     
  31. C        IMODEL   pointeur sur le modele                              
  32. C   ENTREES/SORTIE :                                                  
  33. C   ________                                                          
  34. C                                                                     
  35. C        Rigidité sortie via ecrobj                                   
  36. C---------------------------------------------------------------------   
  37. C***********************************************************************
  38.  
  39.       IMPLICIT INTEGER(I-N)
  40.       IMPLICIT REAL*8 (A-H,O-Z)
  41.  
  42. -INC PPARAM
  43. -INC CCOPTIO
  44. -INC CCGEOME
  45. -INC SMELEME
  46. -INC SMCOORD
  47. -INC SMRIGID
  48. -INC SMMODEL
  49. -INC CCHAMP
  50.  
  51.       MMODEL= IMODEL
  52.       SEGACT MMODEL
  53.  
  54.       NSOUS = MMODEL.KMODEL(/1)
  55. C--------------------------------------------------------------------*
  56. C       initialisation de la rigidité a sortir
  57. C--------------------------------------------------------------------*
  58.       NRIGEL= 0
  59.       SEGINI RI1
  60.       RI1.IFORIG = IFOUR
  61.       RI1.MTYMAT = 'RIGIDITE'
  62.  
  63. C--------------------------------------------------------------------*
  64. C
  65. C       BOUCLE SUR LES ZONES ELEMENTAIRES ( MEME TYPE D'EF )
  66. C
  67. C--------------------------------------------------------------------*
  68.  
  69.       DO 1 ISOUS=1,NSOUS
  70.             IMODE1=MMODEL.KMODEL(ISOUS)
  71.             SEGACT IMODE1
  72. C
  73. C    INITIALISATIONS
  74. C
  75.             MELE   = IMODE1.NEFMOD        
  76.             IPT1 = IMODE1.IMAMOD
  77.             SEGACT IPT1
  78.             ITYP1 = IPT1.ITYPEL
  79. c            write(*,*) 'MELE', MELE
  80. c            write(*,*) 'ITYP1', ITYP1
  81.  
  82.             IF (MELE.eq.263.or.MELE.eq.264) then
  83.  
  84. c          ajout d'une matrice élémentaire de rigidité pour les ddl xfem
  85. c          a désactiver.
  86.  
  87.               CALL RELXFE(IMODE1 , RI1 )
  88.             ENDIF
  89.             IF (MELE.eq.259) then
  90.  
  91. c          ajout d'une matrice élémentaire de rigidité pour chaque 
  92. c          composantes du modele de SURE
  93.  
  94.               CALL RELSUR(IMODE1 , RI1 )
  95.             ENDIF
  96.             SEGSUP IMODE1
  97.   1   CONTINUE
  98.  
  99.       CALL ECROBJ('RIGIDITE', RI1 )
  100. c      SEGSUP RI1  
  101.  
  102.       END
  103.  
  104.  
  105.  
  106.  

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