Lois de matériau orthotrope

La loi de Hooke, avec $\mathrm u\;=\;\frac{\mathrm E}{2\;\mathrm G}\;\;-\;1\;\leq\;0.5$ and $-1\;<\;\mathrm u\;<\;0.5$ n’est en général pas valide pour les matériaux orthotropes.

Allongement transveral

Les paramètres de matériau suivants sont relatifs à des rigidités bidirectionnelles et, à moins que le contraire ne soit précisé, de matériau bois. Comme affiché dans la Figure 01, nous utilisons un système d’axe local.

• E_x = rigidité dans la direction x locale de la surface
• E_y = rigidité dans la direction y locale de la surface
• G_xz = rigidité de cisaillement dans la direction x locale de la surface (direction d’épaisseur de la plaque)
• G_yz = rigidité de cisaillement dans la direction y locale de la surface (direction d’épaisseur de la plaque)
• G_xy = rigidité de cisaillement dans la zone dans le plan
• ν_xy = déformation transversale dans la direction x
• ν_yx = déformation transversale dans la direction y

Les contraintes dans la Figure 02 sont relatives aux rigidités mentionnées ici.

contraintes de cisaillement

La loi de matériau est sujette aux règles suivantes.

Équation 1:

Équation 2:
${\mathrm u}_\mathrm{xy}\;\geqslant\;{\mathrm u}_\mathrm{yx}$

Équation 3:

Équation 4 (rigidités dans la zone du plan):

Le ratio de déformations dans les équations ci-dessus souligne les relations dans la Figure 01.

Les rigidités dans la zone dans le plan sont calculées comme suit.

Équation 5:

Déformation transversale v

Comme expliqué dans la Figure 01, les déformations etcontraintes modifiées dans cette direction résultent d’un comportement dematériau lissé dans la direction respective.

Ratio des déformations :

Équation 6:

Équation 7:

Pour ${\mathrm\sigma}_\mathrm x\; eq\;0,\;{\mathrm\sigma}_\mathrm y\;=\;0,\;{\mathrm\sigma}_\mathrm{xy}\;=\;0$, les équations suivantes sont données avec la loi de Hooke.

Équation 8:

Équation 9:

Équation 10:

Équation 11:

Équation 12:

Équation 13:

Matrice de rigidité

Calcul de la matrice de rigidité globale de la plaque.

Équation 14:

Composants du fléchissement:

Équation 15:
${\mathrm D}_{11}\;=\;\frac{\mathrm d^3}{12}\;\cdot\;{\mathrm d}_{11}\;\widehat{=\;}\frac{{\mathrm E}_\mathrm x\;\mathrm d^3}{12\;\left(1\;-\;{\mathrm\nu}_\mathrm{xy}\;{\mathrm\nu}_\mathrm{yx}\right)}$

Équation 16:
${\mathrm D}_{12}\;=\;\frac{\mathrm d^3}{12}\;\cdot\;{\mathrm d}_{12}\;\widehat{=\;}\mathrm{sgn}\left({\mathrm\nu}_\mathrm{xy}\right)\;\sqrt{{\mathrm\nu}_\mathrm{xy}\;{\mathrm\nu}_\mathrm{yx}\;{\mathrm D}_{11}\;{\mathrm D}_{22}}$

Équation 17:
${\mathrm D}_{22}\;=\;\frac{\mathrm d^3}{12}\;\cdot\;{\mathrm d}_{22}\;\widehat{=\;}\frac{{\mathrm E}_\mathrm y\;\mathrm d^3}{12\;\left(1\;-\;{\mathrm\nu}_\mathrm{xy}\;{\mathrm\nu}_\mathrm{yx}\right)}$

Équation 18:
${\mathrm D}_{33}\;=\;\frac{\mathrm d^3}{12}\;\cdot\;{\mathrm d}_{33}\;\widehat{=\;}{\mathrm G}_\mathrm{xy}\;\frac{\mathrm d^3}{12}$

Composants de membrane:

Équation 19:
${\mathrm D}_{66}\;=\;\mathrm d\;\cdot\;{\mathrm d}_{11}\;\widehat{=\;}\frac{{\mathrm E}_\mathrm x\;\mathrm d}{1\;-\;{\mathrm u}_\mathrm{xy}\;{\mathrm u}_\mathrm{yx}}$

Équation 20:
${\mathrm D}_{67}\;=\;\mathrm d\;\cdot\;{\mathrm d}_{12}\;\widehat{=\;}\mathrm{sgn}\left({\mathrm\nu}_\mathrm{xy}\right)\;\sqrt{{\mathrm\nu}_\mathrm{xy}\;{\mathrm\nu}_\mathrm{yx}\;{\mathrm D}_{66}\;{\mathrm D}_{77}}$

Équation 21:
${\mathrm D}_{77}\;=\;\mathrm d\;\cdot\;{\mathrm d}_{22}\;\widehat{=\;}\frac{{\mathrm E}_\mathrm y\;\mathrm d}{1\;-\;{\mathrm u}_\mathrm{xy}\;{\mathrm u}_\mathrm{yx}}$

Équation 22:
${\mathrm D}_{88}\;=\;\mathrm d\;\cdot\;{\mathrm d}_{33}\;\widehat{=\;}{\mathrm G}_\mathrm{xy}\;\mathrm d$

Composants du cisaillement:

Équation 23:

Équation 24:

La matrice de rigidité doit être définie positive pour ces équations, de sorte que toutes les valeurs propres de la matrice soient positives.

C’est pour cette raison que RFEM contrôle, entre autres, la définition de la déformation transversale relative à l’équation suivante.

Équation 25 :

Exemple

L’exemple suivant (Figure 03) explique le comportement de matériau orthotrope. Un matériau orthotrope sera comparé à un matériau isotrope. De plus, la rigidité d’une plaque orthotrope sera définie avec une rigidité élevée en direction x et y.

• Visionneuse 3D | Google
• Visionneuse 3D | Babylon
• Réalité virtuelle

Système

Structure:

• Épaisseur de plaque 200 mm
• Matériau C 24
• Rigidités orthotropes
  
  $\begin{array}{l}{\mathrm{E}}_{\mathrm{x}} = 1,100.0 \mathrm{kN}/{\mathrm{cm}}^2\\ {\mathrm{E}}_{\mathrm{y}} = 37.0 \mathrm{kN}/{\mathrm{cm}}^2\\ {\mathrm{G}}_{\mathrm{y}} = 6.9 \mathrm{kN}/{\mathrm{cm}}^2\\ {\mathrm{G}}_{\mathrm{x}} = 69.0 \mathrm{kN}/{\mathrm{cm}}^2\\ {\mathrm{G}}_{\mathrm{xy}} = 69.0 \mathrm{kN}/{\mathrm{cm}}^2\\ {\mathrm{u}}_{\mathrm{xy}} = 2.52\\ {\mathrm{u}}_{\mathrm{yx}} = 0.085\end{array}$

  Formule 29
• Rigidités isotropes
  
  $\begin{array}{l}\mathrm{E} = 1,100 \mathrm{kN}/{\mathrm{cm}}^2\\ \mathrm{G} = 500 \mathrm{kN}/{\mathrm{cm}}^2\\ \mathrm{u} = 0.1\end{array}$

  Formule 30
• Dimension l = 2.0 m, L = 4.0 m
• Charge: 20 kN/m²
• Taille du maillage EF: 50 cm

La structure est définie comme encastrée dans la direction z verticale. Les conditions d’appui dans les directions x et y ont été sélectionnées de sorte qu’aucun effet dû aux maintiens n’ait lieu.

Le calcul est réalisé selon l’analyse statique linéaire avec des conditions d’appui et un comportement de matériau linéaire élastique.

La déformation transversale suivante est un résultat de la loi de Hooke appliqué aux valeurs données.

Équation 26 :
$\mathrm\nu\;=\;\left(\sqrt{{\mathrm E}_\mathrm x\;\cdot\;{\mathrm E}_\mathrm y}\right)\;/\;\left(2\;\mathrm G\right)\;-\;1\;=\;6.97$

Cette déformation transversale élevée ne peut pas être atteinte avec le modèle sélectionné. Avec les équations de [1], les valeurs peuvent toutefois être ajustées.

Équation 27 :

Équation 28 :

Équation 29 :

Équation 30 :

• Visionneuse 3D | Google
• Visionneuse 3D | Babylon
• Réalité virtuelle

Surfaces et rigidités

Résultat
Comme prévu, les plus grandes déformations ont lieu avec l’orientation des rigidités dans la direction y (Figure 06). Les réactions d’appui et le moment de la plaque isotrope sont affichés dans la Figure 05.

• Visionneuse 3D | Google
• Visionneuse 3D | Babylon
• Réalité virtuelle

Résultats de matériau isotrope

La plaque avec une rigidité importante dans la direction y (Ey = 1 100 kN/cm^{2) ayant une résistance importante dans cette direction, les réactions d’appui sont également plus élevées (125,4 kNm contre 58,3 kNm).}

• Visionneuse 3D | Google
• Visionneuse 3D | Babylon
• Réalité virtuelle

Déformations et réactions d'appui des plaques orthotropes

Les moments de flexion maximum résultants pour les plaques orthotropes sont égaux à m_x avec la rigidité en direction x et pour m_y avec une rigidité importante en direction y.

La plaque à la rigidité élevée dans la direction-y a son moment fléchissant m_y quasiment en son centre (Figure 07).

Biegemomente der orthotropen Platte mx (links) und my (rechts)

Variation de la déformation transversale

La déformation transversale relative aux diagrammes de déformation peut atteindre les valeurs maximale et minimale listées dans ce tableau.

cette opération est possible.

La plaque à rigidité importante vue précédemment (E_x = 11,000) sera définie avec ces déformations transversales élevées. Les autres rigidités de plaque restent toutefois identiques.

La Figure 08 affiche les résultats de la variation de ν_xy = 5,44 à -5,44.

Pour ν_xy = 5,44 les réactions d’appui sont qualitativement identiques à celles du comportement de matériau isotrope. Le moment fléchissant passe de m_x = 18.1 kNm/m (plaque isotrope) à m_x = 34.9 kNm/m (plaque orthotrope).

Le moment fléchissant est légèrement réduit en comparaison aux déformations transversales habituelles (ν_xy = 2.5).

Ergebnisse der Variation gemäß Tabelle (links: νxy = 5,44, Milieu: νxy = 0, à droite: νxy = -5,44)

Avec ν_xy = 0, la grande amplitude de valeurs des réactions d’appui de la fin libre de plaque est remplacée par une valeur constante de 43 kN/m.

Le moment m_x augmente à 38.1 kNm/m. Relativement au résultat précédent (ν_xy =5.44), l’influence de la déformation transversale est perceptible visuellement. Pour v = 0, aucune déformation ou distorsion n’est provoquée par la déformation transversale.

Pour ν_xy = -5,44 un échec post-critique est affiché en fin de plaque et les réactions d’appui deviennent alors négatives. Le moment maximum a lieu dans le centre de la plaque avec 59,5 kNm/m.

La plaque a maintenant un comportement plus proche de celui d’une plaque sous contraintes uniaxales sans troisième appui dans sa direction longitudinale.

Ce comportement peut être expliqué avec la Figure 01 et la relation qui y est listée.

La déformation transversale négative élevée (ν_xy = -5,44), provoque un échec sous charge de la bordure libre et ne peut ainsi pas être déformée.

L’influence de l’orthotropie dans la direction y est ici quasi-nulle (E_y ≈ 0).

Résumé

Le type de matériau orthotrope permet la définition de presque tous les paramètres de matériau dans RFEM. La variation des déformations transversales permet d’obtenir des résultats très différents. Une déformation transversale suite à la modification des valeurs selon [1] résulte de valeurs proches de la solution pour une poutre à travée simple.

Équation 31:

Des déformations transversales négatives trop élevées affichent un système structural modifié qui ne correspond pas à la modélisation.