Télécharger iddepl.eso

Retour à la liste

Numérotation des lignes :

  1. C IDDEPL SOURCE MB234859 16/12/14 21:15:11 9252
  2.  
  3. C=======================================================================
  4. C= RECHERCHE DES NOMS DE COMPOSANTES DE DEPLACEMENTS =
  5. C= ------------------------------------------------- =
  6. C= =
  7. C= Entrees : =
  8. C= MFR Numero de la FORMULATION element fini =
  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 IDDEPL (MFR,IFOU,IPNOMC,NBROBL,NBRFAC)
  19.  
  20. IMPLICIT INTEGER(I-N)
  21. IMPLICIT REAL*8(A-H,O-Z)
  22.  
  23. -INC CCOPTIO
  24.  
  25. SEGMENT NOMID
  26. CHARACTER*8 LESOBL(NBROBL),LESFAC(NBRFAC)
  27. ENDSEGMENT
  28.  
  29. NBROBL=0
  30. NBRFAC=0
  31. NOMID=0
  32.  
  33. IF (MFR.EQ.65) THEN
  34. NBROBL=3
  35. SEGINI,NOMID
  36. LESOBL(1)='PG '
  37. LESOBL(2)='PC '
  38. LESOBL(3)='T '
  39. *
  40. * 0/ MACRO ELEMENT
  41. *
  42. ELSEIF (MFR.EQ.61)THEN
  43. C
  44. IF (IFOU.EQ.-1.OR.IFOU.EQ.-2) THEN
  45. *
  46. * 0-A/ CONTRAINTES PLANES - DEFORMATIONS PLANES
  47. *
  48. NBROBL=5
  49. SEGINI NOMID
  50. LESOBL(1)='UX '
  51. LESOBL(2)='UY '
  52. LESOBL(3)='RZ '
  53. LESOBL(4)='UM '
  54. LESOBL(5)='RM '
  55. ENDIF
  56.  
  57. C 1 - Elements COQUE (3), COQUE EPAISSE (5), POUTRE (7), COQUE en
  58. C CISAILLEMENT TRANSVERSE (9), TUYAU (13), LINESPRING (15), TUYAU
  59. C FISSURE (17), Barre excentree BAEX (49), LIA2 (51), JOI1(75)
  60. C ====================================================================
  61. ELSE IF (MFR.EQ. 3.OR.MFR.EQ. 5.OR.MFR.EQ. 7.OR.MFR.EQ. 9.OR.
  62. . MFR.EQ.13.OR.MFR.EQ.15.OR.MFR.EQ.17.OR.MFR.EQ.49.OR.
  63. . MFR.EQ.51.OR.MFR.EQ.75.OR.MFR.EQ.74) THEN
  64. C =====
  65. C 1.1 - Tridimensionnel
  66. C =====
  67. IF (IFOU.EQ.2) THEN
  68. NBROBL=6
  69. SEGINI,NOMID
  70. LESOBL(1)='UX '
  71. LESOBL(2)='UY '
  72. LESOBL(3)='UZ '
  73. LESOBL(4)='RX '
  74. LESOBL(5)='RY '
  75. LESOBL(6)='RZ '
  76. C =====
  77. C 1.2 - Fourier
  78. C =====
  79. ELSE IF (IFOU.EQ.1) THEN
  80. NBROBL=4
  81. SEGINI,NOMID
  82. LESOBL(1)='UR '
  83. LESOBL(2)='UZ '
  84. LESOBL(3)='UT '
  85. LESOBL(4)='RT '
  86. C =====
  87. C 1.3 - Axisymetrie
  88. C =====
  89. ELSE IF (IFOU.EQ.0) THEN
  90. NBROBL=3
  91. SEGINI,NOMID
  92. LESOBL(1)='UR '
  93. LESOBL(2)='UZ '
  94. LESOBL(3)='RT '
  95. C =====
  96. C 1.4 - Bidimensionnel PLAN (CP/DP)
  97. C =====
  98. ELSE IF (IFOU.EQ.-1.OR.IFOU.EQ.-2) THEN
  99. NBROBL=3
  100. SEGINI,NOMID
  101. LESOBL(1)='UX '
  102. LESOBL(2)='UY '
  103. LESOBL(3)='RZ '
  104. C =====
  105. C 1.5 - Bidimensionnel PLAN GENE
  106. C =====
  107. ELSE IF (IFOU.EQ.-3) THEN
  108. C Ici il faut distinguer les formulations :
  109. IF (MFR.EQ.03 .OR. MFR.EQ.05) THEN
  110. NBROBL=6
  111. SEGINI,NOMID
  112. LESOBL(1)='UX '
  113. LESOBL(2)='UY '
  114. LESOBL(3)='RZ '
  115. LESOBL(4)='UZ '
  116. LESOBL(5)='RY '
  117. LESOBL(6)='RX '
  118. ELSE
  119. NBROBL=3
  120. SEGINI,NOMID
  121. LESOBL(1)='UX '
  122. LESOBL(2)='UY '
  123. LESOBL(3)='RZ '
  124. ENDIF
  125. ENDIF
  126.  
  127. C 2 - Elements LIQUIDE
  128. C ======================
  129. ELSE IF (MFR.EQ.11) THEN
  130. NBROBL=2
  131. SEGINI,NOMID
  132. LESOBL(1)='P '
  133. LESOBL(2)='PI '
  134.  
  135. C 3 - Elements TUYAU ACOUSTIQUE PUR
  136. C ===================================
  137. ELSE IF (MFR.EQ.41) THEN
  138. NBROBL=2
  139. SEGINI,NOMID
  140. LESOBL(1)='PI '
  141. LESOBL(2)='P '
  142.  
  143. C 4 - Element de RACCORD LITU
  144. C =============================
  145. ELSE IF (MFR.EQ.43) THEN
  146. NBROBL=5
  147. SEGINI,NOMID
  148. LESOBL(1)='UX '
  149. LESOBL(2)='UY '
  150. LESOBL(3)='UZ '
  151. LESOBL(4)='PI '
  152. LESOBL(5)='P '
  153.  
  154. C 5 - Elements HOMOGENEISE
  155. C ==========================
  156. ELSE IF (MFR.EQ.37) THEN
  157. C =====
  158. C 5.1 - Fourier
  159. C =====
  160. IF (IFOU.EQ.1) THEN
  161. NBROBL=6
  162. SEGINI,NOMID
  163. LESOBL(1)='P '
  164. LESOBL(2)='PI '
  165. LESOBL(3)='UR '
  166. LESOBL(4)='RT '
  167. LESOBL(5)='UT '
  168. LESOBL(6)='RR '
  169. C =====
  170. C 5.2 - Axisymetrie
  171. C =====
  172. ELSE IF (IFOU.EQ.0) THEN
  173. NBROBL=4
  174. SEGINI,NOMID
  175. LESOBL(1)='P '
  176. LESOBL(2)='PI '
  177. LESOBL(3)='UR '
  178. LESOBL(4)='RT '
  179. C =====
  180. C 5.3 - Tridimensionnel
  181. C =====
  182. ELSE IF (IFOU.EQ.2) THEN
  183. NBROBL=6
  184. SEGINI,NOMID
  185. LESOBL(1)='P '
  186. LESOBL(2)='PI '
  187. LESOBL(3)='UX '
  188. LESOBL(4)='RY '
  189. LESOBL(5)='UY '
  190. LESOBL(6)='RX '
  191. C =====
  192. C 5.4 - Bidimensionnel PLAN (DP/CP/DPGE)
  193. C =====
  194. ELSE
  195. NBROBL=4
  196. SEGINI,NOMID
  197. LESOBL(1)='P '
  198. LESOBL(2)='PI '
  199. LESOBL(3)='UX '
  200. LESOBL(4)='UY '
  201. ENDIF
  202.  
  203. C 6 - Element de SURFACE LIBRE
  204. C ==============================
  205. ELSE IF (MFR.EQ.23) THEN
  206. NBROBL=3
  207. SEGINI,NOMID
  208. LESOBL(1)='P '
  209. LESOBL(2)='PI '
  210. LESOBL(3)='UZ '
  211.  
  212. C 7 - Element JOINT (35),JOINT CISAILLEMENT (53),JOINT GENERALISE (55)
  213. C COS2 (78)
  214. C =====================================================================
  215. ELSE IF (MFR.EQ.35.OR.MFR.EQ.53.OR.MFR.EQ.55.OR.MFR.EQ.78) THEN
  216. C =====
  217. C 7.1 - Tridimensionnel
  218. C =====
  219. IF (IFOU.EQ.2) THEN
  220. NBROBL=3
  221. SEGINI,NOMID
  222. LESOBL(1)='UX '
  223. LESOBL(2)='UY '
  224. LESOBL(3)='UZ '
  225. C =====
  226. C 7.2 - Bidimensionnel PLAN (CP/DP/DPGE)
  227. C =====
  228. ELSE IF (IFOU.EQ.-1.OR.IFOU.EQ.-2.OR.IFOU.EQ.-3) THEN
  229. NBROBL=2
  230. SEGINI,NOMID
  231. LESOBL(1)='UX '
  232. LESOBL(2)='UY '
  233. C =====
  234. C 7.3 - Axisymetrie
  235. C =====
  236. ELSE IF (IFOU.EQ.0) THEN
  237. NBROBL=2
  238. SEGINI,NOMID
  239. LESOBL(1)='UR '
  240. LESOBL(2)='UZ '
  241. ENDIF
  242.  
  243. C 8 - Elements MASSIFS, de MEMBRANE (25), UNIAXIALE (27)
  244. C ========================================================
  245. ELSE IF (MFR.EQ.1.OR.MFR.EQ.25.OR.MFR.EQ.27.OR.MFR.EQ.31.OR.
  246. $ MFR.EQ.72) THEN
  247. C =====
  248. C 8.1 - Bidimensionnel PLAN (CP/DP)
  249. C =====
  250. IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
  251. NBROBL=2
  252. SEGINI,NOMID
  253. LESOBL(1)='UX '
  254. LESOBL(2)='UY '
  255.  
  256. C =====
  257. C 8.2 - Axisymetrie
  258. C =====
  259. ELSE IF (IFOU.EQ.0) THEN
  260. NBROBL=2
  261. SEGINI,NOMID
  262. LESOBL(1)='UR '
  263. LESOBL(2)='UZ '
  264. C =====
  265. C 8.3 - Fourier
  266. C =====
  267. ELSE IF (IFOU.EQ.1) THEN
  268. NBROBL=3
  269. SEGINI,NOMID
  270. LESOBL(1)='UR '
  271. LESOBL(2)='UZ '
  272. LESOBL(3)='UT '
  273. C =====
  274. C 8.4 - Tridimensionnel
  275. C =====
  276. ELSE IF (IFOU.EQ.2) THEN
  277. NBROBL=3
  278. SEGINI,NOMID
  279. LESOBL(1)='UX '
  280. LESOBL(2)='UY '
  281. LESOBL(3)='UZ '
  282. C =====
  283. C 8.5 - Bidimensionnel PLAN DPGE
  284. C =====
  285. ELSE IF (IFOU.EQ.-3) THEN
  286. NBROBL=5
  287. SEGINI,NOMID
  288. LESOBL(1)='UX '
  289. LESOBL(2)='UY '
  290. LESOBL(3)='UZ '
  291. LESOBL(4)='RY '
  292. LESOBL(5)='RX '
  293. C =====
  294. C 8.6 - Unidimensionnel (1D)
  295. C =====
  296. ELSE IF (IFOU.GE.3.AND.IFOU.LE.15.AND.MFR.EQ.1) THEN
  297. IF (IFOU.LE.6) THEN
  298. NBROBL=1
  299. SEGINI,NOMID
  300. LESOBL(1)='UX '
  301. ELSE IF (IFOU.EQ.7.OR.IFOU.EQ.8) THEN
  302. NBROBL=2
  303. SEGINI,NOMID
  304. LESOBL(1)='UX '
  305. LESOBL(2)='UY '
  306. ELSE IF (IFOU.EQ.9.OR.IFOU.EQ.10) THEN
  307. NBROBL=2
  308. SEGINI,NOMID
  309. LESOBL(1)='UX '
  310. LESOBL(2)='UZ '
  311. ELSE IF (IFOU.EQ.11) THEN
  312. NBROBL=3
  313. SEGINI,NOMID
  314. LESOBL(1)='UX '
  315. LESOBL(2)='UY '
  316. LESOBL(3)='UZ '
  317. ELSE IF (IFOU.EQ.12.OR.IFOU.EQ.13.OR.IFOU.EQ.15) THEN
  318. NBROBL=1
  319. SEGINI,NOMID
  320. LESOBL(1)='UR '
  321. ELSE IF (IFOU.EQ.14) THEN
  322. NBROBL=2
  323. SEGINI,NOMID
  324. LESOBL(1)='UR '
  325. LESOBL(2)='UZ '
  326. ENDIF
  327. ENDIF
  328.  
  329. C 9 - Formulation THERMIQUE
  330. C ===========================
  331. ELSE IF (MFR.EQ.29) THEN
  332. NBROBL=1
  333. SEGINI,NOMID
  334. LESOBL(1)='T '
  335.  
  336. C 10 - Formulation POREUX
  337. C =========================
  338. ELSE IF (MFR.EQ.33) THEN
  339. C ======
  340. C 10.1 - Bidimensionnel PLAN (CP/DP)
  341. C ======
  342. IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
  343. NBROBL=3
  344. SEGINI,NOMID
  345. LESOBL(1)='UX '
  346. LESOBL(2)='UY '
  347. LESOBL(3)='P '
  348. C ======
  349. C 10.2 - Axisymetrie
  350. C ======
  351. ELSE IF (IFOU.EQ.0) THEN
  352. NBROBL=3
  353. SEGINI,NOMID
  354. LESOBL(1)='UR '
  355. LESOBL(2)='UZ '
  356. LESOBL(3)='P '
  357. C ======
  358. C 10.3 - Fourier
  359. C ======
  360. ELSE IF (IFOU.EQ.1) THEN
  361. NBROBL=4
  362. SEGINI,NOMID
  363. LESOBL(1)='UR '
  364. LESOBL(2)='UZ '
  365. LESOBL(3)='UT '
  366. LESOBL(4)='P '
  367. C ======
  368. C 10.4 - Tridimensionnel et bidimensionnel DEFO PLAN GENE
  369. C ======
  370. ELSE IF (IFOU.EQ.2.OR.IFOU.EQ.-3) THEN
  371. NBROBL=4
  372. SEGINI,NOMID
  373. LESOBL(1)='UX '
  374. LESOBL(2)='UY '
  375. LESOBL(3)='UZ '
  376. LESOBL(4)='P '
  377. ENDIF
  378.  
  379. C 11 - Formulation POREUX type Q
  380. C ================================
  381. ELSE IF (MFR.EQ.57) THEN
  382. C ======
  383. C 11.1 - Bidimensionnel PLAN (CP/DP)
  384. C ======
  385. IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
  386. NBROBL=4
  387. SEGINI,NOMID
  388. LESOBL(1)='UX '
  389. LESOBL(2)='UY '
  390. LESOBL(3)='P '
  391. LESOBL(4)='PQ '
  392. C ======
  393. C 11.2 - Axisymetrie
  394. C ======
  395. ELSE IF (IFOU.EQ.0) THEN
  396. NBROBL=4
  397. SEGINI,NOMID
  398. LESOBL(1)='UR '
  399. LESOBL(2)='UZ '
  400. LESOBL(3)='P '
  401. LESOBL(4)='PQ '
  402. C ======
  403. C 11.3 - Fourier
  404. C ======
  405. ELSE IF (IFOU.EQ.1) THEN
  406. NBROBL=5
  407. SEGINI,NOMID
  408. LESOBL(1)='UR '
  409. LESOBL(2)='UZ '
  410. LESOBL(3)='UT '
  411. LESOBL(4)='P '
  412. LESOBL(5)='PQ '
  413. C ======
  414. C 11.4 - Tridimensionnel et bidimensionnel DEFO PLAN GENE
  415. C ======
  416. ELSE IF (IFOU.EQ.2.OR.IFOU.EQ.-3) THEN
  417. NBROBL=5
  418. SEGINI,NOMID
  419. LESOBL(1)='UX '
  420. LESOBL(2)='UY '
  421. LESOBL(3)='UZ '
  422. LESOBL(4)='P '
  423. LESOBL(5)='PQ '
  424. ENDIF
  425.  
  426. C 12 - Formulation POREUX type R
  427. C ================================
  428. ELSE IF (MFR.EQ.59) THEN
  429. C ======
  430. C 12.1 - Bidimensionnel PLAN (CP/DP)
  431. C ======
  432. IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
  433. NBROBL=5
  434. SEGINI,NOMID
  435. LESOBL(1)='UX '
  436. LESOBL(2)='UY '
  437. LESOBL(3)='P '
  438. LESOBL(4)='PQ '
  439. LESOBL(5)='TP '
  440. C ======
  441. C 12.2 - Axisymetrie
  442. C ======
  443. ELSE IF (IFOU.EQ.0) THEN
  444. NBROBL=5
  445. SEGINI,NOMID
  446. LESOBL(1)='UR '
  447. LESOBL(2)='UZ '
  448. LESOBL(3)='P '
  449. LESOBL(4)='PQ '
  450. LESOBL(5)='TP '
  451. C ======
  452. C 12.3 - Fourier
  453. C ======
  454. ELSE IF (IFOU.EQ.1) THEN
  455. NBROBL=6
  456. SEGINI,NOMID
  457. LESOBL(1)='UR '
  458. LESOBL(2)='UZ '
  459. LESOBL(3)='UT '
  460. LESOBL(4)='P '
  461. LESOBL(5)='PQ '
  462. LESOBL(6)='TP '
  463. C ======
  464. C 12.4 - Tridimensionnel et bidimensionnel DEFO PLAN GENE
  465. C ======
  466. ELSE IF(IFOU.EQ.2.OR.IFOU.EQ.-3) THEN
  467. NBROBL=6
  468. SEGINI,NOMID
  469. LESOBL(1)='UX '
  470. LESOBL(2)='UY '
  471. LESOBL(3)='UZ '
  472. LESOBL(4)='P '
  473. LESOBL(5)='PQ '
  474. LESOBL(6)='TP '
  475. ENDIF
  476.  
  477. C 13 - Elements de RACCORD
  478. C ==========================
  479. ELSE IF (MFR.EQ.19.OR.MFR.EQ.21) THEN
  480. NBROBL=2
  481. SEGINI,NOMID
  482. LESOBL(1)='P '
  483. LESOBL(2)='PI '
  484.  
  485. C 14 - Element de RACCORD MASSIF (2e serie de composantes)
  486. C ==========================================================
  487. ELSE IF (MFR.EQ.1019) THEN
  488. C ======
  489. C 14.1 - Bidimensionnel PLAN (CP/DP)
  490. C ======
  491. IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
  492. NBROBL=2
  493. SEGINI,NOMID
  494. LESOBL(1)='UX '
  495. LESOBL(2)='UY '
  496. C ======
  497. C 14.2 - Axisymetrie
  498. C ======
  499. ELSE IF (IFOU.EQ.0) THEN
  500. NBROBL=2
  501. SEGINI,NOMID
  502. LESOBL(1)='UR '
  503. LESOBL(2)='UZ '
  504. C ======
  505. C 14.3 - Fourier
  506. C ======
  507. ELSE IF (IFOU.EQ.1) THEN
  508. NBROBL=3
  509. SEGINI,NOMID
  510. LESOBL(1)='UR '
  511. LESOBL(2)='UZ '
  512. LESOBL(3)='UT '
  513. C ======
  514. C 14.4 - Tridimensionnel et bidimensionnel DEFO PLAN GENE
  515. C ======
  516. ELSE IF (IFOU.EQ.2.OR.IFOU.EQ.-3) THEN
  517. NBROBL=3
  518. SEGINI,NOMID
  519. LESOBL(1)='UX '
  520. LESOBL(2)='UY '
  521. LESOBL(3)='UZ '
  522. ENDIF
  523.  
  524. C 15 - Element de RACCORD COQUE (2e serie de composantes)
  525. C =========================================================
  526. ELSE IF (MFR.EQ.1021) THEN
  527. C ======
  528. C 15.1 - Bidimensionnel PLAN (CP/DP)
  529. C ======
  530. IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
  531. NBROBL=3
  532. SEGINI,NOMID
  533. LESOBL(1)='UX '
  534. LESOBL(2)='UY '
  535. LESOBL(3)='RZ '
  536. C ======
  537. C 15.2 - Axisymetrie
  538. C ======
  539. ELSE IF (IFOU.EQ.0) THEN
  540. NBROBL=3
  541. SEGINI,NOMID
  542. LESOBL(1)='UR '
  543. LESOBL(2)='UZ '
  544. LESOBL(3)='RT '
  545. C ======
  546. C 15.3 - Fourier
  547. C ======
  548. ELSE IF (IFOU.EQ.1) THEN
  549. NBROBL=4
  550. SEGINI,NOMID
  551. LESOBL(1)='UR '
  552. LESOBL(2)='UZ '
  553. LESOBL(3)='UT '
  554. LESOBL(4)='RT '
  555. C ======
  556. C 15.4 - Tridimensionnel et bidimensionnel DEFO PLAN GENE
  557. C ======
  558. ELSE IF (IFOU.EQ.2.OR.IFOU.EQ.-3) THEN
  559. NBROBL=6
  560. SEGINI,NOMID
  561. LESOBL(1)='UX '
  562. LESOBL(2)='UY '
  563. LESOBL(3)='UZ '
  564. LESOBL(4)='RX '
  565. LESOBL(5)='RY '
  566. LESOBL(6)='RZ '
  567. ENDIF
  568.  
  569. C 16 - Element TUYO
  570. C ===================
  571. ELSE IF (MFR.EQ.39) THEN
  572. NBROBL=6+9
  573. SEGINI,NOMID
  574. LESOBL(1)='UX '
  575. LESOBL(2)='UY '
  576. LESOBL(3)='UZ '
  577. LESOBL(4)='RX '
  578. LESOBL(5)='RY '
  579. LESOBL(6)='RZ '
  580. LESOBL(7)='W0 '
  581. LESOBL(8)='U1 '
  582. LESOBL(9)='U2 '
  583. LESOBL(10)='U3 '
  584. LESOBL(11)='U4 '
  585. LESOBL(12)='W1 '
  586. LESOBL(13)='W2 '
  587. LESOBL(14)='W3 '
  588. LESOBL(15)='W4 '
  589. C 17 - Element POI1 materiau MODAL
  590. C ===================
  591. ELSE IF (MFR.EQ.26) THEN
  592. NBROBL=1
  593. SEGINI,NOMID
  594. LESOBL(1)='ALFA'
  595. C 18 - Element POI1 materiau STATIQUE
  596. C ===================
  597. ELSE IF (MFR.EQ.28) THEN
  598. NBROBL=1
  599. SEGINI,NOMID
  600. LESOBL(1)='BETA'
  601. C 19 - Element XFEM (xfem meca rupture)
  602. C ===================
  603. ELSEIF(MFR.EQ.63) THEN
  604. CTY - Element XQ4R (xfem meca rupture en 2D)
  605. IF(IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
  606. NBROBL=2
  607. NBRFAC=18
  608. SEGINI,NOMID
  609. LESOBL(1)='UX '
  610. LESOBL(2)='UY '
  611. LESFAC(1)='AX '
  612. LESFAC(2)='AY '
  613. LESFAC(3)='B1X '
  614. LESFAC(4)='B1Y '
  615. LESFAC(5)='C1X '
  616. LESFAC(6)='C1Y '
  617. LESFAC(7)='D1X '
  618. LESFAC(8)='D1Y '
  619. LESFAC(9)='E1X '
  620. LESFAC(10)='E1Y '
  621. LESFAC(11)='B2X '
  622. LESFAC(12)='B2Y '
  623. LESFAC(13)='C2X '
  624. LESFAC(14)='C2Y '
  625. LESFAC(15)='D2X '
  626. LESFAC(16)='D2Y '
  627. LESFAC(17)='E2X '
  628. LESFAC(18)='E2Y '
  629. CTY - Element XC8R (xfem meca rupture en 3D)
  630. ELSE IF (IFOU.EQ.2) THEN
  631. NBROBL=3
  632. NBRFAC=27
  633. SEGINI,NOMID
  634. LESOBL(1)='UX '
  635. LESOBL(2)='UY '
  636. LESOBL(3)='UZ '
  637. LESFAC(1)='AX '
  638. LESFAC(2)='AY '
  639. LESFAC(3)='AZ '
  640. LESFAC(4)='B1X '
  641. LESFAC(5)='B1Y '
  642. LESFAC(6)='B1Z '
  643. LESFAC(7)='C1X '
  644. LESFAC(8)='C1Y '
  645. LESFAC(9)='C1Z '
  646. LESFAC(10)='D1X '
  647. LESFAC(11)='D1Y '
  648. LESFAC(12)='D1Z '
  649. LESFAC(13)='E1X '
  650. LESFAC(14)='E1Y '
  651. LESFAC(15)='E1Z '
  652. LESFAC(16)='B2X '
  653. LESFAC(17)='B2Y '
  654. LESFAC(18)='B2Z '
  655. LESFAC(19)='C2X '
  656. LESFAC(20)='C2Y '
  657. LESFAC(21)='C2Z '
  658. LESFAC(22)='D2X '
  659. LESFAC(23)='D2Y '
  660. LESFAC(24)='D2Z '
  661. LESFAC(25)='E2X '
  662. LESFAC(26)='E2Y '
  663. LESFAC(27)='E2Z '
  664. ENDIF
  665.  
  666. C 20 - Elements de zones cohesives
  667. C ==================================
  668. ELSEIF(MFR.EQ.77) THEN
  669. C - Element ZCO2 (xfem meca rupture en 2D)
  670. IF(IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
  671. NBROBL=2
  672. NBRFAC=18
  673. SEGINI,NOMID
  674. LESOBL(1)='AX '
  675. LESOBL(2)='AY '
  676. LESFAC(1)='UX '
  677. LESFAC(2)='UY '
  678. LESFAC(3)='B1X '
  679. LESFAC(4)='B1Y '
  680. LESFAC(5)='C1X '
  681. LESFAC(6)='C1Y '
  682. LESFAC(7)='D1X '
  683. LESFAC(8)='D1Y '
  684. LESFAC(9)='E1X '
  685. LESFAC(10)='E1Y '
  686. LESFAC(11)='B2X '
  687. LESFAC(12)='B2Y '
  688. LESFAC(13)='C2X '
  689. LESFAC(14)='C2Y '
  690. LESFAC(15)='D2X '
  691. LESFAC(16)='D2Y '
  692. LESFAC(17)='E2X '
  693. LESFAC(18)='E2Y '
  694. C - Elements ZCO3 et ZCO4(xfem meca rupture en 3D)
  695. ELSE IF (IFOU.EQ.2) THEN
  696. NBROBL=3
  697. NBRFAC=27
  698. SEGINI,NOMID
  699. LESOBL(1)='AX '
  700. LESOBL(2)='AY '
  701. LESOBL(3)='AZ '
  702. LESFAC(1)='UX '
  703. LESFAC(2)='UY '
  704. LESFAC(3)='UZ '
  705. LESFAC(4)='B1X '
  706. LESFAC(5)='B1Y '
  707. LESFAC(6)='B1Z '
  708. LESFAC(7)='C1X '
  709. LESFAC(8)='C1Y '
  710. LESFAC(9)='C1Z '
  711. LESFAC(10)='D1X '
  712. LESFAC(11)='D1Y '
  713. LESFAC(12)='D1Z '
  714. LESFAC(13)='E1X '
  715. LESFAC(14)='E1Y '
  716. LESFAC(15)='E1Z '
  717. LESFAC(16)='B2X '
  718. LESFAC(17)='B2Y '
  719. LESFAC(18)='B2Z '
  720. LESFAC(19)='C2X '
  721. LESFAC(20)='C2Y '
  722. LESFAC(21)='C2Z '
  723. LESFAC(22)='D2X '
  724. LESFAC(23)='D2Y '
  725. LESFAC(24)='D2Z '
  726. LESFAC(25)='E2X '
  727. LESFAC(26)='E2Y '
  728. LESFAC(27)='E2Z '
  729. ENDIF
  730.  
  731. C 21 - Formulation ELECTROSTATIQUE (base MASSIF)
  732. C ==================================
  733. ELSE IF (MFR.EQ.71) THEN
  734. NBROBL=1
  735. SEGINI,NOMID
  736. LESOBL(1)='VEL '
  737.  
  738. C 22 - Formulation DIFFUSION (base MASSIF)
  739. C ============================
  740. ELSE IF (MFR.EQ.73) THEN
  741. C* L'INCOnnue (DDL) contenue dans TYMODE(1) est ajoutee ailleurs
  742. NBROBL=1
  743. SEGINI,NOMID
  744. LESOBL(1)='.. '
  745. ENDIF
  746.  
  747. C Par DEFAUT : segment VIDE
  748. C ===========================
  749. IF (NOMID.EQ.0) THEN
  750. SEGINI,NOMID
  751. ELSE
  752. if (ifomod.eq.6) then
  753. nbrfa0 = nbrfac
  754. NBRFAC = NBROBL + (nbrfa0*2)
  755. segadj nomid
  756. do imo = 1,nbrobl
  757. lesfac(nbrfa0 + imo)(2:4) = lesobl(imo)(1:3)
  758. lesfac(nbrfa0 + imo)(1:1) = 'I'
  759. enddo
  760. do imo = 1,nbrfa0
  761. lesfac(nbrfa0+nbrobl+imo)(2:4) = lesfac(imo)(1:3)
  762. lesfac(nbrfa0+nbrobl+imo)(1:1) = 'I'
  763. enddo
  764. endif
  765. ENDIF
  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