Télécharger la plaquette Cast3M


fr en

Télécharger cficmo.eso

Retour à la liste

Numérotation des lignes : 

cficmo

  1. C CFICMO    SOURCE    PV        17/12/08    21:15:38     9660           
  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.  
  41. -INC PPARAM
  42. -INC CCOPTIO
  43. -INC DECHE
  44. *
  45.       SEGMENT WR12
  46.         REAL*8  EM0(2,NWA(1)),EM1(2,NWA(2)),EM2(2,NWA(3))
  47.         REAL*8  EM3(2,NWA(4)),EM4(2,NWA(5)),EM5(2,NWA(6))
  48.         REAL*8  EM6(2,NWA(7)),EM7(2,NWA(8)),EM8(2,NWA(9))
  49.         REAL*8  SM0(NSTRS),SM1(NSTRS),SM2(NSTRS),SM3(NSTRS)
  50.         REAL*8  SM4(NSTRS),SM5(NSTRS),SM6(NSTRS),SM7(NSTRS)
  51.         REAL*8  SM8(NSTRS)
  52.       ENDSEGMENT
  53. *
  54.        DIMENSION NWA(9)
  55. *
  56. * offsets pour caracteristiques temps de relaxation
  57. *
  58. * deformations planes / axisymetrique
  59. *
  60.       IF ((IFOMOD.EQ.-1.AND.IFOUR.NE.-2).OR.
  61.      &        (IFOMOD.EQ.0.OR.IFOMOD.EQ.1)) THEN
  62.         IPREMO = 10
  63.         IPREMF = 51
  64. *
  65. * contraintes planes / 3D coques minces
  66. *
  67.       ELSE IF((IFOMOD.EQ.-1.AND.IFOUR.EQ.-2).OR.
  68.      &        (IFOMOD.EQ.2.AND.(MFR.EQ.3.OR.MFR.EQ.9))) THEN
  69.         IPREMO = 10
  70.         IPREMF = 46
  71. *
  72. * tridimensionnelle massive
  73. *
  74.       ELSE
  75.         IPREMO = 15
  76.         IPREMF = 57
  77.       ENDIF
  78. *
  79. *     test sur la branche
  80. *
  81.        IF (NBR.EQ.0) THEN
  82.          CALL MAXINT(EM0,NWA(1),TPS1,FTPS1,IRET1)
  83.          CALL MAXINT(EM0,NWA(1),TPS2,FTPS2,IRET2)
  84.          EFIC=FTPS2
  85.          TR=0.
  86.          E1=FTPS1
  87.          E2=FTPS2
  88.          GOTO 10
  89.        ELSE IF (NBR.EQ.1) THEN
  90.          CALL MAXINT(EM1,NWA(2),TPS1,FTPS1,IRET1)
  91.          CALL MAXINT(EM1,NWA(2),TPS2,FTPS2,IRET2)
  92.          TR=XMAT(IPREMO+2)
  93.        ELSE IF (NBR.EQ.2) THEN
  94.          CALL MAXINT(EM2,NWA(3),TPS1,FTPS1,IRET1)
  95.          CALL MAXINT(EM2,NWA(3),TPS2,FTPS2,IRET2)
  96.          TR=XMAT(IPREMO+4)
  97.        ELSE IF (NBR.EQ.3) THEN
  98.          CALL MAXINT(EM3,NWA(4),TPS1,FTPS1,IRET1)
  99.          CALL MAXINT(EM3,NWA(4),TPS2,FTPS2,IRET2)
  100.          TR=XMAT(IPREMO+6)
  101.        ELSE IF (NBR.EQ.4) THEN
  102.          CALL MAXINT(EM4,NWA(5),TPS1,FTPS1,IRET1)
  103.          CALL MAXINT(EM4,NWA(5),TPS2,FTPS2,IRET2)
  104.          TR=XMAT(IPREMO+8)
  105. *
  106. *    5 et plus
  107. *
  108.        ELSE IF (NBR.EQ.5) THEN
  109.          CALL MAXINT(EM5,NWA(6),TPS1,FTPS1,IRET1)
  110.          CALL MAXINT(EM5,NWA(6),TPS2,FTPS2,IRET2)
  111.          TR=XMAT(IPREMF+1)
  112.        ELSE IF (NBR.EQ.6) THEN
  113.          CALL MAXINT(EM6,NWA(7),TPS1,FTPS1,IRET1)
  114.          CALL MAXINT(EM6,NWA(7),TPS2,FTPS2,IRET2)
  115.          TR=XMAT(IPREMF+3)
  116.        ELSE IF (NBR.EQ.7) THEN
  117.          CALL MAXINT(EM7,NWA(8),TPS1,FTPS1,IRET1)
  118.          CALL MAXINT(EM7,NWA(8),TPS2,FTPS2,IRET2)
  119.          TR=XMAT(IPREMF+5)
  120.        ELSE IF (NBR.EQ.8) THEN
  121.          CALL MAXINT(EM8,NWA(9),TPS1,FTPS1,IRET1)
  122.          CALL MAXINT(EM8,NWA(9),TPS2,FTPS2,IRET2)
  123.          TR=XMAT(IPREMF+7)
  124.        ENDIF
  125. *
  126. *
  127.        E1=FTPS1
  128.        E2=FTPS2
  129.        IF(TPS2-TPS1.EQ.0.D0) THEN
  130.          EFIC = FTPS2
  131.        ELSE
  132.          EFIC = FTPS2*(1.D0-EXP(-TR*(TPS2-TPS1)))/(TR*(TPS2-TPS1))
  133.        ENDIF
  134. *
  135. *
  136. 10     IF (IRET1.EQ.0) THEN
  137.          CALL ERREUR(854)
  138.          RETURN
  139.        ELSE IF (IRET2.EQ.0) THEN
  140.          CALL ERREUR(854)
  141.          RETURN
  142.        ENDIF
  143.        END
  144.  
  145.  
  146.  
  147.  
  148.  
  149.  
  150.  
  151.  
  152.  
  153.  
  154.  
  155.  
  156.  
  157.  
  158.  
  159.  
  160.  
  161.  

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