finish ! von Präprozessor-Level auf Begin-Level wechseln /clear ! eventuell vorhandene Datenbasis löschen /NERR,0,99999999,,0,99999999 ! Programmabbruch von ANSYS nach 99999999 Warnungen/Fehlermeldungen /UIS,msgpop,3 ! display only error messages /PNUM,LINE,1 ! Einschalten der Elementnummerierung und Farbkodierung im Plot (hier für Line) /PNUM,KPOI,1 ! Einschalten der Elementnummerierung und Farbkodierung im Plot (hier für Keypoint) /PNUM,area,1 ! Einschalten der Elementnummerierung und Farbkodierung im Plot (hier für Area) /PNUM,volu,1 ! Einschalten der Elementnummerierung und Farbkodierung im Plot (hier für Volumen) /PLOPTS,INFO,auto ! iniate Auto-legend mode, turn legend column on / off /PLOPTS,LEG1,on ! Header portion of legend column on /PLOPTS,LEG2,on ! View portion of legend column on (defaults to ON (except off with contour displays)) /PLOPTS,LEG3,on ! View the contour section of the legend column on /PLOPTS,FRAME,on ! Frame border lines around windows on /PLOPTS,TITLE,on ! Title (bottom left text) on /PLOPTS,MINM,on ! Min-Max symbols on contour displays on /PLOPTS,FILE,on ! filename in legend /PLOPTS,LOGO,off ! no graphical logo in legend /PLOPTS,WINS,on ! Controls whether graphics windows automatically stretch or shrink to adjust to screen size as the legend column is turned off or on /PLOPTS,WP,0 ! Working plane (defaults to OFF) off /PLOPTS,DATE,2 ! Controls the display of the date and time in your legend, 2 (default) - both the date and time are shown /psymb,esys,1 ! Elemtentkoordinatensystem anzeigen /VSCALE,1,1,0 ! Vektoren, z. B. triad, skalieren /TRIAD,ORIG ! Shows the global XYZ coordinate triad on displays, ORIG — Display triad at global origin /pbc,forc,,2 ! plot force value next to symbol /REPLOT /units,si ! SI-Einheiten (1N/mm^2=1e6N/m^2) *abbr,Hertz_Linienkontakt,/input,Hertz_Linienkontakt.inp !Abkürzung anlegen /title, Hertz_Linienkontakt !-Variablen------------- Zylinderradius=0.1 ! [m] Zylinderlaenge=1 ! [m] Plattendicke=0.01 ! [m] Plattenbreite=0.25*Zylinderradius Innenradius_rel=0.95 ! relativer Innenradius [% Zylinderradius] Reibwert=0 El_Z_Axis_Zyl=Zylinderradius/3 El_Z_Axis_Pla=Zylinderradius/3 ElNormZyl=Zylinderradius/10 ElFeinZyl=Zylinderradius/25 ElKontaktZyl=Zylinderradius/75 ElNormPla=Plattendicke/5 ElFeinPla=Zylinderradius/20 esize_verhaeltnis_Zyl_IR=5 esize_verhaeltnis_Zyl_AR=20 esize_verhaeltnis_Pla=30 E1=2.1E11 ! [N/m^2] NUE1=0.3 ! [-] E2=2.1E11 ! [N/m^2] NUE2=0.3 ! [-] E3=2.1E11 ! [N/m^2] NUE3=0.3 ! [-] Interferenz=0 ! => !KEYOPT,3,5,1, CLOSE GAP on CONTA17x !Interferenz=Zylinderradius/1e6 ! => flächenhafte Berührung, Initial size of gap. A zero (or positive) value assumes an initially open gap. ! A negative value defines an interference condition Innenradius=Zylinderradius*Innenradius_rel S=Innenradius_rel Anfangsverschiebung=Zylinderradius/1e4 ! => Kontakt sicherstellen SelectTolerance=0.000000000001 PI=2*ASIN(1) K1=0.9*Innenradius*sin(PI/4) !-Preprocesseing-------- /prep7 ET,1,SOLID45 ! 3-D SOLID ELEMENTS, auch: SOLID185 ET,2,TARGE170 ! 3-D TARGET ELEMENTS ET,3,CONTA174 ! 2-D CONTACT ELEMENTS !ET,3,CONTA175 ! 2-D CONTACT ELEMENTS, 2-D/3-D Node-to-Surface Contact ET,4,SOLID92 ! 3-D SOLID ELEMENTS ET,5,SOLID95 ! 3-D SOLID ELEMENTS *if,Interferenz,eq,0,then KEYOPT,3,5,1 ! CLOSE GAP *else *endif MP,EX,1,E1 ! E-Modul [N/m^2], CYLINDER PROPERTIES MP,NUXY,1,NUE1 ! Querkontraktionszahl MP,MU,1,Reibwert ! Reibwert für Material Nr.1 MP,EX,2,E2 ! E-Modul [N/m^2], BLOCK PROPERTIES MP,NUXY,2,NUE2 ! Querkontraktionszahl MP,MU,2,Reibwert ! Reibwert für Material Nr.2 MP,EX,3,E2 ! E-Modul [N/m^2], CYLINDER PROPERTIES MP,NUXY,3,NUE2 ! Querkontraktionszahl MP,MU,3,Reibwert ! Reibwert für Material Nr.3 /status,global ! hier erkennt ANSYS noch 3 Materialien !-Geometrie------------- ! Die Rollenaussenprofile werden in verschiedenen x-z-Ebenen (im Abstand "Ebenenabstand=Zylinderradius+radiusZ2-Interferenz") gezeichnet ! Symmetrie zur y-z-Ebene => nur Profil für halbe Zylinderhöhe (x>0) zeichnen k,1,Innenradius,Zylinderradius-Interferenz,0 ! Zylinder oben k,2,Innenradius,Zylinderradius-Interferenz,Zylinderlaenge k,3,Zylinderradius,Zylinderradius-Interferenz,Zylinderlaenge k,4,Zylinderradius,Zylinderradius-Interferenz,0 k,5,0,Zylinderradius-Interferenz,Zylinderlaenge k,6,0,Zylinderradius-Interferenz,0 a,1,2,3,4 ! area 1 a,1,2,5,6 ! area 2 vrot,1,,,,,,5,6,180,2 vrot,2,,,,,,5,6,180,2 block,0,K1,Zylinderradius-Interferenz-(K1),Zylinderradius,0,Zylinderlaenge ! V5 block,-K1,0,Zylinderradius-Interferenz-(K1),Zylinderradius,0,Zylinderlaenge ! V6 vsbv,3,5 vsbv,4,6 allsel,all,all block,0,K1,Zylinderradius-Interferenz-(K1),Zylinderradius,0,Zylinderlaenge ! V block,-K1,0,Zylinderradius-Interferenz-(K1),Zylinderradius,0,Zylinderlaenge ! V allsel,all,all vsel,s,volu,,all,,1 ! alle Zylindervolumina mergen nummrg,all numcmp,all block,-Plattenbreite,0,-Plattendicke,0,0,Zylinderlaenge ! Platte ! V7 block,0,Plattenbreite,-Plattendicke,0,0,Zylinderlaenge ! Platte ! V8 seltol,SelectTolerance ! Selektiertoleranz bei Selektion nach loc => If VMIN = VMAX, Toler = SelectTolerance x VMIN allsel,all,all vsel,s,volu,,7,8,1,1 ! alle Plattenvolumina mergen nummrg,all seltol, allsel,all,all numcmp,all allsel,all,all /status,global ! hier erkennt ANSYS nur noch 1 Material !-Vernetzung------------ allsel,all,all lplot seltol,SelectTolerance ! Selektiertoleranz bei Selektion nach loc => If VMIN = VMAX, Toler = SelectTolerance x VMIN lsel,s,loc,z,zylinderlaenge/2,,1 lsel,r,loc,y,-Plattendicke,0 *GET,L_min,line,,NUM,MIN ! GET LOWEST LINE NUMBER lsel,u,line,,L_min lesize,all,El_Z_Axis_Pla allsel,all,all seltol, local,100,1,0,Zylinderradius-Interferenz,0 ! Zylinderkoordinatensystem in Zylinderachse seltol,SelectTolerance ! Selektiertoleranz bei Selektion nach loc => If VMIN = VMAX, Toler = SelectTolerance x VMIN lsel,s,loc,z,zylinderlaenge/2,,1 lsel,r,loc,x,0,Zylinderradius *GET,L_max,line,,NUM,MAX ! GET HIGHEST LINE NUMBER lsel,u,line,,L_max lesize,all,El_Z_Axis_Zyl allsel,all,all seltol, ! Zylinder: allsel,all,all lsel,s,loc,z,0,,1 lsel,r,loc,x,zylinderradius lsel,r,loc,y,180,270 lesize,all,ElFeinZyl,,,esize_verhaeltnis_Zyl_AR allsel,all,all allsel,all,all lsel,s,loc,z,0,,1 lsel,r,loc,x,Innenradius lsel,r,loc,y,180,270 lesize,all,ElFeinZyl*S,,,esize_verhaeltnis_Zyl_IR allsel,all,all allsel,all,all lsel,s,loc,z,0,,1 lsel,r,loc,x,zylinderradius lsel,r,loc,y,270,360 lesize,all,ElFeinZyl,,,1/esize_verhaeltnis_Zyl_AR allsel,all,all allsel,all,all lsel,s,loc,z,0,,1 lsel,r,loc,x,Innenradius lsel,r,loc,y,270,360 lesize,all,ElFeinZyl*S,,,1/esize_verhaeltnis_Zyl_IR allsel,all,all allsel,all,all lsel,s,loc,z,Zylinderlaenge,,1 lsel,r,loc,x,zylinderradius lsel,r,loc,y,180,270 lesize,all,ElFeinZyl,,,esize_verhaeltnis_Zyl_AR allsel,all,all allsel,all,all lsel,s,loc,z,Zylinderlaenge,,1 lsel,r,loc,x,Innenradius lsel,r,loc,y,180,270 lesize,all,ElFeinZyl*S,,,esize_verhaeltnis_Zyl_IR allsel,all,all allsel,all,all lsel,s,loc,z,Zylinderlaenge,,1 lsel,r,loc,x,zylinderradius lsel,r,loc,y,270,360 lesize,all,ElFeinZyl,,,1/esize_verhaeltnis_Zyl_AR allsel,all,all allsel,all,all lsel,s,loc,z,Zylinderlaenge,,1 lsel,r,loc,x,Innenradius lsel,r,loc,y,270,360 lesize,all,ElFeinZyl*S,,,1/esize_verhaeltnis_Zyl_IR allsel,all,all csys,0 ! Platte: allsel,all,all lsel,s,loc,z,0,,1 lsel,r,loc,x,0,Plattenbreite lsel,r,loc,y,0 lesize,all,ElFeinPla,,,1/esize_verhaeltnis_Pla allsel,all,all allsel,all,all lsel,s,loc,z,0,,1 lsel,r,loc,x,0,Plattenbreite lsel,r,loc,y,-Plattendicke lesize,all,ElFeinPla,,,esize_verhaeltnis_Pla allsel,all,all allsel,all,all lsel,s,loc,z,0,,1 lsel,r,loc,x,-Plattenbreite,0 lsel,r,loc,y,0 lesize,all,ElFeinPla,,,esize_verhaeltnis_Pla allsel,all,all allsel,all,all lsel,s,loc,z,0,,1 lsel,r,loc,x,-Plattenbreite,0 lsel,r,loc,y,-Plattendicke lesize,all,ElFeinPla,,,1/esize_verhaeltnis_Pla allsel,all,all allsel,all,all lsel,s,loc,z,zylinderlaenge,,1 lsel,r,loc,x,0,Plattenbreite lsel,r,loc,y,0 lesize,all,ElFeinPla,,,1/esize_verhaeltnis_Pla allsel,all,all allsel,all,all lsel,s,loc,z,zylinderlaenge,,1 lsel,r,loc,x,0,Plattenbreite lsel,r,loc,y,-Plattendicke lesize,all,ElFeinPla,,,esize_verhaeltnis_Pla allsel,all,all allsel,all,all lsel,s,loc,z,zylinderlaenge,,1 lsel,r,loc,x,-Plattenbreite,0 lsel,r,loc,y,0 lesize,all,ElFeinPla,,,esize_verhaeltnis_Pla allsel,all,all allsel,all,all lsel,s,loc,z,zylinderlaenge,,1 lsel,r,loc,x,-Plattenbreite,0 lsel,r,loc,y,-Plattendicke lesize,all,ElFeinPla,,,1/esize_verhaeltnis_Pla allsel,all,all /status,global allsel,all,all esize,ElKontaktZyl !type,1 type,5 mshape,0,3d ! 0 = brick meshing mshkey,1 ! 1 = mapped meshing mat,1 vmesh,1 ! Mesh Aussenring vmesh,2 ! Mesh Aussenring allsel,all,all esize,ElNormZyl type,1 mshape,0,3d ! 0 = brick meshing mshkey,1 ! 1 = mapped meshing mat,1 vmesh,4 vmesh,5 esize,ElFeinZyl type,4 mshape,1,3d ! 0 = brick meshing mshkey,0 ! 1 = mapped meshing mat,1 vmesh,3 vmesh,6 allsel,all,all esize,ElNormPla type,5 mshape,0,3d ! 0 = brick meshing mshkey,1 ! 1 = mapped meshing mat,2 vmesh,7 vmesh,8 /view,1,1,1,1 ! /view,wn,xv,yv,zv -> Defines the viewing direction in window nr. for the display, view to origin triade /angle,1 ! Isometrische Ansicht in Window Nr. (1) /replot,fast allsel,all,all !-Kontakte-------------- allsel,all,all /view,1,1,0,1 ! /view,wn,xv,yv,zv -> Defines the viewing direction in window nr. for the display, view to origin triade /angle,1 /replot,fast eplot asel,s,area,,5,10,5 NSLA,s,1 ! SELECT CONTACT NODES ON CYLINDER nsel,r,loc,x,-Plattenbreite,Plattenbreite,2*Plattenbreite REAL,1 ! share real constant set 1 TYPE,3 MAT,3 ESURF ! GENERATE 3-D CONTA174 ELEMENTS allsel,all,all asel,s,area,,33,39,6 NSLA,,1 ! SELECT TARGET NODES ON BLOCK nsel,r,loc,z,0,Zylinderlaenge REAL,1 ! share real constant set 1 TYPE,2 MAT,3 ESURF ! GENERATE 3-D TARGE170 ELEMENTS allsel,all,all CNCHECK,all ! check contact status !-Randbedingungen------- allsel,all,all NSEL,S,LOC,Y,-Plattendicke ! Platte unten festsetzen D,ALL,UY allsel,all,all NSEL,S,node,,all D,ALL,UX ! Hertz: nur Verschiebung in Kraftwirkungsrichtung möglich D,ALL,UZ allsel,all,all /status,all !-Lasten---------------- NSEL,S,LOC,Y,Zylinderradius ! SELECT TOP EDGE OF MODEL CP,1,UY,ALL ! COUPLE NODES ON TOP EDGE *GET,NC,NODE,,NUM,MIN ! GET LOWEST NODE NUMBER (MASTER) NSEL,ALL FINISH SAVE,MODEL3D /SOLU D,NC,UY,-Anfangsverschiebung ! APPLY SMALL DISPLACEMENT TO ENGAGE CONTACT nsub,2,10,1 ncnv,0 neqit,100 nlgeom,on SOLVE ! SOLVE FIRST LOAD STEP DDELE,NC,UY ! DELETE IMPOSED DISPLACEMENT F,NC,FY,-1e7 ! APPLY LOAD ON MODEL (= -1e4 N/mm) nsub,2,10,1 ncnv,0 neqit,100 nlgeom,on SOLVE ! SOLVE SECOND LOAD STEP FINISH