Télécharger elquoi.eso

Retour à la liste

Numérotation des lignes :

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

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