Code:
/prep7
et, 1, 42 ! 4 Element Flächenelement
et, 2, link1 ! Stabelement
et, 3, beam3 ! Balken
et, 4, 45 ! Volumenelement
mp, ex, 1, 2.10e5 ! E-Modul Stahl
mp, ex, 2, 0.75e5 ! E-Modul Aluminium
mp, ex, 3, 0.10e5 ! E-Modul Kunststoff
mp, nuxy, 3, 0.45 ! Querkontraktionszahl Kunststoff
pi = 3.141592 ! Pi
d1 = 40 ! Durchmesser 40
d2 = 30 ! Durchmesser 30
d3 = 60 ! Durchmesser 60
d4 = 50 ! Durchmesser 50
a1 = (d1/2)**2 * pi ! Querschnittsfläche 1
a2 = (d2/2)**2 * pi ! Querschnittsfläche 2
a3 = (d3/2)**2 * pi ! Querschnittsfläche 3
a4 = (d4/2)**2 * pi ! Querschnittsfläche 4
i1 = pi * ((d1/2)**4) / 4 ! Trägheitsmoment 1
i2 = pi * ((d2/2)**4) / 4 ! Trägheitsmoment 2
i3 = pi * ((d3/2)**4) / 4 ! Trägheitsmoment 3
i4 = pi * ((d4/2)**4) / 4 ! Trägheitsmoment 4
r, 1, a1, i1, d1
r, 2, a2, i2, d2
r, 3, a3, i3, d3
r, 4, a4, i4, d4
a5 = (0.5e5*50) / (2.10e5) ! Berechnung Fläche des Federstabelements: A = (c * l) / E
r, 5, a5 ! Federsteifigkeit [N/mm]
! BALKEN 1:
type, 3 ! Balkenelemente auswählen
mat, 1 ! Stahl Materialgrößen auswählen
! Knoten:
n, 6000, 0, 50, 25
n, 6001, 40, 50, 25
n, 6002, 60, 50, 25
n, 6003, 660, 50, 25
n, 6004, 685, 50, 25
! Balkenelemente:
real, 1
e, 6000, 6001
real, 2
e, 6001, 6002
real, 3
e, 6002, 6003
real, 2
e, 6003, 6004
! BALKEN 2:
type, 3 ! Balkenelemente auswählen
mat, 2 ! Aluminium Materialgrößen auswählen
! Knoten:
n, 6005, 685, 50, 25 ! Selber Knoten wie 6004 für Gelenk
n, 6006, 710, 50, 25
n, 6007, 760, 50, 25
n, 6008, 860, 50, 25
! Balkenelemente:
real, 2
e, 6005, 6006
real, 4
e, 6006, 6007
real, 3
e, 6007, 6008
! FEDERSTAB:
type, 2 ! Stabelemente auswählen
real, 5 ! Querschnittsfläche auswählen
mat, 1 ! Stahl auswählen
n, 6009, 685, 0, 25
e, 6009, 6005
! GELENK:
cp, 1, ux, 6004, 6005
cp, 2, uz, 6004, 6005
cp, 3, uy, 6004, 6005
! LAGER:
! Einspannung rechte Seite
d, 6008, all, 0
! Federlager unten fest einspannen
d, 6009, all, 0
! Einspannung linke Seite
d, 6000, all, 0
! KRÄFTE:
/psf, pres, , 2
sfbeam, 642, 1, pres, 200
sfbeam, 646, 1, pres, 150
/pbc,all,1
/solu ! Lösen
solve