Télécharger gauv1.eso

Retour à la liste

Numérotation des lignes :

  1. C GAUV1 SOURCE BP208322 17/03/01 21:17:30 9325
  2. SUBROUTINE GAUV1(DDAUX,CMATE,VALMAT,VALCAR,N2EL,N2PTEL,
  3. 1 MFR,IFOU,IB,IGAU,EPAIST,MELE,NPINT,NBGMAT,NELMAT,SECT,
  4. 2 LHOOK,TXR,XLOC,XGLOB,D1HOOK,ROTHOO,DDHOMU,CRIGI,
  5. 3 SIG0,NSTRS,DEPST,VAR0,XMAT,NCOMAT,XCAR,TRAC,NCOURB,
  6. 4 NBPGAU,DSIGT,SIGF,VARF,DEFP,KERRE)
  7. *_______________________________________________________________________
  8. *
  9. IMPLICIT INTEGER(I-N)
  10. IMPLICIT REAL*8(A-H,O-Z)
  11. -INC CCOPTIO
  12. *
  13. DIMENSION SIG0(*),DEPST(*),VAR0(*),XMAT(*),XCAR(*),TRAC(*),
  14. . SIGF(*),VARF(*) ,DEFP(*)
  15. DIMENSION VALMAT(*),VALCAR(*),TXR(IDIM,*),CRIGI(*),
  16. . DDAUX(LHOOK,*),DDHOMU(LHOOK,*),XLOC(3,3),
  17. . XGLOB(3,3),D1HOOK(LHOOK,*),ROTHOO(LHOOK,*)
  18. DIMENSION E(4),CONT(4),DCONT(4),NUMETA(2),EPSPLA(4),EPS0(2),
  19. . DSIGT(*),SIGPA(2),SIGPB(2),SIGNA(2),SIGNB(2),SIGSUP(2),SIGINF(2),
  20. . EPSPA(2),EPSPB(2),EPSNA(2),EPSNB(2),EPSUP(2),EPINF(2),
  21. . SIGND(2),SIGNC(2),SIGPC(2),SIGPD(2),EPSND(2),EPSNC(2),EPSPC(2),
  22. . EPSPD(2),SIGBET(2,8),EPSBET(2,8)
  23. CHARACTER*8 CMATE
  24. *
  25. * INITIALISATIONS
  26. *
  27. KERRE=0
  28. CALL ZDANUL(SIGBET,16)
  29. CALL ZDANUL(EPSBET,16)
  30. CALL ZDANUL(EPSPLA, 4)
  31. IPRECI=0
  32. IF(NCOURB.NE.10.AND.NCOURB.NE.18) THEN
  33. KERRE=31
  34. RETURN
  35. ENDIF
  36. IF(NCOURB.EQ.10) IECRAS=0
  37. IF(NCOURB.EQ.18) IECRAS=1
  38. *
  39. * ON RECUPERE LES DONNEES
  40. *
  41. * D'ABORD LES CARACTERISTIQUES
  42. *
  43. TORS = XCAR(1)
  44. XINY = XCAR(2)
  45. XINZ = XCAR(3)
  46. SB = XCAR(4)
  47. SCY = XCAR(5)
  48. SCZ = XCAR(6)
  49. VX = XCAR(7)
  50. VY = XCAR(8)
  51. VZ = XCAR(9)
  52. DX = XCAR(10)
  53. DY = XCAR(11)
  54. DZ = XCAR(12)
  55. *
  56. * ENSUITE LES MATERIAUX
  57. *
  58. YOUNG=XMAT(1)
  59. XNU =XMAT(2)
  60. IF(IECRAS.EQ.1) THEN
  61. * SELON Y
  62. SIGPA(1)=TRAC(11)
  63. SIGPB(1)=TRAC(13)
  64. SIGPC(1)=TRAC(15)
  65. SIGPD(1)=TRAC(17)
  66. EPSPA(1)=TRAC(12)
  67. EPSPB(1)=TRAC(14)
  68. EPSPC(1)=TRAC(16)
  69. EPSPD(1)=TRAC(18)
  70. SIGNA(1)=TRAC( 7)
  71. SIGNB(1)=TRAC( 5)
  72. SIGNC(1)=TRAC( 3)
  73. SIGND(1)=TRAC( 1)
  74. EPSNA(1)=TRAC( 8)
  75. EPSNB(1)=TRAC( 6)
  76. EPSNC(1)=TRAC( 4)
  77. EPSND(1)=TRAC( 2)
  78. * SELON Z
  79. SIGPA(2)=TRAC(29)
  80. SIGPB(2)=TRAC(31)
  81. SIGPC(2)=TRAC(33)
  82. SIGPD(2)=TRAC(35)
  83. EPSPA(2)=TRAC(30)
  84. EPSPB(2)=TRAC(32)
  85. EPSPC(2)=TRAC(34)
  86. EPSPD(2)=TRAC(36)
  87. SIGNA(2)=TRAC(25)
  88. SIGNB(2)=TRAC(23)
  89. SIGNC(2)=TRAC(21)
  90. SIGND(2)=TRAC(19)
  91. EPSNA(2)=TRAC(26)
  92. EPSNB(2)=TRAC(24)
  93. EPSNC(2)=TRAC(22)
  94. EPSND(2)=TRAC(20)
  95. ELSE
  96. * SELON Y
  97. SIGPA(1)=TRAC( 7)
  98. SIGPB(1)=TRAC( 9)
  99. EPSPA(1)=TRAC( 8)
  100. EPSPB(1)=TRAC(10)
  101. SIGNA(1)=TRAC( 3)
  102. SIGNB(1)=TRAC( 1)
  103. EPSNA(1)=TRAC( 4)
  104. EPSNB(1)=TRAC( 2)
  105. * SELON Z
  106. SIGPA(2)=TRAC(17)
  107. SIGPB(2)=TRAC(19)
  108. EPSPA(2)=TRAC(18)
  109. EPSPB(2)=TRAC(20)
  110. SIGNA(2)=TRAC(13)
  111. SIGNB(2)=TRAC(11)
  112. EPSNA(2)=TRAC(14)
  113. EPSNB(2)=TRAC(12)
  114. ENDIF
  115. SIGTO =XMAT(6)
  116. SIGCO =XMAT(7)
  117. *
  118. * PUIS LES EFFORTS
  119. *
  120. * EFFX = SIG0(1)
  121. * EFFY = SIG0(2)
  122. * EFFY = SIG0(3)
  123. * XMOMX= SIG0(4)
  124. * XMOMY= SIG0(5)
  125. * XMOMZ= SIG0(6)
  126. *
  127. * ENFIN LES VARIABLES INTERNES
  128. *
  129. EPSEQ =VAR0(1)
  130. EPS0(1) =VAR0(2)
  131. SIGINF(1)=VAR0(3)
  132. EPINF(1) =VAR0(4)
  133. SIGSUP(1)=VAR0(5)
  134. EPSUP(1) =VAR0(6)
  135. NUMETA(1)=INT(VAR0(7))
  136. EPS0(2) =VAR0(8)
  137. SIGINF(2)=VAR0(9)
  138. EPINF(2) =VAR0(10)
  139. SIGSUP(2)=VAR0(11)
  140. EPSUP(2) =VAR0(12)
  141. NUMETA(2)=INT(VAR0(13))
  142. *
  143. * CALCUL DE L INCREMENT DE CONTRAINTE
  144. *
  145. CALL CALSIG(DEPST,DDAUX,NSTRS,CMATE,VALMAT,VALCAR,
  146. 1 N2EL,N2PTEL,MFR,IFOU,IB,IGAU,EPAIST,NBPGAU,
  147. 2 MELE,NPINT,NBGMAT,NELMAT,SECT,LHOOK,TXR,XLOC,
  148. 3 XGLOB,D1HOOK,ROTHOO,DDHOMU,CRIGI,DSIGT,IRTD)
  149. *
  150. IF(IRTD.NE.1) THEN
  151. KERRE=69
  152. GOTO 998
  153. ENDIF
  154. *
  155. *
  156. *
  157. * CALCUL DES CONTRAINTES INITIALES ET DES INCREMENTS
  158. *
  159. DIVX=XCAR(10)/XCAR(1)
  160. DIVY=XCAR(11)/XCAR(2)
  161. DIVZ=XCAR(12)/XCAR(3)
  162. IF(DIVX.EQ.0.) DIVX=1.D-10/SQRT(XCAR(1)*XCAR(4))
  163. IF(DIVY.EQ.0.) DIVY=1.D-10/SQRT(XCAR(2)*XCAR(4))
  164. IF(DIVZ.EQ.0.) DIVZ=1.D-10/SQRT(XCAR(3)*XCAR(4))
  165. CONT(1) = SIG0(5) * DIVY
  166. CONT(2) = SIG0(6) * DIVZ
  167. CONT(3) = SIG0(4) * DIVX
  168. CONT(4) = SIG0(1) / XCAR(4)
  169. DCONT(1) = DSIGT(5) * DIVY
  170. DCONT(2) = DSIGT(6) * DIVZ
  171. DCONT(3) = DSIGT(4) * DIVX
  172. DCONT(4) = DSIGT(1) / XCAR(4)
  173. *
  174. E(1)= YOUNG
  175. E(2)= YOUNG
  176. E(3)= YOUNG/2.D0/(1.D0+XNU)
  177. E(4)= YOUNG
  178. *
  179. CALL CASTET(SIGPA,SIGPB,SIGNA,SIGNB,SIGTO,SIGCO,EPSPA,EPSPB,
  180. . EPSNA,EPSNB,SIGSUP,SIGINF,EPSUP,EPINF,EPS0,E,CONT,DCONT,
  181. . NUMETA,EPSEQ,IECRAS,SIGND,EPSND,SIGNC,EPSNC,SIGPC,EPSPC,
  182. . SIGPD,EPSPD,SIGBET,EPSBET,IPRECI,IIMPI,KERRE)
  183. *
  184. IF(KERRE.NE.0) RETURN
  185. *
  186. * ON REMET LES EFFORTS
  187. *
  188. SIGF(1) = DCONT(4) * XCAR(4)
  189. SIGF(2) = SIG0(2)+DSIGT(2)
  190. SIGF(3) = SIG0(3)+DSIGT(3)
  191. SIGF(4) = DCONT(3)/ DIVX
  192. SIGF(5) = DCONT(1)/ DIVY
  193. SIGF(6) = DCONT(2)/ DIVZ
  194. *
  195. * PUIS LES VARIABLES INTERNES
  196. *
  197. VARF(1) = EPSEQ
  198. VARF(2) = EPS0(1)
  199. VARF(3) = SIGINF(1)
  200. VARF(4) = EPINF(1)
  201. VARF(5) = SIGSUP(1)
  202. VARF(6) = EPSUP(1)
  203. VARF(7) = REAL(NUMETA(1))
  204. VARF(8) = EPS0(2)
  205. VARF(9) = SIGINF(2)
  206. VARF(10)= EPINF(2)
  207. VARF(11)= SIGSUP(2)
  208. VARF(12)= EPSUP(2)
  209. VARF(13)= REAL(NUMETA(2))
  210. *
  211. * ENFIN LES DEFORMATIONS PLASTIQUES
  212. *
  213. DEFP(1) = EPSPLA(4)
  214. DEFP(2) = 0.D0
  215. DEFP(3) = 0.D0
  216. DEFP(4) = EPSPLA(3) / DX
  217. DEFP(5) = EPSPLA(1) / DY
  218. DEFP(6) = EPSPLA(2) / DZ
  219. *
  220. 998 RETURN
  221. END
  222.  
  223.  
  224.  

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