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

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