Télécharger la plaquette Cast3M


fr en

Télécharger mshea1.eso

Retour à la liste

Numérotation des lignes : 

mshea1

  1. C MSHEA1    SOURCE    CHAT      05/01/13    01:52:55     5004
  2.       SUBROUTINE MSHEA1(DEP,ICAS,FOR,DPLUS,DMOIN, DTRANP,DTRANM,
  3.      >           APLUS,CURFP,NCURFP,CURKP,NCURKP,CURLP,NCURLP,
  4.      >           AMOIN,CURFM,NCURFM,CURKM,NCURKM,CURLM,NCURLM,
  5.      >           FCINI,FCAMP,KCINE)
  6.       IMPLICIT INTEGER(I-N)
  7.       IMPLICIT REAL*8(A-H,O-Z)
  8. C
  9. C======================================================================
  10. C       calcul des caracteristiques de la surface
  11. C       cinematique en tenant compte de la zone de transition
  12. C       WARNING: moyenne arithmetique pour l'amplitude
  13. C                et geometrique pour la raideur
  14. C2003???    WARNING: il faut iterer parceque KCINE0 depend de ALPHP et que
  15. C2003???             ALPHP depend de KCINE0
  16. C======================================================================
  17. C=======================================================================
  18. C
  19. C LISTE D'ECHANGE
  20. C ---------------
  21. C
  22. C DEP          = Deplacement
  23. C ICAS         = Numero de la branche
  24. C FOR          = Effort
  25. C DPLUS        = Deplacement de reference +
  26. C DMOIN        = Deplacement de reference -
  27. C DTRANP       = zone de transition  (sens positif)
  28. C DTRANM       =                     (sens négatif)
  29. C
  30. C CURFP,NCURFP = courbe de charge  +           (x>0,y>0)
  31. C CURKP,NCURKP = courbe de raideur +           (x>0,y>0)
  32. C CURLP,NCURLP = courbe de domaine elastique + (x>0,y>0)
  33. C CURFM,NCURFM = courbe de charge  -           (x<0,y<0)
  34. C CURKM,NCURKM = courbe de raideur -           (x<0,y>0)
  35. C CURLM,NCURLM = courbe de domaine elastique - (x<0,y>0)
  36. C
  37. C FCINI        = Force de reference pour la surface cinematique
  38. C FCAMP        = Taille de la surface cinematique
  39. C KCINE        = Raideur de la surface cinematique
  40. C2003 ALPHP        = Valeur de alpha +
  41. C
  42. C=======================================================================
  43. C
  44.       REAL*8 KCINE,KCINEP,KCINEM
  45.       PARAMETER (XZER=0.D0,UN=1.D0)
  46.       DIMENSION CURFP(2,NCURFP),CURKP(2,NCURKP),CURLP(2,NCURLP)
  47.       DIMENSION CURFM(2,NCURFM),CURKM(2,NCURKM),CURLM(2,NCURLM)
  48. C
  49.       FCINI=FOR
  50.       CALL YOFXCU(DPLUS,CURFP,NCURFP,  FPLUS ,DYSDX,KERRE)
  51.       FPLUS=FPLUS*(UN-APLUS)
  52.       CALL YOFXCU(DPLUS,CURLP,NCURLP,  FCAMPP,DYSDX,KERRE)
  53.       FCAMPP=FCAMPP*(UN-APLUS)
  54.       CALL YOFXCU(DPLUS,CURKP,NCURKP,  KCINEP,DYSDX,KERRE)
  55.       CALL YOFXCU(DMOIN,CURFM,NCURFM,  FMOIN ,DYSDX,KERRE)
  56.       FMOIN=FMOIN*(UN-AMOIN)
  57.       CALL YOFXCU(DMOIN,CURLM,NCURLM,  FCAMPM,DYSDX,KERRE)
  58.       FCAMPM=FCAMPM*(UN-AMOIN)
  59.       CALL YOFXCU(DMOIN,CURKM,NCURKM,  KCINEM,DYSDX,KERRE)
  60.       DDPLUS=DPLUS-FCAMPP/KCINEP
  61.       DDMOIN=DMOIN+FCAMPM/KCINEM
  62. C
  63.       IF (ICAS.EQ.1)THEN
  64. C     CHARGEMENT MONOTONE DIRECTION +
  65.         FCAMP=FCAMPP
  66.         KCINE=KCINEP
  67.       ELSEIF(ICAS.EQ.3)THEN
  68. C     BRANCHE DE UNLOADING ANELASTIQUE -
  69.         ALPHP=(DEP-DMOIN)/(DDPLUS-DMOIN)
  70.         FCAMP=ALPHP*FCAMPP+(UN-ALPHP)*FCAMPM
  71.         KCINE=FCAMP/(ALPHP*FCAMPP/KCINEP+(UN-ALPHP)*FCAMPM/KCINEM)
  72.       ELSEIF(ICAS.EQ.4)THEN
  73. C     CHARGEMENT MONOTONE DIRECTION -
  74.         FCAMP=FCAMPM
  75.         KCINE=KCINEM
  76.       ELSEIF(ICAS.EQ.6)THEN
  77. C     BRANCHE DE UNLOADING ANELASTIQUE +
  78.         ALPHP=(DEP-DDMOIN)/(DPLUS-DDMOIN)
  79.         FCAMP=ALPHP*FCAMPP+(UN-ALPHP)*FCAMPM
  80.         KCINE=FCAMP/(ALPHP*FCAMPP/KCINEP+(UN-ALPHP)*FCAMPM/KCINEM)
  81.       ENDIF
  82. C
  83. C 2003DEPP=DEP-FCINI/KCINE
  84. C 2003ALPHP=(DEPP-DTRANM)/(DTRANP-DTRANM)
  85. C 2003ALPHP=MAX(XZER,ALPHP)
  86. C 2003ALPHP=MIN(UN  ,ALPHP)
  87. C
  88.       RETURN
  89.       END
  90.  
  91.  
  92.  

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