Code:
fini
/clear/filenam, struct ! Name der Substruktur waehrend der Analyse
/prep7
! Berechnungskonstanten
Pi=4*atan(1)
!Anzahl Balkenelemente
Ba_n = 6
!Parameter Linse
Li_Rad = 0.1
Li_h = 0.05
! Zylinder
! die Linien ebenso wie die kp 1, 2, 3, 4 liegen in der xy-Ebene
cylind,,Li_Rad,0,Li_h
lcomb,1,2
lcomb,3,4 ! es bleiben Linien 1 und 3
ldiv,1,,,Ba_n/2,0 ! aufteilen der Linie 1 in Ba_n/2 Teile
ldiv,3,,,Ba_n/2,0 ! aufteilen der Linie 2 in Ba_n/2 Teile
! Materialparameter Balken
et, 2, beam4 ! 3D elastisches Balkenelement
mp, ex, 2, 206e9 ! wikipedia Baustahl [N/m2]
mp, nuxy, 2, 0.28 ! wikipedia Baustahl Poisson Ratio [-]
mp, dens, 2, 7850 ! wikipedia Baustahl [kg/m3]
Ba_b=5e-3 ! breite [m]
Ba_h=1e-3 ! hoehe [m]
Ba_l=1e-2 ! laenge [m]
Ba_A = Ba_h*Ba_b ! Querschnittsflaeche
Ba_Iy = 1/12*(Ba_b*Ba_h*Ba_h*Ba_h) ! Flaechentragheitsmoment
Ba_Ip = Ba_b*Ba_h/12*(Ba_h*Ba_h+Ba_b*Ba_b) ! polares Flaechentraegheitsmoment
r,2,A,Iy,Iy,b,h
rmore,,Ip ! Bernoulli-Balken ohne Scherung
! Balkenlinien generieren
local,11,1 ! Lokales ZylinderKOS
ksel,s,loc,z,0 ! selektiere alle KP mit z=0
kgen,2,all,,,Ba_l, 0, 0, 100 ! generiere KP mit Abstand Ba_l
ksel,s,loc,z,0 ! selektiere alle KP mit z=0
l,1,101
l,2,102
l,3,103
l,4,104
*do,i,9,Ba_n+4
l,i,i+100
*enddo
! Balken in Elemente aufteilen
ksel,s,kp,,100,200
lslk,s
esize,,5 ! Aufteilung in die Elemente
real,2
lmesh,all
! Materialparameter Linse
!et,1,solid95
mp,ex,1,120e9
mp,ey,1,120e9
mp,ez,1,120e9
mp,nuxy,1,0.3
mp,nuyz,1,0.3
mp,nuxz,1,0.3
mp,dens,1,2210
vsel,all
mshkey,0
vmesh,all
gplot