Differences

This shows you the differences between two versions of the page.

--- doc:user:elements:boundaries:pressure_shear [2013/07/12 19:27] – external edit 127.0.0.1
+++ doc:user:elements:boundaries:pressure_shear [2021/09/22 11:43] – boman
@@ Line 1: / Line 1: @@
-====== Pression et cisaillement ======
+====== Pressure and shear ======
-===== Matériaux =====
+Pressure and shear boundary conditions can be applied using special finite elements defined on the boundary of the mesh. The geometry of these elements are lines in 2D corresponding to the boundary edges of the elements of the 2D mesh.
+In 3D, pressure elements are triangles or quads, depending of the shape of the facets of the elements of the mesh near the boundary.
-Comme les interactions de pression/cisaillement sont des interactions de conditions limites (''LoadingInteraction''), il ne faut pas y associer de matériau.
+As any finite element, their definition requires an ''ElementProperty'' object and an ''Interaction'' (''LoadingInteraction'').
+Unlike contact finite element, there is no "material" associated to this kind of element.
-===== Elément =====
+The ''ElementProperty'' contains the type of element and the values for the pressure and the shear, as well as functions defining the time evolution of these values.
-On commence donc par définir un ''[[doc:user:elements:general:def_element_properties|ElementProperties]]'' selon
+These elements should not be mixed with [[doc:user:conditions:dead_loads|dead loads]]. In most cases, pressure elements should be chosen over dead loads because:
+    * Pressure and shear are always defined in local axes, so the direction of the resulting force depends on the orientation of the boundary edge/face on which it is applied.
+    * Pressure is distributed properly over all edges/faces on which the ''LoadingInteraction'' is applied. Consequently, the total applied force does not depend on the mesh.
+{{ doc:user:pressure.gif |Pressure Elements}}
+===== Element =====
+Therefore, the first step consists in defining an ''[[doc:user:elements:general:def_element_properties|ElementProperties]]'', as
   prp = ElementProperties(typeEl)
-  prp.put(param1, valeur1)
+  prp.put(param, value)
-  prp.depend(param1, fct1, Lock1)) #facultatif
+  prp.depend(param, fct, lock)) # optional
   ...
-où
+where
+| ''typeEl'' | desired element (for example ''Pressure[2|3]DElement'' - see below)  |
+| ''param'' | name of the property associated to the element (e.g. ''PRESSURE'') |
+| ''value'' | value of the corresponding property (e.g. value for the pressure  |
+| ''fct''   | function which characterizes the dependency of the property (optional: no fct if the value is constant) |
+| ''lock''  | [[doc:user:general:locks|lock]] which defines the dependency variable of the property (e.g. time-dependence = ''TM'')  |
-|''typeEl'' | type d'élément souhaité (par exemple ''Traction[2|3]DElement'')|
-|''param1''| nom de la propriété relative à l'élément (p.expl. ''PRESSURE'')|
-|''valeur1''| "valeur" de la propriété |
-|''fct1''| numéro de la [[doc:user:general:fonctions|fonction]] dont la propriété dépend (facultatif: si aucune dépendance pas de ''fct'') |
-|''Lock1''| [[doc:user:general:locks|Lock]] définissant la variable de dépendance de la propriété (obligatoire si on a une dépendance) |
-==== Traction[2|3]DElement ====
+==== Pressure2DElement, Pressure3DElement ====
-Elément de traction/cisaillement pour un maillage quadrangulaire (2D) ou hexahédrique (3D).
+Pressure/shear element for a mesh made of triangles/quadrangles in 2D or hexahedra in 3D.
-=== Paramètres ===
+=== Parameters ===
+^   Name                                          ^     Description       ^  Dependency ^
+| ''STIFFMETHOD'' | Method used to compute the stiffness matrix\\ = ''STIFF_ANALYTIC'' : analytic matrix (default)\\ = ''STIFF_NUMERIC''  : numerical matrix  |  -  |
+| ''PRESSURE'' | Pressure on the entity on which the propriety is applied |  time  |
+| ''SHEAR_KSI'' / ''SHEAR_ETA''  | Shear along $\xi$ or $\eta$ on the entity on which the propriety is applied |  time  |
+| ''NPG''    | Number of integration points (default : 2)    |  -  |
-^  Code Metafor  ^  Description  ^  Type de dépendance  ^
+<note tip>Pressure2DElement / Pressure2DElement were named Traction2DElement / Traction3DElement before v3484. The sign of the pressure has also been changed in 2D.  </note>
-| ''STIFFMETHOD'' | Méthode pour calculer la matrice de raideur\\ = ''STIFF_ANALYTIC'' : matrice analytique (défaut)\\ = ''STIFF_NUMERIC''  : matrice numérique           |  -  |
-| ''PRESSURE'' | Pression sur l'entité à laquelle s'applique la propriété |  temps  |
-| ''SHEAR_KSI'' / ''SHEAR_ETA''  | Cisaillement selon $\xi$ ou $\eta$ sur l'entité à laquelle s'applique la propriété |  temps  |
-| ''NPG''    | Nombre de points d'intégration (défaut : 2)    |  -  |
-==== TriangleTraction3DElement ====
+==== TrianglePressure3DElement ====
-Elément de traction/cisaillement pour un maillage tétrahédrique.
+Pressure/shear element for a mesh made of tetrahedra
+=== Parameters ===
+^   Name                                          ^     Description       ^  Dependency ^
+| ''STIFFMETHOD'' | Method used to compute the stiffness matrix\\= ''STIFF_ANALYTIC'' : analytic matrix (default)\\ = ''STIFF_NUMERIC''  : numerical matrix \\ Not applicable, only the numerical one exists ! |  -  |
+| ''PRESSURE'' | Pressure on the entity on which the propriety is applied |  time  |
+| ''SHEAR_KSI'' / ''SHEAR_ETA''  | Shear along $\xi$ or $\eta$ on the entity on which the propriety is applied |  time  |
+| ''NPG''    | Number of integration points (default : 1)    |  -  |
-=== Paramètres ===
-^  Code Metafor  ^  Description  ^  Type de dépendance  ^
+==== ContactPressure2DElement, ContactPressure3DElement ====
-| ''STIFFMETHOD'' | Méthode pour calculer la matrice de raideur\\ = ''STIFF_ANALYTIC'' : matrice analytique (défaut)\\ = ''STIFF_NUMERIC''  : matrice numérique           |  -  |
-| ''PRESSURE'' | Pression sur l'entité à laquelle s'applique la propriété |  temps  |
+Derivations of ''Pressure2DElement'' and ''Pressure3DElement'' with the possibility to set the pressure/shear force of a node to zero, when this node is in contact within a specified contact interaction.
-| ''SHEAR_KSI'' / ''SHEAR_ETA''  | Cisaillement selon $\xi$ ou $\eta$ sur l'entité à laquelle s'applique la propriété |  temps  |
+To determine the contact status, the contact flag of the previous time step is used.
-| ''NPG''    | Nombre de points d'intégration (défaut : 1)    |  -  |
+=== Parameters ===
+The same than in Traction2DElement and Traction3DElement, with the following additional parameters:
+^  Name  ^  Description  ^  Dependency  ^
+| ''NUM_CONTACT_INTERACTION'' | Number of the contact interaction, which is tested for the contact status |  -  |
+| ''TRACTION_AT_TRANSITION_ON'' | If ''True'' (default), the traction of a node belonging to an element at the transition from contact to no contact is not set equal to zero, although the node is in contact  |  -  |
+{{ :doc:user:elements:boundaries:tractionattransition.png?600 |}}
 ===== Interaction =====
-L'interaction se définit de la manière suivante:
+The interaction used for pressure/shear element is a ''LoadingInteraction''. It is defined as:
   load = LoadingInteraction(no)
@@ Line 60: / Line 87: @@
   interactionset.add(load)
-où
+where
-| ''no''       | Numéro de l'interaction |
-| ''gObject1'', ''gObject2''  | Entités géométriques maillées sur lesquelles sont appliquées les conditions aux limites |
-| ''prp''      | [[doc:user:elements:general:def_element_properties|Propriétés]] des [[#Elément|éléments de conditions aux limites]] à générer |
-==== Remarque ====
-Les tractions, pressions et cisaillement sont gérés par des "éléments de traction". Leur création se fait grâce à une ''[[doc:user:general:glossaire#Interaction]]'' spécifique (''LoadingInteraction'') à laquelle on associe un [[doc:user:elements:general:def_element_properties|ElementProperties]]. Il est important de ne pas confondre ces éléments avec les [[doc:user:conditions:dead_loads|charges mortes]]. Dans la plupart des cas, les éléments de tractions sont préférables aux charges mortes:
-    * la traction (ou pression/cisailement) est toujours définie dans des axes locaux si bien que la direction de la force est fonction de l'orientation de l'arête (2D) ou facette (3D) sur laquelle elle s'applique.
+| ''no''       | number of the ''Interaction''  |
-    * La traction est répartie correctement sur toutes les arêtes/facettes de l'entité sur laquelle s'applique l' ''Interaction''. La force résultante obtenue est donc indépendante du maillage utilisé.
+| ''gObject1'', ''gObject2''  | meshed geometric entity where the boundary conditions are applied (one part of the boundary of the mesh)  |
+| ''prp''      | [[doc:user:elements:general:def_element_properties|Properties]] of [[#Element|boundary condition elements]] to be generated  |
-{{ doc:user:pressure.gif |Elements de pression}}
-:!: l'élément de traction est utilisé "à l'envers" en 3D, si bien qu'il faut appliquer une valeur positive pour imposer une pression !!