Télécharger la plaquette Cast3M


fr en

Télécharger reelle.eso

Retour à la liste

Numérotation des lignes : 

reelle

  1. C REELLE    SOURCE    CB215821  16/04/21    21:18:10     8920
  2. C     REELLE    SOURCE   NECHNECH   99/06/08
  3.       SUBROUTINE REELLE(S1,DEP,DEPSITF1,DEPSITF2,DEPE,S2,ENDM,
  4.      &            NSTRS,IFIS,IPLA,BETJEF,NECH0,NECH1)
  5. C
  6.       IMPLICIT INTEGER(I-N)
  7.       IMPLICIT REAL*8(A-H,O-Z)
  8.       DIMENSION S1(4), S2(4), V1(4)
  9.       DIMENSION DEP(4,4),DEPE(4,4)
  10. *      COMMON /DBETJEF/AA,BETA,RB,ALFA,YOUN,XNU,GFC,GFT,CAR,ETA,
  11. *     &          TDEF,TCON,DPSTF1,DPSTF2,TETA,PDT,ICT,ICC,IMOD,IVISS,
  12. *     &                ITER,ISIM,IBB,IGAU1,IZON
  13. *      COMMON /DBETNECH0/DT,DC,ALFG,S0
  14. *      COMMON /DBETNECH1/ENDL
  15.       SEGMENT BETJEF
  16.        REAL*8 AA,BETA,RB,ALPHA,EX,XNU,GFC,GFT,CAR,ETA,TDEF,
  17.      &          TCON,DPSTF1,DPSTF2,TETA,PDT,TP0
  18.        INTEGER ICT,ICC,IMOD,IVIS,ITER,
  19.      &                ISIM,IBB,IGAU,IZON
  20.       ENDSEGMENT
  21.       SEGMENT NECH0
  22.        REAL*8 DT,DC,ALFG,S0,ENDO
  23.       ENDSEGMENT
  24.       SEGMENT NECH1
  25.        REAL*8 ENDL
  26.       ENDSEGMENT
  27. C
  28. *--------------------------------------------------------------------*
  29. *    CALCUL  LA MATRICE ELASTO-PLASTIQUE ENDOMMAGEABLE ET            *
  30. *    DE LA CONTRAINTE REELLE                                         *
  31. *                                                                    *
  32. *                                                                    *
  33. *    ENTREES:                                                        *
  34. *                                                                    *
  35. *        S1        CONTRAINTE EFFECTIVE                              *
  36. *        NSTRS     NOMBRE DE CONTRAINTES                             *
  37. *        DPSTF1    DEFORMATION EQUIVALENTE EN TRACTION               *
  38. *        DPSTF2    DEFORMATION EQUIVALENTE EN COMPRESSION            *
  39. *        DEP       MATRICE ELSATO-PLASTIQUE                          *
  40. *                                                                    *
  41. *    SORTIES:                                                        *
  42. *                                                                    *
  43. *        S2        CONTRAINTE REELLE                                 *
  44. *        DEPE      MATRICE ELASTO-PLASTIQUE ENDOMMAGEABLE            *
  45. *        ENDO      ENDOMMAGEMENT                                     *
  46. *--------------------------------------------------------------------*
  47. C       INITIALISATION
  48. C
  49. C      S0 = 0.25
  50.       ENDT = 0.D0
  51.       ENDC = 0.D0
  52.       ENDM = 0.D0
  53.            NN=NSTRS-1
  54.         CALL ZERO(DEPE,4,4)
  55.         CALL ZERO(S2,4,1)
  56.         SOM1 = 0.D0
  57.         SOM2 = 0.D0
  58. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  59. C       CALCUL DE L'ENDOMMAGEMENT
  60. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  61. C       WRITE(*,*)'K1=',DEPSITF1
  62. C       WRITE(*,*)'K2=',DEPSITF2
  63. C
  64.            CALL DAMAG1(DEPSITF1,DEPSITF2,ENDC,
  65.      &          ENDT,IFIS,IPLA,3,BETJEF,NECH0,NECH1)
  66. C
  67. C
  68. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  69. C       CALCUL DES VALEURS PROPRES DE LA MATRICE DE VECTEUR CONTRAINTE
  70. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  71. C
  72.            CALL PRINC(S1,V1,NSTRS)
  73. C
  74.            DO 10 I=1, NN
  75.             SOM1=SOM1+(V1(I)+ABS(V1(I)))/2.D0
  76.             SOM2=SOM2+ABS(V1(I))
  77.    10 CONTINUE
  78.         IF(SOM2.EQ.0.D0) THEN
  79.               SOM2=SOM2+1.D-10
  80.              ENDIF
  81.            RSIG=SOM1/SOM2
  82.            RSIG=ABS(RSIG)
  83.            S=S0+(1-S0)*RSIG
  84. C      S = 1.D0
  85. C
  86. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  87. C       CALCUL DE L'ENDOMMAGEMENT TOTAL
  88. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  89. C
  90.            ENDM=1.D0-((1.D0-ENDC)*(1.D0-S*ENDT))
  91. C
  92.            IF((ENDM.LT.0.D0).OR.(ENDM.GT.1D0))THEN
  93.              WRITE(*,*)'ENDO=',ENDM
  94. C            STOP
  95.            ENDIF
  96. C
  97. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  98. C
  99.         DO 20 I=1, NSTRS
  100.            S2(I)=(1.D0-ENDM)*S1(I)
  101.    20 CONTINUE
  102. C
  103. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  104. C       CALCUL DE LA MATRICE ELASTO-PLASTIQUE ENDOMMAGEABLE
  105. C~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
  106. C
  107.       DO 30 I=1, NSTRS
  108.       DO 40 J=1, NSTRS
  109.        DEPE(I,J)=(1.D0-ENDM)*DEP(I,J)
  110.    40 CONTINUE
  111.    30 CONTINUE
  112. C
  113.       RETURN
  114.       END
  115.  
  116.  
  117.  
  118.  
  119.  
  120.  

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