finish /clear /prep7 /eshape,0 pi=2*asin(1) !!********Beginn*Eingabebereich*********!! !! Geometrie f_laenge=150 ! freie Laenge e_laenge=400 ! Einspannlaenge breite=150 durchmesser=20 ! Kerbdurchmesser bolzen_d=19 ! Bolzendurchmesser !!Last kraft=10 !! Material Scheibe -> CFK ex=120000 ! Emodul ey=6000 ez=6000 gxy=3200 ! Schubmodul gxz=3200 gyz=2200 nuxy=0.26 ! Querkontraktionszahl nuxz=0.26 nuyz=0.26 dicke=2.1 !! Material Bolzen -> Titan ex2=105000 gx2=41400 nux2=0.34 !! Netzfeinheit teil=40 ! Teilung Radius/4! !!*********Ende*Eingabebereich**********!! mp,ex,1,ex ! Material zuweisen Scheibe mp,ey,1,ey mp,ez,1,ez mp,gxy,1,gxy mp,gxz,1,gxz mp,gyz,1,gyz mp,prxy,1,nuxy mp,prxz,1,nuxz mp,pryz,1,nuyz mp,ex,2,ex2 ! Material zuweisen Bolzen !mp,gxy,2,gx2 mp,prxy,2,nux2 et,1,solid45 ! Elementdefinition Scheibe et,2,targe170 ! Elementdefinition Bolzen r,2,bolzen_d/2 et,3,conta173 ! Kontaktelemente Scheibe keyopt,3,2,1 !keyopt,3,5,1 !keyopt,3,12,6 !keyopt,3,15,2 k,1,durchmesser/2,0,0 ! Geometrie erzeugen k,2,durchmesser,0,0 k,3,f_laenge,0,0 k,4,0,durchmesser/2,0 k,5,0,durchmesser,0 k,6,durchmesser,durchmesser,0 k,7,f_laenge,durchmesser,0 k,8,0,breite/2,0 k,9,durchmesser,breite/2,0 k,10,f_laenge,breite/2,0 k,11,durchmesser/2,,-dicke ! Horizontale Linien l,1,2,teil,5 $ l,2,3,teil/2 $ l,5,6,teil/2 l,6,7,teil/2 $ l,8,9 $ l,9,10,teil/2 ! Vertikale Linien l,4,5,teil,5 $ l,2,6,teil/2 $ l,3,7 $ l,5,8 $ l,6,9,teil/2 $ l,7,10 l,1,11,3 ! Dickenlinie !Linie des Kerbrandes csys,1 $ l,4,1,teil $ csys,0 local,12,0 ! lokales Koordinatensystem ! Flächen erzeugen a,1,2,6,5,4 $ a,2,3,7,6 $ a,5,6,9,8 $ a,6,7,10,9 ! Viertelscheibe zur Vollscheibe spiegeln arsym,x,all $ arsym,y,all nummrg,kp ! doppelte Knoten verschmelzen ! Prüfkörpergeometrie erzeugen k,101,e_laenge,-breite/2 k,102,e_laenge,-durchmesser k,103,e_laenge k,104,e_laenge,durchmesser k,105,e_laenge,breite/2 l,101,102 $ l,102,103 $ l,103,104 l,104,105 $ l,31,101,teil/2 $ l,28,102 l,3,103,teil/2 $ l,7,104 $ l,10,105,teil/2 a,31,101,102,28 $ a,28,102,103,3 a,3,103,104,7 $ a,7,104,105,10 ! Volumen über Dickenlinie austragen vdrag,1,2,3,4,5,6,13 vdrag,7,8,9,10,11,12,13 vdrag,13,14,15,16,17,18,13 vdrag,19,20,,,,,13 nummrg,kp ! Flächen des Kerbrandes für mapped ! mesh verbinden accat,58,59 $ accat,22,23 accat,39,40 $ accat,74,75 allsel vatt,1,1,1,12,1 allsel $ type,1 $ real,1 $ mat,1 mshape,0,3d mshkey,1 vmesh,all ! Scheibe vernetzen !! Bolzen hinzufuegen k,201,,, k,202,,,-dicke k,203,1 circle,201,bolzen_d/2,202,203,360,4 a,201,15,45 a,201,45,46 a,201,46,48 a,201,48,15 l,201,202 vdrag,64,73,76,79,,,132 type,1 real,1 mat,2 ! Material Bolzen vsel,s,,,21,24 ! Bolzen vernetzen mshape,0,3d mshkey,1 vmesh,all type,2 real,2 mat,1 ! Material Kontakt asel,s,,,93 ! Kontaktelemente auf Bolzen asel,a,,,102 asel,a,,,109 asel,a,,,112 areverse,all,0 amesh,all type,3 $ real,2 ! Kontaktelemente Kerbrand Scheibe asel,s,area,,25 $ asel,a,area,,38 asel,a,area,,56 $ asel,a,area,,77 amesh,all ! Kontaktelemete am areverse,all,0 ! Kerbrand aufbringen alls csys,0 ! Scheibe einspannen nsel,s,loc,x,e_laenge d,all,uy,0 d,all,uz,0 d,all,ux,-((durchmesser-bolzen_d)/2)-0.01 /solu nlgeom,on $ rescontrol,,all,1,20 outres,all,all $ nsubst,5,100,2 allsel,all ! Berechnung des solve ! ersten Lastschrittes nlgeom,on autots,on rescontrol,,all,1,20 outres,all,all nsubst,5,10,5 ! Scheibeneinspannung lösen nsel,s,loc,x,e_laenge ddele,all,all d,all,uy,0 ! Scheibe einspannen d,all,uz,0 *get,knoanz,node,all,count ! alle Knoten am Scheibenende zählen f,all,fx,kraft/knoanz ! und Kraft auf die Scheibe aufbringen !ncnv,,1e20 nsel,s,loc,x,0 nsel,r,loc,y,0 d,all,all,0 allsel,all solve ! Zweiter Lastschritt /eof /post1 ! Ausgabe der Knoten- rsys,0 ! verschiebung shell,mid plnsol,u,x, 0,1.0