Télécharger iddein.eso

Retour à la liste

Numérotation des lignes :

  1. C IDDEIN SOURCE FANDEUR 16/12/05 21:39:39 9211
  2.  
  3. C=======================================================================
  4. C= RECHERCHE DES NOMS DE COMPOSANTES DE DEFORMATIONS INELASTIQUES =
  5. C= -------------------------------------------------------------- =
  6. C= =
  7. C= Entrees : =
  8. C= IPMODL pointeur ACTIF sur objet IMODEL =
  9. C= IFOU valeur de IFOUR de CCOPTIO =
  10. C= =
  11. C= Sorties : =
  12. C= IPNOMC pointeur de type NOMID sur les listes de noms de =
  13. C= composantes OBLigatoires et FACultatives =
  14. C= NBROBL nombre de composantes OBLigatoires =
  15. C= NBRFAC nombre de composantes FACultatives =
  16. C=======================================================================
  17.  
  18. SUBROUTINE IDDEIN (IPMODL,IFOU,IPNOMC,NBROBL,NBRFAC)
  19.  
  20. IMPLICIT INTEGER(I-N)
  21. IMPLICIT REAL*8(A-H,O-Z)
  22.  
  23. -INC CCOPTIO
  24. -INC SMMODEL
  25.  
  26. IPNOMC=0
  27. NBROBL=0
  28. NBRFAC=0
  29. NOMID =0
  30.  
  31. IMODEL=IPMODL
  32.  
  33. MELE=NEFMOD
  34. MFR=NUMMFR(MELE)
  35.  
  36. NFOR=FORMOD(/2)
  37. CALL PLACE(FORMOD,NFOR,ITHHY,'THERMOHYDRIQUE')
  38. CALL PLACE(FORMOD,NFOR,ITHER,'THERMIQUE')
  39. CALL PLACE(FORMOD,NFOR,IMAGN,'MAGNETODYNAMIQUE')
  40. CALL PLACE(FORMOD,NFOR,IELEC,'ELECTROSTATIQUE')
  41. CALL PLACE(FORMOD,NFOR,IDIFF,'DIFFUSION')
  42.  
  43. mfr2=0
  44. IF (ITHHY.NE.0) mfr2=65
  45. if (ITHER.NE.0) mfr2=29
  46. if (IMAGN.NE.0) mfr2=69
  47. IF (IELEC.NE.0) mfr2=71
  48. IF (IDIFF.NE.0) mfr2=73
  49. IF (mfr2.NE.0) GOTO 100
  50.  
  51. C Coque integree ou pas ?
  52. NPINT=0
  53. IF (INFMOD(/1).NE.0) NPINT=INFMOD(1)
  54. *
  55. * -1/ MACRO ELEMENT
  56. *
  57. IF (MFR.EQ.61)THEN
  58. C
  59. IF (IFOU.EQ.-1.OR.IFOU.EQ.-2) THEN
  60. *
  61. * 0-A/ CONTRAINTES PLANES - DEFORMATIONS PLANES
  62. *
  63. NBROBL=5
  64. SEGINI NOMID
  65. LESOBL(1)='EIS1'
  66. LESOBL(2)='CIZ1'
  67. LESOBL(3)='GIXY'
  68. LESOBL(4)='EIS2'
  69. LESOBL(5)='CIZ2'
  70. ENDIF
  71.  
  72. C 1 - Element JOINT
  73. C ===================
  74. ELSE IF (MFR.EQ.35) THEN
  75. C =====
  76. C 1.1 - Tridimensionnel
  77. C =====
  78. IF (IFOU.EQ.2) THEN
  79. NBROBL=3
  80. SEGINI,NOMID
  81. LESOBL(1)='FRS1'
  82. LESOBL(2)='FRS2'
  83. LESOBL(3)='FRN '
  84. C =====
  85. C 1.2 - Bidimensionnel PLAN (CP/DP/DPGE) et Axisymetrie
  86. C =====
  87. ELSE IF ((IFOU.EQ.-3).OR.(IFOU.EQ.-2).OR.(IFOU.EQ.-1).OR.
  88. . (IFOU.EQ.0)) THEN
  89. NBROBL=2
  90. SEGINI,NOMID
  91. LESOBL(1)='FRSN'
  92. LESOBL(2)='FRN '
  93. ENDIF
  94.  
  95. C 2 - Element JOINT CISAILLEMENT (2D)
  96. C Bidimensionnel PLAN (CP/DP/DPGE)
  97. C ======================================
  98. ELSE IF (MFR.EQ.53) THEN
  99. IF ((IFOU.EQ.-3).OR.(IFOU.EQ.-2).OR.(IFOU.EQ.-1)) THEN
  100. NBROBL=2
  101. SEGINI,NOMID
  102. LESOBL(1)='FRS1'
  103. LESOBL(2)='FRS2'
  104. ENDIF
  105.  
  106. C 3 - Element JOINT GENERALISE
  107. C ==============================
  108. ELSE IF (MFR.EQ.55) THEN
  109. C =====
  110. C 3.1 - Bidimensionnel PLAN (CP/DP/DPGE)
  111. C =====
  112. IF ((IFOU.EQ.-2).OR.(IFOU.EQ.-1).OR.(IFOU.EQ.-3)) THEN
  113. NBROBL=4
  114. SEGINI,NOMID
  115. LESOBL(1)='FRSN'
  116. LESOBL(2)='FRN '
  117. LESOBL(3)='EITT'
  118. LESOBL(4)='EIZZ'
  119. C =====
  120. C 3.2 - Tridimensionnel
  121. C =====
  122. ELSE IF (IFOU.EQ.2) THEN
  123. NBROBL=6
  124. SEGINI,NOMID
  125. LESOBL(1)='FRS1'
  126. LESOBL(2)='FRS2'
  127. LESOBL(3)='FRN '
  128. LESOBL(4)='EITT'
  129. LESOBL(5)='EIUU'
  130. LESOBL(6)='GITU'
  131. ENDIF
  132.  
  133. C 4 - Elements MASSIFS
  134. C ======================
  135. ELSE IF (MFR.EQ.1.OR.MFR.EQ.31.or.MFR.EQ.63) THEN
  136. C =====
  137. C 4.1 - Bidimensionnel PLAN (CP/DP/DPGE)
  138. C =====
  139. IF (IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
  140. NBROBL=4
  141. SEGINI,NOMID
  142. LESOBL(1)='EIXX'
  143. LESOBL(2)='EIYY'
  144. LESOBL(3)='EIZZ'
  145. LESOBL(4)='GIXY'
  146. C =====
  147. C 4.2 - Axisymetrie
  148. C =====
  149. ELSE IF (IFOU.EQ.0) THEN
  150. NBROBL=4
  151. SEGINI,NOMID
  152. LESOBL(1)='EIRR'
  153. LESOBL(2)='EIZZ'
  154. LESOBL(3)='EITT'
  155. LESOBL(4)='GIRZ'
  156. C =====
  157. C 4.3 - Fourier
  158. C =====
  159. ELSE IF (IFOU.EQ.1) THEN
  160. NBROBL=6
  161. SEGINI,NOMID
  162. LESOBL(1)='EIRR'
  163. LESOBL(2)='EIZZ'
  164. LESOBL(3)='EITT'
  165. LESOBL(4)='GIRZ'
  166. LESOBL(5)='GIRT'
  167. LESOBL(6)='GIZT'
  168. C =====
  169. C 4.4 - Tridimensionnel
  170. C =====
  171. ELSE IF (IFOU.EQ.2) THEN
  172. NBROBL=6
  173. SEGINI,NOMID
  174. LESOBL(1)='EIXX'
  175. LESOBL(2)='EIYY'
  176. LESOBL(3)='EIZZ'
  177. LESOBL(4)='GIXY'
  178. LESOBL(5)='GIXZ'
  179. LESOBL(6)='GIYZ'
  180. C =====
  181. C 4.5 - Unidimensionnel (1D)
  182. C =====
  183. ELSE IF (IFOU.GE.3.AND.IFOU.LE.15.AND.MFR.EQ.1) THEN
  184. NBROBL=3
  185. SEGINI,NOMID
  186. IF (IFOU.LE.11) THEN
  187. LESOBL(1)='EIXX'
  188. LESOBL(2)='EIYY'
  189. LESOBL(3)='EIZZ'
  190. ELSE IF (IFOU.GE.12.AND.IFOU.LE.14) THEN
  191. LESOBL(1)='EIRR'
  192. LESOBL(2)='EIZZ'
  193. LESOBL(3)='EITT'
  194. ELSE IF (IFOU.EQ.15) THEN
  195. LESOBL(1)='EIRR'
  196. C Voir TC pour LESOBL(2)='EIFF'
  197. LESOBL(2)='EIZZ'
  198. LESOBL(3)='EITT'
  199. ENDIF
  200. ENDIF
  201.  
  202. C 5 - Elements COQUES
  203. C =====================
  204. ELSE IF (MFR.EQ.3.OR.MFR.EQ.39) THEN
  205. IF (NPINT.EQ.0) THEN
  206. C =====
  207. C 5.1 - Tridimensionnel et Fourier
  208. C =====
  209. IF (IFOU.EQ.2.OR.IFOU.EQ.1) THEN
  210. NBROBL=6
  211. SEGINI,NOMID
  212. LESOBL(1)='EISS'
  213. LESOBL(2)='EITT'
  214. LESOBL(3)='GIST'
  215. LESOBL(4)='RISS'
  216. LESOBL(5)='RITT'
  217. LESOBL(6)='RIST'
  218. C =====
  219. C 5.2 - Axisymetrie
  220. C =====
  221. ELSE IF (IFOU.EQ.0) THEN
  222. NBROBL=4
  223. SEGINI,NOMID
  224. LESOBL(1)='EISS'
  225. LESOBL(2)='EITT'
  226. LESOBL(3)='RISS'
  227. LESOBL(4)='RITT'
  228. C =====
  229. C 5.3 - Bidimensionnel PLAN (CP/DP/DPGE)
  230. C =====
  231. ELSE IF (IFOU.EQ.-1.OR.IFOU.EQ.-2.OR.IFOU.EQ.-3) THEN
  232. NBROBL=4
  233. SEGINI,NOMID
  234. LESOBL(1)='EISS'
  235. LESOBL(2)='EIZZ'
  236. LESOBL(3)='RISS'
  237. LESOBL(4)='RIZZ'
  238. ENDIF
  239. ELSE
  240. C =====
  241. C 5.4 - Tridimensionnel - Cas particulier
  242. C =====
  243. IF (MFR.EQ.3.AND.IFOU.EQ.2) THEN
  244. NBROBL=4
  245. SEGINI,NOMID
  246. LESOBL(1)='EISS'
  247. LESOBL(2)='EITT'
  248. LESOBL(3)='EINN'
  249. LESOBL(4)='GIST'
  250. ENDIF
  251. ENDIF
  252.  
  253. C 6 - Elements de COQUE EPAISSE
  254. C Tridimensionnel, axisymetrie et Fourier
  255. C =============================================
  256. ELSE IF (MFR.EQ.5) THEN
  257. NBROBL=5
  258. SEGINI,NOMID
  259. LESOBL(1)='EISS'
  260. LESOBL(2)='EITT'
  261. LESOBL(3)='GIST'
  262. LESOBL(4)='GISN'
  263. LESOBL(5)='GITN'
  264.  
  265. C 7 - LIA2 : Element 3D de liaison a 2 noeuds
  266. C =============================================
  267. ELSE IF (MFR.EQ.51) THEN
  268. NBROBL=6
  269. SEGINI,NOMID
  270. LESOBL(1)='FEPS'
  271. LESOBL(2)='FGXY'
  272. LESOBL(3)='FGXZ'
  273. LESOBL(4)='FFCX'
  274. LESOBL(5)='FFCY'
  275. LESOBL(6)='FFCZ'
  276. C 7bis - JOI1: Element 3D de liaison a 2 noeuds
  277. C =============================================
  278. ELSE IF (MFR.EQ.75) THEN
  279. IF (IFOU.EQ.2) THEN
  280. NBROBL=6
  281. SEGINI,NOMID
  282. LESOBL(1)='FEXX '
  283. LESOBL(2)='FGXY '
  284. LESOBL(3)='FGXZ '
  285. LESOBL(4)='FCXX '
  286. LESOBL(5)='FCXY '
  287. LESOBL(6)='FCXZ '
  288. ELSE IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
  289. NBROBL=3
  290. SEGINI,NOMID
  291. LESOBL(1)='FEXX '
  292. LESOBL(2)='FGXY '
  293. LESOBL(3)='FCXZ '
  294. ENDIF
  295.  
  296. C 8 - BAEX : BARRE EXCENTREE
  297. C ============================
  298. ELSE IF (MFR.EQ.49) THEN
  299. NBROBL=1
  300. SEGINI,NOMID
  301. LESOBL(1)='EISS'
  302.  
  303. C 9 - Elements POUTRE et TUYAU
  304. C ==============================
  305. ELSE IF (MFR.EQ.7.OR.MFR.EQ.13) THEN
  306. C =====
  307. C 9.1 - Tridimensionnel
  308. C =====
  309. IF (IFOU.EQ.2) THEN
  310. NBROBL=6
  311. SEGINI,NOMID
  312. LESOBL(1)='FEPS'
  313. LESOBL(2)='FGXY'
  314. LESOBL(3)='FGXZ'
  315. LESOBL(4)='FFCX'
  316. LESOBL(5)='FFCY'
  317. LESOBL(6)='FFCZ'
  318. C =====
  319. C 9.2 - Bidimensionnel PLAN (CP/DP/DPGE)
  320. C =====
  321. ELSE IF (IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
  322. NBROBL=3
  323. SEGINI,NOMID
  324. LESOBL(1)='FEPS'
  325. LESOBL(2)='FGXY'
  326. LESOBL(3)='FFCZ'
  327. ENDIF
  328.  
  329. C 10 - Element de COQUE en CISAILLEMENT TRANSVERSE
  330. C ==================================================
  331. ELSE IF (MFR.EQ.9) THEN
  332. IF (NPINT.EQ.0) THEN
  333. C ======
  334. C 10.1 - Tridimensionnel
  335. C ======
  336. IF (IFOU.EQ.2) THEN
  337. NBROBL=8
  338. SEGINI,NOMID
  339. LESOBL(1)='EISS'
  340. LESOBL(2)='EITT'
  341. LESOBL(3)='GIST'
  342. LESOBL(4)='RISS'
  343. LESOBL(5)='RITT'
  344. LESOBL(6)='RIST'
  345. LESOBL(7)='GISN'
  346. LESOBL(8)='GITN'
  347. ENDIF
  348. ENDIF
  349.  
  350. C 11 - Elements LINESPRING (tridimensionnel)
  351. C ============================================
  352. ELSE IF (MFR.EQ.15) THEN
  353. IF (IFOU.EQ.2) THEN
  354. NBROBL=6
  355. SEGINI,NOMID
  356. LESOBL(1)='EIZZ'
  357. LESOBL(2)='GIXZ'
  358. LESOBL(3)='GIYZ'
  359. LESOBL(4)='RIXX'
  360. LESOBL(5)='RIZZ'
  361. LESOBL(6)='FDJP'
  362. ENDIF
  363.  
  364. C 12 - Elements de MEMBRANE
  365. C ===========================
  366. ELSE IF (MFR.EQ.25) THEN
  367. NBROBL=3
  368. SEGINI,NOMID
  369. LESOBL(1)='EISS'
  370. LESOBL(2)='EITT'
  371. LESOBL(3)='GIST'
  372.  
  373. C 13 - Formulation UNIAXIALE
  374. C =============================
  375. ELSE IF (MFR.EQ.27) THEN
  376. NBROBL=1
  377. SEGINI,NOMID
  378. LESOBL(1)='EISS'
  379.  
  380. C 14 - Formulation THERMIQUE
  381. C ============================
  382. ELSE IF (MFR.EQ.29) THEN
  383. C ======
  384. C 14.1 - Tridimensionnel
  385. C ======
  386. IF (IFOU.EQ.2) THEN
  387. NBROBL=3
  388. SEGINI,NOMID
  389. LESOBL(1)='FRDX'
  390. LESOBL(2)='FRDY'
  391. LESOBL(3)='FRDZ'
  392. C ======
  393. C 14.2 - Bidimensionnel PLAN (CP/DP/DPGE)
  394. C ======
  395. ELSE IF (IFOU.EQ.-1.OR.IFOU.EQ.-2.OR.IFOU.EQ.-3) THEN
  396. NBROBL=2
  397. SEGINI,NOMID
  398. LESOBL(1)='FRDX'
  399. LESOBL(2)='FRDY'
  400. C ======
  401. C 14.3 - Axisymetrie et Fourier
  402. C ======
  403. ELSE IF (IFOU.EQ.1.OR.IFOU.EQ.0) THEN
  404. NBROBL=2
  405. SEGINI,NOMID
  406. LESOBL(1)='FRDX'
  407. LESOBL(2)='FRDY'
  408. C*OF LESOBL(1)='FRDR'
  409. C*OF LESOBL(2)='FRDZ'
  410. C ======
  411. C 14.4 - Unidimensionnel (1D)
  412. C ======
  413. ELSE IF (IFOU.GE.3.AND.IFOU.LE.15) THEN
  414. NBROBL=1
  415. SEGINI,NOMID
  416. C*OF IF (IFOU.LE.11) THEN
  417. LESOBL(1)='FRDX'
  418. C*OF ELSE
  419. C*OF LESOBL(1)='FRDR'
  420. C*OF ENDIF
  421. ENDIF
  422.  
  423. C 15 - Formulation POREUX
  424. C =========================
  425. ELSE IF (MFR.EQ.33) THEN
  426. C =====
  427. C 15.A - Elements POREUX
  428. C ======
  429. IF (MELE.GE.79.AND.MELE.LE.83) THEN
  430. C ========
  431. C 15.A.1 - Bidimensionnel PLAN (CP/DP)
  432. C ========
  433. IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
  434. NBROBL=5
  435. SEGINI,NOMID
  436. LESOBL(1)='EIXX'
  437. LESOBL(2)='EIYY'
  438. LESOBL(3)='EIZZ'
  439. LESOBL(4)='GIXY'
  440. LESOBL(5)='PIDE'
  441. C ========
  442. C 15.A.2 - Axisymetrie
  443. C ========
  444. ELSE IF (IFOU.EQ.0) THEN
  445. NBROBL=5
  446. SEGINI,NOMID
  447. LESOBL(1)='EIRR'
  448. LESOBL(2)='EIZZ'
  449. LESOBL(3)='EITT'
  450. LESOBL(4)='GIRZ'
  451. LESOBL(5)='PIDE'
  452. C ========
  453. C 15.A.3 - Fourier
  454. C ========
  455. ELSE IF (IFOU.EQ.1) THEN
  456. NBROBL=7
  457. SEGINI,NOMID
  458. LESOBL(1)='EIRR'
  459. LESOBL(2)='EIZZ'
  460. LESOBL(3)='EITT'
  461. LESOBL(4)='GIRZ'
  462. LESOBL(5)='GIRT'
  463. LESOBL(6)='GIZT'
  464. LESOBL(7)='PIDE'
  465. C ========
  466. C 15.A.4 - Tridimensionnel et bidimensionnel DEFO PLAN GENE
  467. C ========
  468. ELSE IF (IFOU.EQ.2.OR.IFOU.EQ.-3) THEN
  469. NBROBL=7
  470. SEGINI,NOMID
  471. LESOBL(1)='EIXX'
  472. LESOBL(2)='EIYY'
  473. LESOBL(3)='EIZZ'
  474. LESOBL(4)='GIXY'
  475. LESOBL(5)='GIXZ'
  476. LESOBL(6)='GIYZ'
  477. LESOBL(7)='PIDE'
  478. ENDIF
  479. C ======
  480. C 15.B - Elements JOINT POREUX
  481. C ======
  482. ELSE IF (MELE.GE.108.AND.MELE.LE.110) THEN
  483. C ========
  484. C 15.B.1 - Tridimensionnel
  485. C ========
  486. IF (IFOU.EQ.2) THEN
  487. NBROBL=4
  488. SEGINI,NOMID
  489. LESOBL(1)='FRS1'
  490. LESOBL(2)='FRS2'
  491. LESOBL(3)='FRN '
  492. LESOBL(4)='PIDE'
  493. C ========
  494. C 15.B.2 - Bidimensionnel PLAN (CP/DP) et Axisymetrie
  495. C ========
  496. ELSE IF (IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.0) THEN
  497. NBROBL=3
  498. SEGINI,NOMID
  499. LESOBL(1)='FRSN'
  500. LESOBL(2)='FRN '
  501. LESOBL(3)='PIDE'
  502. ENDIF
  503. ENDIF
  504.  
  505. C 16 - Formulation POREUX type Q
  506. C ================================
  507. ELSE IF (MFR.EQ.57) THEN
  508. C =====
  509. C 16.A - Elements POREUX Q
  510. C ======
  511. IF (MELE.GE.173.AND.MELE.LE.177) THEN
  512. C ========
  513. C 16.A.1 - Bidimensionnel PLAN (CP/DP)
  514. C ========
  515. IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
  516. NBROBL=6
  517. SEGINI,NOMID
  518. LESOBL(1)='EIXX'
  519. LESOBL(2)='EIYY'
  520. LESOBL(3)='EIZZ'
  521. LESOBL(4)='GIXY'
  522. LESOBL(5)='PIDE'
  523. LESOBL(6)='QIDE'
  524. C ========
  525. C 16.A.2 - Axisymetrie
  526. C ========
  527. ELSE IF (IFOU.EQ.0) THEN
  528. NBROBL=6
  529. SEGINI,NOMID
  530. LESOBL(1)='EIRR'
  531. LESOBL(2)='EIZZ'
  532. LESOBL(3)='EITT'
  533. LESOBL(4)='GIRZ'
  534. LESOBL(5)='PIDE'
  535. LESOBL(6)='QIDE'
  536. C ========
  537. C 16.A.3 - Fourier
  538. C ========
  539. ELSE IF (IFOU.EQ.1) THEN
  540. NBROBL=8
  541. SEGINI,NOMID
  542. LESOBL(1)='EIRR'
  543. LESOBL(2)='EIZZ'
  544. LESOBL(3)='EITT'
  545. LESOBL(4)='GIRZ'
  546. LESOBL(5)='GIRT'
  547. LESOBL(6)='GIZT'
  548. LESOBL(7)='PIDE'
  549. LESOBL(8)='QIDE'
  550. C ========
  551. C 16.A.4 - Tridimensionnel et bidimensionnel DEFO PLAN GENE
  552. C ========
  553. ELSE IF (IFOU.EQ.2.OR.IFOU.EQ.-3) THEN
  554. NBROBL=8
  555. SEGINI,NOMID
  556. LESOBL(1)='EIXX'
  557. LESOBL(2)='EIYY'
  558. LESOBL(3)='EIZZ'
  559. LESOBL(4)='GIXY'
  560. LESOBL(5)='GIXZ'
  561. LESOBL(6)='GIYZ'
  562. LESOBL(7)='PIDE'
  563. LESOBL(8)='QIDE'
  564. ENDIF
  565. C ======
  566. C 16.B - Elements JOINT POREUX Q
  567. C ======
  568. ELSE IF (MELE.GE.185.AND.MELE.LE.187) THEN
  569. C ========
  570. C 16.B.1 - Tridimensionnel
  571. C ========
  572. IF (IFOU.EQ.2) THEN
  573. NBROBL=5
  574. SEGINI,NOMID
  575. LESOBL(1)='FRS1'
  576. LESOBL(2)='FRS2'
  577. LESOBL(3)='FRN '
  578. LESOBL(4)='PIDE'
  579. LESOBL(5)='QIDE'
  580. C ========
  581. C 16.B.2 - Bidimensionnel PLAN (CP/DP) et Axisymetrie
  582. C ========
  583. ELSE IF (IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.0) THEN
  584. NBROBL=4
  585. SEGINI,NOMID
  586. LESOBL(1)='FRSN'
  587. LESOBL(2)='FRN '
  588. LESOBL(3)='PIDE'
  589. LESOBL(4)='QIDE'
  590. ENDIF
  591. ENDIF
  592.  
  593. C 17 - Formulation POREUX type R
  594. C ================================
  595. ELSE IF (MFR.EQ.59) THEN
  596. C =====
  597. C 17.A - Elements POREUX R
  598. C ======
  599. IF (MELE.GE.178.AND.MELE.LE.182) THEN
  600. C ========
  601. C 17.A.1 - Bidimensionnel PLAN (CP/DP)
  602. C ========
  603. IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
  604. NBROBL=7
  605. SEGINI,NOMID
  606. LESOBL(1)='EIXX'
  607. LESOBL(2)='EIYY'
  608. LESOBL(3)='EIZZ'
  609. LESOBL(4)='GIXY'
  610. LESOBL(5)='PIDE'
  611. LESOBL(6)='QIDE'
  612. LESOBL(7)='TIDE'
  613. C ========
  614. C 17.A.2 - Axisymetrie
  615. C ========
  616. ELSE IF (IFOU.EQ.0) THEN
  617. NBROBL=7
  618. SEGINI,NOMID
  619. LESOBL(1)='EIRR'
  620. LESOBL(2)='EIZZ'
  621. LESOBL(3)='EITT'
  622. LESOBL(4)='GIRZ'
  623. LESOBL(5)='PIDE'
  624. LESOBL(6)='QIDE'
  625. LESOBL(7)='TIDE'
  626. C ========
  627. C 17.A.3 - Fourier
  628. C ========
  629. ELSE IF (IFOU.EQ.1) THEN
  630. NBROBL=9
  631. SEGINI,NOMID
  632. LESOBL(1)='EIRR'
  633. LESOBL(2)='EIZZ'
  634. LESOBL(3)='EITT'
  635. LESOBL(4)='GIRZ'
  636. LESOBL(5)='GIRT'
  637. LESOBL(6)='GIZT'
  638. LESOBL(7)='PIDE'
  639. LESOBL(8)='QIDE'
  640. LESOBL(9)='TIDE'
  641. C ========
  642. C 17.A.4 - Tridimensionnel et bidimensionnel DEFO PLAN GENE
  643. C ========
  644. ELSE IF (IFOU.EQ.2.OR.IFOU.EQ.-3) THEN
  645. NBROBL=9
  646. SEGINI,NOMID
  647. LESOBL(1)='EIXX'
  648. LESOBL(2)='EIYY'
  649. LESOBL(3)='EIZZ'
  650. LESOBL(4)='GIXY'
  651. LESOBL(5)='GIXZ'
  652. LESOBL(6)='GIYZ'
  653. LESOBL(7)='PIDE'
  654. LESOBL(8)='QIDE'
  655. LESOBL(9)='TIDE'
  656. ENDIF
  657. C ======
  658. C 17.B - Elements JOINT POREUX R
  659. C ======
  660. ELSE IF (MELE.GE.188.AND.MELE.LE.190) THEN
  661. C ========
  662. C 17.B.1 - Tridimensionnel
  663. C ========
  664. IF (IFOU.EQ.2) THEN
  665. NBROBL=6
  666. SEGINI,NOMID
  667. LESOBL(1)='FRS1'
  668. LESOBL(2)='FRS2'
  669. LESOBL(3)='FRN '
  670. LESOBL(4)='PIDE'
  671. LESOBL(5)='QIDE'
  672. LESOBL(6)='TIDE'
  673. C ========
  674. C 17.B.2 - Bidimensionnel PLAN (CP/DP) et Axisymetrie
  675. C ========
  676. ELSE IF (IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.0) THEN
  677. NBROBL=5
  678. SEGINI,NOMID
  679. LESOBL(1)='FRSN'
  680. LESOBL(2)='FRN '
  681. LESOBL(3)='PIDE'
  682. LESOBL(4)='QIDE'
  683. LESOBL(5)='TIDE'
  684. ENDIF
  685. ENDIF
  686.  
  687. C 18 - Elements TUYAU FISSURE (tridimensionnel)
  688. C ===============================================
  689. ELSE IF (MFR.EQ.17) THEN
  690. IF (IFOU.EQ.2) THEN
  691. NBROBL=8
  692. SEGINI,NOMID
  693. LESOBL(1)='FEPS'
  694. LESOBL(2)='FGXY'
  695. LESOBL(3)='FGXZ'
  696. LESOBL(4)='FFCX'
  697. LESOBL(5)='FFCY'
  698. LESOBL(6)='FFCZ'
  699. LESOBL(7)='FPS7'
  700. LESOBL(8)='FPS8'
  701. ENDIF
  702.  
  703. C 19 - Elements de SECTION/FIBRE
  704. C ================================
  705. ELSE IF (MFR.EQ.47) THEN
  706. C ======
  707. C 19.1 - Bidimensionnel PLAN (CP/DP/DPGE)
  708. C ======
  709. IF (IFOU.EQ.-1.OR.IFOU.EQ.-2.OR.IFOU.EQ.-3) THEN
  710. NBROBL=2
  711. SEGINI,NOMID
  712. LESOBL(1)='EIXX'
  713. LESOBL(2)='GIXY'
  714. C*OF LESOBL(3)='GIXZ'
  715. C ======
  716. C 19.2 - Tridimensionnel
  717. C ======
  718. ELSE IF (IFOU.EQ.2) THEN
  719. NBROBL=3
  720. SEGINI,NOMID
  721. LESOBL(1)='EIXX'
  722. LESOBL(2)='GIXY'
  723. LESOBL(3)='GIXZ'
  724. ENDIF
  725.  
  726. C 20 - Elements de zones cohesives
  727. C ==================================
  728. ELSE IF (MFR.EQ.77) THEN
  729. C =====
  730. C 20.1 - Tridimensionnel
  731. C =====
  732. IF (IFOU.GT.0) THEN
  733. NBROBL=3
  734. SEGINI,NOMID
  735. LESOBL(1)='FRS1'
  736. LESOBL(2)='FRS2'
  737. LESOBL(3)='FRN '
  738. C =====
  739. C 20.2 - Bidimensionnel PLAN (CP/DP/DPGE) et Axisymetrie
  740. C =====
  741. ELSE IF (IFOU.LE.0) THEN
  742. NBROBL=2
  743. SEGINI,NOMID
  744. LESOBL(1)='FRSN'
  745. LESOBL(2)='FRN '
  746. ENDIF
  747. c
  748. c element coaxial COS2 (3D pour liaison acier-beton)
  749. c
  750. ELSEIF ( MFR.EQ.78) then
  751. NBROBL= 3
  752. SEGINI NOMID
  753. LESOBL(1)='GINT'
  754. LESOBL(2)='GINN'
  755. LESOBL(3)='GINM'
  756. ENDIF
  757. C ===========================
  758. C Par DEFAUT : segment VIDE
  759. C ===========================
  760. 100 continue
  761. IF (NOMID.EQ.0) SEGINI,NOMID
  762.  
  763. SEGDES,NOMID
  764. IPNOMC=NOMID
  765.  
  766. RETURN
  767. END
  768.  
  769.  
  770.  

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