Télécharger idvisc.eso

Retour à la liste

Numérotation des lignes :

  1. C IDVISC SOURCE YURI 11/10/28 22:15:28 7181
  2.  
  3. C=======================================================================
  4. C= Noms de composantes de materiaux viscoplastiques =
  5. C= =
  6. C= ENTREES : =
  7. C= --------- =
  8. C= MFR Numero de formulation =
  9. C= MOOBL Pointeur sur les mots OBLIGATOIRES =
  10. C= MOFAC Pointeur sur les mots FACULTATIFS =
  11. C= IPLAC Numero du materiau dans MODVIS =
  12. C= =
  13. C= SORTIE : =
  14. C= -------- =
  15. C= IRET =0 si probleme =
  16. C=======================================================================
  17.  
  18. SUBROUTINE IDVISC (MFR,MOOBL,MOFAC,IPLAC,IRET)
  19.  
  20. IMPLICIT INTEGER(I-N)
  21. IMPLICIT REAL*8 (A-H,O-Z)
  22.  
  23. -INC CCOPTIO
  24. -INC SMLMOTS
  25. JGN=8
  26. IRET=1
  27. C ===
  28. C 1 - Modele viscoplastique de GUIONNET
  29. C ===
  30. IF (IPLAC.EQ.1) THEN
  31. MLMOTS=MOOBL
  32. JGA=MOTS(/2)
  33. JGM=JGA+20
  34. SEGADJ,MLMOTS
  35. MOTS(JGA+ 1)='N '
  36. MOTS(JGA+ 2)='K '
  37. MOTS(JGA+ 3)='CK '
  38. MOTS(JGA+ 4)='R0 '
  39. MOTS(JGA+ 5)='CD '
  40. MOTS(JGA+ 6)='M '
  41. MOTS(JGA+ 7)='A '
  42. MOTS(JGA+ 8)='A1 '
  43. MOTS(JGA+ 9)='C '
  44. MOTS(JGA+10)='C1 '
  45. MOTS(JGA+11)='G '
  46. MOTS(JGA+12)='R '
  47. MOTS(JGA+13)='NN '
  48. MOTS(JGA+14)='C2 '
  49. MOTS(JGA+15)='Q '
  50. MOTS(JGA+16)='G1 '
  51. MOTS(JGA+17)='R1 '
  52. MOTS(JGA+18)='BETA'
  53. MOTS(JGA+19)='C0 '
  54. MOTS(JGA+20)='P1M0'
  55. MOOBL=MLMOTS
  56. RETURN
  57. C ===
  58. C 2 - Modele viscoplastique de CHABOCHE
  59. C ===
  60. ELSE IF (IPLAC.EQ.2) THEN
  61. MLMOTS=MOOBL
  62. JGA=MOTS(/2)
  63. JGM=JGA+22
  64. SEGADJ,MLMOTS
  65. MOTS(JGA+ 1)='ALF '
  66. MOTS(JGA+ 2)='N '
  67. MOTS(JGA+ 3)='KK '
  68. MOTS(JGA+ 4)='K0 '
  69. MOTS(JGA+ 5)='ALFK'
  70. MOTS(JGA+ 6)='ALFR'
  71. MOTS(JGA+ 7)='C1 '
  72. MOTS(JGA+ 8)='A1 '
  73. MOTS(JGA+ 9)='PHI '
  74. MOTS(JGA+10)='C2 '
  75. MOTS(JGA+11)='A2 '
  76. MOTS(JGA+12)='B '
  77. MOTS(JGA+13)='BET1'
  78. MOTS(JGA+14)='BET2'
  79. MOTS(JGA+15)='R1 '
  80. MOTS(JGA+16)='R2 '
  81. MOTS(JGA+17)='GAMA'
  82. MOTS(JGA+18)='M '
  83. MOTS(JGA+19)='ETA '
  84. MOTS(JGA+20)='MU '
  85. MOTS(JGA+21)='QMAX'
  86. MOTS(JGA+22)='QSTA'
  87. MOOBL=MLMOTS
  88. RETURN
  89. C ===
  90. C 3 - Modele d'ecrouissage-VISCOsite multiplicative avec enDOMMAGEment
  91. C ===
  92. ELSE IF (IPLAC.EQ.3) THEN
  93. MLMOTS=MOOBL
  94. JGA=MOTS(/2)
  95. JGM=JGA+8
  96. SEGADJ,MLMOTS
  97. MOTS(JGA+1)='N '
  98. MOTS(JGA+2)='M '
  99. MOTS(JGA+3)='KK '
  100. MOTS(JGA+4)='ALP1'
  101. MOTS(JGA+5)='BLP1'
  102. MOTS(JGA+6)='R '
  103. MOTS(JGA+7)='A '
  104. MOTS(JGA+8)='EVOL'
  105. MOOBL=MLMOTS
  106. MLMOTS=MOFAC
  107. JGA=MOTS(/2)
  108. JGM=JGA+1
  109. SEGADJ,MLMOTS
  110. MOTS(JGA+1)='SMAX'
  111. MOFAC=MLMOTS
  112. RETURN
  113. C ===
  114. C 4 - Modele viscoplastique PARFAIT
  115. C ===
  116. ELSE IF (IPLAC.EQ.4) THEN
  117. MLMOTS=MOOBL
  118. JGA=MOTS(/2)
  119. JGM=JGA+3
  120. SEGADJ,MLMOTS
  121. MOTS(JGA+1)='SIGY'
  122. MOTS(JGA+2)='N '
  123. MOTS(JGA+3)='K '
  124. MOOBL=MLMOTS
  125. RETURN
  126. C ===
  127. C 5 - Modele viscoplastique POUDRE_A
  128. C ===
  129. ELSE IF (IPLAC.EQ.5) THEN
  130. MLMOTS=MOOBL
  131. JGA=MOTS(/2)
  132. JGM=JGA+16
  133. SEGADJ,MLMOTS
  134. MOTS(JGA+ 1)='N '
  135. MOTS(JGA+ 2)='A '
  136. MOTS(JGA+ 3)='QSRT'
  137. MOTS(JGA+ 4)='F0 '
  138. MOTS(JGA+ 5)='F1 '
  139. MOTS(JGA+ 6)='F2 '
  140. MOTS(JGA+ 7)='F3 '
  141. MOTS(JGA+ 8)='F4 '
  142. MOTS(JGA+ 9)='F5 '
  143. MOTS(JGA+10)='C0 '
  144. MOTS(JGA+11)='C1 '
  145. MOTS(JGA+12)='C2 '
  146. MOTS(JGA+13)='C3 '
  147. MOTS(JGA+14)='C4 '
  148. MOTS(JGA+15)='C5 '
  149. MOTS(JGA+16)='RHOR'
  150. MOOBL=MLMOTS
  151. RETURN
  152. C ===
  153. C 6 - Modele viscoplastique POUDRE_B
  154. C ===
  155. ELSE IF (IPLAC.EQ.6) THEN
  156. MLMOTS=MOOBL
  157. JGA=MOTS(/2)
  158. JGM=JGA+23
  159. SEGADJ,MLMOTS
  160. MOTS(JGA+ 1)='N '
  161. MOTS(JGA+ 2)='A '
  162. MOTS(JGA+ 3)='QSRT'
  163. MOTS(JGA+ 4)='F0 '
  164. MOTS(JGA+ 5)='F1 '
  165. MOTS(JGA+ 6)='F2 '
  166. MOTS(JGA+ 7)='F3 '
  167. MOTS(JGA+ 8)='F4 '
  168. MOTS(JGA+ 9)='F5 '
  169. MOTS(JGA+10)='C0 '
  170. MOTS(JGA+11)='C1 '
  171. MOTS(JGA+12)='C2 '
  172. MOTS(JGA+13)='C3 '
  173. MOTS(JGA+14)='C4 '
  174. MOTS(JGA+15)='C5 '
  175. MOTS(JGA+16)='G0 '
  176. MOTS(JGA+17)='G1 '
  177. MOTS(JGA+18)='G2 '
  178. MOTS(JGA+19)='G3 '
  179. MOTS(JGA+20)='G4 '
  180. MOTS(JGA+21)='G5 '
  181. MOTS(JGA+22)='G6 '
  182. MOTS(JGA+23)='RHOR'
  183. MOOBL=MLMOTS
  184. RETURN
  185. C ===
  186. C 7 - Modele viscoplastique de OHNO
  187. C ===
  188. ELSE IF (IPLAC.EQ.7) THEN
  189. MLMOTS=MOOBL
  190. JGA=MOTS(/2)
  191. JGM=JGA+24
  192. SEGADJ,MLMOTS
  193. MOTS(JGA+ 1)='ALF '
  194. MOTS(JGA+ 2)='N '
  195. MOTS(JGA+ 3)='KK '
  196. MOTS(JGA+ 4)='K0 '
  197. MOTS(JGA+ 5)='ALFK'
  198. MOTS(JGA+ 6)='ALFR'
  199. MOTS(JGA+ 7)='C1 '
  200. MOTS(JGA+ 8)='A1 '
  201. MOTS(JGA+ 9)='PHI '
  202. MOTS(JGA+10)='C2 '
  203. MOTS(JGA+11)='A2 '
  204. MOTS(JGA+12)='B '
  205. MOTS(JGA+13)='BET1'
  206. MOTS(JGA+14)='BET2'
  207. MOTS(JGA+15)='R1 '
  208. MOTS(JGA+16)='R2 '
  209. MOTS(JGA+17)='GAMA'
  210. MOTS(JGA+18)='M '
  211. MOTS(JGA+19)='ETA '
  212. MOTS(JGA+20)='MU '
  213. MOTS(JGA+21)='QMAX'
  214. MOTS(JGA+22)='QSTA'
  215. MOTS(JGA+23)='EXP1'
  216. MOTS(JGA+24)='EXP2'
  217. MOOBL=MLMOTS
  218. RETURN
  219. C ===
  220. C 8 - Modele viscoplastique DDI
  221. C ===
  222. ELSE IF (IPLAC.EQ.8) THEN
  223. MLMOTS=MOOBL
  224. JGA=MOTS(/2)
  225. JGM=JGA+18
  226. SEGADJ,MLMOTS
  227. MOTS(JGA+ 1)='RP0 '
  228. MOTS(JGA+ 2)='QP '
  229. MOTS(JGA+ 3)='BP '
  230. MOTS(JGA+ 4)='CP1 '
  231. MOTS(JGA+ 5)='CP2 '
  232. MOTS(JGA+ 6)='DP1 '
  233. MOTS(JGA+ 7)='DP2 '
  234. MOTS(JGA+ 8)='KS '
  235. MOTS(JGA+ 9)='N '
  236. MOTS(JGA+10)='RV0 '
  237. MOTS(JGA+11)='QV '
  238. MOTS(JGA+12)='BV '
  239. MOTS(JGA+13)='CV1 '
  240. MOTS(JGA+14)='CV2 '
  241. MOTS(JGA+15)='DV1 '
  242. MOTS(JGA+16)='DV2 '
  243. MOTS(JGA+17)='CVP1'
  244. MOTS(JGA+18)='CVP2'
  245. MOOBL=MLMOTS
  246. MLMOTS=MOFAC
  247. JGA=MOTS(/2)
  248. JGM=JGA+10
  249. SEGADJ,MLMOTS
  250. MOTS(JGA+ 1)='ALP7'
  251. MOTS(JGA+ 2)='BETA'
  252. MOTS(JGA+ 3)='M '
  253. MOTS(JGA+ 4)='SU '
  254. MOTS(JGA+ 5)='SL '
  255. MOTS(JGA+ 6)='B '
  256. MOTS(JGA+ 7)='C '
  257. MOTS(JGA+ 8)='A '
  258. MOTS(JGA+ 9)='R '
  259. MOTS(JGA+10)='KK '
  260. MOFAC=MLMOTS
  261. RETURN
  262. C ===
  263. C 9 - Modele viscoplastique de KOCK
  264. C ===
  265. ELSE IF (IPLAC.EQ.9) THEN
  266. MLMOTS=MOOBL
  267. JGA=MOTS(/2)
  268. JGM=JGA+10
  269. SEGADJ,MLMOTS
  270. MOTS(JGA+ 1)='A '
  271. MOTS(JGA+ 2)='B '
  272. MOTS(JGA+ 3)='M '
  273. MOTS(JGA+ 4)='Q '
  274. MOTS(JGA+ 5)='R '
  275. MOTS(JGA+ 6)='H0 '
  276. MOTS(JGA+ 7)='AP '
  277. MOTS(JGA+ 8)='SB '
  278. MOTS(JGA+ 9)='N '
  279. MOTS(JGA+10)='S0 '
  280. MOOBL=MLMOTS
  281. RETURN
  282. C ====
  283. C 10 - Modele viscoplastique NOUAILLAS_A
  284. C ====
  285. ELSE IF (IPLAC.EQ.10) THEN
  286. MLMOTS=MOOBL
  287. JGA=MOTS(/2)
  288. JGM=JGA+22
  289. SEGADJ,MLMOTS
  290. MOTS(JGA+ 1)='ALF '
  291. MOTS(JGA+ 2)='N '
  292. MOTS(JGA+ 3)='KK '
  293. MOTS(JGA+ 4)='K0 '
  294. MOTS(JGA+ 5)='ALFK'
  295. MOTS(JGA+ 6)='ALFR'
  296. MOTS(JGA+ 7)='CL1 '
  297. MOTS(JGA+ 8)='DNL1'
  298. MOTS(JGA+ 9)='PHI '
  299. MOTS(JGA+10)='CL2 '
  300. MOTS(JGA+11)='DNL2'
  301. MOTS(JGA+12)='B '
  302. MOTS(JGA+13)='GDM1'
  303. MOTS(JGA+14)='GDM2'
  304. MOTS(JGA+15)='PTM1'
  305. MOTS(JGA+16)='PTM2'
  306. MOTS(JGA+17)='GAMA'
  307. MOTS(JGA+18)='M '
  308. MOTS(JGA+19)='ETA '
  309. MOTS(JGA+20)='MU '
  310. MOTS(JGA+21)='QMAX'
  311. MOTS(JGA+22)='QSTA'
  312. MOOBL=MLMOTS
  313. RETURN
  314. C ====
  315. C 11 - Modele viscoplastique NOUAILLAS_B
  316. C ====
  317. ELSE IF (IPLAC.EQ.11) THEN
  318. MLMOTS=MOOBL
  319. JGA=MOTS(/2)
  320. JGM=JGA+14
  321. SEGADJ,MLMOTS
  322. MOTS(JGA+ 1)='ALF '
  323. MOTS(JGA+ 2)='N '
  324. MOTS(JGA+ 3)='KK '
  325. MOTS(JGA+ 4)='K0 '
  326. MOTS(JGA+ 5)='CL1 '
  327. MOTS(JGA+ 6)='DNL1'
  328. MOTS(JGA+ 7)='PHI '
  329. MOTS(JGA+ 8)='CL2 '
  330. MOTS(JGA+ 9)='DNL2'
  331. MOTS(JGA+10)='B '
  332. MOTS(JGA+11)='GDM1'
  333. MOTS(JGA+12)='PTM1'
  334. MOTS(JGA+13)='RMAX'
  335. MOTS(JGA+14)='BR '
  336. MOOBL=MLMOTS
  337. RETURN
  338. C ====
  339. C 12 - Modele viscoplastique VISK2
  340. C ====
  341. ELSE IF (IPLAC.EQ.12) THEN
  342. MLMOTS=MOOBL
  343. JGA=MOTS(/2)
  344. JGM=JGA+5
  345. SEGADJ,MLMOTS
  346. MOTS(JGA+1)='SIGY'
  347. MOTS(JGA+2)='H '
  348. MOTS(JGA+3)='ETA '
  349. MOTS(JGA+4)='HVIS'
  350. MOTS(JGA+5)='N '
  351. MOOBL=MLMOTS
  352. RETURN
  353. C ====
  354. C 13 - Modele VISCOHINTE
  355. C ====
  356. ELSE IF (IPLAC.EQ.13) THEN
  357. MLMOTS=MOOBL
  358. JGA=MOTS(/2)
  359. JGM=JGA+7
  360. SEGADJ,MLMOTS
  361. MOTS(JGA+1)='GAM1'
  362. MOTS(JGA+2)='AL '
  363. MOTS(JGA+3)='KK '
  364. MOTS(JGA+4)='NN '
  365. MOTS(JGA+5)='YO '
  366. MOTS(JGA+6)='YC '
  367. MOTS(JGA+7)='MM '
  368. MOOBL=MLMOTS
  369. RETURN
  370. C ====
  371. C 14 - Modele viscoplastique MISTRAL
  372. C ====
  373. ELSE IF (IPLAC.EQ.14) THEN
  374. MLMOTS=MOOBL
  375. JGA=MOTS(/2)
  376. JGM=JGA+13
  377. SEGADJ,MLMOTS
  378. MOTS(JGA+ 1)='SIP1'
  379. MOTS(JGA+ 2)='SIP2'
  380. MOTS(JGA+ 3)='IBAS'
  381. MOTS(JGA+ 4)='DILT'
  382. MOTS(JGA+ 5)='NBRE'
  383. MOTS(JGA+ 6)='COHI'
  384. MOTS(JGA+ 7)='ACOU'
  385. MOTS(JGA+ 8)='EDIR'
  386. MOTS(JGA+ 9)='RVCE'
  387. MOTS(JGA+10)='ECRX'
  388. MOTS(JGA+11)='DVDI'
  389. MOTS(JGA+12)='CROI'
  390. MOTS(JGA+13)='INCR'
  391. MOOBL=MLMOTS
  392. RETURN
  393. C ====
  394. C 15 - Modele viscoplastique BETON_URGC
  395. C ====
  396. ELSE IF (IPLAC.EQ.15) THEN
  397. MLMOTS=MOOBL
  398. JGA=MOTS(/2)
  399. JGM=JGA+10
  400. SEGADJ,MLMOTS
  401. MOTS(JGA+ 1)='ALFA'
  402. MOTS(JGA+ 2)='BETA'
  403. MOTS(JGA+ 3)='LCS '
  404. MOTS(JGA+ 4)='GFC'
  405. MOTS(JGA+ 5)='GFT'
  406. MOTS(JGA+ 6)='ITON'
  407. MOTS(JGA+ 7)='IMOD'
  408. MOTS(JGA+ 8)='ICT '
  409. MOTS(JGA+ 9)='ICC'
  410. MOTS(JGA+10)='ETA'
  411. MOOBL=MLMOTS
  412. RETURN
  413. C ====
  414. C 16 - Modele LEMENDO (fluage de type LEMaitre avec ENDOmmagement)
  415. C ====
  416. ELSE IF (IPLAC.EQ.16) THEN
  417. MLMOTS=MOOBL
  418. JGA=MOTS(/2)
  419. JGM=JGA+9
  420. SEGADJ,MLMOTS
  421. MOTS(JGA+1)='SIGY'
  422. MOTS(JGA+2)='N '
  423. MOTS(JGA+3)='M '
  424. MOTS(JGA+4)='KK '
  425. MOTS(JGA+5)='ALP1'
  426. MOTS(JGA+6)='BLP1'
  427. MOTS(JGA+7)='R '
  428. MOTS(JGA+8)='A '
  429. MOTS(JGA+9)='EXK '
  430. MOOBL=MLMOTS
  431. RETURN
  432. C ====
  433. C 17 - Modele FLUNOR2 (fluage de type Norton)
  434. C ====
  435. ELSE IF (IPLAC.EQ.17) THEN
  436. MLMOTS=MOOBL
  437. JGA=MOTS(/2)
  438. JGM=JGA+6
  439. SEGADJ,MLMOTS
  440. MOTS(JGA+1)='SIGY'
  441. MOTS(JGA+2)='AF1 '
  442. MOTS(JGA+3)='AF2 '
  443. MOTS(JGA+4)='AF3 '
  444. MOTS(JGA+5)='SFFS'
  445. MOTS(JGA+6)='EFFS'
  446. MOOBL=MLMOTS
  447. RETURN
  448. C ====
  449. C 18 - Modele GATT_MONERIE
  450. C ====
  451. ELSE IF (IPLAC.EQ.18) THEN
  452. MLMOTS=MOOBL
  453. JGA=MOTS(/2)
  454. JGM=JGA+32
  455. SEGADJ,MLMOTS
  456. MOTS(JGA+ 1)='R '
  457. MOTS(JGA+ 2)='DG '
  458. MOTS(JGA+ 3)='K1 '
  459. MOTS(JGA+ 4)='M1 '
  460. MOTS(JGA+ 5)='Q1 '
  461. MOTS(JGA+ 6)='N1 '
  462. MOTS(JGA+ 7)='K2 '
  463. MOTS(JGA+ 8)='M2 '
  464. MOTS(JGA+ 9)='Q2 '
  465. MOTS(JGA+10)='N2 '
  466. MOTS(JGA+11)='OMEG'
  467. MOTS(JGA+12)='Q '
  468. MOTS(JGA+13)='H '
  469. MOTS(JGA+14)='BETA'
  470. MOTS(JGA+15)='K '
  471. MOTS(JGA+16)='A '
  472. MOTS(JGA+17)='Q3 '
  473. MOTS(JGA+18)='N3 '
  474. MOTS(JGA+19)='ADEN'
  475. MOTS(JGA+20)='KGON'
  476. MOTS(JGA+21)='POR0'
  477. MOTS(JGA+22)='BUMI'
  478. MOTS(JGA+23)='EFIS'
  479. MOTS(JGA+24)='DG0 '
  480. MOTS(JGA+25)='CR1 '
  481. MOTS(JGA+26)='CR2 '
  482. MOTS(JGA+27)='CR3 '
  483. MOTS(JGA+28)='CR '
  484. MOTS(JGA+29)='KP '
  485. MOTS(JGA+30)='AP '
  486. MOTS(JGA+31)='BP '
  487. MOTS(JGA+32)='QP '
  488. MOOBL=MLMOTS
  489. MLMOTS=MOFAC
  490. JGA=MOTS(/2)
  491. JGM=JGA+6
  492. SEGADJ,MLMOTS
  493. MOTS(JGA+ 1)='TYPE'
  494. MOTS(JGA+ 2)='COMP'
  495. MOTS(JGA+ 3)='DYN '
  496. MOTS(JGA+ 4)='DYN1'
  497. MOTS(JGA+ 5)='DYN2'
  498. MOTS(JGA+ 6)='DYN3'
  499. MOFAC=MLMOTS
  500. RETURN
  501. C ====
  502. C 19 - Modele UO2 = OTTOSEN + GATT_MONERIE
  503. C ====
  504. ELSE IF (IPLAC.EQ.19) THEN
  505. C= 19.1 - Modes de calcul PLAN DEFO, PLAN GENE, Axisymetrie, Fourier
  506. IF (IFOUR.EQ.-3.OR.IFOUR.EQ.-1.OR.
  507. . IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN
  508. MLMOTS=MOOBL
  509. JGA=MOTS(/2)
  510. JGM=JGA+39
  511. SEGADJ,MLMOTS
  512. MOTS(JGA+ 1)='R '
  513. MOTS(JGA+ 2)='DG '
  514. MOTS(JGA+ 3)='K1 '
  515. MOTS(JGA+ 4)='M1 '
  516. MOTS(JGA+ 5)='Q1 '
  517. MOTS(JGA+ 6)='N1 '
  518. MOTS(JGA+ 7)='K2 '
  519. MOTS(JGA+ 8)='M2 '
  520. MOTS(JGA+ 9)='Q2 '
  521. MOTS(JGA+10)='N2 '
  522. MOTS(JGA+11)='OMEG'
  523. MOTS(JGA+12)='Q '
  524. MOTS(JGA+13)='H '
  525. MOTS(JGA+14)='BETA'
  526. MOTS(JGA+15)='K '
  527. MOTS(JGA+16)='A '
  528. MOTS(JGA+17)='Q3 '
  529. MOTS(JGA+18)='N3 '
  530. MOTS(JGA+19)='ADEN'
  531. MOTS(JGA+20)='KGON'
  532. MOTS(JGA+21)='POR0'
  533. MOTS(JGA+22)='BUMI'
  534. MOTS(JGA+23)='EFIS'
  535. MOTS(JGA+24)='DG0 '
  536. MOTS(JGA+25)='CR1 '
  537. MOTS(JGA+26)='CR2 '
  538. MOTS(JGA+27)='CR3 '
  539. MOTS(JGA+28)='CR '
  540. MOTS(JGA+29)='KP '
  541. MOTS(JGA+30)='AP '
  542. MOTS(JGA+31)='BP '
  543. MOTS(JGA+32)='QP '
  544. IF (IFOUR.EQ.-3.OR.IFOUR.EQ.-1) THEN
  545. MOTS(JGA+33)='LXX '
  546. MOTS(JGA+34)='LYY '
  547. MOTS(JGA+35)='LXY '
  548. MOTS(JGA+36)='LZZ '
  549. MOTS(JGA+37)='PXX '
  550. MOTS(JGA+38)='PYY '
  551. MOTS(JGA+39)='PXY '
  552. C* ELSE IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN
  553. ELSE
  554. MOTS(JGA+33)='LRR '
  555. MOTS(JGA+34)='LZZ '
  556. MOTS(JGA+35)='LRZ '
  557. MOTS(JGA+36)='LOO '
  558. MOTS(JGA+37)='PRR '
  559. MOTS(JGA+38)='PZZ '
  560. MOTS(JGA+39)='PRZ '
  561. ENDIF
  562. MOOBL=MLMOTS
  563. MLMOTS=MOFAC
  564. JGA=MOTS(/2)
  565. JGM=JGA+28
  566. SEGADJ,MLMOTS
  567. MOTS(JGA+ 1)='GFTR'
  568. MOTS(JGA+ 2)='LTR '
  569. MOTS(JGA+ 3)='GS '
  570. MOTS(JGA+ 4)='BTR '
  571. MOTS(JGA+ 5)='WRUP'
  572. MOTS(JGA+ 6)='BILI'
  573. MOTS(JGA+ 7)='EPSR'
  574. MOTS(JGA+ 8)='GFT1'
  575. MOTS(JGA+ 9)='GFT2'
  576. MOTS(JGA+10)='GFT3'
  577. MOTS(JGA+11)='GS1 '
  578. MOTS(JGA+12)='GS2 '
  579. MOTS(JGA+13)='GS3 '
  580. MOTS(JGA+14)='LTR1'
  581. MOTS(JGA+15)='LTR2'
  582. MOTS(JGA+16)='LTR3'
  583. MOTS(JGA+17)='WRU1'
  584. MOTS(JGA+18)='WRU2'
  585. MOTS(JGA+19)='BIL1'
  586. MOTS(JGA+20)='BIL2'
  587. MOTS(JGA+21)='EPSB'
  588. MOTS(JGA+22)='SIMP'
  589. MOTS(JGA+23)='TYPE'
  590. MOTS(JGA+24)='COMP'
  591. MOTS(JGA+25)='DYN '
  592. MOTS(JGA+26)='DYN1'
  593. MOTS(JGA+27)='DYN2'
  594. MOTS(JGA+28)='DYN3'
  595. MOFAC=MLMOTS
  596. RETURN
  597. C= 19.2 - Modes de calcul PLAN CONT et TRIDimensionnel COQUEs minces
  598. ELSE IF ((IFOUR.EQ.2.AND.(MFR.EQ.3.OR.MFR.EQ.9)).OR.
  599. . IFOUR.EQ.-2) THEN
  600. MLMOTS=MOOBL
  601. JGA=MOTS(/2)
  602. JGM=JGA+39
  603. SEGADJ,MLMOTS
  604. MOOBL=MLMOTS
  605. MOTS(JGA+ 1)='R '
  606. MOTS(JGA+ 2)='DG '
  607. MOTS(JGA+ 3)='K1 '
  608. MOTS(JGA+ 4)='M1 '
  609. MOTS(JGA+ 5)='Q1 '
  610. MOTS(JGA+ 6)='N1 '
  611. MOTS(JGA+ 7)='K2 '
  612. MOTS(JGA+ 8)='M2 '
  613. MOTS(JGA+ 9)='Q2 '
  614. MOTS(JGA+10)='N2 '
  615. MOTS(JGA+11)='OMEG'
  616. MOTS(JGA+12)='Q '
  617. MOTS(JGA+13)='H '
  618. MOTS(JGA+14)='BETA'
  619. MOTS(JGA+15)='K '
  620. MOTS(JGA+16)='A '
  621. MOTS(JGA+17)='Q3 '
  622. MOTS(JGA+18)='N3 '
  623. MOTS(JGA+19)='ADEN'
  624. MOTS(JGA+20)='KGON'
  625. MOTS(JGA+21)='POR0'
  626. MOTS(JGA+22)='BUMI'
  627. MOTS(JGA+23)='EFIS'
  628. MOTS(JGA+24)='DG0 '
  629. MOTS(JGA+25)='CR1 '
  630. MOTS(JGA+26)='CR2 '
  631. MOTS(JGA+27)='CR3 '
  632. MOTS(JGA+28)='CR '
  633. MOTS(JGA+29)='KP '
  634. MOTS(JGA+30)='AP '
  635. MOTS(JGA+31)='BP '
  636. MOTS(JGA+32)='QP '
  637. IF (IFOUR.EQ.-2) THEN
  638. MOTS(JGA+33)='LXX '
  639. MOTS(JGA+34)='LYY '
  640. MOTS(JGA+35)='LXY '
  641. MOTS(JGA+36)='LZZ '
  642. MOTS(JGA+37)='PXX '
  643. MOTS(JGA+38)='PYY '
  644. MOTS(JGA+39)='PXY '
  645. C* ELSE IF (IFOUR.EQ.2.AND.(MFR.EQ.3.OR.MFR.EQ.9)) THEN
  646. ELSE
  647. MOTS(JGA+33)='LSS '
  648. MOTS(JGA+34)='LTT '
  649. MOTS(JGA+35)='LST '
  650. MOTS(JGA+36)='LNN '
  651. MOTS(JGA+37)='PSS '
  652. MOTS(JGA+38)='PTT '
  653. MOTS(JGA+39)='PST '
  654. ENDIF
  655. MLMOTS=MOFAC
  656. JGA=MOTS(/2)
  657. JGM=JGA+23
  658. SEGADJ,MLMOTS
  659. MOFAC=MLMOTS
  660. MOTS(JGA+ 1)='GFTR'
  661. MOTS(JGA+ 2)='LTR '
  662. MOTS(JGA+ 3)='GS '
  663. MOTS(JGA+ 4)='BTR '
  664. MOTS(JGA+ 5)='WRUP'
  665. MOTS(JGA+ 6)='BILI'
  666. MOTS(JGA+ 7)='GFT1'
  667. MOTS(JGA+ 8)='GFT2'
  668. MOTS(JGA+ 9)='GS1 '
  669. MOTS(JGA+10)='GS2 '
  670. MOTS(JGA+11)='LTR1'
  671. MOTS(JGA+12)='LTR2'
  672. MOTS(JGA+13)='WRU1'
  673. MOTS(JGA+14)='WRU2'
  674. MOTS(JGA+15)='BIL1'
  675. MOTS(JGA+16)='BIL2'
  676. MOTS(JGA+17)='SIMP'
  677. MOTS(JGA+18)='TYPE'
  678. MOTS(JGA+19)='COMP'
  679. MOTS(JGA+20)='DYN '
  680. MOTS(JGA+21)='DYN1'
  681. MOTS(JGA+22)='DYN2'
  682. MOTS(JGA+23)='DYN3'
  683. RETURN
  684. C= 19.3 - Formulation TRIDIMENSIONNELLE MASSIVE
  685. ELSE IF (IFOUR.EQ.2.AND.MFR.EQ.1) THEN
  686. MLMOTS=MOOBL
  687. JGA=MOTS(/2)
  688. JGM=JGA+44
  689. SEGADJ,MLMOTS
  690. MOTS(JGA+ 1)='R '
  691. MOTS(JGA+ 2)='DG '
  692. MOTS(JGA+ 3)='K1 '
  693. MOTS(JGA+ 4)='M1 '
  694. MOTS(JGA+ 5)='Q1 '
  695. MOTS(JGA+ 6)='N1 '
  696. MOTS(JGA+ 7)='K2 '
  697. MOTS(JGA+ 8)='M2 '
  698. MOTS(JGA+ 9)='Q2 '
  699. MOTS(JGA+10)='N2 '
  700. MOTS(JGA+11)='OMEG'
  701. MOTS(JGA+12)='Q '
  702. MOTS(JGA+13)='H '
  703. MOTS(JGA+14)='BETA'
  704. MOTS(JGA+15)='K '
  705. MOTS(JGA+16)='A '
  706. MOTS(JGA+17)='Q3 '
  707. MOTS(JGA+18)='N3 '
  708. MOTS(JGA+19)='ADEN'
  709. MOTS(JGA+20)='KGON'
  710. MOTS(JGA+21)='POR0'
  711. MOTS(JGA+22)='BUMI'
  712. MOTS(JGA+23)='EFIS'
  713. MOTS(JGA+24)='DG0 '
  714. MOTS(JGA+25)='CR1 '
  715. MOTS(JGA+26)='CR2 '
  716. MOTS(JGA+27)='CR3 '
  717. MOTS(JGA+28)='CR '
  718. MOTS(JGA+29)='KP '
  719. MOTS(JGA+30)='AP '
  720. MOTS(JGA+31)='BP '
  721. MOTS(JGA+32)='QP '
  722. MOTS(JGA+33)='LXX '
  723. MOTS(JGA+34)='LYY '
  724. MOTS(JGA+35)='LZZ '
  725. MOTS(JGA+36)='LXY '
  726. MOTS(JGA+37)='LXZ '
  727. MOTS(JGA+38)='LYZ '
  728. MOTS(JGA+39)='PXX '
  729. MOTS(JGA+40)='PYY '
  730. MOTS(JGA+41)='PZZ '
  731. MOTS(JGA+42)='PXY '
  732. MOTS(JGA+43)='PXZ '
  733. MOTS(JGA+44)='PYZ '
  734. MOOBL=MLMOTS
  735. MLMOTS=MOFAC
  736. JGA=MOTS(/2)
  737. JGM=JGA+28
  738. SEGADJ,MLMOTS
  739. MOTS(JGA+ 1)='GFTR'
  740. MOTS(JGA+ 2)='LTR '
  741. MOTS(JGA+ 3)='GS '
  742. MOTS(JGA+ 4)='BTR '
  743. MOTS(JGA+ 5)='WRUP'
  744. MOTS(JGA+ 6)='BILI'
  745. MOTS(JGA+ 7)='GFT1'
  746. MOTS(JGA+ 8)='GFT2'
  747. MOTS(JGA+ 9)='GFT3'
  748. MOTS(JGA+10)='GS1 '
  749. MOTS(JGA+11)='GS2 '
  750. MOTS(JGA+12)='GS3 '
  751. MOTS(JGA+13)='LTR1'
  752. MOTS(JGA+14)='LTR2'
  753. MOTS(JGA+15)='LTR3'
  754. MOTS(JGA+16)='WRU1'
  755. MOTS(JGA+17)='WRU2'
  756. MOTS(JGA+18)='WRU3'
  757. MOTS(JGA+19)='BIL1'
  758. MOTS(JGA+20)='BIL2'
  759. MOTS(JGA+21)='BIL3'
  760. MOTS(JGA+22)='SIMP'
  761. MOTS(JGA+23)='TYPE'
  762. MOTS(JGA+24)='COMP'
  763. MOTS(JGA+25)='DYN '
  764. MOTS(JGA+26)='DYN1'
  765. MOTS(JGA+27)='DYN2'
  766. MOTS(JGA+28)='DYN3'
  767. MOFAC=MLMOTS
  768. RETURN
  769. C= 19.4 - Formulation UNIDIMENSIONNELLE (1D) MASSIVE
  770. ELSE IF ((IFOUR.GE.3.AND.IFOUR.LE.15).AND.MFR.EQ.1) THEN
  771. MLMOTS=MOOBL
  772. JGA=MOTS(/2)
  773. JGM=JGA+37
  774. SEGADJ,MLMOTS
  775. MOTS(JGA+ 1)='R '
  776. MOTS(JGA+ 2)='DG '
  777. MOTS(JGA+ 3)='K1 '
  778. MOTS(JGA+ 4)='M1 '
  779. MOTS(JGA+ 5)='Q1 '
  780. MOTS(JGA+ 6)='N1 '
  781. MOTS(JGA+ 7)='K2 '
  782. MOTS(JGA+ 8)='M2 '
  783. MOTS(JGA+ 9)='Q2 '
  784. MOTS(JGA+10)='N2 '
  785. MOTS(JGA+11)='OMEG'
  786. MOTS(JGA+12)='Q '
  787. MOTS(JGA+13)='H '
  788. MOTS(JGA+14)='BETA'
  789. MOTS(JGA+15)='K '
  790. MOTS(JGA+16)='A '
  791. MOTS(JGA+17)='Q3 '
  792. MOTS(JGA+18)='N3 '
  793. MOTS(JGA+19)='ADEN'
  794. MOTS(JGA+20)='KGON'
  795. MOTS(JGA+21)='POR0'
  796. MOTS(JGA+22)='BUMI'
  797. MOTS(JGA+23)='EFIS'
  798. MOTS(JGA+24)='DG0 '
  799. MOTS(JGA+25)='CR1 '
  800. MOTS(JGA+26)='CR2 '
  801. MOTS(JGA+27)='CR3 '
  802. MOTS(JGA+28)='CR '
  803. MOTS(JGA+29)='KP '
  804. MOTS(JGA+30)='AP '
  805. MOTS(JGA+31)='BP '
  806. MOTS(JGA+32)='QP '
  807. C*OF A finir en fonction des definitions de ces grandeurs
  808. C*OF et en fonction de IFOUR
  809. IF (IFOUR.LE.11) THEN
  810. MOTS(JGA+33)='LXX '
  811. MOTS(JGA+34)='LYY '
  812. MOTS(JGA+35)='LZZ '
  813. MOTS(JGA+36)='PXX '
  814. MOTS(JGA+37)='PYY '
  815. ELSE
  816. MOTS(JGA+33)='LRR '
  817. MOTS(JGA+34)='LZZ '
  818. MOTS(JGA+35)='LOO '
  819. MOTS(JGA+36)='PRR '
  820. MOTS(JGA+37)='PZZ '
  821. ENDIF
  822. MOOBL=MLMOTS
  823. MLMOTS=MOFAC
  824. JGA=MOTS(/2)
  825. JGM=JGA+28
  826. C* A finir en fonction de la def. des parametres !
  827. SEGADJ,MLMOTS
  828. MOTS(JGA+ 1)='GFTR'
  829. MOTS(JGA+ 2)='LTR '
  830. MOTS(JGA+ 3)='GS '
  831. MOTS(JGA+ 4)='BTR '
  832. MOTS(JGA+ 5)='WRUP'
  833. MOTS(JGA+ 6)='BILI'
  834. MOTS(JGA+ 7)='GFT1'
  835. MOTS(JGA+ 8)='GFT2'
  836. MOTS(JGA+ 9)='GFT3'
  837. MOTS(JGA+10)='GS1 '
  838. MOTS(JGA+11)='GS2 '
  839. MOTS(JGA+12)='GS3 '
  840. MOTS(JGA+13)='LTR1'
  841. MOTS(JGA+14)='LTR2'
  842. MOTS(JGA+15)='LTR3'
  843. MOTS(JGA+16)='WRU1'
  844. MOTS(JGA+17)='WRU2'
  845. MOTS(JGA+18)='WRU3'
  846. MOTS(JGA+19)='BIL1'
  847. MOTS(JGA+20)='BIL2'
  848. MOTS(JGA+21)='BIL3'
  849. MOTS(JGA+22)='SIMP'
  850. MOTS(JGA+23)='TYPE'
  851. MOTS(JGA+24)='COMP'
  852. MOTS(JGA+25)='DYN '
  853. MOTS(JGA+26)='DYN1'
  854. MOTS(JGA+27)='DYN2'
  855. MOTS(JGA+28)='DYN3'
  856. MOFAC=MLMOTS
  857. RETURN
  858. ENDIF
  859. C ====
  860. C 20 - Modele viscoplastique BETON_URGC_ENDO
  861. C ====
  862. ELSE IF (IPLAC.EQ.20) THEN
  863. MLMOTS=MOOBL
  864. JGA=MOTS(/2)
  865. JGM=JGA+12
  866. SEGADJ,MLMOTS
  867. MOTS(JGA+ 1)='ALFA'
  868. MOTS(JGA+ 2)='BETA'
  869. MOTS(JGA+ 3)='LCS '
  870. MOTS(JGA+ 4)='GFC '
  871. MOTS(JGA+ 5)='GFT '
  872. MOTS(JGA+ 6)='ITON'
  873. MOTS(JGA+ 7)='IMOD'
  874. MOTS(JGA+ 8)='DC '
  875. MOTS(JGA+ 9)='DT '
  876. MOTS(JGA+10)='IASO'
  877. MOTS(JGA+11)='ICOU'
  878. MOTS(JGA+12)='ETA '
  879. MOOBL=MLMOTS
  880. RETURN
  881. C ====
  882. C 21 - Modele viscoplastique double endommagement
  883. C ====
  884. ELSE IF (IPLAC.EQ.21) THEN
  885. MLMOTS=MOOBL
  886. JGA=MOTS(/2)
  887. JGM=JGA+11
  888. SEGADJ,MLMOTS
  889. MOTS(JGA+1)='SIGY'
  890. MOTS(JGA+2)='B'
  891. MOTS(JGA+3)='RI'
  892. MOTS(JGA+4)='K'
  893. MOTS(JGA+5)='N'
  894. MOTS(JGA+6)='SD'
  895. MOTS(JGA+7)='RD'
  896. MOTS(JGA+8)='PD'
  897. MOTS(JGA+9)='SC'
  898. MOTS(JGA+10)='RC'
  899. MOTS(JGA+11)='PC'
  900. MOOBL=MLMOTS
  901. RETURN
  902. C ====
  903. C 22 - Modele viscoplastique CHAB_SINH_R
  904. C ====
  905. ELSE IF (IPLAC.EQ.22) THEN
  906. MLMOTS=MOOBL
  907. JGA=MOTS(/2)
  908. JGM=JGA+6
  909. SEGADJ,MLMOTS
  910. MOTS(JGA+1)='SIGY'
  911. MOTS(JGA+2)='B'
  912. MOTS(JGA+3)='RI'
  913. MOTS(JGA+4)='K'
  914. MOTS(JGA+5)='N'
  915. MOTS(JGA+6)='epso'
  916. MOOBL=MLMOTS
  917. RETURN
  918. C ====
  919. C 23 - Modele viscoplastique CHAB_SINH_X
  920. C ====
  921. ELSE IF (IPLAC.EQ.23) THEN
  922. MLMOTS=MOOBL
  923. JGA=MOTS(/2)
  924. JGM=JGA+9
  925. SEGADJ,MLMOTS
  926. MOTS(JGA+1)='SIGY'
  927. MOTS(JGA+2)='B'
  928. MOTS(JGA+3)='RI'
  929. MOTS(JGA+4)='K'
  930. MOTS(JGA+5)='N'
  931. MOTS(JGA+6)='epso'
  932. MOTS(JGA+7)='C1'
  933. MOTS(JGA+8)='gam1'
  934. MOTS(JGA+9)='C2'
  935. MOOBL=MLMOTS
  936. RETURN
  937. C ====
  938. C 24 - Modele viscoplastique CHAB_NOR_R
  939. C ====
  940. ELSE IF (IPLAC.EQ.24) THEN
  941. MLMOTS=MOOBL
  942. JGA=MOTS(/2)
  943. JGM=JGA+5
  944. SEGADJ,MLMOTS
  945. MOTS(JGA+1)='SIGY'
  946. MOTS(JGA+2)='B'
  947. MOTS(JGA+3)='RI'
  948. MOTS(JGA+4)='K'
  949. MOTS(JGA+5)='N'
  950. MOOBL=MLMOTS
  951. RETURN
  952. C ====
  953. C 25 - Modele viscoplastique CHAB_NOR_X
  954. C ====
  955. ELSE IF (IPLAC.EQ.25) THEN
  956. MLMOTS=MOOBL
  957. JGA=MOTS(/2)
  958. JGM=JGA+8
  959. SEGADJ,MLMOTS
  960. MOTS(JGA+1)='SIGY'
  961. MOTS(JGA+2)='B'
  962. MOTS(JGA+3)='RI'
  963. MOTS(JGA+4)='K'
  964. MOTS(JGA+5)='N'
  965. MOTS(JGA+6)='C1'
  966. MOTS(JGA+7)='gam1'
  967. MOTS(JGA+8)='C2'
  968. MOOBL=MLMOTS
  969. RETURN
  970. C ===
  971. C 26 - Modele d'ecrouissage-VISCOsite multiplicative avec enDOMMAGEment anisotrope
  972. C ===
  973. ELSE IF (IPLAC.EQ.26) THEN
  974. MLMOTS=MOOBL
  975. JGA=MOTS(/2)
  976. JGM=JGA+8
  977. SEGADJ,MLMOTS
  978. MOTS(JGA+1)='N '
  979. MOTS(JGA+2)='M '
  980. MOTS(JGA+3)='KK '
  981. MOTS(JGA+4)='ALF2'
  982. MOTS(JGA+5)='BET2'
  983. MOTS(JGA+6)='R '
  984. MOTS(JGA+7)='A '
  985. MOTS(JGA+8)='Q '
  986. MOOBL=MLMOTS
  987. MLMOTS=MOFAC
  988. JGA=MOTS(/2)
  989. JGM=JGA+1
  990. SEGADJ,MLMOTS
  991. MOTS(JGA+1)='SMAX'
  992. MOFAC=MLMOTS
  993. RETURN
  994. C ===
  995. C 27 - Modele FLUTRA
  996. C ===
  997. ELSE IF (IPLAC.EQ.27) THEN
  998. MLMOTS=MOOBL
  999. JGA=MOTS(/2)
  1000. JGM=JGA+24
  1001. SEGADJ,MLMOTS
  1002. MOTS(JGA+1)= 'KTR0'
  1003. MOTS(JGA+2)= 'ATRA'
  1004. MOTS(JGA+3)= 'BTRA'
  1005. MOTS(JGA+4)= 'ACOM'
  1006. MOTS(JGA+5)= 'BCOM'
  1007. MOTS(JGA+6)= 'EKV1'
  1008. MOTS(JGA+7)= 'TKV1'
  1009. MOTS(JGA+8)= 'EKV2'
  1010. MOTS(JGA+9)= 'TKV2'
  1011. MOTS(JGA+10)= 'EKV3'
  1012. MOTS(JGA+11)= 'TKV3'
  1013. MOTS(JGA+12)= 'EKV4'
  1014. MOTS(JGA+13)= 'TKV4'
  1015. MOTS(JGA+14)= 'EKV5'
  1016. MOTS(JGA+15)= 'TKV5'
  1017. MOTS(JGA+16)= 'EKV6'
  1018. MOTS(JGA+17)= 'TKV6'
  1019. MOTS(JGA+18)= 'BETA'
  1020. MOTS(JGA+19)= 'ETA1'
  1021. MOTS(JGA+20)= 'NUPO'
  1022. MOTS(JGA+21)= 'XFDE'
  1023. MOTS(JGA+22)= 'SLMA'
  1024. MOTS(JGA+23)= 'VGEA'
  1025. MOTS(JGA+24)= 'VGEB'
  1026. MOOBL=MLMOTS
  1027. RETURN
  1028. C ===
  1029. C 28 - Modele SYMONDS & COWPER SYCO1
  1030. C ===
  1031. ELSE IF (IPLAC.EQ.28) THEN
  1032. MLMOTS=MOOBL
  1033. JGA=MOTS(/2)
  1034. JGM=JGA+2
  1035. SEGADJ MLMOTS
  1036. MOTS(JGA+1)='PSYC'
  1037. MOTS(JGA+2)='DSYC'
  1038. MOOBL=MLMOTS
  1039. RETURN
  1040. C ===
  1041. C 29 - Modele SYMONDS & COWPER SYCO2
  1042. C ===
  1043. ELSE IF (IPLAC.EQ.29) THEN
  1044. MLMOTS=MOOBL
  1045. JGA=MOTS(/2)
  1046. JGM=JGA+4
  1047. SEGADJ MLMOTS
  1048. MOTS(JGA+1)='PSYC'
  1049. MOTS(JGA+2)='ASYC'
  1050. MOTS(JGA+3)='BSYC'
  1051. MOTS(JGA+4)='CSYC'
  1052. MOOBL=MLMOTS
  1053. RETURN
  1054. ENDIF
  1055. C ==========
  1056. C ERREUR : modele non reference (IRET=0)
  1057. C ==========
  1058. IRET=0
  1059. CALL ERREUR(261)
  1060.  
  1061. RETURN
  1062. END
  1063.  
  1064.  
  1065.  
  1066.  
  1067.  

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