Télécharger elquoi.eso

Retour à la liste

Numérotation des lignes :

  1. C ELQUOI SOURCE CB215821 19/08/20 21:16:58 10287
  2.  
  3. C=======================================================================
  4. C= Sous-programme ELQUOI =
  5. C= --------------------- =
  6. C= =
  7. C= Ce sous programme fournit les caracteristiques des elements finis. =
  8. C= =
  9. C= Entrees : =
  10. C= --------- =
  11. C= NELE : numero de l'element fini =
  12. C= NBPOIN : nombre de points d'integration contraintes =
  13. C= multicouches ( COQ6 - COQ8 ) =
  14. C= INTTYP : type d'integration utilisee par l'operateur =
  15. C= Si l'on desire un segment d'integration pour un champ aux : =
  16. C= 1 = NOEUDS ( appel a RENOEU ) =
  17. C= 2 = points de GAUSS pour centre de GRAVITE et CHAMP CONSTANT =
  18. C= 3 = points de GAUSS pour la RIGIDITE =
  19. C= 4 = points de GAUSS pour la MASSE =
  20. C= 5 = points de GAUSS pour les CONTRAINTES =
  21. C= 6 = on ne calcule pas de segment d'integration =
  22. C= 7 = FACE =
  23. C= 8 = CENTREP1 =
  24. C= 9 = MSOMMET =
  25. C -10= en mecanique on veut calculer les 5 premiers et on les met =
  26. C dans infmod(3 ...) =
  27. C= Dans un autre cas, infell(11) renverra un pointeur sur un segment =
  28. C= d'integration de type RIGIDITE. =
  29. C= IMODEL : pointeur sur une zone elementaire d'un modele =
  30. C= SEGMENT ACTIF en entree et en sortie =
  31. C= Sorties : =
  32. C= --------- =
  33. C= retourne IPTR pointeur sur un segment INFO ACTIF contenant : =
  34. C= La dimension de infell est fixee par NINFOS =
  35. C= infell( 1) numero de l'element fini =
  36. C= infell( 2) nombre de points d'integration en contraintes =
  37. C= multicouche =
  38. C= infell( 3) nombre de points d'integration pour la masse =
  39. C= infell( 4) nombre de points d'integration pour SIGMA =
  40. C= BSIGMA et KSIGMA =
  41. C= infell( 5) nombre de caracteristiques =
  42. C= infell( 6) nombre de points d'integration pour la RIGIDITE =
  43. C= infell( 7) longueur d'un tableau de travail pour l'element =
  44. C= infell( 8) nombre de fonctions de forme =
  45. C= infell( 9) nbre de d.d.l. dans la matrice de RIGIDITE =
  46. C= infell(10) taille de la matrice de Hooke =
  47. C= infell(11) pointeur sur le segment d'integration =
  48. C= infell(12) pointeur sur le 2nd segment d'integration des =
  49. C= elements homogeneises ou fluides (COQ6 ou COQ8) =
  50. C= infell(13) numero de la formulation de l'element fini =
  51. C= = 1 MASSIF =
  52. C= voir aussi = 3 COQUE =
  53. C= NUMMFR.eso = 5 COQUE EPAISSE ( COQ6 , COQ8 ) =
  54. C= = 7 POUTRE =
  55. C= = 9 COQUE AVEC CISAILLEMENT TRANSVERSE =
  56. C= = 11 LIQUIDE =
  57. C= = 13 TUYAU =
  58. C= = 15 LINESPRING =
  59. C= = 17 TUYAU FISSURE =
  60. C= = 19 RACCORD MASSIF =
  61. C= = 21 RACCORD COQUE =
  62. C= = 23 SURFACE LIBRE =
  63. C= = 25 MEMBRANE =
  64. C= = 26 point support MODAL =
  65. C= = 27 UNIAXIALE =
  66. C= = 28 point support STATIQUE =
  67. C= = 29 THERMIQUE =
  68. C= = 31 INCOMPRESSIBLES =
  69. C= = 33 POREUX =
  70. C= = 35 JOINT =
  71. C= = 37 HOMOGENEISE =
  72. C= = 39 TUYO =
  73. C= = 41 TUYAU ACOUSTIQUE PURE =
  74. C= = 43 RACCORD TUYAU FLUIDE =
  75. C= = 45 HYBRIDE =
  76. C= = 47 FIBRE =
  77. C= = 49 BARRE 3D EXCENTREE =
  78. C= = 51 LIASON A 2 NOEUDS =
  79. C= = 53 JOINT CISAILLEMENT (2D) =
  80. C= = 55 JOINT GENERALISE =
  81. C= = 57 POREUX TYPE Q =
  82. C= = 59 POREUX TYPE R =
  83. C= = 61 MACRO ELEMENT =
  84. C= = 63 XFEM meca rupt (2D/3D) =
  85. C= = 65 THERMOHYDRIQUE =
  86. C= = 69 MAGNETODYNAMIQUE =
  87. C= = 71 ELECTROSTATIQUE =
  88. C= = 73 DIFFUSION =
  89. C= = 75 JOINT UNIDIMENSIONNEL =
  90. C= = 77 ZONE COHESIVE =
  91. c = 79 TUYAU 1D pour advection
  92. C* = 78 COAXIAL
  93. c cccccc
  94. C= infell(14) numero de l'element geometrique associe (NUMGEO) =
  95. C= infell(15) nombre maximal de d.d.l. par noeud =
  96. C= infell(16) nombre de composantes de contraintes ou de deform. =
  97. C=======================================================================
  98.  
  99. SUBROUTINE ELQUOI(NELE,NBPOIN,INTTYP,IPTR,IPMODL)
  100.  
  101. IMPLICIT INTEGER(I-N)
  102. IMPLICIT REAL*8(A-H,O-Z)
  103.  
  104. -INC CCOPTIO
  105. -INC CCGEOME
  106. -INC CCHAMP
  107.  
  108. -INC SMMODEL
  109. -INC SMELEME
  110. -INC SMLMOTS
  111. -INC SMINTE
  112.  
  113. C Segment (type LISTENTI) contenant les informations sur un element
  114. SEGMENT INFO
  115. INTEGER INFELL(JG)
  116. ENDSEGMENT
  117.  
  118. C=======================================================================
  119. C= INTEGR : SAVE local des pointeurs sur les segments d'integration =
  120. C= CAS GENERAL =
  121. C= NELFI : Nombre de types d'element FINI (= LNOMTP de bdata.eso) =
  122. C= NBMODE : Nombre de MODEs de calcul (= nb valeurs possibles IFOMOD) =
  123. C= TYINTE : Type de l'INTEGRATION utilise par le calcul =
  124. C= (correspond a la valeur de INTTYP=infell(11) en entree) =
  125. C= NTYNTE : Nombre de TYPEs de segment D'INTEGRATION possibles =
  126. C= (voir la description de INTTYP ci-dessus) =
  127. C= NINFOS : Dimension (JG) du tableau infell (segment INFO) = 16 =
  128. C= On sauvegarde dans INTEGR(I,J) le pointeur sur un segment d'inte- =
  129. C= gration correspondant a l'element fini de type I et une integration =
  130. C= de type J. =
  131. C=======================================================================
  132. C Valeurs de IFOMOD : -1 0 1 2 3 4 5 6
  133. C PLAN AXIS FOUR TRID UNID UNID UNID FREQ
  134. C PLAN AXIS SPHE
  135. PARAMETER ( NELFI = 300 , NTYNTE = 9 , NBMODE = 8 )
  136. PARAMETER ( NINTEG = NELFI * NTYNTE * NBMODE , NINFOS = 16 )
  137.  
  138. INTEGER INTEGR(NELFI,NTYNTE,NBMODE),TYINTE
  139.  
  140. logical dcmate,dcmat2
  141.  
  142. SAVE INTEGR
  143. DATA INTEGR / NINTEG*0 /
  144.  
  145. IPTR = 0
  146.  
  147. C Tests de verification :
  148. IF (NELFI.LT.LNOMTP) THEN
  149. WRITE(IOIMP,*)
  150. WRITE(IOIMP,*) 'ELQUOI : valeur de NELFI incorrecte !'
  151. WRITE(IOIMP,*) ' NELFI = ',NELFI,LNOMTP
  152. CALL ERREUR(5)
  153. RETURN
  154. ENDIF
  155. C On se place sur la bonne "tranche" du tableau INTEGR :
  156. IBMODE = IFOMOD + 2
  157. IF (NBMODE.LT.IBMODE) THEN
  158. WRITE(IOIMP,*)
  159. WRITE(IOIMP,*) 'ELQUOI : valeur de NBMODE incorrecte !'
  160. WRITE(IOIMP,*) ' IBMODE = ',IBMODE,NBMODE,IFOMOD
  161. CALL ERREUR(5)
  162. RETURN
  163. ENDIF
  164. C Quel(s) type(s) de segment est(sont) a traiter ?
  165. TYINTE = INTTYP
  166. IF (INTTYP.EQ.-10) TYINTE = 0
  167.  
  168. IMODEL = IPMODL
  169.  
  170. dcmat2 = .false.
  171. C Recuperation des donnees de infell en entree
  172. MELE = NELE
  173. NBGS = NBPOIN
  174. IELE = NUMGEO(MELE)
  175. MFR = NUMMFR(MELE)
  176. C
  177. dcmate = .false.
  178. do im = 1, imodel.matmod(/2)
  179. if (imodel.matmod(im).eq.'IMPEDANCE') then
  180. dcmate = .true.
  181. meleme = imodel.imamod
  182. ityp1 = meleme.itypel
  183. mfr = nummfr(mele)
  184. if (ityp1.eq.1) mele = 45
  185. iele = numgeo(mele)
  186. endif
  187. enddo
  188. C Nombre de points d'integration dans l'epaisseur (stocke dans IMODEL)
  189. NPINT=0
  190. IF (imodel.INFMOD(/1).NE.0) NPINT = imodel.INFMOD(1)
  191. c IF(MELE.EQ.263) NPINT=0
  192. C Initialisation
  193. IPT2 = 0
  194. JG = NINFOS
  195. SEGINI,INFO
  196.  
  197. C=---------------------------------------------------------------------=
  198. C= REMPLISSAGE DU TABLEAU infell
  199. C=---------------------------------------------------------------------=
  200. C Remplissage de infell(1) : numero de l'element fini
  201. C --------------------------
  202. INFELL(1)=NELE
  203.  
  204. C Remplissage de infell(8) :
  205. C --------------------------
  206. C Sauf cas particuliers, le nombre de fonctions de forme de l'element
  207. C fini est egal au nombre de noeuds de l'element geometrique support.
  208. INFELL(8) = NBNNE(IELE)
  209. C Massifs : TRIx-QUAx 4 a 11, CUxx a PRxx 14 a 17, TExx-PYxx 23 a 26
  210. IF (MELE.GE.4.AND.MELE.LE.11) THEN
  211. INFELL(8) = 3
  212. ELSE IF (MELE.GE.14.AND.MELE.LE.17) THEN
  213. infell(8) = 3
  214. ELSE IF (MELE.GE.23.AND.MELE.LE.26) THEN
  215. infell(8) = 3
  216. C Polygones : POxx 111 a 122
  217. ELSE IF (MELE.GE.111.AND.MELE.LE.122) THEN
  218. IF (MELE.EQ.111) infell(8) = 3
  219. IF (MELE.EQ.112) infell(8) = 4
  220. IF (MELE.EQ.113) infell(8) = 5
  221. IF (MELE.EQ.114) infell(8) = 6
  222. IF (MELE.EQ.115) infell(8) = 7
  223. IF (MELE.EQ.116) infell(8) = 8
  224. IF (MELE.EQ.117) infell(8) = 9
  225. IF (MELE.EQ.118) infell(8) = 10
  226. IF (MELE.EQ.119) infell(8) = 11
  227. IF (MELE.EQ.120) infell(8) = 12
  228. IF (MELE.EQ.121) infell(8) = 13
  229. IF (MELE.EQ.122) infell(8) = 14
  230. C Elements QxRI 183 a 184
  231. ELSE IF (MELE.EQ.183.OR.MELE.EQ.184) THEN
  232. infell(8) = 3
  233. CTY Element XQ4R (xfem meca rupture en 2D)
  234. ELSEIF(MELE.eq.263) THEN
  235. infell(8)=40
  236. CTY Element XC8R (xfem meca rupture en 3D)
  237. ELSE IF(MELE.eq.264) THEN
  238. infell(8)=80
  239. ENDIF
  240.  
  241. C Elements Navier Stokes : cas P1CENTRE
  242. IF(TYINTE.EQ.8) THEN
  243. C cas MACRO 2D QUAF 2D --> TRI3
  244. IF(MELE.EQ.217.OR.MELE.EQ.224.OR.MELE.EQ.231
  245. & .OR.MELE.EQ.238.OR.MELE.EQ.245.OR.MELE.EQ.252)
  246. & infell(8) = 3
  247. C cas MACRO 2D QUAF 2D --> QUA4
  248. IF(MELE.EQ.218.OR.MELE.EQ.225.OR.MELE.EQ.232
  249. & .OR.MELE.EQ.239.OR.MELE.EQ.246.OR.MELE.EQ.253)
  250. & infell(8) = 3
  251. C cas MACRO 3D QUAF 3D
  252. C IF((MELE.GE.226.AND.MELE.LE.229).OR.
  253. C & (MELE.GE.247.AND.MELE.LE.250)) infell(8) = 4
  254. ENDIF
  255. C Elements Navier Stokes : cas MSOMMET
  256. IF(TYINTE.EQ.9) THEN
  257. C cas MACRO QUAF 1D --> SEG2
  258. IF(MELE.EQ.216.OR.MELE.EQ.223.OR.MELE.EQ.230
  259. & .OR.MELE.EQ.237.OR.MELE.EQ.244.OR.MELE.EQ.251) IELE = 2
  260. C cas MACRO QUAF 2D --> TRI3
  261. IF(MELE.EQ.217.OR.MELE.EQ.224.OR.MELE.EQ.231
  262. & .OR.MELE.EQ.238.OR.MELE.EQ.245.OR.MELE.EQ.252) IELE = 4
  263. C cas MACRO QUAF 2D --> QUA4
  264. IF(MELE.EQ.218.OR.MELE.EQ.225.OR.MELE.EQ.232
  265. & .OR.MELE.EQ.239.OR.MELE.EQ.246.OR.MELE.EQ.253) IELE = 8
  266. C cas MACRO QUAF 3D --> CUB8
  267. IF(MELE.EQ.233.OR.MELE.EQ.254) IELE = 14
  268. C cas MACRO QUAF 3D --> PRI6
  269. IF(MELE.EQ.234.OR.MELE.EQ.255) IELE = 16
  270. C cas MACRO QUAF 3D --> TET4
  271. IF(MELE.EQ.235.OR.MELE.EQ.256) IELE = 23
  272. C cas MACRO QUAF 3D --> PYR5
  273. IF(MELE.EQ.236.OR.MELE.EQ.257) IELE = 25
  274. infell(8) = NBNNE(IELE)
  275. infell(14)= IELE
  276. ENDIF
  277.  
  278. C Pour les elements xxxP (79 a 83), TUYO (96), JOPx (108 a 110), xxxQ
  279. C (173 a 177), xxxR (178 a 182), JORx (185 a 187) et JOQx(188 a 190),
  280. C infell(8) est affecte lors du remplissage complet de infell.
  281.  
  282. C Remplissage de infell(13) et infell(14) :
  283. C -----------------------------------------
  284. C infell(13) : numero de la FORMULATION - MFR=NUMMFR(NELE)
  285. infell(13) = MFR
  286. do im = 1,matmod(/2)
  287. if (MFR.eq.27.and.matmod(im).eq.'MODAL') infell(13) = 26
  288. if (MFR.eq.27.and.matmod(im).eq.'STATIQUE') infell(13) = 28
  289. enddo
  290. IF (FORMOD(1).EQ.'CHARGEMENT') THEN
  291. IF (MFR.EQ.1) INFELL(13)=72
  292. IF (MFR.EQ.3.OR.MFR.EQ.5.OR.MFR.EQ.7.OR.MFR.EQ.9.OR.
  293. & MFR.EQ.13) INFELL(13) = 74
  294. ENDIF
  295. MFR = infell(13)
  296. mfr3 = nummfr(nele)
  297. C mfr2 non nul en cas de formulation massive non mecanique
  298. mfr2 = 0
  299. if (MFR.EQ.26.OR.MFR.EQ.28) mfr3 = mfr
  300. C Cas particuliers : elements HYxx (99 a 103)
  301. C MFR = NUMMFR(mele) = 45 mais infell(13)=1 (MASSIF)
  302. IF (MELE.GE.99.AND.MELE.LE.103) infell(13)=1
  303.  
  304. C ELECTROSTATIQUE : MFR=71=infell(13) et NUMMFR(mele)=1 (MASSIF)
  305. IF (FORMOD(1).EQ.'ELECTROSTATIQUE') mfr2=71
  306.  
  307. C DIFFUSION : MFR=73=infell(13) et NUMMFR(mele)=1 (MASSIF)
  308. IF (FORMOD(1).EQ.'DIFFUSION') THEN
  309. IF (MFR .EQ. 1) THEN
  310. infell(13)=73
  311. ELSE
  312. infell(13)=MFR
  313. ENDIF
  314. ENDIF
  315.  
  316. IF (mfr2.NE.0) infell(13) = mfr2
  317. C
  318. C infell(14) : numero de l'element GEOMETRIQUE - IELE=NUMGEO(NELE)
  319. infell(14) = IELE
  320.  
  321. C Remplissage de infell(10), infell(15) et infell(16) :
  322. C -----------------------------------------------------
  323. C infell(10) : dimension de la matrice de Hooke
  324. C infell(15) : nombre maximal de ddl par noeud
  325. C infell(16) : nombre de composantes de contraintes et deformations
  326.  
  327. C MASSIFs (2D) : TRIx-QUAx 4 a 11 + POLYgones : POLY 111 a 122 +
  328. C Elements QxRI 183 a 184
  329. IF ((MELE.GE.4.AND.MELE.LE.11).OR.(MELE.GE.111.AND.MELE.LE.122)
  330. . .OR.MELE.EQ.183.OR.MELE.EQ.184) THEN
  331. IF (IFOUR.EQ.1) THEN
  332. infell(10)=6
  333. infell(15)=3
  334. infell(16)=6
  335. ELSE
  336. infell(10)=4
  337. infell(15)=2
  338. infell(16)=4
  339. ENDIF
  340. C MASSIFs (3D) : CUxx-PRxx 14 a 17, TExx-PYxx 23 a 26
  341. ELSE IF ((MELE.GE.14.AND.MELE.LE.17).OR.
  342. . (MELE.GE.23.AND.MELE.LE.26)) THEN
  343. infell(10)=6
  344. infell(15)=3
  345. infell(16)=6
  346. C Elements LIQUIDE : Lxxx 35 a 40
  347. ELSE IF (MELE.GE.35.AND.MELE.LE.40) THEN
  348. infell(10)=5
  349. infell(15)=2
  350. infell(16)=0
  351. C Elements de RACCORD FLUIDE-STRUCTURE :
  352. C RAC2 12, LIAx 18 a 19, RACO 47, LICO 55, LIC4 94
  353. ELSE IF (MELE.EQ.12.OR.MELE.EQ.18.OR.MELE.EQ.19.OR.MELE.EQ.47.OR.
  354. . MELE.EQ.55.OR.MELE.EQ.94) THEN
  355. infell(10)=6
  356. C Elements de SURFACE LIBRE : LSU2 48, LSU3 53, LSU4 54
  357. ELSE IF (MELE.EQ.48.OR.MELE.EQ.53.OR.MELE.EQ.54) THEN
  358. infell(10)=6
  359. infell(15)=3
  360. C Elements THERMIQUEs : CVxx 57 a 62, THxx 63 a 68
  361. ELSE IF (MELE.GE.57.AND.MELE.LE.68) THEN
  362. IF (MELE.LE.60) THEN
  363. infell(10)=2
  364. infell(15)=1
  365. infell(16)=2
  366. ELSE IF (MELE.GE.61) THEN
  367. infell(10)=3
  368. infell(15)=1
  369. infell(16)=3
  370. ENDIF
  371. C Elements MASSIFs IMCOMPRESSIBLES : ICxx 69 a 78
  372. ELSE IF (MELE.GE.69.AND.MELE.LE.78) THEN
  373. IF (MELE.LE.72) THEN
  374. infell(10)=4
  375. C** IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(10)=4
  376. infell(15)=2
  377. IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(15)=3
  378. infell(16)=4
  379. C** IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(16)=4
  380. ELSE IF (MELE.GE.73) THEN
  381. infell(10)=6
  382. infell(15)=3
  383. infell(16)=6
  384. ENDIF
  385. C Elements MASSIFs POREUX : xxxP 79 a 83
  386. ELSE IF (MELE.GE.79.AND.MELE.LE.83) THEN
  387. IF (MELE.LE.80) THEN
  388. IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) THEN
  389. infell(10)=6
  390. infell(15)=4
  391. infell(16)=7
  392. ELSE
  393. infell(10)=4
  394. infell(15)=3
  395. infell(16)=5
  396. ENDIF
  397. ELSE IF (MELE.GE.81) THEN
  398. infell(10)=6
  399. infell(15)=4
  400. infell(16)=7
  401. ENDIF
  402. C Elements MASSIFs POREUX de type Q : xxxQ 173 a 177
  403. ELSE IF (MELE.GE.173.AND.MELE.LE.177) THEN
  404. IF (MELE.LE.174) THEN
  405. IF (IFOUR.EQ.1) THEN
  406. infell(10)=8
  407. infell(15)=5
  408. infell(16)=8
  409. ELSE
  410. infell(10)=6
  411. infell(15)=4
  412. infell(16)=6
  413. ENDIF
  414. ELSE IF (MELE.GE.175) THEN
  415. infell(10)=8
  416. infell(15)=5
  417. infell(16)=8
  418. ENDIF
  419. C Elements MASSIFs POREUX de type R : xxxR 178 a 182
  420. ELSE IF(MELE.GE.178.AND.MELE.LE.182) THEN
  421. IF (MELE.LE.179) THEN
  422. IF (IFOUR.EQ.1) THEN
  423. infell(10)=9
  424. infell(15)=6
  425. infell(16)=9
  426. ELSE
  427. infell(10)=7
  428. infell(15)=5
  429. infell(16)=7
  430. ENDIF
  431. ELSE IF (MELE.GE.180) THEN
  432. infell(10)=9
  433. infell(15)=6
  434. infell(16)=9
  435. ENDIF
  436. C Elements JOINTS : JOxx 85 a 90
  437. ELSE IF (MELE.GE.85.AND.MELE.LE.90) THEN
  438. IF (MELE.LE.86) THEN
  439. infell(10)=2
  440. infell(16)=2
  441. ELSE IF (MELE.GE.87) THEN
  442. infell(10)=3
  443. infell(16)=3
  444. ENDIF
  445. C
  446. C MACRO ELEMENT CISAILLEMENT FLEXION
  447. C
  448. ELSE IF(MELE.EQ.191) THEN
  449. infell(5)=3
  450. infell(7)=300
  451. GO TO 258
  452.  
  453. C Elements JOINTs CISAILLEMENT et GENERALISES : JCxx-JGxx 168 a 172
  454. ELSE IF (MELE.GE.168.AND.MELE.LE.172) THEN
  455. infell(7)=300
  456. infell(9)=infell(8)*IDIM
  457. infell(15)=IDIM
  458. IF (MELE.EQ.168.OR.MELE.EQ.169) THEN
  459. infell(10)=2
  460. infell(16)=2
  461. ELSE IF (MELE.EQ.170) THEN
  462. infell(10)=4
  463. infell(16)=4
  464. ELSE
  465. infell(10)=6
  466. infell(16)=6
  467. ENDIF
  468. C Elements HOMOGENEISES : TRIH 92, QUAH 126, CUBH 127, TRH6 157
  469. ELSE IF (MELE.EQ.92.OR.MELE.EQ.126.OR.MELE.EQ.127.OR.
  470. . MELE.EQ.157) THEN
  471. infell(10)=10
  472. infell(16)=0
  473. C Elements JOINTS POREUX : JOPx 108 a 110
  474. ELSE IF (MELE.GE.108.AND.MELE.LE.110) THEN
  475. infell(7)=100
  476. C Elements JOINTS POREUX de type Q : JOQx 185 a 187
  477. ELSE IF (MELE.GE.185.AND.MELE.LE.187) THEN
  478. infell(7)=150
  479. C Elements JOINTS POREUX de type R : JORx 188 a 190
  480. ELSE IF (MELE.GE.188.AND.MELE.LE.190) THEN
  481. infell(7)=200
  482. C Elements MASSIFS THERMIQUE/MECANIQUE en 1D : T1Dx-M1Dx 191 a 194
  483. ELSE IF (MELE.GE.191.AND.MELE.LE.194) THEN
  484. infell(10)=3
  485. infell(15)=1
  486. infell(16)=3
  487. ENDIF
  488. C Cas particuliers des autres formulations massives (mfr2 > 0)
  489. IF (mfr2.NE.0) THEN
  490. infell(10)=IDIM
  491. infell(15)=1
  492. infell(16)=IDIM
  493. ENDIF
  494. C Pour les autres elements, les valeurs de infell(10,15,16) sont
  495. C eventuellement affectees (0 par defaut) lors du branchement (GOTO)
  496. C ci-dessous qui effectue le remplissage complet du tableau infell.
  497.  
  498. C Remplissage complet du tableau infell suivant l'element fini
  499. C --------------------------------------------------------------
  500. IF(MELE.GE.1.AND.MELE.LE.194) THEN
  501. C CABL SEG2 SEG3 TRI3 TRI4 TRI6 TRI7 QUA4 QUA5 QUA8 QUA9,RAC2
  502. GOTO ( 300, 2, 3, 4, 300, 6, 300, 8, 300, 10, 300, 12
  503. C RAC3 CUB8 CU20 PRI6 PR15 LIA3 LIA4 LIA6 LIA8 MULT TET4 TE10
  504. . , 300, 14, 15, 16, 17, 18, 19, 300, 300, 300, 23, 24
  505. C PYR5 PY13 COQ3 DKT POUT LISP FAC3 FAC4 FAC6 FAC8 LTR3 LQU4
  506. . , 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36
  507. C LCU8 LPR6 LTE4 LPY5 COQ8 TUYA TUFI COQ2 POI1 BARR RACO LSU2
  508. . , 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48
  509. C COQ4 LISM COF3 RES2 LSU3 LSU4 LICO COQ6 CVS2 CVS3 CVT3 CVT6
  510. . , 49, 50, 51, 300, 53, 54, 55, 56, 57, 58, 59, 60
  511. C CVQ4 CVQ8 THP5 TH13 THP6 TH15 THC8 TH20 ICT3 ICQ4 ICT6 ICQ8
  512. . , 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71, 72
  513. C ICC8 ICT4 ICP6 IC20 IC10 IC15 TRIP QUAP CUBP TETP PRIP TIMO
  514. . , 73, 74, 75, 76, 77, 78, 79, 80, 81, 82, 83, 84
  515. C JOI2 JOI3 JOT3 JOI4 JOI6 JOI8 LISC TRIH DST LIC4 CERC TUYO
  516. . , 85, 86, 87, 88, 89, 90, 91, 92, 93, 94, 95, 96
  517. C LSE2 LITU HYT3 HYQ4 HYT4 HYP6 HYC8 TRIS QUAS POIS FOR3 JOP3
  518. . , 97, 98, 99, 100, 101, 102, 103, 104, 105, 106, 300, 108
  519. C JOP6 JOP8 POL3 POL4 POL5 POL6 POL7 POL8 POL9 PO10 PO11 PO12
  520. . , 109, 110, 111, 111, 111, 111, 111, 111, 111, 111, 111, 111
  521. C PO13 PO14 BAR3 BAEX LIA2 QUAH CUBH ROT3 SEF2 TRF3 QUF4 CUF8
  522. . , 111, 111, 123, 124, 125, 126, 127, 128, 300, 300, 300, 300
  523. C PRF6 TEF4 PYF5 MSE3 MTR6 MQU9 MC27 MP18 MT10 MP14 SEF3 TRF7
  524. . , 300, 300, 300, 300, 300, 300, 300, 300, 300, 300, 300, 300
  525. C QUF9 CF27 PF21 TF15 PF19 SEG6 TR21 QU36 C216 P126 TE56 PY91
  526. . , 300, 300, 300, 300, 300, 300, 300, 300, 300, 300, 300, 300
  527. C TRH6 BSE2 BTR4 BQU5 BCU9 BPR7 BTE5 BPY6 FRO4 SEGS POJS JCT3
  528. . , 157, 300, 300, 300, 300, 300, 300, 300, 300, 166, 167, 168
  529. C JCI4 JGI2 JGT3 JGI4 TRIQ QUAQ CUBQ TETQ PRIQ TRIR QUAR CUBR
  530. . , 169, 170, 171, 172, 173, 174, 175, 176, 177, 178, 179, 180
  531. C TETR PRIR Q4RI Q8RI JOQ3 JOQ6 JOQ8 JOR3 JOR6 JOR8 T1D2 T1D3
  532. . , 181, 182, 183, 184, 185, 186, 187, 188, 189, 190, 191, 192
  533. C M1D2 M1D3
  534. . , 191, 192 ), MELE
  535. ELSEIF(MELE.GE.195.AND.MELE.LE.272) THEN
  536. C LC03 LC07 LC09 LC27 LC21 LC15 LC19
  537. GOTO ( 195, 196, 197, 300, 300, 300, 300
  538. C LS03 LS07 LS09 LS27 LS21 LS15 LS19
  539. . , 195, 196, 197, 300, 300, 300, 300
  540. C BS03 BS07 BS09 BS27 BS21 BS15 BS19
  541. . , 195, 196, 197, 300, 300, 300, 300
  542. C MC03 MC07 MC09 MC27 MC21 MC15 MC19
  543. . , 216, 217, 218, 300, 300, 300, 300
  544. C M103 M107 M109 M127 M121 M115 M119
  545. . , 223, 224, 225, 300, 300, 300, 300
  546. C MS03 MS07 MS09 MS27 MS21 MS15 MS19
  547. . , 223, 224, 225, 300, 300, 300, 300
  548. C QC03 QC07 QC09 QC27 QC21 QC15 QC19
  549. . , 237, 238, 239, 300, 300, 300, 300
  550. C Q103 Q107 Q109 Q127 Q121 Q115 Q119
  551. . , 237, 238, 239, 300, 300, 300, 300
  552. C QS03 QS07 QS09 QS27 QS21 QS15 QS19 CIFL SURE
  553. . , 237, 238, 239, 300, 300, 300, 300, 300, 300
  554. C SHB8 caf2 caf3 xq4r xc8r joi1 ZCO2 ZCO3 ZCO4 TUY2 TUY3 COS2
  555. . , 260, 300, 300, 263, 264, 265, 266, 267, 268, 269, 270, 271
  556. C COA2 ICY5 IC13
  557. . , 272, 25, 26)
  558. c cccccc
  559. . , (MELE-194)
  560. ENDIF
  561. C
  562. C Remplissage de infell pour chaque element
  563. C -------------------------------------------
  564. C Element SEG2
  565. 2 infell(2)=2
  566. infell(4)=2
  567. infell(6)=2
  568. C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=3
  569. infell(7)=0
  570. infell(9)=4
  571. IF (IFOUR.EQ.1) infell(9)=6
  572. infell(10)=1
  573. infell(15)=2
  574. IF (IFOUR.EQ.1) infell(15)=3
  575. infell(16)=1
  576. GOTO 300
  577. C Element SEG3
  578. 3 infell(2)=3
  579. IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(2)=4
  580. infell(4)=3
  581. IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(4)=4
  582. infell(6)=3
  583. IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=4
  584. infell(7)=0
  585. infell(9)=6
  586. IF (IFOUR.EQ.1) infell(9)=9
  587. infell(10)=1
  588. infell(15)=2
  589. IF (IFOUR.EQ.1) infell(15)=3
  590. infell(16)=1
  591. GOTO 300
  592. C Element TRI3
  593. 4 infell(2)=1
  594. infell(3)=4
  595. IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(3)=7
  596. infell(4)=1
  597. infell(6)=1
  598. C IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=4
  599. infell(9)=6
  600. IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(9)=9
  601. GOTO 300
  602. C Element TRI6
  603. 6 IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN
  604. infell(2)=7
  605. infell(4)=7
  606. infell(6)=7
  607. ELSE
  608. infell(2)=4
  609. infell(4)=4
  610. infell(6)=4
  611. ENDIF
  612. infell(3)=7
  613. infell(9)=12
  614. IF (IFOUR.EQ.-3) infell(9)=15
  615. IF (IFOUR.EQ.1) infell(9)=18
  616. GOTO 300
  617. C Element QUA4
  618. 8 infell(2)=4
  619. infell(3)=4
  620. IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(3)=9
  621. infell(4)=4
  622. C** infell(4)=5 CORRECTION MILL LE 8/8/90
  623. infell(6)=4
  624. C** IF(IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=9
  625. infell(9)=8
  626. IF (IFOUR.EQ.-3) infell(9)=11
  627. IF (IFOUR.EQ.1) infell(9)=12
  628. GOTO 300
  629. C Element QUA8
  630. 10 infell(2)=9
  631. infell(3)=9
  632. IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(3)=16
  633. infell(4)=9
  634. infell(6)=9
  635. C IF(IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=16
  636. infell(9)=16
  637. IF (IFOUR.EQ.-3) infell(9)=19
  638. IF (IFOUR.EQ.1) infell(9)=24
  639. GOTO 300
  640. C Element RAC2 (raccord a 4 noeuds en 2 dimensions)
  641. 12 infell(2)=2
  642. infell(4)=2
  643. infell(5)=2
  644. infell(6)=2
  645. infell(9)=8
  646. IF (IFOUR.EQ.1) infell(9)=10
  647. infell(15)=3
  648. GOTO 300
  649. C Element CUB8
  650. 14 infell(2)=8
  651. infell(3)=8
  652. infell(4)=8
  653. infell(6)=8
  654. infell(9)=24
  655. GOTO 300
  656. C Element CU20
  657. 15 infell(2)=27
  658. infell(3)=27
  659. infell(4)=27
  660. infell(6)=27
  661. infell(9)=60
  662. GOTO 300
  663. C Element PRI6
  664. 16 infell(2)=6
  665. cbp infell(3)=6
  666. infell(3)=8
  667. infell(4)=6
  668. infell(6)=6
  669. infell(9)=18
  670. GOTO 300
  671. C Element PR15
  672. 17 infell(2)=12
  673. cbp infell(3)=12
  674. infell(3)=21
  675. infell(4)=12
  676. infell(6)=12
  677. infell(9)=45
  678. GOTO 300
  679. C Element LIA3 (raccord a 6 noeudes en 3 dimensions)
  680. 18 infell(2)=4
  681. infell(4)=4
  682. infell(5)=3
  683. infell(6)=4
  684. infell(9)=15
  685. infell(13)=19
  686. infell(15)=3
  687. GOTO 300
  688. C Element LIA4 (raccord a 8 noeuds en 3 dimensions)
  689. 19 infell(2)=4
  690. infell(4)=4
  691. infell(5)=3
  692. infell(6)=4
  693. infell(9)=20
  694. infell(15)=3
  695. GOTO 300
  696. C Element TET4
  697. 23 infell(2)=1
  698. infell(3)=4
  699. infell(4)=1
  700. infell(6)=1
  701. infell(9)=12
  702. GOTO 300
  703. C Element TE10
  704. 24 infell(2)=5
  705. cbp infell(3)=5
  706. infell(3)=15
  707. infell(4)=5
  708. infell(6)=5
  709. infell(9)=30
  710. GOTO 300
  711. C Element PYR5
  712. 25 infell(2)=5
  713. infell(3)=5
  714. infell(4)=5
  715. infell(6)=5
  716. infell(9)=15
  717. GOTO 300
  718. C Element PY13
  719. cbp 26 infell(2)=27
  720. cbp infell(3)=27
  721. cbp infell(4)=27
  722. cbp infell(6)=27
  723. 26 infell(2)=10
  724. infell(3)=17
  725. infell(4)=10
  726. infell(6)=10
  727. infell(9)=39
  728. GOTO 300
  729. C Element COQ3
  730. 27 infell(2)=3
  731. infell(4)=1
  732. infell(5)=3
  733. infell(6)=1
  734. infell(7)=600
  735. infell(9)=18
  736. IF (FORMOD(1) .EQ. 'DIFFUSION') infell(9)=9
  737. infell(10)=6
  738. infell(15)=6
  739. infell(16)=6
  740. GOTO 300
  741. C Element DKT
  742. 28 infell(2)=3
  743. infell(3)=3
  744. infell(5)=3
  745. infell(7)=600
  746. infell(9)=18
  747. infell(15)=6
  748. IF (NPINT.EQ.0) THEN
  749. C Element DKT COQUE GLOBALE
  750. infell(4)=3
  751. infell(6)=3
  752. infell(10)=6
  753. infell(16)=6
  754. ELSE
  755. C Element DKT COQUE INTEGREE dans l'epaisseur
  756. infell(4)=3*NPINT
  757. infell(6)=3*NPINT
  758. infell(10)=6
  759. infell(16)=4
  760. ENDIF
  761. GOTO 300
  762. C Element POUTRE
  763. 29 infell(2)=2
  764. infell(4)=2
  765. infell(6)=2
  766. infell(7)=200
  767. IF (IFOUR.EQ.-3.OR.IFOUR.EQ.-2.OR.IFOUR.EQ.-1) THEN
  768. C Element POUTRE 2D plan
  769. infell(5)=6
  770. infell(9)=6
  771. infell(10)=3
  772. infell(15)=3
  773. infell(16)=3
  774. ELSE
  775. C Element POUTRE 3D
  776. infell(5)=12
  777. infell(9)=12
  778. infell(10)=6
  779. infell(15)=6
  780. infell(16)=6
  781. ENDIF
  782. GOTO 300
  783. C Element LINE SPRING LISP
  784. 30 infell(2)=3
  785. infell(13)=15
  786. infell(6)=3
  787. infell(4)=3
  788. infell(16)=6
  789. infell(9)=24
  790. infell(7)=100
  791. infell(10)=6
  792. infell(15)=6
  793. infell(5)=5
  794. GOTO 300
  795. C Element FAC3 pour les PRESSIONS uniquement
  796. 31 infell(2)=1
  797. infell(4)=1
  798. infell(6)=1
  799. infell(9)=9
  800. infell(15)=3
  801. GOTO 300
  802. C Element FAC4 pour les PRESSIONS uniquement
  803. 32 infell(2)=4
  804. infell(4)=4
  805. infell(6)=4
  806. infell(9)=12
  807. infell(15)=3
  808. GOTO 300
  809. C Element FAC6 pour les PRESSIONS uniquement
  810. 33 infell(2)=4
  811. infell(4)=4
  812. infell(6)=4
  813. infell(9)=18
  814. infell(15)=3
  815. GOTO 300
  816. C Element FAC8 pour les PRESSIONS uniquement
  817. 34 infell(2)=9
  818. infell(4)=9
  819. infell(6)=9
  820. infell(9)=24
  821. infell(15)=3
  822. GOTO 300
  823. C Element LTR3 : element LIQUIDE triangle a 3 noeuds
  824. 35 infell(2)=4
  825. infell(3)=4
  826. IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(3)=7
  827. infell(4)=4
  828. infell(6)=4
  829. infell(9)=6
  830. GOTO 300
  831. C Element LQU4 : element LIQUIDE quadrangle a 4 noeuds
  832. 36 infell(2)=4
  833. infell(3)=4
  834. IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(3)=9
  835. infell(4)=4
  836. infell(6)=4
  837. infell(9)=8
  838. GOTO 300
  839. C Element LCU8 : Element LIQUIDE CUBE a 8 noeuds
  840. 37 infell(2)=8
  841. infell(3)=8
  842. infell(4)=8
  843. infell(6)=8
  844. infell(9)=16
  845. GOTO 300
  846. C Element LPR6 : Element LIQUIDE PRISME a 6 noeuds
  847. 38 infell(2)=6
  848. infell(3)=6
  849. infell(4)=6
  850. infell(6)=6
  851. infell(9)=12
  852. GOTO 300
  853. C Element LTE4 : Element LIQUIDE TETRAEDRE a 4 noeuds
  854. 39 infell(2)=4
  855. infell(3)=4
  856. infell(4)=4
  857. infell(6)=4
  858. infell(9)=8
  859. GOTO 300
  860. C Element LPY5 : Element LIQUIDE PYRAMIDE a 5 noeuds
  861. 40 infell(2)=5
  862. infell(3)=5
  863. infell(4)=5
  864. infell(6)=5
  865. infell(9)=10
  866. GOTO 300
  867. C Element COQ8 : Element COQUE EPAISSE a 8 noeuds
  868. 41 infell(2)=NBGS
  869. IF (NBGS.EQ.0) infell(2)=8
  870. infell(3)=8
  871. infell(4)=NBGS
  872. IF (NBGS.EQ.0) infell(4)=8
  873. infell(5)=2
  874. infell(6)=8
  875. infell(7)=1700
  876. infell(9)=48
  877. IF (FORMOD(1) .EQ. 'DIFFUSION') infell(9)=24
  878. infell(10)=5
  879. infell(15)=6
  880. infell(16)=5
  881. GOTO 300
  882. C Element TUYA : Element TUYAU
  883. 42 infell(2)=2
  884. infell(4)=2
  885. infell(5)=7
  886. infell(6)=2
  887. infell(7)=200
  888. infell(9)=12
  889. infell(10)=6
  890. infell(15)=6
  891. infell(16)=6
  892. GOTO 300
  893. C Element TUFI : Element TUYAU FISSURE
  894. 43 infell(2)=1
  895. infell(3)=0
  896. infell(4)=1
  897. infell(5)=9
  898. infell(6)=1
  899. infell(7)=200
  900. infell(9)=12
  901. infell(10)=6
  902. infell(15)=6
  903. infell(16)=8
  904. GOTO 300
  905. C Element COQ2
  906. 44 infell(2)=3
  907. infell(3)=3
  908. infell(4)=3
  909. infell(5)=3
  910. infell(6)=3
  911. infell(7)=368
  912. infell(9)=8
  913. IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.0) infell(9)=6
  914. IF (IFOUR.EQ.-3) infell(9)=9
  915. infell(10)=6
  916. IF (IFOUR.LE.0) infell(10)=4
  917. infell(15)=4
  918. IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.0) infell(15)=3
  919. infell(16)=6
  920. IF (IFOUR.LE.0) infell(16)=4
  921. GOTO 300
  922. C Element POI1 (MODEs en DEFOrmations PLANes GENEralisees)
  923. C Valeurs a adapter pour certains MODE UNID (1D)
  924. 45 infell(2)=1
  925. infell(3)=1
  926. infell(4)=1
  927. infell(5)=1
  928. infell(6)=1
  929. infell(7)=50
  930. infell(15)=2
  931. IF (IFOUR.EQ.11) THEN
  932. infell(9)=3
  933. infell(10)=2
  934. infell(16)=2
  935. ELSE
  936. infell(9)=5
  937. IF (IFOUR.GE.3.AND.IFOUR.LE.15) infell(9)=2
  938. infell(10)=1
  939. infell(16)=1
  940. ENDIF
  941. if (mfr.eq.26.or.mfr.eq.28) infell(9) = 1
  942. GOTO 300
  943. C Element BARRE
  944. 46 infell(2)=2
  945. infell(4)=2
  946. infell(5)=1
  947. infell(6)=2
  948. infell(7)=200
  949. infell(9)=6
  950. IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.0) infell(9)=4
  951. IF (IFOUR.EQ.-3) infell(9)=7
  952. IF (FORMOD(1) .EQ. 'DIFFUSION') infell(9)=2
  953. infell(10)=1
  954. infell(15)=3
  955. IF (FORMOD(1) .EQ. 'DIFFUSION') infell(15)=1
  956. IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.0) infell(15)=2
  957. infell(16)=1
  958. GOTO 300
  959. C Element RACO : Element de RACCORD LIQUIDE-COQUE a 4 noeuds en 2D
  960. 47 infell(2)=1
  961. infell(4)=1
  962. infell(5)=2
  963. infell(6)=1
  964. infell(9)=10
  965. IF (IFOUR.EQ.1) infell(9)=12
  966. infell(15)=4
  967. GOTO 300
  968. C Element LSU2 : Element de SURFACE LIBRE a 2 noeuds en 2D
  969. 48 infell(2)=2
  970. infell(4)=2
  971. infell(6)=2
  972. infell(9)=6
  973. GOTO 300
  974. C Element COQ4
  975. 49 infell(2)=5
  976. infell(3)=5
  977. infell(4)=5
  978. infell(5)=3
  979. infell(6)=5
  980. infell(7)=100
  981. infell(9)=24
  982. IF (FORMOD(1) .EQ. 'DIFFUSION') infell(9)=12
  983. infell(10)=8
  984. infell(15)=6
  985. infell(16)=8
  986. GOTO 300
  987. C Element LINE SPRING LISM
  988. 50 infell(2)=1
  989. infell(4)=1
  990. infell(5)=5
  991. infell(6)=3
  992. infell(7)=100
  993. infell(9)=24
  994. infell(10)=6
  995. infell(15)=6
  996. infell(16)=6
  997. GOTO 300
  998. C Element COF3
  999. 51 infell(2)=3
  1000. infell(3)=3
  1001. infell(4)=3
  1002. infell(5)=3
  1003. infell(6)=3
  1004. infell(9)=12
  1005. infell(7)=368
  1006. infell(10)=6
  1007. infell(15)=4
  1008. infell(16)=6
  1009. GOTO 300
  1010. C Element LSU3 : Element de SURFACE LIBRE a 3 noeuds en 3D
  1011. 53 infell(2)=4
  1012. infell(4)=4
  1013. infell(6)=4
  1014. infell(9)=9
  1015. GOTO 300
  1016. C Element LSU4 : Element de SURFACE LIBRE a 4 noeuds en 3D
  1017. 54 infell(2)=4
  1018. infell(4)=4
  1019. infell(6)=4
  1020. infell(9)=12
  1021. GOTO 300
  1022. C Element LICO : Element de RACCORD LIQUIDE-COQUE a 6 noeuds en 3D
  1023. 55 infell(2)=3
  1024. infell(4)=3
  1025. infell(5)=3
  1026. infell(6)=3
  1027. infell(9)=24
  1028. infell(15)=6
  1029. GOTO 300
  1030. C Element COQ6 : Element COQUE EPAISSE a 6 noeuds
  1031. 56 infell(2)=NBGS
  1032. IF(NBGS.EQ.0) infell(2)=6
  1033. infell(3)=6
  1034. infell(4)=NBGS
  1035. IF (NBGS.EQ.0) infell(4)=6
  1036. infell(5)=2
  1037. infell(6)=6
  1038. infell(9)=36
  1039. IF (FORMOD(1) .EQ. 'DIFFUSION') infell(9)=18
  1040. infell(7)=1000
  1041. infell(10)=5
  1042. infell(15)=6
  1043. infell(16)=5
  1044. GOTO 300
  1045. C Element THT3 (thermique)
  1046. 57 infell(2)=1
  1047. infell(6)=1
  1048. C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=4
  1049. infell(9)=3
  1050. IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(9)=9
  1051. GOTO 300
  1052. C Element THT6
  1053. 58 infell(2)=4
  1054. infell(4)=4
  1055. infell(6)=4
  1056. IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=7
  1057. infell(9)=6
  1058. IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(9)=18
  1059. GOTO 300
  1060. C Element THQ4
  1061. 59 infell(2)=4
  1062. infell(6)=4
  1063. C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=9
  1064. infell(9)=4
  1065. IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(9)=12
  1066. GOTO 300
  1067. C Element THQ8
  1068. 60 infell(2)=9
  1069. infell(6)=9
  1070. C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=16
  1071. infell(9)=8
  1072. IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(9)=24
  1073. GOTO 300
  1074. C Element THT4
  1075. 61 infell(2)=1
  1076. infell(6)=1
  1077. infell(9)=4
  1078. GOTO 300
  1079. C Element TH10
  1080. 62 infell(2)=5
  1081. infell(6)=5
  1082. infell(9)=10
  1083. GOTO 300
  1084. C Element THP5
  1085. 63 infell(2)=5
  1086. infell(6)=5
  1087. infell(9)=5
  1088. GOTO 300
  1089. C Element TH13
  1090. 64 infell(2)=27
  1091. infell(6)=27
  1092. infell(9)=13
  1093. GOTO 300
  1094. C Element THP6
  1095. 65 infell(2)=2
  1096. infell(6)=2
  1097. infell(9)=6
  1098. GOTO 300
  1099. C Element TH15
  1100. 66 infell(2)=12
  1101. infell(6)=12
  1102. infell(9)=15
  1103. GOTO 300
  1104. C Element THC8
  1105. 67 infell(2)=8
  1106. infell(6)=8
  1107. infell(9)=8
  1108. GOTO 300
  1109. C Element TH20
  1110. 68 infell(2)=27
  1111. infell(6)=27
  1112. infell(9)=20
  1113. GOTO 300
  1114. C Element ICT3
  1115. 69 infell(2)=1
  1116. infell(3)=4
  1117. IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(3)=7
  1118. infell(4)=1
  1119. infell(6)=1
  1120. C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=4
  1121. infell(9)=6
  1122. IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(9)=9
  1123. GOTO 300
  1124. C Element ICQ4
  1125. 70 infell(2)=4
  1126. infell(3)=4
  1127. IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(3)=9
  1128. infell(4)=4
  1129. C** infell(4)=5 CORRECTION MILL LE 8/8/90
  1130. infell(6)=4
  1131. C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=9
  1132. infell(9)=8
  1133. IF (IFOUR.EQ.-3) infell(9)=11
  1134. IF (IFOUR.EQ.1) infell(9)=12
  1135. GOTO 300
  1136. C Element ICT6
  1137. 71 IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN
  1138. infell(2)=7
  1139. infell(4)=7
  1140. infell(6)=7
  1141. ELSE
  1142. infell(2)=4
  1143. infell(4)=4
  1144. infell(6)=4
  1145. ENDIF
  1146. infell(3)=7
  1147. infell(9)=12
  1148. IF (IFOUR.EQ.-3) infell(9)=15
  1149. IF (IFOUR.EQ.1) infell(9)=18
  1150. GOTO 300
  1151. C Element ICQ8
  1152. 72 infell(2)=9
  1153. infell(3)=9
  1154. IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(3)=16
  1155. infell(4)=9
  1156. infell(6)=9
  1157. C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=16
  1158. infell(9)=16
  1159. IF (IFOUR.EQ.-3) infell(9)=19
  1160. IF (IFOUR.EQ.1) infell(9)=24
  1161. GOTO 300
  1162. C Element ICC8
  1163. 73 infell(2)=8
  1164. infell(3)=8
  1165. infell(4)=8
  1166. infell(6)=8
  1167. infell(9)=24
  1168. GOTO 300
  1169. C Element ICT4
  1170. 74 infell(2)=4
  1171. infell(3)=5
  1172. infell(4)=4
  1173. infell(6)=4
  1174. infell(9)=12
  1175. GOTO 300
  1176. C Element ICP6
  1177. 75 infell(2)=2
  1178. infell(3)=6
  1179. infell(4)=2
  1180. infell(6)=2
  1181. infell(9)=18
  1182. GOTO 300
  1183. C Element IC20
  1184. 76 infell(2)=27
  1185. infell(3)=27
  1186. infell(4)=27
  1187. infell(6)=27
  1188. infell(9)=60
  1189. GOTO 300
  1190. C Element IC10
  1191. 77 infell(2)=5
  1192. infell(3)=5
  1193. infell(4)=5
  1194. infell(6)=5
  1195. infell(9)=30
  1196. GOTO 300
  1197. C Element IC15
  1198. 78 infell(2)=12
  1199. infell(3)=12
  1200. infell(4)=12
  1201. infell(6)=12
  1202. infell(9)=45
  1203. GOTO 300
  1204. C Element TRIP
  1205. 79 IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN
  1206. infell(2)=7
  1207. infell(4)=7
  1208. infell(6)=7
  1209. ELSE
  1210. infell(2)=4
  1211. infell(4)=4
  1212. infell(6)=4
  1213. ENDIF
  1214. infell(3)=7
  1215. infell(8)=9
  1216. infell(9)=15
  1217. IF (IFOUR.EQ.-3.OR.IFOUR.EQ.1) infell(9)=21
  1218. GOTO 300
  1219. C Element QUAP
  1220. 80 infell(2)=9
  1221. infell(3)=9
  1222. IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(3)=9
  1223. infell(4)=9
  1224. infell(6)=9
  1225. C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=16
  1226. infell(8)=12
  1227. infell(9)=20
  1228. IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(9)=28
  1229. GOTO 300
  1230. C Element CUBP
  1231. 81 infell(2)=27
  1232. infell(3)=27
  1233. infell(4)=27
  1234. infell(6)=27
  1235. infell(8)=28
  1236. infell(9)=68
  1237. GOTO 300
  1238. C Element TETP
  1239. 82 infell(2)=5
  1240. infell(3)=5
  1241. infell(4)=5
  1242. infell(6)=5
  1243. infell(8)=14
  1244. infell(9)=34
  1245. GOTO 300
  1246. C Element PRIP
  1247. 83 infell(2)=12
  1248. infell(3)=12
  1249. infell(4)=12
  1250. infell(6)=12
  1251. infell(8)=21
  1252. infell(9)=51
  1253. GOTO 300
  1254. C Element TIMO : Element POUTRE de TIMOSHENCKO
  1255. 84 infell(2)=1
  1256. infell(4)=1
  1257. infell(6)=1
  1258. infell(7)=200
  1259. IF (IFOUR.EQ.-3.OR.IFOUR.EQ.-2.OR.IFOUR.EQ.-1) THEN
  1260. C Element TIMO 2D plan
  1261. infell(5)=6
  1262. infell(9)=6
  1263. infell(10)=3
  1264. infell(15)=3
  1265. infell(16)=3
  1266. ELSE
  1267. C Element TIMO 3D
  1268. infell(5)=12
  1269. infell(9)=12
  1270. infell(10)=6
  1271. infell(15)=6
  1272. infell(16)=6
  1273. ENDIF
  1274. GOTO 300
  1275. C Element JOI2
  1276. 85 infell(2)=4
  1277. infell(3)=2
  1278. infell(4)=2
  1279. infell(5)=0
  1280. infell(6)=2
  1281. infell(7)=300
  1282. infell(9)=8
  1283. GOTO 300
  1284. C Element JOI3
  1285. 86 infell(2)=3
  1286. infell(3)=2
  1287. infell(4)=3
  1288. infell(5)=0
  1289. infell(6)=3
  1290. infell(7)=100
  1291. infell(9)=12
  1292. GOTO 300
  1293. C Element JOT3
  1294. 87 infell(2)=3
  1295. infell(3)=3
  1296. infell(4)=3
  1297. infell(5)=0
  1298. infell(6)=3
  1299. infell(7)=100
  1300. infell(9)=18
  1301. GOTO 300
  1302. C Element JOI4
  1303. 88 infell(2)=4
  1304. infell(3)=4
  1305. infell(4)=4
  1306. infell(5)=0
  1307. infell(6)=4
  1308. infell(7)=100
  1309. infell(9)=24
  1310. GOTO 300
  1311. C Element JOI6
  1312. 89 infell(2)=5
  1313. infell(4)=5
  1314. infell(6)=4
  1315. infell(9)=36
  1316. GOTO 300
  1317. C Element JOI8
  1318. 90 infell(2)=9
  1319. infell(4)=9
  1320. infell(6)=9
  1321. infell(9)=48
  1322. GOTO 300
  1323. C Element LINE SPRING LISC
  1324. 91 infell(2)=3
  1325. infell(4)=3
  1326. infell(5)=5
  1327. infell(6)=5
  1328. infell(7)=100
  1329. infell(9)=36
  1330. infell(10)=2
  1331. infell(15)=6
  1332. infell(16)=6
  1333. GOTO 300
  1334. C Element TRIH
  1335. 92 infell(2)=4
  1336. infell(6)=4
  1337. infell(7)= 0
  1338. IF (IFOUR.EQ.1.OR.IFOUR.EQ.0) THEN
  1339. infell(3)=7
  1340. infell(5)=5
  1341. IF (IFOUR.EQ.0) infell(9)=12
  1342. IF (IFOUR.EQ.1) infell(9)=18
  1343. infell(15)=6
  1344. ELSE
  1345. infell(3)=4
  1346. infell(5)=5
  1347. infell(9)=12
  1348. infell(15)=4
  1349. ENDIF
  1350. GOTO 300
  1351. C Element DST
  1352. 93 infell(2)=3
  1353. infell(3)=3
  1354. infell(4)=3
  1355. infell(5)=3
  1356. infell(6)=3
  1357. infell(7)=600
  1358. infell(9)=18
  1359. infell(10)=8
  1360. infell(15)=6
  1361. infell(16)=8
  1362. GOTO 300
  1363. C Element LIC4 : Element de RACCORD LIQUIDE-COQUE a 8 noeuds en 3D
  1364. 94 infell(2)=4
  1365. infell(3)=4
  1366. infell(4)=4
  1367. infell(5)=3
  1368. infell(6)=4
  1369. infell(9)=32
  1370. infell(15)=6
  1371. GOTO 300
  1372. C Element CERCE
  1373. 95 infell(2)=1
  1374. infell(3)=1
  1375. infell(4)=1
  1376. infell(5)=1
  1377. infell(6)=1
  1378. infell(7)=50
  1379. infell(9)=2
  1380. IF (IFOUR.EQ.1) infell(9)=3
  1381. infell(10)=1
  1382. infell(15)=2
  1383. IF (IFOUR.EQ.1) infell(15)=3
  1384. infell(16)=1
  1385. GOTO 300
  1386. C Element TUYO
  1387. 96 NTETA=40
  1388. infell(2)=2*NTETA
  1389. infell(3)=2*NTETA
  1390. infell(4)=2*NTETA
  1391. infell(5)=7
  1392. infell(6)=2*NTETA
  1393. infell(7)=500
  1394. infell(8)=2+4
  1395. infell(9)=12+2*9
  1396. infell(10)=6
  1397. infell(15)=6+9
  1398. infell(16)=6
  1399. GOTO 300
  1400. C Element LSE2 : TUYAU ACOUSTIQUE PURE
  1401. 97 infell(2)=2
  1402. infell(3)=2
  1403. infell(4)=2
  1404. infell(5)=7
  1405. infell(6)=2
  1406. infell(7)=200
  1407. infell(9)=4
  1408. infell(10)=2
  1409. infell(15)=2
  1410. infell(16)=2
  1411. GOTO 300
  1412. C Element LITU : COUPLAGE LIQUIDE TUYAU
  1413. 98 infell(2)=2
  1414. infell(3)=2
  1415. infell(4)=2
  1416. infell(5)=6
  1417. infell(6)=2
  1418. infell(7)=200
  1419. infell(9)=10
  1420. infell(10)=2
  1421. infell(15)=5
  1422. infell(16)=5
  1423. GOTO 300
  1424. C Element HYT3 : Element HYBRIDE TRI3
  1425. 99 infell(3)=4
  1426. infell(6)=4
  1427. GOTO 300
  1428. C Element HYQ4 : Element HYBRIDE QUA4
  1429. 100 infell(3)=4
  1430. infell(6)=4
  1431. GOTO 300
  1432. C Element HYBRIDE TET4
  1433. 101 infell(3)=4
  1434. infell(6)=4
  1435. GOTO 300
  1436. C Element HYBRIDE PRI6
  1437. 102 infell(3)=6
  1438. infell(6)=6
  1439. GOTO 300
  1440. C Element HYBRIDE CUB8
  1441. 103 infell(3)=8
  1442. infell(6)=8
  1443. GOTO 300
  1444. C Element TRIS (TRI3-SECTION)
  1445. 104 infell(2)=4
  1446. infell(3)=4
  1447. infell(4)=4
  1448. infell(5)=2
  1449. infell(6)=4
  1450. IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.-3) THEN
  1451. infell(10)=2
  1452. infell(16)=2
  1453. ELSE IF (IFOUR.EQ.2) THEN
  1454. infell(10)=3
  1455. infell(16)=3
  1456. ENDIF
  1457. GOTO 300
  1458. C Element QUAS (QUA4-SECTION)
  1459. 105 infell(2)=4
  1460. infell(3)=4
  1461. infell(4)=4
  1462. infell(5)=2
  1463. infell(6)=4
  1464. IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.-3) THEN
  1465. infell(10)=2
  1466. infell(16)=2
  1467. ELSE IF (IFOUR.EQ.2) THEN
  1468. infell(10)=3
  1469. infell(16)=3
  1470. ENDIF
  1471. GOTO 300
  1472. C Element POIS (POI1-SECTION)
  1473. 106 infell(2)=1
  1474. infell(3)=1
  1475. infell(4)=1
  1476. infell(5)=1
  1477. CMILL???? infell(5)=3
  1478. infell(6)=1
  1479. infell(10)=1
  1480. infell(16)=1
  1481. GOTO 300
  1482. C Element JOP3 : JOINTS POREUX
  1483. 108 infell(2)=3
  1484. infell(3)=3
  1485. infell(4)=3
  1486. infell(6)=3
  1487. infell(8)=12
  1488. infell(9)=18
  1489. infell(10)=2
  1490. infell(15)=3
  1491. infell(16)=3
  1492. GOTO 300
  1493. C Element JOP6 : JOINTS POREUX
  1494. 109 infell(2)=7
  1495. infell(3)=7
  1496. infell(4)=7
  1497. infell(6)=7
  1498. infell(8)=21
  1499. infell(9)=45
  1500. infell(10)=3
  1501. infell(15)=4
  1502. infell(16)=4
  1503. GOTO 300
  1504. C Element JOP8 : JOINTS POREUX
  1505. 110 infell(2)=9
  1506. infell(3)=9
  1507. infell(4)=9
  1508. infell(6)=9
  1509. infell(8)=28
  1510. infell(9)=60
  1511. infell(10)=3
  1512. infell(15)=4
  1513. infell(16)=4
  1514. GOTO 300
  1515. C Elements POLYGONE a N cotes
  1516. 111 NBSH=infell(8)
  1517. infell(2)=NBSH
  1518. infell(3)=NBSH
  1519. infell(4)=NBSH
  1520. infell(6)=NBSH
  1521. infell(9)=2*NBSH
  1522. IF (IFOUR.EQ.1) infell(9)=3*NBSH
  1523. IF (IFOUR.EQ.-3) infell(9)=2*NBSH+3
  1524. GOTO 300
  1525. C Element BAR3 (barre a 3 noeuds)
  1526. 123 NPOINT=3
  1527. infell(2)=NPOINT
  1528. infell(4)=NPOINT
  1529. infell(5)=1
  1530. infell(6)=NPOINT
  1531. infell(7)=200
  1532. infell(9)=9
  1533. IF (IFOUR.EQ.-1.OR.IFOUR.EQ.-2.OR.IFOUR.EQ.0) infell(9)=6
  1534. infell(10)=1
  1535. infell(15)=3
  1536. IF (IFOUR.EQ.-1.OR.IFOUR.EQ.-2.OR.IFOUR.EQ.0) infell(15)=2
  1537. infell(16)=1
  1538. GOTO 300
  1539. C Element BAEX (barre excentree 3D a 2 noeuds)
  1540. 124 infell(2)=2
  1541. infell(3)=2
  1542. infell(4)=2
  1543. infell(5)=6
  1544. infell(6)=2
  1545. infell(7)=700
  1546. infell(9)=12
  1547. infell(10)=1
  1548. infell(15)=6
  1549. infell(16)=1
  1550. GOTO 300
  1551. C Element LIA2 : element de liaison a 2 noeuds (6 ddl par noeuds)
  1552. 125 infell(2)=2
  1553. infell(4)=2
  1554. infell(5)=9
  1555. infell(6)=2
  1556. infell(7)=200
  1557. infell(9)=12
  1558. infell(10)=6
  1559. infell(15)=6
  1560. infell(16)=6
  1561. GOTO 300
  1562. C Element QUAH
  1563. 126 infell(2)=4
  1564. infell(3)=4
  1565. IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(3)=9
  1566. infell(5)=5
  1567. infell(6)=4
  1568. infell(7)=0
  1569. IF (IFOUR.EQ.0) infell(9)=16
  1570. IF (IFOUR.EQ.1) infell(9)=24
  1571. infell(15)=6
  1572. GOTO 300
  1573. C Element CUBH
  1574. 127 infell(2)=8
  1575. infell(3)=27
  1576. infell(5)=5
  1577. infell(6)=27
  1578. infell(7)=0
  1579. infell(9)=48
  1580. infell(15)=6
  1581. GOTO 300
  1582. C Element ROT3 (1 pt GAUSS pr MASSE, 7 pts GAUSS pr RIGIDITE)
  1583. 128 infell(3)=1
  1584. infell(6)=7
  1585. infell(9)=3
  1586. infell(15)=1
  1587. GOTO 300
  1588. C Element TR6H
  1589. 157 infell(2)=4
  1590. infell(6)=6
  1591. infell(7)=0
  1592. IF (IFOUR.EQ.1.OR.IFOUR.EQ.0) THEN
  1593. infell(3)=7
  1594. infell(5)=5
  1595. IF (IFOUR.EQ.0) infell(9)=12
  1596. IF (IFOUR.EQ.1) infell(9)=18
  1597. infell(15)=6
  1598. ELSE
  1599. infell(3)=6
  1600. infell(5)=5
  1601. infell(9)=24
  1602. infell(15)=4
  1603. ENDIF
  1604. GOTO 300
  1605. C Element SEGS (SEG2-SECTION 2D)
  1606. 166 infell(2)=2
  1607. infell(3)=2
  1608. infell(4)=2
  1609. infell(5)=3
  1610. infell(6)=2
  1611. IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.-3) THEN
  1612. infell(10)=2
  1613. infell(16)=2
  1614. ELSE IF (IFOUR.EQ.2) THEN
  1615. infell(10)=3
  1616. infell(16)=3
  1617. ENDIF
  1618. GOTO 300
  1619. C Element POJS (POI1-SECTION)
  1620. 167 infell(2)=1
  1621. infell(3)=1
  1622. infell(4)=1
  1623. infell(5)=3
  1624. infell(6)=1
  1625. IF (IFOUR.EQ.-2.OR.IFOUR.EQ.-1.OR.IFOUR.EQ.-3) THEN
  1626. infell(10)=2
  1627. infell(16)=2
  1628. ELSE IF (IFOUR.EQ.2) THEN
  1629. infell(10)=3
  1630. infell(16)=3
  1631. ENDIF
  1632. GOTO 300
  1633. C Element JCT3
  1634. 168 infell(2)=3
  1635. infell(3)=3
  1636. infell(4)=3
  1637. infell(6)=3
  1638. GOTO 300
  1639. C Element JCI4
  1640. 169 infell(2)=4
  1641. infell(3)=4
  1642. infell(4)=4
  1643. infell(6)=4
  1644. GOTO 300
  1645. C Element JGI2
  1646. 170 infell(2)=4
  1647. infell(3)=2
  1648. infell(4)=2
  1649. infell(6)=2
  1650. IF (IFOUR.EQ.-3) infell(9)=infell(9)+3
  1651. infell(15)=2
  1652. GOTO 300
  1653. C Element JGT3
  1654. 171 infell(2)=3
  1655. infell(3)=3
  1656. infell(4)=3
  1657. infell(6)=3
  1658. GOTO 300
  1659. C Element JGI4
  1660. 172 infell(2)=4
  1661. infell(3)=4
  1662. infell(4)=4
  1663. infell(6)=4
  1664. GOTO 300
  1665. C Element TRIQ
  1666. 173 infell(3)=7
  1667. IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN
  1668. infell(2)=7
  1669. infell(4)=7
  1670. infell(6)=7
  1671. ELSE
  1672. infell(2)=4
  1673. infell(4)=4
  1674. infell(6)=4
  1675. ENDIF
  1676. infell(8)=9
  1677. infell(9)=18
  1678. IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(9)=24
  1679. GOTO 300
  1680. C Element QUAQ
  1681. 174 infell(2)=9
  1682. infell(3)=9
  1683. IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(3)=9
  1684. infell(4)=9
  1685. infell(6)=9
  1686. C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=16
  1687. infell(8)=12
  1688. infell(9)=24
  1689. IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(9)=32
  1690. GOTO 300
  1691. C Element CUBQ
  1692. 175 infell(2)=27
  1693. infell(3)=27
  1694. infell(4)=27
  1695. infell(6)=27
  1696. infell(8)=28
  1697. infell(9)=76
  1698. GOTO 300
  1699. C Element TETQ
  1700. 176 infell(2)=5
  1701. infell(3)=5
  1702. infell(4)=5
  1703. infell(6)=5
  1704. infell(8)=14
  1705. infell(9)=38
  1706. GOTO 300
  1707. C Element PRIQ
  1708. 177 infell(2)=12
  1709. infell(3)=12
  1710. infell(4)=12
  1711. infell(6)=12
  1712. infell(8)=21
  1713. infell(9)=57
  1714. GOTO 300
  1715. C Element TRIR
  1716. 178 IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) THEN
  1717. infell(2)=7
  1718. infell(6)=7
  1719. infell(4)=7
  1720. ELSE
  1721. infell(2)=4
  1722. infell(4)=4
  1723. infell(6)=4
  1724. ENDIF
  1725. infell(3)=7
  1726. infell(8)=9
  1727. infell(9)=21
  1728. IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(9)=27
  1729. GOTO 300
  1730. C Element QUAR
  1731. 179 infell(2)=9
  1732. infell(3)=9
  1733. IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(3)=9
  1734. infell(4)=9
  1735. infell(6)=9
  1736. C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=16
  1737. infell(8)=12
  1738. infell(9)=28
  1739. IF (IFOUR.EQ.1.OR.IFOUR.EQ.-3) infell(9)=36
  1740. GOTO 300
  1741. C Element CUBR
  1742. 180 infell(2)=27
  1743. infell(3)=27
  1744. infell(4)=27
  1745. infell(6)=27
  1746. infell(8)=28
  1747. infell(9)=84
  1748. GOTO 300
  1749. C Element TETR
  1750. 181 infell(2)=5
  1751. infell(3)=5
  1752. infell(4)=5
  1753. infell(6)=5
  1754. infell(8)=14
  1755. infell(9)=42
  1756. GOTO 300
  1757. C Element PRIR
  1758. 182 infell(2)=12
  1759. infell(3)=12
  1760. infell(4)=12
  1761. infell(6)=12
  1762. infell(8)=21
  1763. infell(9)=61
  1764. GOTO 300
  1765. C Element Q4RI (QUA4 with 1x1 Gauss points)
  1766. 183 infell(2)=4
  1767. infell(3)=1
  1768. IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(3)=9
  1769. infell(4)=1
  1770. C** infell(4)=5 CORRECTION MILL LE 8/8/90
  1771. infell(6)=1
  1772. C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=9
  1773. infell(9)=8
  1774. IF (IFOUR.EQ.-3) infell(9)=11
  1775. IF (IFOUR.EQ.1) infell(9)=12
  1776. GOTO 300
  1777. C Element Q8RI (QUA8 with 2x2 Gauss points)
  1778. 184 infell(2)=9
  1779. infell(3)=4
  1780. IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(3)=16
  1781. infell(4)=4
  1782. infell(6)=4
  1783. C** IF (IFOUR.EQ.0.OR.IFOUR.EQ.1) infell(6)=16
  1784. infell(9)=16
  1785. IF (IFOUR.EQ.-3) infell(9)=19
  1786. IF (IFOUR.EQ.1) infell(9)=24
  1787. GOTO 300
  1788. C Element JOQ3
  1789. 185 infell(2)=3
  1790. infell(3)=3
  1791. infell(4)=3
  1792. infell(6)=3
  1793. infell(8)=12
  1794. infell(9)=24
  1795. infell(10)=4
  1796. infell(15)=4
  1797. infell(16)=4
  1798. GOTO 300
  1799. C Element JOQ6
  1800. 186 infell(2)=7
  1801. infell(3)=7
  1802. infell(4)=7
  1803. infell(6)=7
  1804. infell(8)=21
  1805. infell(9)=54
  1806. infell(10)=5
  1807. infell(15)=5
  1808. infell(16)=5
  1809. GOTO 300
  1810. C Element JOQ8
  1811. 187 infell(2)=9
  1812. infell(3)=9
  1813. infell(4)=9
  1814. infell(6)=9
  1815. infell(8)=28
  1816. infell(9)=72
  1817. infell(10)=5
  1818. infell(15)=5
  1819. infell(16)=5
  1820. GOTO 300
  1821. C Element JOR3
  1822. 188 infell(2)=3
  1823. infell(3)=3
  1824. infell(4)=3
  1825. infell(6)=3
  1826. infell(8)=12
  1827. infell(9)=30
  1828. infell(10)=5
  1829. infell(15)=5
  1830. infell(16)=5
  1831. GOTO 300
  1832. C Element JOR6
  1833. 189 infell(2)=7
  1834. infell(3)=7
  1835. infell(4)=7
  1836. infell(6)=7
  1837. infell(8)=21
  1838. infell(9)=63
  1839. infell(10)=6
  1840. infell(15)=6
  1841. infell(16)=6
  1842. GOTO 300
  1843. C Element JOR8
  1844. 190 infell(2)=9
  1845. infell(3)=9
  1846. infell(4)=9
  1847. infell(6)=9
  1848. infell(8)=28
  1849. infell(9)=84
  1850. infell(10)=6
  1851. infell(15)=6
  1852. infell(16)=6
  1853. GOTO 300
  1854. C Element T1D2/M1D2 : Element THERMIQUE/MECANIQUE MASSIF 1D a 2 noeuds
  1855. C On suppose que IFOMOD= 3 ou 4 ou 5 (IDIM=1)
  1856. 191 infell(2)=2
  1857. infell(3)=2
  1858. infell(4)=2
  1859. infell(6)=2
  1860. infell(9)=2
  1861. IF ((IFOUR.GE.7.AND.IFOUR.LE.10).OR.IFOUR.EQ.14) infell(9)=3
  1862. IF (IFOUR.EQ.11) infell(9)=4
  1863. GOTO 300
  1864. C Element T1D3/M1D3 : Element THERMIQUE/MECANIQUE MASSIF 1D a 3 noeuds
  1865. C On suppose que IFOMOD= 3 ou 4 ou 5 (IDIM=1)
  1866. 192 infell(2)=3
  1867. infell(3)=3
  1868. infell(4)=3
  1869. infell(6)=3
  1870. infell(9)=3
  1871. IF ((IFOUR.GE.7.AND.IFOUR.LE.10).OR.IFOUR.EQ.14) infell(9)=4
  1872. IF (IFOUR.EQ.11) infell(9)=5
  1873. GOTO 300
  1874. C Element fluide NAVIER_STOKES
  1875. C LC03 LS03 BS03
  1876. 195 infell(6)=2
  1877. GOTO 300
  1878. C LC07 LS07 BS07
  1879. 196 infell(6)=7
  1880. GOTO 300
  1881. C LC09 LS09 BS09
  1882. 197 infell(6)=4
  1883. GOTO 300
  1884. C MC03
  1885. 216 infell(6)=2
  1886. GOTO 300
  1887. C MC07
  1888. 217 infell(6)=7
  1889. GOTO 300
  1890. C MC09
  1891. 218 infell(6)=4
  1892. GOTO 300
  1893. C M103 MS03
  1894. 223 infell(6)=3
  1895. GOTO 300
  1896. C M107 MS07
  1897. 224 infell(6)=7
  1898. GOTO 300
  1899. C M109 MS09
  1900. 225 infell(6)=16
  1901. GOTO 300
  1902. C QC03 Q103 QS03
  1903. 237 infell(6)=3
  1904. GOTO 300
  1905. C QC07 Q107 QS07
  1906. 238 infell(6)=7
  1907. GOTO 300
  1908. C QC09 Q109 QS09
  1909. 239 infell(6)=16
  1910. GOTO 300
  1911. C
  1912. C CIFL
  1913. C
  1914. 258 CONTINUE
  1915. infell(2)=1
  1916. infell(3)=1
  1917. infell(4)=1
  1918. infell(6)=1
  1919. infell(8)=0
  1920. infell(9)=8
  1921. infell(10)=5
  1922. infell(15)=3
  1923. infell(16)=5
  1924. GO TO 300
  1925. C
  1926. C SHB8
  1927. C
  1928. 260 CONTINUE
  1929. infell(2)=5
  1930. infell(3)=5
  1931. infell(4)=5
  1932. infell(6)=5
  1933. infell(8)=0
  1934. infell(9)=24
  1935. infell(10)=6
  1936. infell(13)= 1
  1937. infell(15)=3
  1938. infell(16)=6
  1939. GO TO 300
  1940. C
  1941. C XQ4R
  1942. C
  1943. 263 INFELl(2)=0
  1944. INFELl(3)=64
  1945. INFELl(4)=64
  1946. INFELl(6)=64
  1947. INFELl(15)=20
  1948. INFELl(16)=4
  1949. INFELl(9) = 4*INFELl(15)
  1950. INFELl(10)=4
  1951. GOTO 300
  1952. C
  1953. C XC8R
  1954. C
  1955. 264 INFELl(2)=0
  1956. INFELl(3)=64
  1957. INFELl(4)=64
  1958. INFELl(6)=64
  1959. INFELl(15)=30
  1960. INFELl(16)=6
  1961. INFELl(9) = 8*INFELl(15)
  1962. INFELl(10)=6
  1963. GOTO 300
  1964. C
  1965. C Element JOI1 : element de liaison a 2 noeuds (6 ddl par noeuds)
  1966. C
  1967. 265 infell(2)=1
  1968. infell(3)=1
  1969. infell(4)=1
  1970. infell(5)=0
  1971. infell(6)=1
  1972. infell(7)=700
  1973. infell(8)=0
  1974. infell(13)=75
  1975. infell(14)=2
  1976. IF(IFOUR.EQ.2) THEN
  1977. infell(9)=12
  1978. infell(10)=6
  1979. infell(15)=6
  1980. infell(16)=6
  1981. ELSE IF(IFOUR.EQ.-2.OR.IFOUR.EQ.-1) THEN
  1982. infell(9)=6
  1983. infell(10)=3
  1984. infell(15)=3
  1985. infell(16)=3
  1986. ENDIF
  1987. GOTO 300
  1988. C
  1989. C Element ZCO2 : zone cohesive a 2 noeuds (2 ddls par noeud)
  1990. C
  1991. 266 infell(2)=2
  1992. infell(4)=2
  1993. infell(6)=2
  1994. infell(7)=100
  1995. infell(9)=4
  1996. infell(10)=2
  1997. infell(15)=2
  1998. infell(16)=2
  1999. GOTO 300
  2000. C
  2001. C Element ZCO3 : zone cohesive a 3 noeuds (2 ddls par noeud)
  2002. C
  2003. 267 infell(2)=1
  2004. infell(4)=1
  2005. infell(6)=1
  2006. infell(7)=100
  2007. infell(9)=9
  2008. infell(10)=3
  2009. infell(15)=3
  2010. infell(16)=3
  2011. GOTO 300
  2012. C
  2013. C Element ZCO4 : zone cohesive a 4 noeuds (2 ddls par noeud)
  2014. C
  2015. 268 infell(2)=4
  2016. infell(4)=4
  2017. infell(6)=4
  2018. infell(7)=100
  2019. infell(9)=12
  2020. infell(10)=3
  2021. infell(15)=3
  2022. infell(16)=3
  2023. GOTO 300
  2024. c
  2025. C Element TUY2 : tuyau pour modele d'advection en thermique
  2026. C
  2027. 269 infell(2)=2
  2028. infell(4)=2
  2029. infell(5)=5
  2030. infell(6)=2
  2031. infell(7)=200
  2032. infell(9)=2
  2033. infell(10)=1
  2034. infell(15)=1
  2035. infell(16)=1
  2036. GOTO 300
  2037. c
  2038. C Element TUY3 : tuyau pour modele d'advection en thermique
  2039. C
  2040. 270 infell(2)=3
  2041. infell(4)=3
  2042. infell(5)=5
  2043. infell(6)=3
  2044. infell(7)=200
  2045. infell(9)=2
  2046. infell(10)=1
  2047. infell(15)=1
  2048. infell(16)=1
  2049. GOTO 300
  2050. c
  2051. c element coaxial COS2 (3D pour liaison acier-beton)
  2052. c
  2053. 271 infell(2)=1
  2054. infell(4)=2
  2055. infell(5)=2
  2056. infell(6)=2
  2057. infell(7)=100
  2058. infell(9)=4*idim
  2059. infell(10)=3
  2060. infell(14)=12
  2061. infell(15)=IDIM
  2062. infell(16)=3
  2063. c
  2064. c element coaxial COA2 (3D pour liaison acier-beton)
  2065. c
  2066. 272 infell(2)=1
  2067. infell(4)=2
  2068. infell(5)=2
  2069. infell(6)=2
  2070. infell(7)=100
  2071. infell(9)=4*idim
  2072. infell(10)=3
  2073. infell(14)=12
  2074. infell(15)=IDIM
  2075. infell(16)=3
  2076. C
  2077. 300 IF (IELE.EQ.0) IELE=infell(14)
  2078. NBGST = infell(4)
  2079. NBG = infell(6)
  2080. NBGM = infell(3)
  2081. IF (NBGM.EQ.0) THEN
  2082. NBGM=NBG
  2083. infell(3)=NBGM
  2084. ENDIF
  2085. NBSH = NBNNE(IELE)
  2086. IF (MELE.GE.79.AND.MELE.LE.83 ) NBSH = infell(8)
  2087. IF (MELE.EQ.96) NBSH = infell(8)
  2088. IF (MELE.GE.108.AND.MELE.LE.110) NBSH = infell(8)
  2089. IF (MELE.GE.111.AND.MELE.LE.122) NBSH = infell(8)
  2090. IF (MELE.GE.173.AND.MELE.LE.182) NBSH = infell(8)
  2091. IF (MELE.GE.185.AND.MELE.LE.190) NBSH = infell(8)
  2092. IF (MELE.GE.195.AND.MELE.LE.257) NBSH = infell(8)
  2093. IF (MELE.EQ.263.OR.MELE.EQ.264 ) NBSH = infell(8)
  2094.  
  2095. C Cas particulier des formulations massives (mfr2 > 0)
  2096. C Ici : infell(9) = 1 * nb.noeuds.EF
  2097. IF (mfr2.NE.0) THEN
  2098. infell(9) = NBSH
  2099. ENDIF
  2100. *
  2101. * nombre composante contraintes/defo, materiau IMPEDANCE
  2102. if (dcmate) then
  2103. call idprim(IMODEL,MFR3,IPNOMC,NBROBL,NBRFAC)
  2104. if(tymode(/2).gt.0) then
  2105. if (tymode(1).eq.'LISTMOTS') then
  2106. dcmat2 = .true.
  2107. mlmot5 = ivamod(1)
  2108. segact mlmot5
  2109. nbrobl = mlmot5.mots(/2)*2
  2110. endif
  2111. endif
  2112. if (dcmat2.and.mele.eq.2) then
  2113. infell(9) = nbrobl
  2114. infell(10) = nbrobl/2
  2115. infell(15) = nbrobl/2
  2116. infell(16) = nbrobl/2
  2117. else
  2118. infell(15) = nbrobl
  2119. if (nbrfac.gt.0) infell(15) = infell(15) + nbrfac
  2120. infell(16) = infell(15)
  2121. if (mele.eq.45) infell(9) = nbrobl
  2122. if (mele.eq.2) infell(9) = nbrobl*2
  2123. endif
  2124. endif
  2125.  
  2126. c* write(6,*)'TYINTE',TYINTE
  2127. c* write(6,*)'MELE:',MELE,'IELE:',IELE,'NBSH:',NBSH,'NBG:',NBG
  2128.  
  2129. C Remplissage de infell(12) si MFR=5 cad MELE=41 ou 56 (COQ8 ou COQ6)
  2130. C -------------------------
  2131. IF (MFR.EQ.5.OR.MFR.EQ.74) THEN
  2132. c* write(6,*) ' elquoi appel a renoeu '
  2133. CALL RENOEU(IELE,MELE,NBSH,IPT2,IRT2)
  2134. c* write(6,*) ' elquoi sortie renoeu ipt2 irt2' ,ipt2,irt2
  2135. if(inttyp.eq.-10) then
  2136. infmod(8)=IPT2
  2137. else
  2138. infell(12)=ipt2
  2139. endif
  2140. ENDIF
  2141. c
  2142. c cas XFEM : le 2nd segment d integration contient le MINTE de l'EF
  2143. c correspondant non enrichi (avec un nombre reduit de point de Gauss)
  2144. IF (MFR.EQ.63) THEN
  2145. if(MELE.eq.263) then
  2146. MEL2 = 8
  2147. NBG2 = 4
  2148. NBSH2 = 4
  2149. else if(MELE.eq.264) then
  2150. MEL2 = 14
  2151. NBG2 = 8
  2152. NBSH2 = 8
  2153. endif
  2154. c pour l'instant on ne distingue pas les cas selon TYINTE
  2155. CALL RESHPT(NBG2,NBSH2,IELE,MEL2,0,IPT2,IRT2)
  2156. if(inttyp.eq.-10) then
  2157. infmod(8)=IPT2
  2158. else
  2159. infell(12)=IPT2
  2160. endif
  2161. ENDIF
  2162.  
  2163. C Remplissage de infell(11) : INTEGRATION DE L'ELEMENT FINI
  2164. C ---------------------------
  2165. IDETYI = 0
  2166. 400 CONTINUE
  2167. IF (INTTYP.EQ.-10) TYINTE = TYINTE + 1
  2168.  
  2169. C Tests de verification :
  2170. IF (MELE.GT.NELFI) THEN
  2171. WRITE(IOIMP,*)
  2172. WRITE(IOIMP,*) 'ELQUOI : incoherence NELFI & MELE ',
  2173. & NELFI,MELE
  2174. CALL ERREUR(5)
  2175. RETURN
  2176. ENDIF
  2177. IF (TYINTE.GT.NTYNTE) THEN
  2178. WRITE(IOIMP,*)
  2179. WRITE(IOIMP,*) 'ELQUOI : incoherence TYINTE & NTYNTE',
  2180. & TYINTE,NTYNTE
  2181. CALL ERREUR(5)
  2182. RETURN
  2183. ENDIF
  2184.  
  2185. C == === ====== ===
  2186. C Si le segment n'a pas deja ete rempli ou si c'est un DKT : iin = 0
  2187. C == === ====== ===
  2188. IF (MELE.EQ.28) THEN
  2189. iin = 0
  2190. ELSE
  2191. iin = INTEGR(MELE,TYINTE,IBMODE)
  2192. ENDIF
  2193. C == ==== ======
  2194. C Si le segment d'integration a deja ete rempli : iin est non nul ici !
  2195. C == ==== ======
  2196.  
  2197. IF (iin.EQ.0) THEN
  2198. IPT1 = 0
  2199. C Cas des elements de MECANIQUE :
  2200. C -------------------------------
  2201. C SI FAUX ? Ne devrait-on pas avoir MELE < 57 OU MELE > 68 ?????
  2202. IF (MELE.LT.57.OR.MELE.GT.64) THEN
  2203. C 1 = Champ aux noeuds
  2204. IF (TYINTE.EQ.1) THEN
  2205. IF(MELE.GE.195.and.mele.ne.260.and.mele.ne.263.and.
  2206. & mele.ne.264) THEN
  2207. CALL RESHPT(NBG,NBSH,IELE,MELE,0,IPT1,IRT1)
  2208. ELSE
  2209. CALL RENOEU(IELE,MELE,NBSH,IPT1,IRT1)
  2210. ENDIF
  2211. C 2 = Point de Gauss, centre de gravite et champ CONSTANT
  2212. ELSE IF (TYINTE.EQ.2) THEN
  2213. CALL RESHPT(1,NBSH,IELE,MELE,NPINT,IPT1,IRT1)
  2214. C 3 = Point de Gauss pour la rigidite
  2215. ELSE IF (TYINTE.EQ.3) THEN
  2216. if(mele.eq.263.or.mele.eq.264) then
  2217. CALL RESHPX(NBG,NBSH,IELE,MELE,NPINT,IPT1,IRT1)
  2218. else
  2219. CALL RESHPT(NBG,NBSH,IELE,MELE,NPINT,IPT1,IRT1)
  2220. endif
  2221. C 4 = Point de Gauss pour la masse
  2222. ELSE IF (TYINTE.EQ.4) THEN
  2223. if(mele.eq.263.or.mele.eq.264) then
  2224. CALL RESHPX(NBGM,NBSH,IELE,MELE,NPINT,IPT1,IRT1)
  2225. else
  2226. CALL RESHPT(NBGM,NBSH,IELE,MELE,NPINT,IPT1,IRT1)
  2227. endif
  2228. C 5 = Point de Gauss - calcul des contraintes
  2229. ELSE IF (TYINTE.EQ.5) THEN
  2230. if(mele.eq.263.or.mele.eq.264) then
  2231. CALL RESHPX(NBG,NBSH,IELE,MELE,NPINT,IPT1,IRT1)
  2232. else
  2233. CALL RESHPT(NBGST,NBSH,IELE,MELE,NPINT,IPT1,IRT1)
  2234. endif
  2235. C 6 = Pas de segment d'integration
  2236. ELSE IF (TYINTE.EQ.6) THEN
  2237. IPT1 = 0
  2238. C 7 = Champ aux faces
  2239. C 8 = Champ aux CENTREP1
  2240. ELSE IF (TYINTE.EQ.8) THEN
  2241. CALL RESHPT(NBG,NBSH,IELE,MELE,1,IPT1,IRT1)
  2242. C 9 = Champ aux MSOMMET
  2243. ELSE IF (TYINTE.EQ.9) THEN
  2244. CALL RESHPT(NBG,NBSH,IELE,MELE,2,IPT1,IRT1)
  2245. C Autre cas = Point de Gauss pour la rigidite
  2246. ELSE
  2247. CALL RESHPT(NBG,NBSH,IELE,MELE,NPINT,IPT1,IRT1)
  2248. ENDIF
  2249. C Cas des elements de THERMIQUE :
  2250. C -------------------------------
  2251. C ELSE
  2252. C CALL TINTR(MELE,IPT1)
  2253. ENDIF
  2254. IF (IPT1.NE.0) CALL SAVSEG(IPT1)
  2255. iin = IPT1
  2256. INTEGR(MELE,TYINTE,IBMODE) = iin
  2257. ENDIF
  2258.  
  2259. IF (INTTYP.EQ.-10) THEN
  2260. imodel.infmod(2+TYINTE+IDETYI) = iin
  2261. IF (TYINTE.LT.5) GOTO 400
  2262. IF (TYINTE.EQ.5.AND.IDETYI.EQ.0) THEN
  2263. IF (mele.EQ.263) THEN
  2264. * pour le XQ4R on recommence pour avoir les fonctions de forme
  2265. * classiques du : qua4
  2266. TYINTE = 0
  2267. IDETYI = 5
  2268. mele = 8
  2269. nbsh = 4
  2270. nbg = 4
  2271. nbgm = 4
  2272. nbgst = 4
  2273. GOTO 400
  2274. ELSE IF (mele.EQ.264) THEN
  2275. * pour le XC8R on recommence pour avoir les fonctions de forme
  2276. * classiques du : cub8
  2277. TYINTE = 0
  2278. IDETYI = 5
  2279. mele = 14
  2280. nbsh = 8
  2281. nbg = 8
  2282. nbgm = 8
  2283. nbgst = 8
  2284. GOTO 400
  2285. ENDIF
  2286. ENDIF
  2287. ENDIF
  2288.  
  2289. MINTE=iin
  2290. IF(MINTE.GT.0) SEGACT,MINTE
  2291. infell(11) = iin
  2292.  
  2293. C Sortie de ELQUOI : le segment IPTR=INFO est ACTIF.
  2294. IPTR = INFO
  2295. c* write(6,*)'-> ELQUOI infell',(infell(iou),iou=1,16)
  2296.  
  2297. END
  2298.  
  2299.  
  2300.  

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