v2.02.010 SDLL10 - Poutre de section rectangulaire variable (encastrée-encastrée)#
Résumé:
Ce problème plan consiste à chercher les fréquences et les modes de vibration d’une structure mécanique composée d’une poutre encastrée-encastrée dont l’aire de la section droite varie de façon exponentielle. Ce test de Mécanique des Structures correspond à une analyse dynamique d’un modèle linéique ayant un comportement linéaire. Il comprend une seule modélisation.
Par l’intermédiaire de ce problème, on teste l’élément de poutre en flexion de Timoshenko de section variable ainsi que le calcul des fréquences et des modes de vibration par la méthode de Lanczos. On teste aussi la fonctionnalité « norme à 1. » au point d’amplitude maximale en translation » des modes de vibration.
En utilisant une discrétisation spatiale fine, les résultats obtenus sont en bon accord avec les résultats analytiques donnés dans le guide VPCS.
Solution de référence#
Méthode de calcul utilisée pour la solution de référence#
La solution de référence est celle donnée dans la fiche SDLL10/89 du guide VPCS qui présente la méthode de calcul de la façon suivante:
La pulsation \({\omega}_{i}\) est donnée par les racines de l’équation :
avec :
Les composantes de translation \(\nu\) du mode \({F}_{i}(x)\) sont alors :
\({\phi }_{i}(x)={e}^{\alpha x}\left[\cos(\mathit{rx})-\mathit{ch}(\mathit{rx})+\frac{\cos(\mathit{rL})-\mathit{ch}(\mathit{sL})}{r\mathit{sh}(\mathit{sL})-s\sin(\mathit{rL})}(s\sin(\mathit{rx})-r\mathit{sh}(\mathit{sx}))\right]\)
Résultats de référence#
4 premières fréquences propres et modes propres normés à 1 pour la plus grande composante en translation.
Incertitude sur la solution#
Solution analytique.
Références bibliographiques#
Groupe de travail Analyse Dynamique. Comité de Validation des Progiciels de Calcul de Structure. Société Française des Mécaniciens (1988).
Modélisation A#
Caractéristiques de la modélisation#
Élément de poutre POU_D_T (Poutre droite de Timoshenko)
Découpage : poutre \(\mathrm{AB}\) : 120 mailles SEG2 de section variable.
Conditions limites: |
|
en tous les nœuds |
|
aux nœuds |
|
extrémités |
`` ( NOEUD: (AB) DX: 0., DY: 0., DRZ: 0. )`` |
Noms des nœuds: |
Point \(A\) |
\(x=0.\) |
\(\mathrm{N1}\) |
\(x=0.1\) |
\(\mathrm{N21}\) |
||
\(x=0.2\) |
\(\mathrm{N41}\) |
||
\(x=0.3\) |
\(\mathrm{N61}\) |
||
\(x=0.4\) |
\(\mathrm{N81}\) |
||
\(x=0.5\) |
\(\mathrm{N101}\) |
||
Point \(B\) |
\(x=0.6\) |
\(\mathrm{N121}\) |
Caractéristiques du maillage#
Nombre de nœuds : |
121 |
Nombre de mailles et types : |
120 SEG2 |
Grandeurs testées et résultats#
Ordre du mode propre |
Fréquence Référence |
% tolérance |
1 |
143.303 |
1.6 |
2 |
396.821 |
0.45 |
3 |
779.425 |
0.3 |
4 |
1289.577 |
0.9 |
Les modes propres de Aster ont été normés à 1. au point d’amplitude maximum en translation comme dans la référence.
Mode propre \({F}_{i}(x)\) normé à 1 au point d’amplitude maximale
\(i\) |
\(x=0.1\) |
\(x=0.2\) |
\(x=0.3\) |
\(x=0.4\) |
\(x=0.5\) |
|
Référence Aster % différence % tolérance |
1 |
0.2349 0.2363 0.583 0.6 |
0.6962 0.6970 0.119 0.15 |
0.98960 0.9895 0. 0.1 |
0.8505 0.8516 0.132 0.15 |
0.3507 0.3529 0.631 0.7 |
Référence Aster % différence % tolérance |
2 |
–0.4653 –0.4670 0.37 0.4 |
–0.7558 –0.7555 –0.041 0.1 |
0. –2.910–4 –2.910–4 1.10–3 |
0.9232 0.9226 –0.063 0.1 |
0.6941 0.6971 0.435 0.45 |
Référence Aster % différence % tolérance |
3 |
0.6278 0.6290 0.192 0.2 |
0.1969 0.1952 –0.89 0.9 |
–0.7783 –0.7782 –0.014 0.1 |
0.2406 0.2377 –1.226 1.23 |
0.9366 0.9387 0.228 0.25 |
Référence Aster % différence % tolérance |
4 |
–0.666 –0.6656 –0.081 0.1 |
0.4832 0.4840 0.18 0.2 |
0. 4.610–4 4.610–4 1.10–3 |
–0.5901 –0.5919 0.31 0.35 |
0.9937 0.9928 –0.089 0.1 |
Remarques#
Calculs effectués par :
CALC_MODES
OPTION = “PLUS_PETITE” CALC_FREQ=_F(NMAX_FREQ = 4) SOLVEUR_MODAL=_F(METHODE = “TRI_DIAG”)
Contenu du fichier résultats :
4 premières fréquences propres, vecteurs propres.
Synthèse des résultats#
Modélisation convenable (fréquences et modes propres à moins de \(\text{2\%}\) ) avec un maillage fin.
Un calcul effectué sur un maillage grossier (12 mailles) montre des écarts plus importants avec la solution de référence. Ceci est surtout dû à la façon dont sont normés les modes.