Material Models
The material models are the basis for composing multilayer surfaces to obtain an effective surface stiffness. The Multilayer Surfaces Add-on allows you to freely combine the material models in the RFEM 6 program. The basis of the material models is described in Chapters Materials and Nonlinear Material Behavior of the RFEM manual.
A selection of the possible combinations of the material models is created in the "Multilayer Models" model (see the right column), which you can download for further study of the combinations.
The following list shows a selection of the possible combinations:
- Isotropic layers (e.g. concrete - steel)
- Orthotropic layers (for example, cross-laminated timber)
- Isotropic - orthotropic (e.g. steel - GFRP)
- Isotropic Plastic - Isotropic (e.g. concrete - steel)
- Isotropic Nonlinear Elastic - Orthotropic (e.g. Concrete - Timber)
- Isotropic - Orthotropic Plastic (e.g. Concrete - Timber)
- Isotropic Damage - Orthotropic (e.g. Concrete - Timber)
Stiffnesses for multilayer surfaces without solids
The simpler calculation option in the Multilayer Surfaces add-on is to define different surface layers in the thickness type 'Layers' without solids. However, you can also freely combine the material models here.
Once the layers are defined, the Multilayer Surfaces add-on creates a global stiffness matrix of the surface. In RFEM, the internal forces and deformations are calculated for this surface. In the respective design add-on, such as Timber Design or Stress-Strain Analysis, these internal forces are then divided into the existing layers. Usually, the internal forces are displayed as three integration points per position.
This article explains how to calculate the stiffness matrix for isotropic and orthotropic materials.
stiffness matrix calculation
The material models are based on the following conditions (see Chapter Materials of the RFEM manual):
- All stiffness values ≥ 0
- Overall stiffness matrix of surface must be positive definite.
- Basic equation isotropic:
E | Modulus of elasticity |
G | Shear modulus |
ν | Poisson's ratio |
- Basic Equation Orthotropic:
Local stiffness matrix of each layer
- Isotropic
- Orthotropic
The shear stiffnesses for orthotropic material are as follows:
Gxy | Shear modulus in the panel plane (e.g. 690 N/mm² for C24) |
Gxz | Shear modulus in the direction x over the thickness (e.g. 690 N/mm² for C24) |
Gyz | Shear modulus in direction y over the thickness (e.g. 690 N/mm² for C24) – also called "rolling shear modulus". |
Furthermore, orthotropic material has the special feature that directional stiffnesses can be defined in a surface. In the default case, the local orientation of the surface or layer in the x-direction corresponds to the stiffness in the x-direction. However, since this can be freely defined by means of the angle β in the 'Layers' thickness type, it is necessary to transform the stiffnesses accordingly.
Summed element of each layer:
Bending and torsional elements [Nm]
The matrix elements for bending and torsion are given in the equations below.
If there is only one layer of the thickness type 'Layers', the calculation is based on the parameters described in the RFEM manual]].
For the shear (Element D44/55) different equations apply in the thickness type 'Layers'. They are described in the Shear in Slab Plane section.
Eccentricity terms [Nm/m]
For asymmetrical plates, eccentricity terms arise. An asymmetrical surface can be, for example, in a fire resistance design due to one-sided charring of a cross‑laminated timber plate. The matrix elements are as follows:
Plate Plane [N/m]
In the "Pane Wall" plane, the normal stiffnesses are represented in the plane of the glass pane. The shear force in the pane is calculated using element D88. The matrix elements are as follows:
Shear in slab plane [N/m]
To determine the shear stiffness for an orthotropic material, you have to rotate the stiffnesses according to their orientation to the local surface axis. This has to be done for each layer of the thickness type 'Layers'. In a simple layer structure with a 0° orientation of the cover layer and a 90° orientation of the underlying layer, there is a high shear stiffness, which has to be considered accordingly for the multilayer model. The following image (Source [1]) shows this using an example of a cross-laminated timber plate.
In the laminate theory, the shear stiffness of a layered structure is calculated by transforming all bending and shear components in the respective directions of each layer. You can find more information about this in the literature listed below.
By using the transformation of the stiffness shown in the figure, the stiffnesses are added up. This summation is also known as the "Grashoff integral."
In order to calculate the stiffness in the x- and y-directions, a centroid of stiffness is calculated for each structure of a multilayer surface.
Center of stiffness in y-direction:
In order to determine the orientation per position in the calculation of the shear stiffness, the stiffnesses are determined according to the following equations.
G stands for the shear stiffness of the layers in order to avoid mistakes for the elements of the stiffness matrix (D).
The shear stiffness of each layer can also be displayed in matrix form as follows:
The eccentric shear stiffness in the following equation would always be zero and thus irrelevant for the symmetrical structure of cross-laminated timber (0°/90°/0°) mentioned above. In the case of diagonally glued cross-laminated timber DLT ( Diagonal Laminated Timber ), for example, this eccentricity element is not zero and therefore plays an important role.
Further information can be found in [4] and in this YouTube video.
Calculation of shear stiffness
The shear stiffness is determined in the following steps:
- First, the angle of maximum stiffness is determined. The angle φ shows the modification of the local coordinate system x of the surface with respect to the oriented direction x''.
- All stiffnesses are rotated in the oriented direction x'' according to the equations presented above.
- The pane stiffness matrix of each position (3 x 3) is transformed from the local coordinate system x', y' to the rotated system x'', y". In addition to calculating the directed-stiffness of each individual layer, this is also done for the moduli of elasticity of each layer.
- The shear stiffness is calculated with the equations (Grashoff integral) described above. The shear stiffness is calculated using the individual parts.
- The calculated stiffnesses for the oriented direction of the entire structure are finally calculated back by means of the angular relations and shown as the original stiffnesses D44, D55, and D45 in the stiffness matrix.
Increasing shear stiffness
Since geometries with very narrow surface strips are also possible by modeling laminates as a surface, the shear stiffness must be increased accordingly when calculating such problematic geometries.
The following equation shows this for the X-direction
The length l in the equation above means the shortest length of a box that can be placed over the corresponding geometry.
In another model, which can be downloaded on the right, a narrow surface of 10 cm in width was compared with an identical surface of 20 cm in width.
The shear stiffness of the narrow surface is D44=15,253 kN/m compared to D44=5,970.8 kN/m of the wider surface. As a result, the deformation of the stiffer surface is smaller and the shear load is greater, despite the identical loading.
Stiffnesses for multilayer surfaces with integrated solids
In the future, it will also be possible to define solids together with surfaces in the Multilayer Surfaces add-on. In the case of this type, a surface is also exported to RFEM. Since the generation of stiffnesses and the decomposition of internal forces is more time-consuming, this is explained separately.