Télécharger maxrot.eso

Retour à la liste

Numérotation des lignes :

  1. C MAXROT SOURCE BP208322 17/03/01 21:17:52 9325
  2. SUBROUTINE MAXROT(WTRAV,MFR,IB,IGAU,NBGMAT,NELMAT,NPINT)
  3. C
  4. C ROTATION DES TENSEURS POUR MAXWEL
  5. C SI MATE UNIDIRECTIONNEL
  6. C INSPIRE DE HOOKUN
  7. C
  8. IMPLICIT INTEGER(I-N)
  9. IMPLICIT REAL*8(A-H,O-Z)
  10. PARAMETER(UN=1.D0,DEUX=2.D0,UNDEMI=.5D0)
  11. C
  12. -INC CCHAMP
  13. -INC CCOPTIO
  14.  
  15.  
  16.  
  17. SEGMENT WTRAV
  18. REAL*8 DDAUX(LHOOK,LHOOK),VALMAT(NUMAT)
  19. REAL*8 VALCAR(NUCAR),DSIGT(NSTRS)
  20. REAL*8 TXR(IDIM,IDIM),DDHOMU(LHOOK,LHOOK)
  21. REAL*8 XLOC(3,3),XGLOB(3,3)
  22. REAL*8 D1HOOK(LHOOK,LHOOK),ROTHOO(LHOOK,LHOOK)
  23. ENDSEGMENT
  24. *
  25.  
  26.  
  27. C
  28. C INITIALISATION
  29. C
  30. CALL ZERO(XGLOB,IDIM,IDIM)
  31.  
  32. ***********************
  33. * CAS MASSIF
  34. ***********************
  35. *
  36. IF (MFR.EQ.1.OR.MFR.EQ.33) THEN
  37.  
  38. IF(IGAU.LE.NBGMAT) THEN
  39.  
  40. IF(IDIM.EQ.2)THEN
  41. C
  42. C DEFINITION DES AXES ORTHO./AXES LOCAUX
  43. C
  44.  
  45. XLOC(1,1)=VALMAT(2)
  46. XLOC(2,1)=VALMAT(3)
  47. XLOC(1,2)=-XLOC(2,1)
  48. XLOC(2,2)=XLOC(1,1)
  49.  
  50. ELSEIF(IDIM.EQ.3)THEN
  51. XLOC(1,1)=VALMAT(2)
  52. XLOC(2,1)=VALMAT(3)
  53. XLOC(3,1)=VALMAT(4)
  54. XLOC(1,2)=VALMAT(5)
  55. XLOC(2,2)=VALMAT(6)
  56. XLOC(3,2)=VALMAT(7)
  57. C
  58. CALL CROSS2(XLOC(1,1),XLOC(1,2),XLOC(1,3),IRR)
  59. IF(IRR.EQ.0) THEN
  60. CALL ERREUR(277)
  61. RETURN
  62. ENDIF
  63. ENDIF
  64.  
  65.  
  66.  
  67. C
  68. C DEFINITION DES AXES ORTHO./AXES GLOBAUX
  69. C
  70. DO 1045 K=1,IDIM
  71. DO 1045 J=1,IDIM
  72. DO 1045 I=1,IDIM
  73. XGLOB(K,J)=TXR(J,I)*XLOC(I,K)+XGLOB(K,J)
  74. 1045 CONTINUE
  75.  
  76. ENDIF
  77. C
  78. C
  79. ELSEIF (IGAU.LE.NBGMAT.AND.
  80. + (IB.LE.NELMAT.OR.NBGMAT.GT.1)) THEN
  81. C
  82. C MLR 12/8/99
  83. C POUR LE MOMENT ON BLOQUE LE RESTE
  84. C A BRANCHER PLUS TARD ...!
  85. C
  86. IZOB=0
  87. IF(IZOB.EQ.0) THEN
  88. CALL ERREUR (251)
  89. RETURN
  90. ENDIF
  91.  
  92.  
  93.  
  94.  
  95. C
  96. C Coques minces
  97. C
  98. IF(MFR.EQ.3) THEN
  99. C
  100. C CAS DKT INTEGRE
  101. C
  102. IF (NPINT.NE.0) THEN
  103. * CAS NON ENCORE IMPLEMENTE
  104. CALL ERREUR (251)
  105. RETURN
  106. ELSE
  107. C
  108. C CAS TRIDIMENSIONNEL ET FOURIER
  109. C
  110. IF(IFOUR.EQ.2.OR.IFOUR.EQ.1) THEN
  111. *PV D3HOO1 n'existe pas !!!
  112. * CALL ZERO(D3HOO1,3,3)
  113. COSA=VALMAT(2)
  114. SINA=VALMAT(3)
  115. C
  116. COS2 = COSA**2
  117. SIN2 = SINA**2
  118. SINCOS = SINA * COSA
  119. C
  120. C
  121. C
  122. C CAS AXISYMETRIQUE ET DEFORMATIONS PLANES
  123. C
  124. ELSE IF(IFOUR.EQ.0.OR.IFOUR.EQ.-1.OR.
  125. & IFOUR.EQ.-3) THEN
  126. COSA=VALMAT(2)
  127. SINA=VALMAT(3)
  128. COS2 = COSA**2
  129. SIN2 = SINA**2
  130. C
  131. C
  132. C CAS CONTRAINTES PLANES
  133. C
  134. ELSE IF(IFOUR.EQ.-2) THEN
  135. ******** DDHOOK(1,1)=YG1
  136. ******** DDHOOK(3,3)=YG1
  137. C
  138. ELSE
  139. IRET=0
  140. ENDIF
  141. ENDIF
  142. C
  143. C Coques epaisses
  144. C
  145. ELSE IF(MFR.EQ.5) THEN
  146. C
  147. C CAS TRIDIMENSIONNEL
  148. C
  149. IF(IFOUR.EQ.2) THEN
  150. COSA=VALMAT(2)
  151. SINA=VALMAT(3)
  152. C
  153. COS2 = COSA**2
  154. SIN2 = SINA**2
  155. SINCOS = SINA * COSA
  156. C
  157. C
  158. ELSE
  159. IRET=0
  160. ENDIF
  161. C
  162. C Coques minces avec cisaillement transverse
  163. C
  164. ELSE IF(MFR.EQ.9) THEN
  165. C
  166. IF(IFOUR.EQ.2) THEN
  167. COSA=VALMAT(2)
  168. SINA=VALMAT(3)
  169. C
  170. COS2 = COSA**2
  171. SIN2 = SINA**2
  172. SINCOS = SINA * COSA
  173. C
  174. C
  175. ELSE
  176. IRET=0
  177. ENDIF
  178. C
  179. C Cas des barres
  180. C
  181. ELSE IF(MFR.EQ.27) THEN
  182. C
  183. C RIEN A FAIRE
  184. C
  185. ELSE
  186. IRET=0
  187. GOTO 2035
  188. ENDIF
  189. ENDIF
  190. C
  191. 2035 RETURN
  192.  
  193. END
  194.  
  195.  
  196.  
  197.  
  198.  
  199.  
  200.  
  201.  
  202.  
  203.  
  204.  
  205.  
  206.  
  207.  
  208.  

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