Télécharger idmatr.eso

Retour à la liste

Numérotation des lignes :

idmatr
  1. C IDMATR SOURCE JK148537 25/02/17 21:15:05 12156
  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. do jm = 1,matmod(/2)
  539. IF (matmod(jm).eq.'IMPEDANCE') THEN
  540. INMAT = 3
  541. goto 19
  542. ENDIF
  543. enddo
  544. endif
  545.  
  546. IF (IPLAC.NE.0) THEN
  547. INMAT=3
  548. IF (IPLAC.EQ.1) THEN
  549. INMAT=3
  550. GOTO 10
  551. ENDIF
  552. IF (IPLAC.EQ.2) THEN
  553. *
  554. * MATERIAU ELASTIQUE ORTHOTROPE
  555. *
  556. IF (MFR.EQ.75) THEN
  557. * JOINT UNIDIMENSIONNEL JOI1
  558. *
  559. IF(IFOUR.EQ.2)THEN
  560. JGOBL=12
  561. TABOBL(1)='V1X '
  562. TABOBL(2)='V1Y '
  563. TABOBL(3)='V1Z '
  564. TABOBL(4)='V2X '
  565. TABOBL(5)='V2Y '
  566. TABOBL(6)='V2Z '
  567. TABOBL(7)='KN '
  568. TABOBL(8)='KS1 '
  569. TABOBL(9)='KS2 '
  570. TABOBL(10)='QN '
  571. TABOBL(11)='QS1 '
  572. TABOBL(12)='QS2 '
  573. *
  574. JGFAC=10
  575. TABFAC(1)='MASS'
  576. TABFAC(2)='JX '
  577. TABFAC(3)='JY '
  578. TABFAC(4)='JZ '
  579. TABFAC(5)='ALPN'
  580. TABFAC(6)='ALP1'
  581. TABFAC(7)='ALP2'
  582. TABFAC(8)='ALQN'
  583. TABFAC(9)='ALQ1'
  584. TABFAC(10)='ALQ2'
  585. *
  586. ELSEIF(IFOUR.EQ.-2.OR.IFOUR.EQ.-1)THEN
  587. JGOBL=5
  588. TABOBL(1)='V1X '
  589. TABOBL(2)='V1Y '
  590. TABOBL(3)='KN '
  591. TABOBL(4)='KS '
  592. TABOBL(5)='QS '
  593. *
  594. JGFAC=6
  595. TABFAC(1)='MASS'
  596. TABFAC(2)='JZ'
  597. TABFAC(4)='ALPN'
  598. TABFAC(5)='ALPS'
  599. TABFAC(6)='ALQS'
  600. ENDIF
  601. *
  602. ELSEIF (MFR.EQ.3) THEN
  603. * COQUES MINCES
  604. *
  605. JGOBL=6
  606. TABOBL(1)='YG1 '
  607. TABOBL(2)='YG2 '
  608. TABOBL(3)='NU12'
  609. TABOBL(4)='G12 '
  610. TABOBL(5)='V1X '
  611. TABOBL(6)='V1Y '
  612. *
  613. IF(IFOUR.EQ.-2) THEN
  614. JGFAC=4
  615. ELSE
  616. JGFAC=3
  617. ENDIF
  618. TABFAC(1)='ALP1'
  619. TABFAC(2)='ALP2'
  620. TABFAC(3)='RHO '
  621. IF(IFOUR.EQ.-2) TABFAC(4)='DIM3'
  622.  
  623. ELSEIF (MFR.EQ.9.OR.MFR.EQ.5) THEN
  624. * COQUES AVEC CISAILLEMENT TRANSVERSE
  625. *
  626. JGOBL=8
  627. TABOBL(1)='YG1 '
  628. TABOBL(2)='YG2 '
  629. TABOBL(3)='NU12'
  630. TABOBL(4)='G12 '
  631. TABOBL(5)='G23 '
  632. TABOBL(6)='G13 '
  633. TABOBL(7)='V1X '
  634. TABOBL(8)='V1Y '
  635. *
  636. JGFAC=3
  637. TABFAC(1)='RHO '
  638. TABFAC(2)='ALP1'
  639. TABFAC(3)='ALP2'
  640.  
  641. ELSEIF (MFR.EQ.1.OR.MFR.EQ.31) THEN
  642. * ELEMENTS MASSIFS
  643. *
  644. IF(IDIM.EQ.3)THEN
  645. * ELEMENTS 3D
  646. JGOBL=15
  647. TABOBL(1)='YG1 '
  648. TABOBL(2)='YG2 '
  649. TABOBL(3)='YG3 '
  650. TABOBL(4)='NU12'
  651. TABOBL(5)='NU23'
  652. TABOBL(6)='NU13'
  653. TABOBL(7)='G12 '
  654. TABOBL(8)='G23 '
  655. TABOBL(9)='G13 '
  656. TABOBL(10)='V1X '
  657. TABOBL(11)='V1Y '
  658. TABOBL(12)='V1Z '
  659. TABOBL(13)='V2X '
  660. TABOBL(14)='V2Y '
  661. TABOBL(15)='V2Z '
  662. *
  663. JGFAC=4
  664. TABFAC(1)='RHO '
  665. TABFAC(2)='ALP1'
  666. TABFAC(3)='ALP2'
  667. TABFAC(4)='ALP3'
  668.  
  669. ELSEIF (IDIM.EQ.2) THEN
  670. IF(IFOUR.EQ.-2)THEN
  671. * CONTRAINTE PLANE
  672. JGOBL=9
  673. TABOBL(1)='YG1 '
  674. TABOBL(2)='YG2 '
  675. TABOBL(3)='NU12'
  676. TABOBL(4)='G12 '
  677. TABOBL(5)='V1X '
  678. TABOBL(6)='V1Y '
  679. TABOBL(7)='YG3 '
  680. TABOBL(8)='NU23'
  681. TABOBL(9)='NU13'
  682. *
  683. JGFAC=4
  684. TABFAC(1)='RHO '
  685. TABFAC(2)='ALP1'
  686. TABFAC(3)='ALP2'
  687. TABFAC(4)='DIM3'
  688.  
  689. ELSEIF (IFOUR.EQ.-1.OR.IFOUR.EQ.0.OR.IFOUR.EQ.-3) THEN
  690. * DEFORMATION PLANE ,AXISYMETRIE
  691. JGOBL=9
  692. TABOBL(1)='YG1 '
  693. TABOBL(2)='YG2 '
  694. TABOBL(3)='YG3 '
  695. TABOBL(4)='NU12'
  696. TABOBL(5)='NU23'
  697. TABOBL(6)='NU13'
  698. TABOBL(7)='G12 '
  699. TABOBL(8)='V1X '
  700. TABOBL(9)='V1Y '
  701. *
  702. JGFAC=4
  703. TABFAC(1)='RHO '
  704. TABFAC(2)='ALP1'
  705. TABFAC(3)='ALP2'
  706. TABFAC(4)='ALP3'
  707.  
  708. ELSEIF (IFOUR.EQ.1) THEN
  709. * AXISYMETRIE DE FOURIER
  710. JGOBL=15
  711. TABOBL(1)='YG1 '
  712. TABOBL(2)='YG2 '
  713. TABOBL(3)='YG3 '
  714. TABOBL(4)='NU12'
  715. TABOBL(5)='NU23'
  716. TABOBL(6)='NU13'
  717. TABOBL(7)='G12 '
  718. TABOBL(8)='G23 '
  719. TABOBL(9)='G13 '
  720. TABOBL(10)='V1X '
  721. TABOBL(11)='V1Y '
  722. TABOBL(12)='V1Z '
  723. TABOBL(13)='V2X '
  724. TABOBL(14)='V2Y '
  725. TABOBL(15)='V2Z '
  726. *
  727. JGFAC=4
  728. TABFAC(1)='RHO '
  729. TABFAC(2)='ALP1'
  730. TABFAC(3)='ALP2'
  731. TABFAC(4)='ALP3'
  732. ENDIF
  733. ELSEIF (IDIM.EQ.1) THEN
  734. C= Dans le cas UNID SPHErique, on doit avoir YG2=YG3 et NU12=NU13
  735. C= et dans le cas thermomecanique ALP2=ALP3.
  736. JGOBL=6
  737. TABOBL(1)='YG1 '
  738. TABOBL(2)='YG2 '
  739. TABOBL(3)='YG3 '
  740. TABOBL(4)='NU12'
  741. TABOBL(5)='NU23'
  742. TABOBL(6)='NU13'
  743.  
  744. IF (IFOUR.EQ.6) THEN
  745. JGFAC=2
  746.  
  747. ELSEIF (IFOUR.EQ.5 .OR. IFOUR.EQ.10) THEN
  748. JGFAC=3
  749. TABFAC(3)='ALP3'
  750. ELSEIF (IFOUR.EQ.4.OR.IFOUR.EQ.8.OR.
  751. & IFOUR.EQ.13) THEN
  752. JGFAC=3
  753. TABFAC(3)='ALP2'
  754. ELSE
  755. JGFAC=4
  756. TABFAC(3)='ALP2'
  757. TABFAC(4)='ALP3'
  758. ENDIF
  759. TABFAC(1)='RHO '
  760. TABFAC(2)='ALP1'
  761. ENDIF
  762.  
  763. ELSEIF (MFR.EQ.35) THEN
  764. * ELEMENTS JOINTS
  765. IF (IFOUR.EQ.2) THEN
  766. JGOBL=5
  767. TABOBL(1)='KS1 '
  768. TABOBL(2)='KS2 '
  769. TABOBL(3)='KN '
  770. TABOBL(4)='V1X '
  771. TABOBL(5)='V1Y '
  772. *
  773. JGFAC=2
  774. TABFAC(1)='RHO '
  775. TABFAC(2)='ALPN'
  776. ENDIF
  777. ENDIF
  778.  
  779. ELSEIF (IPLAC.EQ.3)THEN
  780. * MATERIAU ANISOTROPE ELASTIQUE
  781. *
  782. IF(MFR.EQ.75)THEN
  783. * JOINT UNIDIMESIONNEL JOI1
  784. *
  785. IF (IFOUR.EQ.2) THEN
  786. JGOBL=27
  787. TABOBL(1)='V1X '
  788. TABOBL(2)='V1Y '
  789. TABOBL(3)='V1Z '
  790. TABOBL(4)='V2X '
  791. TABOBL(5)='V2Y '
  792. TABOBL(6)='V2Z '
  793. TABOBL(7)='D11 '
  794. TABOBL(8)='D22 '
  795. TABOBL(9)='D33 '
  796. TABOBL(10)='D44 '
  797. TABOBL(11)='D55 '
  798. TABOBL(12)='D66 '
  799. TABOBL(13)='D21 '
  800. TABOBL(14)='D31 '
  801. TABOBL(15)='D32 '
  802. TABOBL(16)='D41 '
  803. TABOBL(17)='D42 '
  804. TABOBL(18)='D43 '
  805. TABOBL(19)='D51 '
  806. TABOBL(20)='D52 '
  807. TABOBL(21)='D53 '
  808. TABOBL(22)='D54 '
  809. TABOBL(23)='D61 '
  810. TABOBL(24)='D62 '
  811. TABOBL(25)='D63 '
  812. TABOBL(26)='D64 '
  813. TABOBL(27)='D65 '
  814. *
  815. JGFAC=10
  816. * MASS: masse totale de l'élément joint
  817. TABFAC(1)='MASS'
  818. TABFAC(2)='JX '
  819. TABFAC(3)='JY '
  820. TABFAC(4)='JZ '
  821. TABFAC(5)='ALP1'
  822. TABFAC(6)='ALP2'
  823. TABFAC(7)='ALP3'
  824. TABFAC(8)='ALQ1'
  825. TABFAC(9)='ALQ2'
  826. TABFAC(10)='ALQ3'
  827. *
  828. ELSEIF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1) THEN
  829. JGOBL=8
  830. TABOBL(1)='V1X '
  831. TABOBL(2)='V1Y '
  832. TABOBL(3)='D11 '
  833. TABOBL(4)='D22 '
  834. TABOBL(5)='D33 '
  835. TABOBL(6)='D21 '
  836. TABOBL(7)='D31 '
  837. TABOBL(8)='D32 '
  838. *
  839. JGFAC=5
  840. * MASS: masse totale de l'élément joint
  841. TABFAC(1)='MASS'
  842. TABFAC(2)='JZ '
  843. TABFAC(3)='ALP1'
  844. TABFAC(4)='ALP2'
  845. TABFAC(5)='ALQ3'
  846. ENDIF
  847. ENDIF
  848. *
  849. * ELEMENTS MASSIFS
  850. *
  851. IF(MFR.EQ.1.OR.MFR.EQ.31)THEN
  852. IF(IDIM.EQ.3)THEN
  853. * ELEMENTS 3D
  854. JGOBL=27
  855. TABOBL(1)='D11 '
  856. TABOBL(2)='D21 '
  857. TABOBL(3)='D22 '
  858. TABOBL(4)='D31 '
  859. TABOBL(5)='D32 '
  860. TABOBL(6)='D33 '
  861. TABOBL(7)='D41 '
  862. TABOBL(8)='D42 '
  863. TABOBL(9)='D43 '
  864. TABOBL(10)='D44 '
  865. TABOBL(11)='D51 '
  866. TABOBL(12)='D52 '
  867. TABOBL(13)='D53 '
  868. TABOBL(14)='D54 '
  869. TABOBL(15)='D55 '
  870. TABOBL(16)='D61 '
  871. TABOBL(17)='D62 '
  872. TABOBL(18)='D63 '
  873. TABOBL(19)='D64 '
  874. TABOBL(20)='D65 '
  875. TABOBL(21)='D66 '
  876. TABOBL(22)='V1X '
  877. TABOBL(23)='V1Y '
  878. TABOBL(24)='V1Z '
  879. TABOBL(25)='V2X '
  880. TABOBL(26)='V2Y '
  881. TABOBL(27)='V2Z '
  882. *
  883. JGFAC=7
  884. TABFAC(1)='RHO '
  885. TABFAC(2)='ALP1'
  886. TABFAC(3)='ALP2'
  887. TABFAC(4)='ALP3'
  888. TABFAC(5)='AL12'
  889. TABFAC(6)='AL13'
  890. TABFAC(7)='AL23'
  891.  
  892. ELSEIF (IDIM.EQ.2) THEN
  893. IF (IFOUR.EQ.-2) THEN
  894. * CONTRAINTE PLANE
  895. JGOBL=12
  896. TABOBL(1)='D11 '
  897. TABOBL(2)='D21 '
  898. TABOBL(3)='D22 '
  899. TABOBL(4)='D41 '
  900. TABOBL(5)='D42 '
  901. TABOBL(6)='D44 '
  902. TABOBL(7)='V1X '
  903. TABOBL(8)='V1Y '
  904. TABOBL(9)='D31 '
  905. TABOBL(10)='D32 '
  906. TABOBL(11)='D33 '
  907. TABOBL(12)='D43 '
  908. *
  909. JGFAC=5
  910. TABFAC(1)='RHO '
  911. TABFAC(2)='ALP1'
  912. TABFAC(3)='ALP2'
  913. TABFAC(4)='AL12'
  914. TABFAC(5)='DIM3'
  915.  
  916. ELSEIF (IFOUR.EQ.-1.OR.IFOUR.EQ.0.OR.IFOUR.EQ.-3) THEN
  917. * DEFORMATION PLANE ,AXISYMETRIE
  918. JGOBL=12
  919. TABOBL(1)='D11 '
  920. TABOBL(2)='D21 '
  921. TABOBL(3)='D22 '
  922. TABOBL(4)='D31 '
  923. TABOBL(5)='D32 '
  924. TABOBL(6)='D33 '
  925. TABOBL(7)='D41 '
  926. TABOBL(8)='D42 '
  927. TABOBL(9)='D43 '
  928. TABOBL(10)='D44 '
  929. TABOBL(11)='V1X '
  930. TABOBL(12)='V1Y '
  931. *
  932. JGFAC=5
  933. TABFAC(1)='RHO '
  934. TABFAC(2)='ALP1'
  935. TABFAC(3)='ALP2'
  936. TABFAC(4)='AL12'
  937. TABFAC(5)='ALP3'
  938.  
  939. ELSEIF (IFOUR.EQ.1) THEN
  940. * AXISYMETRIE DE FOURIER
  941. JGOBL=15
  942. TABOBL(1)='D11 '
  943. TABOBL(2)='D21 '
  944. TABOBL(3)='D22 '
  945. TABOBL(4)='D31 '
  946. TABOBL(5)='D32 '
  947. TABOBL(6)='D33 '
  948. TABOBL(7)='D41 '
  949. TABOBL(8)='D42 '
  950. TABOBL(9)='D43 '
  951. TABOBL(10)='D44 '
  952. TABOBL(11)='D55 '
  953. TABOBL(12)='D65 '
  954. TABOBL(13)='D66 '
  955. TABOBL(14)='V1X '
  956. TABOBL(15)='V1Y '
  957. *
  958. JGFAC=5
  959. TABFAC(1)='RHO '
  960. TABFAC(2)='ALP1'
  961. TABFAC(3)='ALP2'
  962. TABFAC(4)='AL12'
  963. TABFAC(5)='ALP3'
  964. ENDIF
  965. ENDIF
  966. ENDIF
  967.  
  968. ELSEIF (IPLAC.EQ.4) THEN
  969. IF (MFR.EQ.33) THEN
  970. * MILIEU POREUX ISOTROPE
  971. *
  972. JGOBL=12
  973. TABOBL(1) ='YOUN'
  974. TABOBL(2) ='NU '
  975. TABOBL(3) ='RHO '
  976. TABOBL(4) ='ALPH'
  977. TABOBL(5) ='COB '
  978. TABOBL(6) ='MOB '
  979. TABOBL(7) ='ALPM'
  980. TABOBL(8) ='PERM'
  981. TABOBL(9) ='VISC'
  982. TABOBL(10)='KF '
  983. TABOBL(11)='RHOF'
  984. TABOBL(12)='ALPF'
  985. *
  986. IF (IFOUR.EQ.-2) THEN
  987. JGFAC=1
  988. TABFAC(1)='DIM3'
  989. ENDIF
  990. ENDIF
  991.  
  992. ELSEIF (IPLAC.EQ.5) THEN
  993. IF (MFR.EQ.37) THEN
  994. * MILIEU HOMOGENEISE
  995. *
  996. JGOBL=16
  997. TABOBL( 1)='B11 '
  998. TABOBL( 2)='B22 '
  999. TABOBL( 3)='B12 '
  1000. TABOBL( 4)='ROF '
  1001. TABOBL( 5)='ROS '
  1002. TABOBL( 6)='YOUN'
  1003. TABOBL( 7)='CSON'
  1004. TABOBL( 8)='RORF'
  1005. TABOBL( 9)='CREF'
  1006. TABOBL(10)='LCAR'
  1007. TABOBL(11)='E111'
  1008. TABOBL(12)='E112'
  1009. TABOBL(13)='E121'
  1010. TABOBL(14)='E122'
  1011. TABOBL(15)='E221'
  1012. TABOBL(16)='E222'
  1013. ENDIF
  1014. ELSEIF (IPLAC.EQ.6) THEN
  1015. C
  1016. C MATERIAU ELASTIQUE UNIDIRECTIONNEL
  1017. C
  1018. IF ((MFR.EQ.1.OR.MFR.EQ.31).AND.IDIM.EQ.3) THEN
  1019. JGOBL=7
  1020. TABOBL(1)='YOUN'
  1021. TABOBL(2)='V1X '
  1022. TABOBL(3)='V1Y '
  1023. TABOBL(4)='V1Z '
  1024. TABOBL(5)='V2X '
  1025. TABOBL(6)='V2Y '
  1026. TABOBL(7)='V2Z '
  1027.  
  1028. ELSE
  1029. JGOBL=3
  1030. TABOBL(1)='YOUN'
  1031. TABOBL(2)='V1X '
  1032. TABOBL(3)='V1Y '
  1033. ENDIF
  1034.  
  1035. IF((MFR.EQ.1.OR.MFR.EQ.3.OR.MFR.EQ.31).AND.
  1036. . IFOUR.EQ.-2) THEN
  1037. JGFAC=3
  1038. TABFAC(3)='DIM3'
  1039. ELSE
  1040. JGFAC=2
  1041. ENDIF
  1042. TABFAC(1)='RHO '
  1043. TABFAC(2)='ALPH'
  1044.  
  1045. ELSEIF (IPLAC.EQ.7) THEN
  1046. C MODELE ET MATERIAU DE LA SECTION DU MODELE A FIBRE
  1047. C
  1048. JGOBL=2
  1049. TABOBL(1)='MODS'
  1050. TABOBL(2)='MATS'
  1051.  
  1052. JGFAC=1
  1053. TABFAC(1)='MANO'
  1054. C ---- IPLAC 8: Ajoute par Jiang, 22/08/1995 et fleuret 28/05/96
  1055.  
  1056. ELSEIF (IPLAC.EQ.8) THEN
  1057. C MODELE ET MATERIAU POUR LE CABLE PRECONTRAINT
  1058. C
  1059. JGOBL=1
  1060. TABOBL(1)='YOUN'
  1061.  
  1062. JGFAC=8
  1063. TABFAC(1)='FF'
  1064. TABFAC(2)='PHIF'
  1065. TABFAC(3)='GANC'
  1066. TABFAC(4)='RMU0'
  1067. TABFAC(5)='FPRG'
  1068. TABFAC(6)='RH10'
  1069. TABFAC(7)='ALPH'
  1070. TABFAC(8)='RHO '
  1071.  
  1072. ELSEIF (IPLAC.EQ.9) THEN
  1073. C MODAL
  1074. C
  1075. JGOBL=3
  1076. TABOBL(1)='FREQ'
  1077. TABOBL(2)='MASS'
  1078. TABOBL(3)='DEFO'
  1079.  
  1080. JGFAC = 9
  1081. TABFAC(1) = 'AMOR'
  1082. TABFAC(2) = 'CGRA'
  1083. TABFAC(3) = 'MADE'
  1084. TABFAC(4) = 'RICR'
  1085. TABFAC(5) = 'MAIB'
  1086. TABFAC(6) = 'MACR'
  1087. TABFAC(7) = 'AMCR'
  1088. TABFAC(8) = 'ALP0'
  1089. TABFAC(9) = 'ECRO'
  1090.  
  1091. ELSEIF (IPLAC.EQ.10) THEN
  1092. C STATIQUE
  1093. C
  1094. JGOBL=3
  1095. TABOBL(1)='DEFO'
  1096. TABOBL(2)='RIDE'
  1097. TABOBL(3)='MADE'
  1098.  
  1099. JGFAC = 8
  1100. TABFAC(1) = 'AMOR'
  1101. TABFAC(2) = 'RICR'
  1102. TABFAC(3) = 'MAIA'
  1103. TABFAC(4) = 'MAIB'
  1104. TABFAC(5) = 'MACR'
  1105. TABFAC(6) = 'AMCR'
  1106. TABFAC(7) = 'BET0'
  1107. TABFAC(8) = 'ECRO'
  1108.  
  1109. ELSEIF (IPLAC.EQ.11) THEN
  1110. C ZONE_COHESIVE
  1111. C
  1112. JGOBL=2
  1113. TABOBL(1)='KS'
  1114. TABOBL(2)='KN'
  1115. ENDIF
  1116.  
  1117. INMAT=3
  1118. GOTO 20
  1119.  
  1120. ELSE
  1121. INMAT=2
  1122. GOTO 10
  1123. ENDIF
  1124. ENDIF
  1125. INMAT=0
  1126. 10 CONTINUE
  1127. *
  1128. * MATERIAU ELASTIQUE ISOTROPE
  1129. *
  1130. IF (MFR.EQ.35) THEN
  1131. * VALABLE EN 2D COMME EN 3D
  1132. JGOBL=2
  1133. TABOBL(1)='KS '
  1134. TABOBL(2)='KN '
  1135.  
  1136. JGFAC=2
  1137. TABFAC(1)='RHO '
  1138. TABFAC(2)='ALPN'
  1139. *
  1140. ELSEIF (MFR.EQ.78) THEN
  1141. * VALABLE EN 2D COMME EN 3D
  1142. JGOBL=2
  1143. TABOBL(1)='KS '
  1144. TABOBL(2)='KN '
  1145. *
  1146. * JOINT CISAILLEMENT (2D)
  1147. *
  1148. ELSEIF (MFR.EQ.53) THEN
  1149. * VALABLE EN 2D "COMME EN 3D"
  1150.  
  1151. JGOBL=1
  1152. TABOBL(1)='KS '
  1153. *
  1154. JGFAC=2
  1155. TABFAC(1)='RHO '
  1156. TABFAC(2)='ALPN'
  1157. *
  1158. * TOUS LES AUTRES CAS
  1159. *
  1160. ELSE
  1161. JGOBL=2
  1162. TABOBL(1)='YOUN'
  1163. TABOBL(2)='NU '
  1164. *
  1165. IF((MFR.EQ.1.OR.MFR.EQ.3.OR.MFR.EQ.31.OR.MFR.EQ.63).AND.
  1166. . IFOUR.EQ.-2) THEN
  1167. JGFAC=4
  1168. TABFAC(3)='DIM3'
  1169. TABFAC(4)='VISQ'
  1170. ELSE
  1171. JGFAC=3
  1172. TABFAC(3)='VISQ'
  1173. ENDIF
  1174. TABFAC(1)='RHO '
  1175. TABFAC(2)='ALPH'
  1176.  
  1177. * AM CAS FORMULATION NON-LOCALE
  1178.  
  1179. INLOC = 0
  1180. MN3 = INFMOD(/1)
  1181. IF (MN3.GE.13) INLOC=-1*INFMOD(13)
  1182. IF (INLOC.GT.0) THEN
  1183. JGOBL0 = JGFAC
  1184. * moyenne
  1185. IF(INLOC.EQ.1) THEN
  1186. JGFAC=JGOBL0+1
  1187. TABFAC(JGOBL0+1)='LCAR'
  1188. * stress-based
  1189. ELSEIF(INLOC.EQ.2) THEN
  1190. JGFAC=JGOBL0+2
  1191. TABFAC(JGOBL0+1)='LCAR'
  1192. TABFAC(JGOBL0+2)='SBFT'
  1193.  
  1194. * helmholtz
  1195. ELSE IF(INLOC.EQ.3) THEN
  1196. * sellier 03/07/20
  1197. * JGFAC=JGOBL0+6
  1198. * TABFAC(JGOBL0+1)='LCAR'
  1199. * TABFAC(JGOBL0+2)='LCF1'
  1200. * TABFAC(JGOBL0+3)='LCF2'
  1201. * TABFAC(JGOBL0+4)='LCF3'
  1202. * TABFAC(JGOBL0+5)='LCF4'
  1203. * TABFAC(JGOBL0+6)='LCF5'
  1204. JGFAC=JGOBL0+NB_PARA_HELM
  1205. c variables de Helmholtz a diffusions ORThotropes
  1206. c modif liste para heilmholtz sellier 10//11/22
  1207. c print*,'Dans IDMATR ', NB_HELM, 'formulationss Helmholtz'
  1208. do I=1,NB_HELM
  1209. c numero de la variable de Helmholtz
  1210. write (motnl1,'(I1)') I
  1211. c Capacité (infinie si BH=1)
  1212. motnl3='CAP'
  1213. motnl4=motnl3//motnl1
  1214. TABFAC(JGOBL0+NB_PARA_PAR_HELM*(I-1)+1)=motnl4
  1215. c print*,motnl4
  1216. c Blocage de Dirichlet si 1 sinon 0
  1217. motnl3='BLO'
  1218. motnl4=motnl3//motnl1
  1219. TABFAC(JGOBL0+NB_PARA_PAR_HELM*(I-1)+2)=motnl4
  1220. c print*,motnl4
  1221. c Valeur imposee a la zone de blocage de Dirichlet
  1222. motnl3='DEP'
  1223. motnl4=motnl3//motnl1
  1224. TABFAC(JGOBL0+NB_PARA_PAR_HELM*(I-1)+3)=motnl4
  1225. c print*,motnl4
  1226. c Valeur initiale de la variable de Helmholtz
  1227. motnl3='INI'
  1228. motnl4=motnl3//motnl1
  1229. TABFAC(JGOBL0+NB_PARA_PAR_HELM*(I-1)+4)=motnl4
  1230. c print*,motnl4
  1231. c indicateur logique de lineratite pour ne pas passer par istep=3 si 1
  1232. motnl3='LIN'
  1233. motnl4=motnl3//motnl1
  1234. TABFAC(JGOBL0+NB_PARA_PAR_HELM*(I-1)+5)=motnl4
  1235. c print*,motnl4
  1236. do J=1,3
  1237. c coefficients de diffusion
  1238. write (motnld1,'(I1)') J
  1239. c diffusion
  1240. motnl2='DH'
  1241. c 1er indice: numero de la variable d helmholtz
  1242. motnl3=motnl2//motnl1
  1243. c 2eme indice : direction
  1244. motnl4=motnl3//motnld1
  1245. TABFAC(JGOBL0+(I-1)*NB_PARA_PAR_HELM+
  1246. # 5+(J-1)*4+1)=motnl4
  1247. c print*,motnl4
  1248. do K=1,3
  1249. write (motnla1,'(I1)') K
  1250. c direction principales de diffusion
  1251. motnlv1='V'
  1252. c indice numero de la variable
  1253. motnl2=motnlv1//motnl1
  1254. c indice direction principale
  1255. motnl3=motnl2//motnld1
  1256. c projection sur axe indice de base fixe
  1257. motnl4=motnl3//motnla1
  1258. TABFAC(JGOBL0+(I-1)*NB_PARA_PAR_HELM+
  1259. # 5+(J-1)*4+1+K)=motnl4
  1260. c print*,motnl4
  1261. end do
  1262. end do
  1263. end do
  1264. Cc on rajoute 2 inutiles pour compenser un decalage constatte dans castem20
  1265. C TABFAC(JGOBL0+NB_PARA_HELM+1)='FAC1'
  1266. C TABFAC(JGOBL0+NB_PARA_HELM+2)='FAC2'
  1267. * fin modif sellier
  1268. * cas non prevu
  1269. ELSE
  1270. IRET = 0
  1271. CALL ERREUR(5)
  1272. ENDIF
  1273. ENDIF
  1274. *
  1275. ENDIF
  1276. *
  1277. IF (INMAT.EQ.0) GOTO 9999
  1278. *
  1279. 19 CONTINUE
  1280. 20 CONTINUE
  1281. DO jm = 1,matmod(/2)
  1282. IF (matmod(jm).eq.'IMPEDANCE') THEN
  1283. imate = imatee
  1284. INMAT = INMAT + 1
  1285. JGOBL = 0
  1286. JGFAC = 0
  1287.  
  1288. IF(CMATEE.EQ.'IMPELAST') THEN
  1289. *IMPE_ELAS
  1290. JGOBL0 = JGOBL
  1291. JGOBL= JGOBL0+1
  1292. TABOBL(JGOBL0 + 1) ='RAID'
  1293. JGOBL0 = JGFAC
  1294.  
  1295. JGFAC = JGOBL0+4
  1296. TABFAC(JGOBL0+1) = 'AMOR'
  1297. TABFAC(JGOBL0+2) = 'ZNU'
  1298. TABFAC(JGOBL0+3) = 'MASS'
  1299. TABFAC(JGOBL0+4) = 'ALPH'
  1300. GOTO 22
  1301. ELSEIF (CMATEE.EQ.'IMPVOIGT'.or.CMATEE.EQ.'IMPREUSS') THEN
  1302. *IMPE_VOIGT ou IMPE_REUSS
  1303. JGOBL0 = JGOBL
  1304. JGOBL= JGOBL0+2
  1305. TABOBL(JGOBL0 + 1) ='RAID'
  1306. TABOBL(JGOBL0 + 2) ='VISC'
  1307.  
  1308. JGOBL0 = JGFAC
  1309. JGFAC = JGOBL0+2
  1310. TABFAC(JGOBL0+1) = 'MASS'
  1311. TABFAC(JGOBL0+2) = 'AMOR'
  1312. GOTO 22
  1313. ELSEIF (CMATEE.EQ.'IMPCOMPL') then
  1314. *IMPE_COMPLEXE
  1315. JGOBL0 = JGOBL
  1316. JGOBL= JGOBL0+1
  1317. TABOBL(JGOBL0 + 1) ='RAID'
  1318.  
  1319. JGOBL0 = JGFAC
  1320. JGFAC = JGFAC+4
  1321. TABFAC(JGOBL0+1) = 'MOCO'
  1322. TABFAC(JGOBL0+2) = 'VISC'
  1323. TABFAC(JGOBL0+3) = 'MASS'
  1324. TABFAC(JGOBL0+4) = 'AMOR'
  1325. GOTO 22
  1326. ELSE
  1327. ENDIF
  1328. 22 CONTINUE
  1329. * if (mele.eq.45) then
  1330. JGOBL0 = JGFAC
  1331. JGFAC = JGOBL0+3
  1332.  
  1333. IF(CMATEE.EQ.'IMPELAST') THEN
  1334. TABFAC(JGOBL0+1) = TABFAC(JGOBL0-2)
  1335. TABFAC(JGOBL0+2) = TABFAC(JGOBL0-1)
  1336. TABFAC(JGOBL0+3) = TABFAC(JGOBL0)
  1337. TABFAC(JGOBL0-2) = 'CPLE'
  1338. TABFAC(JGOBL0-1) = 'INER'
  1339. TABFAC(JGOBL0) = 'AROT'
  1340. ELSE
  1341. TABFAC(JGOBL0+1) = 'CPLE'
  1342. TABFAC(JGOBL0+2) = 'INER'
  1343. TABFAC(JGOBL0+3) = 'AROT'
  1344. ENDIF
  1345. * endif
  1346. IF (NMAT.GE.INMAT) inmat = inmat+1
  1347. * GOTO 9999
  1348. ENDIF
  1349. ENDDO
  1350. *
  1351. IF (NMAT.GE.INMAT) THEN
  1352. CALL MODNLI(MOMODL,NMOD)
  1353. CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  1354. IF (IPLAC.EQ.1) THEN
  1355. INMAT=INMAT+1
  1356. CALL MODPLA(MOMODL,NMOD)
  1357. CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  1358. IF (IPLAC.EQ.0) THEN
  1359. * ISOTROPE
  1360. IPLAC=1
  1361. * ELSE
  1362. * INMAT=INMAT+1
  1363. ENDIF
  1364. CALL IDPLAS(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  1365. ELSEIF (IPLAC.EQ.2) THEN
  1366. INMAT=INMAT+1
  1367. CALL MODFLU(MOMODL,NMOD)
  1368. CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  1369. IF (IPLAC.EQ.0) THEN
  1370. * NORTON
  1371. IPLAC=1
  1372. * ELSE
  1373. * INMAT=INMAT+1
  1374. ENDIF
  1375. CALL IDFLUA(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  1376. ELSEIF (IPLAC.EQ.3) THEN
  1377. INMAT=INMAT+1
  1378. CALL MODVIS(MOMODL,NMOD)
  1379. CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  1380. IF (IPLAC.EQ.0) THEN
  1381. * ONERA
  1382. IPLAC=2
  1383. * ELSE
  1384. * INMAT=INMAT+1
  1385. ENDIF
  1386. C
  1387. C Cas particulier des modeles GATT_MONERIE et UO2
  1388. C 'RHO ' et 'ALPH' sont obligatoires
  1389. IF (IPLAC.EQ.18.OR.IPLAC.EQ.19) THEN
  1390. JGOBL0 = JGOBL
  1391. JGOBL=JGOBL0+2
  1392. TABOBL(JGOBL0+1)='RHO '
  1393. TABOBL(JGOBL0+2)='ALPH'
  1394.  
  1395. JGOBL0=JGFAC
  1396. JGFAC=JGOBL0-2
  1397. IF (JGFAC.GT.0) THEN
  1398. DO 200 I=1,JGFAC
  1399. TABFAC(I)=TABFAC(I+2)
  1400. 200 CONTINUE
  1401. ENDIF
  1402.  
  1403. C ===
  1404. C Modeles SYMONDS & COWPER SYCO1 et SYCO2
  1405. C ===
  1406. ELSEIF (IPLAC.EQ.28) THEN
  1407. JGOBL0=JGOBL
  1408. JGOBL=JGOBL0+1
  1409. TABOBL(JGOBL0+1)='ECRO'
  1410.  
  1411. ELSEIF (IPLAC.EQ.29) THEN
  1412. JGOBL0=JGOBL
  1413. JGOBL=JGOBL0+1
  1414. TABOBL(JGOBL0+1)='ECRO'
  1415. ENDIF
  1416. CALL IDVISC(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  1417. C
  1418. ELSEIF (IPLAC.EQ.4) THEN
  1419. INMAT=INMAT+1
  1420. CALL MODEND(MOMODL,NMOD)
  1421. CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  1422. IF (IPLAC.EQ.0) THEN
  1423. * MAZARS
  1424. IPLAC=1
  1425. * ELSE
  1426. * INMAT=INMAT+1
  1427. ENDIF
  1428. CALL IDENDO(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  1429. ELSEIF (IPLAC.EQ.5) THEN
  1430. INMAT=INMAT+1
  1431. CALL MODPLE(MOMODL,NMOD)
  1432. CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  1433. IF (IPLAC.EQ.0) THEN
  1434. * TRIAXIAL P/Y
  1435. IPLAC=1
  1436. * ELSE
  1437. * INMAT=INMAT+1
  1438. ENDIF
  1439. CALL IDPLEN(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  1440. ELSEIF (IPLAC.EQ.6) THEN
  1441. INMAT=INMAT+1
  1442. CALL MODENL(MOMODL,NMOD)
  1443. CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  1444. IF (IPLAC.EQ.0) THEN
  1445. * ISOTROPE
  1446. IPLAC=1
  1447. * ELSE
  1448. * INMAT=INMAT+1
  1449. ENDIF
  1450. CALL IDELNL(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  1451. ELSEIF (IPLAC.EQ.7) THEN
  1452. C Pas de composantes supplementaires pour une loi 'VISCO_EXTERNE'
  1453. IRET = 1
  1454. ELSEIF(IPLAC.EQ.0.AND.IMATEE.EQ.31) THEN
  1455. GOTO 9999
  1456. ELSE
  1457. GOTO 9999
  1458. ENDIF
  1459. C IRET code retour de IDPLAS,IDFLUA,IDVISC,IDENDO,IDPLEN,IDELNL
  1460. IF (IRET.EQ.0) GOTO 9999
  1461. ENDIF
  1462.  
  1463. * Parametres CRIP et FUSION :
  1464. do jma=1,matmod(/2)
  1465. if(matmod(jma).eq.'CRIP') then
  1466. JGOBL = JGOBL + 1
  1467. TABOBL(JGOBL) = 'LIMP'
  1468.  
  1469. elseif(matmod(jma).eq.'FUSION') then
  1470. JGOBL = JGOBL + 1
  1471. TABOBL(JGOBL) = 'TFUS'
  1472.  
  1473. endif
  1474. enddo
  1475. *
  1476. GOTO 9999
  1477. ENDIF
  1478.  
  1479. *--------------------------------------------------------------------
  1480. * CAS DE LA FORMULATION POREUX
  1481. *--------------------------------------------------------------------
  1482. CALL PLACE(FORMOD,NFOR,IPORE,'POREUX')
  1483. IF (IPORE.NE.0) THEN
  1484. *
  1485. CALL MODLIN(MOMODL,NMOD)
  1486. CALL PLACE(MOMODL,NMOD,IRET,MATMOD(1))
  1487. IF (IRET.EQ.0) GOTO 9999
  1488. *
  1489. * D'ABORD : CAS NON ISOTROPE
  1490. *
  1491. IF (NMAT.GE.2) THEN
  1492. CALL MODELA(MOMODL,NMOD)
  1493. CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(2))
  1494. IF (IPLAC.NE.0) THEN
  1495. INMAT=3
  1496. IF (IPLAC.EQ.1) THEN
  1497. INMAT=3
  1498. GOTO 30
  1499. ENDIF
  1500. IF (IPLAC.EQ.2) THEN
  1501. *
  1502. * MATERIAU ELASTIQUE ORTHOTROPE
  1503. *
  1504. IF (MFR.EQ.33) THEN
  1505. *
  1506. * ELEMENTS MASSIFS
  1507. *
  1508. IF(IDIM.EQ.3)THEN
  1509. * ELEMENTS 3D
  1510. JGOBL=19
  1511. TABOBL(1)='YG1 '
  1512. TABOBL(2)='YG2 '
  1513. TABOBL(3)='YG3 '
  1514. TABOBL(4)='NU12'
  1515. TABOBL(5)='NU23'
  1516. TABOBL(6)='NU13'
  1517. TABOBL(7)='G12 '
  1518. TABOBL(8)='G23 '
  1519. TABOBL(9)='G13 '
  1520. TABOBL(10)='V1X '
  1521. TABOBL(11)='V1Y '
  1522. TABOBL(12)='V1Z '
  1523. TABOBL(13)='V2X '
  1524. TABOBL(14)='V2Y '
  1525. TABOBL(15)='V2Z '
  1526. TABOBL(16)='COB1'
  1527. TABOBL(17)='COB2'
  1528. TABOBL(18)='COB3'
  1529. TABOBL(19)='MOB '
  1530. *
  1531. JGFAC=12
  1532. TABFAC(1)='RHO '
  1533. TABFAC(2)='ALP1'
  1534. TABFAC(3)='ALP2'
  1535. TABFAC(4)='ALP3'
  1536. TABFAC(5)='ALPM'
  1537. TABFAC(6)='PER1'
  1538. TABFAC(7)='PER2'
  1539. TABFAC(8)='PER3'
  1540. TABFAC(9)='VISC'
  1541. TABFAC(10)='KF '
  1542. TABFAC(11)='RHOF'
  1543. TABFAC(12)='ALPF'
  1544. ELSEIF (IDIM.EQ.2) THEN
  1545. IF(IFOUR.EQ.-2)THEN
  1546. * CONTRAINTE PLANE
  1547. JGOBL=12
  1548. TABOBL(1)='YG1 '
  1549. TABOBL(2)='YG2 '
  1550. TABOBL(3)='NU12'
  1551. TABOBL(4)='G12 '
  1552. TABOBL(5)='V1X '
  1553. TABOBL(6)='V1Y '
  1554. TABOBL(7)='COB1'
  1555. TABOBL(8)='COB2'
  1556. TABOBL(9)='MOB '
  1557. TABOBL(10)='YG3 '
  1558. TABOBL(11)='NU23'
  1559. TABOBL(12)='NU13'
  1560. *
  1561. JGFAC=11
  1562. TABFAC(1)='RHO '
  1563. TABFAC(2)='ALP1'
  1564. TABFAC(3)='ALP2'
  1565. TABFAC(4)='ALPM'
  1566. TABFAC(5)='PER1'
  1567. TABFAC(6)='PER2'
  1568. TABFAC(7)='VISC'
  1569. TABFAC(8)='KF '
  1570. TABFAC(9)='RHOF'
  1571. TABFAC(10)='ALPF'
  1572. TABFAC(11)='DIM3'
  1573. ELSEIF (IFOUR.EQ.-1.OR.IFOUR.EQ.0.OR.IFOUR.EQ.-3) THEN
  1574. * DEFORMATION PLANE ,AXISYMETRIE
  1575. JGOBL=13
  1576. TABOBL(1)='YG1 '
  1577. TABOBL(2)='YG2 '
  1578. TABOBL(3)='YG3 '
  1579. TABOBL(4)='NU12'
  1580. TABOBL(5)='NU23'
  1581. TABOBL(6)='NU13'
  1582. TABOBL(7)='G12 '
  1583. TABOBL(8)='V1X '
  1584. TABOBL(9)='V1Y '
  1585. TABOBL(10)='COB1'
  1586. TABOBL(11)='COB2'
  1587. TABOBL(12)='COB3'
  1588. TABOBL(13)='MOB '
  1589. *
  1590. JGFAC=11
  1591. TABFAC(1)='RHO '
  1592. TABFAC(2)='ALP1'
  1593. TABFAC(3)='ALP2'
  1594. TABFAC(4)='ALP3'
  1595. TABFAC(5)='ALPM'
  1596. TABFAC(6)='PER1'
  1597. TABFAC(7)='PER2'
  1598. TABFAC(8)='VISC'
  1599. TABFAC(9)='KF '
  1600. TABFAC(10)='RHOF'
  1601. TABFAC(11)='ALPF'
  1602. *
  1603. ELSEIF (IFOUR.EQ.1) THEN
  1604. * AXISYMETRIE DE FOURIER
  1605. JGOBL=15
  1606. TABOBL(1)='YG1 '
  1607. TABOBL(2)='YG2 '
  1608. TABOBL(3)='YG3 '
  1609. TABOBL(4)='NU12'
  1610. TABOBL(5)='NU23'
  1611. TABOBL(6)='NU13'
  1612. TABOBL(7)='G12 '
  1613. TABOBL(8)='G23 '
  1614. TABOBL(9)='G13 '
  1615. TABOBL(10)='V1X '
  1616. TABOBL(11)='V1Y '
  1617. TABOBL(12)='COB1'
  1618. TABOBL(13)='COB2'
  1619. TABOBL(14)='COB3'
  1620. TABOBL(15)='MOB '
  1621. *
  1622. JGFAC=12
  1623. TABFAC(1)='RHO '
  1624. TABFAC(2)='ALP1'
  1625. TABFAC(3)='ALP2'
  1626. TABFAC(4)='ALP3'
  1627. TABFAC(5)='ALPM'
  1628. TABFAC(6)='PER1'
  1629. TABFAC(7)='PER2'
  1630. TABFAC(8)='PER3'
  1631. TABFAC(9)='VISC'
  1632. TABFAC(10)='KF '
  1633. TABFAC(11)='RHOF'
  1634. TABFAC(12)='ALPF'
  1635. ENDIF
  1636. ENDIF
  1637. ELSE
  1638. *
  1639. * CAS NON PREVU
  1640. *
  1641. IRET = 0
  1642. GOTO 9999
  1643. ENDIF
  1644. ELSEIF (IPLAC.EQ.3)THEN
  1645. *
  1646. * MATERIAU ANISOTROPE ELASTIQUE
  1647. *
  1648. IF(MFR.EQ.33)THEN
  1649. *
  1650. * ELEMENTS MASSIFS
  1651. *
  1652. IF(IDIM.EQ.3)THEN
  1653. * ELEMENTS 3D
  1654. JGOBL=34
  1655. TABOBL(1)='D11 '
  1656. TABOBL(2)='D21 '
  1657. TABOBL(3)='D22 '
  1658. TABOBL(4)='D31 '
  1659. TABOBL(5)='D32 '
  1660. TABOBL(6)='D33 '
  1661. TABOBL(7)='D41 '
  1662. TABOBL(8)='D42 '
  1663. TABOBL(9)='D43 '
  1664. TABOBL(10)='D44 '
  1665. TABOBL(11)='D51 '
  1666. TABOBL(12)='D52 '
  1667. TABOBL(13)='D53 '
  1668. TABOBL(14)='D54 '
  1669. TABOBL(15)='D55 '
  1670. TABOBL(16)='D61 '
  1671. TABOBL(17)='D62 '
  1672. TABOBL(18)='D63 '
  1673. TABOBL(19)='D64 '
  1674. TABOBL(20)='D65 '
  1675. TABOBL(21)='D66 '
  1676. TABOBL(22)='V1X '
  1677. TABOBL(23)='V1Y '
  1678. TABOBL(24)='V1Z '
  1679. TABOBL(25)='V2X '
  1680. TABOBL(26)='V2Y '
  1681. TABOBL(27)='V2Z '
  1682. TABOBL(28)='COB1'
  1683. TABOBL(29)='COB2'
  1684. TABOBL(30)='COB3'
  1685. TABOBL(31)='CO12'
  1686. TABOBL(32)='CO13'
  1687. TABOBL(33)='CO23'
  1688. TABOBL(34)='MOB '
  1689. *
  1690. JGFAC=18
  1691. TABFAC(1)='RHO '
  1692. TABFAC(2)='ALP1'
  1693. TABFAC(3)='ALP2'
  1694. TABFAC(4)='ALP3'
  1695. TABFAC(5)='AL12'
  1696. TABFAC(6)='AL13'
  1697. TABFAC(7)='AL23'
  1698. TABFAC(8)='ALPM'
  1699. TABFAC(9 )='PER1'
  1700. TABFAC(10)='PER2'
  1701. TABFAC(11)='PER3'
  1702. TABFAC(12)='PE12'
  1703. TABFAC(13)='PE13'
  1704. TABFAC(14)='PE23'
  1705. TABFAC(15)='VISC'
  1706. TABFAC(16)='KF '
  1707. TABFAC(17)='RHOF'
  1708. TABFAC(18)='ALPF'
  1709. ELSEIF (IDIM.EQ.2) THEN
  1710. IF (IFOUR.EQ.-2) THEN
  1711. * CONTRAINTE PLANE
  1712. JGOBL=16
  1713. TABOBL(1)='D11 '
  1714. TABOBL(2)='D21 '
  1715. TABOBL(3)='D22 '
  1716. TABOBL(4)='D41 '
  1717. TABOBL(5)='D42 '
  1718. TABOBL(6)='D44 '
  1719. TABOBL(7)='V1X '
  1720. TABOBL(8)='V1Y '
  1721. TABOBL(9 )='COB1'
  1722. TABOBL(10)='COB2'
  1723. TABOBL(11)='CO12'
  1724. TABOBL(12)='MOB '
  1725. TABOBL(13)='D31 '
  1726. TABOBL(14)='D32 '
  1727. TABOBL(15)='D33 '
  1728. TABOBL(16)='D43 '
  1729. *
  1730. JGFAC=13
  1731. TABFAC(1)='RHO '
  1732. TABFAC(2)='ALP1'
  1733. TABFAC(3)='ALP2'
  1734. TABFAC(4)='AL12'
  1735. TABFAC(5)='ALPM'
  1736. TABFAC(6)='PER1'
  1737. TABFAC(7)='PER2'
  1738. TABFAC(8)='PE12'
  1739. TABFAC(9)='VISC'
  1740. TABFAC(10)='KF '
  1741. TABFAC(11)='RHOF'
  1742. TABFAC(12)='ALPF'
  1743. TABFAC(13)='DIM3'
  1744. *
  1745. ELSEIF (IFOUR.EQ.-1.OR.IFOUR.EQ.0.OR.IFOUR.EQ.-3) THEN
  1746. * DEFORMATION PLANE ,AXISYMETRIE
  1747. JGOBL=17
  1748. TABOBL(1)='D11 '
  1749. TABOBL(2)='D21 '
  1750. TABOBL(3)='D22 '
  1751. TABOBL(4)='D31 '
  1752. TABOBL(5)='D32 '
  1753. TABOBL(6)='D33 '
  1754. TABOBL(7)='D41 '
  1755. TABOBL(8)='D42 '
  1756. TABOBL(9)='D43 '
  1757. TABOBL(10)='D44 '
  1758. TABOBL(11)='V1X '
  1759. TABOBL(12)='V1Y '
  1760. TABOBL(13)='COB1'
  1761. TABOBL(14)='COB2'
  1762. TABOBL(15)='CO12'
  1763. TABOBL(16)='COB3'
  1764. TABOBL(17)='MOB '
  1765. *
  1766. JGFAC=13
  1767. TABFAC(1)='RHO '
  1768. TABFAC(2)='ALP1'
  1769. TABFAC(3)='ALP2'
  1770. TABFAC(4)='AL12'
  1771. TABFAC(5)='ALP3'
  1772. TABFAC(6)='ALPM'
  1773. TABFAC(7)='PER1'
  1774. TABFAC(8)='PER2'
  1775. TABFAC(9)='PE12'
  1776. TABFAC(10)='VISC'
  1777. TABFAC(11)='KF '
  1778. TABFAC(12)='RHOF'
  1779. TABFAC(13)='ALPF'
  1780. *
  1781. ELSEIF (IFOUR.EQ.1) THEN
  1782. * AXISYMETRIE DE FOURIER
  1783. JGOBL=20
  1784. TABOBL(1)='D11 '
  1785. TABOBL(2)='D21 '
  1786. TABOBL(3)='D22 '
  1787. TABOBL(4)='D31 '
  1788. TABOBL(5)='D32 '
  1789. TABOBL(6)='D33 '
  1790. TABOBL(7)='D41 '
  1791. TABOBL(8)='D42 '
  1792. TABOBL(9)='D43 '
  1793. TABOBL(10)='D44 '
  1794. TABOBL(11)='D55 '
  1795. TABOBL(12)='D65 '
  1796. TABOBL(13)='D66 '
  1797. TABOBL(14)='V1X '
  1798. TABOBL(15)='V1Y '
  1799. TABOBL(16)='COB1'
  1800. TABOBL(17)='COB2'
  1801. TABOBL(18)='CO12'
  1802. TABOBL(19)='COB3'
  1803. TABOBL(20)='MOB '
  1804. *
  1805. JGFAC=14
  1806. TABFAC(1)='RHO '
  1807. TABFAC(2)='ALP1'
  1808. TABFAC(3)='ALP2'
  1809. TABFAC(4)='AL12'
  1810. TABFAC(5)='ALP3'
  1811. TABFAC(6)='ALPM'
  1812. TABFAC(7)='PER1'
  1813. TABFAC(8)='PER2'
  1814. TABFAC(9)='PE12'
  1815. TABFAC(10)='PER3'
  1816. TABFAC(11)='VISC'
  1817. TABFAC(12)='KF '
  1818. TABFAC(13)='RHOF'
  1819. TABFAC(14)='ALPF'
  1820. ENDIF
  1821. ENDIF
  1822. ELSE
  1823. *
  1824. * CAS NON PREVU
  1825. *
  1826. IRET = 0
  1827. GOTO 9999
  1828. ENDIF
  1829.  
  1830. ELSEIF (IPLAC.EQ.6) THEN
  1831. C
  1832. C MATERIAU ELASTIQUE UNIDIRECTIONNEL
  1833. C
  1834. IF(MFR.EQ.33)THEN
  1835. IF (IDIM.EQ.3) THEN
  1836. JGOBL=9
  1837. TABOBL(1)='YOUN'
  1838. TABOBL(2)='V1X '
  1839. TABOBL(3)='V1Y '
  1840. TABOBL(4)='V1Z '
  1841. TABOBL(5)='V2X '
  1842. TABOBL(6)='V2Y '
  1843. TABOBL(7)='V2Z '
  1844. TABOBL(8)='COB '
  1845. TABOBL(9)='MOB '
  1846. ELSE
  1847. JGOBL=5
  1848. TABOBL(1)='YOUN'
  1849. TABOBL(2)='V1X '
  1850. TABOBL(3)='V1Y '
  1851. TABOBL(4)='COB '
  1852. TABOBL(5)='MOB '
  1853. ENDIF
  1854. *
  1855. IF(IFOUR.EQ.-2) THEN
  1856. JGFAC=9
  1857. TABFAC(9)='DIM3'
  1858. ELSE
  1859. JGFAC=8
  1860. ENDIF
  1861. TABFAC(1)='RHO '
  1862. TABFAC(2)='ALPH'
  1863. TABFAC(3)='ALPM'
  1864. TABFAC(4)='PERM'
  1865. TABFAC(5)='VISC'
  1866. TABFAC(6)='KF '
  1867. TABFAC(7)='RHOF'
  1868. TABFAC(8)='ALPF'
  1869. ELSE
  1870. *
  1871. * CAS NON PREVU
  1872. *
  1873. IRET = 0
  1874. GO TO 9999
  1875. ENDIF
  1876. *
  1877. ENDIF
  1878. INMAT=3
  1879. GOTO 40
  1880. ELSE
  1881. INMAT=2
  1882. GOTO 30
  1883. ENDIF
  1884. ENDIF
  1885. INMAT=0
  1886. 30 CONTINUE
  1887. *
  1888. * CAS MATERIAU POREUX ELASTIQUE ISOTROPE
  1889. *
  1890. IF(MELE.GE.79.AND.MELE.LE.83)THEN
  1891. *
  1892. JGOBL=4
  1893. TABOBL(1) ='YOUN'
  1894. TABOBL(2) ='NU '
  1895. TABOBL(3) ='COB '
  1896. TABOBL(4) ='MOB '
  1897. *
  1898. IF(IFOUR.EQ.-2) THEN
  1899. JGFAC=9
  1900. TABFAC(9)='DIM3'
  1901. ELSE
  1902. JGFAC=8
  1903. ENDIF
  1904. TABFAC(1)='RHOF'
  1905. TABFAC(2)='ALPF'
  1906. TABFAC(3)='ALPM'
  1907. TABFAC(4)='PERM'
  1908. TABFAC(5)='VISC'
  1909. TABFAC(6)='KF '
  1910. TABFAC(7)='RHO '
  1911. TABFAC(8)='ALPH'
  1912. *
  1913. ELSEIF(MELE.GE.108.AND.MELE.LE.110)THEN
  1914. *
  1915. * CAS DES JOINTS POREUX ISOTROPES (VALABLE EN 2D COMME EN 3D)
  1916. *
  1917. JGOBL=4
  1918. TABOBL(1)='KS '
  1919. TABOBL(2)='KN '
  1920. TABOBL(3)='COB '
  1921. TABOBL(4)='MOB '
  1922. *
  1923. JGFAC=4
  1924. TABFAC(1)='PERT'
  1925. TABFAC(2)='PERH'
  1926. TABFAC(3)='PERB'
  1927. TABFAC(4)='VISC'
  1928. *
  1929. ELSEIF(MELE.GE.173.AND.MELE.LE.177)THEN
  1930. *
  1931. JGOBL=10
  1932. TABOBL(1) ='YOUN'
  1933. TABOBL(2) ='NU '
  1934. TABOBL(3) ='COP1'
  1935. TABOBL(4) ='COP2'
  1936. TABOBL(5) ='CPP1'
  1937. TABOBL(6) ='CPP2'
  1938. TABOBL(7) ='KK11'
  1939. TABOBL(8) ='KK12'
  1940. TABOBL(9) ='KK21'
  1941. TABOBL(10)='KK22'
  1942. *
  1943. IF(IFOUR.EQ.-2) THEN
  1944. JGFAC=9
  1945. TABFAC(9)='DIM3'
  1946. ELSE
  1947. JGFAC=8
  1948. ENDIF
  1949. TABFAC(1)='RHOF'
  1950. TABFAC(2)='ALPF'
  1951. TABFAC(3)='ALPM'
  1952. TABFAC(4)='PK11'
  1953. TABFAC(5)='PK12'
  1954. TABFAC(6)='PK21'
  1955. TABFAC(7)='PK22'
  1956. TABFAC(8)='ALPH'
  1957. *
  1958. ELSEIF(MELE.GE.185.AND.MELE.LE.187)THEN
  1959. *
  1960. * CAS DES JOINTS POREUX ISOTROPES (VALABLE EN 2D COMME EN 3D)
  1961. *
  1962. JGOBL=10
  1963. TABOBL(1)='KS '
  1964. TABOBL(2)='KN '
  1965. TABOBL(3)='COP1'
  1966. TABOBL(4)='COP2'
  1967. TABOBL(5)='CPP1'
  1968. TABOBL(6)='CPP2'
  1969. TABOBL(7)='KK11'
  1970. TABOBL(8)='KK12'
  1971. TABOBL(9)='KK21'
  1972. TABOBL(10)='KK22'
  1973. *
  1974. IF(IFOUR.EQ.-2) THEN
  1975. JGFAC=18
  1976. TABFAC(18)='DIM3'
  1977. ELSE
  1978. JGFAC=17
  1979. ENDIF
  1980. TABFAC(1)='RHOF'
  1981. TABFAC(2)='ALPF'
  1982. TABFAC(3)='ALPM'
  1983. TABFAC(4)='PT11'
  1984. TABFAC(5)='PH11'
  1985. TABFAC(6)='PB11'
  1986. TABFAC(7)='PT12'
  1987. TABFAC(8)='PH12'
  1988. TABFAC(9)='PB12'
  1989. TABFAC(10)='PT21'
  1990. TABFAC(11)='PH21'
  1991. TABFAC(12)='PB21'
  1992. TABFAC(13)='PT22'
  1993. TABFAC(14)='PH22'
  1994. TABFAC(15)='PB22'
  1995. TABFAC(16)='RHO '
  1996. TABFAC(17)='ALPH'
  1997.  
  1998. ELSEIF(MELE.GE.178.AND.MELE.LE.182)THEN
  1999. *
  2000. JGOBL=17
  2001. TABOBL(1)='YOUN'
  2002. TABOBL(2)='NU '
  2003. TABOBL(3)='COP1'
  2004. TABOBL(4)='COP2'
  2005. TABOBL(5)='COP3'
  2006. TABOBL(6)='CPP1'
  2007. TABOBL(7)='CPP2'
  2008. TABOBL(8)='CPP3'
  2009. TABOBL(9)='KK11'
  2010. TABOBL(10)='KK12'
  2011. TABOBL(11)='KK13'
  2012. TABOBL(12)='KK21'
  2013. TABOBL(13)='KK22'
  2014. TABOBL(14)='KK23'
  2015. TABOBL(15)='KK31'
  2016. TABOBL(16)='KK32'
  2017. TABOBL(17)='KK33'
  2018. *
  2019. IF(IFOUR.EQ.-2) THEN
  2020. JGFAC=15
  2021. TABFAC(15)='DIM3'
  2022. ELSE
  2023. JGFAC=14
  2024. ENDIF
  2025. TABFAC(1)='RHOF'
  2026. TABFAC(2)='ALPF'
  2027. TABFAC(3)='ALPM'
  2028. TABFAC(4)='PK11'
  2029. TABFAC(5)='PK12'
  2030. TABFAC(6)='PK13'
  2031. TABFAC(7)='PK21'
  2032. TABFAC(8)='PK22'
  2033. TABFAC(9)='PK23'
  2034. TABFAC(10)='PK31'
  2035. TABFAC(11)='PK32'
  2036. TABFAC(12)='PK33'
  2037. TABFAC(13)='RHO '
  2038. TABFAC(14)='ALPH'
  2039.  
  2040. ELSEIF(MELE.GE.188.AND.MELE.LE.190)THEN
  2041. *
  2042. * CAS DES JOINTS POREUX ISOTROPES (VALABLE EN 2D COMME EN 3D)
  2043. *
  2044. JGOBL=17
  2045. TABOBL(1)='KS '
  2046. TABOBL(2)='KN '
  2047. TABOBL(3)='COP1'
  2048. TABOBL(4)='COP2'
  2049. TABOBL(5)='COP3'
  2050. TABOBL(6)='CPP1'
  2051. TABOBL(7)='CPP2'
  2052. TABOBL(8)='CPP3'
  2053. TABOBL(9)='KK11'
  2054. TABOBL(10)='KK12'
  2055. TABOBL(11)='KK13'
  2056. TABOBL(12)='KK21'
  2057. TABOBL(13)='KK22'
  2058. TABOBL(14)='KK23'
  2059. TABOBL(15)='KK31'
  2060. TABOBL(16)='KK32'
  2061. TABOBL(17)='KK33'
  2062. *
  2063. IF(IFOUR.EQ.-2) THEN
  2064. JGFAC=33
  2065. TABFAC(33)='DIM3'
  2066. ELSE
  2067. JGFAC=32
  2068. ENDIF
  2069. TABFAC(1)='RHOF'
  2070. TABFAC(2)='ALPF'
  2071. TABFAC(3)='ALPM'
  2072. TABFAC(4)='PT11'
  2073. TABFAC(5)='PH11'
  2074. TABFAC(6)='PB11'
  2075. TABFAC(7)='PT12'
  2076. TABFAC(8)='PH12'
  2077. TABFAC(9)='PB12'
  2078. TABFAC(10)='PT13'
  2079. TABFAC(11)='PH13'
  2080. TABFAC(12)='PB13'
  2081. TABFAC(13)='PT21'
  2082. TABFAC(14)='PH21'
  2083. TABFAC(15)='PB21'
  2084. TABFAC(16)='PT22'
  2085. TABFAC(17)='PH22'
  2086. TABFAC(18)='PB22'
  2087. TABFAC(19)='PT23'
  2088. TABFAC(20)='PH23'
  2089. TABFAC(21)='PB23'
  2090. TABFAC(22)='PT31'
  2091. TABFAC(23)='PH31'
  2092. TABFAC(24)='PB31'
  2093. TABFAC(25)='PT32'
  2094. TABFAC(26)='PH32'
  2095. TABFAC(27)='PB32'
  2096. TABFAC(28)='PT33'
  2097. TABFAC(29)='PH33'
  2098. TABFAC(30)='PB33'
  2099. TABFAC(31)='RHO '
  2100. TABFAC(32)='ALPH'
  2101.  
  2102. ENDIF
  2103. *
  2104. IF (INMAT.EQ.0) THEN
  2105. ** IRET = 0
  2106. GOTO 9999
  2107. ENDIF
  2108. *
  2109. 40 CONTINUE
  2110. IF (NMAT.GE.INMAT) THEN
  2111. CALL MODNLI(MOMODL,NMOD)
  2112. CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  2113. IF (IPLAC.EQ.1) THEN
  2114. INMAT=INMAT+1
  2115. CALL MODPLA(MOMODL,NMOD)
  2116. CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  2117. IF (IPLAC.EQ.0) THEN
  2118. * ISOTROPE
  2119. IPLAC=1
  2120. * ELSE
  2121. * INMAT=INMAT+1
  2122. ENDIF
  2123. CALL IDPLAS(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  2124. IF (IRET.EQ.0) GOTO 9999
  2125. ELSEIF (IPLAC.EQ.2) THEN
  2126. INMAT=INMAT+1
  2127. CALL MODFLU(MOMODL,NMOD)
  2128. CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  2129. IF (IPLAC.EQ.0) THEN
  2130. * NORTON
  2131. IPLAC=1
  2132. * ELSE
  2133. * INMAT=INMAT+1
  2134. ENDIF
  2135. CALL IDFLUA(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  2136. IF (IRET.EQ.0) GOTO 9999
  2137. ELSEIF (IPLAC.EQ.3) THEN
  2138. INMAT=INMAT+1
  2139. CALL MODVIS(MOMODL,NMOD)
  2140. CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  2141. IF (IPLAC.EQ.0) THEN
  2142. * ONERA
  2143. IPLAC=2
  2144. * ELSE
  2145. * INMAT=INMAT+1
  2146. ENDIF
  2147. CALL IDVISC(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  2148. IF (IRET.EQ.0) GOTO 9999
  2149. ELSEIF (IPLAC.EQ.4) THEN
  2150. INMAT=INMAT+1
  2151. CALL MODEND(MOMODL,NMOD)
  2152. CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(INMAT))
  2153. IF (IPLAC.EQ.0) THEN
  2154. * MAZARS
  2155. IPLAC=1
  2156. * ELSE
  2157. * INMAT=INMAT+1
  2158. ENDIF
  2159. CALL IDENDO(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  2160. IF (IRET.EQ.0) GOTO 9999
  2161. ELSE
  2162. * GOTO 9999
  2163. ENDIF
  2164. *
  2165. * En cas de creation de materiaux combinant plusieurs materiaux
  2166. * deja existant Ex ELASTIQUE ISOTROPE PLASTIQUE PARFAIT FLUAGE N
  2167. * GOTO 40
  2168. *
  2169. * GOTO 9999
  2170. ENDIF
  2171. GOTO 9999
  2172. ENDIF
  2173.  
  2174. *--------------------------------------------------------------------
  2175. * CAS DE LA FORMULATION CONTACT
  2176. *--------------------------------------------------------------------
  2177. CALL PLACE(FORMOD,NFOR,ICONT,'CONTACT')
  2178. IF (ICONT.NE.0) THEN
  2179. CALL MODFRO(MOMODL,NMOD)
  2180. iplla =0
  2181. do iou = 1,NMAT
  2182. CALL PLACE(MOMODL(1),5,IPLAC,MATMOD(iou))
  2183. iplla = max(iplac,iplla)
  2184. enddo
  2185. IF (iplla.EQ.0) THEN
  2186. IRET = 0
  2187. GOTO 9999
  2188. ENDIF
  2189.  
  2190. IF ((IPLLA.EQ.1) .OR. (IPLLA.EQ.2) .OR. (IPLLA.EQ.3)) THEN
  2191. * CONTACT SIMPLE
  2192. * --------------
  2193. INMAT=1
  2194. JGOBL=0
  2195.  
  2196. JGFAC=2
  2197. TABFAC(1)='JEU'
  2198. TABFAC(2)='ADHE'
  2199.  
  2200. ELSEIF (IPLLA .EQ. 4) THEN
  2201. * FROTTEMENT DE COULOMB
  2202. * ---------------------
  2203. INMAT=1
  2204. JGOBL=1
  2205. TABOBL(1)='MU '
  2206.  
  2207. JGFAC=3
  2208. TABFAC(1)='COHE'
  2209. TABFAC(2)='ADHE'
  2210. TABFAC(3)='JEU'
  2211.  
  2212. ELSEIF (IPLLA .EQ. 5) THEN
  2213. * FROTTEMENTS DE CABLES
  2214. * ---------------------
  2215. INMAT=1
  2216. JGOBL=2
  2217. TABOBL(1)='FF '
  2218. TABOBL(2)='PHIF'
  2219.  
  2220. JGFAC=0
  2221.  
  2222. ELSE
  2223. IRET=0
  2224. CALL ERREUR (261)
  2225. RETURN
  2226. ENDIF
  2227.  
  2228. GOTO 9999
  2229. ENDIF
  2230.  
  2231. *--------------------------------------------------------------------
  2232. * CAS DE LA FORMULATION CONTRAINTE
  2233. *--------------------------------------------------------------------
  2234. CALL PLACE(FORMOD,NFOR,ICNTR,'CONTRAINTE')
  2235. IF (ICNTR.NE.0) THEN
  2236. CALL MODCLI(MOMODL,NMOD)
  2237. iplla =0
  2238. do iou = 1,NMAT
  2239. CALL PLACE(MOMODL(1),NMOD,IPLAC,MATMOD(iou))
  2240. iplla = max(iplac,iplla)
  2241. enddo
  2242. IF (iplla.EQ.0) THEN
  2243. IRET = 0
  2244. GOTO 9999
  2245. ENDIF
  2246.  
  2247.  
  2248. IF (IPLLA .EQ. 1) THEN
  2249. * ROTATION
  2250. * ---------------------
  2251. INMAT=1
  2252. JGOBL=1
  2253. TABOBL(1)='ANGL'
  2254. ELSEIF (IPLLA .EQ. 2) THEN
  2255. * DEPLACEMENT
  2256. * ---------------------
  2257. INMAT=1
  2258. JGOBL=1
  2259. TABOBL(1)='AMPL'
  2260. ENDIF
  2261. ENDIF
  2262. *--------------------------------------------------------------------
  2263. * CAS DE LA FORMULATION MAGNETODYNAMIQUE
  2264. *--------------------------------------------------------------------
  2265. CALL PLACE(FORMOD,NFOR,ICONV,'MAGNETODYNAMIQUE')
  2266. IF (ICONV.NE.0) THEN
  2267. *
  2268. * FORMULATION EN COQUES
  2269. IF(MFR.EQ.3) THEN
  2270. CALL PLACE(MATMOD,NMAT,ISOT,'ISOTROPE')
  2271. IF(ISOT.NE.0) THEN
  2272. JGOBL=3
  2273.  
  2274. TABOBL(1)='ETA'
  2275. TABOBL(2)='PERM'
  2276. TABOBL(3)='EPAI'
  2277. ELSE
  2278. CALL PLACE(MATMOD,NMAT,IORTH,'ORTHOTROPE')
  2279. IF(IORTH.NE.0) THEN
  2280. JGOBL=4
  2281.  
  2282. TABOBL(1)='ETA1'
  2283. TABOBL(2)='ETA2'
  2284. TABOBL(3)='PERM'
  2285. TABOBL(4)='EPAI'
  2286. ENDIF
  2287. ENDIF
  2288. ENDIF
  2289. *
  2290. GOTO 9999
  2291. ENDIF
  2292.  
  2293. *--------------------------------------------------------------------
  2294. * CAS DE LA FORMULATION FISSURE
  2295. *--------------------------------------------------------------------
  2296. CALL PLACE(FORMOD,NFOR,ICONV,'FISSURE')
  2297. IF (ICONV.NE.0) THEN
  2298. *
  2299. * si POISEU_BLASIUS ou POISEU_COLEBROOK ou par defaut
  2300. JGOBL=1
  2301.  
  2302. TABOBL(1)='RUGO'
  2303. * si FROTTEMENT1 ou FROTTEMENT2
  2304. CALL PLACE(MATMOD,NMAT,IFT1,'FROTTEMENT1')
  2305. CALL PLACE(MATMOD,NMAT,IFT2,'FROTTEMENT2')
  2306. IF(IFT1.NE.0.OR.IFT2.NE.0) THEN
  2307. JGOBL=7
  2308.  
  2309. TABOBL(2)='REC'
  2310. TABOBL(3)='FK'
  2311. TABOBL(4)='FA'
  2312. TABOBL(5)='FB'
  2313. TABOBL(6)='FC'
  2314. TABOBL(7)='FD'
  2315. ELSE
  2316. * si FROTTEMENT3 ou FROTTEMENT4
  2317. CALL PLACE(MATMOD,NMAT,IFT1,'FROTTEMENT3')
  2318. CALL PLACE(MATMOD,NMAT,IFT2,'FROTTEMENT4')
  2319. IF(IFT1.NE.0.OR.IFT2.NE.0) THEN
  2320. JGOBL=2
  2321.  
  2322. TABOBL(2)='FK'
  2323. ENDIF
  2324. ENDIF
  2325. GOTO 9999
  2326. ENDIF
  2327.  
  2328. *--------------------------------------------------------------------
  2329. * CAS DE LA FORMULATION NAVIER_STOKES
  2330. *--------------------------------------------------------------------
  2331. CALL PLACE(FORMOD,NFOR,INAST,'NAVIER_STOKES')
  2332. IF (INAST.GT.0) THEN
  2333. iplac = imatee
  2334. if (iplac.eq.0) then
  2335. iret = 0
  2336. goto 9999
  2337. endif
  2338. *
  2339. JGOBL = 0
  2340. JGFAC = 0
  2341.  
  2342. if (iplac.eq.4) then
  2343. *NLIN
  2344. JGOBL0 = JGOBL
  2345. JGOBL = JGOBL0+2
  2346. TABOBL(JGOBL0 + 1) ='RHO '
  2347. TABOBL(JGOBL0 + 2) ='MU '
  2348. JGOBL0 = JGFAC
  2349. JGFAC = JGOBL0+1
  2350. TABFAC(JGOBL0 + 1) ='FREQ'
  2351. endif
  2352. GOTO 9999
  2353. ENDIF
  2354.  
  2355. *--------------------------------------------------------------------
  2356. * CAS DE LA FORMULATION MELANGE
  2357. *--------------------------------------------------------------------
  2358. CALL PLACE(FORMOD,NFOR,ICONV,'MELANGE')
  2359. IF (ICONV.NE.0) THEN
  2360. CALL MODMEL(MOMODL,NMOD)
  2361. CALL PLACE(MOMODL,NMOD,IPLAC,MATMOD(1))
  2362. if (iplac.eq.0) then
  2363. goto 9999
  2364. endif
  2365. INMAT=1
  2366. *
  2367. * a priori elements massifs ou coques : pas de verif
  2368. *
  2369. IF (IPLAC.NE.3.AND.IPLAC.NE.4) THEN
  2370. JGOBL = 0
  2371. JGFAC = 0
  2372. CALL IDMETA(MFR,TABOBL,TABFAC,JGOBL,JGFAC,ITA,IPLAC,IRET)
  2373. IF (IRET.EQ.0) GOTO 9999
  2374. ELSE
  2375. ENDIF
  2376. GOTO 9999
  2377. ENDIF
  2378.  
  2379. *--------------------------------------------------------------------
  2380. * CAS DE LA FORMULATION LIAISON
  2381. *--------------------------------------------------------------------
  2382. CALL PLACE(FORMOD,NFOR,ICONV,'LIAISON')
  2383. IF (ICONV.NE.0) THEN
  2384. iplac = imatee
  2385. if (iplac.eq.0) then
  2386. iret = 0
  2387. goto 9999
  2388. endif
  2389. *
  2390. JGOBL = 0
  2391. JGFAC = 0
  2392.  
  2393. * 'SORT' facultatif dans tous les cas
  2394. IF (iplac.EQ.1) THEN
  2395. *PO_PL_FL
  2396. JGOBL0 = JGOBL
  2397. JGOBL = JGOBL0+7
  2398. TABOBL(JGOBL0 + 1) ='NORM'
  2399. TABOBL(JGOBL0 + 2) ='INER'
  2400. TABOBL(JGOBL0 + 3) ='CONV'
  2401. TABOBL(JGOBL0 + 4) ='VISC'
  2402. TABOBL(JGOBL0 + 5) ='PELO'
  2403. TABOBL(JGOBL0 + 6) ='PRAP'
  2404. TABOBL(JGOBL0 + 7) ='JFLU'
  2405.  
  2406. ELSEIF (iplac.eq.2) then
  2407. *PO_PL_FR
  2408. JGOBL0 = JGOBL
  2409. JGOBL = JGOBL0+7
  2410. TABOBL(JGOBL0 + 1) ='NORM'
  2411. TABOBL(JGOBL0 + 2) ='RAID'
  2412. TABOBL(JGOBL0 + 3) ='JEU'
  2413. TABOBL(JGOBL0 + 4) ='GLIS'
  2414. TABOBL(JGOBL0 + 5) ='ADHE'
  2415. TABOBL(JGOBL0 + 6) ='RTAN'
  2416. TABOBL(JGOBL0 + 7) ='ATAN'
  2417.  
  2418. JGOBL0 = JGFAC
  2419. JGFAC = JGOBL0 +2
  2420. TABFAC(JGOBL0+1) = 'AMOR'
  2421. TABFAC(JGOBL0+2) = 'LOIC'
  2422.  
  2423. ELSEIF (iplac.eq.3) then
  2424. *PO_PL
  2425. JGOBL0 = JGOBL
  2426. JGOBL= JGOBL0+3
  2427. TABOBL(JGOBL0 + 1) ='NORM'
  2428. TABOBL(JGOBL0 + 2) ='RAID'
  2429. TABOBL(JGOBL0 + 3) ='JEU'
  2430.  
  2431. JGOBL0 = JGFAC
  2432. JGFAC = JGOBL0+4
  2433. TABFAC(JGOBL0+1) = 'LOIC'
  2434. TABFAC(JGOBL0+2) = 'PERM'
  2435. TABFAC(JGOBL0 + 3) ='SPLA'
  2436. TABFAC(JGOBL0 + 4) ='AMOR'
  2437.  
  2438. ELSEIF (iplac.eq.4) then
  2439. *PO_PO_FR
  2440. JGOBL0 = JGOBL
  2441. JGOBL= JGOBL0 + 8
  2442. TABOBL(JGOBL0 + 1) ='NORM'
  2443. TABOBL(JGOBL0 + 2) ='RAID'
  2444. TABOBL(JGOBL0 + 3) ='JEU'
  2445. TABOBL(JGOBL0 + 4) ='POIB'
  2446. TABOBL(JGOBL0 + 5) ='ADHE'
  2447. TABOBL(JGOBL0 + 6) ='RTAN'
  2448. TABOBL(JGOBL0 + 7) ='ATAN'
  2449. TABOBL(JGOBL0 + 8) ='GLIS'
  2450.  
  2451. JGOBL0 = JGFAC
  2452. JGFAC = JGOBL0 + 3
  2453. TABFAC(JGOBL0+1) = 'AMOR'
  2454. TABFAC(JGOBL0+2) = 'LOIC'
  2455. TABFAC(JGOBL0+3) = 'MODE'
  2456.  
  2457. ELSEIF (iplac.eq.5) then
  2458. *PO_PO_DP
  2459. JGOBL0 = JGOBL
  2460. JGOBL= JGOBL0 + 6
  2461. TABOBL(JGOBL0 + 1) ='NORM'
  2462. TABOBL(JGOBL0 + 2) ='ECRO'
  2463. TABOBL(JGOBL0 + 3) ='JEU'
  2464. TABOBL(JGOBL0 + 4) ='POIB'
  2465. TABOBL(JGOBL0 + 5) ='PERM'
  2466. TABOBL(JGOBL0 + 6) ='LOIC'
  2467.  
  2468. JGOBL0 = JGFAC
  2469. JGFAC = JGOBL0 + 1
  2470. TABFAC(JGOBL0+1) = 'AMOR'
  2471.  
  2472. ELSEIF (iplac.eq.6) then
  2473. *PO_PO_RP
  2474. JGOBL0 = JGOBL
  2475. JGOBL= JGOBL0 + 6
  2476. TABOBL(JGOBL0 + 1) ='AXRO'
  2477. TABOBL(JGOBL0 + 2) ='ECRO'
  2478. TABOBL(JGOBL0 + 3) ='JEU'
  2479. TABOBL(JGOBL0 + 4) ='POIB'
  2480. TABOBL(JGOBL0 + 5) ='PERM'
  2481. TABOBL(JGOBL0 + 6) ='LOIC'
  2482.  
  2483. JGOBL0 = JGFAC
  2484. JGFAC = JGOBL0 + 2
  2485. TABFAC(JGOBL0+1) = 'AMOR'
  2486. TABFAC(JGOBL0+2) = 'ELAS'
  2487.  
  2488. ELSEIF (iplac.eq.7) then
  2489. *PO_PO
  2490. JGOBL0 = JGOBL
  2491. JGOBL= JGOBL0 + 5
  2492. TABOBL(JGOBL0 + 1) ='NORM'
  2493. TABOBL(JGOBL0 + 2) ='RAID'
  2494. TABOBL(JGOBL0 + 3) ='JEU'
  2495. TABOBL(JGOBL0 + 4) ='POIB'
  2496. TABOBL(JGOBL0 + 5) ='PERM'
  2497.  
  2498. JGOBL0 = JGFAC
  2499. JGFAC = JGOBL0 + 2
  2500. TABFAC(JGOBL0+1) = 'AMOR'
  2501. TABFAC(JGOBL0+2) = 'LOIC'
  2502.  
  2503. ELSEIF (iplac.eq.8) then
  2504. *PO_CE_MO
  2505. JGOBL0 = JGOBL
  2506. JGOBL= JGOBL0 + 8
  2507. TABOBL(JGOBL0 + 1) ='NORM'
  2508. TABOBL(JGOBL0 + 2) ='RAID'
  2509. TABOBL(JGOBL0 + 3) ='PCER'
  2510. TABOBL(JGOBL0 + 4) ='RAYO'
  2511. TABOBL(JGOBL0 + 5) ='GLIS'
  2512. TABOBL(JGOBL0 + 6) ='ADHE'
  2513. TABOBL(JGOBL0 + 7) ='RTAN'
  2514. TABOBL(JGOBL0 + 8) ='ATAN'
  2515.  
  2516. JGOBL0 = JGFAC
  2517. JGFAC= JGOBL0 + 2
  2518. TABFAC(JGOBL0+1) = 'CINT'
  2519. TABFAC(JGOBL0+2) = 'AMOR'
  2520.  
  2521. ELSEIF (iplac.eq.9) then
  2522. *PO_CE_FR
  2523. JGOBL0 = JGOBL
  2524. JGOBL= JGOBL0 + 8
  2525. TABOBL(JGOBL0 + 1) ='NORM'
  2526. TABOBL(JGOBL0 + 2) ='RAID'
  2527. TABOBL(JGOBL0 + 3) ='EXCE'
  2528. TABOBL(JGOBL0 + 4) ='RAYO'
  2529. TABOBL(JGOBL0 + 5) ='GLIS'
  2530. TABOBL(JGOBL0 + 6) ='ADHE'
  2531. TABOBL(JGOBL0 + 7) ='RTAN'
  2532. TABOBL(JGOBL0 + 8) ='ATAN'
  2533.  
  2534. JGOBL0 = JGFAC
  2535. JGFAC = JGOBL0 + 2
  2536. TABFAC(JGOBL0+1) = 'CINT'
  2537. TABFAC(JGOBL0+2) = 'AMOR'
  2538.  
  2539. ELSEIF (iplac.eq.10) then
  2540. *PO_CE
  2541. JGOBL0 = JGOBL
  2542. JGOBL= JGOBL0 + 4
  2543. TABOBL(JGOBL0 + 1) ='NORM'
  2544. TABOBL(JGOBL0 + 2) ='RAID'
  2545. TABOBL(JGOBL0 + 3) ='EXCE'
  2546. TABOBL(JGOBL0 + 4) ='RAYO'
  2547.  
  2548. JGOBL0 = JGFAC
  2549. JGFAC = JGOBL0 + 1
  2550. TABFAC(JGOBL0+1) = 'AMOR'
  2551.  
  2552. ELSEIF (iplac.eq.11) then
  2553. *CE_PL_FR
  2554. JGOBL0 = JGOBL
  2555. JGOBL= JGOBL0 + 8
  2556. TABOBL(JGOBL0 + 1) ='NORM'
  2557. TABOBL(JGOBL0 + 2) ='RAID'
  2558. TABOBL(JGOBL0 + 3) ='JEU'
  2559. TABOBL(JGOBL0 + 4) ='RAYS'
  2560. TABOBL(JGOBL0 + 5) ='GLIS'
  2561. TABOBL(JGOBL0 + 6) ='ADHE'
  2562. TABOBL(JGOBL0 + 7) ='RTAN'
  2563. TABOBL(JGOBL0 + 8) ='ATAN'
  2564.  
  2565. JGOBL0 = JGFAC
  2566. JGFAC = JGOBL0+1
  2567. TABFAC(JGOBL0+1) = 'AMOR'
  2568.  
  2569. ELSEIF (iplac.eq.12) then
  2570. *CE_CE_FR
  2571. JGOBL0 = JGOBL
  2572. JGOBL= JGOBL0 + 9
  2573. TABOBL(JGOBL0 + 1) ='NORM'
  2574. TABOBL(JGOBL0 + 2) ='RAID'
  2575. TABOBL(JGOBL0 + 3) ='EXCE'
  2576. TABOBL(JGOBL0 + 4) ='RAYS'
  2577. TABOBL(JGOBL0 + 5) ='GLIS'
  2578. TABOBL(JGOBL0 + 6) ='ADHE'
  2579. TABOBL(JGOBL0 + 7) ='RTAN'
  2580. TABOBL(JGOBL0 + 8) ='ATAN'
  2581. TABOBL(JGOBL0 + 9) ='RAYB'
  2582.  
  2583. JGOBL0 = JGFAC
  2584. JGFAC = JGOBL0 + 2
  2585. TABFAC(JGOBL0+1) = 'AMOR'
  2586. TABFAC(JGOBL0+2) = 'CINT'
  2587.  
  2588. ELSEIF (iplac.eq.13.or.iplac.eq.14) then
  2589. *PR_PR_IN ou PR_PR_EX
  2590. JGOBL0 = JGOBL
  2591. JGOBL= JGOBL0 + 5
  2592. TABOBL(JGOBL0 + 1) ='NORM'
  2593. TABOBL(JGOBL0 + 2) ='RAID'
  2594. TABOBL(JGOBL0 + 3) ='PFIX'
  2595. TABOBL(JGOBL0 + 4) ='PMOB'
  2596. TABOBL(JGOBL0 + 5) ='ERAI'
  2597.  
  2598. ELSEIF (iplac.eq.15) then
  2599. *LI_LI_FR
  2600. JGOBL0 = JGOBL
  2601. JGOBL= JGOBL0 + 8
  2602. TABOBL(JGOBL0 + 1) ='NORM'
  2603. TABOBL(JGOBL0 + 2) ='LIMA'
  2604. TABOBL(JGOBL0 + 3) ='LIES'
  2605. TABOBL(JGOBL0 + 4) ='RAID'
  2606. TABOBL(JGOBL0 + 5) ='GLIS'
  2607. TABOBL(JGOBL0 + 6) ='ADHE'
  2608. TABOBL(JGOBL0 + 7) ='RTAN'
  2609. TABOBL(JGOBL0 + 8) ='ATAN'
  2610.  
  2611. JGOBL0 = JGFAC
  2612. JGFAC = JGOBL0 + 4
  2613. TABFAC(JGOBL0+1) = 'JEU'
  2614. TABFAC(JGOBL0+2) = 'AMOR'
  2615. TABFAC(JGOBL0+3) = 'RECH'
  2616. TABFAC(JGOBL0+4) = 'SYME'
  2617.  
  2618. ELSEIF (iplac.eq.16) then
  2619. *LI_CE_FR
  2620. JGOBL0 = JGOBL
  2621. JGOBL= JGOBL0 + 8
  2622. TABOBL(JGOBL0 + 1) ='NORM'
  2623. TABOBL(JGOBL0 + 2) ='LIMA'
  2624. TABOBL(JGOBL0 + 3) ='LIES'
  2625. TABOBL(JGOBL0 + 4) ='RAID'
  2626. TABOBL(JGOBL0 + 5) ='GLIS'
  2627. TABOBL(JGOBL0 + 6) ='ADHE'
  2628. TABOBL(JGOBL0 + 7) ='RTAN'
  2629. TABOBL(JGOBL0 + 8) ='ATAN'
  2630.  
  2631. JGOBL0 = JGFAC
  2632. JGFAC = JGOBL0 + 5
  2633. TABFAC(JGOBL0+1) = 'AMOR'
  2634. TABFAC(JGOBL0+2) = 'RECH'
  2635. TABFAC(JGOBL0+3) = 'RAYO'
  2636. TABFAC(JGOBL0+4) = 'ACTN'
  2637. TABFAC(JGOBL0+5) = 'INVE'
  2638.  
  2639. ELSEIF (iplac.eq.17) then
  2640. *PA_FL_RO
  2641. JGOBL0 = JGOBL
  2642. JGOBL= JGOBL0+10
  2643. TABOBL(JGOBL0 + 1) ='LONG'
  2644. TABOBL(JGOBL0 + 2) ='RAYO'
  2645. TABOBL(JGOBL0 + 3) ='VISC'
  2646. TABOBL(JGOBL0 + 4) ='RHOF'
  2647. TABOBL(JGOBL0 + 5) ='PADM'
  2648. TABOBL(JGOBL0 + 6) ='VROT'
  2649. TABOBL(JGOBL0 + 7) ='EPSI'
  2650. TABOBL(JGOBL0 + 8) ='PHII'
  2651. TABOBL(JGOBL0 + 9) ='AFFI'
  2652. TABOBL(JGOBL0 + 10)='TLOB'
  2653.  
  2654. JGOBL0 = JGFAC
  2655. JGFAC = JGOBL0 + 1
  2656. TABFAC(JGOBL0+1) = 'AMOR'
  2657.  
  2658. ELSEIF (iplac.eq.23) then
  2659. *NEWMARK MODAL
  2660. JGOBL0 = JGOBL
  2661. JGOBL= JGOBL0+3
  2662. TABOBL(JGOBL0 + 1) ='JEU'
  2663. TABOBL(JGOBL0 + 2) ='MASS'
  2664. TABOBL(JGOBL0 + 3) ='FREQ'
  2665.  
  2666. JGOBL0 = JGFAC
  2667. JGFAC = JGOBL0 + 3
  2668. TABFAC(JGOBL0+1) = 'EXCE'
  2669. TABFAC(JGOBL0+2) = 'FROT'
  2670. TABFAC(JGOBL0+3) = 'MOFR'
  2671. ENDIF
  2672.  
  2673. * 'SORT' facultatif dans tous les cas
  2674. JGOBL0 = JGFAC
  2675. JGFAC= JGOBL0+1
  2676. TABFAC(JGOBL0+1) = 'SORT'
  2677. GOTO 9999
  2678. ENDIF
  2679.  
  2680. *--------------------------------------------------------------------
  2681. * CAS DE LA FORMULATION ELECTROSTATIQUE
  2682. *--------------------------------------------------------------------
  2683. CALL PLACE(FORMOD,NFOR,IELEC,'ELECTROSTATIQUE')
  2684. IF (IELEC.NE.0) THEN
  2685. C -- Permittivite isotrope
  2686. IF (IMATEE.EQ.1) THEN
  2687. C* IF (CMATEE.EQ.'ISOTROPE') THEN
  2688. C* IF (MATMOD(1).EQ.'ISOTROPE ') THEN
  2689. JGOBL = 1
  2690. TABOBL(1)='PEL '
  2691. C -- Permittivite orthotrope
  2692. ELSEIF (IMATEE.EQ.2) THEN
  2693. C* ELSEIF (CMATEE.EQ.'ORTHOTRO') THEN
  2694. C* ELSEIF (MATMOD(1).EQ.'ORTHOTROPE ') THEN
  2695. C ---- Elements massifs bidimensionnels PLAN et AXISYMETRIQUE
  2696. IF (IDIM.EQ.2) THEN
  2697. IF (IFOMOD.NE.1) THEN
  2698. JGOBL = 4
  2699. TABOBL(1) = 'PE1 '
  2700. TABOBL(2) = 'PE2 '
  2701. TABOBL(3) = 'V1X '
  2702. TABOBL(4) = 'V1Y '
  2703.  
  2704. C ---- Elements massifs bidimensionnels FOURIER
  2705. ELSE
  2706. JGOBL = 5
  2707. TABOBL(1) = 'PE1 '
  2708. TABOBL(2) = 'PE2 '
  2709. TABOBL(3) = 'PE3 '
  2710. TABOBL(4) = 'V1X '
  2711. TABOBL(5) = 'V1Y '
  2712. ENDIF
  2713. C ---- Elements massifs TRIDimensionnels
  2714. ELSEIF (IDIM.EQ.3) THEN
  2715. JGOBL = 9
  2716. TABOBL(1) = 'PE1 '
  2717. TABOBL(2) = 'PE2 '
  2718. TABOBL(3) = 'PE3 '
  2719. TABOBL(4) = 'V1X '
  2720. TABOBL(5) = 'V1Y '
  2721. TABOBL(6) = 'V1Z '
  2722. TABOBL(7) = 'V2X '
  2723. TABOBL(8) = 'V2Y '
  2724. TABOBL(9) = 'V2Z '
  2725. ENDIF
  2726.  
  2727. C -- Permittivite anisotrope
  2728. ELSEIF (IMATEE.EQ.3) THEN
  2729. C* ELSEIF (CMATEE.EQ.'ANISOTRO') THEN
  2730. C* ELSEIF (MATMOD(1).EQ.'ANISOTROPE ') THEN
  2731. C ---- Elements massifs bidimensionnels PLAN et AXISYMETRIQUE
  2732. IF (IDIM.EQ.2) THEN
  2733. IF (IFOMOD.NE.1) THEN
  2734. JGOBL = 5
  2735. TABOBL(1) = 'PE11 '
  2736. TABOBL(2) = 'PE22 '
  2737. TABOBL(3) = 'PE21 '
  2738. TABOBL(4) = 'V1X '
  2739. TABOBL(5) = 'V1Y '
  2740.  
  2741. C ---- Elements massifs bidimensionnels FOURIER
  2742. ELSE
  2743. JGOBL = 6
  2744. TABOBL(1) = 'PE11 '
  2745. TABOBL(2) = 'PE22 '
  2746. TABOBL(3) = 'PE21 '
  2747. TABOBL(4) = 'PE33 '
  2748. TABOBL(5) = 'V1X '
  2749. TABOBL(6) = 'V1Y '
  2750. ENDIF
  2751.  
  2752. C ---- Elements massifs TRIDimensionnels
  2753. ELSEIF (IDIM.EQ.3) THEN
  2754. JGOBL = 12
  2755. TABOBL( 1) = 'PE11 '
  2756. TABOBL( 2) = 'PE22 '
  2757. TABOBL( 3) = 'PE33 '
  2758. TABOBL( 4) = 'PE21 '
  2759. TABOBL( 5) = 'PE31 '
  2760. TABOBL( 6) = 'PE32 '
  2761. TABOBL( 7) = 'V1X '
  2762. TABOBL( 8) = 'V1Y '
  2763. TABOBL( 9) = 'V1Z '
  2764. TABOBL(10) = 'V2X '
  2765. TABOBL(11) = 'V2Y '
  2766. TABOBL(12) = 'V2Z '
  2767. ENDIF
  2768.  
  2769. ELSE
  2770. IRET = 0
  2771. CALL ERREUR(5)
  2772. ENDIF
  2773. GOTO 9999
  2774. ENDIF
  2775.  
  2776. *-----------------------------------------------------------------------
  2777. * CAS DE LA FORMULATION 'DIFFUSION'
  2778. *-----------------------------------------------------------------------
  2779. CALL PLACE(FORMOD,NFOR,IDIFF,'DIFFUSION')
  2780. IF (IDIFF.NE.0) THEN
  2781. JGFAC = 1
  2782. TABFAC(1)='CDIF'
  2783.  
  2784. C- Diffusion isotrope
  2785. IF (IMATEE .EQ. 1) THEN
  2786. JGOBL = 1
  2787. C*8 TABOBL(1) = 'KDIFF '
  2788. TABOBL(1) = 'KD '
  2789.  
  2790. C- Diffusion orthotrope
  2791. ELSEIF (IMATEE .EQ. 2) THEN
  2792. C --- Elements massifs bidimensionnels PLAN et AXISYMETRIQUE
  2793. IF (IDIM.EQ.2) THEN
  2794. IF (IFOMOD.NE.1) THEN
  2795. JGOBL = 4
  2796. TABOBL(1) = 'KD1 '
  2797. TABOBL(2) = 'KD2 '
  2798. TABOBL(3) = 'V1X '
  2799. TABOBL(4) = 'V1Y '
  2800.  
  2801. C --- Elements massifs bidimensionnels FOURIER
  2802. ELSE
  2803. JGOBL = 5
  2804. TABOBL(1) = 'KD1 '
  2805. TABOBL(2) = 'KD2 '
  2806. TABOBL(3) = 'KD3 '
  2807. TABOBL(4) = 'V1X '
  2808. TABOBL(5) = 'V1Y '
  2809. ENDIF
  2810.  
  2811. C --- Elements massifs TRIDimensionnels
  2812. ELSEIF (IDIM.EQ.3) THEN
  2813. JGOBL = 9
  2814. TABOBL(1) = 'KD1 '
  2815. TABOBL(2) = 'KD2 '
  2816. TABOBL(3) = 'KD3 '
  2817. TABOBL(4) = 'V1X '
  2818. TABOBL(5) = 'V1Y '
  2819. TABOBL(6) = 'V1Z '
  2820. TABOBL(7) = 'V2X '
  2821. TABOBL(8) = 'V2Y '
  2822. TABOBL(9) = 'V2Z '
  2823. ENDIF
  2824.  
  2825. C - Diffusion anisotrope
  2826. ELSEIF (IMATEE.EQ.3) THEN
  2827. C --- Elements massifs bidimensionnels PLAN et AXISYMETRIQUE
  2828. IF (IDIM.EQ.2) THEN
  2829. IF (IFOMOD.NE.1) THEN
  2830. JGOBL = 5
  2831. TABOBL(1) = 'KD11 '
  2832. TABOBL(2) = 'KD22 '
  2833. TABOBL(3) = 'KD21 '
  2834. TABOBL(4) = 'V1X '
  2835. TABOBL(5) = 'V1Y '
  2836.  
  2837. C --- Elements massifs bidimensionnels FOURIER
  2838. ELSE
  2839. JGOBL = 6
  2840. TABOBL(1) = 'KD11 '
  2841. TABOBL(2) = 'KD22 '
  2842. TABOBL(3) = 'KD21 '
  2843. TABOBL(4) = 'KD33 '
  2844. TABOBL(5) = 'V1X '
  2845. TABOBL(6) = 'V1Y '
  2846. ENDIF
  2847.  
  2848. C --- Elements massifs TRIDimensionnels
  2849. ELSEIF (IDIM.EQ.3) THEN
  2850. JGOBL = 12
  2851. TABOBL( 1) = 'KD11 '
  2852. TABOBL( 2) = 'KD22 '
  2853. TABOBL( 3) = 'KD33 '
  2854. TABOBL( 4) = 'KD21 '
  2855. TABOBL( 5) = 'KD31 '
  2856. TABOBL( 6) = 'KD32 '
  2857. TABOBL( 7) = 'V1X '
  2858. TABOBL( 8) = 'V1Y '
  2859. TABOBL( 9) = 'V1Z '
  2860. TABOBL(10) = 'V2X '
  2861. TABOBL(11) = 'V2Y '
  2862. TABOBL(12) = 'V2Z '
  2863. ENDIF
  2864. ENDIF
  2865.  
  2866. JGOBL0 = JGOBL
  2867. JGFAC0 = JGFAC
  2868.  
  2869. * Loi FICK :
  2870. * ------------
  2871. IF (INATUU.EQ.0) THEN
  2872. * Rien a ajouter
  2873.  
  2874. * Loi SORET :
  2875. * -------------
  2876. ELSEIF (INATUU.EQ.1) THEN
  2877. C -- Diffusion isotrope
  2878. IF (IMATEE.EQ.1) THEN
  2879. JGOBL = JGOBL+1
  2880. C*8 TABOBL(JGOBL0+1) = 'KSORET '
  2881. TABOBL(JGOBL0+1) = 'KS '
  2882.  
  2883. C -- Diffusion orthotrope
  2884. ELSEIF (IMATEE.EQ.2) THEN
  2885. C ---- Elements massifs bidimensionnels PLAN et AXISYMETRIQUE
  2886. IF (IFOMOD.LE.0) THEN
  2887. JGOBL = JGOBL + 2
  2888. C*8 TABOBL(JGOBL0+1) = 'KSORET1 '
  2889. C*8 TABOBL(JGOBL0+2) = 'KSORET2 '
  2890. TABOBL(JGOBL0+1) = 'KS1 '
  2891. TABOBL(JGOBL0+2) = 'KS2 '
  2892.  
  2893. C ---- Elements massifs bidimensionnels FOURIER
  2894. C ---- Elements massifs TRIDimensionnels
  2895. ELSE
  2896. JGOBL = JGOBL + 3
  2897. C*8 TABOBL(JGOBL0+1) = 'KSORET1 '
  2898. C*8 TABOBL(JGOBL0+2) = 'KSORET2 '
  2899. C*8 TABOBL(JGOBL0+3) = 'KSORET3 '
  2900. TABOBL(JGOBL0+1) = 'KS1 '
  2901. TABOBL(JGOBL0+2) = 'KS2 '
  2902. TABOBL(JGOBL0+3) = 'KS3 '
  2903. ENDIF
  2904.  
  2905. C -- Diffusion anisotrope
  2906. ELSEIF (IMATEE.EQ.3) THEN
  2907. C ---- Elements massifs bidimensionnels PLAN et AXISYMETRIQUE
  2908. IF (IFOMOD.LE.0) THEN
  2909. JGOBL = JGOBL + 3
  2910. C*8 TABOBL(JGOBL0+1) = 'KSORET11'
  2911. C*8 TABOBL(JGOBL0+2) = 'KSORET22'
  2912. C*8 TABOBL(JGOBL0+3) = 'KSORET21'
  2913. TABOBL(JGOBL0+1) = 'KS11 '
  2914. TABOBL(JGOBL0+2) = 'KS22 '
  2915. TABOBL(JGOBL0+3) = 'KS21 '
  2916.  
  2917. C ---- Elements massifs bidimensionnels FOURIER
  2918. ELSEIF (IFOMOD.EQ.1) THEN
  2919. JGOBL = JGOBL + 4
  2920. C*8 TABOBL(JGOBL0+1) = 'KSORET11'
  2921. C*8 TABOBL(JGOBL0+2) = 'KSORET22'
  2922. C*8 TABOBL(JGOBL0+3) = 'KSORET21'
  2923. C*8 TABOBL(JGOBL0+4) = 'KSORET33'
  2924. TABOBL(JGOBL0+1) = 'KS11 '
  2925. TABOBL(JGOBL0+2) = 'KS22 '
  2926. TABOBL(JGOBL0+3) = 'KS21 '
  2927. TABOBL(JGOBL0+4) = 'KS33 '
  2928.  
  2929. C ---- Elements massifs TRIDimensionnels
  2930. ELSEIF (IFOMOD.EQ.2) THEN
  2931. JGOBL = JGOBL + 6
  2932. C*8 TABOBL(JGOBL0+1) = 'KSORET11'
  2933. C*8 TABOBL(JGOBL0+2) = 'KSORET22'
  2934. C*8 TABOBL(JGOBL0+3) = 'KSORET33'
  2935. C*8 TABOBL(JGOBL0+4) = 'KSORET21'
  2936. C*8 TABOBL(JGOBL0+5) = 'KSORET31'
  2937. C*8 TABOBL(JGOBL0+6) = 'KSORET32'
  2938. TABOBL(JGOBL0+1) = 'KS11 '
  2939. TABOBL(JGOBL0+2) = 'KS22 '
  2940. TABOBL(JGOBL0+3) = 'KS33 '
  2941. TABOBL(JGOBL0+4) = 'KS21 '
  2942. TABOBL(JGOBL0+5) = 'KS31 '
  2943. TABOBL(JGOBL0+6) = 'KS32 '
  2944. ENDIF
  2945.  
  2946. ELSE
  2947. C -- Erreur ne devant pas survenir
  2948. CALL ERREUR(261)
  2949. RETURN
  2950. ENDIF
  2951.  
  2952. * Loi UTILISATEUR :
  2953. * -------------------
  2954. * ELSEIF (INATUU.EQ.2) THEN
  2955. * Par defaut pas de composantes autres que celles de l'utilisateur
  2956.  
  2957. ELSEIF (INATUU.EQ.5) THEN
  2958. * diffusion avec dissipation affine
  2959. JGOBL = JGOBL+2
  2960. TABOBL(JGOBL0+1)='DAF0'
  2961. TABOBL(JGOBL0+2)='DAF1'
  2962.  
  2963. ELSEIF (INATUU.EQ.6) THEN
  2964. * diffusion avec dissipation visqueuse
  2965. JGOBL = JGOBL+2
  2966. TABOBL(JGOBL0+1)='DVQ0'
  2967. TABOBL(JGOBL0+2)='DVQ1'
  2968.  
  2969. * Loi inconnue : ERREUR
  2970. * ----------------
  2971. ELSE
  2972. CALL ERREUR(261)
  2973. RETURN
  2974. ENDIF
  2975.  
  2976. * Composantes ADVECTION :
  2977. * -------------------
  2978. CALL PLACE(MATMOD,NMAT,IADVE,'ADVECTION')
  2979. if (iadve .ne. 0) then
  2980.  
  2981. C Cas des Tuyaux 1D (MFR=79)
  2982. C ===========================
  2983. IF (MFR .EQ. 79) THEN
  2984. JGOBL = JGOBL + 1
  2985. TABOBL(JGOBL0 + 1)='VITE'
  2986.  
  2987. JGFAC = JGFAC + 1
  2988. TABFAC(JGFAC0 + 1)='DL'
  2989.  
  2990. C Cas Diffusion Massif (MFR=73)
  2991. C ============================
  2992. ELSEIF (MFR .EQ. 73) THEN
  2993. C Bidimensionnel PLAN DPGE
  2994. C ===========================
  2995. IF (IFOUR.EQ.-3) THEN
  2996. JGOBL = JGOBL + 3
  2997. TABOBL(JGOBL0 + 1)='VITX'
  2998. TABOBL(JGOBL0 + 2)='VITY'
  2999. TABOBL(JGOBL0 + 3)='VITZ'
  3000.  
  3001. JGFAC = JGFAC + 3
  3002. TABFAC(JGFAC0 + 1)='DX'
  3003. TABFAC(JGFAC0 + 2)='DY'
  3004. TABFAC(JGFAC0 + 3)='DZ'
  3005.  
  3006. C Bidimensionnel PLAN (CP/DP)
  3007. C ===========================
  3008. ELSEIF (IFOUR.EQ.-2 .OR. IFOUR.EQ.-1) THEN
  3009. JGOBL = JGOBL + 2
  3010. TABOBL(JGOBL0 + 1)='VITX'
  3011. TABOBL(JGOBL0 + 2)='VITY'
  3012.  
  3013. JGFAC = JGFAC + 2
  3014. TABFAC(JGFAC0 + 1)='DX'
  3015. TABFAC(JGFAC0 + 2)='DY'
  3016.  
  3017. C Axisymetrie
  3018. C ===========================
  3019. ELSEIF (IFOUR .EQ. 0) THEN
  3020. JGOBL = JGOBL + 2
  3021. TABOBL(JGOBL0 + 1)='VITR'
  3022. TABOBL(JGOBL0 + 2)='VITZ'
  3023.  
  3024. JGFAC = JGFAC + 2
  3025. TABFAC(JGFAC0 + 1)='DR'
  3026. TABFAC(JGFAC0 + 2)='DZ'
  3027.  
  3028. C Fourier
  3029. C ===========================
  3030. ELSEIF (IFOUR .EQ. 1) THEN
  3031. CALL ERREUR(21)
  3032. RETURN
  3033. C JGOBL = JGOBL + 3
  3034. C TABOBL(JGOBL0 + 1)='VITR'
  3035. C TABOBL(JGOBL0 + 2)='VITZ'
  3036. C TABOBL(JGOBL0 + 3)='VITT'
  3037. C
  3038. C JGFAC = JGFAC + 3
  3039. C TABFAC(JGFAC0 + 1)='DR'
  3040. C TABFAC(JGFAC0 + 2)='DZ'
  3041. C TABFAC(JGFAC0 + 3)='DT'
  3042.  
  3043. C Tridimensionnel
  3044. C ===========================
  3045. ELSEIF (IFOUR .EQ. 2) THEN
  3046. JGOBL = JGOBL + 3
  3047. TABOBL(JGOBL0 + 1)='VITX'
  3048. TABOBL(JGOBL0 + 2)='VITY'
  3049. TABOBL(JGOBL0 + 3)='VITZ'
  3050.  
  3051. JGFAC = JGFAC + 3
  3052. TABFAC(JGFAC0 + 1)='DX'
  3053. TABFAC(JGFAC0 + 2)='DY'
  3054. TABFAC(JGFAC0 + 3)='DZ'
  3055.  
  3056. C Unidimensionnel (1D)
  3057. C ===========================
  3058. ELSEIF (IFOUR.GE.3 .AND. IFOUR.LE.15) THEN
  3059. IF (IFOUR.LE.6) THEN
  3060. JGOBL = JGOBL + 1
  3061. TABOBL(JGOBL0 + 1)='VITX'
  3062.  
  3063. JGFAC = JGFAC + 1
  3064. TABFAC(JGFAC0 + 1)='DX'
  3065.  
  3066. ELSEIF (IFOUR.EQ.7 .OR. IFOUR.EQ.8) THEN
  3067. JGOBL = JGOBL + 2
  3068. TABOBL(JGOBL0 + 1)='VITX'
  3069. TABOBL(JGOBL0 + 2)='VITY'
  3070.  
  3071. JGFAC = JGFAC + 2
  3072. TABFAC(JGFAC0 + 1)='DX'
  3073. TABFAC(JGFAC0 + 2)='DY'
  3074.  
  3075. ELSEIF (IFOUR.EQ.9 .OR. IFOUR.EQ.10) THEN
  3076. JGOBL = JGOBL + 2
  3077. TABOBL(JGOBL0 + 1)='VITX'
  3078. TABOBL(JGOBL0 + 2)='VITZ'
  3079.  
  3080. JGFAC = JGFAC + 2
  3081. TABFAC(JGFAC0 + 1)='DX'
  3082. TABFAC(JGFAC0 + 2)='DZ'
  3083.  
  3084. ELSEIF (IFOUR.EQ.11) THEN
  3085. JGOBL = JGOBL + 3
  3086. TABOBL(JGOBL0 + 1)='VITX'
  3087. TABOBL(JGOBL0 + 2)='VITY'
  3088. TABOBL(JGOBL0 + 3)='VITZ'
  3089.  
  3090. JGFAC = JGFAC + 3
  3091. TABFAC(JGFAC0 + 1)='DX'
  3092. TABFAC(JGFAC0 + 2)='DY'
  3093. TABFAC(JGFAC0 + 3)='DZ'
  3094.  
  3095. ELSEIF (IFOUR.EQ.12.OR.IFOUR.EQ.13.OR.IFOUR.EQ.15) THEN
  3096. JGOBL = JGOBL + 1
  3097. TABOBL(JGOBL0 + 1)='VITR'
  3098.  
  3099. JGFAC = JGFAC + 1
  3100. TABFAC(JGFAC0 + 1)='DR'
  3101.  
  3102. ELSEIF (IFOUR.EQ.14) THEN
  3103. JGOBL = JGOBL + 2
  3104. TABOBL(JGOBL0 + 1)='VITR'
  3105. TABOBL(JGOBL0 + 2)='VITZ'
  3106.  
  3107. JGFAC = JGFAC + 2
  3108. TABFAC(JGFAC0 + 1)='DR'
  3109. TABFAC(JGFAC0 + 2)='DZ'
  3110. ENDIF
  3111. ENDIF
  3112. ENDIF
  3113. ENDIF
  3114.  
  3115. * WRITE(6,*) 'TABOBL',JGOBL,(TABOBL(i),i=1,JGOBL)
  3116. * WRITE(6,*) 'TABFAC',JGFAC,(TABFAC(i),i=1,JGFAC)
  3117. GOTO 9999
  3118. ENDIF
  3119.  
  3120. *-----------------------------------------------------------------------
  3121. * CAS DE LA FORMULATION 'CHARGEMENT'
  3122. *-----------------------------------------------------------------------
  3123. CALL PLACE(FORMOD,NFOR,ICHAR,'CHARGEMENT')
  3124. IF (ICHAR.NE.0) THEN
  3125. IPLAC = IMATEE
  3126. IF (IPLAC.EQ.1) THEN
  3127. C RAJOUTER DIFFERENTIATION EN FONCIONS DE LA FORMULATION MASSIF -> P
  3128. C COQUE PINF PSUP
  3129. JGOBL=1
  3130. TABOBL(1)='PR '
  3131. ELSE
  3132. IRET = 0
  3133. CALL ERREUR(5)
  3134. ENDIF
  3135. GOTO 9999
  3136. ENDIF
  3137.  
  3138. *--------------------------------------------------------------------
  3139. 9999 CONTINUE
  3140.  
  3141. IF(IMECA.GT.0 .OR. IPORE.GT.0)THEN
  3142. C CB215821 : Ajout de TREF et TALP à la fin des compostantes facultatives
  3143. C pour les FORMULATION MECANIQUE, POREUX
  3144. C Pour ne pas casser l'ordre dans COMP (com2/coml6/comval/comara/...)
  3145. JGFAC=JGFAC+2
  3146. TABFAC(JGFAC-1) ='TREF'
  3147. TABFAC(JGFAC) ='TALP'
  3148. ENDIF
  3149.  
  3150. C Erreur si JGOBL ou JGFAC sont superieurs a ITA
  3151. C (Passage en FORTRAN 77 car la compilation depasse la memoire sur
  3152. C Windows-32bits)
  3153. IF ((JGOBL .GT. ITA) .OR. (JGFAC .GT. ITA)) THEN
  3154. IRET = 0
  3155. CALL ERREUR(5)
  3156. RETURN
  3157. ENDIF
  3158.  
  3159. NBROBL = JGOBL
  3160. NBRFAC = JGFAC
  3161. NOMID = 0
  3162.  
  3163. * SI PROBLEME (IRET = 0), ON SORT AVEC IPNOMC A 0
  3164. *
  3165. IF (IRET.NE.0) THEN
  3166. SEGINI,NOMID
  3167. DO 100 IO = 1,NBROBL
  3168. NOMID.LESOBL(IO) = TABOBL(IO)
  3169. 100 CONTINUE
  3170.  
  3171. DO 110 IO=1,NBRFAC
  3172. NOMID.LESFAC(IO) = TABFAC(IO)
  3173. 110 CONTINUE
  3174. SEGACT,NOMID*NOMOD
  3175. ENDIF
  3176.  
  3177. IPNOMC = NOMID
  3178.  
  3179. END
  3180.  
  3181.  
  3182.  
  3183.  

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