Télécharger takep2.eso

Retour à la liste

Numérotation des lignes :

takep2
  1. C TAKEP2 SOURCE CHAT 05/01/13 03:31:53 5004
  2. C TAKEPO SOURCE AM 95/12/27 22:09:50 1962
  3. SUBROUTINE TAKEP2(SIG0,NSTRS,DEPST,VAR0,XMAT,NCOMAT,XCAR,TRAC,
  4. . NCOURB,SIGF,VARF,DEFP,KERRE)
  5. C-----------------------------------------------------------------------
  6. C
  7. C ROUTINE DE CALCUL DES CONTRAINTES DE LA ROTULE
  8. C POUTRE 2D
  9. C CAS DU TAKEDA TRILINEAIRE SUR LE MOMENT
  10. C
  11. C ENTREE
  12. C SIG0 CONTRAINTES INITIALES
  13. C DEPST INCREMENT OF GENERALIZED STRAIN
  14. C XMAT PROPRIETES DU MATERIAU
  15. C XCAR CARACTERISTIQUE DU MATERIAU
  16. C VAR0 VARFABLES INTERNES AU DEBUT
  17. C TRAC COURBE DE CHARGEMENT
  18. C SORTIE
  19. C SIGF CONTRAINTES A LA FIN
  20. C VARF VARFABLES INTERNES A LA FIN
  21. C DEFP DEFORMATION PLASTIQUE
  22. C
  23. C materiau takeda PP 17/9/92, AVRIL 92 A.P. et P.P.
  24. C PP renewed new ECOU 24/11/95
  25. C-----------------------------------------------------------------------
  26. IMPLICIT INTEGER(I-N)
  27. IMPLICIT REAL*8(A-H,O-Z)
  28. *
  29. DIMENSION SIG0(*),DEPST(*),VAR0(*),XMAT(*),XCAR(*),TRAC(*),
  30. . SIGF(*),VARF(*) ,DEFP(*)
  31. C
  32. C-----------------------------------------------------------------------
  33. C CARACTERISTIQUES
  34. C ================
  35. C
  36. C SECT = XCAR(1) SECTION
  37. C INRZ = XCAR(2) Z-INERTIAL BENDING
  38. C SRDY = XCAR(3) TRANSVERSE SHEAR'S REDUCED SECTION
  39. C-----------------------------------------------------------------------
  40. C VARIABLES INTERNES
  41. C ==================
  42. C
  43. C IFR1 = VARF(1)
  44. C EITA = VARF(2)
  45. C UEP1 = VARF(3)
  46. C PEP1 = VARF(4)
  47. C UEN1 = VARF(5)
  48. C PEN1 = VARF(6)
  49. C UCP1 = VARF(7)
  50. C UCN1 = VARF(8)
  51. C E1P1 = VARF(9)
  52. C E1N1 = VARF(10)
  53. C E2P1 = VARF(11)
  54. C E2N1 = VARF(12)
  55. C USP1 = VARF(13)
  56. C USN1 = VARF(14)
  57. C UTP1 = VARF(15)
  58. C PTP1 = VARF(16)
  59. C UTN1 = VARF(17)
  60. C PTN1 = VARF(18)
  61. C EP11 = VARF(19)
  62. C EN11 = VARF(20)
  63. C DES0 = VARF(21)
  64. C-----------------------------------------------------------------------
  65. C ON RECUPERE LES PARAMETRES DU MODELE
  66. C ====================================
  67. C
  68. C E0 = XMAT(1)
  69. C NU = XMAT(2)
  70. C SFDP = XMAT(6)
  71. C SFDN = XMAT(7)
  72. C PINP = XMAT(8)
  73. C PINN = XMAT(9)
  74. C SRDP = XMAT(10)
  75. C SRDN = XMAT(11)
  76. C
  77. C
  78. IF(NCOURB.EQ.4)THEN
  79. C
  80. C LA COURBE EST SYMETRIQUE
  81. C
  82. EIY=TRAC(3)/TRAC(4)
  83. CMOP=TRAC(3)
  84. CMON=-CMOP
  85. YCUP=TRAC(6)
  86. YCUN=-YCUP
  87. YMOP=TRAC(5)
  88. YMON=-YMOP
  89. EAYP=(TRAC(7)-TRAC(5))/(TRAC(8)-TRAC(6))
  90. EAYN=EAYP
  91. C
  92. XMAT(7)=XMAT(6)
  93. XMAT(9)=-XMAT(8)
  94. XMAT(11)=XMAT(10)
  95. C
  96. ELSEIF(NCOURB.EQ.7)THEN
  97. C
  98. C LA COURBE EST NON SYMETRIQUE
  99. C
  100. EIY=(TRAC(9)-TRAC(5)) / (TRAC(10)-TRAC(6))
  101. CMOP=TRAC(9)
  102. CMON=TRAC(5)
  103. YCUP=TRAC(12)
  104. YCUN=TRAC(4)
  105. YMOP=TRAC(11)
  106. YMON=TRAC(3)
  107. EAYP=(TRAC(13)-TRAC(11))/(TRAC(14)-TRAC(12))
  108. EAYN=(TRAC(3) -TRAC(1) )/(TRAC(4) -TRAC(2))
  109. C
  110. ELSE
  111. C
  112. C LA COURBE EST ERRONEE
  113. C
  114. KERRE=31
  115. RETURN
  116. C
  117. ENDIF
  118. C
  119. C-----------------------------------------------------------------------
  120. C
  121. C TEST ON ELASTIC STIFFNESS
  122. C
  123. EIYY0 = XMAT(1)*XCAR(2)
  124. IF(ABS(EIY-EIYY0)/EIY.GT.1.D-6)THEN
  125. KERRE=36
  126. RETURN
  127. ENDIF
  128. C
  129. XMOMY = SIG0(3)
  130. C
  131. C APPELE AU PROGRAME QUI CALCULE LE MOMENT ( XMOMY )
  132. C CORRESPONDENT A L' INCREMENT LINEAIRE DE MOMENT - DP1 -
  133. C >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>
  134. DO 1 I=1,21
  135. VARF(I)=VAR0(I)
  136. 1 CONTINUE
  137. C
  138. C CONVERSION ENTIER
  139. C
  140. IFR1=INT(VARF(1))
  141. C
  142. C INITIALISATION DE IFR1 ET DE VARF(.,2)
  143. C
  144. IFC1=0
  145. IF(IFR1.EQ.0) THEN
  146. IFR1 = 1
  147. VARF(2)=EIYY0
  148. ENDIF
  149. C
  150. C DEFORMATIONS GENERALISEES
  151. C =========================
  152. C
  153. C INCREMENT DE MOMENT AUTOUR DE YY
  154. C
  155. C DP1 = D_SIGMA / LINEAR_STIFFNESS * CURRENT STIFFNESS
  156. C
  157. C DP1 = DSIGT(3) / EIYY0 * VARF(2)
  158. DP1 = DEPST(3) * VARF(2)
  159. C
  160. C
  161. C===================================================================
  162. CALL DDNSTH(IFR1,IFC1,DP1,XMOMY,VARF(21),
  163. * VARF(2),EIYY0,EAYP,EAYN,CMOP,CMON,YMOP,
  164. * YMON,YCUP,YCUN,XMAT(6),XMAT(7),XMAT(8),XMAT(9),
  165. * XMAT(10),XMAT(11),VARF(3),VARF(4),VARF(5),
  166. * VARF(6),VARF(7),VARF(8),VARF(9),
  167. * VARF(10),VARF(11),VARF(12),VARF(13),
  168. * VARF(14),VARF(15),VARF(16),VARF(17),
  169. * VARF(18),VARF(19),VARF(20))
  170. C===================================================================
  171. C
  172. C CONVERSION ENTIER
  173. C
  174. VARF(1)=DBLE(IFR1)
  175. C
  176. C REMPLISSAGE DU TABLEAU CONT (PARTIE ELASTIQUE)
  177. C
  178. CISA=0.5D0*XMAT(1)/(1.D0+XMAT(2))
  179. C SIGF(1)=SIG0(1) + DSIGT(1)
  180. SIGF(1)=SIG0(1) + XMAT(1)*XCAR(1)*DEPST(1)
  181. C SIGF(2)=SIG0(2) + DSIGT(2)
  182. SIGF(2)=SIG0(2) + CISA*XCAR(3)*DEPST(2)
  183. C
  184. SIGF(3)=XMOMY
  185. C
  186. C LES NOUVAUX VALEURS DES VARIABLES INTERNES
  187. C
  188. C ----- NOTE: VARF(2) CONTAINS THE CURRENT (EI) THAT WILL BE USED
  189. C TO COMPUTE THE TANGENT STIFFNESS MATRIX
  190. C
  191. RETURN
  192. END
  193.  
  194.  
  195.  
  196.  

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