Télécharger plasec.eso

Retour à la liste

Numérotation des lignes :

  1. C PLASEC SOURCE CB215821 16/04/21 21:18:02 8920
  2. C PLASEC SOURCE INSL 24/10/96
  3. SUBROUTINE PLASEC(SIG0,VAR0,DEPST,SIGF,VARF,XMAT,NSTRSS,NMATT,
  4. 1 TRAC,LTRA,MFR,NVARI,CMATE,XCAR,DDHOOK,NCARR,IFOUR)
  5. C----------------------------------------------------------------------
  6. C PLASTICITE MODELE ISOTROPE
  7. C
  8. C ENTREES
  9. C SIG0(NSTRS) = CONTRAINTES INITIALES
  10. C NSTRS = NOMBRE DE CONTRAINTES
  11. C DSIGT(NSTRS) = INCREMENT DE CONTRAINTES TOTALES
  12. C VAR0(NVARI) = VARIABLES INTERNES DEBUT
  13. C VAR0( 1 ) = EPSP(1)
  14. C VAR0( 2 ) = EPSP(2)
  15. C VAR0( 3 ) = EPSP(3)
  16. C XMAT(NMATT) = COMPOSANTES DE MATERIAU
  17. C IVAL(NMATT) = INDICE DES COMPOSANTES DE MATERIAU
  18. C NMATT = NOMBRE DE COMPOSANTES DE MATERIAU
  19. C SORTIES
  20. C SIGF(NSTRS) = CONTRAINTES FINALES
  21. C VARF(NVARI) = VARIABLES INTERNES FINALES
  22. C-----------------------------------------------------------------------
  23. C
  24. IMPLICIT INTEGER(I-N)
  25. IMPLICIT REAL*8(A-H,O-Z)
  26. *
  27. CHARACTER*8 CMATE
  28. DIMENSION SIG0(NSTRSS),VAR0(NVARI),SIGF(NSTRSS),VARF(NVARI)
  29. DIMENSION TRAC(LTRA),XMAT(NMATT),XCAR(NCARR),DEPST(NSTRSS)
  30. DIMENSION EPSLO(3),EPSL1(3),DEPSL(3),S1(3),SIGR(3),YOU(2)
  31. DIMENSION DDHOOK(NSTRSS,NSTRSS),DDHOMU(6,6),DSIGT(6),SIGL0(3)
  32. IRTD=1
  33. C-----------------------------------------------------------------------
  34. IF(CMATE.EQ.'ORTHOTRO'.OR.CMATE.EQ.'UNIDIREC') GOTO 4
  35. GOTO 46
  36. C
  37. 4 CONTINUE
  38. C
  39. C PLASTICITE (DECOUPLE) ORTHOTROPE OU UNIDIRECTIONNELLE
  40. C
  41. CALL ZERO(EPSLO,3,1)
  42. CALL ZERO(DEPSL,3,1)
  43. CALL ZERO(EPSL1,3,1)
  44. CALL ZERO(S1,3,1)
  45. CALL ZERO(SIGR,3,1)
  46. CALL ZERO(SIGL0,3,1)
  47. CALL ZERO(YOU,2,1)
  48. *
  49. NSTRS=3
  50. EPAIST=1.D0
  51. IF(MFR.EQ.3.AND.NSTRSS.GT.3) THEN
  52. EPAIST=XCAR(1)
  53. ENDIF
  54. *
  55. DO I=1,NSTRS
  56. SIGR(I)=SIG0(I)/EPAIST
  57. END DO
  58. C
  59. C DEFINITION DES AXES ORTHO./AXES LOCAUX
  60. C
  61. IF(CMATE.EQ.'ORTHOTRO') THEN
  62. YOU(1) =XMAT(1)
  63. YOU(2) =XMAT(2)
  64. A =XMAT(5)
  65. B =XMAT(6)
  66. ELSE
  67. YOU(1) =XMAT(1)
  68. A =XMAT(2)
  69. B =XMAT(3)
  70. ENDIF
  71. C--------------------------------------------------------------------
  72. C ROTATION DES CONTRAINTES ET DES DEFORMATIONS
  73. C DANS LE REPERE D'ORTHOTROPIE
  74. C
  75. PI=4.D0*ATAN(1.D0)
  76. PHI=(ATAN(B/A)*180.D0/PI)-90.D0
  77. CALL DDRF(DEPST,PHI,DEPSL)
  78. CALL CDRF(SIGR,PHI,SIGL0)
  79. C
  80. IF(CMATE.EQ.'UNIDIREC'.OR.CMATE.EQ.'ORTHOTRO') THEN
  81. EPSL1(1)=VAR0(1)+SIGL0(1)/YOU(1)+DEPSL(1)
  82. CALL COUPLA(EPSL1(1),TRAC,LTRA,YOU(1),S1(1),PENTE)
  83. VARF(1)=EPSL1(1)-S1(1)/YOU(1)
  84. ENDIF
  85. *
  86. IF(CMATE.EQ.'ORTHOTRO'.AND.YOU(2).NE.0.D0) THEN
  87. EPSL1(2)=VAR0(2)+SIGL0(2)/YOU(2)+DEPSL(2)
  88. CALL COUPLA(EPSL1(2),TRAC,LTRA,YOU(2),S1(2),PENTE)
  89. VARF(2)=EPSL1(2)-S1(2)/YOU(2)
  90. ENDIF
  91. *
  92. CALL CDRG(S1,PHI,SIGR)
  93. DO I=1,NSTRS
  94. SIGF(I)=SIGR(I)*EPAIST
  95. END DO
  96. IF(NSTRSS.EQ.4) THEN
  97. SIGF(3)=0.D0
  98. SIGF(4)=SIGR(3)*EPAIST
  99. ENDIF
  100. *
  101. IF(NSTRSS.EQ.6) THEN
  102. CALL ZERO(DDHOMU,6,6)
  103. CALL ZERO(DDHOOK,NSTRSS,NSTRSS)
  104. CALL ZERO(DSIGT,6,1)
  105. *
  106. CALL DOHCOM(XMAT,NMATT,CMATE,IFOUR,NSTRSS,DDHOOK,IRTD)
  107. CALL HOOKMU(EPAIST,0.D0,NSTRSS,DDHOOK,DDHOMU)
  108. CALL BST(DDHOMU,DEPST,NSTRSS,NSTRSS,DSIGT)
  109. *
  110. DO I = 4,NSTRSS
  111. SIGF(I)=SIG0(I)+DSIGT(I)
  112. END DO
  113. ENDIF
  114. GOTO 510
  115. C-----------------------------------------------------------------------
  116. C
  117. C PLASTICITE ISOTROPE UNIAXIALE
  118. C
  119. 46 CONTINUE
  120. SECT=XCAR(1)
  121. YOU1=XMAT(1)
  122. C
  123. EPSX=VAR0(1)+SIG0(1)/YOU1/SECT+DEPST(1)
  124. CALL COUPLA(EPSX,TRAC,LTRA,YOU1,SIGRX,PENTE)
  125. VARF(1)=EPSX-SIGRX/YOU1
  126. SIGF(1)=SIGRX*SECT
  127. C
  128. 510 CONTINUE
  129. C-----------------------------------------------------------------------
  130. 1991 FORMAT(18(1X,E10.3))
  131. 2007 FORMAT(7(6(1X,E12.5),/))
  132. 2008 FORMAT(7(6(1X,I12),/))
  133. 1000 FORMAT(5X,' *** ATTENTION ELEMENT NON COMPATIBLE AVEC CE MODELE'
  134. *,/,10X,'ROUTINE = PLASEC : ELEMENT = ',I3)
  135. RETURN
  136. END
  137.  
  138.  
  139.  
  140.  
  141.  
  142.  
  143.  

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