7.4.4.4 Données des éléments

lecture dans S_PT2

READ (LECELT,1002) M,IEL,IPS,MTYP,EPAIS,KG,(NODB(I), I= 1,IEL/2)

FORMAT(4I5,F10.0,4I5)

READ (LECELT,1003) (NOD(I),I=1,IEL )

FORMAT(6I5)

note	colonnes	variable	définition
1ère ligne:
1	1-5	M	numéro de l'élément $1 \leq \texttt{M} \leq \texttt{NPAR(2)}$
2	6-10	IEL	nombre des n?uds sur le bord (n?uds esclaves) de cet élément = 4, pour les éléments linéaires = 6, pour les éléments quadratiques
	11-15	IPS	code pour l'impression des contraintes = 0, pas d'impression
	16-20	MTYP	numéro du matériau de l'élément $1 \leq \mbox{\texttt{MTYP}} \leq \mbox{\texttt{NPAR(16)}}$
	21-30	EPASI	épaisseur de l'élément dans la direction globale $x$ (direction locale $e_3$)
	31-35	KG	incrément pour la génération automatique
3	36-40	NODB(1)	noeud $P_1$ sur la fibre moyenne (local)
	41-45	NODB(2)	noeud $P_2$ sur la fibre moyenne
	46-50	NODB(2)	noeud $P_3$ sur la fibre moyenne
2ème ligne:
4	1-5	NOD(1)	noeud esclave 1 (local)
	1-5	NOD(1)	noeud esclave 1 (local)
	6-10	NOD(2)	noeud esclave 2 (local)
	11-15	NOD(3)	noeud esclave 3 (local)
	16-20	NOD(4)	noeud esclave 4 (local)
	21-25	NOD(5)	noeud esclave 5 (local)
	26-30	NOD(6)	noeud esclave 6 (local)

notes

1. Il faut absolument donner le dernier élément du groupe, c'est-à-dire :
   M=NUME + INCELT
   Génération des éléments:
   Si les lignes des éléments ($m+1$,$m+2$,.....,$m+j$) manquent, les $j$ éléments sont générés en prenant IPS et MTYP de l'élément M et leurs noeuds sont obtenus à partir des noeuds donnés sur la ligne de l'élément M incrémentés de KG. Seuls les IEL noeuds non nuls sont incrémentés.
2. On donne ici le nombre des n?uds esclaves (sur le bord de l'élément). Cette valeur est égale a deux fois le nombre des n?uds sur la fibre moyenne
3. La numérotation des n?uds sur la fibre moyenne suit la convention de la figure. Pour l'élément linéaire donner 0 ou mettre rien pour le noeud $P_3$.
4. Aussi la numérotation des n?uds esclaves suit la convention de la figure. Pour l'élément linéaire donner 0 ou mettre rien pour les noeuds 5 et 6.
   Attention: Dans la définition globale des n?uds les degrés de liberté $u_y$ et $u_z$ des n?uds esclaves doivent être imposés (code -1), leurs rotations bloquées (code 1).