Télécharger idmatr.eso

Retour à la liste

Numérotation des lignes :

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