Verification Examples FEM-Design 16.

FEM-Design

Verification Examples
version 1.1
2016

1
Verification Examples FEM-Design 16.0

StruSoft AB
Visit the StruSoft website for company and FEM-Design information at
www.strusoft.com

Verification Examples
Copyright © 2016 by StruSoft, all rights reserved.

Trademarks
FEM-Design is a registered trademark of StruSoft.

2
Verification Examples FEM-Design 16.0

Contents
1. Linear static calculations ...........................................................................................................4
1.1 Beam with two point loading at one-third of its span............................................................4
1.2 Calculation of a circular plate with concentrated force at its center......................................6
1.3 A simply supported square plate with uniform load............................................................10
2. Second order analysis...............................................................................................................12
2.1 A column with vertical and horizontal loads........................................................................12
2.2 A plate with in-plane and out-of-plane loads.......................................................................15
3. Stability analysis.......................................................................................................................19
3.1 Flexural buckling analysis of a beam modell with different boundary conditions..............19
3.2 Buckling analysis of a plate with shell modell....................................................................21
3.3 Lateral torsional buckling of an I section with shell modell................................................24
3.4 Lateral torsional buckling of a cantilever with elongated rectangle section........................26
4. Calculation of eigenfrequencies with linear dynamic theory...................................................28
4.1 Continuous mass distribution on a cantilever column.........................................................28
4.2 Free vibration shapes of a clamped circular plate due to its self-weight.............................31
5. Seismic calculation...................................................................................................................33
5.1 Lateral force method with linear shape distribution on a cantilever....................................33
5.2 Lateral force method with fundamental mode shape distribution on a cantilever...............36
5.3 Modal analysis of a concrete frame building.......................................................................38
6. Calculation considering diaphragms.........................................................................................47
6.1. A simple calculation with diaphragms................................................................................47
6.2. The calculation of the shear center.....................................................................................47
7. Calculations considering nonlinear effects...............................................................................52
7.1 Uplift calculation.................................................................................................................52
7.1.1 A trusses with limited compression members................................................................52
7.1.2 A continuous beam with three supports ........................................................................58
7.2 Cracked section analysis by reinforced concrete elements..................................................69
7.2.1 Cracked deflection of a simply supported beam............................................................69
7.2.2 Cracked deflection of a cantilever beam........................................................................74
7.3 Nonlinear soil calculation....................................................................................................80
8. Cross section editor...................................................................................................................81
8.1 Calculation of a compound cross section.............................................................................81
9. Design calculations...................................................................................................................82
9.1 Required reinforcement calculation for a slab.....................................................................82
References...................................................................................................................................103
Notes...........................................................................................................................................104

In this verification handbook we highlighted the analytical results with green and the finite
element results with blue background for better comparison. The analytical closed formulas
are highlighted with a black frame.
If the finite element mesh is not mentioned during the example it means that the
automatically generated mesh was used.

3
Verification Examples FEM-Design 16.0

1. Linear static calculations

1.1 Beam with two point loading at one-third of its span

Fig. 1.1.1 left side shows the simple supported problem. The loads, the geometric and material
properties are as follows:

Force F = 150 kN
Length L=6m
Cross section Steel I beam HEA 300
The second moment of inertia in the relevant direction I1 = 1.8264ˑ10-4 m4
The shear correction factor in the relevant direction ρ2 = 0.21597
The area of the cross section A = 112.53 cm2
Young's modulus E = 210 GPa
Shear modulus G = 80.769 GPa

F F δF=1

L/3 L/3 L/3 L/2 L/2

-F -0.5
V δV

+F +0.5
M δM

M=+FL/3
+L/4

Figure 1.1.1 – The beam theory and the application of a virtual force

The deflection of the mid-span based on the hand calculation (based on virtaul force theorem
[1], see Fig. 1.1.1 right side also):

23 FL3
e= [
2M L 2 2 L 1 L 2 L L 1 L 1
+ +
EI 3 3 3 4 2 6 3 4 6 3 4 2
+
2F
ρ GA
L
0.5 =
3
+ ]
FL
648 EI 3 ρ GA [ ]
23 150⋅63 150⋅6
e= + =0.03151 m=31.51 mm
648 210000000⋅1.8264ˑ10 −4
3⋅0.21597⋅80769000⋅0.011253

4
Verification Examples FEM-Design 16.0

The first part of this equation comes from the bending deformation and the second part comes
from the consideration of the shear deformation as well, because FEM-Design uses Timoshenko
beam theory (see the Scientific Manual).

The deflection and the bending moment at the mid-span based on the linear static calculation
with three 2-noded beam elements (Fig. 1.1.2 and Fig. 1.1.3):

e FEM =31.51 mm and the bending moment M FEM =300 kNm

The theoretical solution in this case (three 2-noded beam elements) must be equal to the finite
element solution because with three beam elements the shape functions order coincides with the
order of the theoretical function of the deflection (the solution of the differtial equation).

Figure 1.1.2 – The finite element model

Figure 1.1.3 – The mid-span deflection [mm]

Therefore the difference between the results of the two calculations is zero.

5
Verification Examples FEM-Design 16.0

1.2 Calculation of a circular plate with concentrated force at its center

In this chapter a circular steel plate with a concentrated force at its center will be analyzed. First
of all the maximum deflection (translation) of the plate will be calculated at its center and then
the bending moments in the plate will be presented.

Two different boundary conditions will be applied at the edge of the plate. In the first case the
edge is clamped (Case I.) and in the second case is simply supported (Case II.), see Fig. 1.2.1.

P P

Figure 1.2.1 – Clamped (Case I.) and simply supported (Case II.) circular plate with concentrated force

The input parameters are as follows:

The concentrated force P = 10 kN

The thickness of the plate h = 0.05 m
The radius of the circular plate R=5m
The elastic modulus E = 210 GPa
The Poisson's ratio ν = 0.3

The ratio between the diameter and the thickness is 2R/h = 200. It means that based on the
geometry the shear deformation only have negligible effects on the maximum deflections. It is
important because FEM-Design uses the Mindlin plate theory (considering the shear
deformation, see Scientific Manual for more details), but in this case the solution of Kirchhoff's
plate theory and the finite element result must be close to each other based on the mentioned
ratio.

The analytical solution of Kirchhoff's plate theory is given in a closed form [2][3].

6
Verification Examples FEM-Design 16.0

Case I.:
For the clamped case the maximum deflection at the center is:
P R2
w cl =
16 π
( E h3
12(1−ν 2) )
The reaction force at the edge:
P
Qr =
2π R

And the bending moment in the radial direction at the edge:

P
M cl =
4π

With the given input parameters the results based on the analytical and the finite element
solutions (with the default finite element mesh size, see Fig. 1.2.2) are:

10⋅52
w cl = =0.002069 m=2.069 mm w clFEM =2.04 mm
16 π
(
210000000⋅0.053
12(1−0.32) )
10 kN kN
Qr = =0.318 QrFEM =0.318
2π 5 m m
P 10 kNm kNm
M cl = = =0.796 M clFEM =0.796
4π 4π m m

Figure 1.2.2 – The clamped (Case I.) and the simply supported (Case II.) plate with the default mesh

7
Verification Examples FEM-Design 16.0

Case II.:
For the simply supported case the maximum deflection in the center is:
P R2
w ss = ( 3+
1+ν )
ν
16 π
( E h3
12(1−ν 2 ) )
The reaction force at the edge:
P
Qr =
2π R

With the given input parameters the results based on the analytical and the finite element
solutions (with the default finite element mesh size, see Fig. 1.2.2) are:

10⋅52
w ss = ( 3+0.3
1+0.3 )
=0.005252 m=5.252 mm w ssFEM =5.00 mm
16 π
(210000000⋅0.053
12( 1−0.32 ) )
10 kN
Qr = =0.318 kN/m Q rFEM =0.318
2π 5 m

Figure 1.2.3 – The deflected shape of Case I. (clamped) and Case II. (simply supported) with the default mesh

Fig. 1.2.3 shows the two deflected shape in side view. The different boundary conditions are
obvious based on the two different displacement shape. The differences between the analytical
solutions and finite element solutions are less than 5% but the results could be more accurate if
the applied mesh is more dense than the default size.

Based on the analytical solution the bending moments in plates under concentrated loads are
infinite. It means that if more and more dense mesh will be applied the bending moment under
the concentrated load will be greater and greater. Thus the following diagram and table (Fig.
1.2.4 and Table 1.2.1) shows the convergence analysis of Case I. respect to the deflection and
bending moment. The deflection converges to the analytical solution (w cl = 2.07 mm) and the
bending moment converges to infinite.

8
Verification Examples FEM-Design 16.0

3,0 16
14
2,5
12
2,0

Bending moment [kNm/m]

10
Deflection [mm]

1,5 8
6
1,0
Deflection 4
0,5
Bending moment 2
0,0 0
0 5000 10000 15000 20000 25000 30000 35000
Number of elements
Figure 1.2.4 – Convergence analysis regarding to deflection and bending moment

Number of elements Deflection Bending moment

Average element size [m]
[pcs.] [mm] [kNm/m]
341 2,034 2,73 0,5
533 2,057 3,91 0,4
957 2,060 3,62 0,3
2035 2,068 4,33 0,2
7994 2,072 5,49 0,1
31719 2,073 6,40 0,05
31772 2,075 10,70 Local refinement 1
31812 2,076 14,30 Local refinement 2
Table 1.2.1 – The convergence analysis

9
Verification Examples FEM-Design 16.0

1.3 A simply supported square plate with uniform load

In this example a simply supported concrete square plate will be analyzed. The external load is a
uniform distributed load (see Fig. 1.3.1). We compare the maximum displacements and
maximum bending moments of the analytical solution of Kirchhoff's plate theory and finite
element results.
The input parameters are in this table:

The intensity of the uniform load p = 40 kN/m2

The thickness of the plate h = 0.25 m
The edge of the square plate a=5m
The elastic modulus E = 30 GPa
Poisson's ratio ν = 0.2

The ratio between the span and the thickness is a/h = 20. It means that based on the geometry
the shear deformation may have effects on the maximum deflection. It is important because
FEM-Design uses the Mindlin plate theory (considering the shear deformation, see Scientific
Manual for more details), therefore in this case the results of Kirchhoff's theory and the finite
element result could be different from each other due to the effect of shear deformations.

Figure 1.3.1 – The square plate with simply supported edges, uniform load and default mesh size

Based on Kirchhoff's plate theory [2][3] the maximum deflection is in the center of the simply
supported square plate and its intensity can be given with the following closed form:

4
pa
w max =0.00416

( )
3
Eh
2
12(1−ν )

10
Verification Examples FEM-Design 16.0

The maximum bending moment in the plate if the Poisson's ratio ν = 0.2:

2
M max =0.0469 p a

According to the input parameters and the analytical solutions the results of this problem are the
following:
The deflection at the center of the plate:
40⋅5 4
w max =0.00416 =0.002556 m=2.556 mm w maxFEM =2.632 mm

(
30000000⋅0.253
12(1−0.22 ) )
The bending moment at the center of the plate:
kNm kNm
M max =0.0469⋅40⋅5 2=46.9 M maxFEM =45.97
m m
Next to the analytical solutions the results of the FE calculations are also indicated (see Fig,
1.3.2 and 1.3.3). The difference is less than 3% and it also comes from the fact that FEM-Design
considers the shear deformation also (Mindlin plate theory).

Figure 1.3.2 – The deflected shape [mm] and the reaction forces [kN/m] with the default mesh

Figure 1.3.3 – The internal forces; Mx – My – Mxy [kNm/m]

11
Verification Examples FEM-Design 16.0

2. Second order analysis

2.1 A column with vertical and horizontal loads

We would like to analyze the following column (see Fig. 2.1.1) with second order theory. First
of all we make a hand calculation with third order theory according to Ref. [6] and [8] with
stability functions. After this step we compare the results with FEM-Design. In this moment we
need to consider that in FEM-Design second order analysis is implemented and the hand
calulation will be based on third order theory therefore the final results won't be exactly the
same. By the FEM-Design calulation we split the column into 3 bar elements thus the finite
element number of the bars was 3 for more precise results.
The input parameters:

Elastic modulus E = 30 GPa

Normal force P = 2468 kN
Horizontal load q = 10 kN/m
Cross section 0.2 m x 0.4 m (rectangle)
Second moment of inertia in the relevant direction I = 0.0002667 m4
Column length L=4m

q
L

Figure 2.1.1 – The column with vertical and horizontal loads

According to Ref. [6] and [8] first of all we need to calculate the following assistant quantities:

P P 2468
ρ= = 2 = 2 =0.500
PE
( )(
π EI
L2
π 30000000⋅0.0002667
42 )
12
Verification Examples FEM-Design 16.0

The constants based on this value for the appropriate stability functions:

s=3.294 ; c=0.666 ; f =1.104

With these values the bending moments and the shear forces based on third order theory and
FEM-Design calculation:

2 2
qL 10⋅4
M clamped = f (1+c) =1.104(1+0.666) =24.52 kNm M 2ndFEMclamped =25.55 kNm
12 12
M roller =0.0 kNm M 2ndFEMroller =0.0 kNm

[
V clamped = 1+
f (1+c )
6 ]( ) [
qL
2
= 1+
1.104(1+0.666)
6 ]( )
10⋅4
2
=26.13 kN

V 2ndFEMclamped =26.38 kN

[
V roller = 1−
f (1+c )
6 ]( ) [
qL
2
= 1−
1.104 (1+0.666)
6 ]( )
10⋅4
2
=13.87 kN

V 2ndFEMroller =13.61 kNm

The differences are less than 5 %.

Figure 2.1.2 – The shear [kN] and bending moment [kNm] diagram with 1st and 2nd order theory

13
Verification Examples FEM-Design 16.0

Figure 2.1.3 – The lateral translations with 1st and 2nd order theory

14
Verification Examples FEM-Design 16.0

2.2 A plate with in-plane and out-of-plane loads

In this chapter we will analyze a rectangular plate with single supported four edges. The load is
a specific normal force at the shorter edge and a lateral distributed total load perpendicular to the
plate (see Fig. 2.2.1). The displacement and the bending moment are the question based on a 2 nd
order analysis. First of all we calculate the results with analytical solution and then we compare
the results with FE calculations.
In this case the material and the geometric properties are the following:

The thickness of the plate h = 0.05 m

The dimensions of the plate a = 8 m; b = 6 m
The elastic modulus E = 210 GPa
Poisson's ratio ν = 0.3
The specific normal force Nx = 1000 kN/m
The lateral distributed load qz = 10 kN/m2

Nx y
qz

Nx
a
b

x
Figure 2.2.1 – The single supported edges, the lateral distributed load and the specific normal force

The maximum displacement and moments based on the 1st order linear calculation:

kNm kNm kNm

w max =35.38 mm , m x , max=18.05 , m y , max =25.62 , m xy , max =13.68
m m m

In Chapter 3.2 the critical specific normal force for this example is:

kN
N cr =2860
m

If the applied specific normal force is not so close to the critical value (now it is lower than the

15
Verification Examples FEM-Design 16.0

half of the critical value) we can assume the second order displacements and internal forces
based on the linear solutions with the following formulas (with blue highlight we indicated the
results of the FE calculation):

1 1
w max ,2nd =w max =35.38 =54.41 mm , w max ,2nd , FEM =54.69 mm

( 1−
Nx
N cr )1−
1000
2860 ( )
1 1 kNm kNm
m x , max ,2nd =mx , max =18.05 =27.76 , m x , max ,2nd , FEM =28.50
m m
( N
1− x
N cr ) 1−
1000
2860( )
1 1 kNm kNm
m y , max ,2 nd =m y , max =25.62 =39.40 , m y , max ,2 nd , FEM =40.30
m m
( N
1− x
N cr ) 1−
1000
2860( )
1 1 kNm kNm
m xy , max ,2nd =mxy , max =13.68 =21.04 m xy , max ,2nd , FEM =20.53
m m
( N
1− x
N cr ) 1−
1000
2860( )
The differences are less than 3 %.

Figure 2.2.2 shows the problem in FEM-Design with the default mesh.

Figure 2.2.2 – The single supported slab with in-plane and out-of-plane loads

The following figures show the moment distribution in the plate and the displacements with
FEM-Design according to 1st and 2nd order theory.

16
Verification Examples FEM-Design 16.0

Figure 2.2.3 – The mx [kNm/m] moment with 1st and 2nd order analysis

Figure 2.2.4 – The my [kNm/m] moment with 1st and 2nd order analysis

Figure 2.2.5 – The mxy [kNm/m] moment with 1st and 2nd order analysis

17
Verification Examples FEM-Design 16.0

Figure 2.2.6 – The vertical translation [mm] with 1st and 2nd order analysis

18
Verification Examples FEM-Design 16.0

3. Stability analysis

3.1 Flexural buckling analysis of a beam modell with different boundary

conditions

The cross section is a rectangular section see Fig. 3.1.1

The material C20/25 concrete
The elastic modulus E = 30000 GPa
The second moment of inertia about the weak axis I2 = 2.667ˑ10-4 m4
The length of the column L=4m
The boundary conditions see Fig. 3.1.1

b=0.4 m
1
L=4 m

a=0.2 m

Figure 3.1.1 – The buckling problem with the different boundary conditions and the cross section

The critical load parameters according to the Euler's theory are as follows and next to the
analytical solution [1] the relevant results of the FEM-Design calculation can be seen.

Pinned-pinned boundary condition:

π 2 E I2
F cr1= =4934.8 kN F crFEM1=4910.9 kN
L2

Fixed-pinned boundary condition:

π2EI2
F cr2= =10094.1 kN F crFEM2 =9974.6 kN
(0.6992 L)2

19
Verification Examples FEM-Design 16.0

Fixed-fixed boundary condition:

π 2E I2
F cr3= =19739.2 kN F crFEM3=19318.7 kN
(0.5 L)2

Fixed-free boundary condition:

π 2E I2
F cr4= =1233.7 kN F crFEM4 =1233.1 kN
(2 L)2

The difference between the two calculations are less than 3% but keep in mind that FEM-Design
considers the shear deformation therefore we can be sure that the Euler's results give higher
critical values. Fig. 3.1.2 shows the first mode shapes of the problem with the different
boundary conditions.

Figure 3.1.2 – The buckling mode shapes for different boundary conditions
pinned-pinned; fixed-pinned; fixed-fixed; fixed-free

20
Verification Examples FEM-Design 16.0

3.2 Buckling analysis of a plate with shell modell

In this chapter we will analyze a rectangular plate with single supported four edges. The load is
a specific normal force at the shorter edge (see Fig. 3.2.1). The critical force parameters are the
questions due to this edge load, therefore it is a stability problem of a plate.
In this case the material and the geometric properties are the following:

The thickness of the plate h = 0.05 m

The dimensions of the plate a = 8 m; b = 6 m
The elastic modulus E = 210 GPa
Poisson's ratio ν = 0.3

The solutions of the differential equation of the plate buckling problem are as follows [6]:

2π
2
( E h3
2 )
(
N cr =
m b n2 a
a
+
mb ) 12(1−ν )
b
2 , m=1, 2, 3 ... , n=1, 2, 3...

y
Nx

Nx
a
b

x
Figure 3.2.1 – The single supported edges and the specific normal force

Figure 3.2.2 shows the problem in FEM-Design with the default mesh.

21
Verification Examples FEM-Design 16.0

Figure 3.2.2 – The stability problem of a plate with single supported edges

According to the analytical solution the first five critical load parameters are:

m=1 , n=1

2π
2
( 210000000⋅0.053
12(1−0.32) )
N cr1=( 1⋅6 12 8
8
+
1⋅6 ) 62
=2860.36
kN
m
N crFEM1 =2862.58
kN
m
m=2 , n=1

2π
2
( 210000000⋅0.053
12(1−0.32) )
N cr2=(2⋅6 1 2 8
+
8 2⋅6 ) 6 2
=3093.77
kN
m
N crFEM2 =3109.96
kN
m

m=3 , n=1

2π
2
( 210000000⋅0.05 3
12(1−0.32 ) )
N cr3=( 3⋅6 12 8
8
+
3⋅6 ) 62
=4784.56
kN
m
N crFEM3 =4884.90
kN
m

m=4 , n=1

2 π2
( 210000000⋅0.05 3
12(1−0.32 ) )
N cr4 = ( 4⋅6 12 8
8
+
4⋅6 ) 6 2
=7322.53
kN
m
N crFEM4 =7655.58
kN
m

m=3 , n=2

2 π 2
( 210000000⋅0.053
2 )
(
N cr5 =
3⋅6 22 8
8
+
3⋅6 ) 12(1−0.3 )
6 2
=10691.41
kN
m
N crFEM5 =10804.62
kN
m

Next to these values we indicated the critical load parameters what were calculated with the
FEM-Design.

22
Verification Examples FEM-Design 16.0

The difference between the calculations less than 5 %.

Figure 3.2.3 shows the first five stability mode shapes of rectangular single supported plate.

Figure 3.2.3 – The first five stability mode shape of the described problem

23
Verification Examples FEM-Design 16.0

3.3 Lateral torsional buckling of an I section with shell modell

The purpose of this example is calculate the lateral torsional critical moment of the following
simple supported beam (see Fig. 3.3.1).

M M

L=10 m

Figure 3.3.1 – The static frame of a simple supported beam loaded with bending moments of both ends

The length of the beam is L = 10 m

The cross section see Fig. 3.3.2
The warping constant of the section Iw = 125841 cm6
The St. Venant torsional inertia It = 15.34 cm4
The minor axis second moment of area Iz = 602.7 cm4
The elastic modulus E = 210 GPa
The shear modulus G = 80.77 GPa

Figure 3.3.2 – The dimensions of the double symmetric cross section

In this case the critical moment can be calculated with the following formula based on the
analytical solution [6]:

π 2E Iz
√
2
I w L G It
M cr = +
L2 Iz π2E Iz

M cr =
10002 √
π 2⋅21000⋅602.7 125841 10002⋅8077⋅15.34
602.7
+ 2
π ⋅21000⋅602.7
=4328 kNcm

24
Verification Examples FEM-Design 16.0

In FEM-Design a shell modell was built to analyze this problem. The bending moments in the
shell model were considered with line loads at the end of the flanges (see Fig. 3.3.3).

The supports provides the simple supported beam effects with a fork support for the shell model
(see Fig. 3.3.3).

Figure 3.3.3 – The FEM model with the supports and the loads (moments) at the ends

From the FEM-Design stability calculation the critical moment value for this lateral torsional
buckling problem is:

M crFEM=4363 kNcm

The critical shape is in Fig 3.3.4. The finite element mesh size was provided based on the
automatic mesh generator of FEM-Design.

Figure 3.3.4 – The critical mode shape of the problem

The difference between the two calculated critical moments is less than 1%.

25
Verification Examples FEM-Design 16.0

3.4 Lateral torsional buckling of a cantilever with elongated rectangle section

The purpose of this example is calculate the critical force at the end of a cantilever beam (see
Fig. 3.4.1). If the load is increasing the state of the cantilever will be unstable due to lateral
torsional buckling.

Pcr

h
L

t
Figure 3.4.1 – The cantilever beam with concentrated load

The input parameters:

The length of the beam is L = 10 m

The cross section t = 40 mm; h = 438 mm; see Fig. 3.4.1
The St. Venant torsional inertia It = 8806246 mm4
The minor axis second moment of area Iz = 2336000 mm4
The elastic modulus E = 210 GPa
The shear modulus G = 80.77 GPa

In this case (elongated rectangle cross section with cantilever boundary condition) the critical
concentrated force at the end can be calculated with the following formula based on analytical
solution:

P cr =
L2 √
4.01 E I z G I t
E Iz

P cr =
4.01⋅210000⋅2336000
100002 √ 80770⋅8806246
210000⋅2336000
=23687 N=23.69 kN

26
Verification Examples FEM-Design 16.0

In FEM-Design a shell modell was built to analyze this problem. The concentrated load at the
end of the cantilever was considered at the top of the beam (see Fig. 3.4.2).

Figure 3.4.2 – The FE model of the cantilever beam with the default mesh

With the FEM-Design stability calculation the critical concentrated force value for this lateral
torsional buckling problem is:
P crFEM=24.00 kN
The critical shape is in Fig 3.4.3. The finite element mesh size was provided based on the
automatic mesh generator of FEM-Design.

Figure 3.4.3 – The critical mode shape of the problem

The difference between the two calculated critical load parameters is less than 2%.

27
Verification Examples FEM-Design 16.0

4. Calculation of eigenfrequencies with linear dynamic theory

4.1 Continuous mass distribution on a cantilever column

Column height H=4m

The cross section square with 0.4 m edge
The second moment of inertia I = 0.002133 m4
The area of the cross section A = 0.16 m2
The shear correction factor ρ = 5/6 = 0.8333
The elastic modulus E = 30 GPa
The shear modulus G = 12.5 GPa
The specific self-weight of the column γ = 25 kN/m3
The mass of the column m = 1.631 t
continuous mass distribution

H=4 m

Figure 4.1.1 – The cantilever with continuous mass distribution

Based on the analytical solution [4] the angular frequencies for this case is:

ω B =μ Bi
√ EI
mH3
; μ B1=3.52; μ B2=22.03; μ B3=61.7

if only the bending deformations are considered.

The angular frequencies are [4]:

28
Verification Examples FEM-Design 16.0

ω S = μ Si
√ ρ GA
mH
; μ S1=0.5 π ; μ S2=1.5 π ; μ S3=2.5 π

if only the shear deformations are considered.

Based on these two equations (considering bending and shear deformation) using the Föppl
theorem the angular frequency for a continuous mass distribution column is:

1 1 1
= +
ω 2n ω 2B ω 2S

Based on the given equations the first three angular frequencies separately for bending and shear
deformations are:

ω B1=3.52
√30000000⋅0.002133
16⋅4 3
=87.16
1
s

ω B2 =22.03
√ 30000000⋅0.002133
1.631⋅4 3
=545.4
1
s

ω B3=61.7
√30000000⋅0.002133
1.631⋅4 3
=1527.7
1
s

ω S1=0.5 π
√ 0.8333⋅12500000⋅0.16
1.631⋅4
=793.9
1
s

ω S2=1.5 π
√ 0.8333⋅12500000⋅0.16
1.631⋅4
=2381.8
1
s

ω S3=2.5 π
√ 0.8333⋅12500000⋅0.16
1.631⋅4
=3969.6
1
s
According to the Föppl theorem the resultant first three angular frequencies of the problem are:

1 1 1
ω n1=86.639 , ω n2=531.64 , ω n3=1425.8
s s s

And based on these results the first three eigenfrequencies are (f = ω/(2π)):

1 1 1
f n1=13.789 , f n2=84.613 , f n3=226.923
s s s

29
Verification Examples FEM-Design 16.0

In FEM-Design to consider the continuous mass distribution 200 beam elements were used for
the cantilever column. The first three planar mode shapes are as follows according to the FE
calculation:

1 1 1
f FEM1 =13.780 , f FEM2 =83.636 , f FEM3 =223.326
s s s

The first three mode shapes can be seen in Fig. 4.1.2.

Figure 4.1.2 – The first three mode shapes for the cantilever with continuous mass distribution

The differences between the analytical and FE solutions are less than 2 %.

30
Verification Examples FEM-Design 16.0

4.2 Free vibration shapes of a clamped circular plate due to its self-weight

In the next example we will analyze a circular clamped plate. The eigenfrequencies are the
question due to the self-weight of the slab.
In this case the material and the geometric properties are the following:

The thickness of the plate h = 0.05 m

The radius of the circular plate R=5m
The elastic modulus E = 210 GPa
Poisson's ratio ν = 0.3
The density ρ = 7.85 t/m3

The solution of the dynamic differential equation for the first two angular frequencies of a
clamped circular plate are [5]:

√ ( )
E h3
2 12 (1−ν 2 )
ω nm= π 2 β nm2 , β 10=1.015 , β 11=1.468
R ρh

Figure 4.2.1 shows the problem in FEM-Design with the clamped edges and with the default
mesh.

Figure 4.2.1 – The clamped circular plate and the default finite element mesh

According to the analytical solution the first two angular frequencies are:

31
Verification Examples FEM-Design 16.0

√( )
210000000⋅0.05 3
2 12(1−0.32 ) 1 1 1
ω 10= π 2 1.0152 =31.83 , f 10=5.066 , f 10FEM =5.129
5 7.85⋅0.05 s s s

√ ( )
210000000⋅0.053
2 12(1−0.3 2) 1 1 1
ω 11= π 2 1.468 2 =66.58 , f 11=10.60 , f 11FEM =10.731
5 7.85⋅0.05 s s s

Based on the angular frequencies we can calculate the eigenfrequencies in a very easy way. Next
to these values we indicated the eigenfrequencies what were calculated with the FEM-Design.
The difference between the calculations less than 2 %.
Figure 4.2.2 shows the first two vibration mode shapes of the circular clamped plate.

Figure 4.2.2 – The first two vibration shape mode of a clamped circular plate

32
Verification Examples FEM-Design 16.0

5. Seismic calculation

5.1 Lateral force method with linear shape distribution on a cantilever

Inputs:
Column height H = 10 m
The cross section square with 0.4 m edge
The second moment of inertia I = 0.002133 m4
The elastic modulus E = 31 GPa
The concentrated mass points 10 pieces of 1.0 t (see Fig. 5.1.1)
The total mass m = 10.0 t

Figure 5.1.1 – The cantilever column with the concentrated mass points, the first vibration shape [T=0.765 s],
the equivalent forces [kN], the shear force diagram [kN] and the bending moment diagram [kNm] with FEM-
Design

First of all based on a hand calculation we determine the first fundamental period:
The first fundamental period of a cantilever column (length H) with a concentrated mass at the
end (m mass) and EI bending stiffness [4]:

2π
T i=

√ 3 EI
mi H i3

33
Verification Examples FEM-Design 16.0

The fundamental period separately for the mass points from bottom to top:

2π 2π
T 1= =0.01411s ; T 2 = =0.03990 s ;

√ 3⋅31000000⋅0.0021333
1⋅13
2π
√
3⋅31000000⋅0.0021333
1⋅2 3
2π
T 3= =0.07330 s ; T 4 = =0.1129 s ;

√ 3⋅31000000⋅0.0021333
1⋅33
2π
√
3⋅31000000⋅0.0021333
1⋅4 3
2π
T 5= =0.1577 s ; T 6 = =0.2073s ;

√ 3⋅31000000⋅0.0021333
1⋅53
2π
√
3⋅31000000⋅0.0021333
1⋅6 3
2π
T 7= =0.2613s ; T 8= =0.3192 s ;

√ 3⋅31000000⋅0.0021333
1⋅7
2π
3
√3⋅31000000⋅0.0021333

2π
1⋅8
3

T 9= =0.3809 s ; T 10= =0.4461 s .

√ 3⋅31000000⋅0.0021333
1⋅9 3 √
3⋅31000000⋅0.0021333
1⋅103

The approximated period based on these values according to the Dunkerley summary and the
result of FE calculation:

√∑
10
T HC = T i2 =0.7758 s T FEM =0.765s
i=1

The difference between the hand calculation and FEM-Design calculation is less than 2%, for
further information on the period calculation see Chapter 4.

The base shear force according to the fundamental period of vibration (see Fig. 5.1.1) and the
response spectrum (see Fig. 5.1.2):

F b=S d (T 1)m λ =0.6588⋅10⋅1.0=6.588 kN

We considered the response acceleration based on the period from FE calculation to get a more
comparable results at the end. Thus the equivalent forces on the different point masses are:
z i mi z i mi zm
F i =F b =6.588 =6.588 i i
∑ z j mj 1⋅1+2⋅1+3⋅1+4⋅1+5⋅1+6⋅1+7⋅1+8⋅1+9⋅1+10⋅1 55

34
Verification Examples FEM-Design 16.0

The equivalent forces from the bottom to the top on each point mass:

F 1 =0.120 kN ; F 2=0.240 kN ; F 3=0.359 kN ; F 4=0.479 kN ; F 5=0.599 kN ;

F 6 =0.719 kN ; F 7 =0.838 kN ; F 8=0.958 kN ; F 9=1.078 kN ; F 10 =1.198 kN

These forces are identical with the FEM-Design calculation and the shear force and bending
moment diagrams are also identical with the hand calculation.

Figure 5.1.2 – The response spectrum [T = 0.765 s; Sd = 0.6588 m/s2]

35
Verification Examples FEM-Design 16.0

5.2 Lateral force method with fundamental mode shape distribution on a

cantilever
Inputs:
Column height H = 10 m
The cross section square with 0.4 m edge
The second moment of inertia I = 0.002133 m4
The elastic modulus E = 31 GPa
The concentrated mass points 10 pieces of 1.0 t (see Fig. 5.1.1)
The total mass m = 10.0 t

Figure 5.2.1 – The cantilever column with the concentrated mass points, the first vibration shape with
the value of the eigenvector [T=0.765 s], the equivalent forces [kN], the shear force diagram [kN] and
the bending moment diagram [kNm] with FEM-Design

The base shear force according to the fundamental period of vibration (see Fig. 5.2.1) and the
response spectrum (see Fig. 5.2.2):

F b=S d (T 1)m λ =0.6588⋅10⋅1.0=6.588 kN

We considered the response acceleration based on the period from FE calculation to get a more
comparable results at the end. Thus the equivalent forces on the different point masses are:

36
Verification Examples FEM-Design 16.0

s i mi
F i =F b =
∑ sjmj
si mi sm
=6.588 =6.588 i i
9.5⋅1+36.2⋅1+77.4⋅1+130.6⋅1+193.4⋅1+ 263.3⋅1+338.4⋅1+416.8⋅1+496.9⋅1+577.6⋅1 2540.1
The equivalent forces from the bottom to the top on each point mass:

F 1 =0.0246 kN ; F 2=0.0939 kN ; F 3=0.201 kN ; F 4=0.339 kN ; F 5=0.502 kN ;

F 6 =0.683 kN ; F 7 =0.878 kN ; F 8=1.081 kN ; F 9=1.289 kN ; F 10 =1.498 kN

These forces are identical with the FEM-Design calculation and the shear force and bending
moment diagrams are also identical with the hand calculation.

Figure 5.2.2 – The response spectrum [T = 0.765 s; Sd = 0.6588 m/s2]

37
Verification Examples FEM-Design 16.0

5.3 Modal analysis of a concrete frame building

In this chapter we show a worked example for modal analysis on a concrete frame building
according to EN 1998-1:2008 with hand calculation and compare the results with FEM-Design.
This example is partly based on [4]. The geometry, the dimensions, the material and the bracing
system are in Fig. 5.3.1-3 and in the following table.
Inputs:
Column height/Total height h = 3.2 m; H = 2·3.2=6.4 m
The cross sections Columns: 30/30 cm; Beams: 30/50 cm
The second moment of inertia Ic = 0.000675 m4; Ib = 0.003125 m4
The elastic modulus E = 28.80 GPa
The concentrated mass points 12 pieces of 13.358 t on 1st storey and
12 pieces of 11.268 t on 2nd storey
(see Fig. 5.3.2)
The total mass 1st storey: m1 = 160.3 t
2nd storey: m2 = 135.2 t
total mass: M = 295.5 t
Reduction factor for elastic modulus α = 0.5
considering the cracking according to EN 1998-
1:2008
Behaviour factors q = 1.5, qd = 1.5
Accidental torsional effect do not considered ξ = 0.05 (damping factor)

Figure 5.3.1 – The concrete frame building with the columns and beams

38
Verification Examples FEM-Design 16.0

The first exercise is the determination of the fundamental periods and mode shapes. There are
several hand calculation modes to get these values but in this chapter the details of the modal
analysis are important therefore we considered the first two fundamental periods based on FEM-
Design calculation (see Fig. 5.3.5). See the details and example on the eigenfrequency
calculation in Chapter 4.
The dead loads and the live loads are considered in the mass points (see Fig. 5.3.2).

Figure 5.3.2 – The frame building with the masses and bracings

Figure 5.3.3 – The side view of the building

39
Verification Examples FEM-Design 16.0

Figure 5.3.4 – The considered design response specra according to EN 1998-1:2008

Figure 5.3.5 – The first two fundamental mode shapes [-], T1 = 0.704 s; T2 = 0.252 s

According to the fundamental periods in Fig. 5.3.5 the response accelerations from Fig. 5.3.4
are:

m
T 1=0.704 s S d1=1.115 ;
s2
m
T 2 =0.252 s S d2=1.57 2
.
s

The second step is to calculate the effective modal masses based on this formula:

2
*
m =
( Φi T m ι )
i
Φ iT m Φ i

40
Verification Examples FEM-Design 16.0

During the hand calculation we assume that the structure is a two degrees of freedom system in
the x direction with the two storeys, because the first two modal shapes are in the same plane
see Fig. 5.3.5. Thus we only consider the seismic loads in one direction because in this way the
hand calculation is more comprehensible.

m*1 =
( 0 [
[40.9 73.6 ] 160.3 0 1
135.2 1 ][ ]) =272.3 t ;
m *1 272.3
= =92.1 %
M 295.5
[40.9 73.6] 160.3
0 [ 0 40.9
135.2 73.6 ][ ]
2

m*2 =
([67.6 −44.5] 160.3
0 [
0 1
135.2 1 ][ ]) =23.23 t ;
m *2 23.23
= =7.9 %
M 295.5
[67.6 −44.5] 160.3
0 [
0 67.6
135.2 −44.5 ][ ]
According to the assumption of a two degrees of freedom system the sum of the effective modal
masses is equal to the total mass:

m*1 m*2 272.3 23.23

+ = + =100.0 %
M M 295.5 295.5

Calculation of the base shear forces:

*
F b1=S d1 m1 =1.115⋅272.3=303.6 kN ; F b2=S d2 m*1=1.570⋅23.23=36.5 kN

The equivalent forces come from this formula:

Φ iT m ι
p i=m Φi S di
Φ iT m Φ i

The equivalent forces at the storeys respect to the mode shapes considering the mentioned two
degrees of freedom model:

0[
[40.9 73.6 ] 160.3 0 1
135.2 1 ][ ]
p1 = [ 160.3 0 40.9
][ ] 1.115= 120.6 kN [ ]
0 135.2 73.6
[40.9 73.6] [
160.3
0
0 40.9
135.2 73.6 ][ ] 183.0

41
Verification Examples FEM-Design 16.0

[67.6 −44.5] 160.3

0
0
[1
135.2 1 ][ ]
p2 = [ 160.3 0
][ 67.6
] [67.6 1.570= 81.98 kN [ ]
0 135.2 −44.5
[
−44.5] 160.3
0
0 67.6
135.2 −44.5 ][
−44.52
]
The equivalent forces on one frame from the six (see Fig. 5.3.1):

[ ][ ]
pf1 = 120.6 /6 = 20.10 kN
183.0 /6 30.50

p =[ 81.98/6 ]=[ 13.66 ] kN

f2
−44.52 /6 −7.420

The shear forces between the storeys respect to the two different mode shapes:

[ ][ ]
V 1= 20.1+30.5 = 50.6 kN
30.5 30.5 [
V 2= 13.66−7.42 = 6.24 kN
−7.42 −7.42 ][ ]
The shear forces in the columns respect to the two different mode shapes:

[ ][ ]
V c1= 50.6 /2 = 25.30 kN
30.5/2 15.25 [
V c2= 6.24 /2 = 3.13 kN
−7.42 /2 −3.71 ][ ]
The bending moments in the columns respect to the two different mode shapes from the relevant
shear forces (by the hand calculation we assumed zero bending moment points in the middle of
the columns):

[ ][ ]
M c1= 25.30⋅3.2/2 = 40.48 kNm
15.25⋅3.2 /2 24.40 [
M c2= 3.13⋅3.2/2 = 5.008 kNm
−3.71⋅3.2/2 −5.936 ][ ]
The bending moments in the beams respect to the two different mode shapes:

[ ][ ]
M b1 = 40.48+24.40 = 64.88 kNm
24.40 24.40 [ ][
M b2 = 5.008−5.936 = −0.928 kNm
−5.936 −5.936 ]
The SRSS summation on the internal forces:

42
Verification Examples FEM-Design 16.0

V c=
[ √ 25.30 2+3.132
√ 15.25 +(−3.71) 15.69
2 2 ][ ]
= 25.49 kN M c=
[ ][ ]
√ 40.482+5.0082 =
√ 24.40 2+5.936 2
40.79 kNm
25.11

M b=
[√√ 24.40 2+(−5.936)2 ][ ]
64.88 2+(−0.928)2 = 64.89 kNm
25.11

The CQC summation on the internal forces:

T 2 0.252
α 12= = =0.358
T 1 0.704
2 3 /2
8 ξ (1+α 12) α 12 8⋅0.05 2 (1+0.358)0.3583/ 2
r 12= 2
= 2
=0.007588
(1−α 122 ) +4 ξ 2 α 12 (1+α 12)2 (1−0.358 2) +4⋅0.052⋅0.358 (1+0.358)2
r= [ 1
0.007588
0.007588
1 ]
And based on these values the results of the CQC summation:

V c=
[√ √ 25.30 2+3.132 +2⋅25.3⋅3.13⋅0.007588
2 2
15.25 +(−3.71) +2⋅15.25⋅(−3.71)⋅0.007588
= 25.52 kN
15.67 ][ ]
M c=
[√√40.482+5.008 2+2⋅40.48⋅5.008⋅0.007588 = 40.83 kNm
24.40 2+5.936 2+2⋅24.40⋅5.936⋅0.007588 25.16 ][ ]
M b=
[√
√ 64.882+(−0.928)2+2⋅64.88⋅(−0.928)⋅0.007588 = 64.88 kNm
2 2
24.40 +(−5.936) +2⋅24.40⋅(−5.936⋅0.007588) 25.07 ][ ]
The following displacements come from the FEM-Design calculation on the complete frame
structure to ensure the comprehensible final results on the P-Δ effect.
The displacements at the storeys respect to the two different mode shapes considering the
displacement behaviour factor:

[ ] [ ][ ]
u1=q d 9.54 =1.5 9.54 = 14.31 mm
17.15 17.15 25.73 [ ] [ ][
u2 =q d 0.818 =1.5 0.818 = 1.227 mm
−0.540 −0.540 −0.810 ]
Based on these values the storey drifting respect to the two different mode shapes:

Δ 1= [ 25.73−14.31
14.31
]=[14.31
11.42 ]
mm [−0.810−1.227
Δ 2= 1.227
]=[−2.037 ]
1.227 mm

43
Verification Examples FEM-Design 16.0

SRSS summation on the story drifting:

Δ=
[√ √14.312+1.227 2
2
11.42 +(−2.027) 11.60
2 ][ ]
= 14.36 mm

P-Δ effect checking on the total building:

[ ][
P tot= (m1+m2 )g = (160.3+135.2)9.81 = 2899 kN
m2 g 135.2⋅9.81 1326 ][ ]
V tot =
[ 6⋅√ 30.5 2+(−7.42)2 ][ ]
6⋅√ 50.62+6.24 2 = 305.9 kN
188.3

P tot1 Δ1 2899⋅14.36 P tot2 Δ 2 1326⋅11.60

θ 1= = =0.0425 θ 2= = =0.0255
V tot1 h 305.9⋅3200 V tot2 h 188.3⋅3200

After the hand calculation let's see the results from the FEM-Design calculation and compare
them to each other. Fig. 5.3.6 shows the effective modal masses from the FE calculation.
Practically these values coincide with the hand calculation.

Figure 5.3.6 – The first two fundamental periods and the effective
modal masses from FEM-Design

Fig. 5.3.7 and the following table shows the equivalent resultant shear forces and the base shear
forces respect to the first two mode shapes. The differences between the two calculations are
less than 2 %.

Storey 1 equivalent Storey 2 equivalent Base shear force

resultant [kN] resultant [kN] [kN]
Hand FEM Hand FEM Hand FEM
Mode shape 1 120.6 121.9 81.98 81.80 303.6 306.9
Mode shape 2 183.0 185.0 – 44.52 – 45.47 36.50 36.33

44
Verification Examples FEM-Design 16.0

Figure 5.3.7 – The equivalent forces respect to the storeys and the base shear forces for the first two mode
shapes [kN]

Fig. 5.3.8-9 and the following table shows the internal forces after the differenr summation
modes (SRSS and CQC). The differences between the two calculations are less than 2 %.

Column shear force Column bending moment Beam bending moment

[kN] [kNm] [kNm]
Storey 1 Storey 2 Storey 1 Storey 2 Storey 1 Storey 2
SRSS
25.49 15.69 40.79 25.11 64.89 25.11
Hand
SRSS (37.18+45.32)/2=
25.78 15.89 27.51 59.33 27.51
FEM 41.25
CQC
25.50 15.67 40.83 25.16 64.88 25.07
Hand
CQC (37.21+45.36)/2=
25.80 15.86 27.46 59.32 27.46
FEM 41.29

Figure 5.3.8 – The shear force [kN] and bending moment diagram [kNm] after the SRSS summation rule

45
Verification Examples FEM-Design 16.0

Figure 5.3.9 – The shear force [kN] and bending moment diagram [kNm] after the CQC summation rule

Fig. 5.3.10 shows the Θ values from FEM-Design. The differences between the hand calculation
and FEM-Design are less than 3 %.

Figure 5.3.10 – The θ values at the different storeys from FEM-Design

46
Verification Examples FEM-Design 16.0

6. Calculation considering diaphragms

6.1. A simple calculation with diaphragms

If we apply two diaphragms on the two storeys of the building from Chapter 5.3 then the
eigenfrequencies and the periods will be the same what we indicated in Chapter 5.3.

6.2. The calculation of the shear center

In this example we show that how can we calculate the shear center of a storey based on the
FEM-Design calculation. We analyzed a bottom fixed cantilever structure made of three
concrete shear walls which are connected to each other at the edges (see Fig. 6.2.1). The
diaphragm is applied at the top plane of the structure (see also Fig. 6.2.1 right side). If the height
of the structure is high enough then the shear center will be on the same geometry point where it
should be when we consider the complete cross section of the shear walls as a “thin-walled” “C”
cross section (see Fig. 6.2.1 left side). Therefore we calculate by hand the shear center of the
“C” profile assumed to be a thin-walled cross section then compare the solution what we can get
from FEM-Design calculation with diaphragms.
Secondly we calculate the idealized bending stiffnesses in the principal rigidity directions by
hand and compare the results what we can calculate with FEM-Design results.
Inputs:
Height of the walls H = 63 m
The thickness of the walls t = 20 cm
The width of wall number 1 and 3 w1 = w3 = 4.0 m
The width of wall number 2 w2 = 6.0 m
The applied Young's modulus of concrete E = 9.396 GPa

Figure 6.2.1 – The geometry of the bracing core and the height of the bottom fixed structure
(the diaphragm is lying on the top plane, see the red line and hatch)

47
Verification Examples FEM-Design 16.0

First of all lets see Fig. 6.2.2. The applied cross section is a symmetric cross section. In the web
the shear stress distribution comes from the shear formula regarding bending (see Fig. 6.2.2).
Therefore it is a second order polynom. In the flanges the shear stress distribution is linear
according to the thin walled theory. With the resultant of these shear stress ditribution (see Fig.
6.2.2, V1, V2 and V3) the position of the shear center can be calculated based on the statical
(equilibrium) equations.

τ
w1
τ
V t V1

V2 G
τmax

w2
xS xG
V3
τ w1

τ
Figure 6.2.2 – The shear stress distribution in a thin-walled cross section if the
shear force acting on the shear center

The shear stress values (see Fig. 6.2.2):

VS 1⋅( 0.2⋅4⋅3) kN
τ= = =0.6665 2
It
( )
3 3 3
0.2⋅6 4⋅6.2 4⋅5.8 m
+ − ⋅0.2
12 12 12
V S max 1⋅(0.2⋅4⋅3+0.2⋅3⋅1.5) kN
τ max = = =0.9164 2
It
( )
3 3 3
0.2⋅6 4⋅6.2 4⋅5.8 m
+ − ⋅0.2
12 12 12

Based on these stresses the resultant in the flanges and in the web:
τ t w 1 0.6665⋅0.2⋅4
V 1=V 3= = =0.2667 kN
2 2
2 2
V 2 = (τ max −τ ) w 2 t+τ w2 t= (0.9164−0.6665) 6⋅0.2+0.6665⋅6⋅0.2=0.9997 kN
3 3

Respect to the equilibrium (sum of the forces):

V =1 kN≈V 2=0.9997 kN

48
Verification Examples FEM-Design 16.0

And also respect to the equilibrium (if the external load is acting on the shear center, see Fig.
6.2.2) the sum of the moments:
V 1 w 2=V 2 x S
V 1 w 2 0.2667⋅6
xS= = =1.601 m
V2 0.9997

Thus the shear center is lying on the symmetry axis and it is x S=1.601 m from the web (see Fig.
6.2.2). In FEM-Design the global coordinate system does not coincide with the symmetry axis
of the structure (see Fig. 6.2.1). Therefore we need no transform the results.
Lets be a selected key node at the diaphragm in the global coordinate system (see Fig. 6.2.1):
x m =0 m ; y m=0 m

Based on the unit forces (1 kN) and moment (1 kNm) on the key node the displacements of the
key node are as follows based on the FEM-Design calculation (see the theory manual
Calculation considering diaphragm chapter also):
According to unit force on key node in X direction:
u xx =1.5852 mm u yx =0.72166 mm φ zx =0.29744⋅10−4 rad

According to unit force on key node in Y direction:

u xy =0.72166 mm u yy =7.3314 mm φ zy =0.10328⋅10−2 rad

According to unit moment on key node around Z direction:

φ zz =0.16283⋅10−3 rad

Based on these finite element results the global coordinates of the shear center of the diaphragm
are:

φ 0.10328⋅10
−2
x S = x m−φ zyzz =0− =−6.343 m
0.16283⋅10−3

φ zx 0.29744⋅10−4
y S = y m +φ zz =0+ =+0.1827 m
0.16283⋅10−3

In FEM-Design the coordinates of the middle point of the web are (see Fig. 6.2.1):
x mid =−4.919 m ; y mid =+0.894 m

49
Verification Examples FEM-Design 16.0

With the distance between these two points we get a comparable solution with the hand
calculation.

√( x −x ) +( y S − y mid ) = √(−6.343−(−4.919)) +( 0.1827−0.894 ) =1.592 m

2 2 2 2
x SFEM = S mid

The difference between FEM and hand calculation is less than 1%.

The gravity center of the cross section (Fig. 6.2.2) can be calculated based on the statical
moments. And of course the gravity center lying on the symmetry axis. The distance of the
gravity center from the web is:
Sy ' 2(0.2⋅4⋅2)
x ' G= = =1.143 m
A 0.2(4+6+4)

With the input Young's modulus and with the second moments of inertia the idealized bending
stiffnesses in the principal directions can be calculated by hand.

( )
3 3 3
0.2⋅6 4⋅6.2 4⋅5.8
E I 1=9396⋅103⋅ + − =1.692⋅108 kNm2
12 12 12

( )
3 3
0.2⋅4 6⋅0.2
E I 2 =9396⋅103⋅ 2 +2(0.2⋅4( 2−1.143)2 )+ +0.2⋅6(1.143)2 =4.585⋅107 kNm 2
12 12

With the finite element results we can calculate the translations of the shear center according to
the unit forces and moment on the key node (see the former calculation method).
The distances between the shear center and the selected key node are:

Δ x=x S− x m=−6.343−0=−6.343 m=−6343 mm

Δ y= y S − y m=+0.1827−0=+0.1827 m=182.7 mm

The translations of the shear center are as follows:

u Sxx =u xx−φ zx Δ y=1.5852−0.29744⋅10−4⋅182.7=1.5798 mm

−4
u Syx =u yx+φ zx Δ x=0.72166+0.29744⋅10 ⋅(−6343)=0.5330 mm

u Sxy =u xy−φ zy Δ y=0.72166−0.10328⋅10−2⋅182.7=0.5330 mm

−2
u Syy =u yy+φ zy Δ x=7.3314+0.10328⋅10 ⋅(−6343)=0.7803 mm

50
Verification Examples FEM-Design 16.0

Based on these values the translations of the shear center in the principal directions:

√( √(
2

)
2
u +u u Sxx −u Syy
u 1= Sxx Syy +
2 2
2
+u Sxy =
1.5798+0.7803
2
+
2 )
1.5798−0.7803 2
+0.5330 =

=1.8463mm

√( √(
2

)
2
u +u uSxx −u Syy
u 2= Sxx Syy −
2 2
2 1.5798+0.7803
+u Sxy =
2
−
2 )
1.5798−0.7803 2
+0.5330 =

=0.5138 mm

According to these values the angles of the principal rigidity directions:

u 1−u Sxx 1.8463−1.5798
α 1FEM =arctan =arctan =26.57o
u Sxy 0.533
u 2−u Sxx 0.5138−1.5798
α 2FEM =arctan =arctan =−63.43o
u Sxy 0.533

The directions coincide with the axes of symmetries (see Fig. 6.2.1-2) which is one of the
principal rigidity direction in this case.

Then with FEM-Design results we can calculate the idealized bending stiffnesses of the
structure:
H3 633
EI 1FEM = = E I 1FEM =1.622⋅108 kNm 2
3 u 2 3⋅( 0.5138/1000)
3 3
H 63 7 2
EI 2FEM = = E I 2FEM =4.514⋅10 kNm
3 u1 3⋅(1.8463/1000)
The difference between FEM and hand calculation is less than 4%.

51
Verification Examples FEM-Design 16.0

7. Calculations considering nonlinear effects

7.1 Uplift calculation

7.1.1 A trusses with limited compression members

In this example a truss will be analyzed. First of all we calculate the normal forces in the truss
members and the maximum deflection for the gived concentrated loads. After this step we
calculate the load multiplier when the vertical truss members reaches its limit compression
bearing capacity what we set. See the inputs in the following table. After the hand calulation we
compare the results with the FEM-Design nonlinear calculation results.
Inputs:
Column height/Span H = 2.0 m; L = 8.0 m
The cross sections KKR 80x80x6
The area of the cross sections A = 1652 mm2
The elastic modulus E = 210 GPa, steel
The concentrated loads F = 40 kN
Limited compression of the vertical truss members Pcr = 700 kN

F F F F F
4 8 12 16
H=2 m

3 7 11 13 15 17
5 9
1

2 6 10 14
L=8 m

δF=1

δF
δF=1

Figure 7.1.1.1 – The truss with the concentrated loads and with the
virtual loads for the translation claculation

52
Verification Examples FEM-Design 16.0

The normal forces in the truss members based on the hand calculation (without further details)
are:
N 1=N 17 =−100 kN ; N 2 =N 14=0 kN ; N 3= N 15 =+84.85 kN ;
N 4 =N 5= N 13= N 16=−60.00 kN ; N 6 =N 10=+60.00 kN ;
N 7 =N 11=+28.28 kN ; N 8 =N 12 =−80.00 kN ; N 9 =−40.00 kN .
The normal forces in the truss members according to the vertical virtual force (see Fig. 7.1.1.1):

N 1,1= N 1,17=−0.5 kN ; N 1,2= N 1,14=N1,9 =0 kN ;

N 1,3= N 1,15=N 1,7 =N 1,11=+0.7071 kN ;
N 1,4 =N 1,5= N 1,13=N 1,16 =−0.5 kN ; N 1,6 =N 1,10 =+0.5 kN ;
N 1,8 =N 1,12=−1.0 kN .
The normal forces in the truss members according to the horizontal virtual force (see Fig.
7.1.1.1):

N 2,2 =N 2,6 = N 2,10= N 2,14=+1.0 kN ;

N 2,1 =N 2,3= N 2,4 =N 2,5 =N 2,7 = N 2,8= N 2,9 =N 2,11= N 2,12=N 2,13=N 2,15 =N 2,16 =N 2,17 =0 kN .
The hand calculation of the vertical translation at the mid-span with the virtual force method:
17
1
e z= ∑ N δ N 1,i li =0.003841 m=3.841 mm
EA i=1 i
The hand calculation of the horizontal translation at right roller with the virtual force method:
17
1
e x= ∑ N δ N 2,i li =0.0006918 m=0.6918 mm
EA i=1 i

Figure 7.1.1.2 – The truss with the concentrated loads in FEM-Design

53
Verification Examples FEM-Design 16.0

Figure 7.1.1.3 – The reaction forces

Figure 7.1.1.4 – The normal forces in the truss members

Figure 7.1.1.5 – The vertical translation at the mid-span and the horizontal translation at the right roller [mm]

54
Verification Examples FEM-Design 16.0

The translations and the normal forces in the truss members based on the hand calculation are
identical with the FEM-Design calculation, see Fig. 7.1.1.2-5.
After this step we would like to know the maximum load multiplier when the vertical truss
members reaches its limit compression bearing capacity what we set, Pcr = 700 kN. The
maximum compression force arises in the side columns, see the hand calculation, N 1 = (–)100
kN. Therefore the load multiplier based on the hand calcualtion is λ = 7.0.
Let's see the FEM-Design uplift calculation considering the limit compression in the vertical
members.

Figure 7.1.1.6 – Large nodal displacements when the side truss

members reached the limit compression value [mm]

Figure 7.1.1.7 – The two different analyzed load multiplier in FEM-Design

With λFEM = 7.00 multiplier the FEM-Design analysis gives the accurate result but with λFEM =
7.01 (see Fig. 7.1.1.7) large nodal displacements occurred, see Fig. 7.1.1.6. Thus by this
structure if we neglect the effect of the side members the complete truss became a statically
over-determinated structure. FEM-Design solve this problem with iterative solver due to the fact
that these kind of problems are nonlinear.

55
Verification Examples FEM-Design 16.0

7.1.2 A continuous beam with three supports

In this example we analyse non-linear supports of a beam. Let's consider a continuous beam
with three supports with the following parameters:
Inputs:
Span length L = 2 m, total length = 2x2 = 4 m
The cross sections Rectangle: 120x150 mm
The elastic modulus E = 30 GPa, concrete C20/25
Intensity of distributed load (total, partial) p = 10 kN/m

In Case I. the distribution of the external load and the nonlinearity of the supports differ from
Case II. See the further details below (Fig. 7.1.2.1 and Fig. 7.1.2.8).
a) Case I.
In this case the distributed load is a total load (Fig. 7.1.2.1). In the first part all of the three
supports behave the same way for compression and tension. In the second part of this case the
middle support only bears tension. We calculate in both cases the deflections, shear forces and
bending moments by hand and compared the results with FEM-Design uplift (nonlinear)
calculations.

Figure 7.1.2.1 – The beam with three supports and uniform distributed load

In first part of this case the maximum deflection comes from the following formula considering
only the bending deformations in the beam:

4 4
2.1 p L 2.1 10⋅2
e max = = =0.0008642 m =0.8642 mm
384 EI 384 30000000⋅0.12⋅0.153 /12

The relevant results with FEM-Design:

Figure 7.1.2.2 – The deflection of the beam with three supports (total load)

56
Verification Examples FEM-Design 16.0

The extremums of the shear force without signs:

3 3 5 5
V 1= p L= 10⋅2=7.5 kN ; V 2 = p L= 10⋅2=12.5 kN
8 8 8 8

The relevant results with FEM-Design:

Figure 7.1.2.3 – The shear diagram of the beam with three supports (total load)

The extremums of the bending moment without signs:

9 2 9 2 1 2 1 2
M midspan = pL= 10⋅2 =2.812 kNm ; M middle = p L = 10⋅2 =5.0 kNm
128 128 8 8

The relevant results with FEM-Design:

Figure 7.1.2.4 – The bending moment diagram of the beam with three supports (total load)

When the middle support only bear tension (second part of this case) basically under the total
vertical load (Fig. 7.1.2.1) the middle support is not active (support nonlinearity). Therefore it
works as a simply supported beam with two supports. The deflection, the shear forces and the
bending moments are the following:
The maximum deflection comes from the following formula considering only the bending
deformations in the beam:

5 p ( L+ L)4 5 10⋅(2+2)4
e max = = =0.03292 m=32.92 mm
384 EI 384 30000000⋅0.12⋅0.153 /12

57
Verification Examples FEM-Design 16.0

The relevant results with FEM-Design:

Figure 7.1.2.5 – The deflection of the beam when the middle support only bear tension (total load)

The maximum of the shear force without sign:

1 1
V= p ( L+L)= 10(2+2)=20 kN
2 2

The relevant results with FEM-Design:

Figure 7.1.2.6 – The shear diagram of the beam when the middle support only bear tension (total load)

The extremum of the bending moment without sign:

1 1
M max = p (L+ L)2= 10⋅(2+2)2=20 kNm
8 8

The relevant results with FEM-Design:

Figure 7.1.2.7 – The bending moment diagram of the beam when the middle support only bear tension (total load)

The differences between the calculated results by hand and by FEM-Design are less than 2%.

58
Verification Examples FEM-Design 16.0

b) Case II.
In this case the distributed load is a partial load (Fig. 7.1.2.8). In the first part all of the three
supports behave the same way for compression and tension. In the second part of this case the
right side support only bears compression. We calculate in both cases the deflections, shear
forces and bending moments by hand and compared the results with FEM-Design calculations.

Figure 7.1.2.8 – The beam with three supports and uniform partial load

The extremums of the deflection come from the following formulas considering only the
bending deformations in the beam (without signs):

2.1 ( p /2)L 4 5 ( p/ 2) L 4 2.1 10/2⋅24 5 10/ 2⋅2 4

e max ≈ + = + =0.001461 m=1.461 mm
384 EI 384 EI 384 EI 384 EI
5 ( p / 2) L4 2 ( p /2) L 4 5 10/2⋅24 2 10 /2⋅2 4
e min ≈ − = − =0.0006173 m=0.6173 mm
384 EI 384 EI 384 EI 384 EI
The relevant results with FEM-Design:

Figure 7.1.2.9 – The deflection of the beam with three supports (partial load)

The extremums of the shear force without signs:

7 7 9 9
V 1= p L= 10⋅2=8.75 kN ; V 2 = p L= 10⋅2=11.25 kN ;
16 16 16 16
1 1
V 3= p L= 10⋅2=1.25 kN
16 16

The relevant results with FEM-Design:

59
Verification Examples FEM-Design 16.0

Figure 7.1.2.10 – The shear diagram of the beam with three supports (partial load)

The extremums of the bending moment without signs:

2 2

M midspan =
7
16 (
pL
=
7
16 ) (
10⋅2 ) 1 2 1 2
=3.828 kNm ; M middle = p L = 10⋅2 =2.5 kNm
2p 2⋅10 16 16

The relevant results with FEM-Design:

Figure 7.1.2.11 – The bending moment diagram of the beam with three supports (partial load)

When the right side support only bear compression (second part of this case) basically under the
partial vertical load (Fig. 7.1.2.8) the right side support is not active (support nonlinearity).
Therefore it works as a simply supported beam with two supports. The deflection, the shear
forces and the bending moments are the following:
The maximum deflection comes from the following formula considering only the bending
deformations in the beam:

5 p L4 5 10⋅24
e midspan = = =0.002058 m =2.058 mm
384 EI 384 EI
1 p L4 1 10⋅24
e right = = =0.006584 m=6.584 mm
24 EI 24 EI

60
Verification Examples FEM-Design 16.0

The relevant results with FEM-Design:

Figure 7.1.2.12 – The deflection of the beam when the right support only bear compression (partial load)

The extremum of the shear force without sign:

1 1
V= p ( L)= 10⋅2=10 kN
2 2

The relevant results with FEM-Design:

Figure 7.1.2.13 – The shear diagram of the beam when the right support only bear compression (partial load)

The extremum of the bending moment without sign:

1 1
M max = p L 2= 10⋅22=5.0 kNm
8 8

The relevant results with FEM-Design:

Figure 7.1.2.14 – The bending moment diagram of the beam when the right support only bear compression
(partial load)

The differences between the calculated results by hand and by FEM-Design are less than 2%.

61
Verification Examples FEM-Design 16.0

7.2 Cracked section analysis by reinforced concrete elements

7.2.1 Cracked deflection of a simply supported beam

Inputs:
Span length Leff = 7.2 m
The cross section Rectangle: b = 300 mm; h = 450 mm
The elastic modulus of concrete Ecm = 31.476 GPa, C25/30
The creep factor φ28 = 2.35
Effective elastic modulus of concrete Eceff = Ecm/(1+φ28) = 9.396 GPa
Mean tensile strength fctm = 2.565 MPa
Elastic modulus of steel bars Es = 200 GPa
Characteristic value of dead load gk = 8.5 kN/m
Characteristic value of live load qk = 12.0 kN/m
Live load combination factor ψ2 = 0.6
Diameter of the longitudinal reinforcement ϕl = 18 mm
Diameter of the stirrup reinforcement ϕs = 8 mm
Area of longitudinal reinforcement Al = 5x182π/4 = 1272.3 m2
Nominal concrete cover cnom = 20 mm

The cross sectional properties without calculation details:

I. stress stadium second moment of inertia II = 3.075x109 mm4
II. stress stadium second moment of inertia III = 2.028x109 mm4
I. stress stadium position of neutral axis xI = 256.4 mm
II. stress stadium position of neutral axis xII = 197.3 mm

ψ2qk
gk

Leff
5ϕ18
M

Mmax

Figure 7.2.1.1 – The simply supported RC beam

62
Verification Examples FEM-Design 16.0

The calculation of deflection according to EN 1992-1-1:

The load value for the quasi-permanent load combination:
kN
p qp=g k +ψ 2 qk =8.5+0.6⋅12=15.7
m
The maximum deflection with cross sectional properties in Stadium I. (uncracked):
4
5 pqp Leff 5 15.7⋅7.2 4
w k.I = = =0.01901 m=19.01 mm
384 E ceff I I 384 9396000⋅0.003075
The maximum deflection with cross sectional properties in Stadium II. (cracked):
4 4
5 p qp L eff 5 15.7⋅7.2
w k.II = = =0.02883 m=28.83 mm
384 E ceff I II 384 9396000⋅0.002028

The maximum bending moment under the quasi-permanent load:

1 1
M max = p qp L eff 2= 15.7⋅7.2 2=101.74 kNm
8 8
The cracking moment with the mean tensile strength:
II 0.003075
M cr = f ctm =2565 =40.74 kNm
h− x I 0.45−0.2564
The interpolation factor considering the mixture of cracked and uncracked behaviour:

[ ( )] [
2

) , 0]=0.9198
2
M cr
ζ =max 1−0.5
M max
, 0 =max 1−0.5
40.74
101.74 (
This value is almost 1.0, it means that the final deflection will be closer to the cracked deflection
than to the uncracked one.
The final deflection with the aim of interpolation factor:
w k =(1−ζ )w k.I +ζ w k.II =(1−0.9198)19.01+0.9198⋅28.83=28.04 mm

63
Verification Examples FEM-Design 16.0

First we modelled the beam with beam elements. In FEM-Design we increased the division
number of the beam finite elements to five to get the more accurate results.

Figure 7.2.1.2 – The cross section of the RC beam in FEM-Design

Fig. 7.2.1.2 shows the applied cross section and reinforcement with the defined input
parameters.

Figure 7.2.1.3 – The deflection of the RC beam in FEM-Design with cracked section analysis

Fig. 7.2.1.3 shows the deflection after the cracked section analysis. The deflection of the beam
model in FEM-Design:
w kFEM =30.59 mm
The difference is less than 9 %.

64
Verification Examples FEM-Design 16.0

Secondly we modelled the beam with shell finite elements. Fig. 7.2.1.4 shows the applied
specific reinforcement with the defined input parameters with slab.

Figure 7.2.1.4 – The specific reinforcement with the shell model in FEM-Design

Figure 7.2.1.5 – The deflection of the RC shell model in FEM-Design with cracked section analysis

Fig. 7.2.1.5 shows the deflection and the finite element mesh after the cracked section analysis.
The deflection of the shell model in FEM-Design:
w kFEM =27.0 mm
The difference is less than 4 %.

65
Verification Examples FEM-Design 16.0

7.2.2 Cracked deflection of a cantilever beam

Inputs:
Span length Leff = 4 m
The cross section Rectangle: b = 300 mm; h = 450 mm
The elastic modulus of concrete Ecm = 31.476 GPa, C25/30
The creep factor φ28 = 2.35
Effective elastic modulus of concrete Eceff = Ecm/(1+φ28) = 9.396 GPa
Mean tensile strength fctm = 2.565 MPa
Elastic modulus of steel bars Es = 200 GPa
Characteristic value of dead load gk = 8.5 kN/m
Characteristic value of live load qk = 12.0 kN/m
Live load combination factor ψ2 = 0.6
Diameter of the longitudinal reinforcement ϕl = 18 mm
Diameter of the stirrup reinforcement ϕs = 8 mm
Area of longitudinal reinforcement Al = 5x182π/4 = 1272.3 m2
Nominal concrete cover cnom = 20 mm

The cross sectional properties without calculation details:

I. stress stadium second moment of inertia II = 3.075x109 mm4
II. stress stadium second moment of inertia III = 2.028x109 mm4
I. stress stadium position of neutral axis xI = 256.4 mm
II. stress stadium position of neutral axis xII = 197.3 mm

ψ2qk 5ϕ18

gk

Leff

Mmax

M
Figure 7.2.2.1 – The cantilever RC beam

66
Verification Examples FEM-Design 16.0

The calculation of deflection according to EN 1992-1-1:

The load value for the quasi-permanent load combination:
kN
p qp=g k +ψ 2 qk =8.5+0.6⋅12=15.7
m
The maximum deflection with cross sectional properties in Stadium I. (uncracked):
4
1 p qp Leff 1 15.7⋅4 4
w k.I = = =0.01739 m=17.39 mm
8 E ceff I I 8 9396000⋅0.003075
The maximum deflection with cross sectional properties in Stadium II. (cracked):
4 4
1 pqp Leff 1 15.7⋅4
w k.II = = =0.02637 m=26.37 mm
8 E ceff I II 8 9396000⋅0.002028
The maximum bending moment under the quasi-permanent load:
1 1
M max = pqp Leff 2= 15.7⋅42 =125.6 kNm
2 2
The cracking moment with the mean tensile strength:
II 0.003075
M cr = f ctm =2565 =40.74 kNm
h− x I 0.45−0.2564
The interpolation factor considering the mixture of cracked and uncracked behaviour:

[ ( )] [
2

) ]
2
M cr
ζ =max 1−0.5
M max
, 0 =max 1−0.5
40.74
125.6 (
, 0 =0.9474

This value is almost 1.0, it means that the final deflection will be closer to the cracked deflection
than to the uncracked one.
The final deflection with the aim of interpolation factor:
w k =(1−ζ )w k.I +ζ w k.II =(1−0.9474) 17.39+0.9474⋅26.37=25.90 mm

67
Verification Examples FEM-Design 16.0

We modelled the beam with beam finite elements. In FEM-Design we increased the division
number of the beam finite elements to five to get the more accurate results.

Figure 7.2.2.2 – The cross section of the RC cantilever in FEM-Design

Fig. 7.2.2.2 shows the applied cross section and reinforcement with the defined input
parameters.

Figure 7.2.2.3 – The deflection of the RC cantilever in FEM-Design with cracked section analysis

Fig. 7.2.2.3 shows the deflection after the cracked section analysis. The deflection of the beam
model in FEM-Design:
w kFEM =27.54 mm
The difference is less than 7 %.

68
Verification Examples FEM-Design 16.0

7.3 Nonlinear soil calculation

This chapter goes beyond the scope of this document, therefore additional informations are
located in:
FEM-Design – Geotechnical modul in 3D, Theoretical background and verification and
validation handbook
http://download.strusoft.com/FEM-Design/inst150x/documents//3dsoilmanual.pdf

69
Verification Examples FEM-Design 16.0

8. Cross section editor

8.1 Calculation of a compound cross section

An example for compound cross section is taken from [7] where the authors calculated the cross
sectional properties with the assumption of thin-walled simplifications. The welded cross section
is consisting of U300 and L160x80x12 (DIN) profiles. In the Section Editor the exact cold
rolled geometry was analyzed as it is seen in Figure 8.1.1.

Figure 8.1.1 – The analyzed cross section

The following table contains the results of the two independent calculations with several cross
sectional properties.

Notation Ref. [1] Section Editor

A [cm2] 86.76 86.30
yG [cm] 1.210 1.442
zG [cm] 19.20 19.22
y'S [cm] 1.39 0.7230
z'S [cm] 10.06 10.36
Iy [cm4] 11379.9 11431.2
Iz [cm4] 4513.3 4372.9
4
Iyz [cm ] 3013.2 3053.5
4
It [cm ] 48.83 52.11
Iw [cm6] - 203082.0
Table 8.1.1 – The results of the example

70
Verification Examples FEM-Design 16.0

9. Design calculations

9.1 Required reinforcement calculation for a slab

In this example we calculate the required reinforcements of a slab due to elliptic and hyperbolic
bending conditions. First of all the applied reinforcement is orthogonal and then the applied
reinforcement is non-orthogonal. We calculate the required reinforcement with hand calculation
and then compare the results with FEM-Design values.
Inputs:
The thickness h = 200 mm
The elastic modulus of concrete Ecm = 33 GPa, C30/37
The Poisson's ratio of concrete ν = 0.2
The design value of compressive strength fcd = 20 MPa
Elastic modulus of steel bars Es = 200 GPa
The design value of yield stress of steel bars fyd = 434.8 MPa
Diameter of the longitudinal reinforcement ϕl = 10 mm
Nominal concrete cover cx = 20 mm; cy = 30 mm
Effective heights dx = 175 mm; dy = 165 mm

I.) Elliptic bending

In the first case the bending condition is an elliptic bending. In FEM-Design the model is a slab
with statically determinant support system and specific moment loads at its edges for the pure
internal force state (see Fig. 9.2.2.1).

Figure 9.2.2.1 – The slab with the edge loads for pure stress state

71
Verification Examples FEM-Design 16.0

Fig. 9.2.2.2 shows the constant internal forces in the slab due to the loads. Fig. 9.2.2.3 shows the
principal moments and their directions based on the FEM-Design calculation. According to the
pure stress state the principal moments and the directions are the same in each elements.

Figure 9.2.2.2 – The mx, my, and mxy internal forces in the slab [kNm]

Figure 9.2.2.3 – The m1 and m2 principal moments and their directions in the slab [kNm]

First of all the reinforcement is orthogonal and the hand calculation and the comparison are the
following:
1. Orthogonal reinforcement (φ=90o)
The reinforcement is orthogonal and their directions concide with the local system (x=ξ,
y=ϑ=η).
The moments in the slab (tensor of the applied moments):
m x =m ξ =+16 kNm /m
m y =mϑ =mη =+8 kNm / m
m xy =mξ ϑ =mξ η =+6 kNm /m

72
Verification Examples FEM-Design 16.0

The first invariant of the tensor: m x +m y =+24 kNm / m

The calculation of the principal moments and their directions:

√( √(
2

)
2
m +m y mx −m y
m 1= x
2
+
2
+mxy 2 =
16+8
2
+
16−8
2 )
+6 2=19.21 kNm / m

√( √(
2

)
2
m +m y m x −m y
2 )
2 16+8 16−8 2
m2 = x − +m xy = − +6 =4.79 kNm / m
2 2 2
m1−mx 19.21−16
α 0 =arctan =arctan =28.15o
m xy 6
Compare these results with Fig. 9.2.2.3. The difference is 0 %.

The design moments (according to [9][10]) if the reinforcement (ξ,η) is orthogonal and their
directions concide with the local co-ordinate system (x,y):
Case a)
cos φ 1−2 cos φ cos 90o 1−2 cos 90 o
mud ξ =mξ −mϑ +mξ ϑ =16−8 +6 =+22 kNm /m
1+cos φ sin φ 1+cos 90 o sin 90 o
1 1 1 1
m ud η =mϑ +m ξ ϑ =8 +6 =+14 kNm /m
1+cos φ sin φ 1+cos 90 o
sin 90o
This is a valid solution! Because mud ξ +mud η =+36 kNm /m>m x +m y =+24 kNm /m
mud ξ =+22 kNm /m mud η =+14 kNm / m

Case b)
cos φ 1+2 cos φ cos 90o 1+2 cos 90 o
mud ξ =mξ +mϑ −mξ ϑ =16+8 −6 =+10 kNm /m
1−cos φ sin φ 1−cos 90 o sin 90 o
1 1 1 1
m ud η =mϑ −m ξ ϑ =8 −6 =+2 kNm /m
1−cos φ sin φ 1−cos 90 o
sin 90o
Invalid solution! Because mud ξ +mud η =+12 kNm /m<m x +m y =+24 kNm /m

Case ξ)

m 2ξ ϑ 62
mud ξ =mξ − =16− =+11.5 kNm/ m
mϑ 8
mud η =0
Invalid solution! Because mud ξ +mud η =+11.5 kNm /m<mx +m y =+24 kNm /m

73
Verification Examples FEM-Design 16.0

Case η)
m ud ξ =0
mξ mϑ −m 2ξ ϑ 16⋅8−62 kNm
m ud η = 2 2
= 2 o 2 o o
=+5.75
m ξ sin φ +m ϑ cos φ −m ξ ϑ sin 2 φ 16⋅sin 90 +8⋅cos 90 −6⋅sin (2⋅90 ) m
Invalid solution! Because mud ξ +mud η =+5.75 kNm /m<m x +m y =+24 kNm / m
The results of the design moments based on FEM-Design are in Fig. 9.2.2.4 and 9.2.2.5. The
difference between the hand and FE calculation is 0%.

Figure 9.2.2.4 – The mud ξ design moment for elliptic bending with orthogonal reinforcement [kNm]

Figure 9.2.2.5 – The mud η design moment for elliptic bending with orthogonal reinforcement [kNm]

74
Verification Examples FEM-Design 16.0

Calculation of the required reinforcement based on the valid design moments:

In x (ξ) direction:
Sum of the moments:

(
mud ξ = f cd x c d x −
xc
2 ) (
; 22000=20 x c 175−
xc
2 ) ; x c =6.403 mm

Sum of the forces:

x c f cd =a sξ f yd ; 6.403⋅20=a sξ 434.8 ; a s ξ =0.2945 mm 2 /mm=294.5 mm 2 /m
Fig. 9.2.2.6 shows the required reinforcement in the relevant direction based on FEM-Design
calculation. The difference is less than 1%.

Figure 9.2.2.6 – The asξ required reinforcement at the bottom for elliptic bending with orthogonal
reinforcement [mm2/m]

In y (η) direction:
Sum of the moments:

(
mud η = f cd x c d y −
xc
2 ) (
; 14000=20 x c 165−
xc
2 ) ; x c =4.298 mm

Sum of the forces:

x c f cd =a sη f yd ; 4.298⋅20=a sη 434.8 ; a s η=0.1977 mm 2 /mm=197.7 mm2 / m
Fig. 9.2.2.7 shows the required reinforcement in the relevant direction based on FEM-Design
calculation. The difference is less than 1%.

75
Verification Examples FEM-Design 16.0

Figure 9.2.2.7 – The asη required reinforcement at the bottom for elliptic bending with orthogonal
reinforcement [mm2/m]

Secondly the reinforcement is non-orthogonal and the hand calculation and the comparison are
the following:
2. Non-orthogonal reinforcement (φ=75o between ξ and η)
The reinforcement is non-orthogonal and the ξ direction concides with the local x direction.
Thus y=ϑ. The angle between the ξ directional reinforcement and η directional reinforcement is
φ=75o.
The moments in the slab (tensor of the applied moments):
m x =m ξ =+16 kNm /m
m y =mϑ =mη =+8 kNm / m
m xy =mξ ϑ =+6 kNm /m
The first invariant of the tensor: m x +m y =+24 kNm / m

The design moments (according to [9][10]) if the reinforcement (ξ,η) is non-orthogonal:

Case a)
cos φ 1−2 cos φ cos 75o 1−2 cos 75o
mud ξ =mξ −mϑ +mξ ϑ =16−8 +6 =+17.35 kNm /m
1+cos φ sin φ 1+cos 75o sin 75o
1 1 1 1
m ud η =mϑ +m ξ ϑ =8 +6 =+12.57 kNm / m
1+cos φ sin φ 1+cos 75o
sin 75o
This is a valid solution! Because mud ξ +mud η =+29.92 kNm / m>mx +m y =+24 kNm /m
mud ξ =+17.35 kNm / m mud η =+12.57 kNm / m

76
Verification Examples FEM-Design 16.0

Case b)
cos φ 1+2 cos φ cos 75o 1+2 cos 75o
mud ξ =mξ +mϑ −mξ ϑ =16+8 −6 =+9.37 kNm /m
1−cos φ sin φ 1−cos 75o sin 75o
1 1 1 1
m ud η =mϑ −m ξ ϑ =8 −6 =+4.58 kNm /m
1−cos φ sin φ 1−cos 75 o
sin 75o
Invalid solution! Because mud ξ +mud η =+13.95 kNm /m<mx +m y =+24 kNm /m

Case ξ)
2
mξ ϑ 62
mud ξ =mξ − =16− =+11.5 kNm/ m
mϑ 8
mud η =0
Invalid solution! Because mud ξ +mud η =+11.5 kNm /m<mx +m y =+24 kNm /m

Case η)
m ud ξ =0
mξ mϑ −m 2ξ ϑ 16⋅8−6 2 kNm
m ud η = 2 2
= 2 o 2 o o
=+7.38
m ξ sin φ +m ϑ cos φ −m ξ ϑ sin 2 φ 16⋅sin 75 +8⋅cos 75 −6⋅sin(2⋅75 ) m
Invalid solution! Because mud ξ +mud η =+7.38 kNm /m <m x +m y =+24 kNm / m
The results of the design moments based on FEM-Design are in Fig. 9.2.2.8 and 9.2.2.9. The
difference between the hand and FE calculation is 0%.

Figure 9.2.2.8 – The mud ξ design moment for elliptic bending with non-orthogonal reinforcement [kNm]

77
Verification Examples FEM-Design 16.0

Figure 9.2.2.9 – The mud η design moment for elliptic bending with non-orthogonal reinforcement [kNm]

Calculation of the required reinforcement based on the valid design moments:

In x (ξ) direction:
Sum of the moments:

(
mud ξ = f cd x c d x −
xc
2 ) (
; 17350=20 x c 175−
xc
2 ) ; x c =5.029 mm

Sum of the forces:

x c f cd =a sξ f yd ; 5.029⋅20=a sξ 434.8 ; a s ξ =0.2313 mm 2 /mm=231.3 mm 2 /m
Fig. 9.2.2.10 shows the required reinforcement in the relevant direction based on FEM-Design
calculation. The difference is less than 1%.

Figure 9.2.2.10 – The asξ required reinforcement at the bottom for elliptic bending with non-orthogonal
reinforcement [mm2/m]

78
Verification Examples FEM-Design 16.0

In η direction:
Sum of the moments:

(
mud η = f cd x c d y −
xc
2 ) (
; 12570=20 x c 165−
xc
2 ) ; x c =3.854 mm

Sum of the forces:

x c f cd =a sη f yd ; 3.854⋅20=a sη 434.8 ; a s η=0.1773 mm 2 /mm=177.3 mm 2 / m
Fig. 9.2.2.11 shows the required reinforcement in the relevant direction based on FEM-Design
calculation. The difference is less than 1%.

Figure 9.2.2.11 – The asη required reinforcement at the bottom for elliptic bending with non-orthogonal
reinforcement [mm2/m]

79
Verification Examples FEM-Design 16.0

II.) Hyperbolic bending

In the second case the bending condition is a hyperbolic bending. In FEM-Design the model is a
slab with statically determinant support system and specific moment loads at its edges for the
pure internal force state (see Fig. 9.2.2.12).

Figure 9.2.2.12 – The slab with the edge loads for pure stress state

Fig. 9.2.2.13 shows the constant internal forces in the slab due to the loads. Fig. 9.2.2.14 shows
the principal moments and their directions based on the FEM-Design calculation. According to
the pure stress state the principal moments and the directions are the same in each elements.

Figure 9.2.2.13 – The mx, my, and mxy internal forces in the slab [kNm]

80
Verification Examples FEM-Design 16.0

Figure 9.2.2.14 – The m1 and m2 principal moments and their directions in the slab [kNm]

Firstly the reinforcement is orthogonal and the hand calculation and the comparison are the
following:
1. Orthogonal reinforcement
The reinforcement is orthogonal and their directions concide with the local system (x=ξ,
y=ϑ=η).
The moments in the slab (tensor of the applied moments):
m x =m ξ =+16 kNm /m
m y =mϑ =mη =−8 kNm / m
m xy =mξ ϑ =mξ η =+6 kNm /m
The first invariant of the tensor: m x +m y =+8 kNm /m

The calculation of the principal moments and their directions:

√( √(
2 2
m 1=
m x +m y
2
+
2 )
mx −m y
+mxy 2 =
16+(−8)
2
+
16−(−8)
2 )
+62=17.42 kNm /m

√( √(
2 2
m +m y
m2 = x
2
−
2 )
m x −m y
+m xy2=
16+(−8)
2
−
16−(−8)
2 )
+6 2=−9.42 kNm /m

17.42−16 19.21−16
α 0 =arctan =arctan =13.32o
6 6
Compare these results with Fig. 9.2.2.14. The difference is 0 %.
The design moments (according to [9][10]) if the reinforcement (ξ,η) is non-orthogonal:

81
Verification Examples FEM-Design 16.0

Case a)
cos φ 1−2 cos φ cos 90o 1−2 cos 90 o
mud ξ =mξ −mϑ +mξ ϑ =16+8 +6 =+22 kNm /m
1+cos φ sin φ 1+cos 90 o sin 90 o
1 1 1 1
m ud η =mϑ +m ξ ϑ =−8 +6 =−2 kNm / m
1+cos φ sin φ 1+cos 90 o
sin 90 o
Invalid solution! Because their have different signs.

Case b)
cos φ 1+2 cos φ cos 90o 1+2 cos 90 o
mud ξ =mξ +mϑ −mξ ϑ =16−8 −6 =+10 kNm /m
1−cos φ sin φ 1−cos 90 o sin 90 o
1 1 1 1
m ud η =mϑ −m ξ ϑ =−8 −6 =−14 kNm /m
1−cos φ sin φ 1−cos 90 o
sin 90 o
Invalid solution! Because their have different signs.

Case ξ)
m 2ξ ϑ 62
mud ξ =mξ − =16− =+20.5 kNm /m
mϑ −8
mud η =0
This is a valid solution at the bottom!
mud ξ =+20.5 kNm /m mud η =0 kNm /m

Case η)
mud ξ =0
mξ mϑ −m 2ξ ϑ 16⋅(−8)−6 2 kNm
m ud η = 2 2
= 2 o 2 o o
=−10.25
m ξ sin φ +m ϑ cos φ −m ξ ϑ sin 2 φ 16⋅sin 90 +(−8)⋅cos 90 −6⋅sin( 2⋅90 ) m
This is a valid solution at the top!
mud ξ =0 kNm /m mud η =−10.25 kNm /m

The results of the design moments based on FEM-Design are in Fig. 9.2.2.15 and 9.2.2.16. The
difference between the hand and FE calculation is 0%.

82
Verification Examples FEM-Design 16.0

Figure 9.2.2.15 – The mud ξ design moment for hyperbolic bending with orthogonal reinforcement [kNm]

Figure 9.2.2.16 – The mud η design moment for hyperbolic bending with orthogonal reinforcement [kNm]

Calculation of the required reinforcement based on the valid design moments:

In x (ξ) direction at the bottom:
Sum of the moments:

(
mud ξ = f cd x c d x −
xc
2 ) (
; 20500=20 x c 175−
xc
2 ) ; x c =5.959 mm

Sum of the forces:

x c f cd =a sξ f yd ; 5.959⋅20=a sξ 434.8 ; a s ξ =0.2741 mm 2 /mm=274.1 mm 2 /m
Fig. 9.2.2.17 shows the required reinforcement in the relevant direction based on FEM-Design
calculation. The difference is less than 1%.

83
Verification Examples FEM-Design 16.0

Figure 9.2.2.17 – The asξ required reinforcement at the bottom for hyperbolic bending with orthogonal
reinforcement [mm2/m]

In y (η) direction at the top:

Sum of the moments:

(
mud η = f cd x c d y −
xc
2 ) (
; 10250=20 x c 165−
xc
2 ) ; x c =3.136 mm

Sum of the forces:

x c f cd =a sη f yd ; 3.136⋅20=a sη 434.8 ; a s η=0.1443 mm 2 / mm=144.3 mm 2 / m
Fig. 9.2.2.18 shows the required reinforcement in the relevant direction based on FEM-Design
calculation. The difference is less than 1%.

Figure 9.2.2.18 – The asη required reinforcement at the top for hyperbolic bending with orthogonal
reinforcement [mm2/m]

84
Verification Examples FEM-Design 16.0

Secondly the reinforcement is non-orthogonal and the hand calculation and the comparison are
the following:
2. Non-orthogonal reinforcement (φ=75o between ξ and η)
The reinforcement is non-orthogonal and the ξ direction concides with the local x direction.
Thus y=ϑ.
The moments in the slab (tensor of the applied moments):
m x =m ξ =+16 kNm /m
m y =mϑ =mη =−8 kNm / m
m xy =mξ ϑ =+6 kNm /m
The first invariant of the tensor: m x +m y =+8 kNm /m

The design moments (according to the theory book) if the reinforcement (ξ,η) is orthogonal and
their directions concide with the local co-ordinate system (x,y):
Case a)
o o
cos φ 1−2 cos φ cos 75 1−2 cos 75
mud ξ =mξ −mϑ +mξ ϑ =16+8 +6 =+20.64 kNm /m
1+cos φ sin φ 1+cos 75o sin 75o
1 1 1 1
m ud η =mϑ +m ξ ϑ =−8 +6 =−0.144 kNm / m
1+cos φ sin φ 1+cos 75
o
sin 75
o

Invalid solution! Because their have different signs.

Case b)
o o
cos φ 1+2 cos φ cos 75 1+2 cos 75
mud ξ =mξ +mϑ −mξ ϑ =16−8 o
−6 o
=+3.78 kNm /m
1−cos φ sin φ 1−cos 75 sin 75
1 1 1 1
m ud η =mϑ −m ξ ϑ =−8 −6 =−17.01 kNm /m
1−cos φ sin φ 1−cos 75 o
sin 75o
Invalid solution! Because their have different signs.

Case ξ)
2
mξ ϑ 62
mud ξ =mξ − =16− =+20.5 kNm /m
mϑ −8
mud η =0
This is a valid solution at the bottom!
mud ξ =+20.5 kNm /m mud η =0 kNm /m

Case η)

85
Verification Examples FEM-Design 16.0

mud ξ =0
mξ mϑ −m 2ξ ϑ 16⋅(−8)−62 kNm
m ud η = 2 2
= 2 o 2 o o
=−14.40
m ξ sin φ +m ϑ cos φ −m ξ ϑ sin 2 φ 16⋅sin 75 +(−8)⋅cos 75 −6⋅sin( 2⋅75 ) m
This is a valid solution at the top!
mud ξ =0 kNm /m mud η =−14.40 kNm /m
The results of the design moments based on FEM-Design are in Fig. 9.2.2.19 and 9.2.2.20. The
difference between the hand and FE calculation is 0%.

Figure 9.2.2.19 – The mud ξ design moment for hyperbolic bending with non-orthogonal reinforcement [kNm]

Figure 9.2.2.20 – The mud η design moment for hyperbolic bending with non-orthogonal reinforcement [kNm]

86
Verification Examples FEM-Design 16.0

Calculation of the required reinforcement based on the valid design moments:

In x (ξ) direction at the bottom:
Sum of the moments:

(
mud ξ = f cd x c d x −
xc
2 ) (
; 20500=20 x c 175−
xc
2 ) ; x c =5.959 mm

Sum of the forces:

x c f cd =a sξ f yd ; 5.959⋅20=a sξ 434.8 ; a s ξ =0.2741 mm 2 /mm=274.1 mm 2 /m
Fig. 9.2.2.21 shows the required reinforcement in the relevant direction based on FEM-Design
calculation. The difference is less than 1%.

Figure 9.2.2.21 – The asξ required reinforcement at the bottom for hyperbolic bending with non-orthogonal
reinforcement [mm2/m]

In η direction at the top:

Sum of the moments:

(
mud η = f cd x c d y −
xc
2 ) (
; 14400=20 x c 165−
xc
2 ) ; x c =4.423 mm

Sum of the forces:

x c f cd =a sη f yd ; 4.423⋅20=a sη 434.8 ; a s η=0.2034 mm 2 /mm=203.4 mm 2 / m
Fig. 9.2.2.22 shows the required reinforcement in the relevant direction based on FEM-Design
calculation. The difference is less than 1%.

87
Verification Examples FEM-Design 16.0

Figure 9.2.2.22 – The asη required reinforcement at the top for hyperbolic bending with non-orthogonal
reinforcement [mm2/m]

88
Verification Examples FEM-Design 16.0

References
[1] Beer F.P., Johnston E.R., DeWolf J.T., Mazurek D.F.: Mechanics of materials, McGraw-Hill,
2012.
[2] Timoshenko S., Woinowsky-Krieger S., Theory of plates and shells, McGraw-Hill, New
York, 1959
[3] Ventsel E., Krauthammer T., Thin plates and shells, Marcel Dekker, New York, 2001
[4] Dulácska E., Joó A., Kollár L., Tartószerkezetek tervezése földrengési hatásokra, Akadémiai
Kiadó, Budapest, 2008.
[5] Sadd M.H., Wave motion and vibration in continuous media, Kingston, Rhode Island, 2009.
[6] Iványi M., Halász O., Stabilitástan, Műegyetemi Kiadó, Budapest, 1995.
[7] Wagner W., Gruttman F., A displacement method for the analysis of flexural shear stresses in
thin walled isotropic composite beams, Computers and Structures, Vol. 80., pp. 1843-1851.,
2002.
[8] Majid K.I., Non-linear structures: matrix methods of analysis and design by computers,
London, Butterworths, 1972.
[9] Németh F., Optimum design of steel bars in reinforced concrete slabs and the criticism of the
Eurocode 2, Vasbetonépítés, 3., pp. 107-114., 2001.
[10] Wood, R.H., The Reinforcement of Slabs in Accordance with a Pre-Determined Field of
Moments, Concrete, Vol. 2, No. 2, pp. 69-76., 1968.

89
Verification Examples FEM-Design 16.0

Notes

90