Télécharger idgrad.eso

Retour à la liste

Numérotation des lignes :

  1. C IDGRAD SOURCE FANDEUR 16/12/05 21:39:42 9211
  2.  
  3. C=======================================================================
  4. C= RECHERCHE DES NOMS DE COMPOSANTES DE GRADIENTS =
  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 IDGRAD (MFR,IFOU,IPNOMC,NBROBL,NBRFAC)
  19.  
  20. IMPLICIT INTEGER(I-N)
  21. IMPLICIT REAL*8(A-H,O-Z)
  22.  
  23. -INC SMMODEL
  24.  
  25. IPNOMC=0
  26. NBROBL=0
  27. NBRFAC=0
  28. NOMID =0
  29.  
  30. C FORMULATION THERMOHYDRIQUE
  31. C ==========================
  32. IF (MFR.EQ.65) THEN
  33. C =====
  34. C Tridimensionnel
  35. C =====
  36. IF (IFOU.EQ.2) THEN
  37. NBROBL=9
  38. SEGINI,NOMID
  39. LESOBL(1)='PG,X '
  40. LESOBL(2)='PG,Y '
  41. LESOBL(3)='PG,Z '
  42. LESOBL(4)='PC,X '
  43. LESOBL(5)='PC,Y '
  44. LESOBL(6)='PC,Z '
  45. LESOBL(7)='T,X '
  46. LESOBL(8)='T,Y '
  47. LESOBL(9)='T,Z '
  48. C =====
  49. C Bidimensionnel PLAN (DP/CP/DPGE)
  50. C =====
  51. ELSE IF (IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
  52. NBROBL=6
  53. SEGINI,NOMID
  54. LESOBL(1)='PG,X '
  55. LESOBL(2)='PG,Y '
  56. LESOBL(3)='PC,X '
  57. LESOBL(4)='PC,Y '
  58. LESOBL(5)='T,X '
  59. LESOBL(6)='T,Y '
  60. C =====
  61. C Axisymetrie et Fourier
  62. C =====
  63. ELSE IF (IFOU.EQ.0.OR.IFOU.EQ.1) THEN
  64. NBROBL=6
  65. SEGINI,NOMID
  66. LESOBL(1)='PG,R '
  67. LESOBL(2)='PG,Z '
  68. LESOBL(3)='PC,R '
  69. LESOBL(4)='PC,Z '
  70. LESOBL(5)='T,R '
  71. LESOBL(6)='T,Z '
  72. ENDIF
  73.  
  74. C 1A- Element JOINT JOI1
  75. C =======================
  76.  
  77. ELSEIF (MFR.EQ.75) THEN
  78. C =====
  79. C 1A.1 - Tridimensionnel et Bidimensionnel PLAN DPGE
  80. C =====
  81. IF (IFOU.EQ.2 .OR. IFOU.EQ.-3) THEN
  82. NBROBL=9
  83. SEGINI,NOMID
  84. LESOBL(1)='UX,X '
  85. LESOBL(2)='UX,Y '
  86. LESOBL(3)='UX,Z '
  87. LESOBL(4)='UY,X '
  88. LESOBL(5)='UY,Y '
  89. LESOBL(6)='UY,Z '
  90. LESOBL(7)='UZ,X '
  91. LESOBL(8)='UZ,Y '
  92. LESOBL(9)='UZ,Z '
  93. C =====
  94. C 1A.2 - Bidimensionnel PLAN (CP/DP) et Axisymetrie
  95. C =====
  96. ELSE IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
  97. NBROBL=4
  98. SEGINI,NOMID
  99. LESOBL(1)='UX,X '
  100. LESOBL(2)='UX,Y '
  101. LESOBL(3)='UY,X '
  102. LESOBL(4)='UY,Y '
  103. C =====
  104. C 1A.3 - Axisymetrie
  105. C =====
  106. ELSE IF (IFOU.EQ.0) THEN
  107. NBROBL=4
  108. SEGINI,NOMID
  109. LESOBL(1)='UR,R '
  110. LESOBL(2)='UR,Z '
  111. LESOBL(3)='UZ,R '
  112. LESOBL(4)='UZ,Z '
  113. ENDIF
  114.  
  115. C 1B- Element JOINT JOI2
  116. C =======================
  117.  
  118. ELSEIF (MFR.EQ.35) THEN
  119. C =====
  120. C 1B.1 - Tridimensionnel et Bidimensionnel PLAN DPGE
  121. C =====
  122. IF (IFOU.EQ.2 .OR. IFOU.EQ.-3) THEN
  123. NBROBL=3
  124. SEGINI,NOMID
  125. LESOBL(1)='U1,N '
  126. LESOBL(2)='U2,N '
  127. LESOBL(3)='UN,N '
  128. C =====
  129. C 1A.2 - Bidimensionnel PLAN (CP/DP) et Axisymetrie
  130. C =====
  131. ELSE IF (IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.0) THEN
  132. NBROBL=2
  133. SEGINI,NOMID
  134. LESOBL(1)='US,N '
  135. LESOBL(2)='UN,N '
  136. ENDIF
  137.  
  138. C 2 - Elements MASSIFS (standard, incompressibles et XFEM)
  139. C ======================
  140. ELSE IF (MFR.EQ.1 .OR. MFR.EQ.31 .OR. MFR.EQ.63) THEN
  141. C =====
  142. C 2.1 - Bidimensionnel PLAN (CP/DP/DPGE) et Tridimensionnel
  143. C =====
  144. IF (IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.IFOU.EQ.-1.OR.IFOU.EQ.2) THEN
  145. NBROBL=9
  146. SEGINI,NOMID
  147. LESOBL(1)='UX,X '
  148. LESOBL(2)='UX,Y '
  149. LESOBL(3)='UX,Z '
  150. LESOBL(4)='UY,X '
  151. LESOBL(5)='UY,Y '
  152. LESOBL(6)='UY,Z '
  153. LESOBL(7)='UZ,X '
  154. LESOBL(8)='UZ,Y '
  155. LESOBL(9)='UZ,Z '
  156. C =====
  157. C 2.2 - 2D Axisymetrie et Fourier
  158. C =====
  159. ELSE IF (IFOU.EQ.0.OR.IFOU.EQ.1) THEN
  160. NBROBL=9
  161. SEGINI,NOMID
  162. LESOBL(1)='UR,R '
  163. LESOBL(2)='UR,Z '
  164. LESOBL(3)='UR,T '
  165. LESOBL(4)='UZ,R '
  166. LESOBL(5)='UZ,Z '
  167. LESOBL(6)='UZ,T '
  168. LESOBL(7)='UT,R '
  169. LESOBL(8)='UT,Z '
  170. LESOBL(9)='UT,T '
  171. C =====
  172. C 2.3 - Massif Unidimensionnel (1D) PLAN et AXIS/SPHE
  173. C =====
  174. ELSE IF (MFR.EQ.1 .AND. (IFOU.GE.3.AND.IFOU.LE.15)) THEN
  175. NBROBL=9
  176. SEGINI,NOMID
  177. IF (IFOU.GE.3.AND.IFOU.LE.11) THEN
  178. LESOBL(1)='UX,X '
  179. LESOBL(2)='UX,Y '
  180. LESOBL(3)='UX,Z '
  181. LESOBL(4)='UY,X '
  182. LESOBL(5)='UY,Y '
  183. LESOBL(6)='UY,Z '
  184. LESOBL(7)='UZ,X '
  185. LESOBL(8)='UZ,Y '
  186. LESOBL(9)='UZ,Z '
  187. ELSE
  188. LESOBL(1)='UR,R '
  189. LESOBL(2)='UR,Z '
  190. LESOBL(3)='UR,T '
  191. LESOBL(4)='UZ,R '
  192. LESOBL(5)='UZ,Z '
  193. LESOBL(6)='UZ,T '
  194. LESOBL(7)='UT,R '
  195. LESOBL(8)='UT,Z '
  196. LESOBL(9)='UT,T '
  197. ENDIF
  198. ENDIF
  199.  
  200. C 3 - Elements COQUES
  201. C =====================
  202. ELSE IF (MFR.EQ.3) THEN
  203. C =====
  204. C 3.1 - Tridimensionnel
  205. C =====
  206. IF (IFOU.EQ.2) THEN
  207. NBROBL=9
  208. SEGINI,NOMID
  209. LESOBL(1)='UX,X'
  210. LESOBL(2)='UX,Y'
  211. LESOBL(3)='UX,Z'
  212. LESOBL(4)='UY,X'
  213. LESOBL(5)='UY,Y'
  214. LESOBL(6)='UY,Z'
  215. LESOBL(7)='UZ,X'
  216. LESOBL(8)='UZ,Y'
  217. LESOBL(9)='UZ,Z'
  218. C =====
  219. C 3.2 - Axisymetrie, Fourier, Bidimensionnel PLAN (CP/DP)
  220. C =====
  221. ELSE IF (IFOU.EQ.1.OR.IFOU.EQ.0.OR.
  222. . IFOU.EQ.-1.OR.IFOU.EQ.-2.or.ifou.eq.-3) THEN
  223. NBROBL=9
  224. SEGINI,NOMID
  225. LESOBL(1)='UU,S'
  226. LESOBL(2)='UU,T'
  227. LESOBL(3)='UU,N'
  228. LESOBL(4)='UV,S'
  229. LESOBL(5)='UV,T'
  230. LESOBL(6)='UV,N'
  231. LESOBL(7)='UW,S'
  232. LESOBL(8)='UW,T'
  233. LESOBL(9)='UW,N'
  234. ENDIF
  235.  
  236. C 4 - Elements COQUES EPAISSES (COQ6,COQ8)
  237. C ==============================
  238. ELSE IF (MFR.EQ.5) THEN
  239. C =====
  240. C 4.1 - Tridimensionnel
  241. C =====
  242. IF (IFOU.EQ.2) THEN
  243. NBROBL=9
  244. SEGINI,NOMID
  245. LESOBL(1)='UX,X'
  246. LESOBL(2)='UX,Y'
  247. LESOBL(3)='UX,Z'
  248. LESOBL(4)='UY,X'
  249. LESOBL(5)='UY,Y'
  250. LESOBL(6)='UY,Z'
  251. LESOBL(7)='UZ,X'
  252. LESOBL(8)='UZ,Y'
  253. LESOBL(9)='UZ,Z'
  254. C =====
  255. C 4.2 - Axisymetrie et Fourier
  256. C =====
  257. ELSE IF (IFOU.EQ.1.OR.IFOU.EQ.0) THEN
  258. NBROBL=9
  259. SEGINI,NOMID
  260. LESOBL(1)='UR,R'
  261. LESOBL(2)='UR,Z'
  262. LESOBL(3)='UR,T'
  263. LESOBL(4)='UZ,R'
  264. LESOBL(5)='UZ,Z'
  265. LESOBL(6)='UZ,T'
  266. LESOBL(7)='UT,R'
  267. LESOBL(8)='UT,Z'
  268. LESOBL(9)='UT,T'
  269. ENDIF
  270.  
  271. C 5 - Element de COQUE en CISAILLEMENT TRANSVERSE (tridimensionnel)
  272. C =================================================
  273. ELSE IF (MFR.EQ.9) THEN
  274. IF (IFOU.EQ.2) THEN
  275. NBROBL=9
  276. SEGINI,NOMID
  277. LESOBL(1)='UX,X'
  278. LESOBL(2)='UX,Y'
  279. LESOBL(3)='UX,Z'
  280. LESOBL(4)='UY,X'
  281. LESOBL(5)='UY,Y'
  282. LESOBL(6)='UY,Z'
  283. LESOBL(7)='UZ,X'
  284. LESOBL(8)='UZ,Y'
  285. LESOBL(9)='UZ,Z'
  286. ENDIF
  287.  
  288. C 6 - Formulation THERMIQUE MASSIF
  289. C ==================================
  290. ELSE IF (MFR.EQ.29) THEN
  291. C =====
  292. C 6.1 - Tridimensionnel
  293. C =====
  294. IF (IFOU.EQ.2) THEN
  295. NBROBL=3
  296. SEGINI,NOMID
  297. LESOBL(1)='T,X '
  298. LESOBL(2)='T,Y '
  299. LESOBL(3)='T,Z '
  300. C =====
  301. C 6.2 - Bidimensionnel PLAN (DP/CP/DPGE)
  302. C =====
  303. ELSE IF (IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
  304. NBROBL=2
  305. SEGINI,NOMID
  306. LESOBL(1)='T,X '
  307. LESOBL(2)='T,Y '
  308. C =====
  309. C 6.3 - Axisymetrie et Fourier
  310. C =====
  311. ELSE IF (IFOU.EQ.0.OR.IFOU.EQ.1) THEN
  312. NBROBL=2
  313. SEGINI,NOMID
  314. LESOBL(1)='T,R '
  315. LESOBL(2)='T,Z '
  316. C =====
  317. C 6.4 - Unidimensionnel (1D)
  318. C =====
  319. ELSE IF (IFOU.GE.3.AND.IFOU.LE.11) THEN
  320. NBROBL=1
  321. SEGINI,NOMID
  322. LESOBL(1)='T,X '
  323. ELSE IF (IFOU.GE.12.AND.IFOU.LE.15) THEN
  324. NBROBL=1
  325. SEGINI,NOMID
  326. LESOBL(1)='T,R '
  327. ENDIF
  328.  
  329. C 7 - Formulation THERMIQUE COQUE
  330. C =================================
  331. ELSE IF (MFR.EQ.39) THEN
  332. C =====
  333. C 7.1 - Tridimensionnel
  334. C =====
  335. IF (IFOU.EQ.2) THEN
  336. NBROBL=2
  337. SEGINI,NOMID
  338. LESOBL(1)='T,X '
  339. LESOBL(2)='T,Y '
  340. C =====
  341. C 7.2 - Axisymetrie et Fourier
  342. C =====
  343. ELSE IF (IFOU.EQ.0) THEN
  344. NBROBL=1
  345. SEGINI,NOMID
  346. LESOBL(1)='T,Z '
  347. ENDIF
  348.  
  349. C 8 - Formulation MAGNETODYNAMIQUE COQUE (tridimensionnel)
  350. C ========================================
  351. ELSE IF (MFR.EQ.69) THEN
  352. NBROBL=3
  353. SEGINI,NOMID
  354. LESOBL(1)='FC,X '
  355. LESOBL(2)='FC,Y '
  356. LESOBL(3)='FC,Z '
  357.  
  358. C 9 - formulation XFEM meca rupt
  359. C =====================================
  360. c ELSE IF (MFR.EQ.63) THEN
  361. CTY - Element XQ4R (xfem meca rupture en 2D)
  362. c IF(IFOU.EQ.-2.OR.ifou.EQ.-1) then
  363. c NBROBL=4
  364. c NBRFAC=36
  365. c NBRFAC=20
  366. c SEGINI,NOMID
  367. c LESOBL(1)='UX,X'
  368. c LESOBL(2)='UX,Y'
  369. c LESOBL(3)='UY,X'
  370. c LESOBL(4)='UY,Y'
  371. c LESFAC(1)='AX,X'
  372. c LESFAC(2)='AX,Y'
  373. c LESFAC(3)='AY,X'
  374. c LESFAC(4)='AY,Y'
  375. c LESFAC(5)='B1XX'
  376. c LESFAC(6)='B1XY'
  377. c LESFAC(7)='B1YX'
  378. c LESFAC(8)='B1YY'
  379. c LESFAC(9)='C1XX'
  380. c LESFAC(10)='C1XY'
  381. c LESFAC(11)='C1YX'
  382. c LESFAC(12)='C1YY'
  383. c LESFAC(13)='D1XX'
  384. c LESFAC(14)='D1XY'
  385. c LESFAC(15)='D1YX'
  386. c LESFAC(16)='D1YY'
  387. c LESFAC(17)='E1XX'
  388. c LESFAC(18)='E1XY'
  389. c LESFAC(19)='E1YX'
  390. c LESFAC(20)='E1YY'
  391. c LESFAC(21)='B2XX'
  392. c LESFAC(22)='B2XY'
  393. c LESFAC(23)='B2YX'
  394. c LESFAC(24)='B2YY'
  395. c LESFAC(25)='C2XX'
  396. c LESFAC(26)='C2XY'
  397. c LESFAC(27)='C2YX'
  398. c LESFAC(28)='C2YY'
  399. c LESFAC(29)='D2XX'
  400. c LESFAC(30)='D2XY'
  401. c LESFAC(31)='D2YX'
  402. c LESFAC(32)='D2YY'
  403. c LESFAC(33)='E2XX'
  404. c LESFAC(34)='E2XY'
  405. c LESFAC(35)='E2YX'
  406. c LESFAC(36)='E2YY'
  407. c ELSE IF (IFOU.EQ.2) THEN
  408. CTY - Element XC8R (xfem meca rupture en 3D)
  409. c NBROBL=9
  410. c NBRFAC=45
  411. c NBRFAC=81
  412. c SEGINI,NOMID
  413. c LESOBL(1)='UX,X'
  414. c LESOBL(2)='UX,Y'
  415. c LESOBL(3)='UX,Z'
  416. c LESOBL(4)='UY,X'
  417. c LESOBL(5)='UY,Y'
  418. c LESOBL(6)='UY,Z'
  419. c LESOBL(7)='UZ,X'
  420. c LESOBL(8)='UZ,Y'
  421. c LESOBL(9)='UZ,Z'
  422. c LESFAC(1)='AX,X'
  423. c LESFAC(2)='AX,Y'
  424. c LESFAC(3)='AX,Z'
  425. c LESFAC(4)='AY,X'
  426. c LESFAC(5)='AY,Y'
  427. c LESFAC(6)='AY,Z'
  428. c LESFAC(7)='AZ,X'
  429. c LESFAC(8)='AZ,Y'
  430. c LESFAC(9)='AZ,Z'
  431. c LESFAC(10)='B1XX'
  432. c LESFAC(11)='B1XY'
  433. c LESFAC(12)='B1XZ'
  434. c LESFAC(13)='B1YX'
  435. c LESFAC(14)='B1YY'
  436. c LESFAC(15)='B1YZ'
  437. c LESFAC(16)='B1ZX'
  438. c LESFAC(17)='B1ZY'
  439. c LESFAC(18)='B1ZZ'
  440. c LESFAC(19)='C1XX'
  441. c LESFAC(20)='C1XY'
  442. c LESFAC(21)='C1XZ'
  443. c LESFAC(22)='C1YX'
  444. c LESFAC(23)='C1YY'
  445. c LESFAC(24)='C1YZ'
  446. c LESFAC(25)='C1ZX'
  447. c LESFAC(26)='C1ZY'
  448. c LESFAC(27)='C1ZZ'
  449. c LESFAC(28)='D1XX'
  450. c LESFAC(29)='D1XY'
  451. c LESFAC(30)='D1XZ'
  452. c LESFAC(31)='D1YX'
  453. c LESFAC(32)='D1YY'
  454. c LESFAC(33)='D1YZ'
  455. c LESFAC(34)='D1ZX'
  456. c LESFAC(35)='D1ZY'
  457. c LESFAC(36)='D1ZZ'
  458. c LESFAC(37)='E1XX'
  459. c LESFAC(38)='E1XY'
  460. c LESFAC(39)='E1XZ'
  461. c LESFAC(40)='E1YX'
  462. c LESFAC(41)='E1YY'
  463. c LESFAC(42)='E1YZ'
  464. c LESFAC(43)='E1ZX'
  465. c LESFAC(44)='E1ZY'
  466. c LESFAC(45)='E1ZZ'
  467. c LESFAC(46)='B2XX'
  468. c LESFAC(47)='B2XY'
  469. c LESFAC(48)='B2XZ'
  470. c LESFAC(49)='B2YX'
  471. c LESFAC(50)='B2YY'
  472. c LESFAC(51)='B2YZ'
  473. c LESFAC(52)='B2ZX'
  474. c LESFAC(53)='B2ZY'
  475. c LESFAC(54)='B2ZZ'
  476. c LESFAC(55)='C2XX'
  477. c LESFAC(56)='C2XY'
  478. c LESFAC(57)='C2XZ'
  479. c LESFAC(58)='C2YX'
  480. c LESFAC(59)='C2YY'
  481. c LESFAC(60)='C2YZ'
  482. c LESFAC(61)='C2ZX'
  483. c LESFAC(62)='C2ZY'
  484. c LESFAC(63)='C2ZZ'
  485. c LESFAC(64)='D2XX'
  486. c LESFAC(65)='D2XY'
  487. c LESFAC(66)='D2XZ'
  488. c LESFAC(67)='D2YX'
  489. c LESFAC(68)='D2YY'
  490. c LESFAC(69)='D2YZ'
  491. c LESFAC(70)='D2ZX'
  492. c LESFAC(71)='D2ZY'
  493. c LESFAC(72)='D2ZZ'
  494. c LESFAC(73)='E2XX'
  495. c LESFAC(74)='E2XY'
  496. c LESFAC(75)='E2XZ'
  497. c LESFAC(76)='E2YX'
  498. c LESFAC(77)='E2YY'
  499. c LESFAC(78)='E2YZ'
  500. c LESFAC(79)='E2ZX'
  501. c LESFAC(80)='E2ZY'
  502. c LESFAC(81)='E2ZZ'
  503. c ENDIF
  504.  
  505. C 10- Formulation MECANIQUE POUTRES ET TUYAUX
  506. C ============================================
  507. ELSE IF (MFR.EQ.7.OR.MFR.EQ.13) THEN
  508. C ======
  509. C 10.1 - Tridimensionnel
  510. C ======
  511. IF (IFOU.EQ.2) THEN
  512. NBROBL=3
  513. SEGINI,NOMID
  514. LESOBL(1)='U,S '
  515. LESOBL(2)='V,S '
  516. LESOBL(3)='W,S '
  517. C ======
  518. C 10.2 - Autres cas
  519. C ======
  520. ELSE
  521. NBROBL=2
  522. SEGINI,NOMID
  523. LESOBL(1)='U,S '
  524. LESOBL(2)='V,S '
  525. ENDIF
  526.  
  527. C 11 - Formulation ELECTROSTATIQUE (base MASSIF)
  528. C ==================================
  529. ELSE IF (MFR.EQ.71) THEN
  530. C ======
  531. C 11.1 - Tridimensionnel
  532. C ======
  533. IF (IFOU.EQ.2) THEN
  534. NBROBL=3
  535. SEGINI,NOMID
  536. LESOBL(1)='VE,X '
  537. LESOBL(2)='VE,Y '
  538. LESOBL(3)='VE,Z '
  539. C ======
  540. C 11.2 - Bidimensionnel PLAN (DP/CP/DPGE)
  541. C ======
  542. ELSE IF (IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
  543. NBROBL=2
  544. SEGINI,NOMID
  545. LESOBL(1)='VE,X '
  546. LESOBL(2)='VE,Y '
  547. C ======
  548. C 11.3 - Axisymetrie et Fourier
  549. C ======
  550. ELSE IF (IFOU.EQ.0.OR.IFOU.EQ.1) THEN
  551. NBROBL=2
  552. SEGINI,NOMID
  553. LESOBL(1)='VE,R '
  554. LESOBL(2)='VE,Z '
  555. C ======
  556. C 11.4 - Unidimensionnel (1D)
  557. C ======
  558. ELSE IF (IFOU.GE.3.AND.IFOU.LE.11) THEN
  559. NBROBL=1
  560. SEGINI,NOMID
  561. LESOBL(1)='VE,X '
  562. ELSE IF (IFOU.GE.12.AND.IFOU.LE.15) THEN
  563. NBROBL=1
  564. SEGINI,NOMID
  565. LESOBL(1)='VE,R '
  566. ENDIF
  567.  
  568. C 12 - Formulation DIFFUSION (base MASSIF)
  569. C ============================
  570. ELSE IF (MFR.EQ.73) THEN
  571. C* L'INCOnnue (DDL) contenue dans TYMODE(1) est ajoutee ailleurs
  572. C* en tete des (noms de) comppsantes (apres l'appel a idgrad)
  573. C ======
  574. C 12.1 - Tridimensionnel
  575. C ======
  576. IF (IFOU.EQ.2) THEN
  577. NBROBL=3
  578. SEGINI,NOMID
  579. LESOBL(1)=',X '
  580. LESOBL(2)=',Y '
  581. LESOBL(3)=',Z '
  582. C ======
  583. C 12.2 - Bidimensionnel PLAN (DP/CP/DPGE)
  584. C ======
  585. ELSE IF (IFOU.EQ.-3.OR.IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
  586. NBROBL=2
  587. SEGINI,NOMID
  588. LESOBL(1)=',X '
  589. LESOBL(2)=',Y '
  590. C ======
  591. C 12.3 - Axisymetrie et Fourier
  592. C ======
  593. ELSE IF (IFOU.EQ.0.OR.IFOU.EQ.1) THEN
  594. NBROBL=2
  595. SEGINI,NOMID
  596. LESOBL(1)=',R '
  597. LESOBL(2)=',Z '
  598. C ======
  599. C 12.4 - Unidimensionnel (1D)
  600. C ======
  601. ELSE IF (IFOU.GE.3.AND.IFOU.LE.11) THEN
  602. NBROBL=1
  603. SEGINI,NOMID
  604. LESOBL(1)=',X '
  605. ELSE IF (IFOU.GE.12.AND.IFOU.LE.15) THEN
  606. NBROBL=1
  607. SEGINI,NOMID
  608. LESOBL(1)=',R '
  609. ENDIF
  610.  
  611. C 13 - Formulation ZONE COHESIVE
  612. C ================================
  613. ELSE IF (MFR.EQ.77) THEN
  614. C =====
  615. C 13.1 - Tridimensionnel
  616. C =====
  617. IF (IFOU.EQ.2) THEN
  618. NBROBL=9
  619. SEGINI,NOMID
  620. LESOBL(1)='AX,X '
  621. LESOBL(2)='AX,Y '
  622. LESOBL(3)='AX,Z '
  623. LESOBL(4)='AY,X '
  624. LESOBL(5)='AY,Y '
  625. LESOBL(6)='AY,Z '
  626. LESOBL(7)='AZ,X '
  627. LESOBL(8)='AZ,Y '
  628. LESOBL(9)='AZ,Z '
  629. C =====
  630. C 13.2 - Bidimensionnel PLAN (CP/DP/DPGE)
  631. C =====
  632. ELSE IF (IFOU.EQ.-2.OR.IFOU.EQ.-1) THEN
  633. NBROBL=4
  634. SEGINI,NOMID
  635. LESOBL(1)='AX,X '
  636. LESOBL(2)='AX,Y '
  637. LESOBL(3)='AY,X '
  638. LESOBL(4)='AY,Y '
  639. ENDIF
  640.  
  641. ENDIF
  642.  
  643. C Par DEFAUT : segment VIDE
  644. C ===========================
  645. IF (NOMID.EQ.0) SEGINI,NOMID
  646.  
  647. SEGDES,NOMID
  648. IPNOMC=NOMID
  649.  
  650. RETURN
  651. END
  652.  
  653.  
  654.  
  655.  
  656.  

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