This is an old revision of the document!

Coques

Matériaux

Le comportement d'une coque se définit de la même façon que les éléments de volume. Dans un premier temps, il faut donc définir 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 (p.expl. `LinearShellMaterial`)
`param`	nom du paramètre du matériau coque (p.expl. `ELASTIC_MODULUS`)
`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)

LinearShellMaterial

Description

Loi élastique simple linéaire (petites déformations/déplacements).

Paramètres

Nom	Code Metafor	Dépendance
Densité	`MASS_DENSITY`
Module de Young	`ELASTIC_MODULUS`
Coefficient de Poisson	`POISSON_RATIO`

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 `LinearShellFirstDegreeElement`)
`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)

LinearShellFirstDegreeElement

Elément de coque de type linéaire (petites deformations/petits déplacements). Ce sont des quadrangles linéaires à quatre noeuds. Par défaut, les contraintes sont intégrées sur dim*2 point d'intégration. Une épaisseur doit être définie, voir paramètres.

LinearShellSecondDegreeElement

Même élément que LinearShellFirstDegreeElement, mais basé sur un maillage du second degré 2D Meshers (Surfaces). La partie membrane peut être intégrée sur dim*3 points d'intégration en considérant des modes EAS, voir paramètres. La partie bending peut être intégrée sur dim*3 points d'intégration en considérant des modes EAS, voir paramètres.

LinearShellNineNodeSecondDegreeElement

Même élément que LinearShellSecondDegreeElement, mais basé sur un maillage du second degré à 9 noeuds 2D Meshers (Surfaces). La partie membrane peut être intégrée sur dim*3 points d'intégration en considérant des modes EAS, voir paramètres. La partie bending peut être intégrée sur dim*3 points d'intégration en considérant des modes EAS, voir paramètres.

LinearShellSixteenNodeThirdDegreeElement

Même élément que LinearShellNineNodeSecondDegreeElement, mais basé sur un maillage du troisième degré à 16 noeuds 2D Meshers (Surfaces). La partie membrane est tjs intégrée sur dim*4 points d'intégration sans modes EAS (pas de locking), voir paramètres. La partie bending est tjs intégrée sur dim*4 points d'intégration sans modes EAS (pas de locking), voir paramètres.

Paramètres des éléments de 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	-
`GRAVITY_X` / `GRAVITY_Y` / `GRAVITY_Z`	Accélération de gravité	temps
`THICKNESS`	épaisseur initiale	-
`MEMBRANE_NPG`	Nombre de points d'intégration selon une direction pour la résolution de la membrane = 2 résolution classique (pour second degré) = 3 résolution EAS avec 11 ou 13 modes ajoutés (pour second degré à 9 ou 8 noeuds) = 4 résolution sans EAS (pour troisième degré à 16 noeuds)	-
`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 de coque se génèrent de la même manière que les éléments de volume, i.e. à l'aide d'une interaction de type FieldApplicator. Il faudra alors les compléter par des éléments d'interface de type DG pour garantir leur consistance.

app = FieldApplicator(no)
app.push(gObject1)
app.push(gObject2)
...
app.addProperty(prp)        # association d'un ElementProperties
interactionset.add(app)             # ajout dans l'InteractionSet

ou

app = interactionset.add(FieldApplicator(no)) #add renvoie une référence
app.push(gObject1)
app.push(gObject2)
...
app.addProperty(prp)        # association d'un ElementProperties

où

`no`	numéro de l' `Interaction` (c-à-d du `FieldApplicator`)
`gObject1`	l'entité géométrique maillée sur laquelle on applique des éléments finis
`prp`	Propriétés des éléments de coque à générer.

Table of Contents