Télécharger curgcs.eso

Retour à la liste

Numérotation des lignes :

  1. C CURGCS SOURCE PV 17/12/08 21:17:10 9660
  2. SUBROUTINE CURGCS(WRK52,WRK53,WRK54,WRK4,NSTRSS,
  3. 1 IFOUR,DDT,IB,IGAU,HCAR,IVIS,T0,TF)
  4. C---------------------------------------------------------------------
  5. C PLASTICITE MODELE BETON
  6. C
  7. C ENTREES
  8. C SIG0(NSTRS) = CONTRAINTES INITIALES
  9. C NSTRS = NOMBRE DE CONTRAINTES
  10. C DEPST(NSTRS) = INCREMENT DE DEFORMATIONS TOTALES
  11. C VAR0(NVARI) = VARIABLES INTERNES DEBUT
  12. C VAR0( 1 ) = IFIS :Indicateur de fissuration (0 1 2)
  13. C VAR0( 2 ) = ANGL :Angle de fissuration
  14. C VAR0( 3 ) = IPLA :Indicateur d'etat en bicompres. (0 1 2 3)
  15. C VAR0( 4 ) = SIG1 :Force d écrouissage de traction
  16. C VAR0( 5 ) = SIG2 :Force d écrouissage de compression
  17. C VAR0( 6 ) = EPS1 :Variable d écrouissage de traction
  18. C VAR0( 7 ) = EPS2 :Variable d écrouissage de compression
  19. C VAR0( 8 ) = TDEF :Taux de déformation
  20. C VAR0( 9 ) = TCON :Taux de contrainte
  21. C VAR0( 10 ) = SIGP(1):Contrainte plastique en mode
  22. C viscoplastique
  23. C VAR0( 11 ) = SIGP(2):
  24. C VAR0( 12 ) = SIGP(3):
  25. C VAR0( 13 ) = SIGP(4):
  26. C VAR0( 14 ) = DPSTV1 :Variable d'écrouissage plastique
  27. C en mode viscoplastique
  28. C VAR0( 15 ) = DPSTV2 :
  29. C VAR0( 16 ) = SIGV1 :Force d ecrouissage de traction
  30. C en mode viscoplastique
  31. C VAR0( 17 ) = SIGV2 :Force d ecrouissage de compression
  32. C en mode viscoplastique
  33. C
  34. C XMAT(NCOMAT) = COMPOSANTES DE MATERIAU
  35. C NCOMAT = NOMBRE DE COMPOSANTES DE MATERIAU
  36. C SORTIES
  37. C SIGF(NSTRS) = CONTRAINTES FINALES
  38. C VARF(NVARI) = VARIABLES INTERNES FINALES
  39. C KERRE = 0 TOUT OK
  40. C---------------------------------------------------------------------
  41. C
  42. C IFOUR INDICE DU TYPE DE PROBLEME
  43. C -2 CONTRAINTES PLANES
  44. C -1 DEFORMATIONS PLANES
  45. C 0 AXISYMETRIQUE
  46. C 1 SERIE DE FOURIER
  47. C 2 TRIDIMENSIONNEL
  48. C---------------------------------------------------------------------
  49. C COMPOSANTES DE MATERIAU
  50. C=====================================================================
  51. C YOUN : Module d'Young
  52. C XNU : Coeficient de Poisson
  53. C RHO : Masse volumique (Facultatif)
  54. C ALPH : Coeficient de dilation thermique (Facultatif)
  55. C ALFA : Resis. tract. simple / resis. compr. simple
  56. C BETA : Resis. compr. biax. / resis. compr. simple (Modele de NADAI)
  57. C RB : Resis. compr. simple
  58. C GFC : Energie de rupt en compression
  59. C GFT : Energie de rupt en traction
  60. C ETA : parametre de viscosite
  61. C ICT : Choix de la courbe de traction (Calibrage)
  62. C ICC : Choix de la courbe de compression (Calibrage)
  63. C
  64. C IVIS : Modele visqueux ou non
  65. C 0 : non visqueux
  66. C 1 : viscoplastique
  67. C 2 : viscoelastoplastique
  68. C 3 : elastoplastique endommageable
  69. C 4 : elasto viscoplastique endommageable
  70. C
  71. C IMOD : Choix du modele Beton
  72. C 1 = Modele BETON_INSA ( Von MISES + Sigma Max en trac.)
  73. C Rheologie 2D
  74. C 2 = Modele BETON_INSA ( Von MISES + Sigma Max en trac.)
  75. C Rheologie 3D
  76. C 3 = Modele BETON_INSA ( Druck Prager + Sigma Max en trac.)
  77. C Rheologie 2D
  78. C 4 = Modele BETON_INSA ( Druck Prager + Sigma Max en trac.)
  79. C Rheologie 3D
  80. C ITER : nombre d'iterations internes admissibles
  81. C Pour modeles fluage
  82. C TZER = DATE MINI POUR APPLICATION DU PREMIER CHARGEMENT
  83. C PGDUR = CHOIX DE PROGRESSION DES PAS DE TEMPS DE DISCRETISATION DES COURBES DE RELAXATION
  84. C PGTZO = CHOIX DE PROGRESSION DES PAS DE TEMPS DE DISCRETISATION DES COURBES DE RELAXATION
  85. C ITYPE = CHOIX DU CODE DE CALCUL DE FLUAGE CEB-FIP, B3, EUROCODE2, ETC
  86. C TAU2 = SECOND TEMPS DE RELAXATION PERMETTANT LE CALCUL DE LA SUITE DES TEMPS DE RELAXATION
  87. C TAU1 = PREMIER TEMPS DE RELAXATION PERMETTANT LA REPONSE ELASTIQUE
  88. C DATSEC = DATE DE DEBUT DE SECHAGE
  89. C DATCUR = DATE DE FIN DE CURE DU BETON
  90. C NBRC = NOMBRE DE BRANCHE DE MAXWELL
  91. C DATCOU = DATE D APPLICATION DU CHARGEMENT
  92. C IFOR = 1 formulation solide
  93. C = 2 liquide
  94. C---------------------------------------------------------------------
  95. C
  96. IMPLICIT INTEGER(I-N)
  97. IMPLICIT REAL*8(A-H,O-Z)
  98. -INC DECHE
  99. CHARACTER*40 FMT,TITRE
  100. DIMENSION SIR(9,4),EPST(4),SIGP(4)
  101. DIMENSION D(6,6),D1(6,6),STRN(4)
  102. DIMENSION SIGR(4),VART(200),VV(36),SIGMF(4)
  103. DIMENSION SIRL(8,4)
  104. C
  105. SEGMENT WRK4
  106. REAL*8 XE(3,NBNN)
  107. ENDSEGMENT
  108.  
  109. SEGMENT BETJEF
  110. REAL*8 AA,BETA,RB,ALFA,YOUN,XNU,GFC,GFT,CAR,ETA,TDEF,
  111. & TCON,DPSTF1,DPSTF2,TETA,PDT,TP00
  112. INTEGER ICT,ICC,IMOD,IVISS,ITER,
  113. & ISIM,IBB,IGAU1,IZON
  114. ENDSEGMENT
  115. SEGMENT VISCO
  116. REAL*8 DPSTV1,DPSTV2,SIGV1,SIGV2,ENDV
  117. ENDSEGMENT
  118. SEGMENT BETFLU
  119. REAL*8 DATCOU,DATCUR,DATSEC,E28,PGTZO,PGDUR,TAU1,TAU2,
  120. & TP0,TZER
  121. INTEGER ITYPE,IMD,NBRC,NCOE,NTZERO,NTPS,IFOR
  122. ENDSEGMENT
  123. SEGMENT NECH0
  124. REAL*8 DT,DC,ALFG,S0,ENDO
  125. ENDSEGMENT
  126. SEGMENT NECH1
  127. REAL*8 DLMT
  128. ENDSEGMENT
  129. SEGMENT NECH2
  130. REAL*8 ATR,GTR,ALPH0
  131. ENDSEGMENT
  132. SEGMENT NECH3
  133. REAL*8 RBT,ALFAT,YOUNT,GFCT,GFTT,ALPH
  134. ENDSEGMENT
  135. C
  136. C NCOMAT = nmatt
  137. C NBNN=XE(/2)
  138. NVARI=VAR0(/1)
  139. SEGINI BETJEF
  140. SEGINI VISCO
  141. SEGINI BETFLU
  142. SEGINI NECH0
  143. SEGINI NECH1
  144. SEGINI NECH2
  145. SEGINI NECH3
  146. PDT = DDT
  147. TP0=TEMP0
  148. TP00=TEMP0
  149. CAR = HCAR
  150. CALL ZERO(SIGP,4,1)
  151. CALL ZERO(SIR,9,4)
  152. CALL ZERO(SIRL,8,4)
  153. C
  154. C TEST DE CONSISTANCE DES DONNEES
  155. C
  156. YOUNT = 0.D0
  157. YOUN = 0.D0
  158. XNU = 0.D0
  159. RHO = 0.D0
  160. ALPH = 0.D0
  161. ALFAT = 0.D0
  162. ALFA = 0.D0
  163. BETA = 0.D0
  164. RBT = 0.D0
  165. RB = 0.D0
  166. GFCT = 0.D0
  167. GFTT = 0.D0
  168. GFC = 0.D0
  169. GFT = 0.D0
  170. ETA = 0.D0
  171. ITER = 0
  172. IMOD = 0
  173. ICT = 0
  174. ICC = 0
  175. DC =0.D0
  176. DT =0.D0
  177. IASO =0
  178. ICOU =0
  179. ENDO = 0.D0
  180. ENDV =0.D0
  181. C
  182. AA=0.D0
  183. BB=0.D0
  184. DK1=0.D0
  185. DK2=0.D0
  186. CNC1=0.D0
  187. CNC1=0.D0
  188. CNT1=0.D0
  189. CNT2=0.D0
  190. C
  191. YOUNT = XMAT( 1)
  192. XNU = XMAT( 2)
  193. RHO = XMAT( 3)
  194. ALPH = XMAT( 4)
  195. ALFAT = XMAT( 5)
  196. BETA = XMAT( 6)
  197. RBT = XMAT( 7)
  198. GFCT = XMAT( 8)
  199. GFTT = XMAT( 9)
  200. ITER = INT(REAL(XMAT(10)))
  201. IMOD = INT(REAL(XMAT(11)))
  202. IF (IVIS.EQ.0) THEN
  203. ICT = INT(REAL(XMAT(12)))
  204. ICC = INT(REAL(XMAT(13)))
  205. ENDIF
  206. IF (IVIS.EQ.1) THEN
  207. ICT = INT(REAL(XMAT(12)))
  208. ICC = INT(REAL(XMAT(13)))
  209. ETA = XMAT(14)
  210. ISIM = 1
  211. ENDIF
  212. IF (IVIS.EQ.2) THEN
  213. ICT = INT(REAL(XMAT(12)))
  214. ICC = INT(REAL(XMAT(13)))
  215. DATCOU = XMAT(14)
  216. NBRC = INT(REAL(XMAT(15)))
  217. IFOR = INT(REAL(XMAT(16)))
  218. ITYPE = INT(REAL(XMAT(17)))
  219. DATCUR = XMAT(18)
  220. DATSEC = XMAT(19)
  221. E28= YOUN
  222. ENDIF
  223. IF (IVIS.EQ.3) THEN
  224. DC = XMAT(12)
  225. DT = XMAT(13)
  226. IASO = INT(REAL(XMAT(14)))
  227. ICOU = INT(REAL(XMAT(15)))
  228. ENDIF
  229. IF (IVIS.EQ.4) THEN
  230. DC = XMAT(12)
  231. DT = XMAT(13)
  232. IASO = INT(REAL(XMAT(14)))
  233. ICOU = INT(REAL(XMAT(15)))
  234. ETA = XMAT(16)
  235. ENDIF
  236. C
  237. C---------------------------------------------------------------------
  238. C INITIALISATION DES VARIABLES DU MODELE D'ENDOMMAGEMENT
  239. C---------------------------------------------------------------------
  240. C A: CALCUL ASSOCIE OU NON
  241. C
  242. IF(IASO.EQ.0) THEN
  243. ALFG = (BETA-1.D0)/(2*BETA-1.D0)
  244. ELSE
  245. ALFG = 0.2
  246. ENDIF
  247. C---------------------------------------------------------------------
  248. C B: CALCUL COUPLE OU NON
  249. C
  250. IF(ICOU.EQ.0) THEN
  251. ATR=0.D0
  252. GTR=0.D0
  253. ELSE
  254. ATR=2.075D0
  255. GTR=0.235D0
  256. ENDIF
  257. C---------------------------------------------------------------------
  258. C
  259. DLMT=0.9999D0
  260. C S0 = 0.25D0
  261. IF(IVIS.EQ.4)THEN
  262. C attention en comportement visqueux refermeture pose probleme
  263. S0 = 1.D0
  264. ELSE
  265. S0 = 0.25D0
  266. C S0 = 1.D0
  267. ENDIF
  268. C
  269. C---------------------------------------------------------------------
  270. C CAS DU COUPLAGE THERMIQUE-MECANIQUE
  271. C---------------------------------------------------------------------
  272. C IF IASO=0 Utilisation des lois f(T) locales
  273. C IF IASO=1 Utilisation des lois du fichier de donnée
  274. C
  275. ITHHER = 1
  276. IF(ITHHER.NE.0)THEN
  277. CALL EVOLTH(T0,TF,BETJEF,NECH0,NECH3)
  278. ELSE
  279. CC write(*,*)'ITHHER=',ITHHER
  280. RB = RBT
  281. ALFA = ALFAT
  282. YOUN = YOUNT
  283. GFC = GFCT
  284. GFT = GFTT
  285. ENDIF
  286. C
  287. C
  288. C---------------------------------------------------------------------
  289. C
  290. C CORRESPONDANCE DES VARIABLES POUR URGCST
  291. C
  292. CALL ZERO(D,6,6)
  293. CALL ZERO(D1,6,6)
  294. CALL ZERO(VART,NVARI,1)
  295. CALL ZERO(SIGR,4,1)
  296. C
  297. IVISS = IVIS
  298. IBB = IB
  299. IGAU1 = IGAU
  300. EPAIST=1.D0
  301. NSTRS=NSTRSS
  302. IFOU=IFOUR
  303. *
  304. DO 1 I=1,NSTRS
  305. SIGR(I)=SIG0(I)/EPAIST
  306. * WRITE(*,*)'SIG0',I,'=',SIG0(I)
  307. STRN(I)=DEPST(I)
  308. * WRITE(*,*)'DEPST',I,'=',DEPST(I)
  309. 1 CONTINUE
  310. *
  311. IF((MELE.EQ.28.OR.MFR.EQ.3).AND.NSTRSS.EQ.4) THEN
  312. NSTRS=3
  313. IFOU=-2
  314. SIGR(3)=SIG0(4)
  315. STRN(3)=DEPST(4)
  316. STRN(4)=DEPST(3)
  317. ENDIF
  318. *
  319. IF((NSTRS.EQ.4).AND.IFOUR.EQ.-2.AND.
  320. * (IMOD.NE.2.AND.IMOD.NE.4)) THEN
  321. NSTRS=3
  322. SIGR(3)=SIG0(4)
  323. STRN(3)=DEPST(4)
  324. STRN(4)=DEPST(3)
  325. ENDIF
  326. *
  327. DO 22 II=1,NVARI
  328. VART(II)=VAR0(II)
  329. 22 CONTINUE
  330. C------------------------------------------------------
  331. C INITIALISATION DES VARIABLES INTERNES
  332. C------------------------------------------------------
  333. C A: MODELE PLASTIQUE
  334. C
  335. IFISU = INT(REAL(VART(1)))
  336. ANGL = VART(2)
  337. IPLA = INT(REAL(VART(3)))
  338. SIG1= VART(4)
  339. SIG2= VART(5)
  340. DPSTF1 = VART(6)
  341. DPSTF2 = VART(7)
  342. C
  343. C------------------------------------------------------
  344. C B: MODELE VISCOPLASTIQUE CLASSIQUE
  345. C
  346. IF (IVIS.EQ.1) THEN
  347. TDEF = VART( 8)
  348. TCON = VART( 9)
  349. SIGP(1) = VART(10)
  350. SIGP(2) = VART(11)
  351. SIGP(3) = VART(12)
  352. SIGP(4) = VART(13)
  353. DPSTV1 = VART(14)
  354. DPSTV2 = VART(15)
  355. SIGV1 = VART(16)
  356. SIGV2 = VART(17)
  357. ENDIF
  358. C
  359. C------------------------------------------------------
  360. C C: MODELE VISCO-ELASTO PLASTIQUE
  361. C
  362. IF (IVIS.EQ.2.AND.IFOR.EQ.1) THEN
  363. C modele visco elasto plastique solide de MAXWELL
  364. C-------------------------------------------------
  365. MC = NBRC + 1
  366. DO 101 I1 = 1,MC
  367. DO 102 J1 = 1,NSTRS
  368. NV = NSTRS * (I1 - 1) + J1
  369. SIR(I1,J1) = VART(7 + NV)
  370. 102 CONTINUE
  371. 101 CONTINUE
  372. C
  373. ELSE IF (IVIS.EQ.2.AND.IFOR.EQ.2) THEN
  374. C modele visco elasto plastique liquide de MAXWELL
  375. C-------------------------------------------------
  376. MC = NBRC + 1
  377. DO 131 I2 = 1,MC
  378. DO 132 J2 = 1,NSTRS
  379. NV = NSTRS * (I2 - 1) + J2
  380. SIRL(I2,J2) = VART(7 + NV)
  381. 132 CONTINUE
  382. 131 CONTINUE
  383. C
  384. ENDIF
  385. C
  386. C----------------------------------------------
  387. C D: MODELE PLASTIQUE ENDOMMAGEABLE
  388. C
  389. IF (IVIS.EQ.3) THEN
  390. ENDO = VART( 8)
  391. DECO = VART( 9)
  392. ENDIF
  393. C
  394. C----------------------------------------------
  395. C D: MODELE VISCOPLASTIQUE ENDOMMAGEABLE
  396. C
  397. IF (IVIS.EQ.4) THEN
  398. TDEF = VART( 8)
  399. TCON = VART( 9)
  400. SIGP(1) = VART(10)
  401. SIGP(2) = VART(11)
  402. SIGP(3) = VART(12)
  403. SIGP(4) = VART(13)
  404. DPSTV1 = VART(14)
  405. DPSTV2 = VART(15)
  406. SIGV1 = VART(16)
  407. SIGV2 = VART(17)
  408. ENDO = VART(18)
  409. ENDV = VART(19)
  410. ENDIF
  411. C -------------------------------------
  412. TETA=ANGL
  413. C---------------------------------------------------------------------
  414. C
  415. GO TO (10,10,30,40),NSTRS
  416. 10 CONTINUE
  417. KERRE=437
  418. WRITE(*,*) '!! ATTENTION DANS BETON NSTRS=',NSTRS
  419. STOP
  420. GO TO 1000
  421. C
  422. 40 CONTINUE
  423. C
  424. 30 CONTINUE
  425. C
  426. C---------------------------------------------------------------------
  427. C
  428. IF(IMOD.GE.1.AND.IMOD.LE.4) THEN
  429. * *-----------------------*
  430. * * MODELE URGC ST *
  431. * *-----------------------*
  432. *
  433. CALL BONE(SIGR,SIGMF,STRN,IPLA,IFISU,SIG1,SIG2
  434. A ,NSTRS,D,D1,IFOU,SIGP,EPST,SIR,SIRL,
  435. B ENDO,ITHHER,T0,TF,BETJEF,VISCO,BETFLU,
  436. C NECH0,NECH1,NECH2,NECH3)
  437. C
  438. ELSE
  439. WRITE(*,*) '!! ATTENTION CE MODELE N EXISTE PAS (URGCST)'
  440. STOP
  441. RETURN
  442. ENDIF
  443. C---------------------------------------------------------------------
  444. C
  445. C
  446. IF((NSTRS.EQ.3).AND.IFOU.EQ.-2.AND.
  447. * (IMOD.NE.2.AND.IMOD.NE.4)) THEN
  448. SIGMF(4)=SIGMF(3)
  449. SIGMF(3)=0.D0
  450. ENDIF
  451. C
  452. IF((MELE.EQ.28.OR.MFR.EQ.3).AND.NSTRSS.EQ.3) THEN
  453. SIGMF(3)=0.D0
  454. SIGMF(4)=SIGMF(3)
  455. ENDIF
  456. C
  457. DO 2 I=1,NSTRSS
  458. SIGF(I)=SIGMF(I)
  459. C WRITE(*,*)'SIGMF(',I,')=',SIGMF(I)
  460. 2 CONTINUE
  461. C
  462. C---------------------------------------------------------------------
  463. C
  464. C---------------------------------------------------------------------
  465. C ACTUALISATION DES VARIABLES INTERNES
  466. C---------------------------------------------------------------------
  467. C A: MODELE PLASTIQUE
  468. C
  469. VART( 1)=FLOAT(IFISU)
  470. VART( 2)=TETA
  471. VART( 3)=FLOAT(IPLA)
  472. VART( 4)=SIG1
  473. VART( 5)=SIG2
  474. VART( 6)=DPSTF1
  475. VART( 7)=DPSTF2
  476. K1 = 7
  477. C
  478. C---------------------------------------------------------------------
  479. C B: MODELE VISCOPLASTIQUE CLASSIQUE
  480. C
  481. IF (IVIS.EQ.1) THEN
  482. VART(8)=TDEF
  483. VART(9)=TCON
  484. VART(10)=SIGP(1)
  485. VART(11)=SIGP(2)
  486. VART(12)=SIGP(3)
  487. VART(13)=SIGP(4)
  488. VART(14)=DPSTV1
  489. VART(15)=DPSTV2
  490. VART(16)=SIGV1
  491. VART(17)=SIGV2
  492. K1 = 17
  493. ENDIF
  494. C---------------------------------------------------------------------
  495. C C: MODELE VISCO-ELASTO PLASTIQUE
  496. C
  497. IF (IVIS.EQ.2.AND.IFOR.EQ.1) THEN
  498. MC = NBRC + 1
  499. DO 103 I1 = 1,MC
  500. DO 104 J1 = 1,NSTRS
  501. NV = NSTRS * (I1 - 1) + J1
  502. VART(7 + NV) = SIR(I1,J1)
  503. 104 CONTINUE
  504. 103 CONTINUE
  505. C DO 123 I1 = 1,MC
  506. C DO 124 J1 = 1,MC
  507. C K = MC * (I1 - 1) + J1
  508. C VART(43 + K) = CODU(I1,J1)
  509. C 124 CONTINUE
  510. C 123 CONTINUE
  511. K1=43
  512. ENDIF
  513. C----------------------------------------
  514. IF (IVIS.EQ.2.AND.IFOR.EQ.2) THEN
  515. DO 108 I1 = 1,NBRC
  516. DO 109 J1 = 1,NSTRS
  517. NV = NSTRS * (I1 - 1) + J1
  518. VART(7 + NV) = SIRL(I1,J1)
  519. 109 CONTINUE
  520. 108 CONTINUE
  521. C DO 143 I1 = 1,NBRC
  522. C DO 144 J1 = 1,NBRC
  523. C K = NBRC * (I1-1) + J1
  524. C VART(43 + K) = CODL(I1,J1)
  525. C 144 CONTINUE
  526. C 143 CONTINUE
  527. K1=43
  528. ENDIF
  529. C
  530. C---------------------------------------------------------------------
  531. C D: MODELE PLASTIQUE ENDOMMAGEABLE
  532. C
  533. IF (IVIS.EQ.3) THEN
  534. VART(8) = ENDO
  535. VART(9) = MAX(VAR0(9),ENDO)
  536. K1= 9
  537. ENDIF
  538. C
  539. C---------------------------------------------------------------------
  540. C B: MODELE VISCOPLASTIQUE ENDOMMAGEABLE
  541. C
  542. IF (IVIS.EQ.4) THEN
  543. VART(8)=TDEF
  544. VART(9)=TCON
  545. VART(10)=SIGP(1)
  546. VART(11)=SIGP(2)
  547. VART(12)=SIGP(3)
  548. VART(13)=SIGP(4)
  549. VART(14)=DPSTV1
  550. VART(15)=DPSTV2
  551. VART(16)=SIGV1
  552. VART(17)=SIGV2
  553. VART(18)=ENDO
  554. VART(19)=ENDV
  555. K1 = 19
  556. ENDIF
  557. C
  558. GO TO (11,11,31,41),NSTRS
  559. 11 CONTINUE
  560. KERRE=437
  561. WRITE(*,*) '!! ATTENTION DANS BETON NSTRS=',NSTRS
  562. STOP
  563. GO TO 1000
  564. C
  565. 41 CONTINUE
  566. C
  567. 31 CONTINUE
  568. C
  569. C CONTRAINTES OU DEFORMATIONS PLANES
  570. C SANS CISAILLEMENTS TRANSVERSAL
  571. C
  572. L1=0
  573. DO 66 I1 = 1,NSTRS
  574. DO 67 J1 = 1,NSTRS
  575. L1=L1+1
  576. VV(L1) = D(J1,I1)
  577. 67 CONTINUE
  578. 66 CONTINUE
  579. *
  580. L1=0
  581. DO 62 I1 = 1,NSTRS
  582. DO 63 J1 = 1,NSTRS
  583. K=K1+(I1-1)*NSTRS+J1
  584. L1=L1+1
  585. VART(K)=VV(L1)
  586. 63 CONTINUE
  587. 62 CONTINUE
  588. C
  589. DO 23 II=1,NVARI
  590. VARF(II)=VART(II)
  591. 23 CONTINUE
  592. C
  593. C---------------------------------------------------------------------
  594. C
  595. 1000 CONTINUE
  596. SEGSUP BETJEF
  597. SEGSUP VISCO
  598. SEGSUP BETFLU
  599. SEGSUP NECH0
  600. SEGSUP NECH1
  601. SEGSUP NECH2
  602. SEGSUP NECH3
  603. C
  604. RETURN
  605. END
  606.  
  607.  
  608.  
  609.  
  610.  
  611.  
  612.  
  613.  
  614.  
  615.  
  616.  
  617.  
  618.  
  619.  
  620.  

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