Télécharger cficmo.eso

Retour à la liste

Numérotation des lignes :

  1. C CFICMO SOURCE BP208322 17/03/01 21:15:31 9325
  2. SUBROUTINE CFICMO(TPS1,TPS2,NBR,WR12,wrk52,wrk53,NWA,EFIC,TR,
  3. & E1,E2)
  4. *
  5. *==============================================================
  6. * calcul du module fictif du modele maxott
  7. *==============================================================
  8. *
  9. *
  10. * entrees
  11. *
  12. * XMAT(NCOMAT) = CARACTERISTIQUES MECANIQUES DU MATERIAU
  13. *
  14. * WRK12
  15. * EMi(2,NWA(i+1)) i=0 a 8
  16. * = EVOLUTION DU MODUE D'YOUNG DE LA BRANCHE i
  17. * SMi(NSTRS) i=1 a 8
  18. * = VARIABLES INTERNES AU DEBUT DU PAS D'INTEGRATION
  19. *
  20. * TPS1 = TEMPS AU DEBUT DU PAS D'INTEGRATION
  21. * TPS2 = TEMPS A LA FIN DU PAS D'INTEGRATION
  22. * NBR = INDICE DE LA CHAINE DE MAXWELL CONCERNEE
  23. * NWA = TABLEAU DES DIMENSIONS DES EVOLUTIONS DES MODULES DE CHAQUE CHAINE
  24. *
  25. *
  26. * sortie
  27. *
  28. * EFIC = MODULE FICTIF DE LA CHAINE DE MAXWELL
  29. * TR = TEMPS DE RELAXATION DE LA CHAINE DE MAXWELL
  30. * E1 = MODULE DE LA CHAINE DE MAXWELL AU DEBUT DU PAS
  31. * E2 = MODULE DE LA CHAINE DE MAXWELL A LA FIN DU PAS
  32. * KERRE = INDICATEUR D'ERREUR
  33. *
  34. *
  35. *
  36. *==============================================================
  37. *
  38. IMPLICIT INTEGER(I-N)
  39. IMPLICIT REAL*8(A-H,O-Z)
  40. -INC CCOPTIO
  41. -INC DECHE
  42. *
  43. SEGMENT WR12
  44. REAL*8 EM0(2,NWA(1)),EM1(2,NWA(2)),EM2(2,NWA(3))
  45. REAL*8 EM3(2,NWA(4)),EM4(2,NWA(5)),EM5(2,NWA(6))
  46. REAL*8 EM6(2,NWA(7)),EM7(2,NWA(8)),EM8(2,NWA(9))
  47. REAL*8 SM0(NSTRS),SM1(NSTRS),SM2(NSTRS),SM3(NSTRS)
  48. REAL*8 SM4(NSTRS),SM5(NSTRS),SM6(NSTRS),SM7(NSTRS)
  49. REAL*8 SM8(NSTRS)
  50. ENDSEGMENT
  51. *
  52. DIMENSION NWA(9)
  53. *
  54. * offsets pour caracteristiques temps de relaxation
  55. *
  56. * deformations planes / axisymetrique
  57. *
  58. IF ((IFOMOD.EQ.-1.AND.IFOUR.NE.-2).OR.
  59. & (IFOMOD.EQ.0.OR.IFOMOD.EQ.1)) THEN
  60. IPREMO = 10
  61. IPREMF = 51
  62. *
  63. * contraintes planes / 3D coques minces
  64. *
  65. ELSE IF((IFOMOD.EQ.-1.AND.IFOUR.EQ.-2).OR.
  66. & (IFOMOD.EQ.2.AND.(MFR.EQ.3.OR.MFR.EQ.9))) THEN
  67. IPREMO = 10
  68. IPREMF = 46
  69. *
  70. * tridimensionnelle massive
  71. *
  72. ELSE
  73. IPREMO = 15
  74. IPREMF = 57
  75. ENDIF
  76. *
  77. * test sur la branche
  78. *
  79. IF (NBR.EQ.0) THEN
  80. CALL MAXINT(EM0,NWA(1),TPS1,FTPS1,IRET1)
  81. CALL MAXINT(EM0,NWA(1),TPS2,FTPS2,IRET2)
  82. EFIC=FTPS2
  83. TR=0.
  84. E1=FTPS1
  85. E2=FTPS2
  86. GOTO 10
  87. ELSE IF (NBR.EQ.1) THEN
  88. CALL MAXINT(EM1,NWA(2),TPS1,FTPS1,IRET1)
  89. CALL MAXINT(EM1,NWA(2),TPS2,FTPS2,IRET2)
  90. TR=XMAT(IPREMO+2)
  91. ELSE IF (NBR.EQ.2) THEN
  92. CALL MAXINT(EM2,NWA(3),TPS1,FTPS1,IRET1)
  93. CALL MAXINT(EM2,NWA(3),TPS2,FTPS2,IRET2)
  94. TR=XMAT(IPREMO+4)
  95. ELSE IF (NBR.EQ.3) THEN
  96. CALL MAXINT(EM3,NWA(4),TPS1,FTPS1,IRET1)
  97. CALL MAXINT(EM3,NWA(4),TPS2,FTPS2,IRET2)
  98. TR=XMAT(IPREMO+6)
  99. ELSE IF (NBR.EQ.4) THEN
  100. CALL MAXINT(EM4,NWA(5),TPS1,FTPS1,IRET1)
  101. CALL MAXINT(EM4,NWA(5),TPS2,FTPS2,IRET2)
  102. TR=XMAT(IPREMO+8)
  103. *
  104. * 5 et plus
  105. *
  106. ELSE IF (NBR.EQ.5) THEN
  107. CALL MAXINT(EM5,NWA(6),TPS1,FTPS1,IRET1)
  108. CALL MAXINT(EM5,NWA(6),TPS2,FTPS2,IRET2)
  109. TR=XMAT(IPREMF+1)
  110. ELSE IF (NBR.EQ.6) THEN
  111. CALL MAXINT(EM6,NWA(7),TPS1,FTPS1,IRET1)
  112. CALL MAXINT(EM6,NWA(7),TPS2,FTPS2,IRET2)
  113. TR=XMAT(IPREMF+3)
  114. ELSE IF (NBR.EQ.7) THEN
  115. CALL MAXINT(EM7,NWA(8),TPS1,FTPS1,IRET1)
  116. CALL MAXINT(EM7,NWA(8),TPS2,FTPS2,IRET2)
  117. TR=XMAT(IPREMF+5)
  118. ELSE IF (NBR.EQ.8) THEN
  119. CALL MAXINT(EM8,NWA(9),TPS1,FTPS1,IRET1)
  120. CALL MAXINT(EM8,NWA(9),TPS2,FTPS2,IRET2)
  121. TR=XMAT(IPREMF+7)
  122. ENDIF
  123. *
  124. *
  125. E1=FTPS1
  126. E2=FTPS2
  127. IF(TPS2-TPS1.EQ.0.D0) THEN
  128. EFIC = FTPS2
  129. ELSE
  130. EFIC = FTPS2*(1.D0-EXP(-TR*(TPS2-TPS1)))/(TR*(TPS2-TPS1))
  131. ENDIF
  132. *
  133. *
  134. 10 IF (IRET1.EQ.0) THEN
  135. CALL ERREUR(854)
  136. RETURN
  137. ELSE IF (IRET2.EQ.0) THEN
  138. CALL ERREUR(854)
  139. RETURN
  140. ENDIF
  141. END
  142.  
  143.  
  144.  
  145.  
  146.  
  147.  
  148.  
  149.  
  150.  
  151.  
  152.  
  153.  
  154.  
  155.  
  156.  
  157.  
  158.  

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