Télécharger idforc.eso

Retour à la liste

Numérotation des lignes :

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

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