Skip to content

Beam elements examples #20

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Open
Fonotec wants to merge 4 commits into
base: main
Choose a base branch
from
Open

Beam elements examples #20

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
4 commits
Select commit Hold shift + click to select a range
b666a31
Add the cantilever beam examples that use beam elements from Nobels (…
Fonotec Oct 22, 2025
d5092a0
Add small update to README
Fonotec Oct 22, 2025
b3be8eb
Update a comment in all the files
Fonotec Oct 22, 2025
ad40123
Update the headers of the files to have the copyright lines in the ra…
Fonotec Oct 22, 2025
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
19 changes: 19 additions & 0 deletions Components/Elements/Beam_Type3_forced_cantilever_beam/README.md
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,19 @@
# Modeling a forced cantilever beam using beam elements

This example shows a simple way to model a cantilever beam using beam elements.

This example has been taken from Nobels (2025), which can be found here: https://repository.tno.nl/SingleDoc?find=UID%20f29c6aa5-7cbd-46b2-8064-305d37841ab4

These examples are the same cantilever beam but with different numerical resolutions.
All the cantilever beams have a length of 4000 mm.
Here the following resolutions are available:
- 4000 mm
- 2000 mm
- 1000 mm
- 500 mm
- 250 mm
- 125 mm
- 62.5 mm
- 31.25 mm

These examples can be compared to analytical solutions.
169 changes: 169 additions & 0 deletions ...eam_Type3_forced_cantilever_beam/resolution_1000mm/cantileverbeam4000mmres1000mm_0000.rad
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,169 @@
#RADIOSS STARTER
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
# Cantilever beam [Beam elements] ("Model")
# Copyright (C) 2025 Nederlandse Organisatie voor toegepast-natuurwetenschappelijk onderzoek (TNO) ("Holder")
# Model is licensed by Holder under CC BY-NC 4.0
# (https://creativecommons.org/licenses/by-nc/4.0/legalcode).
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
##
##--------------------------------------------------------------------------------------------------
## UNIT section
##--------------------------------------------------------------------------------------------------
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/BEGIN
# Runname
OpenRadioss test case
# Input version
2022 0
# input mass unit | input length unit | input time unit |
kg mm ms
# work mass unit | work length unit | work time unit |
kg mm ms
##
##
##--------------------------------------------------------------------------------------------------
## Material Law No 1. Purely linearly elastic model
##--------------------------------------------------------------------------------------------------
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/MAT/LAW1/2
# material title
Steel_DP600
# density |
7.80000000000000E-06
# Young's modulus | Poisson's ratio |
210.00000000000000 0.30000000000000
##
##--------------------------------------------------------------------------------------------------
## NODES
##--------------------------------------------------------------------------------------------------
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/NODE
# node_ID| X coordinate | Y coordinate | Z coordinate |
1 0.000000 4000.000000 0.000000
2 0.000000 0.000000 0.000000
3 1000.000000 0.000000 0.000000
4 2000.000000 0.000000 0.000000
5 3000.000000 0.000000 0.000000
6 4000.000000 0.000000 0.000000
##
##--------------------------------------------------------------------------------------------------
## Boundary Conditions 1
##--------------------------------------------------------------------------------------------------
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/BCS/1
# boundary condition title
constraint1
#Trarot | skew_ID | grnd_ID |
111 111 0 2
##
##--------------------------------------------------------------------------------------------------
## Corresponding boundary nodes
##--------------------------------------------------------------------------------------------------
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/GRNOD/NODE/2
# node group for example for bounary condtions
#item_ID1| item_ID2| item_ID3| item_ID4| item_ID5| item_ID6| item_ID7| item_ID8| item_ID9|item_ID10|
grnodnode
2
##
##--------------------------------------------------------------------------------------------------
## Boundary Conditions 2
##--------------------------------------------------------------------------------------------------
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/BCS/2
# boundary condition title
constraint2
#Trarot | skew_ID | grnd_ID |
000 000 0 3
##
##--------------------------------------------------------------------------------------------------
## Corresponding boundary nodes
##--------------------------------------------------------------------------------------------------
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/GRNOD/NODE/3
# node group for example for bounary condtions
#item_ID1| item_ID2| item_ID3| item_ID4| item_ID5| item_ID6| item_ID7| item_ID8| item_ID9|item_ID10|
grnodnode
6
##
##--------------------------------------------------------------------------------------------------
## Beam Elements
##--------------------------------------------------------------------------------------------------
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/BEAM/3
# beam ID| node_ID1| node_ID2| node_ID3|
1 2 3 1
2 3 4 1
3 4 5 1
4 5 6 1
##
##--------------------------------------------------------------------------------------------------
## Beam Property Set (pid 1)
##--------------------------------------------------------------------------------------------------
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/PROP/BEAM/1
line_1.7
# | I_SMSTR |
4
#beam membrane damp|beam flexural damp |
0.0000000000000 0.0100000000000
# Area | I_YY | I_ZZ | I_XX |
64.0000000000000 341.3333333333333 341.3333333333333 42666.6666666666642
#OmegaDOF| I_shear|
000 000 0
##
##--------------------------------------------------------------------------------------------------
## Function 1
##--------------------------------------------------------------------------------------------------
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/FUNCT/1
Load
# X Y
0.0000000000000e+00 1.0000000000000e+00
1.0000000000000e+30 1.0000000000000e+00
##
##--------------------------------------------------------------------------------------------------
## Line loading for test
##--------------------------------------------------------------------------------------------------
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/CLOAD/3
# title of the imposed load
imposed_load
#fct_IDT | Dir |skew_ID | sens_ID | grnd_ID | | Ascale_x | Fscale_y |
1 Z 1 1.0000000000000000 -0.0001280000000000
##
##--------------------------------------------------------------------------------------------------
## Particles for load in the cente at the boundary
##--------------------------------------------------------------------------------------------------
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/GRNOD/NODE/1
# node group for example for bounary condtions
#item_ID1| item_ID2| item_ID3| item_ID4| item_ID5| item_ID6| item_ID7| item_ID8| item_ID9|item_ID10|
grnodnode
6
##
##--------------------------------------------------------------------------------------------------
## Parts
##--------------------------------------------------------------------------------------------------
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/PART/3
# part title
Our_part
#prop_ID | mat_ID |subset_ID| virtual thickness|
1 2
##
##--------------------------------------------------------------------------------------------------
## Input the desired information for the time file output
##--------------------------------------------------------------------------------------------------
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/TH/NODE/2
# title of nodes to follow in depth
TH_Measuring_Nodes
# var_ID1| var_ID2| var_ID3| var_ID4| var_ID5| var_ID6| var_ID7| var_ID8| var_ID9| var_ID10|
DEF XYZ A DRX DRY DRZ
# node ID| skew_ID or frame_ID
2 0
3 0
4 0
5 0
6 0
69 changes: 69 additions & 0 deletions ...eam_Type3_forced_cantilever_beam/resolution_1000mm/cantileverbeam4000mmres1000mm_0001.rad
View file
Open in desktop
Original file line number Diff line number Diff line change
@@ -0,0 +1,69 @@
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
# Cantilever beam [Beam elements] ("Model")
# Copyright (C) 2025 Nederlandse Organisatie voor toegepast-natuurwetenschappelijk onderzoek (TNO) ("Holder")
# Model is licensed by Holder under CC BY-NC 4.0
# (https://creativecommons.org/licenses/by-nc/4.0/legalcode).
#---1----|----2----|----3----|----4----|----5----|----6----|----7----|----8----|----9----|---10----|
/VERS/2017
/RUN/cantileverbeam/1/
200000.0000000000000
/H3D/NODA/VEL
/H3D/SHELL/TENS/STRESS/NPT=ALL
/H3D/SHELL/TENS/STRAIN/NPT=LOWER
/H3D/SHELL/TENS/STRAIN/NPT=UPPER
/H3D/ELEM/EPSP/NPT=UPPER
/H3D/ELEM/EPSP/NPT=LOWER
/H3D/COMPRESS
0.01
/H3D/DT
0.000000000000000 10000.000000000
/PRINT/-500
/TFILE/0
4.000000000000000
/STOP
0.000000000000000 0.000000000000000 0.000000000000000 1 1
/MON/ON
/DT/NODA/CST/0
0.450000000000000 0.000000000000000

/ANIM/DT
0.000000000000000 10000.000000000
/ANIM/VECT/VEL
/ANIM/VECT/DISP
/ANIM/VECT/ACC
/ANIM/VECT/CONT
/ANIM/VECT/CONT2
/ANIM/VECT/CONT/MAX
/ANIM/VECT/FINT
/ANIM/VECT/FEXT
/ANIM/VECT/VROT

/ANIM/NODA/DT
/ANIM/NODA/DMAS
/ANIM/NODA/DINER
/ANIM/NODA/DENS
/ANIM/NODA/P
/ANIM/NODA/TEMP
/ANIM/NODA/SSP

/ANIM/ELEM/EPSP/ALL
/ANIM/ELEM/DENS
/ANIM/ELEM/DT
/ANIM/ELEM/HOURG
/ANIM/ELEM/P
/ANIM/ELEM/SIGX
/ANIM/ELEM/SIGY
/ANIM/ELEM/SIGZ
/ANIM/ELEM/SIGXY
/ANIM/ELEM/SIGXZ
/ANIM/ELEM/SIGYZ
/ANIM/ELEM/TEMP
/ANIM/ELEM/THIC
/ANIM/ELEM/VONM
/ANIM/ELEM/MACH


/ANIM/COMPRESS
0.01

/TH/TITLE
Loading