Télécharger cuniac.eso

Retour à la liste

Numérotation des lignes :

  1. C CUNIAC SOURCE PV 17/12/08 21:17:09 9660
  2. C UNIACI SOURCE AM 96/12/24 21:33:08 2448
  3. SUBROUTINE CUNIAC(wrk52,wrk53,wrk54,NSTRS1,MFR1)
  4. IMPLICIT INTEGER(I-N)
  5. IMPLICIT REAL*8(A-H,O-Z)
  6. C----------------------------------------------------------------------
  7. C
  8. C COMPORTEMENT UNIDIRECTIONNEL ACIER_UNI (2D ou coque)
  9. C (D'APRES FIBSTE DU MODELE A FIBRE)
  10. C
  11. C----------------------------------------------------------------------
  12. C Pierre Pegon (ISPRA) Juillet 1994
  13. C----------------------------------------------------------------------
  14. C XMAT( 1) YOUNG
  15. C XMAT( 2) NU
  16. C XMAT( 3) V1X
  17. C XMAT( 4) V1Y
  18. C XMAT( 5) MAX. ELAST. STRESS
  19. C XMAT( 6) ULTIMATE STRESS
  20. C XMAT( 7) HARDENING STRAIN
  21. C XMAT( 8) ULTIMATE STRAIN
  22. C XMAT( 9) Ro FACTOR
  23. C XMAT(10) b FACTOR
  24. C XMAT(11) a1 FACTOR
  25. C XMAT(12) a2 FACTOR
  26. C XMAT(13) FACTOR L/D
  27. C XMAT(14) a6 FACTOR
  28. C XMAT(15) c FACTOR
  29. C XMAT(16) a FACTOR
  30. C----------------------------------------------------------------------
  31. C VAR0( 1) VARF( 1) FLAG
  32. C VAR0( 2) VARF( 2) EPSO
  33. C VAR0( 3) VARF( 3) R FACTOR
  34. C VAR0( 4) VARF( 4) TANGENT MODULUS
  35. C VAR0( 5) VARF( 5) LAST REV. STRAIN
  36. C VAR0( 6) VARF( 6) INTERSEPTION POINT
  37. C VAR0( 7) VARF( 7) MAX. REV. POST. STRAIN
  38. C VAR0( 8) VARF( 8) MIN. REV. POST. STRAIN
  39. C VAR0( 9) VARF( 9) LAST REV. STRESS
  40. C VAR0(10) VARF(10) INTERSEPTION POINT
  41. C VAR0(11) VARF(11) MAX. REV. POST. STRESS
  42. C VAR0(12) VARF(12) MIN. REV. POST. STRESS
  43. C----------------------------------------------------------------------
  44. C xcarb( 1) EPAIS (if any)
  45. C----------------------------------------------------------------------
  46. -INC CCOPTIO
  47. -INC DECHE
  48.  
  49. DIMENSION EPSON(5),SIGMA(5)
  50. DIMENSION COBMA(10)
  51. *
  52. KERRE=0
  53. *
  54. C
  55. C Calcul des contraintes unidirectionnelles effectives
  56. C
  57. YOUN=XMAT(1)
  58. V1X=XMAT(3)
  59. V1Y=XMAT(4)
  60. IF (IFOUR.EQ.2.AND.(MFR1.EQ.3.OR.MFR1.EQ.9))THEN
  61. EPAI=xcarb(1)
  62. S0IUNI=(SIG0(1)*V1X**2+SIG0(2)*V1Y**2+V1X*V1Y*SIG0(3))/EPAI
  63. DSIUNI=(DEPST(1)*V1X**2+DEPST(2)*V1Y**2+V1X*V1Y*DEPST(3))
  64. DSIUNI=YOUN*DSIUNI
  65. ELSEIF (IFOUR.LE.0.AND.MFR1.EQ.1)THEN
  66. S0IUNI=SIG0(1)*V1X**2+SIG0(2)*V1Y**2+V1X*V1Y*SIG0(4)
  67. DSIUNI=DEPST(1)*V1X**2+DEPST(2)*V1Y**2+V1X*V1Y*DEPST(4)
  68. DSIUNI=YOUN*DSIUNI
  69. ELSE
  70. KERRE=99
  71. RETURN
  72. ENDIF
  73. C
  74. EYOUN = YOUN
  75. FSY__ = XMAT( 5)
  76. FSU__ = XMAT( 6)
  77. EPSH_ = XMAT( 7)
  78. EPSU_ = XMAT( 8)
  79. ROFAC = XMAT( 9)
  80. BFACT = XMAT(10)
  81. A1FAC = XMAT(11)
  82. A2FAC = XMAT(12)
  83. FACLD = XMAT(13)
  84. A6FAC = XMAT(14)
  85. CFACT = XMAT(15)
  86. AFACT = XMAT(16)
  87. C
  88. IFLAG = NINT(VAR0(1))
  89. STRAI = VAR0(2) + DSIUNI/EYOUN
  90. RFACT = VAR0(3)
  91. TANGE = VAR0(4)
  92. C
  93. DO IFACT = 1, 2
  94. EPSON(IFACT) = VAR0(4+IFACT)
  95. SIGMA(IFACT) = VAR0(8+IFACT)
  96. END DO
  97. C
  98. EPSON(3) = VAR0(2)
  99. SIGMA(3) = S0IUNI
  100. C
  101. DO IFACT = 4, 5
  102. EPSON(IFACT) = VAR0(3+IFACT)
  103. SIGMA(IFACT) = VAR0(7+IFACT)
  104. END DO
  105. C
  106. CALL STEEL1
  107. . ( EYOUN, FSY__, FSU__, EPSH_, EPSU_, ROFAC, BFACT,
  108. . A1FAC, A2FAC, FACLD, A6FAC, CFACT, AFACT, IFLAG,
  109. . STRAI, RFACT, EPSON, SIGMA, TANGE, STRES )
  110. C
  111. VARF(1) = DBLE(IFLAG)
  112. VARF(2) = STRAI
  113. VARF(3) = RFACT
  114. VARF(4) = TANGE
  115. C
  116. DO IFACT = 1, 2
  117. VARF(4+IFACT) = EPSON(IFACT)
  118. VARF(8+IFACT) = SIGMA(IFACT)
  119. END DO
  120. C
  121. DO IFACT = 4, 5
  122. VARF(3+IFACT) = EPSON(IFACT)
  123. VARF(7+IFACT) = SIGMA(IFACT)
  124. END DO
  125. C
  126. IF (IFOUR.EQ.2.AND.(MFR1.EQ.3.OR.MFR1.EQ.9))THEN
  127. EPAI=xcarb(1)
  128. SIGF(1)=EPAI*STRES*V1X**2
  129. SIGF(2)=EPAI*STRES*V1Y**2
  130. SIGF(3)=EPAI*STRES*V1X*V1Y
  131. DO ISTRS=4,NSTRS1
  132. C PPDC
  133. C SIGF(ISTRS)=SIG0(ISTRS)+DSIGT(ISTRS)
  134. SIGF(ISTRS)=SIG0(ISTRS)
  135. C PPDC
  136. ENDDO
  137. ELSEIF (IFOUR.LE.0.AND.MFR1.EQ.1)THEN
  138. SIGF(1)=STRES*V1X**2
  139. SIGF(2)=STRES*V1Y**2
  140. SIGF(3)=0.D0
  141. SIGF(4)=STRES*V1X*V1Y
  142. ENDIF
  143. C
  144. DO ISTRS=1,NSTRS1
  145. DEFP(ISTRS)=DEPST(ISTRS)-((SIGF(ISTRS)-SIG0(ISTRS)))/EYOUN
  146. END DO
  147. C
  148. 1000 RETURN
  149. END
  150.  
  151.  
  152.  
  153.  
  154.  
  155.  
  156.  
  157.  
  158.  
  159.  
  160.  

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