fini
/clear
/prep7

! Alle Maße in mm
! --------------------------------------------------------------

pi=3.141592654
l=250				! Länge Platte (x-Länge)
b=120				! Breite Platte (y-Länge)
W=200				! Position Bohrung in X-Richtung
h=55				! Position Bohrung in Y-Richtung
d=50 				! Durchmesser Bohrung
delta=100			! Dicke der Probe
kl=5				! Kantenlänge
r=l-W+(W*0.25)		! Risslänge
!r=l-W+(W*0.3)
!r=l-W+(W*0.4)
!r=l-W+(W*0.5)
!r=l-W+(W*0.6)
!r=l-W+(W*0.7)

! Element Plane 42 (Scheibe)
! 0 - Elementkoordinatensystem ist parallel zum glob. Koord.sys.
! 0 - Zusätzliche Formen
! 3 - Scheibenbelastung mit Angabe der Scheibendicke (R,1,dicke)
! --------------------------------------------------------------
ET,1,42,0,0,3
R,1,delta

! Material der Scheibe
! --------------------------------------------------------------
!mp,ex,1,31200000	! E-Modul in x-Richtung (in kg/(s²*mm))
!mp,ey,1,31200000	! E-Modul in y-Richtung (in kg/(s²*mm))
!mp,ez,1,31200000	! E-Modul in z-Richtung (in kg/(s²*mm))
!mp,nuxy,1,0.3		! Querkontraktionszahl in x-y-Richtung
!mp,gxy,1,12000000	! Schubmodul in x-y-Richtung

!mp,ex,1,40950000	! E-Modul in x-Richtung (in kg/(s²*mm))
!mp,ey,1,13000000	! E-Modul in y-Richtung (in kg/(s²*mm))
!mp,ez,1,40950000	! E-Modul in z-Richtung (in kg/(s²*mm))
!mp,nuxy,1,0.3		! Querkontraktionszahl in x-y-Richtung
!mp,gxy,1,5900000	! Schubmodul in x-y-Richtung

!mp,ex,1,46100000	! E-Modul in x-Richtung (in kg/(s²*mm))
!mp,ey,1,18800000	! E-Modul in y-Richtung (in kg/(s²*mm))
!mp,ez,1,46100000	! E-Modul in z-Richtung (in kg/(s²*mm))
!mp,nuxy,1,0.3		! Querkontraktionszahl in x-y-Richtung
!mp,gxy,1,9350000	! Schubmodul in x-y-Richtung

mp,ex,1,51300000	! E-Modul in x-Richtung (in kg/(s²*mm))
mp,ey,1,24500000	! E-Modul in y-Richtung (in kg/(s²*mm))
mp,ez,1,51300000	! E-Modul in z-Richtung (in kg/(s²*mm))
mp,nuxy,1,0.3		! Querkontraktionszahl in x-y-Richtung
mp,gxy,1,8800000	! Schubmodul in x-y-Richtung

! Keypoints
! --------------------------------------------------------------
k,1,0,0
k,2,0,b
k,3,l,b
k,4,l,0
k,5,l-r,0
k,6,W-(d/2),h
k,7,W+(d/2),h

! Linien
! --------------------------------------------------------------
l,1,2
l,2,3
l,3,4
l,4,5
l,5,1

! Kreisbogen
! --------------------------------------------------------------
larc,6,7,2,(d/2)
larc,6,7,4,(d/2)

! Fläche
! --------------------------------------------------------------
AL,1,2,3,4,5
AL,6,7
ASBA,1,2

! Vernetzen
! --------------------------------------------------------------
AESIZE,3,kl
AMESH,3

! Zusatzpunkte eingeben (Mittelpunktelement)
! --------------------------------------------------------------
N,100002,l-r-(kl/4),0
N,100003,l-r+(kl/4),0
N,100004,l-r,(kl/4)
N,100005,l-r-kl,kl/2
N,100006,l-r+kl,kl/2
N,100007,l-r-(kl/2),kl
N,100008,l-r+(kl/2),kl

! Einspannungen (unbedingt Knotennummerierung anschauen!)
! ---------------------------------------------------------------
D,100,ALL
D,100002,UY,0
D,120,UY,0,,148,1 ! a/W = 0,25
!D,122,UY,0,,148,1 ! a/W = 0,3
!D,126,UY,0,,148,1 ! a/W = 0,4
!D,130,UY,0,,148,1 ! a/W = 0,5
!D,134,UY,0,,148,1 ! a/W = 0,6
!D,138,UY,0,,148,1 ! a/W = 0,7
D,1,UY,0

! Kraft (unbedingt Knotennummerierung anschauen!)
! ---------------------------------------------------------------
F,173,FY,1000

! Lösen
! ---------------------------------------------------------------
/solu
solve
/post1
pldisp,1