1. Structure de données sd_char_meca#
Un concept de type sd_char_meca contient : un ou plusieurs chargements et/ou une ou plusieurs conditions aux limites affectés sur un modèle mécanique.
Ce concept est produit par l’un des trois opérateurs : AFFE_CHAR_MECA, AFFE_CHAR_MECA_F ou AFFE_CHAR_MECA_C.
Un concept sd_char_meca est toujours associé à un maillage, par l’intermédiaire des entités NOEUD, GROUP_NO, MAILLE, GROUP_MA sur lesquelles sont définis les chargements.
Sur ce maillage doit avoir été affecté un modèle mécanique.
Un concept sd_char_meca peut contenir un sd_ligrel que l’on appellera “ligrel de charge” et des sd_carte.
1.1. Arborescences de la sd_char_meca#
sd_char_meca (K8) :: = record
◊'.CHME' : sd_char_chme
◊”.DUAL” : sd_char_dual
◊”.TYPE” : S E K8
◊”.LISMA01” : S V I
◊”.LISMA02” : S V I
◊”.TRANS01” : S V R
◊”.TRANS02” : S V R
◊”.POIDS_MAILLE” : S V R
sd_char_chme (K13) :: = record
♦”.MODEL.NOMO” : S E K8 long=1
◊”.VEASS” : S E K8 long=1
◊”.EVOL.CHAR” : S E K8 long=1
◊”.VEISS” : S E K24 long=6
◊”.LIGRE” : sd_ligrel
◊”.CIMPO” : sd_carte
◊”.CMULT” : sd_carte
◊”.DPGEN” : sd_carte
◊”.EPSIN” : sd_carte
◊”.F1D2D” : sd_carte
◊”.F1D3D” : sd_carte
◊”.F2D3D” : sd_carte
◊”.FCO2D” : sd_carte
◊”.FCO3D” : sd_carte
◊”.FELEC” : sd_carte
◊”.FL101” : sd_carte
◊”.FL102” : sd_carte
◊”.FORNO” : sd_carte
◊”.IMPE” : sd_carte
◊”.PESAN” : sd_carte
◊”.PRESS” : sd_carte
◊”.ROTAT” : sd_carte
◊”.SIGIN” : sd_carte
◊”.SIINT” : sd_carte
◊”.VNOR” : sd_carte
◊♦”.ONDPL” : sd_carte
♦”.ONDPR” : sd_carte
sd_char_dual (K13) :: = record
♦”.PRDK” : S V K8 lonmax >= nbpaq
♦”.PRDI” : S V I lonmax >= 3*nbpaq
♦”.PRDSO” : S V K8 lonmax >= 4*nbpaq
1.2. Contenu des objets JEVEUX#
“.CHME.MODEL.NOMO”: nom du modèle associé à la charge (K8)
“.CHME.TEMPE.TEMP”: nom du champ de température dans le cas d’un chargement thermique (K8)
“.CHME.VEASS”: nom du vecteur assemblé dont on impose les composantes comme force au second membre (K8)
“.CHME.VEISS”: S E K24 (dim = 6)
L’objet contient 6 chaînes de caractères K(6) définissant les données nécessaires au chargement FORCE_SOL:
K(1): numéro d’unité de l’évolution temporelle de la contribution en rigidité de l’impédance de sol.
K(2): numéro d’unité de l’évolution temporelle de la contribution en masse de l’impédance de sol.
K(3): numéro d’unité de l’évolution temporelle de la contribution en amortissement de l’impédance de sol.
K(4): numéro d’unité de l’évolution temporelle des forces sismiques imposées à l’interface sol-structure.
K(5): lieu d’interface sol-structure défini par GROUP_NO_INTERF.
K(6): lieu d’interface sol-structure défini par SUPER_MAILLE.
“.TYPE”: type de la charge (K8)
Il contient l’une des 2 chaines de caractères :
“MECA_RE” –> réel pour opérateur AFFE_CHAR_MECA
“MECA_FO” –> fonction pour opérateur AFFE_CHAR_MECA_F
“.LISMA01” : S V I (dim = 2* NBMAIL1)
NBMAIL1 = nombre de mailles introduites derrière les mots-clé TOUTou MAILLEoudans le GROUP_MA mot-clé facteur INTE_ELEC.
Cet objet contient la liste des nœuds des éléments linéiques définissant le conducteur principal
“.LISMA02” : S V I (dim = 2* NBMAIL2)
Même chose pour le conducteur secondaire (s’il existe)
“.TRANS01” : S V I (dim = 6)
/ cas où le mot-clé TRANS de FORCE_ELEC est présent :
tx, ty, tz, 0., 0., 0.
(tx, ty, tz) sont les composantes d’une translation du conducteur principal au
conducteur secondaire
/ cas où le mot-clé SYME de FORCE_ELEC est présent :
x0, y0, z0 , nx, ny, nz
(x0, y0, z0) : coordonnées d’un point
(nx, ny, nz) : composantes de la normale commune au conducteur principal et au
conducteur secondaire
1.2.1. sd_char_dual#
Cette structure de donnée est utilisée pour les relations cinématiques dualisées non-linéaires (à la date du 29/9/2015, cela ne concerne que les relations dues au mot clé LIAISON_SOLIDE).
Le problème est de stocker suffisamment d’informations pour pouvoir réactualiser les relations au cours du calcul non linéaire (si TYPE_CHARGE=”SUIV”). Cette structure de données est amenée à évoluer au fur et à mesure que les possibilités de relations non-linéaires seront développées dans le code.
On définit la notion de «paquet» de relations cinématiques. Un paquet correspond à la liste des relations que l’on affecte dans une charge lors d’un appel à la routine aflrch.F90.
Un paquet peut être linéaire (“LIN”) ou non (“NLIN”, “3D3”, …).
Un paquet est constitué de nbrela relations.
Une relation (quand elle est linéarisée) est constituée de nbterm (à gauche du signe «=») et d’un second membre.
Les nbterm coefficients de la relation sont stockés dans la carte .CHME.CMULT.
Le second membre est stocké dans la carte .CHME.CIMPO.
Il est important de savoir que chaque terme d’une relation cinématique dualisée est associé à un élément fini de Dirichlet «tardif» associé à une maille de type SEG3 «tardive» stockée dans l’objet .NEMA du ligrel “.CHME.LIGRE”. La numérotation de ces mailles tardives est «naturelle»: les relations d’un paquet sont stockées bout à bout et les termes d’une relation sont stockés les uns après les autres. C’est ce qui explique le contenu de l’objet .PRDI: on y stocke des informations permettant de se décaler de paquet en paquet.
Objet .PRDK:
Cet objet est de longueur >= nbpaqu si nbpaqu est le nombre de «paquets» de relations contenus dans la charge. L’attribut LONUTI de cet objet est égal à nbpaqu.
Pour k=1,nbpaqu:
.PRDK(k) = code
Avec code=
“LIN” |
Les relations sont toujours linéaires |
“NLIN” |
Les relations sont a priori non linéaires (mais on n’en sait pas plus). Si ces relations sont ensuite utilisées avec TYPE_CHARGE=”SUIV”, le code doit s’arrêter en erreur. |
“3D3” |
Les relations correspondent à une occurrence du mot clé facteur LIAISON_SOLIDE. Les nœuds «solidifiés» sont dans un espace de dimension 3 et ils forment un nuage de points «volumique» (6 mouvements de corps solide). Remarque importante: Les nœuds solidifiés portent des ddls de translation (DX, DYet DZ) mais ils ne portent pas de ddls de rotation (DRX, DRYet DRZ). |
“3D2” |
Idem “3D3” sauf que le nuage des points est «plan» (6 mouvements de corps solide) |
“3D1” |
Idem “3D3” sauf que le nuage des points est «droit» (5 mouvements de corps solide) |
“2D2” |
Les relations correspondent à une occurrence du mot clé facteur LIAISON_SOLIDE. Les nœuds «solidifiés» sont dans un espace de dimension 2 et ils forment un nuage de points «surfacique» (3 mouvements de corps solide). Remarque importante: Les nœuds solidifiés portent des ddls de translation (DXet DY) mais pas de ddl de rotation (DRZ). |
“2D1” |
Idem “2D2” sauf que le nuage des points est «droit» (3 mouvements de corps solide) |
“ROTA3D” |
Les relations correspondent à une occurrence du mot clé facteur LIAISON_SOLIDE. Les nœuds «solidifiés» sont dans un espace de dimension 3. Il existe au moins un nœud solidifié ne portant les 3 ddls de rotation (DRX, DRYet DRZ). |
“ROTA2D” |
Les relations correspondent à une occurrence du mot clé facteur LIAISON_SOLIDE. Les nœuds «solidifiés» sont dans un espace de dimension 2. Il existe au moins un nœud solidifié ne portant le ddl de rotation DRZ. |
Objet .PRDI:
Cet objet est de longueur >= 3* nbpaqu
.PRDI(3*(k-1)+1): nbrela
.PRDI(3*(k-1)+2): nbterm
.PRDI(3*(k-1)+3): term1
nbrela est le nombre de relations du paquet k.
nbterm est le nombre total de termes du paquet k.
term1 est le numéro de la maille tardive associée au 1er terme de la 1ère relation du paquet k.
Objet .PRDSO:
Cet objet est de longueur >= 4* nbpaqu
Il est dédié au mot clé LIAISON_SOLIDE.
.PRDSO(4*(k-1)+1): no1
.PRDSO(4*(k-1)+2): no2
.PRDSO(4*(k-1)+3): no3
.PRDSO(4*(k-1)+4): no4
no1, no2, no3 et no4 sont les noms de 4 nœuds (au plus) qui permettent de définir le corps solide.
4 nœuds sont nécessaires par exemple pour le cas “3D3”, mais pour le cas “2D1”, on ne stocke que 2 nœuds.
1.3. Description des cartes d’une CHAR_MECA#
Nom carte |
Nom grandeur |
Description |
.CIMPO |
DDLI_R, _C, _F |
second membre des équations de conditions aux limites cinématiques |
.CMULT |
DDLM_R, _C |
coefficients des équations de conditions aux limites cinématiques |
.DPGEN |
NEUT_R |
efforts généralisés |
.EPSIN |
EPSI_R, _F |
déformation initialePRE_EPSI |
.F1D1D |
FORC_R, _C, _F |
force répartie linéique en 1D |
.F1D2D |
FORC_R, _C, _F |
force répartie linéique en 2D |
.F1D3D |
FORC_R, _C, _F |
force répartie linéique en 3D |
.F2D3D |
FORC_R, _C, _F |
force répartie surfacique en 3D |
.FCO2D |
FORC_R, _C, _F |
force répartie pour les coques“2D” |
.FCO3D |
FORC_R, _C, _F |
force répartie pour les coques“3D” |
.FELEC |
FELECR |
paramètres positionnels des conducteurs |
.FL101 |
LISTMA |
forces de Laplace |
.FL102 |
LISTMA |
forces de Laplace |
.FLUX |
FTHM_R, _F |
flux “THM” |
.FORNO |
FORC_R, _F |
forces nodales |
.IMPE |
IMPE_R, _F |
impédance (acoustique) |
.ONDE |
ONDE_R, _F |
amplitude de pression d’onde incidente (acoustique) |
.PESAN |
PESA_R |
vient de PESANTEUR |
.PRESS |
PRES_R, PRES_F |
pression répartie |
.SIGIN |
SIEF_R |
vient de RELA_CINE_BP |
.SIINT |
SIEF_R |
Contrainte initialePRE_SIGM |
.ROTAT |
ROTA_R |
vient de ROTATION |
.VNOR |
SOUR_R, SOUR_F |
vitesse normale d’une maille (acoustique) |
.ONDPL |
NEUT_K24 |
vient de ONDE_PLANE |
.ONDPR |
NEUT_R |
vient de ONDE_PLANE |
1.4. Exemple#
CH = AFFE_CHAR_MECA ( modele : mo
DDL_IMPO : ( GROUP_NO : (A B) dy : 0. )
FACE_IMPO : ( MAILLE : M266 dnor: 0. )
PRES_REP : ( GROUP_MA : GRMA13 pres: 60.));
produit :
====> IMPR_CO DE LA STRUCTURE DE DONNEE : CH ????????????????
ATTRIBUT : F CONTENU : T BASE : >G<
NOMBRE D'OBJETS (OU COLLECTIONS) TROUVES :27
===============================================================================
IMPRESSION DU CONTENU DES OBJETS TROUVES :
-------------------------------------------------------------------------------
IMPRESSION SEGMENT DE VALEURS >CH .CHME.CIMPO.DESC <
1 - 24 6 6 -3 1
6 - -3 2 -3 3 -3
11 - 4 -3 5 -3 6
16 - 2 2 2 2 2
21 - 2
-------------------------------------------------------------------------------
IMPRESSION DE LA COLLECTION : CH .CHME.CIMPO.LIMA
IMPRESSION OBJET DE COLLECTION >CH .CHME.CIMPO.LIMA< OC : 1
1 - -1
IMPRESSION OBJET DE COLLECTION >CH .CHME.CIMPO.LIMA< OC : 2
1 - -2
IMPRESSION OBJET DE COLLECTION >CH .CHME.CIMPO.LIMA< OC : 3
1 - -3
IMPRESSION OBJET DE COLLECTION >CH .CHME.CIMPO.LIMA< OC : 4
1 - -4
IMPRESSION OBJET DE COLLECTION >CH .CHME.CIMPO.LIMA< OC : 5
1 - -5
IMPRESSION OBJET DE COLLECTION >CH .CHME.CIMPO.LIMA< OC : 6
1 - -6
-------------------------------------------------------------------------------
IMPRESSION SEGMENT DE VALEURS >CH .CHME.CIMPO.NOLI <
1 - >CH .CHME.LIGRE.LIEL<>CH .CHME.LIGRE.LIEL<
3 - >CH .CHME.LIGRE.LIEL<>CH .CHME.LIGRE.LIEL<
5 - >CH .CHME.LIGRE.LIEL<>CH .CHME.LIGRE.LIEL<
-------------------------------------------------------------------------------
IMPRESSION SEGMENT DE VALEURS >CH .CHME.CIMPO.NOMA <
1 - >M <
-------------------------------------------------------------------------------
IMPRESSION SEGMENT DE VALEURS >CH .CHME.CIMPO.VALE <
1 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
6 - 0.00000E+00
-------------------------------------------------------------------------------
IMPRESSION SEGMENT DE VALEURS >CH .CHME.CMULT.DESC <
1 - 26 6 6 -3 1
6 - -3 2 -3 3 -3
11 - 4 -3 5 -3 6
16 - 2 2 2 2 2
21 - 2
-------------------------------------------------------------------------------
IMPRESSION DE LA COLLECTION : CH .CHME.CMULT.LIMA
IMPRESSION OBJET DE COLLECTION >CH .CHME.CMULT.LIMA< OC : 1
1 - -1
IMPRESSION OBJET DE COLLECTION >CH .CHME.CMULT.LIMA< OC : 2
1 - -2
IMPRESSION OBJET DE COLLECTION >CH .CHME.CMULT.LIMA< OC : 3
1 - -3
IMPRESSION OBJET DE COLLECTION >CH .CHME.CMULT.LIMA< OC : 4
1 - -4
IMPRESSION OBJET DE COLLECTION >CH .CHME.CMULT.LIMA< OC : 5
1 - -5
IMPRESSION OBJET DE COLLECTION >CH .CHME.CMULT.LIMA< OC : 6
1 - -6
-------------------------------------------------------------------------------
IMPRESSION SEGMENT DE VALEURS >CH .CHME.CMULT.NOLI <
1 - >CH .CHME.LIGRE.LIEL<>CH .CHME.LIGRE.LIEL<
3 - >CH .CHME.LIGRE.LIEL<>CH .CHME.LIGRE.LIEL<
5 - >CH .CHME.LIGRE.LIEL<>CH .CHME.LIGRE.LIEL<
-------------------------------------------------------------------------------
IMPRESSION SEGMENT DE VALEURS >CH .CHME.CMULT.NOMA <
1 - >M <
-------------------------------------------------------------------------------
IMPRESSION SEGMENT DE VALEURS >CH .CHME.CMULT.VALE <
1 - 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
6 - 0.00000E+00 1.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
11 - 0.00000E+00 0.00000E+00 -7.07107E-01 0.00000E+00 0.00000E+00
16 - 0.00000E+00 0.00000E+00 0.00000E+00 7.07107E-01 0.00000E+00
21 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 -7.07107E-01
26 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
31 - 7.07107E-01 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
36 - 0.00000E+00
-------------------------------------------------------------------------------
IMPRESSION SEGMENT DE VALEURS >CH .CHME.DPGEN.DESC <
1 - 58 1 1 1 9999
6 - 8190
-------------------------------------------------------------------------------
IMPRESSION DE LA COLLECTION : CH .CHME.DPGEN.LIMA
IMPRESSION OBJET DE COLLECTION >CH .CHME.DPGEN.LIMA< OC : 1
1 - 0
-------------------------------------------------------------------------------
IMPRESSION SEGMENT DE VALEURS >CH .CHME.DPGEN.NOLI <
1 - > <
-------------------------------------------------------------------------------
IMPRESSION SEGMENT DE VALEURS >CH .CHME.DPGEN.NOMA <
1 - >M <
-------------------------------------------------------------------------------
IMPRESSION SEGMENT DE VALEURS >CH .CHME.DPGEN.VALE <
1 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
6 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 1.00000E+00
11 - 1.00000E+00 1.00000E+00
-------------------------------------------------------------------------------
IMPRESSION DE LA COLLECTION : CH .CHME.LIGRE.LIEL
IMPRESSION OBJET DE COLLECTION >CH .CHME.LIGRE.LIEL< OC : 1
1 - -1 -2 -4 -6 39
IMPRESSION OBJET DE COLLECTION >CH .CHME.LIGRE.LIEL< OC : 2
1 - -3 -5 38
-------------------------------------------------------------------------------
IMPRESSION SEGMENT DE VALEURS >CH .CHME.LIGRE.NBNO <
1 - 8
-------------------------------------------------------------------------------
IMPRESSION DE LA COLLECTION : CH .CHME.LIGRE.NEMA
IMPRESSION OBJET DE COLLECTION >CH .CHME.LIGRE.NEMA< OC : 1
1 - 1 -1 -2 4
IMPRESSION OBJET DE COLLECTION >CH .CHME.LIGRE.NEMA< OC : 2
1 - 119 -3 -4 4
IMPRESSION OBJET DE COLLECTION >CH .CHME.LIGRE.NEMA< OC : 3
1 - 57 -5 -6 4
IMPRESSION OBJET DE COLLECTION >CH .CHME.LIGRE.NEMA< OC : 4
1 - 57 -5 -6 4
IMPRESSION OBJET DE COLLECTION >CH .CHME.LIGRE.NEMA< OC : 5
1 - 41 -7 -8 4
IMPRESSION OBJET DE COLLECTION >CH .CHME.LIGRE.NEMA< OC : 6
1 - 41 -7 -8 4
-------------------------------------------------------------------------------
IMPRESSION SEGMENT DE VALEURS >CH .CHME.LIGRE.NOMA <
1 - >M <
-------------------------------------------------------------------------------
IMPRESSION SEGMENT DE VALEURS >CH .CHME.LIGRE.PRNS <
1 - 4096 4096 4096 4096 4096
6 - 4096 4096 4096
-------------------------------------------------------------------------------
IMPRESSION SEGMENT DE VALEURS >CH .CHME.MODEL.NOMO <
1 - >MO <
-------------------------------------------------------------------------------
IMPRESSION SEGMENT DE VALEURS >CH .CHME.PRESS.DESC <
1 - 63 2 2 1 9999
6 - 2 3 6 6
-------------------------------------------------------------------------------
IMPRESSION DE LA COLLECTION : CH .CHME.PRESS.LIMA
IMPRESSION OBJET DE COLLECTION >CH .CHME.PRESS.LIMA< OC : 1
1 - 0
IMPRESSION OBJET DE COLLECTION >CH .CHME.PRESS.LIMA< OC : 2
1 - 0
-------------------------------------------------------------------------------
IMPRESSION SEGMENT DE VALEURS >CH .CHME.PRESS.NOLI <
1 - > <> <
-------------------------------------------------------------------------------
IMPRESSION SEGMENT DE VALEURS >CH .CHME.PRESS.NOMA <
1 - >M <
-------------------------------------------------------------------------------
IMPRESSION SEGMENT DE VALEURS >CH .CHME.PRESS.VALE <
1 - 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00 0.00000E+00
6 - 0.00000E+00 6.00000E+01 0.00000E+00 0.00000E+00 0.00000E+00
11 - 0.00000E+00 0.00000E+00
-------------------------------------------------------------------------------
IMPRESSION SEGMENT DE VALEURS >CH .TYPE <
1 - >MECA_RE <
====> FIN IMPR_CO DE DE STRUCTURE DE DONNEE : CH ????????????????