Télécharger rigixr.eso

Retour à la liste

Numérotation des lignes :

  1. C RIGIXR SOURCE GG250959 17/09/20 21:16:44 9554
  2. subroutine RIGIXR(ISOU,IPOI6,IMODEL,IPINF,
  3. $ IVAMAT,IVACAR,NMATT,CMATE,NCAR1,NBGMAT,NELMAT,IMAT,IRETER)
  4. * as 2009/09/03 : ajout de IMAT en entrée de RIGIXR
  5. c
  6. C Les sous programmes affectés à un type d'élément sont chargés
  7. C de faire le tri des éléments suivant le type d'enrichissement.
  8. c
  9. C Par exemple le XQ4R peut sortir :
  10. c - des matrices ne comportant que des ddl standard du QUA4,
  11. c - des matrices enrichies par le saut du à la fissure,
  12. c - des matrices enrichies par le saut et par les fonctions de
  13. c la mécanique de la rupture
  14. c
  15. C Donc pour 1 type d'EF, on crée 1+NENR objet elementaire IMATTT
  16. C Pour y parvenir, on utilise le tableau LOCIRI (=LOCal IRIgel)
  17. c qui est l'équivalent local de IRIGEL et qui doit etre rempli
  18. C par les sous programmes élémentaires (comme RIGIX.eso).
  19. c Il est dimensionné en dur à (10,6) au lieu de (8,1+NENR)
  20. c pour etre large.
  21. C
  22. C Finalement, RIGIXR.eso recupere LOCIRI, ajuste et remplit IRIGEL
  23. c
  24. IMPLICIT REAL*8 (A-H,O-Z)
  25. IMPLICIT INTEGER (I-N)
  26. C
  27. -INC CCOPTIO
  28. -INC SMRIGID
  29. -INC SMINTE
  30. -INC SMMODEL
  31. -INC SMELEME
  32. C
  33. CHARACTER*8 CMATE
  34. PARAMETER (NBENRMAX=30)
  35. INTEGER LOCIRI(10,(1+NBENRMAX))
  36. c
  37. SEGMENT MPTVAL
  38. INTEGER IPOS(NS),NSOF(NS)
  39. INTEGER IVAL(NCOSOU)
  40. CHARACTER*16 TYVAL(NCOSOU)
  41. ENDSEGMENT
  42. c
  43. c write(6,*) '##### entree dans rigixr #####'
  44. IRETER=0
  45. C
  46. Ccccc on active le modele, les caracteristiques materiau
  47. c* segact,imodel
  48.  
  49. mele = nefmod
  50. C element XQ4R (2D rupture) ou element XC8R (3D rupture)
  51. IF (mele.NE.263 .and. mele.NE.264) then
  52. call erreur (21)
  53. RETURN
  54. ENDIF
  55.  
  56. c* mptval=ivacar
  57. c* segact mptval
  58.  
  59. Ccccc on initialise LOCIRI
  60. do i1=1,10
  61. do i2=1,(1+NBENRMAX)
  62. LOCIRI(i1,i2) = 0
  63. enddo
  64. enddo
  65. C
  66. ccccc Appel au calcul des rigidites elementaires
  67. call rigix(ivamat,ivacar,NMATT,CMATE,
  68. & imodel,IPINF,LOCIRI,NBGMAT,NELMAT,IMAT)
  69. * as 2009/09/03 : ajout de IMAT en entrée de RIGIX
  70. ccccc de combien faut-il augmenter INFELE ?
  71. c write(6,*) 'LOCIRI=',(LOCIRI(1,iou),iou=1,5)
  72. nrigsup = -1
  73. do i=1,(1+NBENRMAX)
  74. if (LOCIRI(1,i).ne.0) nrigsup=nrigsup+1
  75. enddo
  76. if (nrigsup.eq.-1) then
  77. write(*,*) 'on n a meme pas reussi a construire des rigidite',
  78. & 'associes aux ddl std'
  79. CALL ERREUR(21)
  80. return
  81. endif
  82.  
  83. c write(*,*) 'ccccc on ouvre en modification MRIGID'
  84. MRIGID=IPOI6
  85. segact,MRIGID*mod
  86.  
  87. NRIGE = IRIGEL(/1)
  88. nrigini = IRIGEL(/2)
  89. NRIGEL = nrigini + nrigsup
  90. if (nrigsup.gt.0) then
  91. C write(ioimp,*) 'c modification de la taille de MRIGID de RIGI'
  92. C & ,isou,'->',isou+nrigsup
  93. c write(ioimp,*) 'on doit augmenter IRIGEL de',nrigsup,
  94. c & ' rigidites elementaires'
  95. segadj,MRIGID
  96. ENDIF
  97. c write(*,*) 'ccccc remplissage de MRIGID'
  98.  
  99. c Cas ou il ya une partie std -----------------------
  100. IF (LOCIRI(1,1).ne.0) THEN
  101.  
  102. c + partie non enrichie (=std)
  103. do i1 = 1, NRIGE
  104. IRIGEL(i1,isou)=LOCIRI(i1,1)
  105. enddo
  106. COERIG(isou)= 1.D0
  107.  
  108. c + partie enrichie
  109. if (nrigsup.gt.0) then
  110. ia = isou
  111. iloc = 1
  112. do i = 1, nrigsup
  113. ia = ia + 1
  114. c petit ajout pour le cas ou on a "sauté" le H-enrichissement
  115. 1000 continue
  116. iloc = iloc + 1
  117. if( (LOCIRI(1,iloc)) .eq. 0) goto 1000
  118. do i1 = 1,NRIGE
  119. IRIGEL(i1,ia) = LOCIRI(i1,iloc)
  120. enddo
  121. COERIG(ia)= 1.D0
  122. enddo
  123. isou = isou + nrigsup
  124. endif
  125.  
  126. c Cas ou il n y a pas (plus) de partie std -----------------------
  127. ELSE
  128.  
  129. c + partie enrichie
  130. ia = isou - 1
  131. iloc = 1
  132. do i=0, nrigsup
  133. ia = ia + 1
  134. c petit ajout pour le cas ou on a "sauté" le H-enrichissement
  135. 2000 continue
  136. iloc = iloc + 1
  137. if( (LOCIRI(1,iloc)) .eq. 0) goto 2000
  138. do i1=1,NRIGE
  139. IRIGEL(i1,ia) = LOCIRI(i1,iloc)
  140. enddo
  141. COERIG(ia)= 1.D0
  142. enddo
  143. isou = isou + nrigsup
  144.  
  145. ENDIF
  146. *
  147. c write(ioimp,*) 'IRIGEL=',(IRIGEL(1,iou),iou=1,NRIGEL)
  148. c
  149. return
  150. end
  151.  
  152.  
  153.  
  154.  

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