! -------------------------------------------------------------------------------------------- ! Ansys-Eingabedatei zur Untersuchung der Lastverteilung und der Verformung der Laufbahn eines ! Außenrings in einem Gehäuse (dreidimensionales Modell) ! ! Ausgabe: - Laufbahnverformung des Außenrings im belasteten Zustand ! - Lastverteilung im Wälzlager ! -------------------------------------------------------------------------------------------- ! -------------------------------------------------------------------------------------------- !22.01.2012 V. Aul ... ! xx.xx.20xx X. Xxx ... ! -------------------------------------------------------------------------------------------- ! -------------------------------------------------------------------------------------------- ! Mit Notepad++ öffnen ! Ansys v11.0 ! -------------------------------------------------------------------------------------------- ! -------------------------------------------------------------------------------------------- ! Einheiten: ! Länge Zeit Masse Kraft ! m s kg N ! -------------------------------------------------------------------------------------------- FINISH /CLEAR,START /FILNAME,Lastverteilung_Welle_starr_Gehaeuse_starr /TITLE,Lastverteilung_Welle_starr_Gehaeuse_starr (3D) ! -------------------------------------------------------------------------------------------- ! Variablen ! -------------------------------------------------------------------------------------------- ! Wälzlagergeometrie ! Wälzlagergeometrie D_AR_Mantel = 0.200 ![m] Außenradius Außenring "D" vgl, Lagerkatalog FAG-INA D_AR_1 = 0.1767 ![m] Innenradius Außenring "D1" vgl. Lagerkatalog FAG-INA D_IR_IFlaeche = 0.130 ![m] Innenradius(Wellendurchmesser) Innenring "d" vgl. Lagerkatalog FAG-INA D_IR_1 = 0.1547 ![m] AUßenradius Innenring "d1" vgl. Lagerkatalog FAG-INA B = 0.022 ![m} Breite der Lagerringe "B" vgl. Lagerkatalog FAG-INA WK_anz = 18 ! Anzahl der Wälzkörper Wk_d = 0.016 ![m] Durchmesser der Wälzkörper Wk_ausr = 1 ! 0 = Scheitelstellung , 1= Mittelstellung Gh_di = 0.194 ![m] Anschlussmaß "Da" vgl. Lagerkatalog FAG-INA D_AR_1 !Reibungbeiwerte mu_WK_IR=0.2 !Innenring-Wälzkörper mu_WK_AR=0.2 !Außenring-Wälzkörper p_Netz = 0.005 !Vernetzungsparameter (je kleiner, desto feiner vernetzt); ~ D/40 !Kontaktparameter k_4 = 0 !0 = standard k_9 = 0 !0 = standard k_5 = 1 !1 = standard k_6 = 1 !1 = standard ! Verschiebung [m] Innenring bzgl. Gehäuse bzw. Belastung delta_y = -0.000015 ! [m] radial delta_x = -0.000000 ! [m] radial delta_z = -0.000015 ! [m] axial belastung_y = -10000 ! [N] radial belastung_x = -00000 ! [N] radial belastung_z = 00000 ! [N] axial Lkey = 1 ! Belastungsschalter : 1 = Kraftsteuerung ; 2 = Verschiebungssteuerung CSkey = 0 ! Schalter: Ausgabe der Knotenpositionen in Zylinderkoordinaten = 1 , kartesische Koordinaten = 0 ! -------------------------------------------------------------------------------------------- !Zwischengrößen erzeugen ! --------------------------------------------------------------------------------------------- R_wk = D_IR_1/2+(D_AR_1-D_IR_1)/4 !Durchmesser zwischen Wälzkörpermittelpunkten /PREP7 ! -------------------------------------------------------------------------------------------- ! Elemente und Materialeigenschaften ! --------------------------------------------------------------------------------------------- !Außenring ET,1,SOLID185 !3D 8-Knoten lineares Volumenelement KEYOPT,1,2,3 !enhanced strain formulation !Innenring ET,2,SOLID185 KEYOPT,1,2,3 !Wälzkörper ET,3,SOLID185 KEYOPT,1,2,3 !WK-IR-KONTAKT !Kontaktelemete TARGE170 und CONTA173 für Kontakt Wälzkörper-Innenring *DO,j,1,WK_anz ET,600+j,TARGE170 ET,500+j,CONTA173 KEYOPT,500+j,6,k_6 !Aktualisierung der Kontaktsteifigkeit bei jeder Iteration KEYOPT,500+j,4,k_4 !Stelle der Kontakterkennung sind die Integrationspunkte KEYOPT,500+j,9,k_9 !Überlappungen und Offsets werden berücksichtigt und nicht korrigiert KEYOPT,500+j,10,2 !Das Programm darf die Kontaktsteifigkeit automatisch aktualisieren KEYOPT,500+j,5,k_5 !Vernetzungsbedingte Lücken werden automatisch geschlossen MP,MU,500+j,mu_WK_IR !Reibungsbeiwerte R,500+j *ENDDO !WK-IR-KONTAKT !Kontaktelemete TARGE170 und CONTA173 für Kontakt Wälzkörper-AUßenring *DO,j,1,WK_anz ET,800+j,TARGE170 ET,700+j,CONTA173 KEYOPT,700+j,6,k_6 !Aktualisierung der Kontaktsteifigkeit bei jeder Iteration KEYOPT,700+j,4,k_4 !Stelle der Kontakterkennung sind die Integrationspunkte KEYOPT,700+j,9,k_9 !Überlappungen und Offsets werden berücksichtigt und nicht korrigiert KEYOPT,700+j,10,2 !Das Programm darf die Kontaktsteifigkeit automatisch aktualisieren KEYOPT,700+j,5,k_5 !Vernetzungsbedingte Lücken werden automatisch geschlossen MP,MU,700+j,mu_WK_IR !Reibungsbeiwerte R,700+j *ENDDO ET,10,MASS21 !Elementdefinition für Master-Node; nötig für Modellierung der starren Welle R,10,0,0,0,0,0,0,0 ! -------------------------------------------------------------------------------------------- ! Materialeigenschaften ! --------------------------------------------------------------------------------------------- !Außenring MP,EX,1,2.08e11 ![N/m^2] E-Modul MP,NUXY,1,0.3 ! Querkontraktionszahl MP,DENS,1,7.85e3 ![kg/m^3] Dichte !Innenring MP,EX,2,2.08e11 MP,NUXY,2,0.3 MP,DENS,2,7.85e3 !Wälzkörper MP,EX,3,2.08e11 MP,NUXY,3,0.3 MP,DENS,3,7.85e3 ! -------------------------------------------------------------------------------------------- ! Geometrie 1 - Innenring+Aussenring ! -------------------------------------------------------------------------------------------- !Lagerquerschnitt modellieren WPLANE,0,0,0,0,0,0,1,0,1,0 !Arbeitsebene drehen RECTNG,-B/2,B/2,D_IR_IFlaeche/2,D_AR_Mantel/2 !Flächen erzeugen, geeignet subtrahieren und teilen CYL4,0,R_wk,WK_d/2 ASBA,1,2 RECTNG,-B/2,B/2,D_IR_1/2,D_AR_1/2 ASBA,3,1 !Flächen teilen um mapped meshing zu ermöglichen KSEL,S,,,16 !zusätzliche Keypoints erzeugen KSEL,A,,,15 KGEN,2,ALL,,,,(D_AR_Mantel-D_AR_1)/2 ALLS KSEL,S,,,6 KGEN,2,ALL,,,,(D_AR_Mantel/2-(R_wk+Wk_d/2)) ALLS KSEL,S,,,13 KSEL,A,,,14 KGEN,2,ALL,,,,-(D_IR_1-D_IR_IFlaeche)/2 ALLS KSEL,S,,,8 KGEN,2,ALL,,,,-((R_wk-Wk_d/2)-D_IR_IFlaeche/2) ALLS !Linien erzeugen L,16,7 !2 L,6,17 !4 L,15,5 !5 L,13,18 L,8,20 L,14,19 !Flächen teilen ASBL,4,2 ASBL,3,4 ASBL,5,5 !Flächen verkleben AGLUE,1,4,3,6 !Flächen teilen ASBL,2,6 ASBL,7,7 ASBL,8,8 !Flächen verkleben AGLUE,5,2,7,9 !Rotationsachse zum austragen definieren K,100,0,0,0 K,101,0,0,1 !Volumen erzeugen , Rotieren CSYS,1 VROTAT,1,4,3,6 , , ,100,101 ,360, , VROTAT,5,2,7,9 , , ,100,101 ,360, , CSYS,0 !Komponenten bilden vSEL,S,,,1,16,1 CM,AR_V,VOLU vSEL,S,,,17,32,1 CM,IR_V,VOLU ! -------------------------------------------------------------------------------------------- ! Geometrie 2 - Wälzkörper ! -------------------------------------------------------------------------------------------- !ersten Wälzkörper im Ursprung erzeugen SPH4,0,0,Wk_d/2 !Volumen durch Arbeitsebenen teilen um mapped meshing zu ermöglichen (Vier Flächen pro Volumen) WPCSYS,-1,0 VSBW,33 WPLANE,0,0,0,0,0,0,1,0,1,0 VSBW,34 VSBW,35 WPLANE,0,0,0,0,0,0,1,1,0,0 VSEL,S,,,33,40,1 VSBW, ALL WPCSYS,-1,0 !Wälzkörper an Ausgangsposition verschieben VSEL,S,,,33,48,1 !Komponente erzeugen CM,WK_1_V,VOLU VGEN,2,WK_1_V,,,,R_wk,,,,1 !Bei Mittelstellung ersten WK weiter drehen *IF,Wk_ausr,EQ,1,THEN CSYS,1 VGEN,2,WK_1_V,,,,360/WK_anz/2,,,,1 CSYS,0 *ENDIF !restliche WK erzeugen CSYS,1 VGEN,WK_anz,WK_1_V,,,,360/WK_anz,,,, CSYS,0 !Komponente bilden VSEL,S,,,33,1000,1 CM,WK_ges,VOLU ! -------------------------------------------------------------------------------------------- ! Vernetzen - Welle, Innenring, Wälzkörper, Außenring ! -------------------------------------------------------------------------------------------- !Aussenring vernetzen ALLS TYPE,1 MAT,1 ESIZE,p_Netz MSHKEY,1 VMESH,AR_V !Innenring vernetzen ALLS TYPE,2 MAT,2 ESIZE,p_Netz MSHKEY,1 VMESH,IR_V !Wälzkörper vernetzen ALLS TYPE,3 MAT,3 ESIZE,p_Netz MSHKEY,1 VMESH,WK_ges ! -------------------------------------------------------------------------------------------- ! Definition Kontakt: Wälzkörper - Aussenring / Wälzkörper -Innenring ! -------------------------------------------------------------------------------------------- !Komponenten vorbereiten !Aussenring,Laufbahnflächen ASEL,S,,,11 ASEL,A,,,14 ASEL,A,,,18 ASEL,A,,,22 ASEL,A,,,39 ASEL,A,,,35 ASEL,A,,,31 ASEL,A,,,28 ASEL,A,,,45 ASEL,A,,,48 ASEL,A,,,52 ASEL,A,,,56 ASEL,A,,,62 ASEL,A,,,64 ASEL,A,,,67 ASEL,A,,,70 NSLA,S,1 CM,ARWK_cont_surf,NODE !Innenring,Laufbahnflächen ALLS ASEL,S,,,92 ASEL,A,,,97 ASEL,A,,,101 ASEL,A,,,105 ASEL,A,,,75 ASEL,A,,,80 ASEL,A,,,84 ASEL,A,,,88 ASEL,A,,,109 ASEL,A,,,114 ASEL,A,,,118 ASEL,A,,,122 ASEL,A,,,126 ASEL,A,,,130 ASEL,A,,,133 ASEL,A,,,136 NSLA,S,1 CM,IRWK_cont_surf,NODE !Wälzkörper ! Anpassen der Wälzkörperanzahl weil nur die Hälfte des Lagers modelliert wird (Symmetrie) !*IF,Wk_ausr,EQ,0,THEN !Schleife: Für jeden Wälzkörper werden die Oberflächenknoten zu Komponenten zusammengefasst *DO,j,1,WK_anz CSYS,1 !lokale Koordinatensysteme definieren um Wälzkörper zu beschreiben !Wahl je nach Wälzkörperausrichtung *IF,Wk_ausr,EQ,0,THEN CLOCAL,10+j,2,R_wk,(j-1)*360/WK_anz+90,0 *ELSEIF,WK_ausr,EQ,1 CLOCAL,10+j,2,R_wk,(j-0.5)*360/WK_anz+90,0 *ENDIF !Koordinatensystem definieren CSYS,10+j !Knoten auf der Oberfläche der Wälzkörper auswählen VSEL,S,LOC,X,,Wk_d/2+1 VSEL,U,,,IR_V VSEL,U,,,AR_V NSLV,S,1 NSEL,R,EXT NSEL,R,LOC,X,Wk_d/2 !Komponenten der Knoten auf WK-Oberfläche bilden CM,WK_%(j)%_targ_surf,NODE *ENDDO ! -------------------------------------------------------------------------------------------- ! Definition Kontakt: Wälzkörper - Laufbahndurchmesser Innenring ! -------------------------------------------------------------------------------------------- !Schleife über die Wälzkörper zum Aufbringen der Kontaktelemente *DO,i,1,WK_anz !Konstantenset 500+i wählen REAL,500+i ALLS !TARGE169 auf Innenring aufbringen CMSEL,S,IRWK_cont_surf,NODE TYPE,600+i ESURF,ALL !CONTA173 auf Wälzkörper aufbringen CMSEL,S,WK_%(i)%_targ_surf,NODE TYPE,500+i MAT,500+i ESURF,ALL *ENDDO ! -------------------------------------------------------------------------------------------- ! Definition Kontakt: Wälzkörper - Laufbahndurchmesser Außenring ! -------------------------------------------------------------------------------------------- !Schleife über die Wälzkörper zum Aufbringen der Kontaktelemente *DO,i,1,WK_anz !Konstantenset 700+i wählen REAL,700+i ALLS !TARGE170 auf Außenring aufbringen CMSEL,S,ARWK_cont_surf,NODE TYPE,800+i ESURF,ALL !CONTA173 auf Wälzkörper aufbringen CMSEL,S,WK_%(i)%_targ_surf,NODE TYPE,700+i MAT,700+i ESURF,ALL *ENDDO ! -------------------------------------------------------------------------------------------- ! Erzeugen der constraint equations ("Spinnennetz") am Innenring -> starre Welle ! -------------------------------------------------------------------------------------------- CSYS,0 ALLS !maximale Knotennummer ermitteln *GET,nd_max,NODE,,NUM,MAX !Masternode erzeugen + Element aufbringen, im Koordinatenursprung mst_nd=nd_max+1 N,mst_nd,0,0,0 TYPE,10 REAL,10 E,mst_nd ! Knoten auf Innenfläche des Innenrings auswählen ALLS ASEL,S,,,73 ASEL,A,,,78 ASEL,A,,,82 ASEL,A,,,86 ASEL,A,,,90 ASEL,A,,,95 ASEL,A,,,99 ASEL,A,,,103 ASEL,A,,,107 ASEL,A,,,112 ASEL,A,,,116 ASEL,A,,,120 ASEL,A,,,124 ASEL,A,,,128 ASEL,A,,,131 ASEL,A,,,134 !Komponente bilden Fläche CM,IR_A,AREA CMSEL,S,IR_A,AREA !Knoten auswählen NSLA,S,1 CSYS,0 NSEL,A,,,mst_nd !Kopplung aller Knoten der Innenfläche des Innenrings mit dem Masternode bzgl. der Freiheitsgrade UX,UY,UZ --> starre Welle CERIG,mst_nd,ALL,UXYZ ! -------------------------------------------------------------------------------------------- ! Koordinatensystem an Mittenknoten der WK drehen um Freiheitsgradkopplung in Umfangrichtung !zu definieren zu können (Käfig) ! -------------------------------------------------------------------------------------------- !Mittenknoten der WK auswählen !ersten Knoten auswählen ALLS CSYS,11 NSEL,S,LOC,X,,0 !Schleife über restliche Knoten *DO,i,1,WK_anz-1 CSYS,11+i NSEL,A,LOC,X,,0 *ENDDO !Komponente bilden (Mittenknoten der WK) CM,WK_MID_NODES,NODE CMSEL,S,WK_MID_NODES,NODE !Knotenkoordinatensystem nach dem zylindrischen globalen Koordinatensystem ausrichten, Y-Richtung = Umfangsrichtung CSYS,1 NROTAT,ALL CSYS,0 ! -------------------------------------------------------------------------------------------- ! Randbedingungen ! -------------------------------------------------------------------------------------------- FINISH /SOLU !statische Berechnung ANTYPE,STATIC !Knoten am Außenring des Gehäuses festhalten, starres Gehäuse, Einbaumaß Da(GH_di) berücksichtigt !Knoten an Mantelfläche des AR auswählen CSYS,1 ALLS NSEL,S,LOC,X,D_AR_Mantel/2 !Komponente bilden CM,AR_RB_1,NODE CSYS,0 !Knoten an Stirnseiten des AR auswählen (definiert durch GH Geometrie (GH_di bzw. Da)) ALLS NSEL,S,LOC,Z,B/2 NSEL,A,LOC,Z,-B/2 CSYS,1 NSEL,R,LOC,X,GH_di/2,D_AR_Mantel/2 CSYS,0 !Komponente bilden CM,AR_RB_2,NODE !Beide Komponenten (Alle relevante Knoten) auswählen und festhalten CMSEL,S,AR_RB_1,NODE CMSEL,A,AR_RB_2,NODE D,ALL,UX,0 D,ALL,UY,0 D,ALL,UZ,0 !Wälzkörper auf Abstand halten (Käfig), Kopplung der Verschiebung in Umfangsrichtung CMSEL,S,WK_MID_NODES,NODE CP,666,UY,WK_MID_NODES !Kräfte Verschiebungen der Welle definieren ALLS !Kraftgesteuert *IF,Lkey ,EQ,1,THEN F,mst_nd,FY,belastung_y F,mst_nd,FX,belastung_x F,mst_nd,FZ,belastung_z D,mst_nd,UX,0 D,mst_nd,UZ,0 !Verschiebungsgesteuert *ELSEIF,Lkey,EQ,2 D,mst_nd,UY,delta_y D,mst_nd,UX,delta_x D,mst_nd,UZ,delta_z *ENDIF ! Anzahl der Elemente im gesamten Modell *GET,no_elem,ELEM,,COUNT ! -------------------------------------------------------------------------------------------- !Solvereinstellungen ! -------------------------------------------------------------------------------------------- ! keine Berücksichtigung großer Verformungen NLGEOM,OFF !Automatische Step Size Control aktiviert AUTOTS,ON !erster Zeitschritt bei T= 0.01 von 1 NSUB,100 !Wähle Sparse direkten Solver EQSLV,PCG,1.0E-8,2.9 ! Activates a line search to be used with Newton-Raphson, automatische Wahl LNSRCH,AUTO *GET,CPU_strt,ACTIVE,,TIME,CPU ! CPU-Zeit vor Simulation [h] *GET,wll_strt,ACTIVE,,TIME,WALL !Lösung starten SOLVE *GET,CPU_end,ACTIVE,,TIME,CPU ! CPU-Zeit nach Simulation [h] *GET,wll_end,ACTIVE,,TIME,WALL CPU_tme=CPU_end-CPU_strt ! CPU-Zeit für Simulation [s] wll_tme=(wll_end-wll_strt)*3600 ! Tatsächliche Zeit in [s] ! -------------------------------------------------------------------------------------------- ! Postprocessing ! -------------------------------------------------------------------------------------------- /POST1 ! ganzes Modell darstellen /EXPAND,2,RECT,HALF,0.00001 /REPLOT ! -------------------------------------------------------------------------------------------- ! Auslesen der ! - Knotennummern, ! - Anfangspositionen, ! - Verschiebungen, ! - Kontaktspannungen und ! - Reaktionskräfte ! an Laufbahn des Außenrings ! -------------------------------------------------------------------------------------------- !Bestimmen der Anzahl der Knoten auf Aussenringlaufbahn !Relevante Flächen auswählen ALLS ASEL,S,,,14 ASEL,A,,,18 ASEL,A,,,35 ASEL,A,,,31 ASEL,A,,,64 ASEL,A,,,67 ASEL,A,,,48 ASEL,A,,,52 !Knoten auswählen NSLA,S,1 !Komponente bilden CM,ARWK_cont_surf1,NODE !Knotenanzahl ermitteln *GET,AR_lfb_no_nd,NODE,,COUNT ALLS !Bestimmen der Anzahl der Knoten für phi = konst in Zylinderkoordinaten auf Aussenringlaufbahn CSYS,1 CMSEL,S,ARWK_cont_surf1,NODE NSEL,R,LOC,Y,180 CSYS,0 NSEL,R,LOC,Y,0 !Knotenanzahl ermitteln *GET,no_nd_phi_const,NODE,,COUNT CSYS,0 !Array zum spreichern der Ergebnisse *DIM,res,ARRAY,AR_lfb_no_nd,11 ! Ausgehend von der "12 Uhr Position" werden die Knotennummern gegen den Uhrzeigersinn ausgelesen ! und in res(j,1) (Array) gespeichert (1. Spalte) ! 2. Quadrant !Knoten auf Außenringlaufbahn auswählen 90-180° CSYS,1 CMSEL,S,ARWK_cont_surf1 NSEL,R,LOC,Y,90,180 CSYS,0 !Knotenanzahl ermitteln *GET,AR_nd_1,NODE,,COUNT !Schleife über alle Knoten *DO,j,1,AR_nd_1 !Knotenpostionen des nächsten Knoten ermitteln *GET,xx,NODE,,MXLOC,X *GET,yy,NODE,,MXLOC,Y *GET,zz,NODE,,MXLOC,Z !Knotennummer dieses Knoten ermitteln nd=NODE(xx,yy,zz) !Ablegen in dem Array res(j,1)=nd !Entfernen dieses Knoten aus dem Set NSEL,U,,,nd *ENDDO ! 3. Quadrant !Knoten auf Außenringlaufbahn auswählen 180-270° CSYS,1 CMSEL,S,ARWK_cont_surf1 NSEL,R,LOC,Y,180,270 !Schon ausgewählte Knoten auslassen NSEL,U,LOC,Y,180 CSYS,0 !Knotenanzahl ermitteln *GET,AR_nd_2,NODE,,COUNT !Schleife über alle Knoten *DO,j,1,AR_nd_2 !Knotenpostionen des nächsten Knoten ermitteln *GET,xx,NODE,,MNLOC,X *GET,yy,NODE,,MXLOC,Y *GET,zz,NODE,,MXLOC,Z !Knotennummer dieses Knoten ermitteln nd=NODE(xx,yy,zz) !Ablegen in dem Array res(j+AR_nd_1,1)=nd !Entfernen dieses Knoten aus dem Set NSEL,U,,,nd *ENDDO ! 4. Quadrant !Knoten auf Außenringlaufbahn auswählen 270-360° CSYS,1 CMSEL,S,ARWK_cont_surf1 NSEL,R,LOC,Y,270,360 !Schon ausgewählte Knoten auslassen NSEL,U,LOC,Y,270 CSYS,0 !Knotenanzahl ermitteln *GET,AR_nd_3,NODE,,COUNT !Schleife über alle Knoten *DO,j,1,AR_nd_3 !Knotenpostionen des nächsten Knoten ermitteln *GET,xx,NODE,,MXLOC,X *GET,yy,NODE,,MXLOC,Y *GET,zz,NODE,,MXLOC,Z !Knotennummer dieses Knoten ermitteln nd=NODE(xx,yy,zz) !Ablegen in dem Array res(j+AR_nd_1+AR_nd_2,1)=nd !Entfernen dieses Knoten aus dem Set NSEL,U,,,nd *ENDDO ! 1. Quadrant !Knoten auf Außenringlaufbahn auswählen 0-90° CSYS,1 CMSEL,S,ARWK_cont_surf1 !Schon ausgewählte Knoten auslassen NSEL,R,LOC,Y,0,90 NSEL,U,LOC,Y,90 NSEL,U,LOC,Y,0 CSYS,0 !Knotenanzahl ermitteln *GET,AR_nd_4,NODE,,COUNT !Schleife über alle Knoten *DO,j,1,AR_nd_4 !Knotenpostionen des nächsten Knoten ermitteln *GET,xx,NODE,,MXLOC,X *GET,yy,NODE,,MXLOC,Y *GET,zz,NODE,,MXLOC,Z !Knotennummer dieses Knoten ermitteln nd=NODE(xx,yy,zz) !Ablegen in dem Array res(j+AR_nd_1+AR_nd_2+AR_nd_3,1)=nd !Entfernen dieses Knoten aus dem Set NSEL,U,,,nd *ENDDO ! -------------------------------------------------------------------------------------------- ! Speichern der Ergebnisse ! -------------------------------------------------------------------------------------------- ALLS !Reaktionskräfte am Lager ermitteln !Knoten auf Laufbahn des Aussenrings auswählen CMSEL,S,ARWK_cont_surf,NODE !Kräfte an diesen Knoten aufsummieren FSUM *GET,FR_X,FSUM,,ITEM,FX *GET,FR_Y,FSUM,,ITEM,FY *GET,FR_Z,FSUM,,ITEM,FZ ! Verschiebungen, Spannungen, Kräfte, ... !Schleife über alle Knoten *DO,j,1,AR_lfb_no_nd ALLS !Ausgabe der Koordinaten in Zylinder- oder in kartesischen Koordinaten CSYS,0 *IF,CSkey,EQ,1,THEN CSYS,1 *ENDIF !Knotenkoordinaten in Array ablegen *GET,res(j,2),NODE,res(j,1),LOC,X *GET,res(j,3),NODE,res(j,1),LOC,Y *GET,res(j,4),NODE,res(j,1),LOC,Z !Knotenverschiebungen in Array ablegen CSYS,0 *GET,res(j,5),NODE,res(j,1),U,X *GET,res(j,6),NODE,res(j,1),U,Y *GET,res(j,7),NODE,res(j,1),U,Z ! *GET,res(j,8),NODE,res(j,1),CONT,PRES !Knoten auswählen und Knotenkräfte für diesen Knoten ermitteln NSEL,S,,,res(j,1) FSUM,,BOTH !Knotenkräfte in Array ablegen *GET,res(j,9),FSUM,0,ITEM,FX *GET,res(j,10),FSUM,0,ITEM,FY *GET,res(j,11),FSUM,0,ITEM,FZ *ENDDO ! -------------------------------------------------------------------------------------------- ! Ausgabedatei ! ----------------------------------------------------------------- *CREATE,MAKRO_lastvert ! Makro erzeugen *CFOPEN,'RES_Lastverteilung','dat' ! Datei "RES_Lastverteilung.dat" erstellen ! Überschrift *VWRITE ('*************************************************************************************') *VWRITE ('Ergebnisse aus Berechnungen am dreidimensionalen Wälzlagermodell:') *VWRITE (' - Reaktionskräfte am Gesamtlager') *VWRITE (' - Positionen und Verschiebungen der Knoten am Laufbahndurchmesser des Außenrings') *VWRITE (' - Kontaktspannungen und Kräfte an WK-AR-Kontakten') *VWRITE ('*************************************************************************************') *VWRITE (' ') *VWRITE (' ') ! Allgemeines *VWRITE ('Anzahl der Elemente im Gesamtmodell:') *VWRITE (' ') *VWRITE,'no_elem [-]' (T17,A) *VWRITE,no_elem (T17,F0.0) *VWRITE (' ') *VWRITE ('CPU-Zeit für Lösung [s]:') *VWRITE (' ') *VWRITE,'CPU_tme' (T17,A) *VWRITE,CPU_tme (T17,F0.4) *VWRITE (' ') *VWRITE ('Tatsächliche Zeit für Lösung [s]:') *VWRITE (' ') *VWRITE,'wll_tme' (T17,A) *VWRITE,wll_tme (T17,F0.4) *VWRITE (' ') ! Reaktionskräfte *VWRITE ('Reaktionskräfte durch Verschiebung des Innenrings:') *VWRITE (' ') *VWRITE,'Fz_axial [N]','Fy_rad [N]' (T17,A,T41,A,T65,A) *VWRITE,FR_Z,FR_Y (T17,F0.2,T41,F0.2) ! Anfangspositionen, Verschiebungen, Kontaktspannungen und Kontaktkräfte über den Umfang *VWRITE (' ') *VWRITE (' ') *VWRITE ('Anfangspositionen, Verschiebungen, Kontaktspannungen und Kontaktkräfte über den') *VWRITE ('Umfang. Start: 12 Uhr Position') *VWRITE (' ') *VWRITE,'Knoten','POS x','POS y','POS z','DISP x','DISP y','DISP z','p_max','Fx','Fy','Fz' (A,T17,A,T41,A,T65,A,T89,A,T113,A,T137,A,T161,A,T185,A,T209,A,T233,A) *VWRITE,' ','[m]','[m]','[m]','[m]','[m]','[m]','[MPa]','[N]','[N]','[N]' (A,T17,A,T41,A,T65,A,T89,A,T113,A,T137,A,T161,A,T185,A,T209,A,T233,A) *VWRITE,res(1,1),res(1,2),res(1,3),res(1,4),res(1,5),res(1,6),res(1,7),res(1,8),res(1,9),res(1,10),res(1,11) (F0.0,T17,E0.10,T41,E0.10,T65,E0.10,T89,E0.10,T113,E0.10,T137,E0.10,T161,E0.10,T185,E0.10,T209,E0.10,T233,E0.10) *VWRITE (' ')! *CFCLOS ! Datei schließen *END ! Makro beenden /INPUT,MAKRO_lastvert