Télécharger relmod.eso

Retour à la liste

Numérotation des lignes :

  1. C RELMOD SOURCE CB215821 19/08/20 21:21:34 10287
  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. -INC CCOPTIO
  42. -INC CCGEOME
  43. -INC SMELEME
  44. -INC SMCOORD
  45. -INC SMRIGID
  46. -INC SMMODEL
  47. -INC CCHAMP
  48.  
  49. MMODEL= IMODEL
  50. SEGACT MMODEL
  51.  
  52. NSOUS = MMODEL.KMODEL(/1)
  53. C--------------------------------------------------------------------*
  54. C initialisation de la rigidité a sortir
  55. C--------------------------------------------------------------------*
  56. NRIGEL= 0
  57. SEGINI RI1
  58.  
  59. C--------------------------------------------------------------------*
  60. C
  61. C BOUCLE SUR LES ZONES ELEMENTAIRES ( MEME TYPE D'EF )
  62. C
  63. C--------------------------------------------------------------------*
  64.  
  65. DO 1 ISOUS=1,NSOUS
  66. IMODE1=MMODEL.KMODEL(ISOUS)
  67. SEGACT IMODE1
  68. C
  69. C INITIALISATIONS
  70. C
  71. MELE = IMODE1.NEFMOD
  72. IPT1 = IMODE1.IMAMOD
  73. SEGACT IPT1
  74. ITYP1 = IPT1.ITYPEL
  75. c write(*,*) 'MELE', MELE
  76. c write(*,*) 'ITYP1', ITYP1
  77.  
  78. IF (MELE.eq.263.or.MELE.eq.264) then
  79.  
  80. c ajout d'une matrice élémentaire de rigidité pour les ddl xfem
  81. c a désactiver.
  82.  
  83. CALL RELXFE(IMODE1 , RI1 )
  84. ENDIF
  85. IF (MELE.eq.259) then
  86.  
  87. c ajout d'une matrice élémentaire de rigidité pour chaque
  88. c composantes du modele de SURE
  89.  
  90. CALL RELSUR(IMODE1 , RI1 )
  91. ENDIF
  92. SEGSUP IMODE1
  93. 1 CONTINUE
  94.  
  95. CALL ECROBJ('RIGIDITE', RI1 )
  96. c SEGSUP RI1
  97.  
  98. END
  99.  
  100.  
  101.  

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