Interface de type DG pour coques

Une fois les interactions de coques définies, il faut leur associer des interactions d'interface de type DG (Discontinuous Galerkin), en suivant le schéma traditionnel matériau-élément-interaction.

Matériaux

On commence par définit un matériau selon

materialno = materialset.define (number, type)
materialno = materialset(number)
materialno.put(param, valeur)
materialno.depend(param, fct, Key(Lock)))
...

`number`	numéro du matériau (doit être unique, parmi tous les éléments dans `materialset` et pas seulement parmi les coques)
`type`	type de coque (par exemple `DgShellMaterial`)
`param`	nom du paramètre du matériau coque (par exemple `STABILIZATION_PARAMETER`)
`valeur`	valeur du paramètre `param` correspondant
`fct`	numéro de la fonction dont la propriété dépend (facultatif: si aucune dépendance pas de `fct`)
`Lock`	Lock définissant la variable de dépendance de la propriété (obligatoire si on a une dépendance)

DgShellMaterial

Description

Loi DG pour les éléments d'interface entre coques.

Paramètres

Nom	Code Metafor	Dépendance
Paramètre de stabilisation (>1)	`STABILIZATION_PARAMETER`

Eléments

Une fois le matériau défini, il faut l'associer à un ElementProperties, selon

prp = ElementProperties(typeEl)
prp.put(param1, valeur1)
prp.depend(param1, fct1, Lock1)) #facultatif
...

où

`typeEl`	type d'élément souhaité (par exemple `DgShellFirstDegreeElement`)
`param1`	nom de la propriété relative à l'élément (p.expl. `MATERIAL` pour associer le matériau souhaité)
`valeur1`	“valeur” de la propriété
`fct1`	numéro de la fonction dont la propriété dépend (facultatif: si aucune dépendance pas de `fct`)
`Lock1`	Lock définissant la variable de dépendance de la propriété (obligatoire si on a une dépendance)

DgShellFirstDegreeElement

Elément d'interface de type DG à insérer entre des coques. Ce sont des lignes linéaires à deux noeuds, mais ils interfèrent sur tous les ddl des éléments adjacents. Par défaut, les contraintes sont intégrées sur dim*2 points d'intégration. Le matériau utilisé est dgShellMaterial.

DgShellSecondDegreeElement

Même élément que DgShellFirstDegreeElement, mais au second degré pour éléments de coque à 8 noeuds.

DgShellNineNodeSecondDegreeElement

Même élément que DgShellSecondDegreeElement, mais au second degré pour éléments de coque à 9 noeuds.

DgShellSixteenNodeThirdDegreeElement

Même élément que DgShellNineNodeSecondDegreeElement, mais au troisième degré pour éléments de coque à 16 noeuds.

Paramètres des éléments d'interface DG pour coques

Code Metafor	Description	Type de dépendance
`MATERIAL`	Numéro du matériau à considérer	-
`STIFFMETHOD`	Méthode pour calculer la matrice de raideur = `STIFF_ANALYTIC` : matrice analytique (défaut) = `STIFF_NUMERIC` : matrice numérique	-
`BENDING_NPG`	Nombre de points d'intégration selon une direction pour la résolution de la flexion = 2 résolution classique (pour second degré) = 3 résolution EAS avec 22 ou 7 modes ajoutés(pour second degré à 9 ou 8 noeuds) = 4 résolution sans EAS (pour troisième degré à 16 noeuds)	-

Interactions

Les éléments d'interface DG pour coques sont générés en début de calcul par une Interaction DgShellInteraction définie sur les entités voulues.

dg = DgShellInteraction(numero)
dg.push(gObject1)
dg.push(gObject2)
...
dg.addProperty(prp)
interactionset.add(dg)                 # ajout dans l'InteractionSet

ou

dg = interactionset.add( DgShellInteraction(numero) )
dg.push(gObject1)
dg.push(gObject2)
...
dg.addProperty(prp)

`numero`	Numéro utilisateur de l'interaction
`gObject1`, `gObject2`	Entité géométrique maillées
`prp`	Propriétés des éléments d'interface DG à générer

Attention pour le moment, si deux entités ont des noeuds en commun, il y aura trop d'éléments générés. C'est à corriger…