Télécharger mshea1.eso

Retour à la liste

Numérotation des lignes :

  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