Autres chargements ================== Mot-clé LIAISON_INTERF ---------------------- .. code-block:: text LIAISON_INTERF = _F( ◆ MACR_ELEM_DYNA = macr_elem_dyna, ◇ TYPE_LIAISON = / "RIGIDE" (par défaut), / "SOUPLE", ) Le mot-clé facteur LIAISON_INTERF permet dedéfinir des relations linéaires entre les degrés de liberté physiques des interfaces de la partie de modèle en éléments finis et les coordonnées généralisées de modes de représentation réduite des mouvements d’interface contenus dans certains macro-éléments de condensation statique. Il est utilisable avec un modèle contenant à la fois des éléments finis et des macro-éléments statiques condensant certains sous-domaines. .. code-block:: text ◆ MACR_ELEM_DYNA Nom du macr_elem_dyna qui sert à définir les relations linéaires entre les degrés de liberté physiques de l’interface entre le domaine non condensé modélisé en éléments finis et un domaine condensé par le macro-élément et les composantes des nœud assimilés à des coordonnées généralisées de modes de mouvements d’interface. Cela est nécessaire seulement quand les modes de mouvements d’interface sont une base réduite de tous les modes contraints correspondant chacun à un mode de déplacement pour chaque degré de liberté physique de l’interface. On génère ainsi des relations de type LIAISON_DDL dont les coefficients sont calculés de façon transparente pour l’utilisateur entre les nœuds de l’interface dynamique du macro-élément et ceux associés à la base de réduction qui a servi à constituer le macro-élément. .. code-block:: text ◇ TYPE_LIAISON = ’RIGIDE’/’SOUPLE’ * Si ’RIGIDE’, on écrit la relation entre les degrés de liberté physiques de l’interface :math:`{U}_{\Sigma}` et les composantes des nœud assimilés à des coordonnées généralisées :math:`q` de modes de mouvements d’interface :math:`\Phi` sous la forme de produit simple: :math:`{U}_{\Sigma}={\Phi}_{q}` . Ce choix permet d’avoir une liaison plus rigide qu’en prenant en compte tous les modes contraints correspondant chacun à un mode de déplacement pour chaque degré de liberté physique de l’interface. * Si ’SOUPLE’, on écrit la relation entre les degrés de liberté physiques de l’interface :math:`{U}_{\Sigma}` et les composantes des nœud assimilés à des coordonnées généralisées :math:`q` de modes de mouvements d’interface :math:`\Phi` sous la forme de produit double: :math:`{\Phi}^{T}{U}_{\Sigma}={\Phi}^{T}\Phi q` . Ce choix permet d’avoir une liaison plus souple qu’en prenant en compte tous les modes contraints correspondant chacun à un mode de déplacement pour chaque degré de liberté physique de l’interface. Mot-clé RELA_CINE_BP -------------------- .. code-block:: text RELA_CINE_BP = _F( ◆ CABLE_BP = cabl_precont, ◇ RELA_CINE = / "OUI" (par défaut), / "NON", # Si: equal_to("RELA_CINE", 'OUI') ◇ SIGM_BPEL = / "OUI", / "NON" (par défaut), # Si: equal_to("RELA_CINE", 'NON') ◇ SIGM_BPEL = "OUI", ◇ TYPE_EPX = / "ADHE" (par défaut), / "GLIS", / "FROT", ◇ DIST_MIN = float, ) Ce type de chargement peut être défini pour un système mécanique comprenant une structure béton et ses câbles de précontrainte. Les profils initiaux de tension dans les câbles, ainsi que les coefficients des relations cinématiques entre les degrés de liberté des nœuds des câbles et les degrés de liberté des nœuds de la structure béton sont déterminés préalablement par l’opérateur DEFI_CABLE_BP [:external:ref:`U4.42.04 `]. Les concepts cabl_precont produits par cet opérateur apportent toutes les informations nécessaires à la définition du chargement. Les occurrences multiples sont autorisées pour le mot-clé facteur RELA_CINE_BP, afin de permettre dans un même appel à l’opérateur AFFE_CHAR_MECA de définir les contributions de chacun des groupes de câbles ayant fait l’objet d’appels distincts à l’opérateur DEFI_CABLE_BP [:external:ref:`U4.42.04 `]. À chaque groupe de câbles considéré, défini par un concept cabl_precont, est associée une occurrence du mot-clé facteur RELA_CINE_BP. Le chargement ainsi défini sert ensuite à calculer l’état d’équilibre de l’ensemble structure béton / câbles de précontrainte. Cependant, la prise en compte de ce type de chargement n’est pas effective dans tous les opérateurs de résolution. Le chargement de type RELA_CINE_BP n’est reconnu pour l’instant que par l’opérateur STAT_NON_LINE [:external:ref:`U4.51.03 `], comportements incrémentaux exclusivement. .. code-block:: text ◆ CABLE_BP Concept de type cabl_precont produit par l’opérateur DEFI_CABLE_BP [:external:ref:`U4.42.04 `]. Ce concept apporte d’une part la carte des contraintes initiales dans les éléments des câbles d’un même groupe, et d’autre part les listes des relations cinématiques entre les degrés de liberté des nœuds de ces câbles et les degrés de liberté des nœuds de la structure béton. .. code-block:: text ◇ SIGM_BPEL = ’OUI’/’NON’ Indicateur de type texte par lequel on spécifie la prise en compte des contraintes initiales dans les câbles; la valeur par défaut est ’NON’. Dans le cas ’NON’, seul le liaisonnement cinématique est pris en compte. C’est utile si on enchaîne des STAT_NON_LINE alors qu’on a des câbles de précontrainte. Pour le premier STAT_NON_LINE il faut avoir mis ’OUI’, de telle sorte que l’on met en place la tension dans les câbles. En revanche, pour les STAT_NON_LINE suivants, il ne faut considérer comme chargement que les liaisons cinématiques et donc définir le chargement avec SIGM_BPEL = ’NON’, sinon la tension est comptée deux fois. Depuis la restitution de la macro pour mettre en tension les câbles, l’utilisateur ne devrait plus avoir besoin de faire un AFFE_CHAR_MECA avec SIGM_BPEL = ’OUI’, cela devrait ainsi éviter les risques d’erreur. .. code-block:: text ◇ RELA_CINE = ’OUI’/’NON’ Indicateur de type texte par lequel on spécifie la prise en compte des relations cinématiques entre les degrés de liberté des nœuds des câbles et les degrés de liberté des nœuds de la structure béton; valeur par défaut ’OUI’. .. code-block:: text ◇ DIST_MIN = float Voir LIAISON_SOLIDE §:ref:`5.16 `. .. code-block:: text ◇ TYPE_EPX = ’ADHE’/’GLIS’/’FROT’ Ce mot-clé n’a d’effet que dans CALC_EUROPLEXUS. Il permet d’indiquer si on souhaite des liaisons câble-béton totales (c’est-à-dire dans les 3 directions de l’espace, correspondant à ce chargement Aster), des liaisons glissantes ou frottantes. Mot-clé VITE_FACE ----------------- .. code-block:: text VITE_FACE = _F( ◇ GROUP_MA = grma, ◆ / VNOR = float, / DIRECTION = float, # Si: exists("DIRECTION") ◆ VITE = float, ) Le mot-clé facteur VITE_FACE permet d’appliquer des vitesses à une face dans le cas de la modélisation 2D_FLUIDE, AXIS_FLUIDE, 3D_FLUIDE ou éléments d'interface fluide/structure ( \*_FLUI_STRU). Affectation topologique: .. code-block:: text ◇ GROUP_MA Le chargement est affecté sur les mailles qui sont nécessairement des segments en 2D et des triangles ou quadrangles en 3D. Il y a deux manières d’affecter une vitesse: 1. Si on renseigne VNOR, alors la vitesse (de norme VNOR) sera **normale** à la face sur laquelle elle est appliquée (définie par GROUP_MA); 2. Si on renseigne DIRECTION qui donne la direction de la vitesse alors VITE sera la norme de cette vitesse. Cette méthode n’est utilisable que si on a choisi une formulation de type U_PHI (mot-clef dans AFFE_MODELE). La direction et la norme (VNOR ou VITE) peuvent être des réels ou des fonctions. .. _RefHeading__26049_722879224: Mot-clé ONDE_PLANE ------------------ .. code-block:: text ONDE_PLANE =_F( ◇ GROUP_MA = grma, ◆ TYPE_ONDE = / "S", / "P", / "SV", / "SH", ◆ DIRECTION = (kx,ky,kz), ◆ FONC_SIGNAL = f, ◇ DEPL_IMPO = fd, ◇ COOR_SOURCE = (x0, y0, z0), # Si: exists("COOR_SOURCE") ◇ COOR_REFLECHI = (x1, y1, z1), ) Le mot-clé facteur ONDE_PLANE (une seule occurrence est permise) permet d’imposer un chargement sismique par onde plane, correspondant aux chargements classiquement rencontrés lors des calculs d’interaction sol-structure par les équations intégrales (voir [:external:ref:`R4.05.01 `] ). En harmonique, une onde plane élastique est caractérisée par sa direction, sa pulsation et son type (onde :math:`P` pour les ondes de compression, ondes :math:`S` , :math:`\mathit{SV}` ou :math:`\mathit{SH}` pour les ondes de cisaillement). En transitoire, la donnée de la pulsation, correspondant à une onde stationnaire en temps, doit être remplacée par la donnée d’un profil de déplacement dont on va prendre en compte la propagation au cours du temps dans la direction de l’onde. Plus précisément, on caractérise: * Une onde :math:`P` par la fonction :math:`u(x,t)=f\left(k\cdot x-{C}_{p}t\right)` ; * Une onde :math:`S` par la fonction :math:`u(x,t)=f\left(k\cdot x-{C}_{s}t\right)\wedge k` ; Avec: * :math:`k` le vecteur unitaire de direction; * :math:`f` le profil de l’onde donné selon la direction :math:`k` ; **Il ne peut y avoir qu’une seule occurrence de ce mot-clef dans** AFFE_CHAR_MECA *\*.** Affectation topologique: .. code-block:: text ◇ GROUP_MA Le chargement est affecté sur les mailles de frontières absorbantes concernées par l’introduction de l’onde incidente. Si rien n’est donné, par défaut, ce sont toutes les mailles de la modélisation ABSO qui sont concernées. .. code-block:: text ◇ COOR_SOURCE Cette liste de paramètres réels permet de déterminer le déphasage en temps dû au passage de l’onde, en indiquant les coordonnées du point d’entrée :math:`{\mathrm{x}}_{0}` du chargement par onde plane dans la structure. Si ce paramètre n'est pas présent, alors on renseignera l'opérande FONC_SIGNAL avec une nappe de signaux temporels. .. code-block:: text ◇ COOR_REFLECHI Cette liste de paramètres réels permet la prise en compte de l'onde réfléchie, activée si on donne les coordonnéesdu point de sortie :math:`{\mathrm{x}}_{1}`. Dans ce cas, l'expression du profil de l'onde en tenant compte du déphasage spatial lié au passage de l'onde :math:`\dot{f}\left(k.(\mathrm{x}-{\mathrm{x}}_{0})-{C}_{m}t\right)` devient alors :math:`\dot{f}\left(k.(\mathrm{x}-{\mathrm{x}}_{0})-{C}_{m}t\right)+\dot{f}\left(k.(2{\mathrm{x}}_{1}-{\mathrm{x}}_{0}-\mathrm{x})-{C}_{m}t\right)` avec :math:`m=S` ou :math:`P` . .. image:: images/10000201000004FD000006BF7BB4B4564A10FC81.png :width: 3.8693in :height: 5.2335in .. _RefImage_10000201000004FD000006BF7BB4B4564A10FC81.png: Figure 9.4-1: Schéma d'une propagation d'onde plane inclinée dans un modèle 2D Pour clarifier l’utilisation des opérandes COOR_SOURCE et COOR_REFLECHI on fait référence au schéma dans la . Si on considère une onde plane inclinée d’un angle :math:`\alpha` avec la verticale, on identifie de façon évidente le premier point qui rentre en contact avec le domaine d’étude, c’est-à-dire le point :math:`{\mathrm{x}}_{0}` . Ce point est considéré comme le point source de l’onde incidente et ses coordonnées représentent la donnée d’entrée de l’opérande COOR_SOURCE. Avec le même raisonnement, si l’on souhaite rajouter la prise en compte de l’onde réfléchie, on peut clairement identifier le premier point de la surface libre qui génère l’onde réfléchie, c’est-à-dire le point :math:`{\mathrm{x}}_{1}` . Ce point est considéré comme le point source de l’onde réfléchie et ses coordonnées représentent la donnée d’entrée de l’opérande COOR_REFLECHI. Remarques concernant l’utilisation des opérandes COOR_SOURCE et COOR_REFLECHI: * Dans le cas d’une modélisation 2D (D_PLAN, par exemple) on travaille dans le repère XY. Dans ce cas il faudra rentrer tout simplement les coordonnés X et Y du point source de l’onde incidente (COOR_SOURCE) et, éventuellement, de l’onde réfléchie (COOR_REFLECHI). La troisième coordonnée, même si elle est rentrée, ne sera pas considérée; * Dans le cas d’une propagation verticale, le déphasage de l’onde est indépendante des coordonnées X et Y en 3D (X en 2D). En effet, le front d’onde est parfaitement horizontal et il est déphasé exclusivement par rapport à la côte verticale (Z en 3D et Y en 2D). Par conséquent, on peut rentrer n’importe quelle coordonnée X et Y en 3D (X en 2D) pour les points source de l’onde incidente et de l’onde réfléchie. Ce qui compte est de bien rentrer la coordonnée Z en 3D (Y en 2D) pour bien calculer le déphasage lié à la propagation verticale. .. code-block:: text ◆ TYPE_ONDE = ’P’/’S’/’SV’/’SH’ Type de l’onde (compression ou cisaillement). .. csv-table:: ":math:`P`", "onde de compression" ":math:`\mathit{SV}`", "onde de cisaillement (uniquement en 3D)" ":math:`\mathit{SH}`", "onde de cisaillement (uniquement en 3D)" ":math:`S`", "onde de cisaillement (uniquement en 2D)" Les directions des ondes :math:`P`, :math:`\mathit{SV}`et :math:`\mathit{SH}`sont déterminées à partir du vecteur :math:`V` renseigné par DIRECTION. A savoir : * :math:`P` est colinéaire à :math:`V` et normé à 1, * :math:`\mathit{SH}` est l'intersection du plan horizontal et du plan normal à :math:`V`, et normé à 1, * :math:`\mathit{SV}` est le produit vectoriel de :math:`\mathit{SH}`et de :math:`P`. Il existe un cas d'indétermination avec cette règle quand le plan horizontal et le plan normal sont confondus. Dans ce cas, si :math:`V`= :math:`Z` purement vertical, on impose :math:`\mathit{SH}`= :math:`Y`, et :math:`\mathit{SV}`= :math:`X`. .. code-block:: text ◆ DIRECTION Direction de propagation de l’onde. .. code-block:: text ◆ FONC_SIGNAL Dérivée du profil de l’onde :math:`f(t)` pour :math:`t\in [0,+\infty [` . **Attention:** c’est la fonction correspondant à la vitesse :math:`v(t)=\dot{u}(t)` que l’utilisateur donne dans FONC_SIGNAL. Si l'opérande COOR_SOURCE est absent, on donne une nappe de signaux temporels en vitesse paramétrés par des valeurs de coordonnées de propagation :math:`k.\mathrm{x}` pour la variable :math:`X` . .. code-block:: text ◇ DEPL_IMPO Profil de l’onde :math:`\mathit{fd}(t)` pour :math:`t\in [0,+\infty [` . **Attention:** c’est la fonction correspondant à l’intégrale de la vitesse :math:`v(t)=\dot{u}(t)` que l’utilisateur donne **obligatoirement** dans DEPL_IMPO. Cette donnée **facultative** n’est activée que si on a des raideurs ajoutées sur la frontière absorbante, activées par la donnée d’une valeur non nulle de l’opérande LONG_CARA du mot clé ELAS de DEFI_MATERIAU. Si l'opérande COOR_SOURCE est absent, on donne une nappe de signaux temporels en déplacement paramétrés par des valeurs de coordonnées de propagation :math:`k.\mathrm{x}` pour la variable :math:`X` . Remarque concernant l’utilisation de cette charge: Il est fortement conseillé d’orienter au moyen des mots clés ORIE_PEAU_2D ou ORIE_PEAU_3D de MODI_MAILLAGE les mailles de peau affectées par une modélisation de frontière absorbante vers l’extérieur du domaine qu’elles délimitent. Mot-clé ONDE_FLUI ----------------- .. code-block:: text ONDE_FLUI = _F( ◇ GROUP_MA = grma, ◆ PRES = float, ) Le mot-clé facteur ONDE_FLUI permet d’appliquer une amplitude de pression d’onde incidente sinusoïdale arrivant normalement à une face sur les modélisations 3D_FLUIDE, 2D_FLUIDE et AXIS_FLUIDE. Affectation topologique: .. code-block:: text ◇ GROUP_MA L’onde est appliquée sur ces mailles . Ce sont nécessairement des mailles de bord (segments en 2D et triangles/quadrangles en 3D). ◆ PRES Amplitude de pression d’onde incidente sinusoïdale arrivant normalement à la face. Mot-clé FLUX_THM_REP -------------------- .. code-block:: text FLUX_THM_REP = _F( ◆ / TOUT = "OUI" (ou non renseigné), / GROUP_MA = grma, ◆ | FLUN = float, | FLUN_HYDR1 = float, | FLUN_HYDR2 = float, | FLUN_FRAC = fonction / formule / nappe, ) Le mot-clé facteur FLUX_THM_REP permet d’appliquer à un domaine de milieu continu 2D ou 3D un flux de chaleur et/ou un apport de masse fluide (flux hydraulique). Les flux hydrauliques (eau et vapeur) sont définis par: .. math:: :label: eq-20 \begin{array}{c}{\varphi}^{e}={\rho}_{e}\left(\nabla {P}_{e}-{\rho}_{e}g\right)\cdot n\\ {\varphi}^{v}={\rho}_{v}\left(\nabla {P}_{v}-{\rho}_{v}g\right)\cdot n\end{array} Avec la masse volumique de l’eau :math:`{\rho}_{e}` , la masse volumique de la vapeur :math:`{\rho}_{\nu}` , la pression de l’eau :math:`{P}_{e}` (degré de liberté PRE1 ) et la pression de la vapeur :math:`{P}_{v}` (degré de liberté PRE2 ). Le flux de chaleur est défini par: .. math:: :label: eq-21 {\varphi}_{T}={\lambda}_{T}\frac{\partial T}{\partial n}+{h}_{m}^{e}{\varphi}^{e}+{h}_{m}^{v}{\varphi}^{v}+{h}_{m}^{a}{\varphi}^{a} Avec les enthalpies massiques de l’eau :math:`{h}_{m}^{l}` , de la vapeur :math:`{h}_{m}^{v}` , de l’air :math:`{h}_{m}^{a}` et le flux d’air :math:`{\varphi}^{a}` . Affectation topologique: .. code-block:: text ◆ / TOUT = "OUI" (ou non renseigné), / GROUP_MA = grma, Le chargement est affecté sur ces mailles Ce sont nécessairement des mailles de bord (segments en 2D et triangles/quadrangles en 3D).Ce chargement n’est utilisable qu’en (thermo)-hydraulique ou en hydraulique pur (THM,THH,THHM,HMou HHM). .. code-block:: text ◆ | FLUN Valeur du flux de chaleur. .. code-block:: text ◆ | FLUN_HYDR1 Valeur du flux hydraulique associé au constituant eau. .. code-block:: text ◆ | FLUN_HYDR2 Valeur du flux hydraulique associé au constituant vapeur. .. code-block:: text ◆ | FLUN_FRAC Valeur du flux hydraulique injecté dans une interface cohésive pour les modèles HM-XFEM. Ce chargement est imposé exclusivement à l'aide d'une fonction de l'espace sur les éléments principaux et non pas les éléments de bord. Pour plus de détails concernant l'utilisation de ce mot-clef, se reporter aux documentations [:external:ref:`R7.02.18 `] et [:external:ref:`U2.05.02 `]. Mot-clé ECHANGE_THM ------------------- .. code-block:: text ECHANGE_THM = _F( ◆ / TOUT = "OUI" (ou non renseigné), / GROUP_MA = grma, ◆ COEF_11 = float, ◇ COEF_12 = float, ◇ COEF_21 = float, ◇ COEF_22 = float, ◆ PRE1_EXT = float, ◇ PRE2_EXT = float, ) Le mot-clé facteur ECHANGE_THM permet pour une modélisation non saturée d’appliquer à un domaine de milieu continu 2D ou 3D une condition d'échange sur les pressions hydrauliques extérieures. Les flux hydrauliques (eau et air) qui en résultent sont définis par: .. math:: :label: eq-22 \begin{array}{c}\left({M}_{w}+{M}_{\text{vp}}\right).n=C11({p}_{1}-{p}_{{}_{1}}^{\text{ext}})-C12({p}_{2}-{p}_{{}_{2}}^{\text{ext}})\\ \left({M}_{\text{ad}}+{M}_{\text{as}}\right).n=C21({p}_{1}-{p}_{{}_{1}}^{\text{ext}})-C22({p}_{2}-{p}_{{}_{2}}^{\text{ext}})\end{array} Avec la pression capillaire :math:`{p}_{1}` (degré de liberté PRE1 ) et la pression gaz :math:`{p}_{2}` (degré de liberté PRE2 ). :math:`{M}_{w}+{M}_{\text{vp}}` correspond au flux d'eau (somme du flux d'eau liquide :math:`{M}_{w}` et du flux de vapeur :math:`{M}_{\text{vp}}` ) et :math:`{M}_{\text{ad}}+{M}_{\text{as}}` correspond au flux d'air (somme du flux d'air dissous :math:`{M}_{\mathit{ad}}` et du flux d'air sec :math:`{M}_{\text{as}}` ). Affectation topologique: .. code-block:: text ◆ / TOUT = "OUI" (ou non renseigné), / GROUP_MA = grma, Le chargement est affecté sur ces mailles Ce sont nécessairement des mailles de bord (segments en 2D et triangles/quadrangles en 3D).Ce chargement n’est utilisable que pour les modélisations non saturés de typeTHH2MS,THH2S,HH2MSou HH2S. ◆ PRE1_EXT Valeur de la pression capillaire extérieure. .. code-block:: text ◇ PRE2_EXT Valeur de la pression de gaz extérieure (0 par défaut). .. code-block:: text ◆ COEF_11 coefficient d'échange C11 (cf. eq. ) reliant le flux d'eau à la pression capillaire. .. code-block:: text ◇ COEF_12 coefficient d'échange C12 (cf. eq. ) reliant le flux d'eau à la pression de gaz (0 par défaut). .. code-block:: text ◇ COEF_21 coefficient d'échange C21 (cf. eq. ) reliant le flux d'air à la pression capillaire (0 par défaut). .. code-block:: text ◇ COEF_22 coefficient d'échange C22 (cf. eq. ) reliant le flux d'air à la pression de gaz (0 par défaut). Mot-clé ECHANGE_THM_HR ---------------------- .. code-block:: text ECHANGE_THM_HR = _F( ◆ / TOUT = "OUI" (ou non renseigné), / GROUP_MA = grma, ◆ HR_EXT = float, ◆ ALPHA = float, ◆ PVAP_SAT = float, ) Le mot-clé facteur ECHANGE_THM_HR permet pour une modélisation non saturée d’appliquer à un domaine de milieu continu 2D ou 3D une condition mixte d'échange sur la densité de vapeur (proportionelle à l'humidité relative). Le flux d'eau s'écrit alors : .. math:: \left({M}_{w}+{M}_{\text{vp}}\right).n= \alpha.\left( \rho_{vp} -{\rho_{vp}}^{ext} \right) avec :math:`\rho_{vp}` la densité de vapeur. Sachant que l'humidité relative HR est fonction de la densité via la densité de vapeur saturante :math:`\rho_{vp}^{sat}` : .. math:: \rho_{vp} = HR.\rho_{vp}^{sat} Le flux s'exprime alors .. math:: \left({M}_{w}+{M}_{\text{vp}}\right).n= \alpha.\rho_{vp}^{sat}\left( HR -{HR}^{ext} \right) On peut exprimer cette équation en fonction de la pression capillaire (Pc qui est une inconnue du modèle) grace à la relation de Kelvin soit : .. math:: HR = exp\left[ \frac{Pc.M_{H2O}}{\rho_{l}.R.T}\right] avec :math:`{M}_{H2O}` la masse molaire de l'eau, R la constant des gaz parfaits et T la température. In fine l'expression du flux devient : .. math:: \left({M}_{w}+{M}_{\text{vp}}\right).n= \alpha.\rho_{vp}.\left( exp\left[ \frac{Pc.M_{H2O}}{\rho_{l}.R.T}\right] -{HR}^{ext} \right) Affectation topologique: .. code-block:: text ◆ / TOUT = "OUI" (ou non renseigné), / GROUP_MA = grma, .. code-block:: text ◆ HR_EXT Humidité relative extérieure. .. code-block:: text ◆ ALPHA Coefficient d'échange ◆ PVAP_SAT Pression de vapeur saturante Mot-clé FORCE_SOL ----------------- .. code-block:: text FORCE_SOL = _F( ◆ / GROUP_NO_INTERF = grno, / SUPER_MAILLE = ma, ◆ | UNITE_RESU_MASS = unit, | UNITE_RESU_RIGI = unit, | UNITE_RESU_AMOR = unit, ◇ UNITE_RESU_FORC = unit, ◇ NB_PAS_TRONCATURE = int, ◇ TYPE = / "BINAIRE", / "ASCII" (par défaut), ) Le mot-clé facteur FORCE_SOL permet de prendre en compte la force interne d’un domaine de sol en utilisant les évolutions temporelles des contributions en rigidité, masse et amortissement de l’impédance de sol. L’impédance de sol extraite à l’instant initial permet de constituer par MACR_ELEM_DYNA un macro-élément représentant le comportement du domaine de sol que l’on ajoute au modèle de structure. L’interface dynamique du macro-élément est décrite soit par une super-maille du modèle contenant à la fois la structure et ce macro-élément, soit par un groupe de nœuds si l’interface physique coïncide avec l’interface dynamique modale. On peut également prendre en compte, si elle existe, l’évolution temporelle des forces sismiques, affectée à cette même interface dynamique sous forme d’unité logique . Ce type de charge est pris en compte dans la commande DYNA_NON_LINE. Un exemple d’utilisation est fourni dans le cas test MISS03B [:external:ref:`V1.10.122 `]. **Il ne peut y avoir qu’une seule occurrence de ce mot-clef dans** AFFE_CHAR_MECA *\*.** Affectation topologique: .. code-block:: text ◆ / GROUP_NO_INTERF = grno, / SUPER_MAILLE = ma, Ces opérandes permettent de décrire l’interface dynamique du macro-élément représentant le comportement du domaine de sol que l’on ajoute au modèle de structure soit par une super-maille du modèle contenant à la fois la structure et ce macro-élément par le mot clé SUPER_MAILLE, soit par un groupe de nœuds par le mot clé GROUP_NO_INTERF si l’interface physique coïncide avec l’interface dynamique modale. .. code-block:: text ◆ | UNITE_RESU_MASS = unit, | UNITE_RESU_RIGI = unit, | UNITE_RESU_AMOR = unit, Ces opérandes permettent d’introduire les évolutions temporelles des contributions en rigidité, masse et amortissement de l’impédance de sol sous forme d’unités logiques. .. code-block:: text ◇ TYPE = / "BINAIRE", / "ASCII" (par défaut) Cet opérande permet d’introduire le format d'écriture commun ('ASCII' ou 'BINAIRE') des évolutions temporelles des contributions en rigidité, masse et amortissement de l’impédance de sol renseignées sous forme d’unités logiques par les opérandes précédents UNITE_RESU_* . La valeur de l'opérande doit être cohérente avec la valeur entrée dans l'opérande TYPE_FICHIER_TEMPS lors du calcul de l'option TYPE_RESU='FICHIER_TEMPS' par CALC_MISS. .. code-block:: text ◇ UNITE_RESU_FORC Cet opérande permet d’introduire, si elle existe et sous forme d’unité logique, l’évolution temporelle des forces sismiques, affectée à l’interface dynamique du macro-élément représentant le comportement du domaine de sol que l’on ajoute au modèle de structure. .. code-block:: text ◇ NB_PAS_TRONCATURE = int, Cet opérande permet d’introduire (s'il est présent) une troncature du calcul de la contribution au second membre de la force de sol dans DYNA_NON_LINE. Cette troncature s'opère sur la sommation à partir des champs de déplacements stockés aux instants précédents. Suite à des tests, on se rend compte qu'il suffit souvent de tronquer la somme sur les 50 à 100 derniers pas de temps au lieu de faire la somme complète sur tous les instants stockés précédents, qui est l'option par défaut. Il est donc conseillé d'utiliser comme valeur pratique la valeur entière du maximum entre 50 et le dixième du nombre total de pas de temps stockés.