Télécharger iddein.eso

Retour à la liste

Numérotation des lignes :

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

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