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

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