7.6.4.1 Caractéristiques du groupe

lecture dans IELMNT

READ (LECG,1000) NPAR

FORMAT(20I4)

note	colonnes	variable	définition
	1-4	NPAR(1)	donner le nombre 9
	5-8	NPAR(2)	(NUME) nombre d'éléments dans ce groupe $\texttt{NUME} \> 0$
1	9-12	NPAR(3)	(INDNL) code indiquant le type d'analyse non linéaire = 0, mis par défaut a "1" = 1, matériau non linéaire seulement
	13-16	NPAR(4)	(IPARX) ordre d'élément paraxial = 0, uniquement
	17-20	NPAR(5)	inutilisé
2	21-24	NPAR(6)	(NDIM) nombre de "ddl locaux" des noeuds des éléments = 3, ddl $u_x$, $u_y$ et $u_z$ (mécanique seulement) = 0, mis par défaut à 3
3	25-28	NPAR(7)	(MXNODS) nombre maximum de noeuds utilisés pour décrire un élément du groupe: =4 à 8 pour les quadrangles =3 à 6 pour les triangles mis par défaut à 8
	29-32	NPAR(8)	(NTYPEL) type d'élément géométrique = 1, ce groupe contient des quadrangles seulement = 2, ce groupe contient des triangles seulement = 0, ce groupe peut contenir les deux types d'éléments
	33-36	NPAR(9)	inutilisé
4	37-40	NPAR(10)	(NINT) ordre d'intégration numérique pas de valeur par défaut $1 \leq \mbox{\texttt{NINT}} \leq 4$
	41-44	NPAR(11)	inutilisé
	45-48	NPAR(12)	inutilisé
	49-52	NPAR(13)	inutilisé
5	53-56	NPAR(14)	(INCTYP) type de champ incident = 0, pas de champ incident = 1, onde plane harmonique = 2, onde plane signal Ricker = 3, onde plane quelconque = 4, champ quelconque
	57-60	NPAR(15)	(MODEL) numéro de modèle choisi pour le matériau = 1, seulement par défaut mis à 1
	61-64	NPAR(16)	(NUMMAT) nombre de "jeux de propriétés" différents pour le modèle de comportement choisi $\texttt{NUMMAT} \geq 1$ pas de valeur par défaut
	65-68	NPAR(17)	inutilisé
6	69-72	NPAR(18)	(NMBNO) nombre de noeuds paraxiaux (seulement si INCTYP=4)
6	73-76	NPAR(19)	(NMBDT) nombre de pas de temps (seulement si INCTYP=4)
	77-80	NPAR(20)	(IEXPL) code permettant de définir si le groupe est traité implicitement ou explicitement dans un calcul dynamique avec un code d'intégration IOPE=10 (cf. paragraphe 1.7). = 1, groupe explicite (seulement)

notes

1. Dans la version actuelle, donner NPAR(3)= 0 ou 1.
   Seules les non linéarités dues aux matériaux peuvent être modélisées: les déplacements et les déformations sont supposés infiniment petits.
2. Les éléments paraxiaux n'existent qu'en mécanique pure.
3. NPAR(7) fixe le nombre maximum de noeuds autorisé pour décrire un quelconque élément du groupe.
   Un minimum de 3 noeuds et un maximum de 8 noeuds sont utilisés pour décrire ces éléments.
4. Un ordre d'intégration de 2 est suffisant.
   Si l'élément est fortement déformé, un ordre d'intégration plus élevé doit être utilisé.
   

   NTYPEL	type d'élément	nombre max. des points d'intégr.(NINTMX)
   		NINT
   		1	2	3	4
   1	quadrangle	1	4	9	16
   2	triangle	1	3	4	7
   0	les deux	1	4	9	16

5. Pour définir un champ incident il y a deux possibilités:
   1. soit une onde plane:
      • harmonique
      • ondulation de Ricker
      • quelconque
   2. soit un champ quelconque
   1. Pour les ondes planes il faut définir les matériaux, les caractéristiques des ondes et les coordonnées du point de contrôle. Le point de contrôle est un point situé sur la surface libre (cf figure ci-dessous).
      Les deux premiers cas correspondent aux courbes analytiques et l'amplitude des déplacements est donnée au point de contrôle.
      Pour les ondes planes quelconques il faut donner les accélérations horizontales et verticales au point de contrôle sous formes des courbes de chargement définies dans la section 7.6.4.2.3.
      Dans les deux cas le programme calcule les champs de contraintes et de vitesses nécessaires au niveaux des paraxiaux, en tenant compte des déphasages de l'arrivée d'ondes incidentes directes ou réfléchies.
      
      Il faut rappeler que l'on suppose que le point de contrôle correspondant à la solution de demi espace, n'est pas affecté par la présence des irrégularités géométriques ou du comportement.
      Attention: Pour INCTYP=1, 2, et 3 c'est la composition des ondes $Sv$ et $P$ qui permet d'obtenir le champ attendu au point de contrôle.
   2. Pour le champ quelconque il faut que les champs de vitesses et de contraintes soient a priori calculés et donnés sur un fichier accès direct (PREFIXE.prxi) dont la structure est la suivante :
      NMBNO,NMBDT,T0PRX,DTPRX
      
      NUM(1),RHO(1),Cs(1),Cp(1)
      NUM(2),RHO(2),Cs(2),Cp(2)
      ..........................................
      ..........................................
      XNUM($n$),RHO($n$),Cs($n$),Cp($n$) $n$=NMBNO
      
      Vx(1),Vx(2),......................Vx(NMBDT) 1er noeud
      Vy(1),Vy(2),.......................Vy(NMBDT)
      Vz(1),Vz(2),.......................Vz(NMBDT)
      SIGxx(1),SIGxx(2),........SIGxx(NMBDT)
      SIGyy(1),SIGyy(2),........SIGyy(NMBDT)
      SIGzz(1),SIGzz(2),..........SIGzz(NMBDT)
      SIGxy(1),SIGxy(2),.........SIGxy(NMBDT)
      SIGxz(1),SIGxz(2),.........SIGxz(NMBDT)
      SIGyz(1),SIGyz(2),.........SIGyz(NMBDT)
      
      Vx(1),Vx(2),......................Vy(NMBDT) 2ème noeud
      Vy(1),Vy(2),.......................Vz(NMBDT)
      Vz(1),Vz(2),.......................Vz(NMBDT)
      SIGxx(1),SIGxx(2),........SIGxx(NMBDT)
      SIGyy(1),SIGyy(2),........SIGyy(NMBDT)
      SIGzz(1),SIGzz(2),..........SIGzz(NMBDT)
      SIGxy(1),SIGxy(2),.........SIGxy(NMBDT)
      SIGxz(1),SIGxz(2),.........SIGxz(NMBDT)
      SIGyz(1),SIGyz(2),.........SIGyz(NMBDT)
      ..........................................
      ..........................................
      
      Vx(1),Vx(2),......................Vx(NMBDT) $n$ème noeud
      Vy(1),Vy(2),.......................Vy(NMBDT)
      Vz(1),Vz(2),.......................Vz(NMBDT)
      SIGxx(1),SIGxx(2),........SIGxx(NMBDT)
      SIGyy(1),SIGyy(2),........SIGyy(NMBDT)
      SIGzz(1),SIGzz(2),..........SIGzz(NMBDT)
      SIGxy(1),SIGxy(2),..........SIGxy(NMBDT)
      SIGxz(1),SIGxz(2),..........SIGxz(NMBDT)
      SIGyz(1),SIGyz(2),.........SIGyz(NMBDT)
      avec

      NMBNO: nombre de noeuds paraxiaux

      NMBDT: nombre de pas de temps

      T0PRX: temps initial

      DTPRX: pas de temps

      NUM($i$): numéro du noeud $i$

      RHO($i$) : masse volumique au noeud $i$

      Cs($i$): vitesse d'onde $S$ au noeud $i$

      Cp($i$): vitesse d'onde $P$ au noeud $i$

      Vx($j$): vitesse direction $x$ du noeud $i$ à l'instant $j$

      Vy($j$): vitesse direction $y$ du noeud $i$ à l'instant $j$

      Vz($j$): vitesse verticale du noeud $i$ à l'instant $j$

      SIGxx($j$): contrainte $xx$ du noeud $i$ à l'instant $j$

      SIGyy($j$): contrainte $yy$ du noeud $i$ à l'instant $j$

      SIGzz($j$): contrainte $zz$ du noeud $i$ à l'instant $j$

      SIGxy($j$): contrainte $xy$ du noeud $i$ à l'instant $j$

      SIGxz($j$): contrainte $xz$ du noeud $i$ à l'instant $j$

      SIGyz($j$): contrainte $yz$ du noeud $i$ à l'instant $j$

6. Un seul modèle de comportement possible: milieu élastique (MODEL=1).