Télécharger creep.eso

Retour à la liste

Numérotation des lignes :

  1. C CREEP SOURCE CB215821 16/04/21 21:16:08 8920
  2. SUBROUTINE CREEP (DECRA,DESWA,STATEV,serd,EC0,ESW0,P,QTILD,
  3. & TEMP,DTEMP,PREDEF,DPRED,TIME,DTIME,
  4. & CMNAME,leximp,LEND,COORDS,NSTATV,NOEL,NPT,
  5. & layer,kspt,KSTEP,KINC)
  6. IMPLICIT INTEGER(I-N)
  7. IMPLICIT REAL*8(A-H,O-Z)
  8. -INC CCREEL
  9. C-----------------------------------------------------------------------
  10. C
  11. C DESCRIPTION FONCTIONNELLE :
  12. C -------------------------
  13. C
  14. C Version du module externe CREEP integree au code pour eviter tout
  15. C probleme d'edition de lien et pour permettre le passage de tests
  16. C de non regression.
  17. C
  18. C En utilisation normale, ce module doit etre surcharge par celui
  19. C de meme nom fourni par l'utilisateur.
  20. C
  21. C Si cette version de CREEP est appelee, cela signifie qu'il n'y a
  22. C pas eu surcharge par le module utilisateur
  23. C ATTENTION : si l'utilisateur oublie de surcharger et affecte a sa
  24. C loi le meme numero que l'un de ceux pris en charge ci-apres, il y
  25. C a risque d'erreur non controlee a l'execution.
  26. C
  27. C CONVENTION : on autorise la programmation de l'utilisateur a se
  28. C servir de KSTEP comme code retour, de la meme maniere que l'on
  29. C s'en sert ci-dessous. Regles a respecter :
  30. C 1. Pas d'initialisation superflue de KSTEP en entrant dans CREEP.
  31. C KSTEP est initialise a 1 par WCREE0, avant l'appel a CREEP.
  32. C 2. En cas d'erreur, KSTEP est affecte d'une valeur differente de 1
  33. C
  34. C-----------------------------------------------------------------------
  35. C
  36. C DESCRIPTION DE L'INTERFACE :
  37. C --------------------------
  38. C
  39. C OUT : DECRA : REAL*8(5)
  40. C DECRA(1) : increment scalaire de deformation de
  41. C fluage, au debut du pas si LEND=0, a la fin du
  42. C pas si LEND=1.
  43. C Les composantes 2 a 5 ne sont pas exploitees
  44. C dans le cas d'une adherence a CAST3M.
  45. C
  46. C OUT : DESWA : REAL*8(5)
  47. C DESWA(1) : increment scalaire de deformation de
  48. C gonflement, au debut du pas si LEND=0, a la fin
  49. C du pas si LEND=1.
  50. C Les composantes 2 a 5 ne sont pas exploitees
  51. C dans le cas d'une adherence a CAST3M.
  52. C
  53. C Remarque : la routine CREEP evalue des VITESSES
  54. C de deformation (fluage et gonflement) au debut
  55. C ou a la fin du pas, suivant la valeur de LEND.
  56. C Les increments DECRA(1) et DESWA(1) sont ensuite
  57. C determines par les produits de ces vitesses par
  58. C le pas de temps DTIME.
  59. C
  60. C IN/OUT : STATEV : REAL*8(*), variables internes supplementaires
  61. C Il s'agit des eventuelles variables internes
  62. C s'ajoutant aux 4 pre-definies 'EC0 ', 'ESW0',
  63. C 'P ' et 'QTLD'.
  64. C Valeurs au debut du pas si LEND=0.
  65. C Peuvent etre mises a jour si LEND=1 : valeurs a
  66. C la fin du pas.
  67. C
  68. C IN/OUT : serd : REAL*8, puissance volumique de deformation
  69. C inelastique (produit contracte du tenseur des
  70. C contraintes et du tenseur des vitesses de
  71. C deformation inelastique).
  72. C Entree/sortie facultative de la routine CREEP,
  73. C non exploitee par CAST3M pour l'instant.
  74. C
  75. C IN : EC0 : REAL*8, deformation de fluage cumulee.
  76. C 1ere variable interne parmi les 4 pre-definies.
  77. C Valeur au debut du pas.
  78. C
  79. C IN : ESW0 : REAL*8, deformation de gonflement cumulee.
  80. C 2eme variable interne parmi les 4 pre-definies.
  81. C Valeur au debut du pas.
  82. C
  83. C IN : P : REAL*8, 1er invariant du tenseur des contraintes
  84. C (1/3 trace(SIGMA)).
  85. C 3eme variable interne parmi les 4 pre-definies.
  86. C Valeur au debut du pas si LEND=0, a la fin du
  87. C pas si LEND=1.
  88. C
  89. C IN : QTILD : REAL*8, contrainte equivalente de Von Mises
  90. C (SQRT(3/2 S:S) ou S deviateur des contraintes).
  91. C 4eme variable interne parmi les 4 pre-definies.
  92. C Valeur au debut du pas si LEND=0, a la fin du
  93. C pas si LEND=1.
  94. C
  95. C IN : TEMP : REAL*8, temperature a la fin du pas.
  96. C DTEMP : REAL*8, increment de temperature sur le pas.
  97. C PREDEF : REAL*8(*), vecteur des parametres externes de la
  98. C loi de comportement, valeurs a la fin du pas.
  99. C DPRED : REAL*8(*), increments des parametres externes
  100. C sur le pas de temps.
  101. C
  102. C IN : TIME : REAL*8(2)
  103. C TIME(1) : duree cumulee des iterations internes
  104. C a la fin du pas en cours.
  105. C TIME(2) : instant intermediaire absolu a la fin
  106. C du pas en cours.
  107. C DTIME : REAL*8, valeur du pas de temps en cours.
  108. C
  109. C Remarque : l'instant absolu correspondant a la
  110. C precedente iteration interne convergee est
  111. C TIME(2)-DTIME.
  112. C
  113. C IN : CMNAME : CHARACTER*16, identifiant de la loi de
  114. C comportement.
  115. C On conserve le type 'chaine de caracteres' pour
  116. C l'identifiant de la loi, afin de preserver la
  117. C compatibilite avec ABAQUS.
  118. C Dans le cas d'une adherence a CAST3M, la loi est
  119. C identifiee par le numero qui lui a ete attribue:
  120. C argument 'NUME_LOI' de l'operateur MODE.
  121. C Par convention, ce numero est encode dans les 4
  122. C derniers caracteres de la chaine, et doit etre
  123. C recupere dans CREEP par une instruction du type
  124. C K4ILOI = CMNAME(13:16)
  125. C avec K4ILOI variable locale de type CHARACTER*4.
  126. C
  127. C IN : leximp : INTEGER, indicateur d'integration par un schema
  128. C Entree non active dans le cas d'une adherence
  129. C a CAST3M.
  130. C
  131. C IN : LEND : INTEGER, indicateur de debut/fin de pas.
  132. C LEND=0 : les entrees P, QTILD sont definies au
  133. C debut du pas ; les sorties DECRA(1), DESWA(1)
  134. C sont calculees au debut du pas.
  135. C LEND=1 : les entrees P, QTILD sont definies a la
  136. C fin du pas ; les sorties DECRA(1), DESWA(1) sont
  137. C calculees a la fin du pas.
  138. C
  139. C IN : COORDS : REAL*8(*), coordonnees cartesiennes du point
  140. C d'integration courant.
  141. C
  142. C IN NSTATV : INTEGER, nombre de variables internes
  143. C supplementaires (en plus des 4 pre-definies).
  144. C
  145. C IN : NOEL : INTEGER, numero de l'element courant.
  146. C NPT : INTEGER, numero du point d'integration courant.
  147. C
  148. C IN : layer : INTEGER, numero de couche pour des coques
  149. C composites ou des solides multi-couches.
  150. C kspt : INTEGER, numero de section dans la couche
  151. C courante.
  152. C Entrees non actives dans le cas d'une adherence
  153. C a CAST3M.
  154. C
  155. C IN/OUT : KSTEP : INTEGER
  156. C Entree n'ayant pas de sens dans le cas d'une
  157. C adherence a CAST3M.
  158. C On autorise la programmation de l'utilisateur a
  159. C se servir de KSTEP comme code retour de CREEP
  160. C (cf. description fonctionnelle ci-dessus).
  161. C
  162. C IN : KINC : INTEGER, compteur d'iterations internes.
  163. C Incrementation geree par l'appelant CCREEP.
  164. C
  165. C-----------------------------------------------------------------------
  166. C
  167. C MODELES PRIS EN CHARGE par le module CREEP 'bouchon' :
  168. C ----------------------------------------------------
  169. C
  170. C ILOI = 21 : Modele 'FLUAGE' 'POLYNOMIAL' de CAST3M, INPLAS=21
  171. C ILOI = 23 RESTRICTION A LA FORMULATION MASSIVE
  172. C Le modele ILOI=23 comprend deux variables internes
  173. C supplementaires artificielles :
  174. C 'PSUP' <=> 'P ' et 'QSUP' <=> 'QTLD'
  175. C
  176. C ILOI = 22 : Modele 'FLUAGE' 'POLYNOMIAL' de CAST3M, INPLAS=21
  177. C ILOI = 24 RESTRICTION A LA FORMULATION MASSIVE
  178. C Instanciation des composantes de materiau au fur et a
  179. C mesure des iterations internes (uniquement les coef.
  180. C du polynome de la loi de fluage).
  181. C Le modele comprend cette fois des parametres externes.
  182. C Les composantes de materiau sont evaluees par les
  183. C memes fonctions externes que celles utilisees par
  184. C l'operateur VARI.
  185. C MODELE 'MAQUETTE'
  186. C Le modele ILOI=24 comprend deux variables internes
  187. C supplementaires artificielles :
  188. C 'PSUP' <=> 'P ' et 'QSUP' <=> 'QTLD'
  189. C
  190. C-----------------------------------------------------------------------
  191. C Arguments de l'interface
  192. C
  193. CHARACTER*16 CMNAME
  194. C
  195. INTEGER leximp, LEND, NSTATV, NOEL, NPT, layer, kspt,
  196. & KSTEP, KINC
  197. C
  198. REAL*8 DECRA(5), DESWA(5), STATEV(*),
  199. & serd, EC0, ESW0, P, QTILD,
  200. & TEMP, DTEMP, PREDEF(*), DPRED(*), TIME(2), DTIME,
  201. & COORDS(*)
  202. C
  203. C Variables locales
  204. C
  205. C Tableaux de travail
  206. C
  207. REAL*8 VALPAR(5)
  208. C
  209. C Variables scalaires
  210. C
  211. INTEGER NBPAR
  212. C
  213. REAL*8 TFIS, TUO2, FACF, DSIU, DGRA,
  214. & XAF0, XAF1, XAF2, XAF3, XAF4, XAF5, XAF6,
  215. & EC0PT
  216. C
  217. CHARACTER*4 K4ILOI
  218. C
  219. C
  220. C
  221. C------------------- Debut du code executable --------------------------
  222. C
  223. C
  224. K4ILOI = CMNAME(13:16)
  225. C
  226. C=======================================================================
  227. C ILOI = 21 :
  228. C Modele 'FLUAGE' 'POLYNOMIAL' de CAST3M, INPLAS=21
  229. C RESTRICTION A LA FORMULATION MASSIVE
  230. C=======================================================================
  231. C
  232. IF (K4ILOI.EQ.' 21') THEN
  233. C
  234. C........Affectation des parametres pour evaluation des coefficients
  235. C du polynome de la loi de fluage
  236. C
  237. TFIS = 1.5D+19
  238. TUO2 = 1.5D+3
  239. FACF = 1.D0
  240. DSIU = 0.95D0
  241. DGRA = 8.0D-6
  242.  
  243. ** pour eviter les overflow
  244. qtild=MAX(min(1D65,QTiLD),-1d65)
  245.  
  246.  
  247. C
  248. C........Evaluation des coefficients du polynome de la loi de fluage
  249. C
  250. XAF0 = 0.D0 * FACF
  251. XAF1 = ( ( EXP(0.183D0*(100.D0-DSIU)) ) *
  252. & ( (0.2031D-13+(0.67861D-32*TFIS)) / (DGRA*DGRA) ) *
  253. & ( EXP(-45294.4D0/TUO2) ) )
  254. & * FACF
  255. XAF2 = 1.D0
  256. XAF3 = ( ( EXP(0.241D0*(100.D0-DSIU)) ) *
  257. & ( 0.1524D-25+(0.50825D-44*TFIS) ) *
  258. & ( EXP(-66431.8D0/TUO2) ) )
  259. & * FACF
  260. XAF4 = 4.5D0
  261. XAF5 = ( (1.447D-34*TFIS) * (EXP(-3624.D0/TUO2)) ) * FACF
  262. XAF6 = 1.D0
  263. C
  264. C........Calcul de la vitesse de deformation de fluage et de l'increment
  265. C........Pas de terme de gonflement
  266. C
  267. EC0PT = XAF0 + XAF1 * (QTILD**XAF2) + XAF3 * (QTILD**XAF4)
  268. & + XAF5 * (QTILD**XAF6)
  269. DECRA(1) = EC0PT * DTIME
  270. C
  271. DESWA(1) = 0.D0
  272. C
  273. C=======================================================================
  274. C ILOI = 22 :
  275. C Modele 'FLUAGE' 'POLYNOMIAL' de CAST3M, INPLAS=21
  276. C RESTRICTION A LA FORMULATION MASSIVE
  277. C Instanciation des composantes de materiau au fur et a mesure des
  278. C iterations internes (uniquement les coef. du polynome de la loi
  279. C de fluage).
  280. C Le modele comprend cette fois des parametres externes.
  281. C Les composantes de materiau sont evaluees par les memes fonctions
  282. C externes que celles utilisees par l'operateur VARI.
  283. C MODELE 'MAQUETTE'
  284. C=======================================================================
  285. C
  286. ELSE IF (K4ILOI.EQ.' 22') THEN
  287. C
  288. IERUT = 0
  289. C
  290. C........Affectation des parametres et evaluation des coefficients
  291. C du polynome de la loi de fluage
  292. C
  293. IF ( LEND.EQ.0 ) THEN
  294. C
  295. NBPAR = 3
  296. VALPAR(1) = PREDEF(1)-DPRED(1)
  297. VALPAR(2) = PREDEF(2)-DPRED(2)
  298. VALPAR(3) = PREDEF(3)-DPRED(3)
  299. XAF0 = FVAF0(VALPAR,NBPAR,IERUT)
  300. IF (IERUT.GT.0) THEN
  301. KSTEP = -22000-IERUT
  302. RETURN
  303. ENDIF
  304. C
  305. XAF2 = FVAF2(VALPAR,NBPAR,IERUT)
  306. IF (IERUT.GT.0) THEN
  307. KSTEP = -22000-IERUT
  308. RETURN
  309. ENDIF
  310. C
  311. XAF4 = FVAF4(VALPAR,NBPAR,IERUT)
  312. IF (IERUT.GT.0) THEN
  313. KSTEP = -22000-IERUT
  314. RETURN
  315. ENDIF
  316. C
  317. XAF5 = FVAF5(VALPAR,NBPAR,IERUT)
  318. IF (IERUT.GT.0) THEN
  319. KSTEP = -22000-IERUT
  320. RETURN
  321. ENDIF
  322. C
  323. XAF6 = FVAF6(VALPAR,NBPAR,IERUT)
  324. IF (IERUT.GT.0) THEN
  325. KSTEP = -22000-IERUT
  326. RETURN
  327. ENDIF
  328. C
  329. NBPAR = 4
  330. VALPAR(4) = PREDEF(4)-DPRED(4)
  331. XAF3 = FVAF3(VALPAR,NBPAR,IERUT)
  332. IF (IERUT.GT.0) THEN
  333. KSTEP = -22000-IERUT
  334. RETURN
  335. ENDIF
  336. C
  337. NBPAR = 5
  338. VALPAR(5) = PREDEF(5)-DPRED(5)
  339. XAF1 = FVAF1(VALPAR,NBPAR,IERUT)
  340. IF (IERUT.GT.0) THEN
  341. KSTEP = -22000-IERUT
  342. RETURN
  343. ENDIF
  344. C
  345. ELSE
  346. C
  347. NBPAR = 3
  348. XAF0 = FVAF0(PREDEF,NBPAR,IERUT)
  349. IF (IERUT.GT.0) THEN
  350. KSTEP = -22000-IERUT
  351. RETURN
  352. ENDIF
  353. C
  354. XAF2 = FVAF2(PREDEF,NBPAR,IERUT)
  355. IF (IERUT.GT.0) THEN
  356. KSTEP = -22000-IERUT
  357. RETURN
  358. ENDIF
  359. C
  360. XAF4 = FVAF4(PREDEF,NBPAR,IERUT)
  361. IF (IERUT.GT.0) THEN
  362. KSTEP = -22000-IERUT
  363. RETURN
  364. ENDIF
  365. C
  366. XAF5 = FVAF5(PREDEF,NBPAR,IERUT)
  367. IF (IERUT.GT.0) THEN
  368. KSTEP = -22000-IERUT
  369. RETURN
  370. ENDIF
  371. C
  372. XAF6 = FVAF6(PREDEF,NBPAR,IERUT)
  373. IF (IERUT.GT.0) THEN
  374. KSTEP = -22000-IERUT
  375. RETURN
  376. ENDIF
  377. C
  378. NBPAR = 4
  379. XAF3 = FVAF3(PREDEF,NBPAR,IERUT)
  380. IF (IERUT.GT.0) THEN
  381. KSTEP = -22000-IERUT
  382. RETURN
  383. ENDIF
  384. C
  385. NBPAR = 5
  386. XAF1 = FVAF1(PREDEF,NBPAR,IERUT)
  387. IF (IERUT.GT.0) THEN
  388. KSTEP = -22000-IERUT
  389. RETURN
  390. ENDIF
  391. C
  392. ENDIF
  393. C
  394. C........Calcul de la vitesse de deformation de fluage et de l'increment
  395. C........Pas de terme de gonflement
  396. C
  397. EC0PT = XAF0 + XAF1 * (QTILD**XAF2) + XAF3 * (QTILD**XAF4)
  398. & + XAF5 * (QTILD**XAF6)
  399. DECRA(1) = EC0PT * DTIME
  400. C
  401. DESWA(1) = 0.D0
  402. C
  403. C=======================================================================
  404. C ILOI = 23 :
  405. C Modele 'FLUAGE' 'POLYNOMIAL' de CAST3M, INPLAS=21
  406. C RESTRICTION A LA FORMULATION MASSIVE
  407. C Equivalent au modele ILOI=21
  408. C Variables internes supplementaires 'PSUP' et 'QSUP'
  409. C=======================================================================
  410. C
  411. ELSE IF (K4ILOI.EQ.' 23') THEN
  412. C
  413. C........Affectation des parametres pour evaluation des coefficients
  414. C du polynome de la loi de fluage
  415. C
  416. TFIS = 1.5D+19
  417. TUO2 = 1.5D+3
  418. FACF = 1.D0
  419. DSIU = 0.95D0
  420. DGRA = 8.0D-6
  421. C
  422. C........Evaluation des coefficients du polynome de la loi de fluage
  423. C
  424. XAF0 = 0.D0 * FACF
  425. = ( ( EXP(0.183D0*(100.D0-DSIU)) ) *
  426. & ( (0.2031D-13+(0.67861D-32*TFIS)) / (DGRA*DGRA) ) *
  427. & ( EXP(-45294.4D0/TUO2) ) )
  428. & * FACF
  429. XAF2 = 1.D0
  430. XAF3 = ( ( EXP(0.241D0*(100.D0-DSIU)) ) *
  431. & ( 0.1524D-25+(0.50825D-44*TFIS) ) *
  432. & ( EXP(-66431.8D0/TUO2) ) )
  433. & * FACF
  434. XAF4 = 4.5D0
  435. XAF5 = ( (1.447D-34*TFIS) * (EXP(-3624.D0/TUO2)) ) * FACF
  436. XAF6 = 1.D0
  437. C
  438. C........Calcul de la vitesse de deformation de fluage et de l'increment
  439. C........Pas de terme de gonflement
  440. C
  441. EC0PT = XAF0 + XAF1 * (QTILD**XAF2) + XAF3 * (QTILD**XAF4)
  442. & + XAF5 * (QTILD**XAF6)
  443. DECRA(1) = EC0PT * DTIME
  444. C
  445. DESWA(1) = 0.D0
  446. C
  447. C........MAJ des variables internes supplementaires le cas echeant
  448. C
  449. IF ( LEND.EQ.1 ) THEN
  450. STATEV(1) = P
  451. STATEV(2) = QTILD
  452. ENDIF
  453. C
  454. C=======================================================================
  455. C ILOI = 24 :
  456. C Modele 'FLUAGE' 'POLYNOMIAL' de CAST3M, INPLAS=21
  457. C RESTRICTION A LA FORMULATION MASSIVE
  458. C Instanciation des composantes de materiau au fur et a mesure des
  459. C iterations internes (uniquement les coef. du polynome de la loi
  460. C de fluage).
  461. C Le modele comprend cette fois des parametres externes.
  462. C Les composantes de materiau sont evaluees par les memes fonctions
  463. C externes que celles utilisees par l'operateur VARI.
  464. C MODELE 'MAQUETTE'
  465. C Equivalent au modele ILOI=22
  466. C Variables internes supplementaires 'PSUP' et 'QSUP'
  467. C=======================================================================
  468. C
  469. ELSE IF (K4ILOI.EQ.' 24') THEN
  470. C
  471. C........Affectation des parametres et evaluation des coefficients
  472. C du polynome de la loi de fluage
  473. C MAJ des variables internes supplementaires le cas echeant
  474. C
  475. IF ( LEND.EQ.0 ) THEN
  476. C
  477. NBPAR = 3
  478. VALPAR(1) = PREDEF(1)-DPRED(1)
  479. VALPAR(2) = PREDEF(2)-DPRED(2)
  480. VALPAR(3) = PREDEF(3)-DPRED(3)
  481. XAF0 = FVAF0(VALPAR,NBPAR,IERUT)
  482. IF (IERUT.GT.0) THEN
  483. KSTEP = -24000-IERUT
  484. RETURN
  485. ENDIF
  486. C
  487. XAF2 = FVAF2(VALPAR,NBPAR,IERUT)
  488. IF (IERUT.GT.0) THEN
  489. KSTEP = -24000-IERUT
  490. RETURN
  491. ENDIF
  492. C
  493. XAF4 = FVAF4(VALPAR,NBPAR,IERUT)
  494. IF (IERUT.GT.0) THEN
  495. KSTEP = -24000-IERUT
  496. RETURN
  497. ENDIF
  498. C
  499. XAF5 = FVAF5(VALPAR,NBPAR,IERUT)
  500. IF (IERUT.GT.0) THEN
  501. KSTEP = -24000-IERUT
  502. RETURN
  503. ENDIF
  504. C
  505. XAF6 = FVAF6(VALPAR,NBPAR,IERUT)
  506. IF (IERUT.GT.0) THEN
  507. KSTEP = -24000-IERUT
  508. RETURN
  509. ENDIF
  510. C
  511. NBPAR = 4
  512. VALPAR(4) = PREDEF(4)-DPRED(4)
  513. XAF3 = FVAF3(VALPAR,NBPAR,IERUT)
  514. IF (IERUT.GT.0) THEN
  515. KSTEP = -24000-IERUT
  516. RETURN
  517. ENDIF
  518. C
  519. NBPAR = 5
  520. VALPAR(5) = PREDEF(5)-DPRED(5)
  521. XAF1 = FVAF1(VALPAR,NBPAR,IERUT)
  522. IF (IERUT.GT.0) THEN
  523. KSTEP = -24000-IERUT
  524. RETURN
  525. ENDIF
  526. C
  527. ELSE
  528. C
  529. NBPAR = 3
  530. XAF0 = FVAF0(PREDEF,NBPAR,IERUT)
  531. IF (IERUT.GT.0) THEN
  532. KSTEP = -24000-IERUT
  533. RETURN
  534. ENDIF
  535. C
  536. XAF2 = FVAF2(PREDEF,NBPAR,IERUT)
  537. IF (IERUT.GT.0) THEN
  538. KSTEP = -24000-IERUT
  539. RETURN
  540. ENDIF
  541. C
  542. XAF4 = FVAF4(PREDEF,NBPAR,IERUT)
  543. IF (IERUT.GT.0) THEN
  544. KSTEP = -24000-IERUT
  545. RETURN
  546. ENDIF
  547. C
  548. XAF5 = FVAF5(PREDEF,NBPAR,IERUT)
  549. IF (IERUT.GT.0) THEN
  550. KSTEP = -24000-IERUT
  551. RETURN
  552. ENDIF
  553. C
  554. XAF6 = FVAF6(PREDEF,NBPAR,IERUT)
  555. IF (IERUT.GT.0) THEN
  556. KSTEP = -24000-IERUT
  557. RETURN
  558. ENDIF
  559. C
  560. NBPAR = 4
  561. XAF3 = FVAF3(PREDEF,NBPAR,IERUT)
  562. IF (IERUT.GT.0) THEN
  563. KSTEP = -24000-IERUT
  564. RETURN
  565. ENDIF
  566. C
  567. NBPAR = 5
  568. XAF1 = FVAF1(PREDEF,NBPAR,IERUT)
  569. IF (IERUT.GT.0) THEN
  570. KSTEP = -24000-IERUT
  571. RETURN
  572. ENDIF
  573. C
  574. STATEV(1) = P
  575. STATEV(2) = QTILD
  576. C
  577. ENDIF
  578. C
  579. C........Calcul de la vitesse de deformation de fluage et de l'increment
  580. C........Pas de terme de gonflement
  581. C
  582. EC0PT = XAF0 + XAF1 * (QTILD**XAF2) + XAF3 * (QTILD**XAF4)
  583. & + XAF5 * (QTILD**XAF6)
  584. DECRA(1) = EC0PT * DTIME
  585. C
  586. DESWA(1) = 0.D0
  587. C
  588. C=======================================================================
  589. C Modele externe non reconnu
  590. C=======================================================================
  591. C
  592. ELSE
  593. C
  594. KSTEP = 0
  595. C
  596. ENDIF
  597. C
  598. RETURN
  599. END
  600.  
  601.  
  602.  
  603.  
  604.  
  605.  
  606.  

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