Télécharger elquoi.eso

Retour à la liste

Numérotation des lignes :

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

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