Télécharger rigjoi.eso

Retour à la liste

Numérotation des lignes :

  1. C RIGJOI SOURCE AM 11/05/26 21:16:38 6982
  2. C RIGLI2 SOURCE CHAT 05/01/13 03:00:02 5004
  3. C RIGLI2 SOURCE DJER 94/12/19 21:15:07 654
  4. SUBROUTINE RIGJOI(NMATT,REL,LRE,WORK,IDIM,CMATE)
  5. C====================================================================
  6. C
  7. C MODIFIE REL (LRE,LRE) MATRICE DE RIGIDITE EN INTRODUISANT DES
  8. C RAIDEUR LOCALES EQUIVALENTES
  9. C====================================================================
  10. IMPLICIT INTEGER(I-N)
  11. IMPLICIT REAL*8(A-H,O-Z)
  12. CHARACTER*8 CMATE
  13. DIMENSION REL(LRE,LRE),WORK(NMATT)
  14. DIMENSION RE2(3,3),RE3(6,6)
  15. C
  16. CALL ZERO(REL,LRE,LRE)
  17. XZER = 0.D0
  18. C
  19. IF (IDIM.EQ.3) THEN
  20. CALL ZERO(RE3,6,6)
  21. DO I =1,6
  22. RE3(I,I)=WORK(I+6)
  23. END DO
  24. * SI anisotrope
  25. IF (CMATE.EQ.'ANISOTRO') THEN
  26. K=0
  27. DO I =2,6
  28. DO J =1,I-1
  29. K=K+1
  30. RE3(I,J)=WORK(K+12)
  31. RE3(J,I)=RE3(I,J)
  32. END DO
  33. END DO
  34. ENDIF
  35. C
  36. DO I =1,6
  37. DO J =1,6
  38. REL(I,J)=RE3(I,J)
  39. REL(I+6,J)=-1*RE3(I,J)
  40. REL(I,J+6)=-1*RE3(I,J)
  41. REL(I+6,J+6)=RE3(I,J)
  42. END DO
  43. END DO
  44. *
  45. *
  46. *
  47. ELSE IF (IDIM.EQ.2) THEN
  48. CALL ZERO(RE2,3,3)
  49. DO I =1,3
  50. RE2(I,I)=WORK(I+2)
  51. END DO
  52. * SI anisotrope
  53. IF (CMATE.EQ.'ANISOTRO') THEN
  54. K=0
  55. DO I =2,3
  56. DO J =1,I-1
  57. K=K+1
  58. RE2(I,J)=WORK(K+5)
  59. RE2(J,I)=RE2(I,J)
  60. END DO
  61. END DO
  62. ENDIF
  63. C
  64. DO I =1,3
  65. DO J =1,3
  66. REL(I,J)=RE2(I,J)
  67. REL(I,J+3)=-1*RE2(I,J)
  68. REL(I+3,J)=-1*RE2(I,J)
  69. REL(I+3,J+3)=RE2(I,J)
  70. END DO
  71. END DO
  72. ENDIF
  73. 5 CONTINUE
  74. *
  75. *-------------------------------------------------------------
  76. * WRITE(*,*) '** REL APRES ( dans RIGJOI) ** '
  77. * DO 112 IA=1,LRE
  78. * WRITE(*,1991) (REL(IA,IC),IC=1,LRE)
  79. * 112 CONTINUE
  80. *1991 FORMAT(18(1X,E12.5))
  81. *-------------------------------------------------------------
  82. RETURN
  83. END
  84.  
  85.  
  86.  
  87.  
  88.  
  89.  

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