User Story
Dans le premier exemple, nous procédons à une conception préliminaire, en nous concentrant particulièrement sur le calcul de la force totale. Il est sélectionné comme un plan carré 2D. Il appartient au Groupe 1 du WTG-Merkblatt-M3 :
- G1 : Valeurs qualitatives avec des exigences de précision faibles pour une utilisation dans l'étude de base ou la conception préliminaire. L'effort et les exigences pour le niveau de détail sont réduits, car souvent toutes les conditions limites ne sont pas entièrement clarifiées.
- R1 : Isolée (sans bâtiments environnants), analyse de directions de vent importantes individuelles.
- Z1 : Valeurs moyennes statistiques, à condition qu'elles concernent des processus d'écoulement stationnaires où les fluctuations (par exemple, dues à la turbulence de l'écoulement approchant) peuvent être suffisamment capturées par d'autres mesures.
- S1 : Effets statiques. Il est suffisant de représenter le modèle structurel avec le détail mécanique nécessaire, mais sans propriétés de masse et d'amortissement.
Les dimensions de l'exemple sont montrées à la Figure 1, et l'hypothèse d'entrée est illustrée dans le Tableau 1 :
Table 1 : Données d'entrée de l'exemple de vérification du plan carré 2D
| Modèle | Plan carré 2D |
|---|---|
| Dimension | a = 1 m |
| Vitesse de base du vent | V = 30 m/s |
| Densité de l'air | ρ = 1,225 kg/m³ |
| Solveur | Basé sur la pression |
| Modèle de turbulence | k-ω SST stable |
| Catégorie de terrain | 2 |
| Type de profil de vitesse de vent dans RFEM | Pic |
| Algorithme numérique | Algorithme SIMPLE |
| Discrétisation | Ordre deux |
| Pression résiduelle | 10⁻⁴ |
| Viscosité cinématique | ν = 1,5 × 10⁻⁵ |
Dans cet exemple, nous allons comparer les valeurs de la force du vent entre EN 1991-1-4 et RWIND. La formule de la force du vent dans la section 5.3 de l'Eurocode est définie :
|
cscd |
Structural factor |
|
cf |
Force coefficient for the structure or structural element |
|
qp(ze ) |
Peak velocity pressure at reference height ze |
|
Aref |
Reference area of the structure or structural element |
Les coefficients de force (d=b=1 → Cf,0=2,10) des sections rectangulaires avec des coins vifs et sans écoulement à l'extrémité libre peuvent être obtenus dans la Figure 7.23 de l'EN 1991-1-4, et le facteur de réduction (ψr) pour une section carrée à coins arrondis (r/b=0 → ψ r=1) peut être obtenu de la Figure 7.24 de l'EN 1991-1-4. Des valeurs indicatives du facteur d'effet d'extrémité ψλ=0,63 en fonction du rapport de compacité φ=1 par rapport à la finesse λ=2 peuvent être obtenues dans la Figure 7.36 de l'EN 1991-1-4.
Le coefficient de force cf des éléments structurels de la section rectangulaire avec le vent soufflant normalement sur une surface plate doit être déterminé par l'Expression (7.9) de l'EN 1991-1-4 :
|
cf,0 |
The force coefficient of rectangular sections with sharp corners and without free-end flow |
|
ψr |
The reduction factor for square sections with rounded corners |
|
ψλ |
The end-effect factor for elements with free-end flow |
- Vitesse moyenne du vent
La vitesse moyenne du vent vm (ze) à la hauteur de référence z_e dépend de la rugosité du terrain, de l'orographie du terrain, et de la vitesse de base du vent vb. Elle est déterminée à l'aide de l'équation (4.3) de l'EN1991-1-4 :
- Turbulence du vent
L'intensité de la turbulence Iv (ze) à la hauteur de référence ze est définie comme l'écart type de la turbulence divisé par la vitesse moyenne du vent. Elle est calculée conformément à l'équation 4.7 de l'EN1991-1-4. Pour le cas examiné, ze est inférieur à zmin :
- Pression de base de la vitesse
La pression de base de la vitesse q_b est la pression correspondant à la quantité de vent déterminée à la vitesse de base du vent vb. La pression de base de la vitesse est calculée selon la relation fondamentale spécifiée dans l'EN1991 -14 §4.5(1) :
où ρ est la densité de l'air conformément à l'EN1991-1-4 §4.5(1). Dans ce calcul, la valeur suivante est considérée ρ=1,225 kg/m3.
- Pression de vitesse de pointe
La pression de vitesse de pointe qp (ze) à la hauteur de référence ze inclut les fluctuations de vitesse moyennes et à court terme. Elle est déterminée selon l'équation 4.8 de l'EN1991-1-4 :
Ensuite, la force du vent peut être calculée :
L'image 2 montre une étude de sensibilité du maillage dans RWIND pour la plaque 2D. À mesure que la densité du maillage augmente de 10 % à 40 %, le coefficient de force Cf diminue et se stabilise à 1,23 à partir de 30 %, indiquant des résultats indépendants du maillage et assurant la précision de la simulation sans raffinement inutile.
De plus, l'étude du maillage numérique doit être effectuée selon le lien suivant :
Le WTG-Merkblatt M3 fournit deux méthodes clés pour valider les résultats de simulation. La méthode de la Taux de Réussite évalue combien des valeurs simulées Pi correspondent correctement aux valeurs de référence Oi dans une tolérance définie, en utilisant une approche de classification binaire (succès ou échec). Cette approche évalue la fiabilité de la simulation en calculant un taux de réussite q, similaire aux fonctions de confiance utilisées en théorie de la fiabilité. En revanche, la méthode Erreur Quadratique Moyenne Normalisée (e2) offre une évaluation de la précision plus détaillée en quantifiant la déviation quadratique moyenne entre les valeurs simulées et de référence, normalisée pour tenir compte des différences d'échelle. Ensemble, ces méthodes fournissent des mesures qualitatives et quantitatives pour la validation de la simulation.
Résultats dans RWIND et Comparaison à l'Eurocode
Dans RWIND, les résultats des forces totales (dans les Figures 3 et 4) sont disponibles dans l'onglet Info du modèle d'édition. La différence entre RWIND et l'Eurocode est d'environ Wrel = 2,45% (inférieure au critère mentionné dans WTG) ; alors le taux de réussite peut être obtenu comme q=100%, ce qui montre un bon accord. La faible erreur quadratique moyenne normalisée e²=0,0005 confirme un fort accord entre la simulation et les mesures, répondant efficacement aux normes de validation.
Voici le modèle 3D de la plaque 2D :