Verification Examples

2008
AxisVM 9 Verification Examples 2



Linear static ............................................................................................................. 3
Supported bar with concentrated loads. ....................................................................................................................... 4
Thermally loaded bar structure. .................................................................................................................................... 5
Continously supported beam with constant distributed load......................................................................................... 6
External prestessed beam........................................................................................................................................... 9
Periodically supported infinite membrane wall with constant distributed load. ........................................................... 11
Clamped beam examination with plane stress elements............................................................................................ 13
Clamped thin square plate.......................................................................................................................................... 16
Plate with fixed support and constant distributed load................................................................................................ 18
Annular plate. ............................................................................................................................................................. 19
All edges simply supported plate with partial distributed load. ................................................................................... 21
Clamped plate with linear distributed load. ................................................................................................................. 23
Hemisphere displacement. ......................................................................................................................................... 25
Nonlinear static...................................................................................................... 27
3D beam structure. ..................................................................................................................................................... 28
Plate with fixed end and bending moment.................................................................................................................. 30
Dynamic................................................................................................................. 33
Deep simply supported beam..................................................................................................................................... 34
Clamped thin rhombic plate........................................................................................................................................ 37
Cantilevered thin square plate.................................................................................................................................... 39
Cantilevered tapered membrane. ............................................................................................................................... 42
Flat grillages. .............................................................................................................................................................. 45
Stability .................................................................................................................. 49
Simply supported beam.............................................................................................................................................. 50
Simply supported beam.............................................................................................................................................. 52
Design ................................................................................................................... 53
N-M interaction curve of cross-section (EN 1992-1-1:2004)....................................................................................... 54
RC beam deflection according to EC2, EN 1992-1-1:2004. ....................................................................................... 55
Required steel reinforcement of RC plate according to EC2, EN 1992-1-1:2004. ...................................................... 57
Earth-quake design using response-spectrum method. ............................................................................................. 59




AxisVM 9 Verification Examples 3

















































































Linear static
AxisVM 9 Verification Examples 4

Software Release Number: R2
Date: 12. 03. 2008.
Tested by: InterCAD
Page number:
File name: [Link]

Thema

Supported bar with concentrated loads.
Analysis
Type
Linear analysis.
Geometry



Side view

Section Area = 1,0 m
2


Loads

Axial direction forces P
1
= -200 N, P
2
= 100 N, P
3
= -40 N
Boundary
Conditions
Fix ends, at R
1
and R
5
.
Material
Properties
E = 20000 kN / cm
2

= 0,3
Element
types
Beam element
Mesh

Target

R
1
, R
5
support forces
Results
Theory AxisVM %
R
1
[N] 118,00 118,00 0,00
R
5
[N] 22,00 22,00 0,00




AxisVM 9 Verification Examples 5

Software Release Number: R2
Date: 12. 03. 2008.
Tested by: InterCAD
Page number:
File name: [Link]

Thema

Thermally loaded bar structure.
Analysis
Type
Linear analysis.
Geometry



Side view
Sections:
Steel: A
S
= x 10
-4
m
2

Copper: A
C
= x 10
-4
m
2


Loads

P = -12 kN (Point load)
Temperature rise of 10 C in the structure after assembly.
Boundary
Conditions
The upper end of bars are fixed.
Material
Properties
Steel: E
S
= 20700 kN / cm
2
, = 0,3 ,
S
= 1,2 x 10
-5
C
-1

Copper: E
C
= 11040 kN / cm
2
, = 0,3 ,
C
= 1,7 x 10
-5
C
-1


Element
types
Beam element
Target

S
max
in the three bars.


Results
Theory AxisVM %
Steel S
max
[MPa] 23824000 23848000 0,10
Cooper S
max
[MPa] 7185300 7198900 0,19



AxisVM 9 Verification Examples 6

Software Release Number: R2
Date: 12. 03. 2008.
Tested by: InterCAD
Page number:
File name: [Link]

Thema

Continously supported beam with point loads.
Analysis
Type
Linear analysis.
Geometry



Side view

(Section width = 1,00 m, height
1
= 0,30 m, height
2
= 0,60 m)

Loads

P
1
= -300 kN, P
2
= -1250 kN, P
3
= -800 kN, P
4
= -450 kN
Boundary
Conditions
Elastic supported.
From A to D is K
z
= 25000 kN/m/m.
From D to F is K
z
= 15000 kN/m/m.
Material
Properties
E = 3000 kN/cm
2

= 0,3
Element
types
Three node beam element. Shear deformation is taken into account.
Target e
z
, M
y
, V
z
, R
z

Results





























Results

Diagram e
z


Diagram M
y




AxisVM 9 Verification Examples 7


































































Diagram V
z


Diagram R



Reference AxisVM e [%]
e
A
[m] 0,006 0,006 0,00
e
B
[m] 0,009 0,009 0,00
e
C
[m] 0,014 0,014 0,00
e
D
[m] 0,015 0,015 0,00
e
E
[m] 0,015 0,015 0,00
e
F
[m] 0,013 0,013 0,00



Reference AxisVM e [%]
M
A
[KNm] 0,0 0,0 0,00
M
B
[KNm] 88,5 87,5 -1,13
M
C
[KNm] 636,2 632,8 -0,53
M
D
[KNm] 332,8 329,7 -0,93
M
E
[KNm] 164,2 163,3 -0,55
M
F
[KNm] 0,0 0,0 0,00









AxisVM 9 Verification Examples 8


Results




Reference AxisVM e [%]
V
A
[KN] 0,0 0,0 0,00
V
B
[KN] 112,1 113,1 0,89
V
C
[KN] 646,8 647,2 0,06
V
D
[KN] 335,0 334,9 -0,03
V
E
[KN] 267,8 267,5 -0,11
V
F
[KN] 0,0 0,0 0,00



Reference AxisVM e [%]
R
A
[KN/m
2
] 145,7 154,0 5,70
R
B
[KN/m
2
] 219,5 219,4 -0,05
R
C
[KN/m
2
] 343,8 346,0 0,64
R
D
[KN/m
2
] 386,9 386,4 -0,13
R
E
[KN/m
2
] 224,5 224,7 0,09
R
F
[KN/m
2
] 201,2 200,8 -0,20



AxisVM 9 Verification Examples 9

Software Release Number: R2
Date: 12. 03. 2008.
Tested by: InterCAD
Page number:
File name: [Link]

Thema

External prestessed beam.
Analysis
Type
Linear analysis.
Geometry




Side view

Loads

p = -50 kN /m distributed load
Length change = -6,52E-3 at beam 5-6

Boundary
Conditions
eY = eZ = = 0 at node 1
eX = eY = eZ = 0 at node 4

Material
Properties
E = 2,1E11 N / m
2

Beam 1-5, 5-6, 6-4 A = 4,5E-3 m
2
I
z
= 0,2E-5 m
4

Truss 2-5, 3-6 A = 3,48E-3 m
2
I
z
= 0,2E-5 m
4

Beam 1-4 A = 1,1516E-2 m
2
I
z
= 2,174E-4 m
4


Mesh

Element
types
Three node beam element, 1-5, 5-6, 6-4, 1-4 (shear deformation is taken into account)
Truss element 2-5, 3-6

Target

N
X
at beam 6-7
M
y,max
at beam 2-3
e
z
at node 2


AxisVM 9 Verification Examples 10


Results

1 2 3 4
5 6
2,000 4,000 2,000
8,000
0
,
6
0
0
X
Z

Diagram e
z






ROBOT V6

AxisVM %
N
x
[kN] 584,56 585,70 0,19
M
y
[kNm] 49,26 49,60 0,68
e
z
[mm] -0,5421 -0,5469 0,89



AxisVM 9 Verification Examples 11

Software Release Number: R2
Date: 12. 03. 2008.
Tested by: InterCAD
Page number:
File name: [Link]

Thema

Periodically supported infinite membrane wall with constant distributed load.
Analysis
Type
Linear analysis.
Geometry


Side view

(thickness = 20,0 cm)
Loads

p = 200 kN / m
Boundary
Conditions
vertical support at every 4,0 m support length is 0,4 m
Material
Properties
E = 880 kN / cm
2

= 0,16
Element
types
Parabolic quadrilateral membrane (plane stress)
Mesh

Target

S
xx
at 1-10 nodes (1-5 at middle, 6-10 at support)
AxisVM 9 Verification Examples 12


Results

Node
Analytical [kN/cm
2
] AxisVM [kN/cm
2
]
%
1 0,1313 0,1310 -0,23
2 0,0399 0,0395 -1,00
3 -0,0093 -0,0095 2,15
4 -0,0412 -0,0410 -0,49
5 -0,1073 -0,1070 -0,28
6 -0,9317 -0,9270 -0,50
7 0,0401 0,0426 6,23
8 0,0465 0,0469 0,86
9 0,0538 0,0540 0,37
10 0,1249 0,1240 -0,72


Reference:

Dr. Blcskey Elemr Dr. Orosz rpd:
Vasbeton szerkezetek Faltartk, Lemezek, Trolk


AxisVM 9 Verification Examples 13

Software Release Number: R2
Date: 12. 03. 2008.
Tested by: InterCAD
Page number:
File name: [Link]

Thema

Clamped beam examination with plane stress elements.
Analysis
Type
Linear analysis.
Geometry




Side view

Loads

p = -25 kN/m
Boundary
Conditions
Both ends built-in.

Material
Properties
E = 880 kN / cm
2

= 0
Element
types
Parabolic quadrilateral membrane (plane stress)
Mesh

C
Clamped
edge
1
3,000 0,250
0
,
5
0
0
0,375
C
Clamped
edge
X
Z

Side view
AxisVM 9 Verification Examples 14


Target

xy, max
at section C
Results

Diagram
xy



5,28
791,56
5,14
5,28
791,56
5,14
5,28
791,56
5,14
5,28
791,56
5,14
Y
Z

Diagram
xy
at section C
AxisVM 9 Verification Examples 15


2
'
4
3 '
/ 5 , 787
00260416 , 0 25 , 0
0078125 , 0 625 , 65
00260416 , 0
25 , 0
0078125 , 0
) ( 625 , 65
m kN
I b
S V
m I
m b
m S
theory beam from kN V
y
y
xy
y
y
=

=
=
=
=
=

AxisVM result
xy
= 791,6 kN / m
2


Difference = +0,52 %

AxisVM result
kN n V
xy
34 , 65 = =

Difference = +0,43 %


AxisVM 9 Verification Examples 16

Software Release Number: R2
Date: 12. 03. 2008.
Tested by: InterCAD
Page number:
File name: [Link]

Thema

Clamped thin square plate.
Analysis
Type
Linear analysis.
Geometry





Top view

(thickness = 5,0 cm)

Loads

P = -10 kN (at the middle of the plate)
Boundary
Conditions
eX = ez = eZ = fiX = fiY = fiZ = 0 along all edges

Material
Properties
E = 20000 kN / cm
2

= 0,3
Element
types
Plate element (Parabolic quadrilateral, heterosis)
Mesh
4,000
4
,
0
0
0
X
Y

Target

Displacement of middle of the plate
AxisVM 9 Verification Examples 17

Results
-0,001
-0,006
-0,012
-0,019
-0,024
-0,026
-0,026
-0,024
-0,019
-0,012
-0,006
-0,001
-0,006
-0,022
-0,043
-0,065
-0,081
-0,087
-0,087
-0,081
-0,065
-0,043
-0,022
-0,006
-0,012
-0,043
-0,084
-0,125
-0,156
-0,168
-0,168
-0,156
-0,125
-0,084
-0,043
-0,012
-0,019
-0,065
-0,125
-0,187
-0,237
-0,257
-0,257
-0,237
-0,187
-0,125
-0,065
-0,019
-0,024
-0,081
-0,156
-0,237
-0,307
-0,337
-0,337
-0,307
-0,237
-0,156
-0,081
-0,024
-0,026
-0,087
-0,168
-0,257
-0,337
-0,383
-0,383
-0,337
-0,257
-0,168
-0,087
-0,026
-0,026
-0,087
-0,168
-0,257
-0,337
-0,383
-0,383
-0,337
-0,257
-0,168
-0,087
-0,026
-0,024
-0,081
-0,156
-0,237
-0,307
-0,337
-0,337
-0,307
-0,237
-0,156
-0,081
-0,024
-0,019
-0,065
-0,125
-0,187
-0,237
-0,257
-0,257
-0,237
-0,187
-0,125
-0,065
-0,019
-0,012
-0,043
-0,084
-0,125
-0,156
-0,168
-0,168
-0,156
-0,125
-0,084
-0,043
-0,012
-0,006
-0,022
-0,043
-0,065
-0,081
-0,087
-0,087
-0,081
-0,065
-0,043
-0,022
-0,006
-0,001
-0,006
-0,012
-0,019
-0,024
-0,026
-0,026
-0,024
-0,019
-0,012
-0,006
-0,001
X
Y
Z


Displacements


Mode Mesh Book
1
Timoshenko
2
AxisVM Diff
1
[%] Diff
2
[%]
1 2x2 0,402 0,420 4,48 10,53
2 4x4 0,416 0,369 -11,30 -2,89
3 8x8 0,394 0,381 -3,30 0,26
4 12x12 0,387 0,383 -1,03 0,79
5 16x16 0,385 0,383 -0,52 0,79
0,38


References:

1.) The Finite Element Method (Fourth Edition) Volume 2.
/O.C. Zienkiewicz and R.L. Taylor/ McGraw-Hill Book Company 1991 London
2.) Result of analytical solution of Timoshenko





Convergency
-15,00
-10,00
-5,00
0,00
5,00
10,00
15,00
1 2 3 4 5
Mesh density
D
i
s
p
l
a
c
e
m
e
n
t
s
Diff1 [%]
Diff2 [%]

AxisVM 9 Verification Examples 18

Software Release Number: R2
Date: 12. 03. 2008.
Tested by: InterCAD
Page number:
File name: plate2_1.axs

Thema

Plate with fixed support and constant distributed load.
Analysis
Type
Linear analysis.
Geometry



Top view

(thickness = 15,0 cm)
Loads P = -5 kN / m
2

Boundary
Conditions
eX = eY = eZ = fiX = fiY = fiZ = 0 along all edges

Material
Properties
E = 990 kN/cm
2

= 0,16
Element
types
Parabolic triangle plate element
Mesh

Target Maximal eX (found at Node1) and maximal m
x
(found at Node2)
Results
Component Nastran AxisVM %
eZ,max [mm] -1,613 -1,593 -1,24
mx,max [kNm/m] 3,060 3,059 -0,03


AxisVM 9 Verification Examples 19

Software Release Number: R2
Date: 12. 03. 2008.
Tested by: InterCAD
Page number:
File name: [Link]

Thema

Annular plate.
Analysis
Type
Linear analysis.
Geometry



Top view

(thickness = 22,0 cm)

Loads

Edge load: Q = 100 kN / m

Distributed load: q = 100 kN / m
2
Boundary
Conditions


Material
Properties
E = 880 kN / cm
2

= 0,3
Element
types
Plate element (parabolic quadrilateral, heterosis)
AxisVM 9 Verification Examples 20

Mesh
4,000
1,000 3,000
X
Y

Target S
max
, e
max

Results
Theory AxisVM
Model
S
max
S
max %
[kN/cm2] [kN/cm2]
a.) 2,82 2,78 -1,42
b.) 6,88 6,76 -1,74
c.) 14,22 14,10 -0,84
d.) 1,33 1,33 0,00
e.) 2,35 2,25 -4,26
f.) 9,88 9,88 0,00
g.) 4,79 4,76 -0,63
h.) 7,86 7,86 0,00


Theory AxisVM
Model
e
max
e
max %
[mm] [mm]
a.) 77,68 76,10 -2,03
b.) 226,76 220,84 -2,61
c.) 355,17 352,89 -0,64
d.) 23,28 23,42 0,60
e.) 44,26 44,50 0,54
f.) 123,19 123,17 -0,02
g.) 112,14 111,94 -0,18
h.) 126,83 126,81 -0,02


Reference:

S. Timoshenko and S. Woinowsky-Krieger: Theory of Plates And Shells


AxisVM 9 Verification Examples 21

Software Release Number: R2
Date: 12. 03. 2008.
Tested by: InterCAD
Page number:
File name: [Link]

Thema

All edges simply supported plate with partial distributed load.
Analysis
Type
Linear analysis.
Geometry



Top view

(thickness = 22,0 cm)

Loads

Distributed load: q = -10 kN / m
2
(middle of the plate at 2,0 x 2,0 m area)
Boundary
Conditions
a.) eX = eY = eZ = 0 along all edges (soft support)
b.) eX = eY = eZ = 0 along all edges = 0 perpendicular the edges (hard support)
Material
Properties
E = 880 kN / cm
2

= 0,3
Element
types
Plate element (Heterosis)
Mesh
5,000
1
0
,
0
0
0
X
Y


AxisVM 9 Verification Examples 22


Target

m
x, max
, m
y, max

Results a.)

Moment Theory AxisVM %
m
x, max
[kNm/m]
7,24 7,34 1,38
m
y, max
[kNm/m]
5,32 5,39 1,32


b.)


Moment Theory AxisVM %
m
x, max
[kNm/m]
7,24 7,28 0,55
m
y, max
[kNm/m]
5,32 5,35 0,56


Reference:

S. Timoshenko and S. Woinowsky-Krieger: Theory of Plates And Shells


AxisVM 9 Verification Examples 23

Software Release Number: R2
Date: 12. 03. 2008.
Tested by: InterCAD
Page number:
File name: [Link]

Thema

Clamped plate with linear distributed load.
Analysis
Type
Linear analysis.
Geometry



Top view

(thickness = 22,0 cm)

Loads

Distributed load: q = -10 kN / m
2

Boundary
Conditions
eX = eY = eZ = fiX = fiY= fiZ = 0 along all edges


Material
Properties
E = 880 kN / cm
2

= 0,3

Element
types
Plate element (Heterosis)
Mesh
1 2 3
4
q
10,000
1
0
,
0
0
0 1 2 3
4
q
X
Y

AxisVM 9 Verification Examples 24


Target

m
x
, m
y

Results

Results Theory AxisVM %
m
x, 1
[kNm/m] 11,50 11,48 -0,17
m
y, 1
[kNm/m] 11,50 11,48 -0,17
m
x, 2
[kNm/m] 33,40 33,23 -0,51
m
x, 3
[kNm/m] 17,90 17,83 -0,39
m
y, 4
[kNm/m] 25,70 25,53 -0,66


Reference:

S. Timoshenko and S. Woinowsky-Krieger: Theory of Plates And Shells


AxisVM 9 Verification Examples 25

Software Release Number: R2
Date: 12. 03. 2008.
Tested by: InterCAD
Page number:
File name: [Link]

Thema

Hemisphere displacement.
Analysis
Type
Linear analysis.
Geometry



Hemisphere (Axonometric view)

t = 0,04 m

Loads

Point load P = 2,0 kN

2,0 kN
2,0 kN
A
B
C
2,0 kN
2,0 kN
A
B
C
X Y
Z



AxisVM 9 Verification Examples 26

Boundary
Conditions
eX = eY = eZ = 0 at A
eX = eY = eZ = 0 at B

Material
Properties
E = 6825 kN / cm
2

= 0,3

Element
types
Shell element
1.) guadrilateral parabolic
2.) triangle parabolic
Target

e
x
at point A
Results





e
x
[m] e

[%]
Theory 0,185
AxisVM quadrilateral 0,186 0,54
AxisVM triangle 0,185 0,00



AxisVM 9 Verification Examples 27

















































































Nonlinear static
AxisVM 9 Verification Examples 28

Software Release Number: R2
Date: 12. 03. 2008.
Tested by: InterCAD
Page number:
File name: [Link]

Thema

3D beam structure.
Analysis
Type
Geometrical nonlinear analysis.
Geometry




F
z
=-600,00 kN
F
y
=-300,00 kN
F
z
=-600,00 kN
F
y
=-300,00 kN
F
z
=-600,00 kN
Node1
Beam1
A
B
C
D
F
z
=-600,00 kN
F
y
=-300,00 kN
F
z
=-600,00 kN
F
y
=-300,00 kN
F
z
=-600,00 kN
Node1
Beam1
A
B
C
D
X
Y
Z
3,000 m
1
,
7
3
2

m
1,732 m 1,732 m
3
,
0
0
0

m
1
,
7
3
2

m
X
Y
4
,
0
0
0

m
X
Z



Loads

Py = -300 kN
Pz = -600 kN
Boundary
Conditions
eX = eY = eZ = 0 at A, B, C and D

Material
Properties
S 275
E = 21000 kN / cm
2

= 0,3
Cross-
Section
Properties
HEA 300
Ax = 112.56 cm
2
; Ix = 85.3 cm
4
; Iy = 18268.0 cm
4
; Iz = 6309.6 cm
4


Element
types
Beam
Target

eX, eY, eZ, at Node1
Nx, Vy, Vz, Tx, My, Mz of Beam1 at Node1











AxisVM 9 Verification Examples 29

Results

Comparison with the results obtained using Nastran V4

Component
Nastran

AxisVM %
eX [mm] 17,898 17,881 -0,09
eY [mm] -75,702 -75,663 -0,05
eZ [mm] -42,623 -42,597 -0,06
Nx [kN] -283,15 -283,25 0,04
Vy [kN] -28,09 -28,10 0,04
Vx [kN] -106,57 -106,48 -0,08
Tx [kNm] -4,57 -4,57 0,00
My [kNm] -519,00 -518,74 -0,05
Mz [kNm] 148,94 148,91 -0,02




AxisVM 9 Verification Examples 30

Software Release Number: R2
Date: 12. 03. 2008.
Tested by: InterCAD
Page number:
File name: [Link]

Thema

Plate with fixed end and bending moment.
Analysis
Type
Geometrical nonlinear analysis.
Geometry




Edge1
Edge2
12,0 m
Edge1
Edge2
1
,
0

m
X
Y
Z



Loads Mz = 2600 kNm (2x1300 Nm) acting on Edge2
Boundary
Conditions
eX = eY = eZ = fiX = fiY = fiZ = 0 along Edge1

Material
Properties
E = 20000 N / mm
2

= 0
Cross
Section
Properties
Plate thickness: 150 mm
Rib on Edge2: circular D = 500 mm (for distributing load to the mid-side-node)

Element
types
Parabolic quadrilateral shell (heterosis)
Rib on Edge2 for distributing load to the mid-side-node
AxisVM 9 Verification Examples 31


Target

Z
at Edge2
Results


Edge1
Edge2
5,5502 rad
12,0 m
Edge1
Edge2
5,5502 rad
1
,
0

m
X
Y
Z



Theoretical results based on the differential equation of the flexible beam:

rad
Nm M
m
m N E
b a
I
E I
M
E I
M
z
plate
plate
plate
plate plate
plate
z
plate z
plate plate
5467 . 5
10 2 10 8125 . 2
12 10 6 . 2
10 6 . 2
12
10 2
10 8125 . 2
12
15 . 0 1
12
10 4
6
6
2 10
4
3 3
=


=
=
=
=
=

= =
=

=
=

l
l
l



Comparison the AxisVM result with the theoretical one:

Component Theory AxisVM %
fiZ [rad] 5,5467 5,5502 0,06



AxisVM 9 Verification Examples 32































BLANK








AxisVM 9 Verification Examples 33
















































































Dynamic
AxisVM 9 Verification Examples 34

Software Release Number: R2
Date: 12. 03. 2008.
Tested by: InterCAD
Page number:
File name: [Link]

Thema

Deep simply supported beam.
Analysis
Type
Dynamic analysis.
Geometry





Beam (Axonometric view)
Cross section (square 2,0 m x 2,0 m)


Loads

Self-weight
Boundary
Conditions
eX = eY = eZ = fiX = 0 at A
eY = eZ = 0 at B

Material
Properties
E = 20000 kN / cm
2

= 0,3
= 8000 kg / m
3


Element
types
Three node beam element (shear deformation is taken into account)
Target

First 7 mode shapes
AxisVM 9 Verification Examples 35

Results

Mode 1: f = 43,16 Hz

Mode 3: f = 124,01 Hz

Mode 5: f = 152,50 Hz

Mode 7: f = 293,55 Hz













Mode 2: f = 43,16 Hz

Mode 4: f = 152,50 Hz

Mode 6: f = 293,55 Hz


AxisVM 9 Verification Examples 36


Results

Comparison with NAFEMS example

Mode NAFEMS (Hz) AxisVM (Hz) %
1 42,65 43,16 -1,20
2 42,65 43,16 -1,20
3 125,00 124,01 0,79
4 148,31 152,50 -2,83
5 148,31 152,50 -2,83
6 284,55 293,55 -3,16
7 284,55 293,55 -3,16



AxisVM 9 Verification Examples 37

Software Release Number: R2
Date: 12. 03. 2008.
Tested by: InterCAD
Page number:
File name: [Link]

Thema

Clamped thin rhombic plate.
Analysis
Type
Dynamic analysis.
Geometry




Top view of plane

(thickness = 5,0 cm)

Loads

Self-weight
Boundary
Conditions
eX = eY = fiZ = 0 at all nodes (ie: eX, eY, fiZ constained at all nodes)
eZ = fiX = fiY = 0 along the 4 edges

Material
Properties
E = 20000 kN / cm
2

= 0,3
= 8000 kg / m
3


Element
types
Parabolic quadrilateral shell element (heterosis)
Mesh
10,000
1
0
,
0
0
0
X
Y


AxisVM 9 Verification Examples 38


Target

First 6 mode shapes
Results
eR
0,506
0,470
0,433
0,397
0,361
0,325
0,289
0,253
0,217
0,181
0,144
0,108
0,072
0,036
0


Mode 1: f = 8,02 Hz
eR
0,486
0,451
0,416
0,382
0,347
0,312
0,278
0,243
0,208
0,174
0,139
0,104
0,069
0,035
0


Mode 3: f = 18,41 Hz
eR
0,498
0,462
0,427
0,391
0,356
0,320
0,284
0,249
0,213
0,178
0,142
0,107
0,071
0,036
0


Mode 5: f = 24,62 Hz

eR
0,463
0,429
0,396
0,363
0,330
0,297
0,264
0,231
0,198
0,165
0,132
0,099
0,066
0,033
0


Mode 2: f = 13,02 Hz
eR
0,520
0,483
0,446
0,409
0,372
0,335
0,297
0,260
0,223
0,186
0,149
0,112
0,074
0,037
0


Mode 4: f = 19,33 Hz
eR
0,449
0,417
0,385
0,353
0,321
0,289
0,257
0,225
0,192
0,160
0,128
0,096
0,064
0,032
0


Mode 6: f = 28,24 Hz


Results

Comparison with NAFEMS example

Mode NAFEMS (Hz) AxisVM (Hz) %
1 7,94 8,02 1,01
2 12,84 13,02 1,40
3 17,94 18,41 2,62
4 19,13 19,33 1,05
5 24,01 24,62 2,54
6 27,92 28,24 1,15



AxisVM 9 Verification Examples 39

Software Release Number: R2
Date: 12. 03. 2008.
Tested by: InterCAD
Page number:
File name: [Link]

Thema

Cantilevered thin square plate.
Analysis
Type
Dynamic analysis.
Geometry



Top view
(thickness = 5,0 cm)

Loads

Self-weight
Boundary
Conditions
eX = eY = eZ = fiX = fiY = fiZ = 0 along y-axis

Material
Properties
E = 20000 kN / cm
2

= 0,3
= 8000 kg / m
3


Element
types
Parabolic quadrilateral shell element (heterosis).
Mesh

AxisVM 9 Verification Examples 40


Target

First 5 mode shapes
Results


Mode 1: f = 0,42 Hz


Mode 3: f = 2,53 Hz


Mode 5: f = 3,68 Hz













AxisVM 9 Verification Examples 41



Mode 2: f = 1,02 Hz


Mode 4: f = 3,22 Hz




Comparison with NAFEMS example

Mode NAFEMS (Hz) AxisVM (Hz) %
1 0,421 0,420 -0,24
2 1,029 1,020 -0,87
3 2,580 2,530 -1,94
4 3,310 3,220 -2,72
5 3,750 3,680 -1,87



AxisVM 9 Verification Examples 42

Software Release Number: R2
Date: 12. 03. 2008.
Tested by: InterCAD
Page number:
File name: [Link]

Thema

Cantilevered tapered membrane.
Analysis
Type
Dynamic analysis.
Geometry



Side view
(thickness = 10,0 cm)

Loads

Self-weight
Boundary
Conditions
eZ = 0 at all nodes (ie: eZ constained at all nodes)
eX = eY = 0 along y-axis

Material
Properties
E = 20000 kN / cm
2

= 0,3
= 8000 kg / m
3


Element
types
Parabolic quadrilateral membrane (plane stress)
Mesh
10,000
1
,
0
0
0
5
,
0
0
0
X
Y

AxisVM 9 Verification Examples 43


Target

First 4 mode shapes
Results


Mode 1: f = 44,33 Hz
10,000
1
,
0
0
0
5
,
0
0
0
X
Y


Mode 2: f = 128,36 Hz









10,000
1,000
5,000
X
Y
AxisVM 9 Verification Examples 44




10,000
1
,
0
0
0
5
,
0
0
0
X
Y

Mode 3: f = 162,48 Hz
10,000
1
,
0
0
0
5
,
0
0
0
X
Y

Mode 4: f = 241,22 Hz

Results

Comparison with NAFEMS example

Mode NAFEMS (Hz) AxisVM (Hz) %
1 44,62 44,33 -0,65
2 130,03 128,36 -1,28
3 162,70 162,48 -0,14
4 246,05 241,22 -1,96


AxisVM 9 Verification Examples 45

Software Release Number: R2
Date: 12. 03. 2008.
Tested by: InterCAD
Page number:
File name: [Link]

Thema

Flat grillages.
Analysis
Type
Dynamic analysis.
Geometry




Top view

Loads

Self-weight
Boundary
Conditions
eX = eY = eZ = 0 at the ends (simple supported beams)

Material
Properties
E = 20000 kN / cm
2

G = 7690 kN / cm
2

= 0,3
= 7860 kg / m
3


Cross
Section
A = 0,004 m
2

Ix = 2,5E-5 m
4

Iy = Iz = 1,25E-5 m
4

Element
types
Three node beam element (shear deformation is taken into account)
Mesh
1,500 1,500 1,500
1
,
0
0
0
2
,
0
0
0
0,500 1,000 4,500
1,000 0,500
X
Y

AxisVM 9 Verification Examples 46


Target

First 3 mode shapes
Results
1
,
8
7
9
1
,
5
8
6
1
,
2
4
1
1
,
1
1
4
1
,
6
0
5
1
,
6
3
8
1
,
0
3
5
1
,
6
7
9
X
Y Z

Mode 1: f = 16,90 Hz
-
2
,
0
6
5
0
,
8
5
6
2
,
2
5
4
2
,
0
4
0
-
1
,
8
1
3
1
,
9
3
8
-
1
,
8
3
7
X
Y Z

Mode 2: f = 20,64 Hz
-
1
,
6
2
0
1
,
5
8
5
-
1
,
9
9
2
-
1
,
6
6
7
-
1
,
5
8
1
1
,
7
2
1
-
1
,
8
4
5
-
1
,
1
3
0
2
,
0
4
0
X
Y Z


Mode 3: f = 51,76 Hz

AxisVM 9 Verification Examples 47




Mode Reference AxisVM (Hz) %
1 16,85 16,90 0,30
2 20,21 20,64 2,13
3 53,30 51,76 -2,89


Reference:

C.T.F. ROSS: Finite Element Methods In Engineering Science


AxisVM 9 Verification Examples 48

































BLANK
AxisVM 9 Verification Examples 49

















































































Stability
AxisVM 9 Verification Examples 50

Software Release Number: R2
Date: 12. 03. 2008.
Tested by: InterCAD
Page number:
File name: [Link]

Thema

Simply supported beam.
Analysis
Type
Buckling analysis.
Geometry



Front view


1
2 3
4 5
6
7
8
9
S
G
1
2
10,0
1
,
0
2
0
,
0
y
z

Cross section (I
z
=168,3 cm
4
, I
t
=12,18 cm
4
, I
w
=16667 cm
6
)

Loads

Bending moment at both ends of beam
M
A
= 1,0 kNm, M
B
= -1,0 kNm
Boundary
Conditions
eX = eY = eZ = 0 at A
eX = eY = eZ = 0 at B
k
z
= k
w
= 1
Material
Properties
E = 20600 kN / cm
2

= 0,3
G = 7923 kg / m
2

Element
types
Parabolic quadrilateral shell element (heterosis)
Mesh

AxisVM 9 Verification Examples 51


Target

M
cr
= ? (for lateral torsional buckling)
Results


Analytical solution

Z
t
Z
W Z
cr
I E
I G L
I
I
L
I E
M


+

=
2
2
2
2

kNm kNcm M
cr
51 , 124 12451
3 , 168 20600
18 , 12 7923 200
3 , 168
16667
200
3 , 168 20600
2
2
2
2
= =


+

=

AxisVM result

M
cr
= 125,3 kNm

Difference

+0,6%


AxisVM 9 Verification Examples 52

Software Release Number: R2
Date: 12. 03. 2008.
Tested by: InterCAD
Page number:
File name: [Link]

Thema

Simply supported beam.
Analysis
Type
Buckling analysis.
Geometry





Front view (L = 1,0 m)



1
2
3 4
5
S
G
1
2
12,0
1
0
,
0
y
z

1
2
3 4
5
S
G
2
30,0
1
0
,
0
y
z
Section A
1
Section A
2

Cross-sections

Loads

P = -1,0 N at point B.

Boundary
Conditions
eX = eY = eZ = 0 at A
eY = eZ = 0 at B

Material
Properties
E = 20000 kN / cm
2

= 0,3

Element
types
Beam element
Target

P
cr
= ? (for inplane buckling)
Results


Theory AxisVM e [%]
P
cr
[N] 3340,0 3337,0 -0,09




AxisVM 9 Verification Examples 53

















































































Design
AxisVM 9 Verification Examples 54

Software Release Number: R2
Date: 12. 03. 2008.
Tested by: InterCAD
Page number:
File name: RC [Link]

Thema

N-M interaction curve of cross-section (EN 1992-1-1:2004).
Analysis
Type
Linear static analysis+design.
Geometry


220








328


Section: 300x400 mm
Covering: 40 mm

Loads


Boundary
Conditions

Material
Properties
Concrete:
f
cd
=14,2 N/mm
2

e
c1
=0,0007 e
cu
=0,0035 (rectangular - diagram)
Steel:
f
sd
=348 N/mm
2

e
su
=0,015

Target

Compare the program results with with hand calculation at keypoints of M-N interaction
curve.
Results

N
1


2

6


5 3

4









Reference: Dr. Kollr L. P., Vasbetonszerkezetek I. Megyetemi kiad

N [kN] M [kNm] M(N) AxisVM e%
1 -2561 +61 +61,4 +0,7
2 -1221 +211 +212,3 +0,6
3 0 +70 +70,5 +0,7
4 +861 -61 -61,4 +0,7
5 0 -190 -191,2 +0,6
6 -362 +211 +212,3 +0,6

AxisVM 9 Verification Examples 55

Software Release Number: R2
Date: 12. 03. 2008.
Tested by: InterCAD
Page number:
File name: [Link]

Thema

RC beam deflection according to EC2, EN 1992-1-1:2004.
Analysis
Type
Material nonlinear analysis.
Geometry




q = 17 kN/m




L = 5,60 m

Side view




220

35 cm covering = 3 cm
= 0,5
420


25 cm

Section

Loads

q = 17 kN /m distributed load
Boundary
Conditions
Simply supported beam.

Material
Properties
Concrete: C25/30, = 2,1
Steel: B500B

Element
types
Parabolic quadrilateral plate element (Heterosis)
Target

e
z, max















AxisVM 9 Verification Examples 56


Results

-
0
,
0
0
2
-
0
,
0
0
2
-
0
,
0
0
2
-
0
,
0
0
2
-
5
,
2
3
9
-
5
,
2
3
9
-
1
0
,
1
0
1
-
1
0
,
1
0
1
-
1
4
,
2
4
2
-
1
4
,
2
4
2
-
1
7
,
3
9
3
-
1
7
,
3
9
3
-
1
9
,
3
6
0
-
1
9
,
3
6
0
-
2
0
,
0
2
9
-
2
0
,
0
2
9
-
1
9
,
3
6
0
-
1
9
,
3
6
0
-
1
7
,
3
9
3
-
1
7
,
3
9
3
-
1
4
,
2
4
2
-
1
4
,
2
4
2
-
1
0
,
1
0
1
-
1
0
,
1
0
1
-
5
,
2
3
9
-
5
,
2
3
9
-
5
,
2
3
9
-
5
,
2
3
9
-
1
0
,
1
0
1
-
1
0
,
1
0
1
-
1
4
,
2
4
2
-
1
4
,
2
4
2
-
1
7
,
3
9
3
-
1
7
,
3
9
3
-
1
9
,
3
6
0
-
1
9
,
3
6
0
-
2
0
,
0
2
9
-
2
0
,
0
2
9
-
2
0
,
0
2
9
-
2
0
,
0
2
9
-
1
9
,
3
6
0
-
1
9
,
3
6
0
-
1
7
,
3
9
3
-
1
7
,
3
9
3
-
1
4
,
2
4
2
-
1
4
,
2
4
2
-
1
0
,
1
0
1
-
1
0
,
1
0
1
-
5
,
2
3
9
-
5
,
2
3
9
-
5
,
2
3
9
-
5
,
2
3
9
-
1
0
,
1
0
1
-
1
0
,
1
0
1
-
1
4
,
2
4
2
-
1
4
,
2
4
2
-
1
7
,
3
9
3
-
1
7
,
3
9
3
-
1
9
,
3
6
0
-
1
9
,
3
6
0
-
2
0
,
0
2
9
-
2
0
,
0
2
9
-
2
0
,
0
2
9
-
2
0
,
0
2
9
-
1
9
,
3
6
0
-
1
9
,
3
6
0
-
1
7
,
3
9
3
-
1
7
,
3
9
3
-
1
4
,
2
4
2
-
1
4
,
2
4
2
-
1
0
,
1
0
1
-
1
0
,
1
0
1
-
5
,
2
3
9
-
5
,
2
3
9
X
Z

Diagram e
z



Aproximate calculation:

mm e e e
I II
_ 06 , 20 ) 1 ( = + =
where,
e
I
is the deflection which was calculated with the uncracked inertia moment
e
II
is the deflection which was calculated with the cracked inertia moment

2
1
|
|

\
|
=
s
sr

Calculation with integral of :

e = 19,82 mm

Calculation with AxisVM:

e = 20,03 mm (different +1,1%)



AxisVM 9 Verification Examples 57

Software Release Number: R2
Date: 12. 03. 2008.
Tested by: InterCAD
Page number:
File name: [Link]

Thema

Required steel reinforcement of RC plate according to EC2, EN 1992-1-1:2004.
Analysis
Type
Linear analysis.
Geometry

50 kN 50 kN
4,0
X
Y
Szabvny : Eurocode
Eset : ST1

Side view

Cross-section

Loads

Pz = -50 kN point load
Boundary
Conditions
Clamped cantilever plate.

Material
Properties
Concrete: C25/30
Steel: B500A
Element
types
Parabolic quadrilateral plate element (heterosis)
Mesh

Clamped
edge
4,0
1
,
0
Clamped
edge
X
Y
Szabvny : Eurocode
Eset : ST1

Top view

AxisVM 9 Verification Examples 58

Target

A
XT
steel reinforcement along x direction at the top of the support

Results


Clamped
edge
ST1, axf: 2093 mm
2
/m
4
,
0
1
,
0
Clamped
edge
ST1, axf: 2093 mm
2
/m
X Y
Z
Lineris szmts
Szabvny : Eurocode
Eset : ST1
E (W) : 1,09E-11
E (P) : 1,09E-11
E (ER) : 8,49E-13
Komp. : axf [mm2 /m]

Diagram A
XT



Calculation according to EC2:

2
/ 6 , 16
5 , 1
25
mm N f
cd
= =

2
/ 435
15 , 1
500
mm N f
yd
= =


54 , 0
435 20000 0035 , 0
20000 0035 , 0 85 , 0
0
=
+

=
+

=
yd S cu
S cu
c
f E
E c

d = 300 53 = 247 mm


kNm
x
d f x b M M
c
cd c Rd sd
200
2
=
|

\
|
= =

)
`

>
=
55
439 h
x
c


54 , 0 22 , 0
247
55
0
= < = = =
c
c
c
d
x

Steel reinforcement is yielding

2
2099
435
6 , 16 1000 55
mm
f
f x b
A
yd
cd c
S
=

=

=


Calculation with AxisVM:

A
XT
= m mm / 2093
2


Different = -0,3 %


AxisVM 9 Verification Examples 59

Software Release Number: R2
Date: 12. 03. 2008.
Tested by: InterCAD
Page number:
File name: [Link]

Thema

Earth-quake design using response-spectrum method.
Analysis
Type
Linear frequency analysis with 5 modes.
Linear static analysis.
Geometry


6
,
0
0
0
7,000
8
,0
0
0 30,0
5
,
1
9
6
5
,0
0
0
90,0
90,0
X
Y
Code : Eurocode
Case : FR +


Top view



4
,
0
0
0
3
,
5
0
0
X
Z
Code : Eurocode
Case : FR +


Front view

AxisVM 9 Verification Examples 60

All columns 60x40 cm
All nodal masses are Mx=My=Mz=100000 kg
Column A
Column B
All beams 60x40 cm
Inertia about vertical axis is multiplied by 1000.
All supports are constrained in all directions.
eX=eY=eZ=fiX=fiY=fiZ=0
Support C
Node D
All columns 60x40 cm
All nodal masses are Mx=My=Mz=100000 kg
Column A
Column B
All beams 60x40 cm
Inertia about vertical axis is multiplied by 1000.
All supports are constrained in all directions.
eX=eY=eZ=fiX=fiY=fiZ=0
Support C
Node D
X
Y
Z
Code : Eurocode
Case : ST1

Perspective view

Section beams: 60x40 cm
Ax=2400 cm2 Ay=2000 cm2 Az=2000 cm2
Ix=751200 cm4 Iy=720000 cm2 Iz=320000000 cm4

Section columns: 60x40 cm
Ax=2400 cm2 Ay=2000 cm2 Az=2000 cm2
Ix=751200 cm4 Iy=720000 cm2 Iz=320000 cm4

Loads

Nodal masses on eight nodes. Mx=My=Mz=100000 kg
Model self-weight is excluded.

Spectrum for X and Y direction of seismic action:
S
d
[m/s
2
]
T[s]

1,150
2,156
0,300
2,0000
0,709
T[s] S
d
1 0 1,150
2 0,2000 2,156
3 0,6000 2,156
4 1,3000 0,995
5 3,0000 0,300
6 4,0000 0,300
. . . . . .


Boundary
Conditions
Nodes at the columns bottom ends are constrained in all directions.
eX=eY=eZ=fiX=fiY=fiZ=0

Material
Properties
C25/30
E=3050 kN/cm2 =0,2 = 0

AxisVM 9 Verification Examples 61


Element
types
Three node straight prismatic beam element. Shear deformation is taken into account.
Target

Compare the model results with SAP2000 v6.13 results.
The results are combined for all modes and all direction of spectral acceleration.
CQC combination are used for modes in each direction of acceleration.
SRSS combination are used for combination of directions.

Results

Period times of first 5 modes
Mode T[s] SAP2000 T[s] AxisVM Difference [%]
1 0,7450 0,7450 0
2 0,7099 0,7099 0
3 0,3601 0,3601 0
4 0,2314 0,2314 0
5 0,2054 0,2054 0

Modal participating mass ratios in X and Y directions
Mode X
SAP2000
X
AxisVM
Difference
%
Y
SAP2000
Y
AxisVM
Difference
%
1 0,5719 0,5719 0 0,3153 0,3154 +0,03
2 0,3650 0,3650 0 0,4761 0,4760 -0,02
3 0 0 0 0,1261 0,1261 0
4 0,0460 0,0460 0 0,0131 0,0131 0
5 0,0170 0,0170 0 0,0562 0,0562 0
Summ 1,0000 1,0000 0 0,9868 0,9868 0

Internal forces at the bottom end of Column A and Column B
Column A
SAP2000
Column A
AxisVM
Difference
%
Column B
SAP2000
Column B
AxisVM
Difference
%
Nx [kN] 315,11 315,15 +0,01 557,26 557,29 +0,005
Vy [kN] 280,34 280,34 0 232,88 232,88 0
Vz [kN] 253,49 253,49 0 412,04 412,04 0
Tx [kNm] 34,42 34,41 -0,03 34,47 34,46 -0,03
My [kNm] 625,13 625,12 -0,002 1038,74 1038,70 -0,004
Mz [kNm] 612,31 612,31 0 553,41 553,41 0

Support forces of Support C
Support C
SAP2000
Support C
AxisVM
Difference
%
Rx [kN] 280,34 280,34 0
Ry [kN] 253,49 253,49 0
Rz [kN] 315,11 315,15 +0,01
Rxx [kNm] 625,13 625,12 -0,002
Ryy [kNm] 612,31 612,31 0
Rzz [kNm] 34,42 34,41 -0,03

Displacements of Node D
Node D
SAP2000
Node D
AxisVM
Difference
%
eX [mm] 33,521 33,521 0
eY [mm] 19,944 19,945 +0,005
eZ [mm] 0,229 0,229 0
X [rad] 0,00133 0,00133 0
Y [rad] 0,00106 0,00106 0
Z [rad] 0,00257 0,00257 0


AxisVM 9 Verification Examples 62

Normal forces:




AxisVM 9 Verification Examples 63

Bending moments:




AxisVM 9 Verification Examples 64




AxisVM 9 Verification Examples 65

Displacements: