Télécharger idmatr.eso

Retour à la liste

Numérotation des lignes :

idmatr
  1. C IDMATR SOURCE CB215821 24/04/12 21:16:21 11897
  2.  
  3. SUBROUTINE IDMATR(MFR,IPMODL,IPNOMC,NBROBL,NBRFAC)
  4.  
  5. *--------------------------------------------------------------------*
  6. * Noms de composantes de materiaux *
  7. *--------------------------------------------------------------------*
  8. * *
  9. * ENTREES: *
  10. * MFR Numero de formulation *
  11. * IPMODL objet modele elementaire ( segment actif ) *
  12. * *
  13. * SORTIES: *
  14. * IPNOMC pointeur sur les listes de noms de composantes *
  15. * obligatoires et facultatives *
  16. * NBROBL nombre de composantes obligatoires *
  17. * NBRFAC nombre de composantes facultatives *
  18. * *
  19. *--------------------------------------------------------------------*
  20. *
  21. IMPLICIT INTEGER(I-N)
  22. IMPLICIT REAL*8(A-H,O-Z)
  23.  
  24.  
  25. -INC PPARAM
  26. -INC CCOPTIO
  27. -INC SMLMOTS
  28. -INC SMMODEL
  29. POINTEUR NOMID1.NOMID
  30. *
  31. ******Modif SELLIER ****************************************************
  32. * augmentation nombre de caracteristiques
  33. * PARAMETER (ITA=100)
  34. * CHARACTER*16 MOMODL(100)
  35. PARAMETER (ITA=500)
  36. CHARACTER*16 MOMODL(500)
  37. integer I,J,K,L
  38. c nombre de variables et de parametres par variables
  39. c include './nombre_helmholtz.h'
  40. -INC HNBRHEL
  41. CHARACTER*1 motnl1,motnlv1,motnld1,motnla1
  42. CHARACTER*2 motnl2
  43. CHARACTER*3 motnl3
  44. CHARACTER*4 motnl4
  45. ******fin modif SELLIER ************************************************
  46.  
  47. CHARACTER*(LOCOMP) TABOBL(ITA),TABFAC(ITA)
  48. *
  49. IMODEL=IPMODL
  50. * Le segment existe-t-il deja?
  51. IF (lnomid(6).NE.0) THEN
  52. nomid = lnomid(6)
  53. SEGACT,nomid
  54. nbrobl = lesobl(/2)
  55. nbrfac = lesfac(/2)
  56. IF (nbrobl+nbrfac.EQ.0) GOTO 765
  57. SEGINI,nomid1=nomid
  58. ipnomc=nomid1
  59. RETURN
  60. ENDIF
  61. 765 CONTINUE
  62.  
  63. JGOBL=0
  64. JGFAC=0
  65. ipnomc=0
  66. IRET = 1
  67. *
  68. NMAT=MATMOD(/2)
  69. NFOR=FORMOD(/2)
  70. MELE=NEFMOD
  71.  
  72. *--------------------------------------------------------------------
  73. * CAS DE LA FORMULATION THERMOHYDRIQUE
  74. *--------------------------------------------------------------------
  75. CALL PLACE(FORMOD,NFOR,ithehy,'THERMOHYDRIQUE')
  76. IF (ithehy.ne.0) then
  77. if ( matmod(1).eq.'SCHREFLER') then
  78. JGOBL = 20
  79. TABOBL(1)='KGG'
  80. TABOBL(2)='KGC'
  81. TABOBL(3)='KGT'
  82. TABOBL(4)='KCG'
  83. TABOBL(5)='KCC'
  84. TABOBL(6)='KCT'
  85. TABOBL(7)='KTG'
  86. TABOBL(8)='KTC'
  87. TABOBL(9)='KTT'
  88. TABOBL(10)='CGG'
  89. TABOBL(11)='CGC'
  90. TABOBL(12)='CGT'
  91. TABOBL(13)='CCG'
  92. TABOBL(14)='CCC'
  93. TABOBL(15)='CCT'
  94. TABOBL(16)='CTG'
  95. TABOBL(17)='CTC'
  96. TABOBL(18)='CTT'
  97. TABOBL(19)='KTGG'
  98. TABOBL(20)='KTCG'
  99. else
  100. IRET = 0
  101. call erreur(5)
  102. endif
  103. GO TO 9999
  104. ENDIF
  105.  
  106. *--------------------------------------------------------------------
  107. * CAS DE LA FORMULATION LIQUIDE ET ELEMENT DE RACCORD LITU
  108. *--------------------------------------------------------------------
  109. CALL PLACE(FORMOD,NFOR,ILIQU,'LIQUIDE')
  110. IF (ILIQU.NE.0) THEN
  111. IF (MFR.NE.41.AND.MFR.NE.43) THEN
  112. JGOBL = 6
  113. TABOBL(1)='RHO '
  114. TABOBL(2)='CSON'
  115. TABOBL(3)='RORF'
  116. TABOBL(4)='CREF'
  117. TABOBL(5)='LCAR'
  118. TABOBL(6)='G '
  119.  
  120. ELSEIF (MFR.EQ.41) THEN
  121. JGOBL = 5
  122. TABOBL(1)='RHO '
  123. TABOBL(2)='CSON'
  124. TABOBL(3)='RORF'
  125. TABOBL(4)='CREF'
  126. TABOBL(5)='LCAR'
  127. c* ELSEIF (MFR.EQ.43) THEN
  128. ELSE
  129. JGOBL = 3
  130. TABOBL(1)='RHO '
  131. TABOBL(2)='LCAR'
  132. TABOBL(3)='RORF'
  133. ENDIF
  134. GOTO 9999
  135. ENDIF
  136.  
  137. *--------------------------------------------------------------------
  138. * CAS DE LA FORMULATION THERMIQUE
  139. * CAS DE LA FORMULATION DARCY
  140. *--------------------------------------------------------------------
  141. CALL PLACE(FORMOD,NFOR,ITHER ,'THERMIQUE')
  142. CALL PLACE(FORMOD,NFOR,IDARCY,'DARCY' )
  143. IF (ITHER.NE.0 .OR. IDARCY.NE.0) THEN
  144.  
  145. IF (ITHER.NE.0) THEN
  146.  
  147. C Cas particuliers de THERMIQUE CONVECTION :
  148. CALL PLACE(MATMOD,NMAT,ICONV,'CONVECTION')
  149. IF (ICONV.NE.0) THEN
  150. C IF (MFR.EQ.3.OR.MFR.EQ.5.OR.MFR.EQ.9) THEN
  151. CC Coques thermiques
  152. C JGOBL =2
  153. C TABOBL(1)='HINF'
  154. C TABOBL(2)='HSUP'
  155. C JGFAC =2
  156. C TABFAC(1)='TCINF'
  157. C TABFAC(2)='TCSUP'
  158. C else
  159. C Elements massifs
  160. JGOBL = 1
  161. TABOBL(1)='H'
  162. JGFAC=1
  163. TABFAC(1)='TC'
  164. C endif
  165. GOTO 9999
  166. ENDIF
  167.  
  168. C Cas particuliers de THERMIQUE RAYONNEMENT :
  169. CALL PLACE(MATMOD,NMAT,IRAYE,'RAYONNEMENT')
  170. IF (IRAYE.NE.0) THEN
  171. C Cas particuliers de THERMIQUE RAYONNEMENT :
  172. CALL PLACE(MATMOD,NMAT,ICAVE,'CAVITE')
  173. CALL PLACE(MATMOD,NMAT,IFACA,'FAC_A_FAC')
  174. CALL PLACE(MATMOD,NMAT,IINFI,'INFINI')
  175.  
  176. IF (ICAVE.NE.0) THEN
  177. C RAYONNEMENT en CAVITE :
  178. IF (MFR.EQ.3.OR.MFR.EQ.5.OR.MFR.EQ.9) THEN
  179. JGOBL=2
  180. TABOBL(1)='EINF'
  181. TABOBL(2)='ESUP'
  182. else
  183. JGOBL=1
  184. TABOBL(1)='EMIS'
  185. endif
  186.  
  187. JGFAC=3
  188. TABFAC(1)='CABS'
  189. TABFAC(2)='TABS'
  190. TABFAC(3)='H'
  191.  
  192. ELSEIF (IFACA.NE.0) THEN
  193. C RAYONNEMENT FACE_A_FACE :
  194. JGOBL=1
  195. TABOBL(1)='EMIS'
  196. JGFAC = 1
  197. TABFAC(1)='H'
  198.  
  199. ELSEIF (IINFI.NE.0) THEN
  200. C RAYONNEMENT a l'INFINI :
  201. IF (MFR.EQ.3.OR.MFR.EQ.5.OR.MFR.EQ.9) THEN
  202. JGOBL=3
  203. TABOBL(1)='EINF'
  204. TABOBL(2)='ESUP'
  205. TABOBL(3)='T_IN'
  206. ELSE
  207. JGOBL=2
  208. TABOBL(1)='EMIS'
  209. TABOBL(2)='T_IN'
  210. ENDIF
  211.  
  212. JGFAC = 2
  213. TABFAC(1)='E_IN'
  214. TABFAC(2)='H'
  215.  
  216. ELSE
  217. IRET = 0
  218. CALL ERREUR(5)
  219. ENDIF
  220. GOTO 9999
  221. ENDIF
  222.  
  223. C Cas particuliers de THERMIQUE SOURCE :
  224. ISRCE = 0
  225. CALL PLACE(MATMOD,NMAT,ISRCE,'SOURCE')
  226. IF (ISRCE.NE.0) THEN
  227. IF (INATUU.EQ.1.AND.IMATEE.EQ.1) THEN
  228. JGOBL = 1
  229. TABOBL(1)='QVOL'
  230. JGFAC = 2
  231. TABFAC(1)='QINF'
  232. TABFAC(2)='QSUP'
  233. GOTO 9999
  234. ELSEIF (INATUU.EQ.2.AND.IMATEE.EQ.1) THEN
  235. JGOBL = 3
  236. TABOBL(1)='QTOT'
  237. TABOBL(2)='ORIG'
  238. TABOBL(3)='RGAU'
  239. GOTO 9999
  240. ELSEIF (INATUU.EQ.2.AND.IMATEE.EQ.2) THEN
  241. JGOBL = 5
  242. TABOBL(1)='QTOT'
  243. TABOBL(2)='ORIG'
  244. TABOBL(3)='RGAU'
  245. TABOBL(4)='DIRE'
  246. TABOBL(5)='ZGAU'
  247. GOTO 9999
  248. ELSEIF (INATUU.EQ.2.AND.IMATEE.EQ.3) THEN
  249. JGOBL = 7
  250. TABOBL(1)='QTOT'
  251. TABOBL(2)='ORIG'
  252. TABOBL(3)='RGAU'
  253. TABOBL(4)='DIRE'
  254. TABOBL(5)='ZGAU'
  255. TABOBL(6)='DIRL'
  256. TABOBL(7)='LGAU'
  257. GOTO 9999
  258. ELSE
  259. c write(6,*) 'INATUU, IMATEE =',INATUU, IMATEE
  260. WRITE(6,*) ' Dans IDMATR : numero IMATEE non prevu'
  261. CALL ERREUR(21)
  262. RETURN
  263. ENDIF
  264. ENDIF
  265.  
  266. ENDIF
  267.  
  268. C* Cas THERMIQUE et DARCY
  269. CALL PLACE(MATMOD,NMAT,IORTH,'ORTHOTROPE')
  270. CALL PLACE(MATMOD,NMAT,IANIS,'ANISOTROPE')
  271. IF (IORTH.EQ.0.AND.IANIS.EQ.0)THEN
  272. JGOBL = 1
  273. IF(MFR .EQ. 75)THEN
  274. C Cas des JOI1 (MFR=75) ==> Ressorts THERMIQUES
  275. C ====================
  276. TABOBL(1)='KT'
  277. ELSE
  278. TABOBL(1)='K '
  279. ENDIF
  280.  
  281. ELSEIF (IORTH.NE.0) THEN
  282. * ELEMENTS COQUES
  283. IF (MFR.EQ.3.OR.MFR.EQ.5.OR.MFR.EQ.9) THEN
  284. JGOBL=5
  285. TABOBL(1)='K1 '
  286. TABOBL(2)='K2 '
  287. TABOBL(3)='K3 '
  288. TABOBL(4)='V1X '
  289. TABOBL(5)='V1Y '
  290.  
  291. ELSEIF (MFR.EQ.1.OR.MFR.EQ.45) THEN
  292. * ELEMENTS MASSIFS et HYBRIDES
  293. IF (IDIM.EQ.2) THEN
  294. IF (IFOMOD.EQ.1) THEN
  295. * ELEMENT MASSIF DE FOURIER
  296. JGOBL=5
  297. TABOBL(1)='K1 '
  298. TABOBL(2)='K2 '
  299. TABOBL(3)='V1X '
  300. TABOBL(4)='V1Y '
  301. TABOBL(5)='K3 '
  302.  
  303. ELSE
  304. JGOBL=4
  305. TABOBL(1)='K1 '
  306. TABOBL(2)='K2 '
  307. TABOBL(3)='V1X '
  308. TABOBL(4)='V1Y '
  309. ENDIF
  310.  
  311. ELSEIF (IDIM.EQ.3) THEN
  312. JGOBL=9
  313. TABOBL(1)='K1 '
  314. TABOBL(2)='K2 '
  315. TABOBL(3)='K3 '
  316. TABOBL(4)='V1X '
  317. TABOBL(5)='V1Y '
  318. TABOBL(6)='V1Z '
  319. TABOBL(7)='V2X '
  320. TABOBL(8)='V2Y '
  321. TABOBL(9)='V2Z '
  322. ENDIF
  323. ENDIF
  324.  
  325. ELSEIF (IANIS.NE.0) THEN
  326. * ELEMENTS MASSIFS
  327. IF (MFR.EQ.1.OR.MFR.EQ.45) THEN
  328. IF (IDIM.EQ.2) THEN
  329. IF (IFOMOD.EQ.1) THEN
  330. * ELEMENT MASSIF DE FOURIER
  331. JGOBL=6
  332. TABOBL(1)='K11 '
  333. TABOBL(2)='K22 '
  334. TABOBL(3)='K21 '
  335. TABOBL(4)='V1X '
  336. TABOBL(5)='V1Y '
  337. TABOBL(6)='K33 '
  338.  
  339. ELSE
  340. JGOBL=5
  341. TABOBL(1)='K11 '
  342. TABOBL(2)='K22 '
  343. TABOBL(3)='K21 '
  344. TABOBL(4)='V1X '
  345. TABOBL(5)='V1Y '
  346. ENDIF
  347.  
  348. ELSEIF (IDIM.EQ.3) THEN
  349. JGOBL=12
  350. TABOBL(1)='K11 '
  351. TABOBL(2)='K22 '
  352. TABOBL(3)='K33 '
  353. TABOBL(4)='K21 '
  354. TABOBL(5)='K31 '
  355. TABOBL(6)='K32 '
  356. TABOBL(7)='V1X '
  357. TABOBL(8)='V1Y '
  358. TABOBL(9)='V1Z '
  359. TABOBL(10)='V2X '
  360. TABOBL(11)='V2Y '
  361. TABOBL(12)='V2Z '
  362. ENDIF
  363. ENDIF
  364. ENDIF
  365. C*
  366. C* Cas THERMIQUE CONDUCTION, THERMIQUE PHASE ou THERMIQUE ADVECTION :
  367. IF (ITHER.NE.0) THEN
  368.  
  369. CALL PLACE(MATMOD,NMAT,IPHA,'PHASE')
  370. IF (IPHA.NE.0) THEN
  371. JGOBL0 = JGOBL
  372. JGOBL= JGOBL0+4
  373. TABOBL(JGOBL0+1)='RHO '
  374. TABOBL(JGOBL0+2)='C '
  375. TABOBL(JGOBL0+3)='QLAT'
  376. TABOBL(JGOBL0+4)='TPHA'
  377.  
  378. JGFAC=1
  379. TABFAC(1)='H '
  380. GOTO 9999
  381. ENDIF
  382.  
  383. CALL PLACE(MATMOD,NMAT,IADVE,'ADVECTION')
  384. if (iadve .ne. 0) then
  385. C Cas des Tuyaux 1D (MFR=79)
  386. C ===========================
  387. IF (MFR .EQ. 79) THEN
  388. JGFAC = 3
  389. TABFAC(1)='RHO '
  390. TABFAC(2)='C '
  391. TABFAC(3)='VITE'
  392.  
  393. C Cas Massif (MFR=1)
  394. C ===========================
  395. ELSEIF (MFR .EQ. 1) THEN
  396. C Bidimensionnel PLAN DPGE
  397. C ===========================
  398. IF (IFOUR.EQ.-3) THEN
  399. JGFAC = 5
  400. TABFAC(1)='RHO '
  401. TABFAC(2)='C '
  402. TABFAC(3)='VITX'
  403. TABFAC(4)='VITY'
  404. TABFAC(5)='VITZ'
  405.  
  406. C Bidimensionnel PLAN (CP/DP)
  407. C ===========================
  408. ELSEIF (IFOUR.EQ.-2 .OR. IFOUR.EQ.-1) THEN
  409. JGFAC = 4
  410. TABFAC(1)='RHO '
  411. TABFAC(2)='C '
  412. TABFAC(3)='VITX'
  413. TABFAC(4)='VITY'
  414.  
  415. C Axisymetrie
  416. C ===========================
  417. ELSEIF (IFOUR .EQ. 0) THEN
  418. JGFAC = 4
  419. TABFAC(1)='RHO '
  420. TABFAC(2)='C '
  421. TABFAC(3)='VITR'
  422. TABFAC(4)='VITZ'
  423.  
  424. C Fourier
  425. C ===========================
  426. ELSEIF (IFOUR .EQ. 1) THEN
  427. CALL ERREUR(21)
  428. RETURN
  429. C JGFAC = 5C
  430. C TABFAC(1)='RHO '
  431. C TABFAC(2)='C '
  432. C TABFAC(3)='VITR'
  433. C TABFAC(4)='VITZ'
  434. C TABFAC(5)='VITT'
  435.  
  436. C Tridimensionnel
  437. C ===========================
  438. ELSEIF (IFOUR .EQ. 2) THEN
  439. JGFAC = 5
  440. TABFAC(1)='RHO '
  441. TABFAC(2)='C '
  442. TABFAC(3)='VITX'
  443. TABFAC(4)='VITY'
  444. TABFAC(5)='VITZ'
  445.  
  446. C Unidimensionnel (1D)
  447. C ===========================
  448. ELSEIF (IFOUR.GE.3 .AND. IFOUR.LE.15) THEN
  449. IF (IFOUR.LE.6) THEN
  450. JGFAC = 3
  451. TABFAC(1)='RHO '
  452. TABFAC(2)='C '
  453. TABFAC(3)='VITX'
  454.  
  455. ELSEIF (IFOUR.EQ.7 .OR. IFOUR.EQ.8) THEN
  456. JGFAC = 4
  457. C Verifier l'utilite des composantes au dela de 4
  458. TABFAC(1)='RHO '
  459. TABFAC(2)='C '
  460. TABFAC(3)='VITX'
  461. TABFAC(4)='VITY'
  462.  
  463. ELSEIF (IFOUR.EQ.9 .OR. IFOUR.EQ.10) THEN
  464. JGFAC = 4
  465. C Verifier l'utilite des composantes au dela de 4
  466. TABFAC(1)='RHO '
  467. TABFAC(2)='C '
  468. TABFAC(3)='VITX'
  469. TABFAC(4)='VITZ'
  470.  
  471. ELSEIF (IFOUR.EQ.11) THEN
  472. JGFAC = 5
  473. C Verifier l'utilite des composantes au dela de 4
  474. TABFAC(1)='RHO '
  475. TABFAC(2)='C '
  476. TABFAC(3)='VITX'
  477. TABFAC(4)='VITY'
  478. TABFAC(5)='VITZ'
  479.  
  480. ELSEIF (IFOUR.EQ.12.OR.IFOUR.EQ.13.OR.IFOUR.EQ.15) THEN
  481. JGFAC = 3
  482. TABFAC(1)='RHO '
  483. TABFAC(2)='C '
  484. TABFAC(3)='VITR'
  485.  
  486. ELSEIF (IFOUR.EQ.14) THEN
  487. JGFAC = 4
  488. C Verifier l'utilite des composantes au dela de 4
  489. TABFAC(1)='RHO '
  490. TABFAC(2)='C '
  491. TABFAC(3)='VITR'
  492. TABFAC(4)='VITZ'
  493. ENDIF
  494. ENDIF
  495. ENDIF
  496. GOTO 9999
  497. endif
  498.  
  499. C composantes facultatives THERMIQUE CONDUCTION seule
  500. IF(MFR .EQ. 75)THEN
  501. C Cas des JOI1 (MFR=75) ==> Ressorts THERMIQUES
  502. C ====================
  503. JGFAC = 4
  504. TABFAC(1)='M'
  505. TABFAC(2)='C'
  506. TABFAC(3)='TINI'
  507. ELSE
  508. JGFAC = 4
  509. TABFAC(1)='RHO '
  510. TABFAC(2)='C '
  511. TABFAC(3)='H '
  512. TABFAC(4)='TINI'
  513. ENDIF
  514. ENDIF
  515. GOTO 9999
  516. ENDIF
  517.  
  518. *--------------------------------------------------------------------
  519. * CAS DE LA FORMULATION MECANIQUE
  520. *--------------------------------------------------------------------
  521. CALL PLACE(FORMOD,NFOR,IMECA,'MECANIQUE')
  522. IF (IMECA.NE.0) THEN
  523. *
  524. CALL MODLIN(MOMODL,NMOD)
  525. CALL PLACE(MOMODL,NMOD,IRET,MATMOD(1))
  526. IF (IRET.EQ.0) GOTO 9999
  527. IF (NMAT.GE.2) THEN
  528.  
  529. CALL MODELA(MOMODL,NMOD)
  530. CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(2))
  531. if ((IPLAC.EQ.9.OR.IPLAC.EQ.10).AND.NMAT.GT.2) then
  532. INMAT = 3
  533. goto 19
  534. endif
  535.  
  536. IF (IPLAC.NE.0) THEN
  537. INMAT=3
  538. IF (IPLAC.EQ.1) THEN
  539. INMAT=3
  540. GOTO 10
  541. ENDIF
  542. IF (IPLAC.EQ.2) THEN
  543. *
  544. * MATERIAU ELASTIQUE ORTHOTROPE
  545. *
  546. IF (MFR.EQ.75) THEN
  547. * JOINT UNIDIMENSIONNEL JOI1
  548. *
  549. IF(IFOUR.EQ.2)THEN
  550. JGOBL=12
  551. TABOBL(1)='V1X '
  552. TABOBL(2)='V1Y '
  553. TABOBL(3)='V1Z '
  554. TABOBL(4)='V2X '
  555. TABOBL(5)='V2Y '
  556. TABOBL(6)='V2Z '
  557. TABOBL(7)='KN '
  558. TABOBL(8)='KS1 '
  559. TABOBL(9)='KS2 '
  560. TABOBL(10)='QN '
  561. TABOBL(11)='QS1 '
  562. TABOBL(12)='QS2 '
  563. *
  564. JGFAC=10
  565. TABFAC(1)='MASS'
  566. TABFAC(2)='JX '
  567. TABFAC(3)='JY '
  568. TABFAC(4)='JZ '
  569. TABFAC(5)='ALPN'
  570. TABFAC(6)='ALP1'
  571. TABFAC(7)='ALP2'
  572. TABFAC(8)='ALQN'
  573. TABFAC(9)='ALQ1'
  574. TABFAC(10)='ALQ2'
  575. *
  576. ELSEIF(IFOUR.EQ.-2.OR.IFOUR.EQ.-1)THEN
  577. JGOBL=5
  578. TABOBL(1)='V1X '
  579. TABOBL(2)='V1Y '
  580. TABOBL(3)='KN '
  581. TABOBL(4)='KS '
  582. TABOBL(5)='QS '
  583. *
  584. JGFAC=6
  585. TABFAC(1)='MASS'
  586. TABFAC(2)='JZ'
  587. TABFAC(4)='ALPN'
  588. TABFAC(5)='ALPS'
  589. TABFAC(6)='ALQS'
  590. ENDIF
  591. *
  592. ELSEIF (MFR.EQ.3) THEN
  593. * COQUES MINCES
  594. *
  595. JGOBL=6
  596. TABOBL(1)='YG1 '
  597. TABOBL(2)='YG2 '
  598. TABOBL(3)='NU12'
  599. TABOBL(4)='G12 '
  600. TABOBL(5)='V1X '
  601. TABOBL(6)='V1Y '
  602. *
  603. IF(IFOUR.EQ.-2) THEN
  604. JGFAC=4
  605. ELSE
  606. JGFAC=3
  607. ENDIF
  608. TABFAC(1)='ALP1'
  609. TABFAC(2)='ALP2'
  610. TABFAC(3)='RHO '
  611. IF(IFOUR.EQ.-2) TABFAC(4)='DIM3'
  612.  
  613. ELSEIF (MFR.EQ.9.OR.MFR.EQ.5) THEN
  614. * COQUES AVEC CISAILLEMENT TRANSVERSE
  615. *
  616. JGOBL=8
  617. TABOBL(1)='YG1 '
  618. TABOBL(2)='YG2 '
  619. TABOBL(3)='NU12'
  620. TABOBL(4)='G12 '
  621. TABOBL(5)='G23 '
  622. TABOBL(6)='G13 '
  623. TABOBL(7)='V1X '
  624. TABOBL(8)='V1Y '
  625. *
  626. JGFAC=3
  627. TABFAC(1)='RHO '
  628. TABFAC(2)='ALP1'
  629. TABFAC(3)='ALP2'
  630.  
  631. ELSEIF (MFR.EQ.1.OR.MFR.EQ.31) THEN
  632. * ELEMENTS MASSIFS
  633. *
  634. IF(IDIM.EQ.3)THEN
  635. * ELEMENTS 3D
  636. JGOBL=15
  637. TABOBL(1)='YG1 '
  638. TABOBL(2)='YG2 '
  639. TABOBL(3)='YG3 '
  640. TABOBL(4)='NU12'
  641. TABOBL(5)='NU23'
  642. TABOBL(6)='NU13'
  643. TABOBL(7)='G12 '
  644. TABOBL(8)='G23 '
  645. TABOBL(9)='G13 '
  646. TABOBL(10)='V1X '
  647. TABOBL(11)='V1Y '
  648. TABOBL(12)='V1Z '
  649. TABOBL(13)='V2X '
  650. TABOBL(14)='V2Y '
  651. TABOBL(15)='V2Z '
  652. *
  653. JGFAC=4
  654. TABFAC(1)='RHO '
  655. TABFAC(2)='ALP1'
  656. TABFAC(3)='ALP2'
  657. TABFAC(4)='ALP3'
  658.  
  659. ELSEIF (IDIM.EQ.2) THEN
  660. IF(IFOUR.EQ.-2)THEN
  661. * CONTRAINTE PLANE
  662. JGOBL=9
  663. TABOBL(1)='YG1 '
  664. TABOBL(2)='YG2 '
  665. TABOBL(3)='NU12'
  666. TABOBL(4)='G12 '
  667. TABOBL(5)='V1X '
  668. TABOBL(6)='V1Y '
  669. TABOBL(7)='YG3 '
  670. TABOBL(8)='NU23'
  671. TABOBL(9)='NU13'
  672. *
  673. JGFAC=4
  674. TABFAC(1)='RHO '
  675. TABFAC(2)='ALP1'
  676. TABFAC(3)='ALP2'
  677. TABFAC(4)='DIM3'
  678.  
  679. ELSEIF (IFOUR.EQ.-1.OR.IFOUR.EQ.0.OR.IFOUR.EQ.-3) THEN
  680. * DEFORMATION PLANE ,AXISYMETRIE
  681. JGOBL=9
  682. TABOBL(1)='YG1 '
  683. TABOBL(2)='YG2 '
  684. TABOBL(3)='YG3 '
  685. TABOBL(4)='NU12'
  686. TABOBL(5)='NU23'
  687. TABOBL(6)='NU13'
  688. TABOBL(7)='G12 '
  689. TABOBL(8)='V1X '
  690. TABOBL(9)='V1Y '
  691. *
  692. JGFAC=4
  693. TABFAC(1)='RHO '
  694. TABFAC(2)='ALP1'
  695. TABFAC(3)='ALP2'
  696. TABFAC(4)='ALP3'
  697.  
  698. ELSEIF (IFOUR.EQ.1) THEN
  699. * AXISYMETRIE DE FOURIER
  700. JGOBL=15
  701. TABOBL(1)='YG1 '
  702. TABOBL(2)='YG2 '
  703. TABOBL(3)='YG3 '
  704. TABOBL(4)='NU12'
  705. TABOBL(5)='NU23'
  706. TABOBL(6)='NU13'
  707. TABOBL(7)='G12 '
  708. TABOBL(8)='G23 '
  709. TABOBL(9)='G13 '
  710. TABOBL(10)='V1X '
  711. TABOBL(11)='V1Y '
  712. TABOBL(12)='V1Z '
  713. TABOBL(13)='V2X '
  714. TABOBL(14)='V2Y '
  715. TABOBL(15)='V2Z '
  716. *
  717. JGFAC=4
  718. TABFAC(1)='RHO '
  719. TABFAC(2)='ALP1'
  720. TABFAC(3)='ALP2'
  721. TABFAC(4)='ALP3'
  722. ENDIF
  723. ELSEIF (IDIM.EQ.1) THEN
  724. C= Dans le cas UNID SPHErique, on doit avoir YG2=YG3 et NU12=NU13
  725. C= et dans le cas thermomecanique ALP2=ALP3.
  726. JGOBL=6
  727. TABOBL(1)='YG1 '
  728. TABOBL(2)='YG2 '
  729. TABOBL(3)='YG3 '
  730. TABOBL(4)='NU12'
  731. TABOBL(5)='NU23'
  732. TABOBL(6)='NU13'
  733.  
  734. IF (IFOUR.EQ.6) THEN
  735. JGFAC=2
  736.  
  737. ELSEIF (IFOUR.EQ.5 .OR. IFOUR.EQ.10) THEN
  738. JGFAC=3
  739. TABFAC(3)='ALP3'
  740. ELSEIF (IFOUR.EQ.4.OR.IFOUR.EQ.8.OR.
  741. . IFOUR.EQ.13) THEN
  742. JGFAC=3
  743. TABFAC(3)='ALP2'
  744. ELSE
  745. JGFAC=4
  746. TABFAC(3)='ALP2'
  747. TABFAC(4)='ALP3'
  748. ENDIF
  749. TABFAC(1)='RHO '
  750. TABFAC(2)='ALP1'
  751. ENDIF
  752.  
  753. ELSEIF (MFR.EQ.35) THEN
  754. * ELEMENTS JOINTS
  755. IF (IFOUR.EQ.2) THEN
  756. JGOBL=5
  757. TABOBL(1)='KS1 '
  758. TABOBL(2)='KS2 '
  759. TABOBL(3)='KN '
  760. TABOBL(4)='V1X '
  761. TABOBL(5)='V1Y '
  762. *
  763. JGFAC=2
  764. TABFAC(1)='RHO '
  765. TABFAC(2)='ALPN'
  766. ENDIF
  767. ENDIF
  768.  
  769. ELSEIF (IPLAC.EQ.3)THEN
  770. * MATERIAU ANISOTROPE ELASTIQUE
  771. *
  772. IF(MFR.EQ.75)THEN
  773. * JOINT UNIDIMESIONNEL JOI1
  774. *
  775. IF (IFOUR.EQ.2) THEN
  776. JGOBL=27
  777. TABOBL(1)='V1X '
  778. TABOBL(2)='V1Y '
  779. TABOBL(3)='V1Z '
  780. TABOBL(4)='V2X '
  781. TABOBL(5)='V2Y '
  782. TABOBL(6)='V2Z '
  783. TABOBL(7)='D11 '
  784. TABOBL(8)='D22 '
  785. TABOBL(9)='D33 '
  786. TABOBL(10)='D44 '
  787. TABOBL(11)='D55 '
  788. TABOBL(12)='D66 '
  789. TABOBL(13)='D21 '
  790. TABOBL(14)='D31 '
  791. TABOBL(15)='D32 '
  792. TABOBL(16)='D41 '
  793. TABOBL(17)='D42 '
  794. TABOBL(18)='D43 '
  795. TABOBL(19)='D51 '
  796. TABOBL(20)='D52 '
  797. TABOBL(21)='D53 '
  798. TABOBL(22)='D54 '
  799. TABOBL(23)='D61 '
  800. TABOBL(24)='D62 '
  801. TABOBL(25)='D63 '
  802. TABOBL(26)='D64 '
  803. TABOBL(27)='D65 '
  804. *
  805. JGFAC=10
  806. * MASS: masse totale de l'élément joint
  807. TABFAC(1)='MASS'
  808. TABFAC(2)='JX '
  809. TABFAC(3)='JY '
  810. TABFAC(4)='JZ '
  811. TABFAC(5)='ALP1'
  812. TABFAC(6)='ALP2'
  813. TABFAC(7)='ALP3'
  814. TABFAC(8)='ALQ1'
  815. TABFAC(9)='ALQ2'
  816. TABFAC(10)='ALQ3'
  817. *
  818. ELSEIF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1) THEN
  819. JGOBL=8
  820. TABOBL(1)='V1X '
  821. TABOBL(2)='V1Y '
  822. TABOBL(3)='D11 '
  823. TABOBL(4)='D22 '
  824. TABOBL(5)='D33 '
  825. TABOBL(6)='D21 '
  826. TABOBL(7)='D31 '
  827. TABOBL(8)='D32 '
  828. *
  829. JGFAC=5
  830. * MASS: masse totale de l'élément joint
  831. TABFAC(1)='MASS'
  832. TABFAC(2)='JZ '
  833. TABFAC(3)='ALP1'
  834. TABFAC(4)='ALP2'
  835. TABFAC(5)='ALQ3'
  836. ENDIF
  837. ENDIF
  838. *
  839. * ELEMENTS MASSIFS
  840. *
  841. IF(MFR.EQ.1.OR.MFR.EQ.31)THEN
  842. IF(IDIM.EQ.3)THEN
  843. * ELEMENTS 3D
  844. JGOBL=27
  845. TABOBL(1)='D11 '
  846. TABOBL(2)='D21 '
  847. TABOBL(3)='D22 '
  848. TABOBL(4)='D31 '
  849. TABOBL(5)='D32 '
  850. TABOBL(6)='D33 '
  851. TABOBL(7)='D41 '
  852. TABOBL(8)='D42 '
  853. TABOBL(9)='D43 '
  854. TABOBL(10)='D44 '
  855. TABOBL(11)='D51 '
  856. TABOBL(12)='D52 '
  857. TABOBL(13)='D53 '
  858. TABOBL(14)='D54 '
  859. TABOBL(15)='D55 '
  860. TABOBL(16)='D61 '
  861. TABOBL(17)='D62 '
  862. TABOBL(18)='D63 '
  863. TABOBL(19)='D64 '
  864. TABOBL(20)='D65 '
  865. TABOBL(21)='D66 '
  866. TABOBL(22)='V1X '
  867. TABOBL(23)='V1Y '
  868. TABOBL(24)='V1Z '
  869. TABOBL(25)='V2X '
  870. TABOBL(26)='V2Y '
  871. TABOBL(27)='V2Z '
  872. *
  873. JGFAC=7
  874. TABFAC(1)='RHO '
  875. TABFAC(2)='ALP1'
  876. TABFAC(3)='ALP2'
  877. TABFAC(4)='ALP3'
  878. TABFAC(5)='AL12'
  879. TABFAC(6)='AL13'
  880. TABFAC(7)='AL23'
  881.  
  882. ELSEIF (IDIM.EQ.2) THEN
  883. IF (IFOUR.EQ.-2) THEN
  884. * CONTRAINTE PLANE
  885. JGOBL=12
  886. TABOBL(1)='D11 '
  887. TABOBL(2)='D21 '
  888. TABOBL(3)='D22 '
  889. TABOBL(4)='D41 '
  890. TABOBL(5)='D42 '
  891. TABOBL(6)='D44 '
  892. TABOBL(7)='V1X '
  893. TABOBL(8)='V1Y '
  894. TABOBL(9)='D31 '
  895. TABOBL(10)='D32 '
  896. TABOBL(11)='D33 '
  897. TABOBL(12)='D43 '
  898. *
  899. JGFAC=5
  900. TABFAC(1)='RHO '
  901. TABFAC(2)='ALP1'
  902. TABFAC(3)='ALP2'
  903. TABFAC(4)='AL12'
  904. TABFAC(5)='DIM3'
  905.  
  906. ELSEIF (IFOUR.EQ.-1.OR.IFOUR.EQ.0.OR.IFOUR.EQ.-3) THEN
  907. * DEFORMATION PLANE ,AXISYMETRIE
  908. JGOBL=12
  909. TABOBL(1)='D11 '
  910. TABOBL(2)='D21 '
  911. TABOBL(3)='D22 '
  912. TABOBL(4)='D31 '
  913. TABOBL(5)='D32 '
  914. TABOBL(6)='D33 '
  915. TABOBL(7)='D41 '
  916. TABOBL(8)='D42 '
  917. TABOBL(9)='D43 '
  918. TABOBL(10)='D44 '
  919. TABOBL(11)='V1X '
  920. TABOBL(12)='V1Y '
  921. *
  922. JGFAC=5
  923. TABFAC(1)='RHO '
  924. TABFAC(2)='ALP1'
  925. TABFAC(3)='ALP2'
  926. TABFAC(4)='AL12'
  927. TABFAC(5)='ALP3'
  928.  
  929. ELSEIF (IFOUR.EQ.1) THEN
  930. * AXISYMETRIE DE FOURIER
  931. JGOBL=15
  932. TABOBL(1)='D11 '
  933. TABOBL(2)='D21 '
  934. TABOBL(3)='D22 '
  935. TABOBL(4)='D31 '
  936. TABOBL(5)='D32 '
  937. TABOBL(6)='D33 '
  938. TABOBL(7)='D41 '
  939. TABOBL(8)='D42 '
  940. TABOBL(9)='D43 '
  941. TABOBL(10)='D44 '
  942. TABOBL(11)='D55 '
  943. TABOBL(12)='D65 '
  944. TABOBL(13)='D66 '
  945. TABOBL(14)='V1X '
  946. TABOBL(15)='V1Y '
  947. *
  948. JGFAC=5
  949. TABFAC(1)='RHO '
  950. TABFAC(2)='ALP1'
  951. TABFAC(3)='ALP2'
  952. TABFAC(4)='AL12'
  953. TABFAC(5)='ALP3'
  954. ENDIF
  955. ENDIF
  956. ENDIF
  957.  
  958. ELSEIF (IPLAC.EQ.4) THEN
  959. IF (MFR.EQ.33) THEN
  960. * MILIEU POREUX ISOTROPE
  961. *
  962. JGOBL=12
  963. TABOBL(1) ='YOUN'
  964. TABOBL(2) ='NU '
  965. TABOBL(3) ='RHO '
  966. TABOBL(4) ='ALPH'
  967. TABOBL(5) ='COB '
  968. TABOBL(6) ='MOB '
  969. TABOBL(7) ='ALPM'
  970. TABOBL(8) ='PERM'
  971. TABOBL(9) ='VISC'
  972. TABOBL(10)='KF '
  973. TABOBL(11)='RHOF'
  974. TABOBL(12)='ALPF'
  975. *
  976. IF (IFOUR.EQ.-2) THEN
  977. JGFAC=1
  978. TABFAC(1)='DIM3'
  979. ENDIF
  980. ENDIF
  981.  
  982. ELSEIF (IPLAC.EQ.5) THEN
  983. IF (MFR.EQ.37) THEN
  984. * MILIEU HOMOGENEISE
  985. *
  986. JGOBL=16
  987. TABOBL( 1)='B11 '
  988. TABOBL( 2)='B22 '
  989. TABOBL( 3)='B12 '
  990. TABOBL( 4)='ROF '
  991. TABOBL( 5)='ROS '
  992. TABOBL( 6)='YOUN'
  993. TABOBL( 7)='CSON'
  994. TABOBL( 8)='RORF'
  995. TABOBL( 9)='CREF'
  996. TABOBL(10)='LCAR'
  997. TABOBL(11)='E111'
  998. TABOBL(12)='E112'
  999. TABOBL(13)='E121'
  1000. TABOBL(14)='E122'
  1001. TABOBL(15)='E221'
  1002. TABOBL(16)='E222'
  1003. ENDIF
  1004. ELSEIF (IPLAC.EQ.6) THEN
  1005. C
  1006. C MATERIAU ELASTIQUE UNIDIRECTIONNEL
  1007. C
  1008. IF ((MFR.EQ.1.OR.MFR.EQ.31).AND.IDIM.EQ.3) THEN
  1009. JGOBL=7
  1010. TABOBL(1)='YOUN'
  1011. TABOBL(2)='V1X '
  1012. TABOBL(3)='V1Y '
  1013. TABOBL(4)='V1Z '
  1014. TABOBL(5)='V2X '
  1015. TABOBL(6)='V2Y '
  1016. TABOBL(7)='V2Z '
  1017.  
  1018. ELSE
  1019. JGOBL=3
  1020. TABOBL(1)='YOUN'
  1021. TABOBL(2)='V1X '
  1022. TABOBL(3)='V1Y '
  1023. ENDIF
  1024.  
  1025. IF((MFR.EQ.1.OR.MFR.EQ.3.OR.MFR.EQ.31).AND.
  1026. . IFOUR.EQ.-2) THEN
  1027. JGFAC=3
  1028. TABFAC(3)='DIM3'
  1029. ELSE
  1030. JGFAC=2
  1031. ENDIF
  1032. TABFAC(1)='RHO '
  1033. TABFAC(2)='ALPH'
  1034.  
  1035. ELSEIF (IPLAC.EQ.7) THEN
  1036. C MODELE ET MATERIAU DE LA SECTION DU MODELE A FIBRE
  1037. C
  1038. JGOBL=2
  1039. TABOBL(1)='MODS'
  1040. TABOBL(2)='MATS'
  1041.  
  1042. JGFAC=1
  1043. TABFAC(1)='MANO'
  1044. C ---- IPLAC 8: Ajoute par Jiang, 22/08/1995 et fleuret 28/05/96
  1045.  
  1046. ELSEIF (IPLAC.EQ.8) THEN
  1047. C MODELE ET MATERIAU POUR LE CABLE PRECONTRAINT
  1048. C
  1049. JGOBL=1
  1050. TABOBL(1)='YOUN'
  1051.  
  1052. JGFAC=8
  1053. TABFAC(1)='FF'
  1054. TABFAC(2)='PHIF'
  1055. TABFAC(3)='GANC'
  1056. TABFAC(4)='RMU0'
  1057. TABFAC(5)='FPRG'
  1058. TABFAC(6)='RH10'
  1059. TABFAC(7)='ALPH'
  1060. TABFAC(8)='RHO '
  1061.  
  1062. ELSEIF (IPLAC.EQ.9) THEN
  1063. C MODAL
  1064. C
  1065. JGOBL=3
  1066. TABOBL(1)='FREQ'
  1067. TABOBL(2)='MASS'
  1068. TABOBL(3)='DEFO'
  1069.  
  1070. JGFAC = 9
  1071. TABFAC(1) = 'AMOR'
  1072. TABFAC(2) = 'CGRA'
  1073. TABFAC(3) = 'MADE'
  1074. TABFAC(4) = 'RICR'
  1075. TABFAC(5) = 'MAIB'
  1076. TABFAC(6) = 'MACR'
  1077. TABFAC(7) = 'AMCR'
  1078. TABFAC(8) = 'ALP0'
  1079. TABFAC(9) = 'ECRO'
  1080.  
  1081. ELSEIF (IPLAC.EQ.10) THEN
  1082. C STATIQUE
  1083. C
  1084. JGOBL=3
  1085. TABOBL(1)='DEFO'
  1086. TABOBL(2)='RIDE'
  1087. TABOBL(3)='MADE'
  1088.  
  1089. JGFAC = 8
  1090. TABFAC(1) = 'AMOR'
  1091. TABFAC(2) = 'RICR'
  1092. TABFAC(3) = 'MAIA'
  1093. TABFAC(4) = 'MAIB'
  1094. TABFAC(5) = 'MACR'
  1095. TABFAC(6) = 'AMCR'
  1096. TABFAC(7) = 'BET0'
  1097. TABFAC(8) = 'ECRO'
  1098.  
  1099. ELSEIF (IPLAC.EQ.11) THEN
  1100. C ZONE_COHESIVE
  1101. C
  1102. JGOBL=2
  1103. TABOBL(1)='KS'
  1104. TABOBL(2)='KN'
  1105. ENDIF
  1106.  
  1107. INMAT=3
  1108. GOTO 20
  1109.  
  1110. ELSE
  1111. INMAT=2
  1112. GOTO 10
  1113. ENDIF
  1114. ENDIF
  1115. INMAT=0
  1116. 10 CONTINUE
  1117. *
  1118. * MATERIAU ELASTIQUE ISOTROPE
  1119. *
  1120. IF (MFR.EQ.35) THEN
  1121. * VALABLE EN 2D COMME EN 3D
  1122. JGOBL=2
  1123. TABOBL(1)='KS '
  1124. TABOBL(2)='KN '
  1125.  
  1126. JGFAC=2
  1127. TABFAC(1)='RHO '
  1128. TABFAC(2)='ALPN'
  1129. *
  1130. ELSEIF (MFR.EQ.78) THEN
  1131. * VALABLE EN 2D COMME EN 3D
  1132. JGOBL=2
  1133. TABOBL(1)='KS '
  1134. TABOBL(2)='KN '
  1135. *
  1136. * JOINT CISAILLEMENT (2D)
  1137. *
  1138. ELSEIF (MFR.EQ.53) THEN
  1139. * VALABLE EN 2D "COMME EN 3D"
  1140.  
  1141. JGOBL=1
  1142. TABOBL(1)='KS '
  1143. *
  1144. JGFAC=2
  1145. TABFAC(1)='RHO '
  1146. TABFAC(2)='ALPN'
  1147. *
  1148. * TOUS LES AUTRES CAS
  1149. *
  1150. ELSE
  1151. JGOBL=2
  1152. TABOBL(1)='YOUN'
  1153. TABOBL(2)='NU '
  1154. *
  1155. IF((MFR.EQ.1.OR.MFR.EQ.3.OR.MFR.EQ.31.OR.MFR.EQ.63).AND.
  1156. . IFOUR.EQ.-2) THEN
  1157. JGFAC=4
  1158. TABFAC(3)='DIM3'
  1159. TABFAC(4)='VISQ'
  1160. ELSE
  1161. JGFAC=3
  1162. TABFAC(3)='VISQ'
  1163. ENDIF
  1164. TABFAC(1)='RHO '
  1165. TABFAC(2)='ALPH'
  1166.  
  1167. * AM CAS FORMULATION NON-LOCALE
  1168.  
  1169. INLOC = 0
  1170. MN3 = INFMOD(/1)
  1171. IF (MN3.GE.13) INLOC=-1*INFMOD(13)
  1172. IF (INLOC.GT.0) THEN
  1173. JGOBL0 = JGFAC
  1174. * moyenne
  1175. IF(INLOC.EQ.1) THEN
  1176. JGFAC=JGOBL0+1
  1177. TABFAC(JGOBL0+1)='LCAR'
  1178. * stress-based
  1179. ELSEIF(INLOC.EQ.2) THEN
  1180. JGFAC=JGOBL0+2
  1181. TABFAC(JGOBL0+1)='LCAR'
  1182. TABFAC(JGOBL0+2)='SBFT'
  1183.  
  1184. * helmholtz
  1185. ELSE IF(INLOC.EQ.3) THEN
  1186. * sellier 03/07/20
  1187. * JGFAC=JGOBL0+6
  1188. * TABFAC(JGOBL0+1)='LCAR'
  1189. * TABFAC(JGOBL0+2)='LCF1'
  1190. * TABFAC(JGOBL0+3)='LCF2'
  1191. * TABFAC(JGOBL0+4)='LCF3'
  1192. * TABFAC(JGOBL0+5)='LCF4'
  1193. * TABFAC(JGOBL0+6)='LCF5'
  1194. JGFAC=JGOBL0+NB_PARA_HELM
  1195. c variables de Helmholtz a diffusions ORThotropes
  1196. c modif liste para heilmholtz sellier 10//11/22
  1197. c print*,'Dans IDMATR ', NB_HELM, 'formulationss Helmholtz'
  1198. do I=1,NB_HELM
  1199. c numero de la variable de Helmholtz
  1200. write (motnl1,'(I1)') I
  1201. c Capacité (infinie si BH=1)
  1202. motnl3='CAP'
  1203. motnl4=motnl3//motnl1
  1204. TABFAC(JGOBL0+NB_PARA_PAR_HELM*(I-1)+1)=motnl4
  1205. c print*,motnl4
  1206. c Blocage de Dirichlet si 1 sinon 0
  1207. motnl3='BLO'
  1208. motnl4=motnl3//motnl1
  1209. TABFAC(JGOBL0+NB_PARA_PAR_HELM*(I-1)+2)=motnl4
  1210. c print*,motnl4
  1211. c Valeur imposee a la zone de blocage de Dirichlet
  1212. motnl3='DEP'
  1213. motnl4=motnl3//motnl1
  1214. TABFAC(JGOBL0+NB_PARA_PAR_HELM*(I-1)+3)=motnl4
  1215. c print*,motnl4
  1216. c Valeur initiale de la variable de Helmholtz
  1217. motnl3='INI'
  1218. motnl4=motnl3//motnl1
  1219. TABFAC(JGOBL0+NB_PARA_PAR_HELM*(I-1)+4)=motnl4
  1220. c print*,motnl4
  1221. c indicateur logique de lineratite pour ne pas passer par istep=3 si 1
  1222. motnl3='LIN'
  1223. motnl4=motnl3//motnl1
  1224. TABFAC(JGOBL0+NB_PARA_PAR_HELM*(I-1)+5)=motnl4
  1225. c print*,motnl4
  1226. do J=1,3
  1227. c coefficients de diffusion
  1228. write (motnld1,'(I1)') J
  1229. c diffusion
  1230. motnl2='DH'
  1231. c 1er indice: numero de la variable d helmholtz
  1232. motnl3=motnl2//motnl1
  1233. c 2eme indice : direction
  1234. motnl4=motnl3//motnld1
  1235. TABFAC(JGOBL0+(I-1)*NB_PARA_PAR_HELM+
  1236. # 5+(J-1)*4+1)=motnl4
  1237. c print*,motnl4
  1238. do K=1,3
  1239. write (motnla1,'(I1)') K
  1240. c direction principales de diffusion
  1241. motnlv1='V'
  1242. c indice numero de la variable
  1243. motnl2=motnlv1//motnl1
  1244. c indice direction principale
  1245. motnl3=motnl2//motnld1
  1246. c projection sur axe indice de base fixe
  1247. motnl4=motnl3//motnla1
  1248. TABFAC(JGOBL0+(I-1)*NB_PARA_PAR_HELM+
  1249. # 5+(J-1)*4+1+K)=motnl4
  1250. c print*,motnl4
  1251. end do
  1252. end do
  1253. end do
  1254. Cc on rajoute 2 inutiles pour compenser un decalage constatte dans castem20
  1255. C TABFAC(JGOBL0+NB_PARA_HELM+1)='FAC1'
  1256. C TABFAC(JGOBL0+NB_PARA_HELM+2)='FAC2'
  1257. * fin modif sellier
  1258. * cas non prevu
  1259. ELSE
  1260. IRET = 0
  1261. CALL ERREUR(5)
  1262. ENDIF
  1263. ENDIF
  1264. *
  1265. ENDIF
  1266. *
  1267. IF (INMAT.EQ.0) GOTO 9999
  1268. *
  1269. 19 CONTINUE
  1270. 20 CONTINUE
  1271. DO jm = 1,matmod(/2)
  1272. IF (matmod(jm).eq.'IMPEDANCE') THEN
  1273. imate = imatee
  1274. INMAT = INMAT + 1
  1275. JGOBL = 0
  1276. JGFAC = 0
  1277.  
  1278. IF(CMATEE.EQ.'IMPELAST') THEN
  1279. *IMPE_ELAS
  1280. JGOBL0 = JGOBL
  1281. JGOBL= JGOBL0+1
  1282. TABOBL(JGOBL0 + 1) ='RAID'
  1283. JGOBL0 = JGFAC
  1284.  
  1285. JGFAC = JGOBL0+4
  1286. TABFAC(JGOBL0+1) = 'AMOR'
  1287. TABFAC(JGOBL0+2) = 'ZNU'
  1288. TABFAC(JGOBL0+3) = 'MASS'
  1289. TABFAC(JGOBL0+4) = 'ALPH'
  1290. GOTO 22
  1291. ELSEIF (CMATEE.EQ.'IMPVOIGT'.or.CMATEE.EQ.'IMPREUSS') THEN
  1292. *IMPE_VOIGT ou IMPE_REUSS
  1293. JGOBL0 = JGOBL
  1294. JGOBL= JGOBL0+2
  1295. TABOBL(JGOBL0 + 1) ='RAID'
  1296. TABOBL(JGOBL0 + 2) ='VISC'
  1297.  
  1298. JGOBL0 = JGFAC
  1299. JGFAC = JGOBL0+2
  1300. TABFAC(JGOBL0+1) = 'MASS'
  1301. TABFAC(JGOBL0+2) = 'AMOR'
  1302. GOTO 22
  1303. ELSEIF (CMATEE.EQ.'IMPCOMPL') then
  1304. *IMPE_COMPLEXE
  1305. JGOBL0 = JGOBL
  1306. JGOBL= JGOBL0+1
  1307. TABOBL(JGOBL0 + 1) ='RAID'
  1308.  
  1309. JGOBL0 = JGFAC
  1310. JGFAC = JGFAC+4
  1311. TABFAC(JGOBL0+1) = 'MOCO'
  1312. TABFAC(JGOBL0+2) = 'VISC'
  1313. TABFAC(JGOBL0+3) = 'MASS'
  1314. TABFAC(JGOBL0+4) = 'AMOR'
  1315. GOTO 22
  1316. ELSE
  1317. ENDIF
  1318. 22 CONTINUE
  1319. * if (mele.eq.45) then
  1320. JGOBL0 = JGFAC
  1321. JGFAC = JGOBL0+3
  1322.  
  1323. IF(CMATEE.EQ.'IMPELAST') THEN
  1324. TABFAC(JGOBL0+1) = TABFAC(JGOBL0-2)
  1325. TABFAC(JGOBL0+2) = TABFAC(JGOBL0-1)
  1326. TABFAC(JGOBL0+3) = TABFAC(JGOBL0)
  1327. TABFAC(JGOBL0-2) = 'CPLE'
  1328. TABFAC(JGOBL0-1) = 'INER'
  1329. TABFAC(JGOBL0) = 'AROT'
  1330. ELSE
  1331. TABFAC(JGOBL0+1) = 'CPLE'
  1332. TABFAC(JGOBL0+2) = 'INER'
  1333. TABFAC(JGOBL0+3) = 'AROT'
  1334. ENDIF
  1335. * endif
  1336. IF (NMAT.GE.INMAT) inmat = inmat+1
  1337. * GOTO 9999
  1338. ENDIF
  1339. ENDDO
  1340. *
  1341. IF (NMAT.GE.INMAT) THEN
  1342. CALL MODNLI(MOMODL,NMOD)
  1343. CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  1344. IF (IPLAC.EQ.1) THEN
  1345. INMAT=INMAT+1
  1346. CALL MODPLA(MOMODL,NMOD)
  1347. CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  1348. IF (IPLAC.EQ.0) THEN
  1349. * ISOTROPE
  1350. IPLAC=1
  1351. * ELSE
  1352. * INMAT=INMAT+1
  1353. ENDIF
  1354. CALL IDPLAS(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  1355. ELSEIF (IPLAC.EQ.2) THEN
  1356. INMAT=INMAT+1
  1357. CALL MODFLU(MOMODL,NMOD)
  1358. CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  1359. IF (IPLAC.EQ.0) THEN
  1360. * NORTON
  1361. IPLAC=1
  1362. * ELSE
  1363. * INMAT=INMAT+1
  1364. ENDIF
  1365. CALL IDFLUA(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  1366. ELSEIF (IPLAC.EQ.3) THEN
  1367. INMAT=INMAT+1
  1368. CALL MODVIS(MOMODL,NMOD)
  1369. CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  1370. IF (IPLAC.EQ.0) THEN
  1371. * ONERA
  1372. IPLAC=2
  1373. * ELSE
  1374. * INMAT=INMAT+1
  1375. ENDIF
  1376. C
  1377. C Cas particulier des modeles GATT_MONERIE et UO2
  1378. C 'RHO ' et 'ALPH' sont obligatoires
  1379. IF (IPLAC.EQ.18.OR.IPLAC.EQ.19) THEN
  1380. JGOBL0 = JGOBL
  1381. JGOBL=JGOBL0+2
  1382. TABOBL(JGOBL0+1)='RHO '
  1383. TABOBL(JGOBL0+2)='ALPH'
  1384.  
  1385. JGOBL0=JGFAC
  1386. JGFAC=JGOBL0-2
  1387. IF (JGFAC.GT.0) THEN
  1388. DO 200 I=1,JGFAC
  1389. TABFAC(I)=TABFAC(I+2)
  1390. 200 CONTINUE
  1391. ENDIF
  1392.  
  1393. C ===
  1394. C Modeles SYMONDS & COWPER SYCO1 et SYCO2
  1395. C ===
  1396. ELSEIF (IPLAC.EQ.28) THEN
  1397. JGOBL0=JGOBL
  1398. JGOBL=JGOBL0+1
  1399. TABOBL(JGOBL0+1)='ECRO'
  1400.  
  1401. ELSEIF (IPLAC.EQ.29) THEN
  1402. JGOBL0=JGOBL
  1403. JGOBL=JGOBL0+1
  1404. TABOBL(JGOBL0+1)='ECRO'
  1405. ENDIF
  1406. CALL IDVISC(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  1407. C
  1408. ELSEIF (IPLAC.EQ.4) THEN
  1409. INMAT=INMAT+1
  1410. CALL MODEND(MOMODL,NMOD)
  1411. CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  1412. IF (IPLAC.EQ.0) THEN
  1413. * MAZARS
  1414. IPLAC=1
  1415. * ELSE
  1416. * INMAT=INMAT+1
  1417. ENDIF
  1418. CALL IDENDO(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  1419. ELSEIF (IPLAC.EQ.5) THEN
  1420. INMAT=INMAT+1
  1421. CALL MODPLE(MOMODL,NMOD)
  1422. CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  1423. IF (IPLAC.EQ.0) THEN
  1424. * TRIAXIAL P/Y
  1425. IPLAC=1
  1426. * ELSE
  1427. * INMAT=INMAT+1
  1428. ENDIF
  1429. CALL IDPLEN(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  1430. ELSEIF (IPLAC.EQ.6) THEN
  1431. INMAT=INMAT+1
  1432. CALL MODENL(MOMODL,NMOD)
  1433. CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  1434. IF (IPLAC.EQ.0) THEN
  1435. * ISOTROPE
  1436. IPLAC=1
  1437. * ELSE
  1438. * INMAT=INMAT+1
  1439. ENDIF
  1440. CALL IDELNL(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  1441. ELSEIF (IPLAC.EQ.7) THEN
  1442. C Pas de composantes supplementaires pour une loi 'VISCO_EXTERNE'
  1443. IRET = 1
  1444. ELSEIF(IPLAC.EQ.0.AND.IMATEE.EQ.31) THEN
  1445. GOTO 9999
  1446. ELSE
  1447. GOTO 9999
  1448. ENDIF
  1449. C IRET code retour de IDPLAS,IDFLUA,IDVISC,IDENDO,IDPLEN,IDELNL
  1450. IF (IRET.EQ.0) GOTO 9999
  1451. ENDIF
  1452.  
  1453. * Parametres CRIP et FUSION :
  1454. do jma=1,matmod(/2)
  1455. if(matmod(jma).eq.'CRIP') then
  1456. JGOBL = JGOBL + 1
  1457. TABOBL(JGOBL) = 'LIMP'
  1458.  
  1459. elseif(matmod(jma).eq.'FUSION') then
  1460. JGOBL = JGOBL + 1
  1461. TABOBL(JGOBL) = 'TFUS'
  1462.  
  1463. endif
  1464. enddo
  1465. *
  1466. GOTO 9999
  1467. ENDIF
  1468.  
  1469. *--------------------------------------------------------------------
  1470. * CAS DE LA FORMULATION POREUX
  1471. *--------------------------------------------------------------------
  1472. CALL PLACE(FORMOD,NFOR,IPORE,'POREUX')
  1473. IF (IPORE.NE.0) THEN
  1474. *
  1475. CALL MODLIN(MOMODL,NMOD)
  1476. CALL PLACE(MOMODL,NMOD,IRET,MATMOD(1))
  1477. IF (IRET.EQ.0) GOTO 9999
  1478. *
  1479. * D'ABORD : CAS NON ISOTROPE
  1480. *
  1481. IF (NMAT.GE.2) THEN
  1482. CALL MODELA(MOMODL,NMOD)
  1483. CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(2))
  1484. IF (IPLAC.NE.0) THEN
  1485. INMAT=3
  1486. IF (IPLAC.EQ.1) THEN
  1487. INMAT=3
  1488. GOTO 30
  1489. ENDIF
  1490. IF (IPLAC.EQ.2) THEN
  1491. *
  1492. * MATERIAU ELASTIQUE ORTHOTROPE
  1493. *
  1494. IF (MFR.EQ.33) THEN
  1495. *
  1496. * ELEMENTS MASSIFS
  1497. *
  1498. IF(IDIM.EQ.3)THEN
  1499. * ELEMENTS 3D
  1500. JGOBL=19
  1501. TABOBL(1)='YG1 '
  1502. TABOBL(2)='YG2 '
  1503. TABOBL(3)='YG3 '
  1504. TABOBL(4)='NU12'
  1505. TABOBL(5)='NU23'
  1506. TABOBL(6)='NU13'
  1507. TABOBL(7)='G12 '
  1508. TABOBL(8)='G23 '
  1509. TABOBL(9)='G13 '
  1510. TABOBL(10)='V1X '
  1511. TABOBL(11)='V1Y '
  1512. TABOBL(12)='V1Z '
  1513. TABOBL(13)='V2X '
  1514. TABOBL(14)='V2Y '
  1515. TABOBL(15)='V2Z '
  1516. TABOBL(16)='COB1'
  1517. TABOBL(17)='COB2'
  1518. TABOBL(18)='COB3'
  1519. TABOBL(19)='MOB '
  1520. *
  1521. JGFAC=12
  1522. TABFAC(1)='RHO '
  1523. TABFAC(2)='ALP1'
  1524. TABFAC(3)='ALP2'
  1525. TABFAC(4)='ALP3'
  1526. TABFAC(5)='ALPM'
  1527. TABFAC(6)='PER1'
  1528. TABFAC(7)='PER2'
  1529. TABFAC(8)='PER3'
  1530. TABFAC(9)='VISC'
  1531. TABFAC(10)='KF '
  1532. TABFAC(11)='RHOF'
  1533. TABFAC(12)='ALPF'
  1534. ELSEIF (IDIM.EQ.2) THEN
  1535. IF(IFOUR.EQ.-2)THEN
  1536. * CONTRAINTE PLANE
  1537. JGOBL=12
  1538. TABOBL(1)='YG1 '
  1539. TABOBL(2)='YG2 '
  1540. TABOBL(3)='NU12'
  1541. TABOBL(4)='G12 '
  1542. TABOBL(5)='V1X '
  1543. TABOBL(6)='V1Y '
  1544. TABOBL(7)='COB1'
  1545. TABOBL(8)='COB2'
  1546. TABOBL(9)='MOB '
  1547. TABOBL(10)='YG3 '
  1548. TABOBL(11)='NU23'
  1549. TABOBL(12)='NU13'
  1550. *
  1551. JGFAC=11
  1552. TABFAC(1)='RHO '
  1553. TABFAC(2)='ALP1'
  1554. TABFAC(3)='ALP2'
  1555. TABFAC(4)='ALPM'
  1556. TABFAC(5)='PER1'
  1557. TABFAC(6)='PER2'
  1558. TABFAC(7)='VISC'
  1559. TABFAC(8)='KF '
  1560. TABFAC(9)='RHOF'
  1561. TABFAC(10)='ALPF'
  1562. TABFAC(11)='DIM3'
  1563. ELSEIF (IFOUR.EQ.-1.OR.IFOUR.EQ.0.OR.IFOUR.EQ.-3) THEN
  1564. * DEFORMATION PLANE ,AXISYMETRIE
  1565. JGOBL=13
  1566. TABOBL(1)='YG1 '
  1567. TABOBL(2)='YG2 '
  1568. TABOBL(3)='YG3 '
  1569. TABOBL(4)='NU12'
  1570. TABOBL(5)='NU23'
  1571. TABOBL(6)='NU13'
  1572. TABOBL(7)='G12 '
  1573. TABOBL(8)='V1X '
  1574. TABOBL(9)='V1Y '
  1575. TABOBL(10)='COB1'
  1576. TABOBL(11)='COB2'
  1577. TABOBL(12)='COB3'
  1578. TABOBL(13)='MOB '
  1579. *
  1580. JGFAC=11
  1581. TABFAC(1)='RHO '
  1582. TABFAC(2)='ALP1'
  1583. TABFAC(3)='ALP2'
  1584. TABFAC(4)='ALP3'
  1585. TABFAC(5)='ALPM'
  1586. TABFAC(6)='PER1'
  1587. TABFAC(7)='PER2'
  1588. TABFAC(8)='VISC'
  1589. TABFAC(9)='KF '
  1590. TABFAC(10)='RHOF'
  1591. TABFAC(11)='ALPF'
  1592. *
  1593. ELSEIF (IFOUR.EQ.1) THEN
  1594. * AXISYMETRIE DE FOURIER
  1595. JGOBL=15
  1596. TABOBL(1)='YG1 '
  1597. TABOBL(2)='YG2 '
  1598. TABOBL(3)='YG3 '
  1599. TABOBL(4)='NU12'
  1600. TABOBL(5)='NU23'
  1601. TABOBL(6)='NU13'
  1602. TABOBL(7)='G12 '
  1603. TABOBL(8)='G23 '
  1604. TABOBL(9)='G13 '
  1605. TABOBL(10)='V1X '
  1606. TABOBL(11)='V1Y '
  1607. TABOBL(12)='COB1'
  1608. TABOBL(13)='COB2'
  1609. TABOBL(14)='COB3'
  1610. TABOBL(15)='MOB '
  1611. *
  1612. JGFAC=12
  1613. TABFAC(1)='RHO '
  1614. TABFAC(2)='ALP1'
  1615. TABFAC(3)='ALP2'
  1616. TABFAC(4)='ALP3'
  1617. TABFAC(5)='ALPM'
  1618. TABFAC(6)='PER1'
  1619. TABFAC(7)='PER2'
  1620. TABFAC(8)='PER3'
  1621. TABFAC(9)='VISC'
  1622. TABFAC(10)='KF '
  1623. TABFAC(11)='RHOF'
  1624. TABFAC(12)='ALPF'
  1625. ENDIF
  1626. ENDIF
  1627. ELSE
  1628. *
  1629. * CAS NON PREVU
  1630. *
  1631. IRET = 0
  1632. GOTO 9999
  1633. ENDIF
  1634. ELSEIF (IPLAC.EQ.3)THEN
  1635. *
  1636. * MATERIAU ANISOTROPE ELASTIQUE
  1637. *
  1638. IF(MFR.EQ.33)THEN
  1639. *
  1640. * ELEMENTS MASSIFS
  1641. *
  1642. IF(IDIM.EQ.3)THEN
  1643. * ELEMENTS 3D
  1644. JGOBL=34
  1645. TABOBL(1)='D11 '
  1646. TABOBL(2)='D21 '
  1647. TABOBL(3)='D22 '
  1648. TABOBL(4)='D31 '
  1649. TABOBL(5)='D32 '
  1650. TABOBL(6)='D33 '
  1651. TABOBL(7)='D41 '
  1652. TABOBL(8)='D42 '
  1653. TABOBL(9)='D43 '
  1654. TABOBL(10)='D44 '
  1655. TABOBL(11)='D51 '
  1656. TABOBL(12)='D52 '
  1657. TABOBL(13)='D53 '
  1658. TABOBL(14)='D54 '
  1659. TABOBL(15)='D55 '
  1660. TABOBL(16)='D61 '
  1661. TABOBL(17)='D62 '
  1662. TABOBL(18)='D63 '
  1663. TABOBL(19)='D64 '
  1664. TABOBL(20)='D65 '
  1665. TABOBL(21)='D66 '
  1666. TABOBL(22)='V1X '
  1667. TABOBL(23)='V1Y '
  1668. TABOBL(24)='V1Z '
  1669. TABOBL(25)='V2X '
  1670. TABOBL(26)='V2Y '
  1671. TABOBL(27)='V2Z '
  1672. TABOBL(28)='COB1'
  1673. TABOBL(29)='COB2'
  1674. TABOBL(30)='COB3'
  1675. TABOBL(31)='CO12'
  1676. TABOBL(32)='CO13'
  1677. TABOBL(33)='CO23'
  1678. TABOBL(34)='MOB '
  1679. *
  1680. JGFAC=18
  1681. TABFAC(1)='RHO '
  1682. TABFAC(2)='ALP1'
  1683. TABFAC(3)='ALP2'
  1684. TABFAC(4)='ALP3'
  1685. TABFAC(5)='AL12'
  1686. TABFAC(6)='AL13'
  1687. TABFAC(7)='AL23'
  1688. TABFAC(8)='ALPM'
  1689. TABFAC(9 )='PER1'
  1690. TABFAC(10)='PER2'
  1691. TABFAC(11)='PER3'
  1692. TABFAC(12)='PE12'
  1693. TABFAC(13)='PE13'
  1694. TABFAC(14)='PE23'
  1695. TABFAC(15)='VISC'
  1696. TABFAC(16)='KF '
  1697. TABFAC(17)='RHOF'
  1698. TABFAC(18)='ALPF'
  1699. ELSEIF (IDIM.EQ.2) THEN
  1700. IF (IFOUR.EQ.-2) THEN
  1701. * CONTRAINTE PLANE
  1702. JGOBL=16
  1703. TABOBL(1)='D11 '
  1704. TABOBL(2)='D21 '
  1705. TABOBL(3)='D22 '
  1706. TABOBL(4)='D41 '
  1707. TABOBL(5)='D42 '
  1708. TABOBL(6)='D44 '
  1709. TABOBL(7)='V1X '
  1710. TABOBL(8)='V1Y '
  1711. TABOBL(9 )='COB1'
  1712. TABOBL(10)='COB2'
  1713. TABOBL(11)='CO12'
  1714. TABOBL(12)='MOB '
  1715. TABOBL(13)='D31 '
  1716. TABOBL(14)='D32 '
  1717. TABOBL(15)='D33 '
  1718. TABOBL(16)='D43 '
  1719. *
  1720. JGFAC=13
  1721. TABFAC(1)='RHO '
  1722. TABFAC(2)='ALP1'
  1723. TABFAC(3)='ALP2'
  1724. TABFAC(4)='AL12'
  1725. TABFAC(5)='ALPM'
  1726. TABFAC(6)='PER1'
  1727. TABFAC(7)='PER2'
  1728. TABFAC(8)='PE12'
  1729. TABFAC(9)='VISC'
  1730. TABFAC(10)='KF '
  1731. TABFAC(11)='RHOF'
  1732. TABFAC(12)='ALPF'
  1733. TABFAC(13)='DIM3'
  1734. *
  1735. ELSEIF (IFOUR.EQ.-1.OR.IFOUR.EQ.0.OR.IFOUR.EQ.-3) THEN
  1736. * DEFORMATION PLANE ,AXISYMETRIE
  1737. JGOBL=17
  1738. TABOBL(1)='D11 '
  1739. TABOBL(2)='D21 '
  1740. TABOBL(3)='D22 '
  1741. TABOBL(4)='D31 '
  1742. TABOBL(5)='D32 '
  1743. TABOBL(6)='D33 '
  1744. TABOBL(7)='D41 '
  1745. TABOBL(8)='D42 '
  1746. TABOBL(9)='D43 '
  1747. TABOBL(10)='D44 '
  1748. TABOBL(11)='V1X '
  1749. TABOBL(12)='V1Y '
  1750. TABOBL(13)='COB1'
  1751. TABOBL(14)='COB2'
  1752. TABOBL(15)='CO12'
  1753. TABOBL(16)='COB3'
  1754. TABOBL(17)='MOB '
  1755. *
  1756. JGFAC=13
  1757. TABFAC(1)='RHO '
  1758. TABFAC(2)='ALP1'
  1759. TABFAC(3)='ALP2'
  1760. TABFAC(4)='AL12'
  1761. TABFAC(5)='ALP3'
  1762. TABFAC(6)='ALPM'
  1763. TABFAC(7)='PER1'
  1764. TABFAC(8)='PER2'
  1765. TABFAC(9)='PE12'
  1766. TABFAC(10)='VISC'
  1767. TABFAC(11)='KF '
  1768. TABFAC(12)='RHOF'
  1769. TABFAC(13)='ALPF'
  1770. *
  1771. ELSEIF (IFOUR.EQ.1) THEN
  1772. * AXISYMETRIE DE FOURIER
  1773. JGOBL=20
  1774. TABOBL(1)='D11 '
  1775. TABOBL(2)='D21 '
  1776. TABOBL(3)='D22 '
  1777. TABOBL(4)='D31 '
  1778. TABOBL(5)='D32 '
  1779. TABOBL(6)='D33 '
  1780. TABOBL(7)='D41 '
  1781. TABOBL(8)='D42 '
  1782. TABOBL(9)='D43 '
  1783. TABOBL(10)='D44 '
  1784. TABOBL(11)='D55 '
  1785. TABOBL(12)='D65 '
  1786. TABOBL(13)='D66 '
  1787. TABOBL(14)='V1X '
  1788. TABOBL(15)='V1Y '
  1789. TABOBL(16)='COB1'
  1790. TABOBL(17)='COB2'
  1791. TABOBL(18)='CO12'
  1792. TABOBL(19)='COB3'
  1793. TABOBL(20)='MOB '
  1794. *
  1795. JGFAC=14
  1796. TABFAC(1)='RHO '
  1797. TABFAC(2)='ALP1'
  1798. TABFAC(3)='ALP2'
  1799. TABFAC(4)='AL12'
  1800. TABFAC(5)='ALP3'
  1801. TABFAC(6)='ALPM'
  1802. TABFAC(7)='PER1'
  1803. TABFAC(8)='PER2'
  1804. TABFAC(9)='PE12'
  1805. TABFAC(10)='PER3'
  1806. TABFAC(11)='VISC'
  1807. TABFAC(12)='KF '
  1808. TABFAC(13)='RHOF'
  1809. TABFAC(14)='ALPF'
  1810. ENDIF
  1811. ENDIF
  1812. ELSE
  1813. *
  1814. * CAS NON PREVU
  1815. *
  1816. IRET = 0
  1817. GOTO 9999
  1818. ENDIF
  1819.  
  1820. ELSEIF (IPLAC.EQ.6) THEN
  1821. C
  1822. C MATERIAU ELASTIQUE UNIDIRECTIONNEL
  1823. C
  1824. IF(MFR.EQ.33)THEN
  1825. IF (IDIM.EQ.3) THEN
  1826. JGOBL=9
  1827. TABOBL(1)='YOUN'
  1828. TABOBL(2)='V1X '
  1829. TABOBL(3)='V1Y '
  1830. TABOBL(4)='V1Z '
  1831. TABOBL(5)='V2X '
  1832. TABOBL(6)='V2Y '
  1833. TABOBL(7)='V2Z '
  1834. TABOBL(8)='COB '
  1835. TABOBL(9)='MOB '
  1836. ELSE
  1837. JGOBL=5
  1838. TABOBL(1)='YOUN'
  1839. TABOBL(2)='V1X '
  1840. TABOBL(3)='V1Y '
  1841. TABOBL(4)='COB '
  1842. TABOBL(5)='MOB '
  1843. ENDIF
  1844. *
  1845. IF(IFOUR.EQ.-2) THEN
  1846. JGFAC=9
  1847. TABFAC(9)='DIM3'
  1848. ELSE
  1849. JGFAC=8
  1850. ENDIF
  1851. TABFAC(1)='RHO '
  1852. TABFAC(2)='ALPH'
  1853. TABFAC(3)='ALPM'
  1854. TABFAC(4)='PERM'
  1855. TABFAC(5)='VISC'
  1856. TABFAC(6)='KF '
  1857. TABFAC(7)='RHOF'
  1858. TABFAC(8)='ALPF'
  1859. ELSE
  1860. *
  1861. * CAS NON PREVU
  1862. *
  1863. IRET = 0
  1864. GO TO 9999
  1865. ENDIF
  1866. *
  1867. ENDIF
  1868. INMAT=3
  1869. GOTO 40
  1870. ELSE
  1871. INMAT=2
  1872. GOTO 30
  1873. ENDIF
  1874. ENDIF
  1875. INMAT=0
  1876. 30 CONTINUE
  1877. *
  1878. * CAS MATERIAU POREUX ELASTIQUE ISOTROPE
  1879. *
  1880. IF(MELE.GE.79.AND.MELE.LE.83)THEN
  1881. *
  1882. JGOBL=4
  1883. TABOBL(1) ='YOUN'
  1884. TABOBL(2) ='NU '
  1885. TABOBL(3) ='COB '
  1886. TABOBL(4) ='MOB '
  1887. *
  1888. IF(IFOUR.EQ.-2) THEN
  1889. JGFAC=9
  1890. TABFAC(9)='DIM3'
  1891. ELSE
  1892. JGFAC=8
  1893. ENDIF
  1894. TABFAC(1)='RHOF'
  1895. TABFAC(2)='ALPF'
  1896. TABFAC(3)='ALPM'
  1897. TABFAC(4)='PERM'
  1898. TABFAC(5)='VISC'
  1899. TABFAC(6)='KF '
  1900. TABFAC(7)='RHO '
  1901. TABFAC(8)='ALPH'
  1902. *
  1903. ELSEIF(MELE.GE.108.AND.MELE.LE.110)THEN
  1904. *
  1905. * CAS DES JOINTS POREUX ISOTROPES (VALABLE EN 2D COMME EN 3D)
  1906. *
  1907. JGOBL=4
  1908. TABOBL(1)='KS '
  1909. TABOBL(2)='KN '
  1910. TABOBL(3)='COB '
  1911. TABOBL(4)='MOB '
  1912. *
  1913. JGFAC=4
  1914. TABFAC(1)='PERT'
  1915. TABFAC(2)='PERH'
  1916. TABFAC(3)='PERB'
  1917. TABFAC(4)='VISC'
  1918. *
  1919. ELSEIF(MELE.GE.173.AND.MELE.LE.177)THEN
  1920. *
  1921. JGOBL=10
  1922. TABOBL(1) ='YOUN'
  1923. TABOBL(2) ='NU '
  1924. TABOBL(3) ='COP1'
  1925. TABOBL(4) ='COP2'
  1926. TABOBL(5) ='CPP1'
  1927. TABOBL(6) ='CPP2'
  1928. TABOBL(7) ='KK11'
  1929. TABOBL(8) ='KK12'
  1930. TABOBL(9) ='KK21'
  1931. TABOBL(10)='KK22'
  1932. *
  1933. IF(IFOUR.EQ.-2) THEN
  1934. JGFAC=9
  1935. TABFAC(9)='DIM3'
  1936. ELSE
  1937. JGFAC=8
  1938. ENDIF
  1939. TABFAC(1)='RHOF'
  1940. TABFAC(2)='ALPF'
  1941. TABFAC(3)='ALPM'
  1942. TABFAC(4)='PK11'
  1943. TABFAC(5)='PK12'
  1944. TABFAC(6)='PK21'
  1945. TABFAC(7)='PK22'
  1946. TABFAC(8)='ALPH'
  1947. *
  1948. ELSEIF(MELE.GE.185.AND.MELE.LE.187)THEN
  1949. *
  1950. * CAS DES JOINTS POREUX ISOTROPES (VALABLE EN 2D COMME EN 3D)
  1951. *
  1952. JGOBL=10
  1953. TABOBL(1)='KS '
  1954. TABOBL(2)='KN '
  1955. TABOBL(3)='COP1'
  1956. TABOBL(4)='COP2'
  1957. TABOBL(5)='CPP1'
  1958. TABOBL(6)='CPP2'
  1959. TABOBL(7)='KK11'
  1960. TABOBL(8)='KK12'
  1961. TABOBL(9)='KK21'
  1962. TABOBL(10)='KK22'
  1963. *
  1964. IF(IFOUR.EQ.-2) THEN
  1965. JGFAC=18
  1966. TABFAC(18)='DIM3'
  1967. ELSE
  1968. JGFAC=17
  1969. ENDIF
  1970. TABFAC(1)='RHOF'
  1971. TABFAC(2)='ALPF'
  1972. TABFAC(3)='ALPM'
  1973. TABFAC(4)='PT11'
  1974. TABFAC(5)='PH11'
  1975. TABFAC(6)='PB11'
  1976. TABFAC(7)='PT12'
  1977. TABFAC(8)='PH12'
  1978. TABFAC(9)='PB12'
  1979. TABFAC(10)='PT21'
  1980. TABFAC(11)='PH21'
  1981. TABFAC(12)='PB21'
  1982. TABFAC(13)='PT22'
  1983. TABFAC(14)='PH22'
  1984. TABFAC(15)='PB22'
  1985. TABFAC(16)='RHO '
  1986. TABFAC(17)='ALPH'
  1987.  
  1988. ELSEIF(MELE.GE.178.AND.MELE.LE.182)THEN
  1989. *
  1990. JGOBL=17
  1991. TABOBL(1)='YOUN'
  1992. TABOBL(2)='NU '
  1993. TABOBL(3)='COP1'
  1994. TABOBL(4)='COP2'
  1995. TABOBL(5)='COP3'
  1996. TABOBL(6)='CPP1'
  1997. TABOBL(7)='CPP2'
  1998. TABOBL(8)='CPP3'
  1999. TABOBL(9)='KK11'
  2000. TABOBL(10)='KK12'
  2001. TABOBL(11)='KK13'
  2002. TABOBL(12)='KK21'
  2003. TABOBL(13)='KK22'
  2004. TABOBL(14)='KK23'
  2005. TABOBL(15)='KK31'
  2006. TABOBL(16)='KK32'
  2007. TABOBL(17)='KK33'
  2008. *
  2009. IF(IFOUR.EQ.-2) THEN
  2010. JGFAC=15
  2011. TABFAC(15)='DIM3'
  2012. ELSE
  2013. JGFAC=14
  2014. ENDIF
  2015. TABFAC(1)='RHOF'
  2016. TABFAC(2)='ALPF'
  2017. TABFAC(3)='ALPM'
  2018. TABFAC(4)='PK11'
  2019. TABFAC(5)='PK12'
  2020. TABFAC(6)='PK13'
  2021. TABFAC(7)='PK21'
  2022. TABFAC(8)='PK22'
  2023. TABFAC(9)='PK23'
  2024. TABFAC(10)='PK31'
  2025. TABFAC(11)='PK32'
  2026. TABFAC(12)='PK33'
  2027. TABFAC(13)='RHO '
  2028. TABFAC(14)='ALPH'
  2029.  
  2030. ELSEIF(MELE.GE.188.AND.MELE.LE.190)THEN
  2031. *
  2032. * CAS DES JOINTS POREUX ISOTROPES (VALABLE EN 2D COMME EN 3D)
  2033. *
  2034. JGOBL=17
  2035. TABOBL(1)='KS '
  2036. TABOBL(2)='KN '
  2037. TABOBL(3)='COP1'
  2038. TABOBL(4)='COP2'
  2039. TABOBL(5)='COP3'
  2040. TABOBL(6)='CPP1'
  2041. TABOBL(7)='CPP2'
  2042. TABOBL(8)='CPP3'
  2043. TABOBL(9)='KK11'
  2044. TABOBL(10)='KK12'
  2045. TABOBL(11)='KK13'
  2046. TABOBL(12)='KK21'
  2047. TABOBL(13)='KK22'
  2048. TABOBL(14)='KK23'
  2049. TABOBL(15)='KK31'
  2050. TABOBL(16)='KK32'
  2051. TABOBL(17)='KK33'
  2052. *
  2053. IF(IFOUR.EQ.-2) THEN
  2054. JGFAC=33
  2055. TABFAC(33)='DIM3'
  2056. ELSE
  2057. JGFAC=32
  2058. ENDIF
  2059. TABFAC(1)='RHOF'
  2060. TABFAC(2)='ALPF'
  2061. TABFAC(3)='ALPM'
  2062. TABFAC(4)='PT11'
  2063. TABFAC(5)='PH11'
  2064. TABFAC(6)='PB11'
  2065. TABFAC(7)='PT12'
  2066. TABFAC(8)='PH12'
  2067. TABFAC(9)='PB12'
  2068. TABFAC(10)='PT13'
  2069. TABFAC(11)='PH13'
  2070. TABFAC(12)='PB13'
  2071. TABFAC(13)='PT21'
  2072. TABFAC(14)='PH21'
  2073. TABFAC(15)='PB21'
  2074. TABFAC(16)='PT22'
  2075. TABFAC(17)='PH22'
  2076. TABFAC(18)='PB22'
  2077. TABFAC(19)='PT23'
  2078. TABFAC(20)='PH23'
  2079. TABFAC(21)='PB23'
  2080. TABFAC(22)='PT31'
  2081. TABFAC(23)='PH31'
  2082. TABFAC(24)='PB31'
  2083. TABFAC(25)='PT32'
  2084. TABFAC(26)='PH32'
  2085. TABFAC(27)='PB32'
  2086. TABFAC(28)='PT33'
  2087. TABFAC(29)='PH33'
  2088. TABFAC(30)='PB33'
  2089. TABFAC(31)='RHO '
  2090. TABFAC(32)='ALPH'
  2091.  
  2092. ENDIF
  2093. *
  2094. IF (INMAT.EQ.0) THEN
  2095. ** IRET = 0
  2096. GOTO 9999
  2097. ENDIF
  2098. *
  2099. 40 CONTINUE
  2100. IF (NMAT.GE.INMAT) THEN
  2101. CALL MODNLI(MOMODL,NMOD)
  2102. CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  2103. IF (IPLAC.EQ.1) THEN
  2104. INMAT=INMAT+1
  2105. CALL MODPLA(MOMODL,NMOD)
  2106. CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  2107. IF (IPLAC.EQ.0) THEN
  2108. * ISOTROPE
  2109. IPLAC=1
  2110. * ELSE
  2111. * INMAT=INMAT+1
  2112. ENDIF
  2113. CALL IDPLAS(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  2114. IF (IRET.EQ.0) GOTO 9999
  2115. ELSEIF (IPLAC.EQ.2) THEN
  2116. INMAT=INMAT+1
  2117. CALL MODFLU(MOMODL,NMOD)
  2118. CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  2119. IF (IPLAC.EQ.0) THEN
  2120. * NORTON
  2121. IPLAC=1
  2122. * ELSE
  2123. * INMAT=INMAT+1
  2124. ENDIF
  2125. CALL IDFLUA(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  2126. IF (IRET.EQ.0) GOTO 9999
  2127. ELSEIF (IPLAC.EQ.3) THEN
  2128. INMAT=INMAT+1
  2129. CALL MODVIS(MOMODL,NMOD)
  2130. CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  2131. IF (IPLAC.EQ.0) THEN
  2132. * ONERA
  2133. IPLAC=2
  2134. * ELSE
  2135. * INMAT=INMAT+1
  2136. ENDIF
  2137. CALL IDVISC(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  2138. IF (IRET.EQ.0) GOTO 9999
  2139. ELSEIF (IPLAC.EQ.4) THEN
  2140. INMAT=INMAT+1
  2141. CALL MODEND(MOMODL,NMOD)
  2142. CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  2143. IF (IPLAC.EQ.0) THEN
  2144. * MAZARS
  2145. IPLAC=1
  2146. * ELSE
  2147. * INMAT=INMAT+1
  2148. ENDIF
  2149. CALL IDENDO(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  2150. IF (IRET.EQ.0) GOTO 9999
  2151. ELSE
  2152. * GOTO 9999
  2153. ENDIF
  2154. *
  2155. * En cas de creation de materiaux combinant plusieurs materiaux
  2156. * deja existant Ex ELASTIQUE ISOTROPE PLASTIQUE PARFAIT FLUAGE N
  2157. * GOTO 40
  2158. *
  2159. * GOTO 9999
  2160. ENDIF
  2161. GOTO 9999
  2162. ENDIF
  2163.  
  2164. *--------------------------------------------------------------------
  2165. * CAS DE LA FORMULATION CONTACT
  2166. *--------------------------------------------------------------------
  2167. CALL PLACE(FORMOD,NFOR,ICONT,'CONTACT')
  2168. IF (ICONT.NE.0) THEN
  2169. CALL MODFRO(MOMODL,NMOD)
  2170. iplla =0
  2171. do iou = 1,NMAT
  2172. CALL PLACE(MOMODL(1),5,IPLAC,MATMOD(iou))
  2173. iplla = max(iplac,iplla)
  2174. enddo
  2175. IF (iplla.EQ.0) THEN
  2176. IRET = 0
  2177. GOTO 9999
  2178. ENDIF
  2179.  
  2180. IF ((IPLLA.EQ.1) .OR. (IPLLA.EQ.2) .OR. (IPLLA.EQ.3)) THEN
  2181. * CONTACT SIMPLE
  2182. * --------------
  2183. INMAT=1
  2184. JGOBL=0
  2185.  
  2186. JGFAC=2
  2187. TABFAC(1)='JEU'
  2188. TABFAC(2)='ADHE'
  2189.  
  2190. ELSEIF (IPLLA .EQ. 4) THEN
  2191. * FROTTEMENT DE COULOMB
  2192. * ---------------------
  2193. INMAT=1
  2194. JGOBL=1
  2195. TABOBL(1)='MU '
  2196.  
  2197. JGFAC=3
  2198. TABFAC(1)='COHE'
  2199. TABFAC(2)='ADHE'
  2200. TABFAC(3)='JEU'
  2201.  
  2202. ELSEIF (IPLLA .EQ. 5) THEN
  2203. * FROTTEMENTS DE CABLES
  2204. * ---------------------
  2205. INMAT=1
  2206. JGOBL=2
  2207. TABOBL(1)='FF '
  2208. TABOBL(2)='PHIF'
  2209.  
  2210. JGFAC=0
  2211.  
  2212. ELSE
  2213. IRET=0
  2214. CALL ERREUR (261)
  2215. RETURN
  2216. ENDIF
  2217.  
  2218. GOTO 9999
  2219. ENDIF
  2220.  
  2221. *--------------------------------------------------------------------
  2222. * CAS DE LA FORMULATION CONTRAINTE
  2223. *--------------------------------------------------------------------
  2224. CALL PLACE(FORMOD,NFOR,ICNTR,'CONTRAINTE')
  2225. IF (ICNTR.NE.0) THEN
  2226. CALL MODCLI(MOMODL,NMOD)
  2227. iplla =0
  2228. do iou = 1,NMAT
  2229. CALL PLACE(MOMODL(1),NMOD,IPLAC,MATMOD(iou))
  2230. iplla = max(iplac,iplla)
  2231. enddo
  2232. IF (iplla.EQ.0) THEN
  2233. IRET = 0
  2234. GOTO 9999
  2235. ENDIF
  2236.  
  2237.  
  2238. IF (IPLLA .EQ. 1) THEN
  2239. * ROTATION
  2240. * ---------------------
  2241. INMAT=1
  2242. JGOBL=1
  2243. TABOBL(1)='ANGL'
  2244. ELSEIF (IPLLA .EQ. 2) THEN
  2245. * DEPLACEMENT
  2246. * ---------------------
  2247. INMAT=1
  2248. JGOBL=1
  2249. TABOBL(1)='AMPL'
  2250. ENDIF
  2251. ENDIF
  2252. *--------------------------------------------------------------------
  2253. * CAS DE LA FORMULATION MAGNETODYNAMIQUE
  2254. *--------------------------------------------------------------------
  2255. CALL PLACE(FORMOD,NFOR,ICONV,'MAGNETODYNAMIQUE')
  2256. IF (ICONV.NE.0) THEN
  2257. *
  2258. * FORMULATION EN COQUES
  2259. IF(MFR.EQ.3) THEN
  2260. CALL PLACE(MATMOD,NMAT,ISOT,'ISOTROPE')
  2261. IF(ISOT.NE.0) THEN
  2262. JGOBL=3
  2263.  
  2264. TABOBL(1)='ETA'
  2265. TABOBL(2)='PERM'
  2266. TABOBL(3)='EPAI'
  2267. ELSE
  2268. CALL PLACE(MATMOD,NMAT,IORTH,'ORTHOTROPE')
  2269. IF(IORTH.NE.0) THEN
  2270. JGOBL=4
  2271.  
  2272. TABOBL(1)='ETA1'
  2273. TABOBL(2)='ETA2'
  2274. TABOBL(3)='PERM'
  2275. TABOBL(4)='EPAI'
  2276. ENDIF
  2277. ENDIF
  2278. ENDIF
  2279. *
  2280. GOTO 9999
  2281. ENDIF
  2282.  
  2283. *--------------------------------------------------------------------
  2284. * CAS DE LA FORMULATION FISSURE
  2285. *--------------------------------------------------------------------
  2286. CALL PLACE(FORMOD,NFOR,ICONV,'FISSURE')
  2287. IF (ICONV.NE.0) THEN
  2288. *
  2289. * si POISEU_BLASIUS ou POISEU_COLEBROOK ou par defaut
  2290. JGOBL=1
  2291.  
  2292. TABOBL(1)='RUGO'
  2293. * si FROTTEMENT1 ou FROTTEMENT2
  2294. CALL PLACE(MATMOD,NMAT,IFT1,'FROTTEMENT1')
  2295. CALL PLACE(MATMOD,NMAT,IFT2,'FROTTEMENT2')
  2296. IF(IFT1.NE.0.OR.IFT2.NE.0) THEN
  2297. JGOBL=7
  2298.  
  2299. TABOBL(2)='REC'
  2300. TABOBL(3)='FK'
  2301. TABOBL(4)='FA'
  2302. TABOBL(5)='FB'
  2303. TABOBL(6)='FC'
  2304. TABOBL(7)='FD'
  2305. ELSE
  2306. * si FROTTEMENT3 ou FROTTEMENT4
  2307. CALL PLACE(MATMOD,NMAT,IFT1,'FROTTEMENT3')
  2308. CALL PLACE(MATMOD,NMAT,IFT2,'FROTTEMENT4')
  2309. IF(IFT1.NE.0.OR.IFT2.NE.0) THEN
  2310. JGOBL=2
  2311.  
  2312. TABOBL(2)='FK'
  2313. ENDIF
  2314. ENDIF
  2315. GOTO 9999
  2316. ENDIF
  2317.  
  2318. *--------------------------------------------------------------------
  2319. * CAS DE LA FORMULATION NAVIER_STOKES
  2320. *--------------------------------------------------------------------
  2321. CALL PLACE(FORMOD,NFOR,INAST,'NAVIER_STOKES')
  2322. IF (INAST.GT.0) THEN
  2323. iplac = imatee
  2324. if (iplac.eq.0) then
  2325. iret = 0
  2326. goto 9999
  2327. endif
  2328. *
  2329. JGOBL = 0
  2330. JGFAC = 0
  2331.  
  2332. if (iplac.eq.4) then
  2333. *NLIN
  2334. JGOBL0 = JGOBL
  2335. JGOBL = JGOBL0+1
  2336. TABOBL(JGOBL0 + 1) ='REYN'
  2337. JGOBL0 = JGFAC
  2338. JGFAC = JGOBL0+1
  2339. TABFAC(JGOBL0 + 1) ='FREQ'
  2340. endif
  2341. GOTO 9999
  2342. ENDIF
  2343.  
  2344. *--------------------------------------------------------------------
  2345. * CAS DE LA FORMULATION MELANGE
  2346. *--------------------------------------------------------------------
  2347. CALL PLACE(FORMOD,NFOR,ICONV,'MELANGE')
  2348. IF (ICONV.NE.0) THEN
  2349. CALL MODMEL(MOMODL,NMOD)
  2350. CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(1))
  2351. if (iplac.eq.0) then
  2352. goto 9999
  2353. endif
  2354. INMAT=1
  2355. *
  2356. * a priori elements massifs ou coques : pas de verif
  2357. *
  2358. IF (IPLAC.NE.3.AND.IPLAC.NE.4) THEN
  2359. JGOBL = 0
  2360. JGFAC = 0
  2361. CALL IDMETA(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  2362. IF (IRET.EQ.0) GOTO 9999
  2363. ELSE
  2364. ENDIF
  2365. GOTO 9999
  2366. ENDIF
  2367.  
  2368. *--------------------------------------------------------------------
  2369. * CAS DE LA FORMULATION LIAISON
  2370. *--------------------------------------------------------------------
  2371. CALL PLACE(FORMOD,NFOR,ICONV,'LIAISON')
  2372. IF (ICONV.NE.0) THEN
  2373. iplac = imatee
  2374. if (iplac.eq.0) then
  2375. iret = 0
  2376. goto 9999
  2377. endif
  2378. *
  2379. JGOBL = 0
  2380. JGFAC = 0
  2381.  
  2382. * 'SORT' facultatif dans tous les cas
  2383. IF (iplac.EQ.1) THEN
  2384. *PO_PL_FL
  2385. JGOBL0 = JGOBL
  2386. JGOBL = JGOBL0+7
  2387. TABOBL(JGOBL0 + 1) ='NORM'
  2388. TABOBL(JGOBL0 + 2) ='INER'
  2389. TABOBL(JGOBL0 + 3) ='CONV'
  2390. TABOBL(JGOBL0 + 4) ='VISC'
  2391. TABOBL(JGOBL0 + 5) ='PELO'
  2392. TABOBL(JGOBL0 + 6) ='PRAP'
  2393. TABOBL(JGOBL0 + 7) ='JFLU'
  2394.  
  2395. ELSEIF (iplac.eq.2) then
  2396. *PO_PL_FR
  2397. JGOBL0 = JGOBL
  2398. JGOBL = JGOBL0+7
  2399. TABOBL(JGOBL0 + 1) ='NORM'
  2400. TABOBL(JGOBL0 + 2) ='RAID'
  2401. TABOBL(JGOBL0 + 3) ='JEU'
  2402. TABOBL(JGOBL0 + 4) ='GLIS'
  2403. TABOBL(JGOBL0 + 5) ='ADHE'
  2404. TABOBL(JGOBL0 + 6) ='RTAN'
  2405. TABOBL(JGOBL0 + 7) ='ATAN'
  2406.  
  2407. JGOBL0 = JGFAC
  2408. JGFAC = JGOBL0 +2
  2409. TABFAC(JGOBL0+1) = 'AMOR'
  2410. TABFAC(JGOBL0+2) = 'LOIC'
  2411.  
  2412. ELSEIF (iplac.eq.3) then
  2413. *PO_PL
  2414. JGOBL0 = JGOBL
  2415. JGOBL= JGOBL0+3
  2416. TABOBL(JGOBL0 + 1) ='NORM'
  2417. TABOBL(JGOBL0 + 2) ='RAID'
  2418. TABOBL(JGOBL0 + 3) ='JEU'
  2419.  
  2420. JGOBL0 = JGFAC
  2421. JGFAC = JGOBL0+4
  2422. TABFAC(JGOBL0+1) = 'LOIC'
  2423. TABFAC(JGOBL0+2) = 'PERM'
  2424. TABFAC(JGOBL0 + 3) ='SPLA'
  2425. TABFAC(JGOBL0 + 4) ='AMOR'
  2426.  
  2427. ELSEIF (iplac.eq.4) then
  2428. *PO_PO_FR
  2429. JGOBL0 = JGOBL
  2430. JGOBL= JGOBL0 + 8
  2431. TABOBL(JGOBL0 + 1) ='NORM'
  2432. TABOBL(JGOBL0 + 2) ='RAID'
  2433. TABOBL(JGOBL0 + 3) ='JEU'
  2434. TABOBL(JGOBL0 + 4) ='POIB'
  2435. TABOBL(JGOBL0 + 5) ='ADHE'
  2436. TABOBL(JGOBL0 + 6) ='RTAN'
  2437. TABOBL(JGOBL0 + 7) ='ATAN'
  2438. TABOBL(JGOBL0 + 8) ='GLIS'
  2439.  
  2440. JGOBL0 = JGFAC
  2441. JGFAC = JGOBL0 + 3
  2442. TABFAC(JGOBL0+1) = 'AMOR'
  2443. TABFAC(JGOBL0+2) = 'LOIC'
  2444. TABFAC(JGOBL0+3) = 'MODE'
  2445.  
  2446. ELSEIF (iplac.eq.5) then
  2447. *PO_PO_DP
  2448. JGOBL0 = JGOBL
  2449. JGOBL= JGOBL0 + 6
  2450. TABOBL(JGOBL0 + 1) ='NORM'
  2451. TABOBL(JGOBL0 + 2) ='ECRO'
  2452. TABOBL(JGOBL0 + 3) ='JEU'
  2453. TABOBL(JGOBL0 + 4) ='POIB'
  2454. TABOBL(JGOBL0 + 5) ='PERM'
  2455. TABOBL(JGOBL0 + 6) ='LOIC'
  2456.  
  2457. JGOBL0 = JGFAC
  2458. JGFAC = JGOBL0 + 1
  2459. TABFAC(JGOBL0+1) = 'AMOR'
  2460.  
  2461. ELSEIF (iplac.eq.6) then
  2462. *PO_PO_RP
  2463. JGOBL0 = JGOBL
  2464. JGOBL= JGOBL0 + 6
  2465. TABOBL(JGOBL0 + 1) ='AXRO'
  2466. TABOBL(JGOBL0 + 2) ='ECRO'
  2467. TABOBL(JGOBL0 + 3) ='JEU'
  2468. TABOBL(JGOBL0 + 4) ='POIB'
  2469. TABOBL(JGOBL0 + 5) ='PERM'
  2470. TABOBL(JGOBL0 + 6) ='LOIC'
  2471.  
  2472. JGOBL0 = JGFAC
  2473. JGFAC = JGOBL0 + 2
  2474. TABFAC(JGOBL0+1) = 'AMOR'
  2475. TABFAC(JGOBL0+2) = 'ELAS'
  2476.  
  2477. ELSEIF (iplac.eq.7) then
  2478. *PO_PO
  2479. JGOBL0 = JGOBL
  2480. JGOBL= JGOBL0 + 5
  2481. TABOBL(JGOBL0 + 1) ='NORM'
  2482. TABOBL(JGOBL0 + 2) ='RAID'
  2483. TABOBL(JGOBL0 + 3) ='JEU'
  2484. TABOBL(JGOBL0 + 4) ='POIB'
  2485. TABOBL(JGOBL0 + 5) ='PERM'
  2486.  
  2487. JGOBL0 = JGFAC
  2488. JGFAC = JGOBL0 + 2
  2489. TABFAC(JGOBL0+1) = 'AMOR'
  2490. TABFAC(JGOBL0+2) = 'LOIC'
  2491.  
  2492. ELSEIF (iplac.eq.8) then
  2493. *PO_CE_MO
  2494. JGOBL0 = JGOBL
  2495. JGOBL= JGOBL0 + 8
  2496. TABOBL(JGOBL0 + 1) ='NORM'
  2497. TABOBL(JGOBL0 + 2) ='RAID'
  2498. TABOBL(JGOBL0 + 3) ='PCER'
  2499. TABOBL(JGOBL0 + 4) ='RAYO'
  2500. TABOBL(JGOBL0 + 5) ='GLIS'
  2501. TABOBL(JGOBL0 + 6) ='ADHE'
  2502. TABOBL(JGOBL0 + 7) ='RTAN'
  2503. TABOBL(JGOBL0 + 8) ='ATAN'
  2504.  
  2505. JGOBL0 = JGFAC
  2506. JGFAC= JGOBL0 + 2
  2507. TABFAC(JGOBL0+1) = 'CINT'
  2508. TABFAC(JGOBL0+2) = 'AMOR'
  2509.  
  2510. ELSEIF (iplac.eq.9) then
  2511. *PO_CE_FR
  2512. JGOBL0 = JGOBL
  2513. JGOBL= JGOBL0 + 8
  2514. TABOBL(JGOBL0 + 1) ='NORM'
  2515. TABOBL(JGOBL0 + 2) ='RAID'
  2516. TABOBL(JGOBL0 + 3) ='EXCE'
  2517. TABOBL(JGOBL0 + 4) ='RAYO'
  2518. TABOBL(JGOBL0 + 5) ='GLIS'
  2519. TABOBL(JGOBL0 + 6) ='ADHE'
  2520. TABOBL(JGOBL0 + 7) ='RTAN'
  2521. TABOBL(JGOBL0 + 8) ='ATAN'
  2522.  
  2523. JGOBL0 = JGFAC
  2524. JGFAC = JGOBL0 + 2
  2525. TABFAC(JGOBL0+1) = 'CINT'
  2526. TABFAC(JGOBL0+2) = 'AMOR'
  2527.  
  2528. ELSEIF (iplac.eq.10) then
  2529. *PO_CE
  2530. JGOBL0 = JGOBL
  2531. JGOBL= JGOBL0 + 4
  2532. TABOBL(JGOBL0 + 1) ='NORM'
  2533. TABOBL(JGOBL0 + 2) ='RAID'
  2534. TABOBL(JGOBL0 + 3) ='EXCE'
  2535. TABOBL(JGOBL0 + 4) ='RAYO'
  2536.  
  2537. JGOBL0 = JGFAC
  2538. JGFAC = JGOBL0 + 1
  2539. TABFAC(JGOBL0+1) = 'AMOR'
  2540.  
  2541. ELSEIF (iplac.eq.11) then
  2542. *CE_PL_FR
  2543. JGOBL0 = JGOBL
  2544. JGOBL= JGOBL0 + 8
  2545. TABOBL(JGOBL0 + 1) ='NORM'
  2546. TABOBL(JGOBL0 + 2) ='RAID'
  2547. TABOBL(JGOBL0 + 3) ='JEU'
  2548. TABOBL(JGOBL0 + 4) ='RAYS'
  2549. TABOBL(JGOBL0 + 5) ='GLIS'
  2550. TABOBL(JGOBL0 + 6) ='ADHE'
  2551. TABOBL(JGOBL0 + 7) ='RTAN'
  2552. TABOBL(JGOBL0 + 8) ='ATAN'
  2553.  
  2554. JGOBL0 = JGFAC
  2555. JGFAC = JGOBL0+1
  2556. TABFAC(JGOBL0+1) = 'AMOR'
  2557.  
  2558. ELSEIF (iplac.eq.12) then
  2559. *CE_CE_FR
  2560. JGOBL0 = JGOBL
  2561. JGOBL= JGOBL0 + 9
  2562. TABOBL(JGOBL0 + 1) ='NORM'
  2563. TABOBL(JGOBL0 + 2) ='RAID'
  2564. TABOBL(JGOBL0 + 3) ='EXCE'
  2565. TABOBL(JGOBL0 + 4) ='RAYS'
  2566. TABOBL(JGOBL0 + 5) ='GLIS'
  2567. TABOBL(JGOBL0 + 6) ='ADHE'
  2568. TABOBL(JGOBL0 + 7) ='RTAN'
  2569. TABOBL(JGOBL0 + 8) ='ATAN'
  2570. TABOBL(JGOBL0 + 9) ='RAYB'
  2571.  
  2572. JGOBL0 = JGFAC
  2573. JGFAC = JGOBL0 + 2
  2574. TABFAC(JGOBL0+1) = 'AMOR'
  2575. TABFAC(JGOBL0+2) = 'CINT'
  2576.  
  2577. ELSEIF (iplac.eq.13.or.iplac.eq.14) then
  2578. *PR_PR_IN ou PR_PR_EX
  2579. JGOBL0 = JGOBL
  2580. JGOBL= JGOBL0 + 5
  2581. TABOBL(JGOBL0 + 1) ='NORM'
  2582. TABOBL(JGOBL0 + 2) ='RAID'
  2583. TABOBL(JGOBL0 + 3) ='PFIX'
  2584. TABOBL(JGOBL0 + 4) ='PMOB'
  2585. TABOBL(JGOBL0 + 5) ='ERAI'
  2586.  
  2587. ELSEIF (iplac.eq.15) then
  2588. *LI_LI_FR
  2589. JGOBL0 = JGOBL
  2590. JGOBL= JGOBL0 + 8
  2591. TABOBL(JGOBL0 + 1) ='NORM'
  2592. TABOBL(JGOBL0 + 2) ='LIMA'
  2593. TABOBL(JGOBL0 + 3) ='LIES'
  2594. TABOBL(JGOBL0 + 4) ='RAID'
  2595. TABOBL(JGOBL0 + 5) ='GLIS'
  2596. TABOBL(JGOBL0 + 6) ='ADHE'
  2597. TABOBL(JGOBL0 + 7) ='RTAN'
  2598. TABOBL(JGOBL0 + 8) ='ATAN'
  2599.  
  2600. JGOBL0 = JGFAC
  2601. JGFAC = JGOBL0 + 4
  2602. TABFAC(JGOBL0+1) = 'JEU'
  2603. TABFAC(JGOBL0+2) = 'AMOR'
  2604. TABFAC(JGOBL0+3) = 'RECH'
  2605. TABFAC(JGOBL0+4) = 'SYME'
  2606.  
  2607. ELSEIF (iplac.eq.16) then
  2608. *LI_CE_FR
  2609. JGOBL0 = JGOBL
  2610. JGOBL= JGOBL0 + 8
  2611. TABOBL(JGOBL0 + 1) ='NORM'
  2612. TABOBL(JGOBL0 + 2) ='LIMA'
  2613. TABOBL(JGOBL0 + 3) ='LIES'
  2614. TABOBL(JGOBL0 + 4) ='RAID'
  2615. TABOBL(JGOBL0 + 5) ='GLIS'
  2616. TABOBL(JGOBL0 + 6) ='ADHE'
  2617. TABOBL(JGOBL0 + 7) ='RTAN'
  2618. TABOBL(JGOBL0 + 8) ='ATAN'
  2619.  
  2620. JGOBL0 = JGFAC
  2621. JGFAC = JGOBL0 + 5
  2622. TABFAC(JGOBL0+1) = 'AMOR'
  2623. TABFAC(JGOBL0+2) = 'RECH'
  2624. TABFAC(JGOBL0+3) = 'RAYO'
  2625. TABFAC(JGOBL0+4) = 'ACTN'
  2626. TABFAC(JGOBL0+5) = 'INVE'
  2627.  
  2628. ELSEIF (iplac.eq.17) then
  2629. *PA_FL_RO
  2630. JGOBL0 = JGOBL
  2631. JGOBL= JGOBL0+10
  2632. TABOBL(JGOBL0 + 1) ='LONG'
  2633. TABOBL(JGOBL0 + 2) ='RAYO'
  2634. TABOBL(JGOBL0 + 3) ='VISC'
  2635. TABOBL(JGOBL0 + 4) ='RHOF'
  2636. TABOBL(JGOBL0 + 5) ='PADM'
  2637. TABOBL(JGOBL0 + 6) ='VROT'
  2638. TABOBL(JGOBL0 + 7) ='EPSI'
  2639. TABOBL(JGOBL0 + 8) ='PHII'
  2640. TABOBL(JGOBL0 + 9) ='AFFI'
  2641. TABOBL(JGOBL0 + 10)='TLOB'
  2642.  
  2643. JGOBL0 = JGFAC
  2644. JGFAC = JGOBL0 + 1
  2645. TABFAC(JGOBL0+1) = 'AMOR'
  2646.  
  2647. ELSEIF (iplac.eq.23) then
  2648. *NEWMARK MODAL
  2649. JGOBL0 = JGOBL
  2650. JGOBL= JGOBL0+3
  2651. TABOBL(JGOBL0 + 1) ='JEU'
  2652. TABOBL(JGOBL0 + 2) ='MASS'
  2653. TABOBL(JGOBL0 + 3) ='FREQ'
  2654.  
  2655. JGOBL0 = JGFAC
  2656. JGFAC = JGOBL0 + 3
  2657. TABFAC(JGOBL0+1) = 'EXCE'
  2658. TABFAC(JGOBL0+2) = 'FROT'
  2659. TABFAC(JGOBL0+3) = 'MOFR'
  2660. ENDIF
  2661.  
  2662. * 'SORT' facultatif dans tous les cas
  2663. JGOBL0 = JGFAC
  2664. JGFAC= JGOBL0+1
  2665. TABFAC(JGOBL0+1) = 'SORT'
  2666. GOTO 9999
  2667. ENDIF
  2668.  
  2669. *--------------------------------------------------------------------
  2670. * CAS DE LA FORMULATION ELECTROSTATIQUE
  2671. *--------------------------------------------------------------------
  2672. CALL PLACE(FORMOD,NFOR,IELEC,'ELECTROSTATIQUE')
  2673. IF (IELEC.NE.0) THEN
  2674. C -- Permittivite isotrope
  2675. IF (IMATEE.EQ.1) THEN
  2676. C* IF (CMATEE.EQ.'ISOTROPE') THEN
  2677. C* IF (MATMOD(1).EQ.'ISOTROPE ') THEN
  2678. JGOBL = 1
  2679. TABOBL(1)='PEL '
  2680. C -- Permittivite orthotrope
  2681. ELSEIF (IMATEE.EQ.2) THEN
  2682. C* ELSEIF (CMATEE.EQ.'ORTHOTRO') THEN
  2683. C* ELSEIF (MATMOD(1).EQ.'ORTHOTROPE ') THEN
  2684. C ---- Elements massifs bidimensionnels PLAN et AXISYMETRIQUE
  2685. IF (IDIM.EQ.2) THEN
  2686. IF (IFOMOD.NE.1) THEN
  2687. JGOBL = 4
  2688. TABOBL(1) = 'PE1 '
  2689. TABOBL(2) = 'PE2 '
  2690. TABOBL(3) = 'V1X '
  2691. TABOBL(4) = 'V1Y '
  2692.  
  2693. C ---- Elements massifs bidimensionnels FOURIER
  2694. ELSE
  2695. JGOBL = 5
  2696. TABOBL(1) = 'PE1 '
  2697. TABOBL(2) = 'PE2 '
  2698. TABOBL(3) = 'PE3 '
  2699. TABOBL(4) = 'V1X '
  2700. TABOBL(5) = 'V1Y '
  2701. ENDIF
  2702. C ---- Elements massifs TRIDimensionnels
  2703. ELSEIF (IDIM.EQ.3) THEN
  2704. JGOBL = 9
  2705. TABOBL(1) = 'PE1 '
  2706. TABOBL(2) = 'PE2 '
  2707. TABOBL(3) = 'PE3 '
  2708. TABOBL(4) = 'V1X '
  2709. TABOBL(5) = 'V1Y '
  2710. TABOBL(6) = 'V1Z '
  2711. TABOBL(7) = 'V2X '
  2712. TABOBL(8) = 'V2Y '
  2713. TABOBL(9) = 'V2Z '
  2714. ENDIF
  2715.  
  2716. C -- Permittivite anisotrope
  2717. ELSEIF (IMATEE.EQ.3) THEN
  2718. C* ELSEIF (CMATEE.EQ.'ANISOTRO') THEN
  2719. C* ELSEIF (MATMOD(1).EQ.'ANISOTROPE ') THEN
  2720. C ---- Elements massifs bidimensionnels PLAN et AXISYMETRIQUE
  2721. IF (IDIM.EQ.2) THEN
  2722. IF (IFOMOD.NE.1) THEN
  2723. JGOBL = 5
  2724. TABOBL(1) = 'PE11 '
  2725. TABOBL(2) = 'PE22 '
  2726. TABOBL(3) = 'PE21 '
  2727. TABOBL(4) = 'V1X '
  2728. TABOBL(5) = 'V1Y '
  2729.  
  2730. C ---- Elements massifs bidimensionnels FOURIER
  2731. ELSE
  2732. JGOBL = 6
  2733. TABOBL(1) = 'PE11 '
  2734. TABOBL(2) = 'PE22 '
  2735. TABOBL(3) = 'PE21 '
  2736. TABOBL(4) = 'PE33 '
  2737. TABOBL(5) = 'V1X '
  2738. TABOBL(6) = 'V1Y '
  2739. ENDIF
  2740.  
  2741. C ---- Elements massifs TRIDimensionnels
  2742. ELSEIF (IDIM.EQ.3) THEN
  2743. JGOBL = 12
  2744. TABOBL( 1) = 'PE11 '
  2745. TABOBL( 2) = 'PE22 '
  2746. TABOBL( 3) = 'PE33 '
  2747. TABOBL( 4) = 'PE21 '
  2748. TABOBL( 5) = 'PE31 '
  2749. TABOBL( 6) = 'PE32 '
  2750. TABOBL( 7) = 'V1X '
  2751. TABOBL( 8) = 'V1Y '
  2752. TABOBL( 9) = 'V1Z '
  2753. TABOBL(10) = 'V2X '
  2754. TABOBL(11) = 'V2Y '
  2755. TABOBL(12) = 'V2Z '
  2756. ENDIF
  2757.  
  2758. ELSE
  2759. IRET = 0
  2760. CALL ERREUR(5)
  2761. ENDIF
  2762. GOTO 9999
  2763. ENDIF
  2764.  
  2765. *-----------------------------------------------------------------------
  2766. * CAS DE LA FORMULATION 'DIFFUSION'
  2767. *-----------------------------------------------------------------------
  2768. CALL PLACE(FORMOD,NFOR,IDIFF,'DIFFUSION')
  2769. IF (IDIFF.NE.0) THEN
  2770. JGFAC = 1
  2771. TABFAC(1)='CDIF'
  2772.  
  2773. C- Diffusion isotrope
  2774. IF (IMATEE .EQ. 1) THEN
  2775. JGOBL = 1
  2776. C*8 TABOBL(1) = 'KDIFF '
  2777. TABOBL(1) = 'KD '
  2778.  
  2779. C- Diffusion orthotrope
  2780. ELSEIF (IMATEE .EQ. 2) THEN
  2781. C --- Elements massifs bidimensionnels PLAN et AXISYMETRIQUE
  2782. IF (IDIM.EQ.2) THEN
  2783. IF (IFOMOD.NE.1) THEN
  2784. JGOBL = 4
  2785. TABOBL(1) = 'KD1 '
  2786. TABOBL(2) = 'KD2 '
  2787. TABOBL(3) = 'V1X '
  2788. TABOBL(4) = 'V1Y '
  2789.  
  2790. C --- Elements massifs bidimensionnels FOURIER
  2791. ELSE
  2792. JGOBL = 5
  2793. TABOBL(1) = 'KD1 '
  2794. TABOBL(2) = 'KD2 '
  2795. TABOBL(3) = 'KD3 '
  2796. TABOBL(4) = 'V1X '
  2797. TABOBL(5) = 'V1Y '
  2798. ENDIF
  2799.  
  2800. C --- Elements massifs TRIDimensionnels
  2801. ELSEIF (IDIM.EQ.3) THEN
  2802. JGOBL = 9
  2803. TABOBL(1) = 'KD1 '
  2804. TABOBL(2) = 'KD2 '
  2805. TABOBL(3) = 'KD3 '
  2806. TABOBL(4) = 'V1X '
  2807. TABOBL(5) = 'V1Y '
  2808. TABOBL(6) = 'V1Z '
  2809. TABOBL(7) = 'V2X '
  2810. TABOBL(8) = 'V2Y '
  2811. TABOBL(9) = 'V2Z '
  2812. ENDIF
  2813.  
  2814. C - Diffusion anisotrope
  2815. ELSEIF (IMATEE.EQ.3) THEN
  2816. C --- Elements massifs bidimensionnels PLAN et AXISYMETRIQUE
  2817. IF (IDIM.EQ.2) THEN
  2818. IF (IFOMOD.NE.1) THEN
  2819. JGOBL = 5
  2820. TABOBL(1) = 'KD11 '
  2821. TABOBL(2) = 'KD22 '
  2822. TABOBL(3) = 'KD21 '
  2823. TABOBL(4) = 'V1X '
  2824. TABOBL(5) = 'V1Y '
  2825.  
  2826. C --- Elements massifs bidimensionnels FOURIER
  2827. ELSE
  2828. JGOBL = 6
  2829. TABOBL(1) = 'KD11 '
  2830. TABOBL(2) = 'KD22 '
  2831. TABOBL(3) = 'KD21 '
  2832. TABOBL(4) = 'KD33 '
  2833. TABOBL(5) = 'V1X '
  2834. TABOBL(6) = 'V1Y '
  2835. ENDIF
  2836.  
  2837. C --- Elements massifs TRIDimensionnels
  2838. ELSEIF (IDIM.EQ.3) THEN
  2839. JGOBL = 12
  2840. TABOBL( 1) = 'KD11 '
  2841. TABOBL( 2) = 'KD22 '
  2842. TABOBL( 3) = 'KD33 '
  2843. TABOBL( 4) = 'KD21 '
  2844. TABOBL( 5) = 'KD31 '
  2845. TABOBL( 6) = 'KD32 '
  2846. TABOBL( 7) = 'V1X '
  2847. TABOBL( 8) = 'V1Y '
  2848. TABOBL( 9) = 'V1Z '
  2849. TABOBL(10) = 'V2X '
  2850. TABOBL(11) = 'V2Y '
  2851. TABOBL(12) = 'V2Z '
  2852. ENDIF
  2853. ENDIF
  2854.  
  2855. JGOBL0 = JGOBL
  2856. JGFAC0 = JGFAC
  2857.  
  2858. * Loi FICK :
  2859. * ------------
  2860. IF (INATUU.EQ.0) THEN
  2861. * Rien a ajouter
  2862.  
  2863. * Loi SORET :
  2864. * -------------
  2865. ELSEIF (INATUU.EQ.1) THEN
  2866. C -- Diffusion isotrope
  2867. IF (IMATEE.EQ.1) THEN
  2868. JGOBL = JGOBL+1
  2869. C*8 TABOBL(JGOBL0+1) = 'KSORET '
  2870. TABOBL(JGOBL0+1) = 'KS '
  2871.  
  2872. C -- Diffusion orthotrope
  2873. ELSEIF (IMATEE.EQ.2) THEN
  2874. C ---- Elements massifs bidimensionnels PLAN et AXISYMETRIQUE
  2875. IF (IFOMOD.LE.0) THEN
  2876. JGOBL = JGOBL + 2
  2877. C*8 TABOBL(JGOBL0+1) = 'KSORET1 '
  2878. C*8 TABOBL(JGOBL0+2) = 'KSORET2 '
  2879. TABOBL(JGOBL0+1) = 'KS1 '
  2880. TABOBL(JGOBL0+2) = 'KS2 '
  2881.  
  2882. C ---- Elements massifs bidimensionnels FOURIER
  2883. C ---- Elements massifs TRIDimensionnels
  2884. ELSE
  2885. JGOBL = JGOBL + 3
  2886. C*8 TABOBL(JGOBL0+1) = 'KSORET1 '
  2887. C*8 TABOBL(JGOBL0+2) = 'KSORET2 '
  2888. C*8 TABOBL(JGOBL0+3) = 'KSORET3 '
  2889. TABOBL(JGOBL0+1) = 'KS1 '
  2890. TABOBL(JGOBL0+2) = 'KS2 '
  2891. TABOBL(JGOBL0+3) = 'KS3 '
  2892. ENDIF
  2893.  
  2894. C -- Diffusion anisotrope
  2895. ELSEIF (IMATEE.EQ.3) THEN
  2896. C ---- Elements massifs bidimensionnels PLAN et AXISYMETRIQUE
  2897. IF (IFOMOD.LE.0) THEN
  2898. JGOBL = JGOBL + 3
  2899. C*8 TABOBL(JGOBL0+1) = 'KSORET11'
  2900. C*8 TABOBL(JGOBL0+2) = 'KSORET22'
  2901. C*8 TABOBL(JGOBL0+3) = 'KSORET21'
  2902. TABOBL(JGOBL0+1) = 'KS11 '
  2903. TABOBL(JGOBL0+2) = 'KS22 '
  2904. TABOBL(JGOBL0+3) = 'KS21 '
  2905.  
  2906. C ---- Elements massifs bidimensionnels FOURIER
  2907. ELSEIF (IFOMOD.EQ.1) THEN
  2908. JGOBL = JGOBL + 4
  2909. C*8 TABOBL(JGOBL0+1) = 'KSORET11'
  2910. C*8 TABOBL(JGOBL0+2) = 'KSORET22'
  2911. C*8 TABOBL(JGOBL0+3) = 'KSORET21'
  2912. C*8 TABOBL(JGOBL0+4) = 'KSORET33'
  2913. TABOBL(JGOBL0+1) = 'KS11 '
  2914. TABOBL(JGOBL0+2) = 'KS22 '
  2915. TABOBL(JGOBL0+3) = 'KS21 '
  2916. TABOBL(JGOBL0+4) = 'KS33 '
  2917.  
  2918. C ---- Elements massifs TRIDimensionnels
  2919. ELSEIF (IFOMOD.EQ.2) THEN
  2920. JGOBL = JGOBL + 6
  2921. C*8 TABOBL(JGOBL0+1) = 'KSORET11'
  2922. C*8 TABOBL(JGOBL0+2) = 'KSORET22'
  2923. C*8 TABOBL(JGOBL0+3) = 'KSORET33'
  2924. C*8 TABOBL(JGOBL0+4) = 'KSORET21'
  2925. C*8 TABOBL(JGOBL0+5) = 'KSORET31'
  2926. C*8 TABOBL(JGOBL0+6) = 'KSORET32'
  2927. TABOBL(JGOBL0+1) = 'KS11 '
  2928. TABOBL(JGOBL0+2) = 'KS22 '
  2929. TABOBL(JGOBL0+3) = 'KS33 '
  2930. TABOBL(JGOBL0+4) = 'KS21 '
  2931. TABOBL(JGOBL0+5) = 'KS31 '
  2932. TABOBL(JGOBL0+6) = 'KS32 '
  2933. ENDIF
  2934.  
  2935. ELSE
  2936. C -- Erreur ne devant pas survenir
  2937. CALL ERREUR(261)
  2938. RETURN
  2939. ENDIF
  2940.  
  2941. * Loi UTILISATEUR :
  2942. * -------------------
  2943. * ELSEIF (INATUU.EQ.2) THEN
  2944. * Par defaut pas de composantes autres que celles de l'utilisateur
  2945.  
  2946. ELSEIF (INATUU.EQ.5) THEN
  2947. * diffusion avec dissipation affine
  2948. JGOBL = JGOBL+2
  2949. TABOBL(JGOBL0+1)='DAF0'
  2950. TABOBL(JGOBL0+2)='DAF1'
  2951.  
  2952. ELSEIF (INATUU.EQ.6) THEN
  2953. * diffusion avec dissipation visqueuse
  2954. JGOBL = JGOBL+2
  2955. TABOBL(JGOBL0+1)='DVQ0'
  2956. TABOBL(JGOBL0+2)