/title, Berechnung der mech. Spannungen am Schneidkeil (2D)


!	****************
!	Modelldefinition
!	****************

/prep7				! Preprozessor starten
shpp,modify,11,100		! Warnung erst bei Parallelenabweichung >100°

!  ---- Definition der Geometrie- und Lastparameter

Mikwin = 3			! Mikrowinkel in Grad
Fase   = 0.1			! Fase in mm
KraftY = 10			! Passivkraft in N
KraftX = 500			! Schnittkraft in N
eingr  = 2			! Eingriffstiefe in mm
b      = 2			! Schnittbreite in mm


!
!  ---- Werkstoffeigenschaften ----
!

mp, ex,   1, 2.11e5		! Elastizitätsmodul in N/mm² (Diamant)
mp, dens, 1, 4.12e-6		! Dichte in Kg/mm³  (Diamant)
mp, nuxy, 1, 0.3   		! Querkontraktionszahl

mp, ex,   2, 8.0e3		! Elastizitätsmodul in N/mm² (Holz)
mp, dens, 2, 0.8e-6		! Dichte in Kg/mm³   (Holz)
mp, nuxy, 2, 0.3   		! Querkontraktionszahl  (Holz)

!
!  ---- Elementtypen ----
!	
et, 1,182,,,3			! Plane 182	(2-D Structural Solid) ebener Spannungszustand
r,1,b				! Realkonstante für Elementtyp Plane 182: Schnittbreite b

et, 2,169			! Target 169	(2-D Target Segment)
et, 3,171			! Contact 171	(2-D Surface-to-Surface Contact)

keyopt,3,5,3			! Close gap/reduce penetration with auto Contact surface offset (CNOF)
keyopt,3,10,4			! Änderung der Kontaktsteifigkeit bei jedem Substep möglich

fkn=0.3				! Skalierungsfakor für die Kontaktsteifigkeit
pmin=0				! Lower limit of initial allowable penetration
pmax=pmin			! Upper limit of initial allowable penetration

r,2,,,fkn			! Realkonstanten für Kontakt zwischen Holz und Diamand
rmore,pmax,pmin

!
!  ---- Geometrie ----
!
!  Schneidkeil

imme,off			! Immediate Printing abstellen

K,1,0,0
K,2,Fase,TAN(Mikwin*3.1416/180)*Fase,
K,3,Fase-(TAN(15*3.1416/180)*(eingr-TAN(Mikwin*3.1416/180)*Fase)),eingr
K,4,Fase-(Sin(15*3.1416/180)*15),(TAN(Mikwin*3.1416/180)*Fase)+COS(15*3.1416/180)*15
K,5,-5.5235,13.6892
K,6,-1.9,0.1662

LSTR,1,2,
LSTR,2,3,
LSTR,3,4,
LSTR,4,5,
LSTR,5,6,
LSTR,6,1,

al,1,2,3,4,5,6			! Fläche 1 (Schneidkeil)

esize,0.5				! Default-Elementlänge 0.5 mm
lesize,1,0.005			! Elementlänge an der Linie 1: 0.005 mm
lesize,2,0.015,,,10		! Elementlänge an der Linie 2: 0.015 mm, Längenverhältnis: 1:10
lesize,6,0.015,,,0.1		! Elementlänge an der Linie 6: 0.015 mm, Längenverhältnis: 10:1
lesize,3,0.15,,,10		! Elementlänge an der Linie 3: 0.15 mm, Längenverhältnis: 1:10
lesize,5,0.15,,,0.1		! Elementlänge an der Linie 5: 0.15 mm, Längenverhältnis: 10:1

mat,1				! Werkstoff 1 = Diamant
asel,s,,,1 $ amesh,all		! Fläche 1 vernetzen (Elementtyp 1, Realkonstantenset 1)

alls
imme,on				! Immediate Printing einschalten
eplo				! Elemente (Netz) darstellen

C***
C***	Werkstoffgenerierung und Vernetzung
C***

!
!  ---- Geometrie ----
!
!  Werkstoff

K,7,0,0
K,8,Fase+0.01,(TAN(Mikwin*3.1416/180)*Fase)-0.01,
K,9,(Fase-(TAN(15*3.1416/180)*(eingr-TAN(Mikwin*3.1416/180)*Fase)))+0.1,eingr
K,10,5,eingr
K,11,5,-2
K,12,-5,-2
K,13,-5,0

LSTR,7,8,
LSTR,8,9,
LSTR,9,10,
LSTR,10,11
LSTR,11,12,
LSTR,12,13,
LSTR,13,7,


al,7,8,8,9,10,11,12,13		! Fläche 2 (Werkstoff)

esize,0.5			! Default-Elementlänge 0.5 mm
lesize,7,0.005			! Elementlänge an der Linie 7: 0.005 mm (Deckungsgleich L 1)
lesize,8,0.015,,,10		! Elementlänge an der Linie 8: 0.015 mm, Längenverhältnis: 1:10 (Deckungsgleich L 2)
lesize,9,0.15,,,0.1		! Elementlänge an der Linie 9: 0.15 mm, Längenverhältnis: 10:1
lesize,10,0.15,,,0.1		! Elementlänge an der Linie 10: 0.15 mm, Längenverhältnis: 10:1
lesize,11,0.15,,,0.1		! Elementlänge an der Linie 11: 0.15 mm, Längenverhältnis: 10:1
lesize,12,0.15,,,0.1		! Elementlänge an der Linie 12: 0.15 mm, Längenverhältnis: 10:1
lesize,13,0.15,,,0.1		! Elementlänge an der Linie 13: 0.15 mm, Längenverhältnis: 10:1

mat,2				! Werkstoff 2 = Holz
asel,s,,,1 $ amesh,all		! Fläche 2 vernetzen (Elementtyp 1, Realkonstantenset 1)

alls
imme,on				! Immediate Printing einschalten
eplo				! Elemente (Netz) darstellen