This is an old revision of the document!
Table of Contents
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…