v6.03.108 SSNP108 - Élément de béton précontraint en compression#
Résumé:
On considère une structure élémentaire constituée d’une plaque de béton carrée traversée par un câble de précontrainte dont la fibre neutre est confondue avec l’axe de symétrie horizontal de la plaque. Le bord vertical gauche de la plaque est fixe. Le câble est mis en traction à ses deux extrémités afin de pré-contraindre la plaque en compression. Les pertes de tension le long du câble sont négligées.
Le but de ce cas-test est de valider, sur une configuration simple, la méthode de calcul de l’état d’équilibre d’une structure de béton précontrainte. Les résultats sont validés par comparaison à une solution analytique.
Solution de référence#
Solution formelle#
Les déplacements dans le câble et dans le béton sont continus et homogènes. Ainsi le déplacement horizontal sur l’intervalle \([0;L]\) vaut \(u(x)=\frac{{u}_{L}}{L}x\)
Les contraintes normales dans la plaque de béton et le câble d’acier s’écrivent respectivement (hypothèse d’élasticité):
\(\lbrace \begin{array}{c}{\sigma}_{b}={E}_{b}\frac{\partial u}{\partial x}={E}_{b}\frac{{u}_{L}}{L}\\ {\sigma}_{a}={E}_{a}\frac{\partial u}{\partial x}+{\sigma}_{0}={E}_{a}\frac{{u}_{L}}{L}+{\sigma}_{0}\end{array}\)
où \({\sigma}_{0}=\frac{{F}_{0}}{{S}_{a}}\) est la précontrainte initiale dans le câble
et \({u}_{L}\) est le déplacement horizontal à l’abscisse \(L\) .
L’équilibre de l’ensemble plaque et câble s’écrit: \({\sigma}_{b}{S}_{b}+{\sigma}_{a}{S}_{a}=0\mathrm{\Rightarrow }{u}_{L}=\frac{-L{F}_{0}}{{E}_{b}eH+{E}_{a}{S}_{a}}\)
L’effort normal dans le câble vaut: \({N}_{a}={E}_{a}{S}_{a}\frac{{u}_{L}}{F}+{F}_{0}={F}_{0}\frac{{E}_{b}eH}{{E}_{b}eH+{E}_{a}{S}_{a}}\)
L’effort total sur la section verticale de la plaque de béton vaut:
\({N}_{b}={E}_{b}eH\frac{{u}_{L}}{L}=-{F}_{0}\frac{{E}_{b}eH}{{E}_{b}eH+{E}_{a}{S}_{a}}\)
On en déduit la densité linéique d’effort normal sur la plaque de béton
\({N}_{xx}=-{F}_{0}\frac{{E}_{b}e}{{E}_{b}eH+{E}_{a}{S}_{a}}\)
Valeurs numériques de référence#
Les valeurs numériques de référence sont:
\({u}_{L}=–1,11013974{.10}^{-5}m\) |
\({N}_{a}=1,99825153{.10}^{5}N\) |
\({N}_{b}=–1,99825153{.10}^{5}N\) |
\({N}_{xx}=–9,99125765{.10}^{4}N/m\) |
Modélisation A#
Caractéristiques de la modélisation#
La figure ci-dessous donne une représentation simplifiée du maillage.
La plaque de béton est représentée par un élément DKT, supporté par une maille quadrangle à 4 nœuds.
Une épaisseur \(e=0,6m\) lui est affectée, ainsi qu’un matériau béton pour lequel sont définis les comportements ELAS (module d’Young \({E}_{b}={3.10}^{10}\mathit{Pa}\) ) et BPEL_BETON: les paramètres caractéristiques de cette relation sont fixés à 0 car on néglige les pertes de tension le long du câble de précontrainte.
Le nœud \(\mathit{NB001001}\) est encastré: DX, DY, DZ, DRX, DRY et DRZ sont bloqués. Le nœud \(\mathit{NB002001}\) est appuyé bi-latéralement: DX et DZ sont bloqués.
Le câble est représenté par 4 éléments MECA_BARRE, supportés par 4 mailles segments à 2 nœuds.
Une aire de section droite \({S}_{a}=1,5{.10}^{-4}{m}^{2}\) leur est affectée, ainsi qu’un matériau acier pour lequel sont définis les comportements ELAS (module d’Young \({E}_{a}=2,1{.10}^{11}\mathit{Pa}\) ) et BPEL_ACIER: les paramètres caractéristiques de cette relation sont fixés à 0 (pertes de tension négligées), à l’exception de la contrainte limite élastique pour laquelle une valeur nulle est illicite (\({f}_{\mathit{prg}}=1,77{.10}^{9}\mathit{Pa}\) ).
La tension \({F}_{0}={2.10}^{5}N\) est appliquée aux nœuds \(\mathit{NC001001}\) et \(\mathit{NC001005}\) . Cette valeur de tension est cohérente avec les valeurs de section et de limite élastique, pour un câble de précontrainte de type toron.
Le calcul de l’état d’équilibre de l’ensemble plaque et câble est effectué en un seul pas, le comportement étant élastique.
Étapes de calcul et fonctionnalités testées#
Les principales étapes de calcul correspondent aux fonctionnalités que l’on souhaite valider:
opérateur DEFI_MATERIAU : définition des relations de comportement BPEL_BETON et BPEL_ACIER, dans le cas particulier où les pertes de tension le long du câble de précontrainte sont négligées (valeurs par défaut des paramètres);
opérateur DEFI_CABLE_BP: détermination d’un profil de tension constant le long du câble de précontrainte, les pertes étant négligées; calcul des coefficients des relations cinématiques entre les degrés de liberté des nœuds du câble et les degrés de liberté des nœuds «voisins» de la plaque de béton;
opérateur AFFE_CHAR_MECA: définition d’un chargement de type RELA_CINE_BP;
opérateur STAT_NON_LINE, option COMPORTEMENT: calcul de l’état d’équilibre en tenant compte du chargement de type RELA_CINE_BP.
Résultats de la modélisation A#
Valeurs testées#
Densité linéique d’effort normal sur la section verticale de la plaque de béton#
On compare les valeurs extraites du champ SIEF_ELNO issu de STAT_NON_LINE aux valeurs théoriques de référence. L’extraction se fait sur la maille \(\mathit{QD001001}\) représentant la plaque de béton.
La composante sur laquelle portent les tests est \(\mathit{NXX}\) .
La tolérance d’écart relatif par rapport à la référence vaut 10–6 %.
Nœud |
Valeur de référence |
Tolérance |
NB001001 |
–9,99125765.104N/m |
–4,43.10–11% |
NB001002 |
–9,99125765.104N/m |
–4,43.10–11% |
NB002001 |
–9,99125765.104N/m |
–4,43.10–11% |
NB002002 |
–9,99125765.104N/m |
–4,43.10–11% |
Déplacement horizontal des nœuds de la plaque de béton#
On compare les valeurs extraites du champ DEPL issu de STAT_NON_LINE aux valeurs théoriques de référence.
La composante sur laquelle portent les tests est \(\mathit{DX}\) .
La tolérance d’écart relatif par rapport à la référence vaut 10–6 %.
Nœud |
Valeur de référence |
Tolérance |
NB001002 |
–1,11013974.10–5m |
–1,05.10–9% |
NB002002 |
–1,11013974.10–5m |
–1,05.10–9% |
Effort normal dans le câble#
On compare les valeurs extraites du champ SIEF_ELNO issu de STAT_NON_LINE aux valeurs théoriques de référence. L’extraction se fait sur les mailles \(\mathit{SG001001}\) pour le nœud \(\mathit{NC001001}\) , \(\mathit{SG001002}\) pour le nœud \(\mathit{NC001002}\) , \(\mathit{SG001003}\) pour le nœud \(\mathit{NC001003}\) , et \(\mathit{SG001004}\) pour les nœuds \(\mathit{NC001004}\) et \(\mathit{NC001005}\) .
La composante sur laquelle portent les tests est \(N\) .
La tolérance d’écart relatif par rapport à la référence vaut 10–6 %.
Nœud |
Valeur de référence |
Tolérance |
\(\mathit{NC001001}\) |
1,99825153.105N |
–4,44.10–11% |
\(\mathit{NC001002}\) |
1,99825153.105N |
–4,44.10–11% |
\(\mathit{NC001003}\) |
1,99825153.105N |
–4,44.10–11% |
\(\mathit{NC001004}\) |
1,99825153.105N |
–4,44.10–11% |
\(\mathit{NC001005}\) |
1,99825153.105N |
–4,44.10–11% |
Déplacement horizontal des nœuds du câble#
On compare les valeurs extraites du champ DEPL issu de STAT_NON_LINE aux valeurs théoriques de référence.
La composante sur laquelle portent les tests est \(\mathit{DX}\) .
La tolérance d’écart relatif par rapport à la référence vaut 10–6 %.
Nœud |
Valeur de référence |
Tolérance |
NC001002 |
–2,77534935.10–6m |
–1,05.10–9% |
NC001003 |
–5,55069869.10–6m |
7,56.10–10% |
NC001004 |
–8,32604804.10–6m |
1,56.10–10% |
NC001005 |
–1,11013974.10–5m |
–1,05.10–9% |
Remarques#
Les valeurs calculées correspondent effectivement à celles théoriquement attendues. On obtient bien un état de compression pour la plaque de béton.
On observe une différence infinitésimale entre le déplacement horizontal au nœud \(\mathit{NC001005}\) appartenant au câble et le déplacement horizontal aux nœuds \(\mathit{NB001002}\) et \(\mathit{NB002002}\) de la plaque de béton. Les valeurs relevées devraient être identiques, mais les erreurs d’arrondis apparaissant dans les coefficients des relations cinématiques expliquent cette différence infinitésimale.
Modélisation B#
Caractéristiques de la modélisation#
Pour cette modélisation, les nœuds «câble» et «béton» sont confondus.
La figure ci-dessous donne une représentation simplifiée du maillage.
La plaque de béton est représentée par 10 éléments DKT, supportés par 10 mailles QUAD4.
Une épaisseur \(e=0,6m\) lui est affectée, ainsi qu’un matériau béton pour lequel sont définis les comportements ELAS (module d’Young \({E}_{b}={3.10}^{10}\mathit{Pa}\) ) et BPEL_BETON: les paramètres caractéristiques de cette relation sont fixés à 0 car on néglige les pertes de tension le long du câble de précontrainte.
Les nœuds \(\mathit{N1}\) , \(\mathit{N5}\) et \(\mathit{N3}\) sont encastrés: \(\mathit{DX},\mathit{DY},\mathit{DZ},\mathit{DRX},\mathit{DRY}\) et \(\mathit{DRZ}\) sont bloqués.
Le câble est représenté par 5 éléments MECA_BARRE, supportés par 5 mailles SEG2.
Une aire de section droite \({S}_{a}=1,5{.10}^{-4}{m}^{2}\) leur est affectée, ainsi qu’un matériau acier pour lequel sont définis les comportements ELAS (module d’Young \({E}_{a}=2,1{.10}^{11}\mathit{Pa}\) ) et BPEL_ACIER: les paramètres caractéristiques de cette relation sont fixés à 0 (pertes de tension négligées), à l’exception de la contrainte limite élastique pour laquelle une valeur nulle est illicite (\({f}_{\mathit{prg}}=1,77{.10}^{9}\mathit{Pa}\) ).
La tension \({F}_{0}={2.10}^{5}N\) est appliquée aux nœuds \(\mathit{N5}\) et \(\mathit{N6}\) . Cette valeur de tension est cohérente avec les valeurs de section et de limite élastique, pour un câble de précontrainte de type toron.
Le calcul de l’état d’équilibre de l’ensemble plaque et câble est effectué en un seul pas, le comportement étant élastique.
Etapes de calcul et fonctionnalités testées#
Les principales étapes de calcul correspondent aux fonctionnalités que l’on souhaite valider:
opérateur DEFI_MATERIAU : définition des relations de comportement BPEL_BETON et BPEL_ACIER, dans le cas particulier où les pertes de tension le long du câble de précontrainte sont négligées (valeurs par défaut des paramètres);
opérateur DEFI_CABLE_BP: détermination d’un profil de tension constant le long du câble de précontrainte, les pertes étant négligées; calcul des coefficients des relations cinématiques entre les degrés de liberté des nœuds du câble et les degrés de liberté des nœuds «voisins» de la plaque de béton;
opérateur AFFE_CHAR_MECA: définition d’un chargement de type RELA_CINE_BP;
opérateur STAT_NON_LINE, option COMPORTEMENT: calcul de l’état d’équilibre en tenant compte du chargement de type RELA_CINE_BP.
opérateur POST_RELEVE_T, NOM_CMP=‘N’ sur tout le champ SIEF_ELNO.
Résultats de la modélisation B#
Valeurs testées#
Densité linéique d’effort normal sur la section verticale de la plaque de béton#
Nœud |
Valeur de référence |
Tolérance |
\(\mathit{N1}\) |
–9,99125765.104N/m |
1.526 % |
\(\mathit{N3}\) |
–9,99125765.104N/m |
1.526 % |
Effort normal dans le câble#
Nœud |
Valeur de référence |
Tolérance |
\(\mathit{N5}\) |
1,99825153.105N |
–0.001 % |
\(\mathit{N6}\) |
1,99825153.105N |
–0.193 % |
Effort normal dans le câble via la commande POST_RELEVE_T#
Nœud |
Valeur de référence |
Valeur calculée |
\(\mathit{N5}\) |
1,998224892.105N |
1,9982248921222.105N |
\(\mathit{N6}\) |
1,994393252.105N |
1,9943932520206.105N |
Remarques#
Les valeurs calculées correspondent effectivement à celles théoriquement attendues. On obtient bien un état de compression pour la plaque de béton.
Synthèse des résultats#
Les résultats obtenus sont validés par comparaison à une solution analytique de référence avec une très bonne précision.
Les fonctionnalités particulières testées sont les suivantes:
opérateur DEFI_MATERIAU: définition des paramètres caractéristiques des matériaux acier et béton permettant le calcul de la tension le long du câble de précontrainte, suivant les règles du BPEL;
opérateur DEFI_CABLE_BP: calcul de la tension le long du câble et des coefficients des relations cinématiques entre les degrés de liberté des nœuds du câble et les degrés de liberté des nœuds «voisins» de la plaque de béton;
opérateur AFFE_CHAR_MECA: définition d’un chargement de type RELA_CINE_BP;
opérateur STAT_NON_LINE, option COMPORTEMENT: calcul de l’état d’équilibre en tenant compte du chargement de type RELA_CINE_BP.