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.  
  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.  
  61. C--------------------------------------------------------------------*
  62. C
  63. C BOUCLE SUR LES ZONES ELEMENTAIRES ( MEME TYPE D'EF )
  64. C
  65. C--------------------------------------------------------------------*
  66.  
  67. DO 1 ISOUS=1,NSOUS
  68. IMODE1=MMODEL.KMODEL(ISOUS)
  69. SEGACT IMODE1
  70. C
  71. C INITIALISATIONS
  72. C
  73. MELE = IMODE1.NEFMOD
  74. IPT1 = IMODE1.IMAMOD
  75. SEGACT IPT1
  76. ITYP1 = IPT1.ITYPEL
  77. c write(*,*) 'MELE', MELE
  78. c write(*,*) 'ITYP1', ITYP1
  79.  
  80. IF (MELE.eq.263.or.MELE.eq.264) then
  81.  
  82. c ajout d'une matrice élémentaire de rigidité pour les ddl xfem
  83. c a désactiver.
  84.  
  85. CALL RELXFE(IMODE1 , RI1 )
  86. ENDIF
  87. IF (MELE.eq.259) then
  88.  
  89. c ajout d'une matrice élémentaire de rigidité pour chaque
  90. c composantes du modele de SURE
  91.  
  92. CALL RELSUR(IMODE1 , RI1 )
  93. ENDIF
  94. SEGSUP IMODE1
  95. 1 CONTINUE
  96.  
  97. CALL ECROBJ('RIGIDITE', RI1 )
  98. c SEGSUP RI1
  99.  
  100. END
  101.  
  102.  
  103.  

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