正交各向异性材料法则

一般来说如果胡克定律 $\mathrm\nu\;=\;\frac{\mathrm E}{2\;\mathrm G}\;\;-\;1\;\leq\;0,5$ und $-1\;<\;\mathrm\nu\;<\;0,5$ ，那么就不适用于各向异性材料。

横向应变

下列材料参数属于二维刚度，一般没有特别说明，指的材料是木材。 基于一个局部坐标系，如图 01 所示。

• E_x = 面上局部 x 方向上的刚度
• E_y = 面上局部 y 方向上的刚度
• G_xz = 面上局部 x 方向上的抗剪刚度（板厚度方向）
• G_yz = 面上局部 y 方向上的抗剪刚度（板厚度方向）
• G_xy = 板平面上的抗剪刚度
• ν_xy = 在 x 方向上的横向应变
• ν_yx = 在 y 方向上的横向应变

在图 02 中的应力和这里的提到的刚度相关。

剪应力

材料法则受如下的限制：

公式 1：

公式 2：
${\mathrm\nu}_\mathrm{xy}\;\geqslant\;{\mathrm\nu}_\mathrm{yx}$

公式 3：

公式 4（在板平面上的刚度）：

上述公式中的应变比强调了图 01 中的关系。

计算在板平面上的刚度如下。

公式 5：

横向应变 ν

如在图 01 中所解释的，在该方向上改变的变形和应力是由相应方向上更平滑的材料行为引起的。

应变比：

公式 6：

公式 7：

如果 ${\mathrm\sigma}_\mathrm x\;\neq\;0,\;{\mathrm\sigma}_\mathrm y\;=\;0,\;{\mathrm\sigma}_\mathrm{xy}\;=\;0$ 从胡克定律得出下面公式。

公式 8：

公式 9：

公式 10：

公式 11：

公式 12：

公式 13：

刚度矩阵

计算板的整体刚度矩阵。

公式 14：

弯曲部分：

公式 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)}$

公式 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}}$

公式 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)}$

公式 18：
${\mathrm D}_{33}\;=\;\frac{\mathrm d^3}{12}\;\cdot\;{\mathrm d}_{33}\;\widehat{=\;}{\mathrm G}_\mathrm{xy}\;\frac{\mathrm d^3}{12}$

板壳部分：

公式 19：
${\mathrm D}_{66}\;=\;\mathrm d\;\cdot\;{\mathrm d}_{11}\;\widehat{=\;}\frac{{\mathrm E}_\mathrm x\;\mathrm d}{1\;-\;{\mathrm\nu}_\mathrm{xy}\;{\mathrm\nu}_\mathrm{yx}}$

公式 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}}$

公式 21：
${\mathrm D}_{77}\;=\;\mathrm d\;\cdot\;{\mathrm d}_{22}\;\widehat{=\;}\frac{{\mathrm E}_\mathrm y\;\mathrm d}{1\;-\;{\mathrm\nu}_\mathrm{xy}\;{\mathrm\nu}_\mathrm{yx}}$

公式 22：
${\mathrm D}_{88}\;=\;\mathrm d\;\cdot\;{\mathrm d}_{33}\;\widehat{=\;}{\mathrm G}_\mathrm{xy}\;\mathrm d$

剪切部分：

公式 23：

公式 24：

当然，这些公式的前提条件一直是刚度矩阵为正，即矩阵的所有特征值也都是正的。

出于这个原因 RFEM 根据下面的公式来检查横向应变的定义。

公式 25：

举例

下面的示例（图 03）展示正交各向异性材料的材料行为。 这里的例子，与最左边的各向异性材料比较， 在右边设置了各向同性材料。 并且，各向异性材料的板上设置在 x 方向上的抗剪刚度远大于 y 方向上的抗剪刚度。

• 3D Viewer | Google
• 3D Viewer | Babylon
• 虚拟现实

系统

结构：

• 板厚 200mm
• 材料 C24
• 正交各向异性刚度
  
  $\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{\nu }}_{\mathrm{xy}} = 2,52\\ {\mathrm{\nu }}_{\mathrm{yx}} = 0,085\end{array}$

  公式 29
• 各向同性刚度
  
  $\begin{array}{l}\mathrm{E} = 1.100 \mathrm{kN}/{\mathrm{cm}}^2\\ \mathrm{G} = 500 \mathrm{kN}/{\mathrm{cm}}^2\\ \mathrm{\nu } = 0,1\end{array}$

  公式 30
• 尺寸 b = 2,0 m，l = 4,0 m
• 荷载 20 kN/m²
• 有限元网格尺寸 50 cm

该结构在垂直的 z 方向是固定不可移动的。 对于在 x 方向和 y 方向上的支座条件，在软件中选择没有约束。

该计算基于一阶理论的线性弹性材料行为和支座条件得出。

通过上面的已知值胡克定律给出了下面的横向应变。

公式 26：
$\mathrm\nu\;=\;\left(\sqrt{{\mathrm E}_\mathrm x\;\cdot\;{\mathrm E}_\mathrm y}\right)\;/\;\left(2\;\mathrm G\right)\;-\;1\;=\;6,97$

选定的材料模型无法实现这种高的横向应变。 但是这些值可以通过 [1] 中的这个公式进行调整。

公式 27：

公式 28：

公式 29：

公式 30：

• 3D Viewer | Google
• 3D Viewer | Babylon
• 虚拟现实

面的刚度

结果
正如预期的那样，在 y 方向上的刚度产生最大的变形（图 06）。 各向同性的板的支座反力和弯矩如图 05 所示。

• 3D Viewer | Google
• 3D Viewer | Babylon
• 虚拟现实

各向同性材料的结果

由于板在 y 方向上具有高刚度 (E_y = 1.100 kN/cm²)，那么在这个方向上就具有高承载力，所以这里的支座反力也很大 (125,4 kNm 到 58,3 kNm)。

• 3D Viewer | Google
• 3D Viewer | Babylon
• 虚拟现实

各向异性板的变形和支座反力

正交各向异性板中的 x 和 y 方向的最大弯矩分别为 m_x 和 m_y，在 y 方向上刚度更大。

板在 y 方向上有更大的刚度，那么最大弯矩 m_y 在板中心的附近（见图 07）。

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

横向应变的变量

横向应变按照应变关系在表格中列出了得出的最大值和最小值。

开始时引入的具有很大刚度 (E_x = 11.000) 的板在这里给它定义一个很大的横向应变。 板的其它刚度保持不变。

图 08 中是由 ν_xy = 5.44 到 -5.44 变量得出的结果。

对于 ν_xy = 5.44 支座反力性质上与各向同性的材料行为是相同的。 弯矩从 m_x = 18.1 kNm/m (isotrope Platte) 增加到 m_x = 34.9 kNm/m（各向异性材料板）。

与具有较常见横向应变 (ν_xy = 2,5) 的正交各向异性材料板相比，弯矩有稍微变小。

Ergebnisse der Variation gemäß Tabelle (links: νxy= 5.44，中间： νxy= 0，右： νxy= -5.44）

当 ν_xy = 0 时板的自由端的支座的偏移值变为近似常量值 43 kN/m。

弯矩 m_x 增大到 38.1 kNm/m。 和前面的结果 (ν_xy =5,44) 相比较这显示了横向应变的影响。 对于 ν =0 横向应变没有引起变形或者应变。

当 ν_xy = -5.44 时在板的自由端出现超临界状况，支座反力为负。 最大弯矩出现在板中心为 59.5 kNm/m。

现在这个板更像一个单轴拉紧的板，在它的纵向上没有第三个支座。

根据图 01 结合这里引入的关系就可以解释这种行为。

由于横向应变值为负且绝对值很大 (ν_xy = -5,44)，在板的自由端完全过压，这里不会变形。

正交各向异性在 y 方向上的影响在这里接近零 (E_y ≈ 0)。

总结

利用 RFEM 中的正交各向材料模型几乎可以对任何材料参数进行定义。 基于横向应变的变量可以得出非常不同的结果。 按照 [1] 中进行数值修改过的横向应变得出一个接近单跨梁的解。

公式 31：

过高的负的横向应变表示一个修改的静力学系统，这不在与建模相关。