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

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