Télécharger la plaquette Cast3M


fr en

Télécharger dogo88.eso

Retour à la liste

Numérotation des lignes : 

dogo88

  1. C DOGO88    SOURCE    CHAT      05/01/12    22:53:12     5004
  2.       SUBROUTINE DOGO88(VELA,MATE,EPAIST,IFOU,LHOOK,DDHOOK,IRET)
  3. C
  4. C=======================================================================
  5. C
  6. C  MATRICE DE HOOKE DES ELEMENTS JOINTS GENERALISE JGI2 (2D),
  7. C                                          JGI4 ET JGT3 (3D)
  8. C
  9. C  ENTREES
  10. C     VELA() = materiau dans un tableau de travail
  11. C     MATE   = Nom du materiau
  12. C     IFOU   = num{ro d'harmonique de  fourier: IFOUR de CCOPTIO
  13. C     LHOOK  = taille de la matrice de hooke
  14. C
  15. C  SORTIES
  16. C     DDHOOK(LHOOK,LHOOK) = matrice de hooke
  17. C     IRET = 1 si option existante 0 SINON
  18. C
  19. C  REMARQUE : ATTENTION : DANS LES CAS  CONTRAINTES PLANES,
  20. C             DEFORMATIONS PLANES ET AXISYMETRIQUE,
  21. C             LA MATRICE DE HOOK N'A PLUS LA DIMENSION
  22. C             (3X3), MAIS (2X2). C'EST DONC LES DIMENSIONS DU CAS
  23. C             BIDIMENSIONNEL.
  24. C
  25. C  CECI EST UN DOUBLE DE DOHO88, MAIS AVEC LA THEORIE DE PARVIS
  26. C  POUR LE CALCUL DE LA MATRICE DE TRANSFORMATION. LA MATRICE DE
  27. C  ROTATION ROTHOO TRANSFORME UN VECTEUR (EN L'OCCURENCE LE VECTEUR
  28. C  DEPLACEMENT) DU REPERE D'ORTHOTROPIE AU REPERE GLOBAL. IL FAUT
  29. C  BIEN VOIR QUE CECI EST UNE TRANSFORMATION D'UN VECTEUR EN UN AUTRE
  30. C  VECTEUR, ET NON PAS LA TRANSFORMATION D'UN TENSEUR (EX: TENSEUR
  31. C  DES CONTRAINTES) EN UN AUTRE TENSEUR. C'EST DONC CETTE SUBROUTINE
  32. C  QUI EST LA BONNE, CAR JE RAISONNE EN TERMES DE DEPLACEMENTS, ET
  33. C  NON EN TERMES DE CONTRAINTES, DANS LE CAS DES JOINTS (cf. LA
  34. C  RELATION FORCES/DEPLACEMENTS ET NON PAS CONTRAINTES/DEFORMATIONS)
  35. C
  36. C          P. PEGON 2001 d'apres DOUO88
  37. C=======================================================================
  38. C
  39.       IMPLICIT INTEGER(I-N)
  40.       IMPLICIT REAL*8(A-H,O-Z)
  41.       CHARACTER*8 MATE
  42. C
  43.       DIMENSION VELA(*),DDHOOK(LHOOK,*)
  44.       REAL*8 D3HOO1(3,3),D3HOO2(3,3),ROTHOO(3,3)
  45.       DATA UN,DEUX,UNDEMI/1.D0,2.0D0,.5D0/
  46. C
  47.       IRET=1
  48.       CALL ZERO(DDHOOK,LHOOK,LHOOK)
  49. C---->attention, LHOOK vaut 6...
  50. C     CALL ZERO(D3HOO1,LHOOK,LHOOK)
  51. C
  52. C     MATERIAU ISOTROPE
  53. C
  54.       IF(MATE.EQ.'ISOTROPE') THEN
  55. C
  56.         YOU  =VELA(1)
  57.         XNU  =VELA(2)
  58. C
  59. C     CAS CONTRAINTES PLANES DANS LE PLAN (S1,N)
  60. C
  61.         IF(IFOU.EQ.-2) THEN
  62.           AUX =YOU/(UN-XNU*XNU)
  63.           AUX1=AUX*XNU
  64. *
  65.           DDHOOK(1,1)=(UN-XNU)*AUX*UNDEMI/EPAIST
  66. *
  67.           DDHOOK(2,2)=   AUX  / EPAIST
  68.           DDHOOK(2,3)=   AUX1
  69. *
  70.           DDHOOK(3,2)=   AUX1
  71.           DDHOOK(3,3)=   AUX  * EPAIST
  72. C
  73. C      CAS DEFORMATIONS PLANES DANS LE PLAN (S1,N)
  74. C
  75.         ELSE IF(IFOU.EQ.-3.OR.IFOU.EQ.-1) THEN
  76.           AUX0=UN/((UN+XNU)*(UN-DEUX*XNU))
  77.           AUX =YOU*AUX0*(UN-XNU)
  78.           AUX1=YOU*AUX0*XNU
  79. *
  80.           DDHOOK(1,1)=YOU*UNDEMI/(UN+XNU)/EPAIST
  81. *
  82.           DDHOOK(2,2)=   AUX  / EPAIST
  83.           DDHOOK(2,3)=   AUX1
  84.           DDHOOK(2,4)=   AUX1
  85. *
  86.           DDHOOK(3,2)=   AUX1
  87.           DDHOOK(3,3)=   AUX  * EPAIST
  88.           DDHOOK(3,4)=   AUX1 * EPAIST
  89. *
  90.           DDHOOK(4,2)=   AUX1
  91.           DDHOOK(4,3)=   AUX1 * EPAIST
  92.           DDHOOK(4,4)=   AUX  * EPAIST
  93. C
  94. C     CAS TRIDIMENSIONNEL
  95. C
  96.         ELSE IF(IFOU.EQ.2) THEN
  97.           AUX0=UN/((UN+XNU)*(UN-DEUX*XNU))
  98.           AUX =YOU*AUX0*(UN-XNU)
  99.           AUX1=YOU*AUX0*XNU
  100. *
  101.           DDHOOK(1,1)=YOU*UNDEMI/(UN+XNU)/EPAIST
  102. *
  103.           DDHOOK(2,2)=DDHOOK(1,1)
  104. *
  105.           DDHOOK(3,3)=   AUX  / EPAIST
  106.           DDHOOK(3,4)=   AUX1
  107.           DDHOOK(3,5)=   AUX1
  108. *
  109.           DDHOOK(4,3)=   AUX1
  110.           DDHOOK(4,4)=   AUX  * EPAIST
  111.           DDHOOK(4,5)=   AUX1 * EPAIST
  112. *
  113.           DDHOOK(5,3)=   AUX1
  114.           DDHOOK(5,4)=   AUX1 * EPAIST
  115.           DDHOOK(5,5)=   AUX  * EPAIST
  116. *
  117.           DDHOOK(6,6)=YOU*UNDEMI/(UN+XNU)*EPAIST
  118. C
  119. C     OPTION NON DEFINIE
  120. C
  121.         ELSE
  122.           IRET=0
  123.         ENDIF
  124. C
  125. C     MATERIAU ORTHOTROPE (Pour le moment rien!)
  126. C
  127. C     ELSE IF(MATE.EQ.'ORTHOTRO') THEN
  128. C
  129. C     CAS TRIDIMENSIONNEL
  130. C
  131. C       IF(IFOU.EQ.2) THEN
  132. C         ZK1 =VELA(1)
  133. C         ZK2 =VELA(2)
  134. C         ZK3 =VELA(3)
  135. C         COSA=VELA(4)
  136. C         SINA=VELA(5)
  137. C
  138. C         D3HOO1(1,1)=ZK1
  139. C         D3HOO1(2,2)=ZK2
  140. C         D3HOO1(3,3)=ZK3
  141. C
  142. C         ROTHOO(1,1) =  COSA
  143. C         ROTHOO(1,2) = -SINA
  144. C         ROTHOO(1,3) =  0.0D0
  145. C         ROTHOO(2,1) =  SINA
  146. C         ROTHOO(2,2) =  COSA
  147. C         ROTHOO(2,3) =  0.0D0
  148. C         ROTHOO(3,1) =  0.0D0
  149. C         ROTHOO(3,2) =  0.0D0
  150. C         ROTHOO(3,3) =  1.0D0
  151. C
  152. C       PASSAGE DANS LE REPERE DE L'ELEMENT:
  153. C
  154. C         CALL PRODT (D3HOO2,D3HOO1,ROTHOO,3,3)
  155. C         DO 400 J=1,3
  156. C           DO 400 I=1,3
  157. C              DDHOOK(I,J) = D3HOO2(I,J)
  158. C 400     CONTINUE
  159. C
  160. C           OPTION NON DEFINIE
  161. C
  162. C         ELSE
  163. C           IRET=0
  164. C         ENDIF
  165. C
  166. C     OPTION NON DEFINIE
  167. C
  168.       ELSE
  169.         IRET=0
  170.       ENDIF
  171. C
  172.       RETURN
  173.       END
  174.  
  175.  
  176.  
  177.  
  178.  

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