Télécharger sigjoi.eso

Retour à la liste

Numérotation des lignes :

sigjoi
  1. C SIGJOI SOURCE AM 11/05/26 21:16:41 6982
  2. C SIGLIA SOURCE CHAT 05/01/13 03:16:39 5004
  3. C SIGLIA SOURCE DJER 23/01/96
  4. SUBROUTINE SIGJOI(NMATT,XDDL,WORK,LRE,XSTRS,IDIM,NSTRS,CMATE)
  5. C====================================================================
  6. C CALCULE DES CONTRAINTES POUR JOI1 :
  7. C element de liaison a 2 noeuds (6 ddl par noeuds)
  8. C (ON CALCUL en fait DES EFFORTs NODAUX )
  9. C====================================================================
  10. IMPLICIT INTEGER(I-N)
  11. IMPLICIT REAL*8(A-H,O-Z)
  12. CHARACTER*8 CMATE
  13. DIMENSION XDDL(LRE),WORK(NMATT),XSTRS(*),DEPL(12)
  14. DIMENSION HK3(6,6),HK2(3,3)
  15. *-------------------------------------------------------------
  16. DO I =1,LRE
  17. DEPL(I) = XDDL(I)
  18. END DO
  19. C
  20. IF (IDIM.EQ.3) THEN
  21. C CAS D'UN MATERIAU ORTHOTROPE
  22. C CALCUL DES EFFORTS NODAUX
  23. IF (CMATE.EQ.'ORTHOTRO') THEN
  24. DO I =1,NSTRS
  25. XSTRS(I)=WORK(I+6)*(DEPL(I+NSTRS)-DEPL(I))
  26. END DO
  27. C CAS D'UN MATERIAU ANISOTROPE
  28. ELSE IF (CMATE.EQ.'ANISOTRO') THEN
  29. C INSERTION DES RAIDEURS MATERIELLES DANS UNE MATRICE 6*6
  30. K = 1
  31. DO I=1,NSTRS
  32. HK3(I,I)=WORK(I+6)
  33. DO J=1,I-1
  34. HK3(I,J)=WORK(K+12)
  35. HK3(J,I)=HK3(I,J)
  36. K = K+1
  37. ENDDO
  38. ENDDO
  39. C CALCUL DES EFFORTS NODAUX
  40. DO I=1,NSTRS
  41. DO J=1,NSTRS
  42. XSTRS(I)=XSTRS(I)+HK3(I,J)*(DEPL(J+NSTRS)-DEPL(J))
  43. ENDDO
  44. ENDDO
  45. ENDIF
  46. C
  47. *
  48. C
  49. ELSE IF (IDIM.EQ.2) THEN
  50. C CAS D'UN MATERIAU ORTHOTROPE
  51. C CALCUL DES EFFORTS NODAUX
  52. IF (CMATE.EQ.'ORTHOTRO') THEN
  53. DO I =1,NSTRS
  54. XSTRS(I)=WORK(I+2)*(DEPL(I+NSTRS)-DEPL(I))
  55. END DO
  56. C CAS D'UN MATERIAU ANISOTROPE
  57. ELSE IF (CMATE.EQ.'ANISOTRO') THEN
  58. C INSERTION DES RAIDEURS MATERIELLES DANS UNE MATRICE 3*3
  59. K = 1
  60. DO I=1,NSTRS
  61. HK2(I,I)=WORK(I+2)
  62. DO J=1,I-1
  63. HK2(I,J)=WORK(K+5)
  64. HK2(J,I)=HK2(I,J)
  65. K = K+1
  66. ENDDO
  67. ENDDO
  68. C CALCUL DES EFFORTS NODAUX
  69. DO I=1,NSTRS
  70. DO J=1,NSTRS
  71. XSTRS(I)=XSTRS(I)+HK2(I,J)*(DEPL(J+NSTRS)-DEPL(J))
  72. ENDDO
  73. ENDDO
  74. ENDIF
  75. ENDIF
  76. *
  77. *
  78. *
  79. *-------------------------------------------------------------
  80. * WRITE(*,*) '** EFFORT ( dans siglia) ** '
  81. * WRITE(*,1992) (XSTRS(IC),IC=1,12)
  82. *1992 FORMAT(6(6(1X,E12.5),/))
  83. *-------------------------------------------------------------
  84. RETURN
  85. END
  86.  
  87.  
  88.  
  89.  
  90.  
  91.  

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