Télécharger flutra.eso

Retour à la liste

Numérotation des lignes :

  1. C FLUTRA SOURCE CB215821 16/04/21 21:16:55 8920
  2. SUBROUTINE FLUTRA(SIG0,NSTRS,DEPST,VAR0,NVARI,XMAT,NMATT,
  3. * IFOUR,TDT,IB,IGAU,TETA1,TETA2,ITHER,NMAT,
  4. * SIGF,VARF,CMAT,DEPSI,SIGI,DEP,EPSILO,EPSEQU,
  5. * EPSFLU,DSIGT,SIGEF,AN,BN,CN,KERRE)
  6. C
  7. C___________________________________________________________________
  8. C
  9. C MODELE N. REVIRON
  10. C
  11. C DEFINITION DES DIFFERENTES VARIABLES ENTREES, SORTIES,
  12. C PARAMETRES
  13. C
  14. C___________________________________________________________________
  15. C
  16. C NSTRS = TAILLE DU VECTEUR DES CONTRAINTES
  17. C
  18. C ENTREES 3D
  19. C
  20. C SIG0(NSTRS) = VECTEUR DES CONTRAINTES INITIALES APPARENTES
  21. C 1 : XX
  22. C 2 : YY
  23. C 3 : ZZ
  24. C 4 : XY
  25. C 5 : ZX
  26. C 6 : YZ
  27. C
  28. C SIGI(NSTRS) = VECTEUR DES CONTRAINTES INITIALES EFFECTIVES
  29. C (SE11,SE22,SE12)
  30. C idem
  31. C
  32. C DEPST(NSTRS)= VECTEUR DES INCREMENTS DE DEFORMATIONS TOTALES
  33. C idem
  34. C
  35. C ENTREES 2D CONTRAINTES PLANES
  36. C
  37. C SIG0(NSTRS) = VECTEUR DES CONTRAINTES INITIALES APPARENTES
  38. C 1 : XX
  39. C 2 : YY
  40. C (3 : ZZ)
  41. C 4 : XY
  42. C
  43. C SIGI(NSTRS) = VECTEUR DES CONTRAINTES INITIALES EFFECTIVES
  44. C (SE11,SE22,SE12)
  45. C idem
  46. C
  47. C DEPST(NSTRS)= VECTEUR DES INCREMENTS DE DEFORMATIONS TOTALES
  48. C 1 : XX
  49. C 2 : YY
  50. C (3 : ZZ)
  51. C 4 : XY
  52. C
  53. C DEPSI (UTILSER DANS LES SUBROUTINES)
  54. C 1 : XX
  55. C 2 : YY
  56. C 3 : XY
  57. C (4 : ZZ)
  58. C
  59. C ENTREES 2D DEFORMATIONS PLANES
  60. C
  61. C SIG0(NSTRS) = VECTEUR DES CONTRAINTES INITIALES APPARENTES
  62. C 1 : XX
  63. C 2 : YY
  64. C 3 : ZZ
  65. C 4 : XY
  66. C
  67. C SIGI(NSTRS) = VECTEUR DES CONTRAINTES INITIALES EFFECTIVES
  68. C (SE11,SE22,SE12)
  69. C idem
  70. C
  71. C DEPST(NSTRS)= VECTEUR DES INCREMENTS DE DEFORMATIONS TOTALES
  72. C 1 : XX
  73. C 2 : YY
  74. C (3 : ZZ = 0)
  75. C 4 : XY
  76. C
  77. C ENTREES AXISYMETRIE
  78. C
  79. C SIG0(NSTRS) = VECTEUR DES CONTRAINTES INITIALES APPARENTES
  80. C 1 : RR
  81. C 2 : ZZ
  82. C 3 : TT
  83. C 4 : RZ
  84. C
  85. C SIGI(NSTRS) = VECTEUR DES CONTRAINTES INITIALES EFFECTIVES
  86. C idem
  87. C
  88. C DEPST(NSTRS)= VECTEUR DES INCREMENTS DE DEFORMATIONS TOTALES
  89. C idem
  90. C___________________________________________________________
  91. C
  92. C VAR0(NVARI) = VECTEUR DES VARIABLES INTERNES INITIALES
  93. C___________________________________________________________
  94. C
  95. C VARIABLES INTERNES MODELE MECANIQUE ISOTROPE
  96. C
  97. C VAR0(1) = EPSTIM
  98. C VAR0(2) = DINI
  99. C VAR0(3) = EF11
  100. C VAR0(4) = EF22
  101. C VAR0(5) = EF33
  102. C VAR0(6) = EF12
  103. C VAR0(7) = EF31
  104. C VAR0(8) = EF23
  105. C VAR0(9) = SE11
  106. C VAR0(10) = SE22
  107. C VAR0(11) = SE12
  108. C__________________________________________________________
  109. C
  110. C NVARI = TAILLE DU VECTEUR DE VARIABLES INTERNES
  111. C__________________________________________________________
  112. C
  113. C XMAT(NVARI) = VECTEUR DES PARAMETRES DU MATERIAU
  114. C
  115. C PARAMETRES STANDARD
  116. C XMAT(1) = YOUN : MODULE D YOUNG
  117. C XMAT(2) = XNU : COEFFICIENT DE POISSON ELASTIQUE
  118. C XMAT(3) = RHO : MASSE VOLUMIQUE
  119. C XMAT(4) = ALPH : COEFFICIENT DE DILATATION THERMIQUE
  120. C
  121. C PARAMETRES MODELE MECANIQUE
  122. C XMAT(5) = KTR0 : KTR0 Mazars
  123. C XMAT(6) = ATRA : ATRA Mazars
  124. C XMAT(7) = BTRA : BTRA Mazars
  125. C XMAT(11) = A : COEF De Schutter Module d YOUNG
  126. C XMAT(12) = B : COEF De Schutter fc : non utilise
  127. C XMAT(13) = C : COEF De Schutter ft
  128. C XMAT(14) = D : COEF De Schutter Gft
  129. C___________________________________________________
  130. C
  131. C NMATT = TAILLE DU VECTEUR DES PARAMETRES
  132. C MELE = NUMERO DE L ELEMENT FINI UTILISE
  133. C XE = VECTEUR DES COORDONNEES DE L ELEMENT
  134. C NBNN = TAILLE DU VECTEUR XE
  135. C IFOUR = INDICATEUR DU TYPE DE CALCUL
  136. C -3 : DEFORMATION PLANE GENERALISEE
  137. C -2 : CONTRAINTES PLANES
  138. C -1 : DEFORMATION PLANE
  139. C 0 : AXISYMETRIQUE
  140. C 1 : SERIE DE FOURIER
  141. C 2 : TRIDIMENSIONNEL
  142. C MFR = NUMERO DE FORMULATION
  143. C 1 : ELEMENTS MASSIFS
  144. C AUTRE: VOIR nummfr.eso
  145. C ITTHER = INDICATEUR DE CHARGEMENT THERMIQUE
  146. C 0 : PAS DE CHARGEMENT THERMIQUE
  147. C 1 : PRESENCE DE CHARGEMENT THERMIQUE
  148. C TINI = TEMPERATURE EN DEBUT DE PAS
  149. C TFIN = TEMPERATURE EN FIN DE PAS
  150. C TDT = PAS DE TEMPS
  151. C IB = NUMERO DE L ELEMENT CALCULE
  152. C IGAU = NUMERO DU POINT DE GAUSS
  153. C DEP(6,6) = MATRICE ELASTO-PLASTIQUE TANGENTE
  154. C DINI = ENDOMMAGEMENT EN DEBUT DE PAS
  155. C D = ENDOMMAGEMENT EN FIN DE PAS
  156. C SIG0 = VECTEUR DES CONTRAINTES INITIALES APPARENTES
  157. C SIGF = VECTEUR DES CONTRAINTES FINALES APPARENTES
  158. C SIGI = VECTEUR DES CONTRAINTES INITIALES EFFECTIVES
  159. C EPSILO = VECTEUR DES DEFORMATIONS ELASTIQUES FIN PAS
  160. C____________________________________________________
  161. C
  162. C VARIABLES EN SORTIES
  163. C____________________________________________________
  164. C
  165. C SIGF(NSTRS) = VECTEUR DES CONTRAINTES FINALES
  166. C REMARQUE :
  167. C SIGF est dans un premier temps les contraintes effectives finales
  168. C puis les contraintes apparentes finales par multiplication par D
  169. C
  170. C DEPSFT = VECTEUR INCREMENT DES DEFORMATIONS DE FLUAGE
  171. C (EF11,EF22,EF33,EF12,EF31,EF23)
  172. C VARF(NVARI) = VECTEUR DES VARIABLES INTERNES FINALES
  173. C KERRE = INDICATEUR D ERREUR
  174. C 0 : PAS D ERREUR
  175. C
  176. C____________________________________________________
  177. C____________________________________________________
  178. C____________________________________________________
  179. C____________________________________________________
  180. C
  181. C
  182. IMPLICIT INTEGER(I-N)
  183. IMPLICIT REAL*8(A-H,O-Z)
  184. DIMENSION XMAT(NMATT),CMAT(NMAT)
  185. DIMENSION VAR0(NVARI),VARF(NVARI)
  186. DIMENSION DEPST(NSTRS),DEPSI(NSTRS)
  187. DIMENSION SIG0(NSTRS),SIGF(NSTRS),SIGI(NSTRS)
  188. DIMENSION DEP(NSTRS,NSTRS)
  189. DIMENSION EPSILO(NSTRS)
  190. DIMENSION EPSEQU(NSTRS),EPSFLU(NSTRS)
  191. DIMENSION DSIGT(NSTRS),SIGEF(NSTRS)
  192. DIMENSION AN(NSTRS),BN(NSTRS,NSTRS),CN(NSTRS,NSTRS)
  193. C
  194. C TEST SUR L OPTION DE CALCUL
  195. C LE CALCUL N EST POSSIBLE QUE DANS LE CAS
  196. C CONTRAINTES PLANES : IFOUR = -2
  197. C DEFORMATION PLANE : IFOUR = -1
  198. C AXISYMETRIQUE : IFOUR = 0
  199. C TRIDIMENSIONEL : IFOUR = 2
  200. C
  201. IF ((IFOUR.NE.-2).AND.(IFOUR.NE.-1).AND.(IFOUR.NE.0)
  202. $ .AND.(IFOUR.NE.2)) THEN
  203. WRITE(*,*) 'L OPTION DE CALCUL N EST PAS IMPLEMENTE'
  204. WRITE(*,*) 'IL FAUT CHOISIR PARMIS LES OPTIONS SUIVANTES :'
  205. WRITE(*,*) ' CONTRAINTES PLANES'
  206. WRITE(*,*) ' DEFORMATIONS PLANES'
  207. WRITE(*,*) ' AXISYMETRIQUE'
  208. WRITE(*,*) ' TRIDIMENSIONEL'
  209. KERRE = 99
  210. RETURN
  211. ENDIF
  212. C
  213. C INITIALISATION DES MATRICES
  214. C
  215. XPRECIS = 1.D-6
  216. CALL ZERO(CMAT,NMAT,1)
  217. CALL ZERO(SIGF,NSTRS,1)
  218. C
  219. C RECUPERATION DES VALEURS DES PARAMETRES MATERIAU
  220. C
  221. C
  222. C REMPLISSAGE D UN NOUVEAU VECTEUR DE MATERIAU
  223. C
  224. DO 300 I =1, NMATT
  225. CMAT(I)= XMAT(I)
  226. 300 CONTINUE
  227. C
  228. C AJOUT DU TEMPS A LA FIN
  229. C
  230. CMAT(NMAT) = TDT
  231. C
  232. C GESTION DES TEMPERATURES
  233. C
  234. CMAT(29) = TETA1
  235. CMAT(30) = TETA2
  236. C
  237. C RECUPERATION DES VARIABLES INTERNES INITIALES
  238. C
  239. C REMPLISSAGE DES VARIABLES INTERNES FINALES
  240. C
  241. DO 200 I =1, NVARI
  242. VARF(I)= VAR0(I)
  243. 200 CONTINUE
  244. C
  245. C_______________________________________________________________________
  246. C
  247. C REDIMENSIONNEMENT DES VECTEURS DANS LE CAS DE CONTRAINTES PLANES
  248. C_______________________________________________________________________
  249. C
  250. C
  251. C
  252. C CHANGEMENT DE NOM POUR TRAITEMENT
  253. C
  254. DO 100 I =1, NSTRS
  255. DEPSI(I)= DEPST(I)
  256. 100 CONTINUE
  257. C
  258. IPLAN = 0
  259. IF ((NSTRS.EQ.4).AND.(IFOUR.EQ.-2)) THEN
  260. IPLAN = 1
  261. DEPSI(3) = DEPST(4)
  262. DEPSI(4) = DEPST(3)
  263. IF (ABS(SIG0(3)).GE.1.D-3) THEN
  264. WRITE(*,*)'PROBLEME AVEC LES CONTRAINTES PLANES'
  265. KERRE=99
  266. RETURN
  267. ENDIF
  268. ENDIF
  269. C_______________________________________________________________________
  270. C
  271. C MODELE MECANIQUE
  272. C_______________________________________________________________________
  273. C
  274. CALL CAMECA(NSTRS,DEPSI,VAR0,NVARI,CMAT,NMAT,
  275. * IFOUR,IB,IGAU,SIGF,VARF,SIGI,DEP,
  276. * EPSILO,EPSEQU,EPSFLU,DSIGT,SIGEF,
  277. * AN,BN,CN,KERRE)
  278. C
  279. C
  280. C_______________________________________________________________________
  281. C
  282. C REDIMENSIONNEMENT DES VECTEURS DANS LE CAS DE CONTRAINTES PLANES
  283. C_______________________________________________________________________
  284. C
  285. IF (IPLAN.EQ.1) THEN
  286. temp = SIGF(4)
  287. SIGF(4)=SIGF(3)
  288. SIGF(3)=temp
  289. ENDIF
  290. C
  291. END
  292.  
  293.  
  294.  
  295.  
  296.  

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