Télécharger cmaxro.eso

Retour à la liste

Numérotation des lignes :

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

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