# -*- coding: mbcs -*- # Do not delete the following import lines from abaqus import * from abaqusConstants import * import __main__ import math import numpy as np import os import section import regionToolset import displayGroupMdbToolset as dgm import part import material import assembly import step import interaction import load import mesh import optimization import job import sketch import visualization import xyPlot import displayGroupOdbToolset as dgo import connectorBehavior ### Benennung ### Modell_Name = ['WNV_B_FE_40_07_0125_12_08_11_75', 'WNV_B_FE_40_07_0125_06_08_11_75', 'WNV_B_FE_20_07_0125_12_08_11_75', 'WNV_B_FE_20_07_0125_06_08_11_75', 'WNV_B_FE_80_07_0125_06_08_11_75', 'WNV_B_FE_80_07_0125_16_08_11_75', 'WNV_B_FE_125_07_0125_12_08_11_75', 'WNV_B_FE_125_07_0125_25_08_11_75', 'WNV_B_FE_125_07_0125_40_08_11_75', 'WNV_Bi_FE_40_07_0125_12_08_11_75_H7n6', 'WNV_Bi_FE_40_07_0125_12_08_11_75_00', 'WNV_Bi_FE_40_07_0125_12_08_11_75_H7h6', 'WNV_Bi_FE_40_07_0125_12_08_11_75_H8e8', 'WNV_B_FE_40_07_0125_12_08_11_20', 'WNV_B_FE_40_07_0125_12_08_11_150', 'WNV_B_FE_40_00_0125_12_08_11_75', 'WNV_B_FE_40_03_0125_12_08_11_75', 'WNV_B_FE_40_085_0125_12_08_11_75', 'WNV_B_FE_40_07_0125_12_05_11_75', 'WNV_B_FE_40_07_0125_12_08_05_75', 'WNV_B_FE_40_07_0125_12_08_15_75', 'WNV_B_FE_40_07_005_12_08_11_75', 'WNV_B_FE_40_07_04_12_08_11_75'] #WNV-Bezeichnungen Wellen_Name = ['W_B_FE_40_07_0107_0125_12', 'W_B_FE_40_07_0107_0125_06', 'W_B_FE_20_07_0053_0125_12', 'W_B_FE_20_07_0053_0125_06', 'W_B_FE_80_07_0214_0125_06', 'W_B_FE_80_07_0214_0125_16', 'W_B_FE_125_07_0334_0125_12', 'W_B_FE_125_07_0334_0125_25', 'W_B_FE_125_07_0334_0125_40', 'W_B_FE_40_07_H7_0107_0125_12', 'W_B_FE_40_07_00_0107_0125_12', 'W_B_FE_40_07_H7_0107_0125_12', 'W_B_FE_40_07_H8_0107_0125_12', 'W_B_FE_40_07_0028_0125_12', 'W_B_FE_40_07_0214_0125_12', 'W_B_FE_40_00_0079_0125_12', 'W_B_FE_40_03_0082_0125_12', 'W_B_FE_40_085_0154_0125_12', 'W_B_FE_40_07_0060_0125_12', 'W_B_FE_40_07_0107_0125_12', 'W_B_FE_40_07_0107_0125_12', 'W_B_FE_40_07_0107_005_12', 'W_B_FE_40_07_0107_04_12'] Naben_Name = ['N_FE_40_08_11', 'N_FE_40_08_11', 'N_FE_20_08_11', 'N_FE_20_08_11', 'N_FE_80_08_11', 'N_FE_80_08_11', 'N_FE_125_08_11', 'N_FE_125_08_11', 'N_FE_125_05_11', 'N_FE_40_08_11', 'N_FE_40_08_11', 'N_FE_40_08_11', 'N_FE_40_08_11', 'N_FE_40_08_11', 'N_FE_40_08_11', 'N_FE_40_08_11', 'N_FE_40_08_11', 'N_FE_40_08_11', 'N_FE_40_05_11', 'N_FE_40_08_05', 'N_FE_40_08_15', 'N_FE_40_08_11', 'N_FE_40_08_11'] Welle_innen_Name = [' ', ' ', ' ', ' ', ' ', ' ', ' ', ' ', ' ', 'Wi_FE_28_n6','Wi_FE_28_00', 'Wi_FE_28_h6', 'Wi_FE_28_e8'] odb_Pfad = ['C:/temp/V-B_B1/WNV_B_FE_40_07_0125_12_08_11_75.odb', 'C:/temp/V-B_B1/WNV_B_FE_40_07_0125_06_08_11_75.odb', 'C:/temp/V-B_B1/WNV_B_FE_20_07_0125_12_08_11_75.odb', 'C:/temp/V-B_B1/WNV_B_FE_20_07_0125_06_08_11_75.odb', 'C:/temp/V-B_B1/WNV_B_FE_80_07_0125_06_08_11_75.odb', 'C:/temp/V-B_B1/WNV_B_FE_80_07_0125_16_08_11_75.odb', 'C:/temp/V-B_B1/WNV_B_FE_125_07_0125_12_08_11_75.odb', 'C:/temp/V-B_B1/WNV_B_FE_125_07_0125_25_08_11_75.odb', 'C:/temp/V-B_B1/WNV_B_FE_125_07_0125_40_08_11_75.odb', 'C:/temp/V-B_B1/WNV_Bi_FE_40_07_0125_12_08_11_75_H7n6.odb', 'C:/temp/V-B_B1/WNV_Bi_FE_40_07_0125_12_08_11_75_00.odb', 'C:/temp/V-B_B1/WNV_Bi_FE_40_07_0125_12_08_11_75_H7h6.odb', 'C:/temp/V-B_B1/WNV_Bi_FE_40_07_0125_12_08_11_75_H8e8.odb', 'C:/temp/V-B_B1/WNV_B_FE_40_07_0125_12_08_11_20.odb', 'C:/temp/V-B_B1/WNV_B_FE_40_07_0125_12_08_11_150.odb', 'C:/temp/V-B_B1/WNV_B_FE_40_00_0125_12_08_11_75.odb', 'C:/temp/V-B_B1/WNV_B_FE_40_03_0125_12_08_11_75.odb', 'C:/temp/V-B_B1/WNV_B_FE_40_085_0125_12_08_11_75.odb', 'C:/temp/V-B_B1/WNV_B_FE_40_07_0125_12_05_11_75.odb', 'C:/temp/V-B_B1/WNV_B_FE_40_07_0125_12_08_05_75.odb', 'C:/temp/V-B_B1/WNV_B_FE_40_07_0125_12_08_15_75.odb', 'C:/temp/V-B_B1/WNV_B_FE_40_07_005_12_08_11_75.odb', 'C:/temp/V-B_B1/WNV_B_FE_40_07_04_12_08_11_75.odb'] #WNV-Bezeichnungen ### Parameter ### # Wellenabmessung list_di = [28.0, 28.0, 14.0, 14.0, 56.0, 56.0, 87.5, 87.5, 87.5, 28.0, 28.0, 28.0, 28.0, 28.0, 28.0, 0.0, 12.0, 34.0, 28.0, 28.0, 28.0, 28.0, 28.0] # Innendurchmesser der Hohlwelle list_da = [40.0, 40.0, 20.0, 20.0, 80.0, 80.0, 125.0, 125.0, 125.0, 40.0, 40.0, 40.0, 40.0, 40.0, 40.0, 40.0, 40.0, 40.0, 40.0, 40.0, 40.0, 40.0, 40.0] # Außendurchmesser der Hohlwelle h_absatz = [5.0, 5.0, 2.5, 2.5, 10.0, 10.0, 15.625, 15.625, 15.625, 5.0, 5.0, 5.0, 5.0, 5.0, 5.0, 5.0, 5.0, 5.0, 5.0, 5.0, 5.0, 2.0, 16.0] # Absatzhöhe # Freistich t1 = [0.2, 0.3, 0.2, 0.3, 0.3, 0.3, 0.4, 0.4, 0.5, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2, 0.2] # Tiefe Freistich f = [2.5, 2.5, 2.5, 2.5, 2.5, 4.0, 4.0, 5.0, 7.0, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5, 2.5] # Länge Freistich radius = [1.2, 0.6, 1.2, 0.6, 0.6, 1.6, 1.2, 2.5, 4.0, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2, 1.2] # Radius im Freistich winkel = 15.0 # Winkel der Freistichschräge # Nabenabmessung list_Di = [40.0, 40.0, 20.0, 20.0, 80.0, 80.0, 125.0, 125.0, 125.0, 40.0, 40.0, 40.0, 40.0, 40.0, 40.0, 40.0, 40.0, 40.0, 40.0, 40.0, 40.0, 40.0, 40.0] # Innendurchmesser der Nabe list_Da = [50.0, 50.0, 25.0, 25.0, 100.0, 100.0, 156.25, 156.25, 156.25, 50.0, 50.0, 50.0, 50.0, 50.0, 50.0, 50.0, 50.0, 50.0, 60.0, 50.0, 50.0, 50.0, 50.0] # Außendurchmesser der Hohlwelle list_lf = [44.0, 44.0, 22.0, 22.0, 88.0, 88.0, 137.5, 137.5, 137.5, 44.0, 44.0, 44.0, 44.0, 44.0, 44.0, 44.0, 44.0, 44.0, 44.0, 20.0, 60.0, 44.0, 44.0] # Fugenlänge # Abmessung der innenliegende Welle da_i = 28.0 # Außendurchmesser der innenliegenden Welle # Materialkennwerte E = 210000.0 quer = 0.3 # Übermaß uebermass = [0.107, 0.107, 0.053, 0.053, 0.214, 0.214, 0.334, 0.334, 0.334, 0.107, 0.107, 0.107, 0.107, 0.028, 0.214, 0.079, 0.082, 0.154, 0.060, 0.107, 0.107, 0.107, 0.107] reibwert = 0.2 uebermass_wi = [0,0,0,0,0,0,0,0,0, 0.007, 0, 0.021, 0.073] # Netzfeinheit z = 0.125 # falls float (Kommazahl) dann folgende Schritte notwendig EI = int(8*z) ZW = int(16*z) DR = int (24*z) VI = int(32*z) FU = int(80*z) SE = int(160*z) for i in range (9,10): mdb.Model(name=Modell_Name[i], modelType=STANDARD_EXPLICIT) # mehrere Modelle erstellen session.viewports['Viewport: 1'].setValues(displayedObject=None) lw = 2*list_da[i]+list_lf[i]+22+f[i] # Länge der Welle l_absatz = list_da[i]+22 # Absatzlänge # Wellenmodellierung s = mdb.models[Modell_Name[i]].ConstrainedSketch(name='__profile__', sheetSize=200.0) g, v, d, c = s.geometry, s.vertices, s.dimensions, s.constraints s.setPrimaryObject(option=STANDALONE) s.ConstructionLine(point1=(0.0, -100.0), point2=(0.0, 100.0)) s.FixedConstraint(entity=g[2]) s.Line(point1=(list_di[i]/2, 0.0), point2=(list_da[i]/2+h_absatz[i], 0.0)) # Startpunkt Innendurchmesser zum Absatzdurchmesser s.HorizontalConstraint(entity=g[3], addUndoState=False) s.Line(point1=(list_da[i]/2+h_absatz[i], 0.0), point2=(list_da[i]/2+h_absatz[i], l_absatz)) # Absatzlinie s.VerticalConstraint(entity=g[4], addUndoState=False) s.PerpendicularConstraint(entity1=g[3], entity2=g[4], addUndoState=False) s.Line(point1=(list_da[i]/2+h_absatz[i], l_absatz), point2=(list_da[i]/2+(radius[i]-t1[i]), l_absatz)) # Linie von Absatz bis Start Radius am Freistich s.HorizontalConstraint(entity=g[5], addUndoState=False) s.PerpendicularConstraint(entity1=g[4], entity2=g[5], addUndoState=False) s.ArcByCenterEnds(center=(list_da[i]/2+(radius[i]-t1[i]), l_absatz+radius[i]), point1=(list_da[i]/2+(radius[i]-t1[i]), l_absatz), point2=(list_da[i]/2-t1[i], l_absatz+radius[i]), direction=CLOCKWISE) # erster Radius s.Line(point1=(list_da[i]/2-t1[i], l_absatz+radius[i]), point2=(list_da[i]/2-t1[i], (l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))-2*radius[i]*sin(winkel/2*pi/180))) # Horizontale Linie s.VerticalConstraint(entity=g[7], addUndoState=False) s.ArcByCenterEnds(center=(list_da[i]/2+(radius[i]-t1[i]), (l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin(winkel/2*pi/180))), point1=(list_da[i]/2-t1[i], (l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))-2*radius[i]*sin((winkel/2)*pi/180)), point2=(list_da[i]/2-(radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i])), (l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))), direction=CLOCKWISE) # zweiter Radius im Freistich s.Line(point1=(list_da[i]/2-(radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i])), (l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))), point2=(list_da[i]/2, l_absatz+f[i])) s.Line(point1=(list_da[i]/2, l_absatz+f[i]), point2=(list_da[i]/2, lw)) s.VerticalConstraint(entity=g[10], addUndoState=False) s.Line(point1=(list_da[i]/2, lw), point2=(list_di[i]/2, lw)) s.HorizontalConstraint(entity=g[11], addUndoState=False) s.PerpendicularConstraint(entity1=g[10], entity2=g[11], addUndoState=False) s.Line(point1=(list_di[i]/2, lw), point2=(list_di[i]/2, 0.0)) s.VerticalConstraint(entity=g[12], addUndoState=False) s.PerpendicularConstraint(entity1=g[11], entity2=g[12], addUndoState=False) if list_da[i] > 120: s.CoincidentConstraint(entity1=v[9], entity2=v[7]) welle = mdb.models[Modell_Name[i]].Part(name=Wellen_Name[i], dimensionality=THREE_D, type=DEFORMABLE_BODY) welle = mdb.models[Modell_Name[i]].parts[Wellen_Name[i]] welle.BaseSolidRevolve(sketch=s, angle=360.0, flipRevolveDirection=OFF) s.unsetPrimaryObject() welle = mdb.models[Modell_Name[i]].parts[Wellen_Name[i]] session.viewports['Viewport: 1'].setValues(displayedObject=welle) del mdb.models[Modell_Name[i]].sketches['__profile__'] # Nabenmodellierung s = mdb.models[Modell_Name[i]].ConstrainedSketch(name='__profile__', sheetSize=200.0) g, v, d, c = s.geometry, s.vertices, s.dimensions, s.constraints s.setPrimaryObject(option=STANDALONE) s.ConstructionLine(point1=(0.0, -100.0), point2=(0.0, 100.0)) s.FixedConstraint(entity=g[2]) s.rectangle(point1=(list_Di[i]/2, (l_absatz+1.5)), point2=(list_Da[i]/2, (l_absatz+f[i]+list_lf[i]))) # Nabenskizze nabe = mdb.models[Modell_Name[i]].Part(name=Naben_Name[i], dimensionality=THREE_D, type=DEFORMABLE_BODY) nabe = mdb.models[Modell_Name[i]].parts[Naben_Name[i]] nabe.BaseSolidRevolve(sketch=s, angle=360.0, flipRevolveDirection=OFF) # um wie viel Grad ich den Körper drehe s.unsetPrimaryObject() nabe = mdb.models[Modell_Name[i]].parts[Naben_Name[i]] session.viewports['Viewport: 1'].setValues(displayedObject=nabe) del mdb.models[Modell_Name[i]].sketches['__profile__'] # Innenliegende Welle if i > 8 and i < 13 : # Bedingung des WNV 10-13 mit innenliegender Welle sind s = mdb.models[Modell_Name[i]].ConstrainedSketch(name='__profile__', sheetSize=200.0) g, v, d, c = s.geometry, s.vertices, s.dimensions, s.constraints s.setPrimaryObject(option=STANDALONE) s.ConstructionLine(point1=(0.0, -100.0), point2=(0.0, 100.0)) s.FixedConstraint(entity=g[2]) s.rectangle(point1=(0.0, 0.0), point2=(da_i/2, lw)) # Abmessung der innenliegenden Welle welle_innen = mdb.models[Modell_Name[i]].Part(name=Welle_innen_Name[i], dimensionality=THREE_D, type=DEFORMABLE_BODY) welle_innen = mdb.models[Modell_Name[i]].parts[Welle_innen_Name[i]] welle_innen.BaseSolidRevolve(sketch=s, angle=360.0, flipRevolveDirection=OFF) s.unsetPrimaryObject() welle_innen = mdb.models[Modell_Name[i]].parts[Welle_innen_Name[i]] session.viewports['Viewport: 1'].setValues(displayedObject=welle_innen) del mdb.models[Modell_Name[i]].sketches['__profile__'] ### Partitionieren ### #Skizzenerstellung für Partitionierung der Stirnfläche#Welle welle = mdb.models[Modell_Name[i]].parts[Wellen_Name[i]] session.viewports['Viewport: 1'].setValues(displayedObject=welle) f_w, e_w, d1_w = welle.faces, welle.edges, welle.datums alpha =pi/4 Ebene_Hilfe_Welle = welle.DatumPlaneByPrincipalPlane(principalPlane=XZPLANE, offset=lw) #Hilfsebene Ebene_Hilfe_Welle = welle.datums[Ebene_Hilfe_Welle.id] Edge_Hilfe_Welle = e_w.findAt(((sin (alpha)*(list_da[i]/2),lw,cos (alpha)*(list_da[i]/2)),),) #Hilskante t = welle.MakeSketchTransform(sketchPlane=Ebene_Hilfe_Welle, sketchUpEdge=Edge_Hilfe_Welle[0], sketchPlaneSide=SIDE1, origin=(0,lw,0)) s = mdb.models[Modell_Name[i]].ConstrainedSketch( name='__profile__', sheetSize=256.32, gridSpacing=6.4, transform=t) g_w, v_w, d_w, c_w = s.geometry, s.vertices, s.dimensions, s.constraints s.setPrimaryObject(option=SUPERIMPOSE) welle = mdb.models[Modell_Name[i]].parts[Wellen_Name[i]] welle.projectReferencesOntoSketch(sketch=s, filter=COPLANAR_EDGES) s.CircleByCenterPerimeter(center=(0, 0), point1=(list_da[i]/2-1, 0)) #Zeichne Kreis welle = mdb.models[Modell_Name[i]].parts[Wellen_Name[i]] f_w = welle.faces pickedFaces = f_w.findAt(((sin (alpha)*(list_da[i]/2-1),lw,cos (alpha)*(list_da[i]/2-1)),),) e_w, d1_w = welle.edges, welle.datums welle.PartitionFaceBySketch(sketchUpEdge=Edge_Hilfe_Welle[0], faces=pickedFaces, sketch=s) s.unsetPrimaryObject() del mdb.models[Modell_Name[i]].sketches['__profile__'] #Sweepen der Skizze# welle = mdb.models[Modell_Name[i]].parts[Wellen_Name[i]] face = e_w.findAt(((sin (alpha)*(list_da[i]/2-1),lw,cos (alpha)*(list_da[i]/2-1)),),) # Auswahl der Faces cell_ges = welle.cells # über gesamtes Bauteil partitionieren e_w, d1_w = welle.edges, welle.datums pickedEdges = face welle.PartitionCellByExtrudeEdge(line=d1_w[1], cells=cell_ges, edges=pickedEdges, sense=FORWARD) # Sweepbefehl # Vierteln durch EbenenErstellung bei Welle # welle = mdb.models[Modell_Name[i]].parts[Wellen_Name[i]] session.viewports['Viewport: 1'].setValues(displayedObject=welle) Ebene_w1 = welle.DatumPlaneByPrincipalPlane(principalPlane=XYPLANE, offset=0.0) Ebene_w1 = welle.datums[Ebene_w1.id] cell_w_ges = welle.cells welle.PartitionCellByDatumPlane(datumPlane=Ebene_w1, cells=cell_w_ges) Ebene_w2 = welle.DatumPlaneByPrincipalPlane(principalPlane=YZPLANE, offset=0.0) Ebene_w2 = welle.datums[Ebene_w2.id] cell_w_ges = welle.cells welle.PartitionCellByDatumPlane(datumPlane=Ebene_w2, cells=cell_w_ges) # Partitionieren der Fugenflächen Ebene_w3 = welle.DatumPlaneByPrincipalPlane(principalPlane=XZPLANE, offset=(l_absatz-2)) # Kante vor dem Absatz Ebene_w3 = welle.datums[Ebene_w3.id] #cell_w_part = welle.cells.findAt(((-list_da[i]/2-0.01,0,0),),((list_da[i]/2-0.01,0,0),),) #welle.PartitionCellByDatumPlane(datumPlane=Ebene_w3, cells=cell_w_part) cell_w_ges = welle.cells welle.PartitionCellByDatumPlane(datumPlane=Ebene_w3, cells=cell_w_ges) welle = mdb.models[Modell_Name[i]].parts[Wellen_Name[i]] Ebene_w4 = welle.DatumPlaneByPrincipalPlane(principalPlane=XZPLANE, offset=(l_absatz+f[i]+list_lf[i])) Ebene_w4 = welle.datums[Ebene_w4.id] cell_w_part = welle.cells.findAt(((-list_di[i]/2-0.01,l_absatz+1.5+list_lf[i],0),),((list_di[i]/2+0.01,l_absatz+1.5+list_lf[i],0),),) welle.PartitionCellByDatumPlane(datumPlane=Ebene_w4, cells=cell_w_part) cell_w_ges = welle.cells welle.PartitionCellByDatumPlane(datumPlane=Ebene_w4, cells=cell_w_ges) welle = mdb.models[Modell_Name[i]].parts[Wellen_Name[i]] del welle.features['Partition cell-5'] # Partitionieren der Freistich Kanten Ebene_w5 = welle.DatumPlaneByPrincipalPlane(principalPlane=XZPLANE, offset=(l_absatz+f[i])) # Kante nach Freistich Ebene_w5 = welle.datums[Ebene_w5.id] cell_w_part = welle.cells.findAt(((-list_di[i]/2-0.01,l_absatz+f[i],0),),) welle.PartitionCellByDatumPlane(datumPlane=Ebene_w5, cells=cell_w_part) cell_w_ges = welle.cells welle.PartitionCellByDatumPlane(datumPlane=Ebene_w5, cells=cell_w_ges) welle = mdb.models[Modell_Name[i]].parts[Wellen_Name[i]] del welle.features['Partition cell-5'] Ebene_w6 = welle.DatumPlaneByPrincipalPlane(principalPlane=XZPLANE, offset=(l_absatz+radius[i])) # Kante nach dem erstem Radius Ebene_w6 = welle.datums[Ebene_w6.id] cell_w_part = welle.cells.findAt(((-list_di[i]/2-0.01,l_absatz+radius[i],0),),((list_di[i]/2+0.01,l_absatz+radius[i],0),),) welle.PartitionCellByDatumPlane(datumPlane=Ebene_w6, cells=cell_w_part) cell_w_ges = welle.cells welle.PartitionCellByDatumPlane(datumPlane=Ebene_w6, cells=cell_w_ges) welle = mdb.models[Modell_Name[i]].parts[Wellen_Name[i]] del welle.features['Partition cell-5'] Ebene_w7 = welle.DatumPlaneByPrincipalPlane(principalPlane=XZPLANE, offset=(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))) # Kante vor der Schräge des Freistich Ebene_w7 = welle.datums[Ebene_w7.id] cell_w_part = welle.cells.findAt(((-list_di[i]/2-0.01,(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),0),),) welle.PartitionCellByDatumPlane(datumPlane=Ebene_w7, cells=cell_w_part) cell_w_ges = welle.cells welle.PartitionCellByDatumPlane(datumPlane=Ebene_w7, cells=cell_w_ges) welle = mdb.models[Modell_Name[i]].parts[Wellen_Name[i]] del welle.features['Partition cell-5'] Ebene_w8 = welle.DatumPlaneByPrincipalPlane(principalPlane=XZPLANE, offset=(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180))) # Kante nach der Horizontalen Linie Ebene_w8 = welle.datums[Ebene_w8.id] cell_w_part = welle.cells.findAt(((-list_di[i]/2-0.01,(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),0),),) welle.PartitionCellByDatumPlane(datumPlane=Ebene_w8, cells=cell_w_part) cell_w_ges = welle.cells welle.PartitionCellByDatumPlane(datumPlane=Ebene_w8, cells=cell_w_ges) welle = mdb.models[Modell_Name[i]].parts[Wellen_Name[i]] del welle.features['Partition cell-5'] Ebene_w9 = welle.DatumPlaneByPrincipalPlane(principalPlane=XZPLANE, offset=(l_absatz-20)) # 2. Kante vor dem Absatz Ebene_w9 = welle.datums[Ebene_w9.id] cell_w_ges = welle.cells welle.PartitionCellByDatumPlane(datumPlane=Ebene_w9, cells=cell_w_ges) Ebene_w10 = welle.DatumPlaneByPrincipalPlane(principalPlane=XZPLANE, offset=(l_absatz+f[i]+4)) # Kante nach Freistich Ebene_w10 = welle.datums[Ebene_w10.id] cell_w_ges = welle.cells welle.PartitionCellByDatumPlane(datumPlane=Ebene_w10, cells=cell_w_ges) Ebene_w11 = welle.DatumPlaneByPrincipalPlane(principalPlane=XZPLANE, offset=(l_absatz+f[i]+list_lf[i]/2)) # Kante nach Freistich Ebene_w11 = welle.datums[Ebene_w11.id] cell_w_ges = welle.cells welle.PartitionCellByDatumPlane(datumPlane=Ebene_w11, cells=cell_w_ges) # Vierteln durch EbenenErstellung bei Nabe # Ebene_n1 = nabe.DatumPlaneByPrincipalPlane(principalPlane=YZPLANE, offset=0.0) Ebene_n1 = nabe.datums[Ebene_n1.id] Ebene_n2 = nabe.DatumPlaneByPrincipalPlane(principalPlane=XYPLANE, offset=0.0) Ebene_n2 = nabe.datums[Ebene_n2.id] cell_n_ges = nabe.cells nabe.PartitionCellByDatumPlane(datumPlane=Ebene_n1, cells=cell_n_ges) cell_n_ges = nabe.cells nabe.PartitionCellByDatumPlane(datumPlane=Ebene_n2, cells=cell_n_ges) Ebene_n3 = nabe.DatumPlaneByPrincipalPlane(principalPlane=XZPLANE, offset=l_absatz+f[i]) Ebene_n3 = nabe.datums[Ebene_n3.id] cell_n_ges = nabe.cells nabe.PartitionCellByDatumPlane(datumPlane=Ebene_n3, cells=cell_n_ges) Ebene_n4 = nabe.DatumPlaneByPrincipalPlane(principalPlane=XZPLANE, offset=l_absatz+f[i]+4) Ebene_n4 = nabe.datums[Ebene_n4.id] cell_n_ges = nabe.cells nabe.PartitionCellByDatumPlane(datumPlane=Ebene_n4, cells=cell_n_ges) Ebene_n5 = nabe.DatumPlaneByPrincipalPlane(principalPlane=XZPLANE, offset=l_absatz+f[i]+list_lf[i]/2) Ebene_n5 = nabe.datums[Ebene_n5.id] cell_n_ges = nabe.cells nabe.PartitionCellByDatumPlane(datumPlane=Ebene_n5, cells=cell_n_ges) # Vierteln durch EbenenErstellung bei innenliegenden Wellen # if i > 8 and i < 13 : # Bedingung des WNV 10-13 mit innenliegender Welle sind Ebene_wi1 = welle_innen.DatumPlaneByPrincipalPlane(principalPlane=YZPLANE, offset=0.0) Ebene_wi1 = welle_innen.datums[Ebene_wi1.id] Ebene_wi2 = welle_innen.DatumPlaneByPrincipalPlane(principalPlane=XYPLANE, offset=0.0) Ebene_wi2 = welle_innen.datums[Ebene_wi2.id] cell_wi_ges = welle_innen.cells welle_innen.PartitionCellByDatumPlane(datumPlane=Ebene_wi1, cells=cell_wi_ges) cell_wi_ges = welle_innen.cells welle_innen.PartitionCellByDatumPlane(datumPlane=Ebene_wi2, cells=cell_wi_ges) # Partitionieren der Fugenflächen Ebene_wi3 = welle_innen.DatumPlaneByPrincipalPlane(principalPlane=XZPLANE, offset=(l_absatz-2)) Ebene_wi3 = welle_innen.datums[Ebene_wi3.id] cell_wi_ges = welle_innen.cells welle_innen.PartitionCellByDatumPlane(datumPlane=Ebene_wi3, cells=cell_wi_ges) Ebene_wi4 = welle_innen.DatumPlaneByPrincipalPlane(principalPlane=XZPLANE, offset=(l_absatz+f[i]+list_lf[i])) Ebene_wi4 = welle_innen.datums[Ebene_wi4.id] cell_wi_ges = welle_innen.cells welle_innen.PartitionCellByDatumPlane(datumPlane=Ebene_wi4, cells=cell_wi_ges) # Partitionieren restlicher Kanten Ebene_wi5 = welle_innen.DatumPlaneByPrincipalPlane(principalPlane=XZPLANE, offset=(l_absatz+f[i])) # Kante nach Freistich Ebene_wi5 = welle_innen.datums[Ebene_wi5.id] cell_wi_ges = welle_innen.cells welle_innen.PartitionCellByDatumPlane(datumPlane=Ebene_wi5, cells=cell_wi_ges) Ebene_wi6 = welle_innen.DatumPlaneByPrincipalPlane(principalPlane=XZPLANE, offset=(l_absatz+radius[i])) # Kante nach dem ersten Radius Ebene_wi6 = welle_innen.datums[Ebene_wi6.id] cell_wi_ges = welle_innen.cells welle_innen.PartitionCellByDatumPlane(datumPlane=Ebene_wi6, cells=cell_wi_ges) Ebene_wi7 = welle_innen.DatumPlaneByPrincipalPlane(principalPlane=XZPLANE, offset=(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))) # Kante vor der Schräge des Freistich Ebene_wi7 = welle_innen.datums[Ebene_wi7.id] cell_wi_ges = welle_innen.cells welle_innen.PartitionCellByDatumPlane(datumPlane=Ebene_wi7, cells=cell_wi_ges) Ebene_wi8 = welle_innen.DatumPlaneByPrincipalPlane(principalPlane=XZPLANE, offset=(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180))) # Kante nach der Horizontalen Linie Ebene_wi8 = welle_innen.datums[Ebene_wi8.id] cell_wi_ges = welle_innen.cells welle_innen.PartitionCellByDatumPlane(datumPlane=Ebene_wi8, cells=cell_wi_ges) Ebene_wi9 = welle_innen.DatumPlaneByPrincipalPlane(principalPlane=XZPLANE, offset=(l_absatz-20)) # 2. Kante vor dem Absatz Ebene_wi9 = welle_innen.datums[Ebene_wi9.id] cell_wi_ges = welle_innen.cells welle_innen.PartitionCellByDatumPlane(datumPlane=Ebene_wi9, cells=cell_wi_ges) Ebene_wi10 = welle_innen.DatumPlaneByPrincipalPlane(principalPlane=XZPLANE, offset=(l_absatz+f[i]+4)) # 2. Kante vor dem Absatz Ebene_wi10 = welle_innen.datums[Ebene_wi10.id] cell_wi_ges = welle_innen.cells welle_innen.PartitionCellByDatumPlane(datumPlane=Ebene_wi10, cells=cell_wi_ges) Ebene_wi11 = welle_innen.DatumPlaneByPrincipalPlane(principalPlane=XZPLANE, offset=(l_absatz+f[i]+list_lf[i]/2)) # 2. Kante vor dem Absatz Ebene_wi11 = welle_innen.datums[Ebene_wi11.id] cell_wi_ges = welle_innen.cells welle_innen.PartitionCellByDatumPlane(datumPlane=Ebene_wi11, cells=cell_wi_ges) #Skizzenerstellung für Partitionierung der Absatzfläche#Welle welle = mdb.models[Modell_Name[i]].parts[Wellen_Name[i]] session.viewports['Viewport: 1'].setValues(displayedObject=welle) f_w, e_w, d1_w = welle.faces, welle.edges, welle.datums alpha =pi/4 Ebene_Hilfe_Welle = welle.DatumPlaneByPrincipalPlane(principalPlane=XZPLANE, offset=l_absatz) #Hilfsebene Ebene_Hilfe_Welle = welle.datums[Ebene_Hilfe_Welle.id] Edge_Hilfe_Welle = e_w.findAt(((sin (alpha)*(list_da[i]/2+h_absatz[i]),l_absatz,cos (alpha)*(list_da[i]/2+h_absatz[i])),),) #Hilskante t = welle.MakeSketchTransform(sketchPlane=Ebene_Hilfe_Welle, sketchUpEdge=Edge_Hilfe_Welle[0], sketchPlaneSide=SIDE1, origin=(0,l_absatz,0)) s = mdb.models[Modell_Name[i]].ConstrainedSketch( name='__profile__', sheetSize=256.32, gridSpacing=6.4, transform=t) g_w, v_w, d_w, c_w = s.geometry, s.vertices, s.dimensions, s.constraints s.setPrimaryObject(option=SUPERIMPOSE) welle = mdb.models[Modell_Name[i]].parts[Wellen_Name[i]] welle.projectReferencesOntoSketch(sketch=s, filter=COPLANAR_EDGES) s.CircleByCenterPerimeter(center=(0, 0), point1=(0, list_da[i]/2+radius[i])) #Zeichne Kreis welle = mdb.models[Modell_Name[i]].parts[Wellen_Name[i]] f_w = welle.faces pickedFaces = [0,0,0,0] for k in range (0,4): # Diese Schleife soll die Faces in jedem Viertel erzeugen welle = mdb.models[Modell_Name[i]].parts[Wellen_Name[i]] Edge_Hilfe_Welle = e_w.findAt(((sin (alpha)*(list_da[i]/2+h_absatz[i]),l_absatz,cos (alpha)*(list_da[i]/2+h_absatz[i])),),) #Hilskante f_w = welle.faces pickedFaces[k] = f_w.findAt(((sin (alpha)*(list_da[i]/2+radius[i]),l_absatz,cos (alpha)*(list_da[i]/2+radius[i])),),) e_w, d1_w = welle.edges, welle.datums welle.PartitionFaceBySketch(sketchUpEdge=Edge_Hilfe_Welle[0], faces=pickedFaces[k], sketch=s) alpha = alpha + pi/2 s.unsetPrimaryObject() del mdb.models[Modell_Name[i]].sketches['__profile__'] #Sweepen der Skizze# alpha = pi/4 face = [0,0,0,0] # ist denke notwendig sonst fehlt Definition in der Schleife for j in range (0,4): welle = mdb.models[Modell_Name[i]].parts[Wellen_Name[i]] face[j] = e_w.findAt(((sin (alpha)*(list_da[i]/2+radius[i]),l_absatz,cos (alpha)*(list_da[i]/2+radius[i])),),) # Auswahl der Faces cell_ges = welle.cells # über gesamtes Bauteil partitionieren e_w, d1_w = welle.edges, welle.datums pickedEdges = face[j] welle.PartitionCellByExtrudeEdge(line=d1_w[1], cells=cell_ges, edges=pickedEdges, sense=FORWARD) # Sweepbefehl alpha = alpha + pi/2 #Skizzenerstellung für Partitionierung der Stirnfläche#Nabe nabe = mdb.models[Modell_Name[i]].parts[Naben_Name[i]] session.viewports['Viewport: 1'].setValues(displayedObject=nabe) f_n, e_n, d1_n = nabe.faces, nabe.edges, nabe.datums beta= pi/4 Ebene_Hilfe_Nabe = nabe.DatumPlaneByPrincipalPlane(principalPlane=XZPLANE, offset=(l_absatz+f[i]+list_lf[i])) #Hilfsebene Ebene_Hilfe_Nabe = nabe.datums[Ebene_Hilfe_Nabe.id] Edge_Hilfe_Nabe = e_n.findAt(((sin (beta)*(list_Di[i]/2),(l_absatz+f[i]+list_lf[i]),cos (beta)*(list_Di[i]/2)),),) #Hilskante t = nabe.MakeSketchTransform(sketchPlane=Ebene_Hilfe_Nabe, sketchUpEdge=Edge_Hilfe_Nabe[0], sketchPlaneSide=SIDE1, origin=(0,l_absatz+f[i]+list_lf[i],0)) s = mdb.models[Modell_Name[i]].ConstrainedSketch( name='__profile__', sheetSize=256.32, gridSpacing=6.4, transform=t) g_n, v_n, d_n, c_n = s.geometry, s.vertices, s.dimensions, s.constraints s.setPrimaryObject(option=SUPERIMPOSE) nabe = mdb.models[Modell_Name[i]].parts[Naben_Name[i]] nabe.projectReferencesOntoSketch(sketch=s, filter=COPLANAR_EDGES) s.CircleByCenterPerimeter(center=(0, 0), point1=(0, list_Di[i]/2+1)) #Zeichnet Kreis pickedFaces = [0,0,0,0] for k in range (0,4): # Diese Schleife soll die Faces in jedem Viertel erzeugen nabe = mdb.models[Modell_Name[i]].parts[Naben_Name[i]] Edge_Hilfe_Nabe = e_n.findAt(((sin (beta)*(list_Di[i]/2),(l_absatz+f[i]+list_lf[i]),cos (beta)*(list_Di[i]/2)),),) #Hilskante f_n = nabe.faces pickedFaces = f_n.findAt(((sin (beta)*(list_Di[i]/2+1),l_absatz+f[i]+list_lf[i],cos (beta)*(list_Di[i]/2+1)),),) e_n, d1_n = nabe.edges, nabe.datums nabe.PartitionFaceBySketch(sketchUpEdge=Edge_Hilfe_Nabe[0], faces=pickedFaces, sketch=s) beta = beta + pi/2 s.unsetPrimaryObject() del mdb.models[Modell_Name[i]].sketches['__profile__'] #Sweepen der Skizze# alpha = pi/4 face = [0,0,0,0] # ist denke notwendig sonst fehlt Definition in der Schleife for j in range (0,4): nabe = mdb.models[Modell_Name[i]].parts[Naben_Name[i]] face[j] = e_n.findAt(((sin (alpha)*(list_Di[i]/2+1),l_absatz+f[i]+list_lf[i],cos (alpha)*(list_Di[i]/2+1)),),) # Auswahl der Faces cell_ges = nabe.cells # über gesamtes Bauteil partitionieren e_n, d1_n = nabe.edges, nabe.datums pickedEdges = face[j] nabe.PartitionCellByExtrudeEdge(line=d1_n[1], cells=cell_ges, edges=pickedEdges, sense=FORWARD) # Sweepbefehl alpha = alpha + pi/2 #Skizzenerstellung für Partitionierung der Stirnfläche#innenliegende Welle if i > 8 and i < 13: # Bedingung des WNV 10-13 mit innenliegender Welle sind welle_innen = mdb.models[Modell_Name[i]].parts[Welle_innen_Name[i]] session.viewports['Viewport: 1'].setValues(displayedObject=welle_innen) f_wi, e_wi, d1_wi = welle_innen.faces, welle_innen.edges, welle_innen.datums beta= pi/4 Ebene_Hilfe_Welle_innen = welle_innen.DatumPlaneByPrincipalPlane(principalPlane=XZPLANE, offset=lw) #Hilfsebene Ebene_Hilfe_Welle_innen = welle_innen.datums[Ebene_Hilfe_Welle_innen.id] Edge_Hilfe_Welle_innen = e_wi.findAt(((sin (beta)*(da_i/2),lw,cos (beta)*(da_i/2)),),) #Hilskante t = welle_innen.MakeSketchTransform(sketchPlane=Ebene_Hilfe_Welle_innen, sketchUpEdge=Edge_Hilfe_Welle_innen[0], sketchPlaneSide=SIDE1, origin=(0,lw,0)) s = mdb.models[Modell_Name[i]].ConstrainedSketch( name='__profile__', sheetSize=256.32, gridSpacing=6.4, transform=t) g_wi, v_wi, d_wi, c_wi = s.geometry, s.vertices, s.dimensions, s.constraints s.setPrimaryObject(option=SUPERIMPOSE) welle_innen = mdb.models[Modell_Name[i]].parts[Welle_innen_Name[i]] welle_innen.projectReferencesOntoSketch(sketch=s, filter=COPLANAR_EDGES) s.CircleByCenterPerimeter(center=(0, 0), point1=(0, da_i/2-1)) #Zeichnet Kreis pickedFaces = [0,0,0,0] for k in range (0,4): # Diese Schleife soll die Faces in jedem Viertel erzeugen welle_innen = mdb.models[Modell_Name[i]].parts[Welle_innen_Name[i]] Edge_Hilfe_Welle_innen = e_wi.findAt(((sin (beta)*(da_i/2),lw,cos (beta)*(da_i/2)),),) #Hilskante f_wi = welle_innen.faces pickedFaces[k] = f_wi.findAt(((sin (beta)*(da_i/2-1),lw,cos (beta)*(da_i/2-1)),),) e_wi, d1_wi = welle_innen.edges, welle_innen.datums welle_innen.PartitionFaceBySketch(sketchUpEdge=Edge_Hilfe_Welle_innen[0], faces=pickedFaces[k], sketch=s) beta = beta + pi/2 s.unsetPrimaryObject() del mdb.models[Modell_Name[i]].sketches['__profile__'] #Sweepen der Skizze# alpha = pi/4 face = [0,0,0,0] # ist denke notwendig sonst fehlt Definition in der Schleife for j in range (0,4): welle_innen = mdb.models[Modell_Name[i]].parts[Welle_innen_Name[i]] face[j] = e_wi.findAt(((sin (alpha)*(da_i/2-1),lw,cos (alpha)*(da_i/2-1)),),) # Auswahl der Faces cell_ges = welle_innen.cells # über gesamtes Bauteil partitionieren e_wi, d1_wi = welle_innen.edges, welle_innen.datums pickedEdges = face[j] welle_innen.PartitionCellByExtrudeEdge(line=d1_wi[1], cells=cell_ges, edges=pickedEdges, sense=FORWARD) # Sweepbefehl alpha = alpha + pi/2 # Materialzuweisung und Erstellen von Section mdb.models[Modell_Name[i]].Material(name='Stahl') # Materialname vergeben mdb.models[Modell_Name[i]].materials['Stahl'].Elastic(table=((E, quer), )) # Materialkennwerte definieren hier: elastisches Verhalten mdb.models[Modell_Name[i]].HomogeneousSolidSection(name='Stahl', material='Stahl', # Section definieren thickness=None) nabe = mdb.models[Modell_Name[i]].parts[Naben_Name[i]] # Nabe Material zuweisen cell_ges = nabe.cells nabe.Set(cells=cell_ges, name='Alles') #Als Setname bietet sich soetwas wie: Alles, Wholepart, Gesamt oder sowas an :) region = nabe.sets['Alles'] nabe.SectionAssignment(region=region, sectionName='Stahl', #SectionAssignment Zuweisung Section der Parts offset=0.0, offsetType=MIDDLE_SURFACE, offsetField='', thicknessAssignment=FROM_SECTION) welle = mdb.models[Modell_Name[i]].parts[Wellen_Name[i]] # Welle Material zuweisen session.viewports['Viewport: 1'].setValues(displayedObject=welle) cell_ges = welle.cells welle.Set(cells=cell_ges, name='Alles') #Als Setname bietet sich soetwas wie: Alles, Wholepart, Gesamt oder sowas an :) region = welle.sets['Alles'] welle.SectionAssignment(region=region, sectionName='Stahl', #SectionAssignment Zuweisung Section der Parts offset=0.0, offsetType=MIDDLE_SURFACE, offsetField='', thicknessAssignment=FROM_SECTION) if i > 8 and i < 13 : # Bedingung des WNV 10-13 mit innenliegender Welle sind welle_innen = mdb.models[Modell_Name[i]].parts[Welle_innen_Name[i]] # innenliegende Welle Material zuweisen session.viewports['Viewport: 1'].setValues(displayedObject=welle_innen) cell_ges = welle_innen.cells welle_innen.Set(cells=cell_ges, name='Alles') #Als Setname bietet sich soetwas wie: Alles, Wholepart, Gesamt oder sowas an :) region = welle_innen.sets['Alles'] welle_innen.SectionAssignment(region=region, sectionName='Stahl', #SectionAssignment Zuweisung Section der Parts offset=0.0, offsetType=MIDDLE_SURFACE, offsetField='', thicknessAssignment=FROM_SECTION) # Baugruppen Erstellung wnv = mdb.models[Modell_Name[i]].rootAssembly #Baugruppe erstellen wnv.DatumCsysByDefault(CARTESIAN) nabe = mdb.models[Modell_Name[i]].parts[Naben_Name[i]] # Nabe auswählen wnv.Instance(name=Naben_Name[i], part=nabe, dependent=ON) welle = mdb.models[Modell_Name[i]].parts[Wellen_Name[i]] # Welle auswählen wnv.Instance(name=Wellen_Name[i], part=welle, dependent=ON) if i > 8 and i < 13 : welle_innen = mdb.models[Modell_Name[i]].parts[Welle_innen_Name[i]] # Welle auswählen wnv.Instance(name=Welle_innen_Name[i], part=welle_innen, dependent=ON) # Step erstellen # Querpressverbund mdb.models[Modell_Name[i]].StaticStep(name='Querpressverbund', previous='Initial', #s.o. für vorherigen Stepname maxNumInc=100, #Maximale rechenincrementanzahl initialInc=0.1, #incrementgröße des ersten Rechenschritts minInc=1e-05, #kleinste Incrementgröße maxInc=1.0, #größte Incrementgröße nlgeom=ON) #Berechnen mit Nichtlinearen Operationen ON = an OFF = aus duh. # Reibung mdb.models[Modell_Name[i]].StaticStep(name='Reibung', previous='Querpressverbund', #s.o. für vorherigen Stepname maxNumInc=100, #Maximale rechenincrementanzahl initialInc=0.1, #incrementgröße des ersten Rechenschritts minInc=1e-05, #kleinste Incrementgröße maxInc=1.0, #größte Incrementgröße nlgeom=ON) #Berechnen mit Nichtlinearen Operationen ON = an OFF = aus duh. # Erstellen von Surface # Welleneinspannung welle = mdb.models[Modell_Name[i]].parts[Wellen_Name[i]] session.viewports['Viewport: 1'].setValues(displayedObject=welle) v_w = welle.vertices Einspannung_Welle = v_w.findAt(((list_da[i]/2, l_absatz+f[i]+list_lf[i]/2, 0),),) welle.Set(vertices=Einspannung_Welle, name='Einspannung Welle') # Erstellen von Surface # Wellenkontaktfläche Kontaktflaeche1_PV = f_w.findAt(((sin(pi/4)*list_da[i]/2,l_absatz+f[i]+2,cos(pi/4)*list_da[i]/2),),) # nicht genau die erzwugte Kante auswählen sonst gibt es Probleme Kontaktflaeche2_PV = f_w.findAt(((sin(3*pi/4)*list_da[i]/2,l_absatz+f[i]+2,cos(3*pi/4)*list_da[i]/2),),) Kontaktflaeche3_PV = f_w.findAt(((sin(5*pi/4)*list_da[i]/2,l_absatz+f[i]+2,cos(5*pi/4)*list_da[i]/2),),) Kontaktflaeche4_PV = f_w.findAt(((sin(7*pi/4)*list_da[i]/2,l_absatz+f[i]+2,cos(7*pi/4)*list_da[i]/2),),) Kontaktflaeche5_PV = f_w.findAt(((sin(pi/4)*list_da[i]/2,l_absatz+f[i]+5,cos(pi/4)*list_da[i]/2),),) # nicht genau die erzwugte Kante auswählen sonst gibt es Probleme Kontaktflaeche6_PV = f_w.findAt(((sin(3*pi/4)*list_da[i]/2,l_absatz+f[i]+5,cos(3*pi/4)*list_da[i]/2),),) Kontaktflaeche7_PV = f_w.findAt(((sin(5*pi/4)*list_da[i]/2,l_absatz+f[i]+5,cos(5*pi/4)*list_da[i]/2),),) Kontaktflaeche8_PV = f_w.findAt(((sin(7*pi/4)*list_da[i]/2,l_absatz+f[i]+5,cos(7*pi/4)*list_da[i]/2),),) Kontaktflaeche_PV_Welle = [Kontaktflaeche1_PV, Kontaktflaeche2_PV, Kontaktflaeche3_PV, Kontaktflaeche4_PV, Kontaktflaeche5_PV, Kontaktflaeche6_PV, Kontaktflaeche7_PV, Kontaktflaeche8_PV] welle.Surface(side1Faces=Kontaktflaeche_PV_Welle, name='Kontaktflaeche PV') # Erstellen von Surface # Wellenkontaktfläche für innenliegende Welle if i > 8 and i < 13: Kontaktflaeche1_inne = f_w.findAt(((sin(pi/4)*list_di[i]/2,l_absatz-3,cos(pi/4)*list_di[i]/2),),) # nicht genau die erzwugte Kante auswählen sonst gibt es Probleme Kontaktflaeche2_inne = f_w.findAt(((sin(3*pi/4)*list_di[i]/2,l_absatz-3,cos(3*pi/4)*list_di[i]/2),),) Kontaktflaeche3_inne = f_w.findAt(((sin(5*pi/4)*list_di[i]/2,l_absatz-3,cos(5*pi/4)*list_di[i]/2),),) Kontaktflaeche4_inne = f_w.findAt(((sin(7*pi/4)*list_di[i]/2,l_absatz-3,cos(7*pi/4)*list_di[i]/2),),) Kontaktflaeche5_inne = f_w.findAt(((sin(pi/4)*list_di[i]/2,l_absatz,cos(pi/4)*list_di[i]/2),),) Kontaktflaeche6_inne = f_w.findAt(((sin(3*pi/4)*list_di[i]/2,l_absatz,cos(3*pi/4)*list_di[i]/2),),) Kontaktflaeche7_inne = f_w.findAt(((sin(5*pi/4)*list_di[i]/2,l_absatz,cos(5*pi/4)*list_di[i]/2),),) Kontaktflaeche8_inne = f_w.findAt(((sin(7*pi/4)*list_di[i]/2,l_absatz,cos(7*pi/4)*list_di[i]/2),),) Kontaktflaeche9_inne = f_w.findAt(((sin(pi/4)*list_di[i]/2,l_absatz+radius[i]+0.1,cos(pi/4)*list_di[i]/2),),) # nicht genau die erzwugte Kante auswählen sonst gibt es Probleme Kontaktflaeche10_inne = f_w.findAt(((sin(3*pi/4)*list_di[i]/2,l_absatz+radius[i]+0.1,cos(3*pi/4)*list_di[i]/2),),) # Flächen der horizontaalen Linie im Freistich Kontaktflaeche11_inne = f_w.findAt(((sin(5*pi/4)*list_di[i]/2,l_absatz+radius[i]+0.1,cos(5*pi/4)*list_di[i]/2),),) Kontaktflaeche12_inne = f_w.findAt(((sin(7*pi/4)*list_di[i]/2,l_absatz+radius[i]+0.1,cos(7*pi/4)*list_di[i]/2),),) Kontaktflaeche13_inne = f_w.findAt(((sin(pi/4)*list_di[i]/2,(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))-0.05,cos(pi/4)*list_di[i]/2),),) # nicht genau die erzwugte Kante auswählen sonst gibt es Probleme Kontaktflaeche14_inne = f_w.findAt(((sin(3*pi/4)*list_di[i]/2,(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))-0.05,cos(3*pi/4)*list_di[i]/2),),) # Fläche bei 2. Rundung im Freistich Kontaktflaeche15_inne = f_w.findAt(((sin(5*pi/4)*list_di[i]/2,(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))-0.05,cos(5*pi/4)*list_di[i]/2),),) Kontaktflaeche16_inne = f_w.findAt(((sin(7*pi/4)*list_di[i]/2,(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))-0.05,cos(7*pi/4)*list_di[i]/2),),) Kontaktflaeche17_inne = f_w.findAt(((sin(pi/4)*list_di[i]/2,l_absatz+f[i]-0.1,cos(pi/4)*list_di[i]/2),),) # nicht genau die erzwugte Kante auswählen sonst gibt es Probleme Kontaktflaeche18_inne = f_w.findAt(((sin(3*pi/4)*list_di[i]/2,l_absatz+f[i]-0.1,cos(3*pi/4)*list_di[i]/2),),) # Flächen der schrägen Kante Kontaktflaeche19_inne = f_w.findAt(((sin(5*pi/4)*list_di[i]/2,l_absatz+f[i]-0.1,cos(5*pi/4)*list_di[i]/2),),) Kontaktflaeche20_inne = f_w.findAt(((sin(7*pi/4)*list_di[i]/2,l_absatz+f[i]-0.1,cos(7*pi/4)*list_di[i]/2),),) Kontaktflaeche21_inne = f_w.findAt(((sin(pi/4)*list_di[i]/2,l_absatz+f[i]+3,cos(pi/4)*list_di[i]/2),),) # nicht genau die erzwugte Kante auswählen sonst gibt es Probleme Kontaktflaeche22_inne = f_w.findAt(((sin(3*pi/4)*list_di[i]/2,l_absatz+f[i]+3,cos(3*pi/4)*list_di[i]/2),),) # Fugenflächen innen anfang Kontaktflaeche23_inne = f_w.findAt(((sin(5*pi/4)*list_di[i]/2,l_absatz+f[i]+3,cos(5*pi/4)*list_di[i]/2),),) Kontaktflaeche24_inne = f_w.findAt(((sin(7*pi/4)*list_di[i]/2,l_absatz+f[i]+3,cos(7*pi/4)*list_di[i]/2),),) Kontaktflaeche25_inne = f_w.findAt(((sin(pi/4)*list_di[i]/2,lw-2,cos(pi/4)*list_di[i]/2),),) # nicht genau die erzwugte Kante auswählen sonst gibt es Probleme Kontaktflaeche26_inne = f_w.findAt(((sin(3*pi/4)*list_di[i]/2,lw-2,cos(3*pi/4)*list_di[i]/2),),) Kontaktflaeche27_inne = f_w.findAt(((sin(5*pi/4)*list_di[i]/2,lw-2,cos(5*pi/4)*list_di[i]/2),),) Kontaktflaeche28_inne = f_w.findAt(((sin(7*pi/4)*list_di[i]/2,lw-2,cos(7*pi/4)*list_di[i]/2),),) Kontaktflaeche29_inne = f_w.findAt(((sin(pi/4)*list_di[i]/2,l_absatz+f[i]+6,cos(pi/4)*list_di[i]/2),),) # nicht genau die erzwugte Kante auswählen sonst gibt es Probleme Kontaktflaeche30_inne = f_w.findAt(((sin(3*pi/4)*list_di[i]/2,l_absatz+f[i]+6,cos(3*pi/4)*list_di[i]/2),),) # Fugenflächen innen Kontaktflaeche31_inne = f_w.findAt(((sin(5*pi/4)*list_di[i]/2,l_absatz+f[i]+6,cos(5*pi/4)*list_di[i]/2),),) Kontaktflaeche32_inne = f_w.findAt(((sin(7*pi/4)*list_di[i]/2,l_absatz+f[i]+6,cos(7*pi/4)*list_di[i]/2),),) Kontaktflaeche33_inne = f_w.findAt(((sin(pi/4)*list_di[i]/2, 2,cos(pi/4)*list_di[i]/2),),) # nicht genau die erzwugte Kante auswählen sonst gibt es Probleme Kontaktflaeche34_inne = f_w.findAt(((sin(3*pi/4)*list_di[i]/2, 2,cos(3*pi/4)*list_di[i]/2),),) Kontaktflaeche35_inne = f_w.findAt(((sin(5*pi/4)*list_di[i]/2, 2,cos(5*pi/4)*list_di[i]/2),),) Kontaktflaeche36_inne = f_w.findAt(((sin(7*pi/4)*list_di[i]/2, 2,cos(7*pi/4)*list_di[i]/2),),) Kontaktflaeche_inne_Welle = [Kontaktflaeche1_inne, Kontaktflaeche2_inne, Kontaktflaeche3_inne, Kontaktflaeche4_inne, Kontaktflaeche5_inne, Kontaktflaeche6_inne, Kontaktflaeche7_inne, Kontaktflaeche8_inne, Kontaktflaeche9_inne, Kontaktflaeche10_inne, Kontaktflaeche11_inne, Kontaktflaeche12_inne, Kontaktflaeche13_inne, Kontaktflaeche14_inne, Kontaktflaeche15_inne, Kontaktflaeche16_inne, Kontaktflaeche17_inne, Kontaktflaeche18_inne, Kontaktflaeche19_inne, Kontaktflaeche20_inne, Kontaktflaeche21_inne, Kontaktflaeche22_inne, Kontaktflaeche23_inne, Kontaktflaeche24_inne, Kontaktflaeche25_inne, Kontaktflaeche26_inne, Kontaktflaeche27_inne, Kontaktflaeche28_inne, Kontaktflaeche29_inne, Kontaktflaeche30_inne, Kontaktflaeche31_inne, Kontaktflaeche32_inne, Kontaktflaeche33_inne, Kontaktflaeche34_inne, Kontaktflaeche35_inne, Kontaktflaeche36_inne] welle.Surface(side1Faces=Kontaktflaeche_inne_Welle, name='Kontaktflaeche Innen') # Erstellen von Surface # Nabeneinspannung nabe = mdb.models[Modell_Name[i]].parts[Naben_Name[i]] session.viewports['Viewport: 1'].setValues(displayedObject=nabe) v_n = nabe.vertices Einspannung_Nabe = v_n.findAt(((list_Da[i]/2, l_absatz+f[i]+list_lf[i]/2, 0),),) nabe.Set(vertices=Einspannung_Nabe, name='Einspannung Nabe') # Erstellen von Surface # Nabenkontaktfläche Kontaktflaeche1_PV = f_n.findAt(((sin(pi/4)*list_Di[i]/2,l_absatz+f[i]+2 ,cos(pi/4)*list_Di[i]/2),),) # nicht genau die erzwugte Kante auswählen sonst gibt es Probleme Kontaktflaeche2_PV = f_n.findAt(((sin(3*pi/4)*list_Di[i]/2,l_absatz+f[i]+2,cos(3*pi/4)*list_Di[i]/2),),) Kontaktflaeche3_PV = f_n.findAt(((sin(5*pi/4)*list_Di[i]/2,l_absatz+f[i]+2,cos(5*pi/4)*list_Di[i]/2),),) Kontaktflaeche4_PV = f_n.findAt(((sin(7*pi/4)*list_Di[i]/2,l_absatz+f[i]+2,cos(7*pi/4)*list_Di[i]/2),),) Kontaktflaeche5_PV = f_n.findAt(((sin(pi/4)*list_Di[i]/2,l_absatz+f[i]+6 ,cos(pi/4)*list_Di[i]/2),),) # nicht genau die erzwugte Kante auswählen sonst gibt es Probleme Kontaktflaeche6_PV = f_n.findAt(((sin(3*pi/4)*list_Di[i]/2,l_absatz+f[i]+6,cos(3*pi/4)*list_Di[i]/2),),) Kontaktflaeche7_PV = f_n.findAt(((sin(5*pi/4)*list_Di[i]/2,l_absatz+f[i]+6,cos(5*pi/4)*list_Di[i]/2),),) Kontaktflaeche8_PV = f_n.findAt(((sin(7*pi/4)*list_Di[i]/2,l_absatz+f[i]+6,cos(7*pi/4)*list_Di[i]/2),),) Kontaktflaeche_PV_Nabe = [Kontaktflaeche1_PV, Kontaktflaeche2_PV, Kontaktflaeche3_PV, Kontaktflaeche4_PV, Kontaktflaeche5_PV, Kontaktflaeche6_PV, Kontaktflaeche7_PV, Kontaktflaeche8_PV] nabe.Surface(side1Faces=Kontaktflaeche_PV_Nabe, name='Kontaktflaeche PV') # Erstellen von Surface # innenliegende Welleneinspannung if i > 8 and i < 13: welle_innen = mdb.models[Modell_Name[i]].parts[Welle_innen_Name[i]] session.viewports['Viewport: 1'].setValues(displayedObject=welle_innen) v_wi = welle_innen.vertices Einspannung_Welle_innen = v_wi.findAt(((da_i/2, l_absatz+f[i]+list_lf[i]/2, 0),),) welle_innen.Set(vertices=Einspannung_Welle_innen, name='Einspannung Welle innen') # Erstellen von Surface # Wellenkontaktfläche für innenliegende Welle Kontaktflaeche1 = f_wi.findAt(((sin(pi/4)*da_i/2,2,cos(pi/4)*da_i/2),),) # nicht genau die erzweugte Kante auswählen sonst gibt es Probleme Kontaktflaeche2 = f_wi.findAt(((sin(3*pi/4)*da_i/2,2,cos(3*pi/4)*da_i/2),),) Kontaktflaeche3 = f_wi.findAt(((sin(5*pi/4)*da_i/2,2,cos(5*pi/4)*da_i/2),),) Kontaktflaeche4 = f_wi.findAt(((sin(7*pi/4)*da_i/2,2,cos(7*pi/4)*da_i/2),),) Kontaktflaeche5 = f_wi.findAt(((sin(pi/4)*da_i/2,l_absatz+0.2,cos(pi/4)*da_i/2),),) # nicht genau die erzweugte Kante auswählen sonst gibt es Probleme Kontaktflaeche6 = f_wi.findAt(((sin(3*pi/4)*da_i/2,l_absatz+0.2,cos(3*pi/4)*da_i/2),),) Kontaktflaeche7 = f_wi.findAt(((sin(5*pi/4)*da_i/2,l_absatz+0.2,cos(5*pi/4)*da_i/2),),) Kontaktflaeche8 = f_wi.findAt(((sin(7*pi/4)*da_i/2,l_absatz+0.2,cos(7*pi/4)*da_i/2),),) Kontaktflaeche9 = f_wi.findAt(((sin(pi/4)*da_i/2,l_absatz+radius[i]+0.1,cos(pi/4)*da_i/2),),) # nicht genau die erzweugte Kante auswählen sonst gibt es Probleme Kontaktflaeche10 = f_wi.findAt(((sin(3*pi/4)*da_i/2,l_absatz+radius[i]+0.1,cos(3*pi/4)*da_i/2),),) Kontaktflaeche11 = f_wi.findAt(((sin(5*pi/4)*da_i/2,l_absatz+radius[i]+0.1,cos(5*pi/4)*da_i/2),),) Kontaktflaeche12 = f_wi.findAt(((sin(7*pi/4)*da_i/2,l_absatz+radius[i]+0.1,cos(7*pi/4)*da_i/2),),) Kontaktflaeche13 = f_wi.findAt(((sin(pi/4)*da_i/2,(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))-0.05,cos(pi/4)*da_i/2),),) # nicht genau die erzweugte Kante auswählen sonst gibt es Probleme Kontaktflaeche14 = f_wi.findAt(((sin(3*pi/4)*da_i/2,(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))-0.05,cos(3*pi/4)*da_i/2),),) Kontaktflaeche15 = f_wi.findAt(((sin(5*pi/4)*da_i/2,(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))-0.05,cos(5*pi/4)*da_i/2),),) Kontaktflaeche16 = f_wi.findAt(((sin(7*pi/4)*da_i/2,(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))-0.05,cos(7*pi/4)*da_i/2),),) Kontaktflaeche17 = f_wi.findAt(((sin(pi/4)*da_i/2,l_absatz+f[i]-0.1,cos(pi/4)*da_i/2),),) # nicht genau die erzweugte Kante auswählen sonst gibt es Probleme Kontaktflaeche18 = f_wi.findAt(((sin(3*pi/4)*da_i/2,l_absatz+f[i]-0.1,cos(3*pi/4)*da_i/2),),) Kontaktflaeche19 = f_wi.findAt(((sin(5*pi/4)*da_i/2,l_absatz+f[i]-0.1,cos(5*pi/4)*da_i/2),),) Kontaktflaeche20 = f_wi.findAt(((sin(7*pi/4)*da_i/2,l_absatz+f[i]-0.1,cos(7*pi/4)*da_i/2),),) Kontaktflaeche21 = f_wi.findAt(((sin(pi/4)*da_i/2,l_absatz+f[i]+2,cos(pi/4)*da_i/2),),) Kontaktflaeche22 = f_wi.findAt(((sin(3*pi/4)*da_i/2,l_absatz+f[i]+2,cos(3*pi/4)*da_i/2),),) Kontaktflaeche23 = f_wi.findAt(((sin(5*pi/4)*da_i/2,l_absatz+f[i]+2,cos(5*pi/4)*da_i/2),),) Kontaktflaeche24 = f_wi.findAt(((sin(7*pi/4)*da_i/2,l_absatz+f[i]+2,cos(7*pi/4)*da_i/2),),) Kontaktflaeche25 = f_wi.findAt(((sin(pi/4)*da_i/2,lw-2,cos(pi/4)*da_i/2),),) # nicht genau die erzweugte Kante auswählen sonst gibt es Probleme Kontaktflaeche26 = f_wi.findAt(((sin(3*pi/4)*da_i/2,lw-2,cos(3*pi/4)*da_i/2),),) Kontaktflaeche27 = f_wi.findAt(((sin(5*pi/4)*da_i/2,lw-2,cos(5*pi/4)*da_i/2),),) Kontaktflaeche28 = f_wi.findAt(((sin(7*pi/4)*da_i/2,lw-2,cos(7*pi/4)*da_i/2),),) Kontaktflaeche29 = f_wi.findAt(((sin(pi/4)*da_i/2,l_absatz-3,cos(pi/4)*da_i/2),),) # nicht genau die erzweugte Kante auswählen sonst gibt es Probleme Kontaktflaeche30 = f_wi.findAt(((sin(3*pi/4)*da_i/2,l_absatz-3,cos(3*pi/4)*da_i/2),),) Kontaktflaeche31 = f_wi.findAt(((sin(5*pi/4)*da_i/2,l_absatz-3,cos(5*pi/4)*da_i/2),),) Kontaktflaeche32 = f_wi.findAt(((sin(7*pi/4)*da_i/2,l_absatz-3,cos(7*pi/4)*da_i/2),),) Kontaktflaeche33 = f_wi.findAt(((sin(pi/4)*da_i/2,l_absatz+f[i]+6,cos(pi/4)*da_i/2),),) Kontaktflaeche34 = f_wi.findAt(((sin(3*pi/4)*da_i/2,l_absatz+f[i]+6,cos(3*pi/4)*da_i/2),),) Kontaktflaeche35 = f_wi.findAt(((sin(5*pi/4)*da_i/2,l_absatz+f[i]+6,cos(5*pi/4)*da_i/2),),) Kontaktflaeche36 = f_wi.findAt(((sin(7*pi/4)*da_i/2,l_absatz+f[i]+6,cos(7*pi/4)*da_i/2),),) Kontaktflaeche_Welle_innen = [Kontaktflaeche1, Kontaktflaeche2, Kontaktflaeche3, Kontaktflaeche4, Kontaktflaeche5, Kontaktflaeche6, Kontaktflaeche7, Kontaktflaeche8, Kontaktflaeche9, Kontaktflaeche10, Kontaktflaeche11, Kontaktflaeche12, Kontaktflaeche13, Kontaktflaeche14, Kontaktflaeche15, Kontaktflaeche15, Kontaktflaeche16, Kontaktflaeche17, Kontaktflaeche18, Kontaktflaeche19, Kontaktflaeche20, Kontaktflaeche21, Kontaktflaeche22, Kontaktflaeche23, Kontaktflaeche24, Kontaktflaeche25, Kontaktflaeche26, Kontaktflaeche27, Kontaktflaeche28, Kontaktflaeche29, Kontaktflaeche30, Kontaktflaeche31, Kontaktflaeche32, Kontaktflaeche33, Kontaktflaeche34, Kontaktflaeche35, Kontaktflaeche36] welle_innen.Surface(side1Faces=Kontaktflaeche_Welle_innen, name='Kontaktflaeche') # Erstellung von Sets # Welle # tangential welle = mdb.models[Modell_Name[i]].parts[Wellen_Name[i]] session.viewports['Viewport: 1'].setValues(displayedObject=welle) e_w = welle.edges if list_di[i] > 1: Tangentialedge_Welle_Absatz= e_w.findAt(((sin (pi/4)*(list_da[i]/2+h_absatz[i]),0,cos (pi/4)*(list_da[i]/2+h_absatz[i])),), ((sin (3*pi/4)*(list_da[i]/2+h_absatz[i]),0,cos (3*pi/4)*(list_da[i]/2+h_absatz[i])),), ((sin (5*pi/4)*(list_da[i]/2+h_absatz[i]),0,cos (5*pi/4)*(list_da[i]/2+h_absatz[i])),), ((sin (7*pi/4)*(list_da[i]/2+h_absatz[i]),0,cos (7*pi/4)*(list_da[i]/2+h_absatz[i])),), ((sin (pi/4)*(list_da[i]/2+radius[i]),0,cos (pi/4)*(list_da[i]/2+radius[i])),), ((sin (3*pi/4)*(list_da[i]/2+radius[i]),0,cos (3*pi/4)*(list_da[i]/2+radius[i])),), ((sin (5*pi/4)*(list_da[i]/2+radius[i]),0,cos (5*pi/4)*(list_da[i]/2+radius[i])),), ((sin (7*pi/4)*(list_da[i]/2+radius[i]),0,cos (7*pi/4)*(list_da[i]/2+radius[i])),), ((sin (pi/4)*(list_da[i]/2-1),0,cos (pi/4)*(list_da[i]/2-1)),), ((sin (3*pi/4)*(list_da[i]/2-1),0,cos (3*pi/4)*(list_da[i]/2-1)),), ((sin (5*pi/4)*(list_da[i]/2-1),0,cos (5*pi/4)*(list_da[i]/2-1)),), ((sin (7*pi/4)*(list_da[i]/2-1),0,cos (7*pi/4)*(list_da[i]/2-1)),), ((sin (pi/4)*(list_di[i]/2),0,cos (pi/4)*(list_di[i]/2)),), ((sin (3*pi/4)*(list_di[i]/2),0,cos (3*pi/4)*(list_di[i]/2)),), ((sin (5*pi/4)*(list_di[i]/2),0,cos (5*pi/4)*(list_di[i]/2)),), ((sin (7*pi/4)*(list_di[i]/2),0,cos (7*pi/4)*(list_di[i]/2)),),) else: Tangentialedge_Welle_Absatz= e_w.findAt(((sin (pi/4)*(list_da[i]/2+h_absatz[i]),0,cos (pi/4)*(list_da[i]/2+h_absatz[i])),), ((sin (3*pi/4)*(list_da[i]/2+h_absatz[i]),0,cos (3*pi/4)*(list_da[i]/2+h_absatz[i])),), ((sin (5*pi/4)*(list_da[i]/2+h_absatz[i]),0,cos (5*pi/4)*(list_da[i]/2+h_absatz[i])),), ((sin (7*pi/4)*(list_da[i]/2+h_absatz[i]),0,cos (7*pi/4)*(list_da[i]/2+h_absatz[i])),), ((sin (pi/4)*(list_da[i]/2+radius[i]),0,cos (pi/4)*(list_da[i]/2+radius[i])),), ((sin (3*pi/4)*(list_da[i]/2+radius[i]),0,cos (3*pi/4)*(list_da[i]/2+radius[i])),), ((sin (5*pi/4)*(list_da[i]/2+radius[i]),0,cos (5*pi/4)*(list_da[i]/2+radius[i])),), ((sin (7*pi/4)*(list_da[i]/2+radius[i]),0,cos (7*pi/4)*(list_da[i]/2+radius[i])),), ((sin (pi/4)*(list_da[i]/2-1),0,cos (pi/4)*(list_da[i]/2-1)),), ((sin (3*pi/4)*(list_da[i]/2-1),0,cos (3*pi/4)*(list_da[i]/2-1)),), ((sin (5*pi/4)*(list_da[i]/2-1),0,cos (5*pi/4)*(list_da[i]/2-1)),), ((sin (7*pi/4)*(list_da[i]/2-1),0,cos (7*pi/4)*(list_da[i]/2-1)),),) welle.Set(edges=Tangentialedge_Welle_Absatz, name='Tangential_01') if list_di[i] > 1: Tangentialedge_Welle= e_w.findAt(((sin (pi/4)*(list_da[i]/2+h_absatz[i]),l_absatz-2,cos (pi/4)*(list_da[i]/2+h_absatz[i])),), # Kante vor Absatz ((sin (3*pi/4)*(list_da[i]/2+h_absatz[i]),l_absatz-2,cos (3*pi/4)*(list_da[i]/2+h_absatz[i])),), ((sin (5*pi/4)*(list_da[i]/2+h_absatz[i]),l_absatz-2,cos (5*pi/4)*(list_da[i]/2+h_absatz[i])),), ((sin (7*pi/4)*(list_da[i]/2+h_absatz[i]),l_absatz-2,cos (7*pi/4)*(list_da[i]/2+h_absatz[i])),), ((sin (pi/4)*(list_da[i]/2+radius[i]),l_absatz-2,cos (pi/4)*(list_da[i]/2+radius[i])),), ((sin (3*pi/4)*(list_da[i]/2+radius[i]),l_absatz-2,cos (3*pi/4)*(list_da[i]/2+radius[i])),), ((sin (5*pi/4)*(list_da[i]/2+radius[i]),l_absatz-2,cos (5*pi/4)*(list_da[i]/2+radius[i])),), ((sin (7*pi/4)*(list_da[i]/2+radius[i]),l_absatz-2,cos (7*pi/4)*(list_da[i]/2+radius[i])),), ((sin (pi/4)*(list_da[i]/2-1),l_absatz-2,cos (pi/4)*(list_da[i]/2-1)),), ((sin (3*pi/4)*(list_da[i]/2-1),l_absatz-2,cos (3*pi/4)*(list_da[i]/2-1)),), ((sin (5*pi/4)*(list_da[i]/2-1),l_absatz-2,cos (5*pi/4)*(list_da[i]/2-1)),), ((sin (7*pi/4)*(list_da[i]/2-1),l_absatz-2,cos (7*pi/4)*(list_da[i]/2-1)),), ((sin (pi/4)*(list_di[i]/2),l_absatz-2,cos (pi/4)*(list_di[i]/2)),), ((sin (3*pi/4)*(list_di[i]/2),l_absatz-2,cos (3*pi/4)*(list_di[i]/2)),), ((sin (5*pi/4)*(list_di[i]/2),l_absatz-2,cos (5*pi/4)*(list_di[i]/2)),), ((sin (7*pi/4)*(list_di[i]/2),l_absatz-2,cos (7*pi/4)*(list_di[i]/2)),), ((sin (pi/4)*(list_da[i]/2+h_absatz[i]),l_absatz,cos (pi/4)*(list_da[i]/2+h_absatz[i])),), # Absatzkante ((sin (3*pi/4)*(list_da[i]/2+h_absatz[i]),l_absatz,cos (3*pi/4)*(list_da[i]/2+h_absatz[i])),), ((sin (5*pi/4)*(list_da[i]/2+h_absatz[i]),l_absatz,cos (5*pi/4)*(list_da[i]/2+h_absatz[i])),), ((sin (7*pi/4)*(list_da[i]/2+h_absatz[i]),l_absatz,cos (7*pi/4)*(list_da[i]/2+h_absatz[i])),), ((sin (pi/4)*(list_da[i]/2+radius[i]),l_absatz,cos (pi/4)*(list_da[i]/2+radius[i])),), ((sin (3*pi/4)*(list_da[i]/2+radius[i]),l_absatz,cos (3*pi/4)*(list_da[i]/2+radius[i])),), ((sin (5*pi/4)*(list_da[i]/2+radius[i]),l_absatz,cos (5*pi/4)*(list_da[i]/2+radius[i])),), ((sin (7*pi/4)*(list_da[i]/2+radius[i]),l_absatz,cos (7*pi/4)*(list_da[i]/2+radius[i])),), ((sin (pi/4)*(list_di[i]/2),l_absatz+radius[i],cos (pi/4)*(list_di[i]/2)),), # nach Erstem Radius Kante ((sin (3*pi/4)*(list_di[i]/2),l_absatz+radius[i],cos (3*pi/4)*(list_di[i]/2)),), ((sin (5*pi/4)*(list_di[i]/2),l_absatz+radius[i],cos (5*pi/4)*(list_di[i]/2)),), ((sin (7*pi/4)*(list_di[i]/2),l_absatz+radius[i],cos (7*pi/4)*(list_di[i]/2)),), ((sin (pi/4)*(list_da[i]/2-t1[i]),l_absatz+radius[i],cos (pi/4)*(list_da[i]/2-t1[i])),), ((sin (3*pi/4)*(list_da[i]/2-t1[i]),l_absatz+radius[i],cos (3*pi/4)*(list_da[i]/2-t1[i])),), ((sin (5*pi/4)*(list_da[i]/2-t1[i]),l_absatz+radius[i],cos (5*pi/4)*(list_da[i]/2-t1[i])),), ((sin (7*pi/4)*(list_da[i]/2-t1[i]),l_absatz+radius[i],cos (7*pi/4)*(list_da[i]/2-t1[i])),), ((sin (pi/4)*(list_da[i]/2-1),l_absatz+radius[i],cos (pi/4)*(list_da[i]/2-1)),), ((sin (3*pi/4)*(list_da[i]/2-1),l_absatz+radius[i],cos (3*pi/4)*(list_da[i]/2-1)),), ((sin (5*pi/4)*(list_da[i]/2-1),l_absatz+radius[i],cos (5*pi/4)*(list_da[i]/2-1)),), ((sin (7*pi/4)*(list_da[i]/2-1),l_absatz+radius[i],cos (7*pi/4)*(list_da[i]/2-1)),), ((sin (pi/4)*(list_di[i]/2),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos (pi/4)*(list_di[i]/2)),), # nach horizontaler Linie Kante ((sin (3*pi/4)*(list_di[i]/2),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos (3*pi/4)*(list_di[i]/2)),), ((sin (5*pi/4)*(list_di[i]/2),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos (5*pi/4)*(list_di[i]/2)),), ((sin (7*pi/4)*(list_di[i]/2),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos (7*pi/4)*(list_di[i]/2)),), ((sin (pi/4)*(list_da[i]/2-t1[i]),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos (pi/4)*(list_da[i]/2-t1[i])),), ((sin (3*pi/4)*(list_da[i]/2-t1[i]),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos (3*pi/4)*(list_da[i]/2-t1[i])),), ((sin (5*pi/4)*(list_da[i]/2-t1[i]),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos (5*pi/4)*(list_da[i]/2-t1[i])),), ((sin (7*pi/4)*(list_da[i]/2-t1[i]),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos (7*pi/4)*(list_da[i]/2-t1[i])),), ((sin (pi/4)*(list_da[i]/2-1),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos (pi/4)*(list_da[i]/2-1)),), ((sin (3*pi/4)*(list_da[i]/2-1),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos (3*pi/4)*(list_da[i]/2-1)),), ((sin (5*pi/4)*(list_da[i]/2-1),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos (5*pi/4)*(list_da[i]/2-1)),), ((sin (7*pi/4)*(list_da[i]/2-1),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos (7*pi/4)*(list_da[i]/2-1)),), ((sin (pi/4)*(list_di[i]/2),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos (pi/4)*(list_di[i]/2)),), # nach Zweitem Radius Kante ((sin (3*pi/4)*(list_di[i]/2),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos (3*pi/4)*(list_di[i]/2)),), ((sin (5*pi/4)*(list_di[i]/2),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos (5*pi/4)*(list_di[i]/2)),), ((sin (7*pi/4)*(list_di[i]/2),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos (7*pi/4)*(list_di[i]/2)),), ((sin (pi/4)*(list_da[i]/2-(radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos (pi/4)*(list_da[i]/2-(radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i])))),), ((sin (3*pi/4)*(list_da[i]/2-(radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos (3*pi/4)*(list_da[i]/2-(radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i])))),), ((sin (5*pi/4)*(list_da[i]/2-(radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos (5*pi/4)*(list_da[i]/2-(radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i])))),), ((sin (7*pi/4)*(list_da[i]/2-(radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos (7*pi/4)*(list_da[i]/2-(radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i])))),), ((sin (pi/4)*(list_da[i]/2-1),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos (pi/4)*(list_da[i]/2-1)),), ((sin (3*pi/4)*(list_da[i]/2-1),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos (3*pi/4)*(list_da[i]/2-1)),), ((sin (5*pi/4)*(list_da[i]/2-1),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos (5*pi/4)*(list_da[i]/2-1)),), ((sin (7*pi/4)*(list_da[i]/2-1),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos (7*pi/4)*(list_da[i]/2-1)),), ((sin (pi/4)*(list_di[i]/2),l_absatz+f[i],cos (pi/4)*(list_di[i]/2)),), # nach Erstem Radius Kante ((sin (3*pi/4)*(list_di[i]/2),l_absatz+f[i],cos (3*pi/4)*(list_di[i]/2)),), ((sin (5*pi/4)*(list_di[i]/2),l_absatz+f[i],cos (5*pi/4)*(list_di[i]/2)),), ((sin (7*pi/4)*(list_di[i]/2),l_absatz+f[i],cos (7*pi/4)*(list_di[i]/2)),), ((sin (pi/4)*(list_da[i]/2),l_absatz+f[i],cos (pi/4)*(list_da[i]/2)),), ((sin (3*pi/4)*(list_da[i]/2),l_absatz+f[i],cos (3*pi/4)*(list_da[i]/2)),), ((sin (5*pi/4)*(list_da[i]/2),l_absatz+f[i],cos (5*pi/4)*(list_da[i]/2)),), ((sin (7*pi/4)*(list_da[i]/2),l_absatz+f[i],cos (7*pi/4)*(list_da[i]/2)),), ((sin (pi/4)*(list_da[i]/2-1),l_absatz+f[i],cos (pi/4)*(list_da[i]/2-1)),), ((sin (3*pi/4)*(list_da[i]/2-1),l_absatz+f[i],cos (3*pi/4)*(list_da[i]/2-1)),), ((sin (5*pi/4)*(list_da[i]/2-1),l_absatz+f[i],cos (5*pi/4)*(list_da[i]/2-1)),), ((sin (7*pi/4)*(list_da[i]/2-1),l_absatz+f[i],cos (7*pi/4)*(list_da[i]/2-1)),), ((sin (pi/4)*(list_di[i]/2),l_absatz+f[i]+list_lf[i],cos (pi/4)*(list_di[i]/2)),), # nach Erstem Radius Kante ((sin (3*pi/4)*(list_di[i]/2),l_absatz+f[i]+list_lf[i],cos (3*pi/4)*(list_di[i]/2)),), ((sin (5*pi/4)*(list_di[i]/2),l_absatz+f[i]+list_lf[i],cos (5*pi/4)*(list_di[i]/2)),), ((sin (7*pi/4)*(list_di[i]/2),l_absatz+f[i]+list_lf[i],cos (7*pi/4)*(list_di[i]/2)),), ((sin (pi/4)*(list_da[i]/2),l_absatz+f[i]+list_lf[i],cos (pi/4)*(list_da[i]/2)),), ((sin (3*pi/4)*(list_da[i]/2),l_absatz+f[i]+list_lf[i],cos (3*pi/4)*(list_da[i]/2)),), ((sin (5*pi/4)*(list_da[i]/2),l_absatz+f[i]+list_lf[i],cos (5*pi/4)*(list_da[i]/2)),), ((sin (7*pi/4)*(list_da[i]/2),l_absatz+f[i]+list_lf[i],cos (7*pi/4)*(list_da[i]/2)),), ((sin (pi/4)*(list_da[i]/2-1),l_absatz+f[i]+list_lf[i],cos (pi/4)*(list_da[i]/2-1)),), ((sin (3*pi/4)*(list_da[i]/2-1),l_absatz+f[i]+list_lf[i],cos (3*pi/4)*(list_da[i]/2-1)),), ((sin (5*pi/4)*(list_da[i]/2-1),l_absatz+f[i]+list_lf[i],cos (5*pi/4)*(list_da[i]/2-1)),), ((sin (7*pi/4)*(list_da[i]/2-1),l_absatz+f[i]+list_lf[i],cos (7*pi/4)*(list_da[i]/2-1)),), ((sin (pi/4)*(list_di[i]/2),l_absatz+f[i]+list_lf[i]/2,cos (pi/4)*(list_di[i]/2)),), # nach Erstem Radius Kante ((sin (3*pi/4)*(list_di[i]/2),l_absatz+f[i]+list_lf[i]/2,cos (3*pi/4)*(list_di[i]/2)),), ((sin (5*pi/4)*(list_di[i]/2),l_absatz+f[i]+list_lf[i]/2,cos (5*pi/4)*(list_di[i]/2)),), ((sin (7*pi/4)*(list_di[i]/2),l_absatz+f[i]+list_lf[i]/2,cos (7*pi/4)*(list_di[i]/2)),), ((sin (pi/4)*(list_da[i]/2),l_absatz+f[i]+list_lf[i]/2,cos (pi/4)*(list_da[i]/2)),), ((sin (3*pi/4)*(list_da[i]/2),l_absatz+f[i]+list_lf[i]/2,cos (3*pi/4)*(list_da[i]/2)),), ((sin (5*pi/4)*(list_da[i]/2),l_absatz+f[i]+list_lf[i]/2,cos (5*pi/4)*(list_da[i]/2)),), ((sin (7*pi/4)*(list_da[i]/2),l_absatz+f[i]+list_lf[i]/2,cos (7*pi/4)*(list_da[i]/2)),), ((sin (pi/4)*(list_da[i]/2-1),l_absatz+f[i]+list_lf[i]/2,cos (pi/4)*(list_da[i]/2-1)),), ((sin (3*pi/4)*(list_da[i]/2-1),l_absatz+f[i]+list_lf[i]/2,cos (3*pi/4)*(list_da[i]/2-1)),), ((sin (5*pi/4)*(list_da[i]/2-1),l_absatz+f[i]+list_lf[i]/2,cos (5*pi/4)*(list_da[i]/2-1)),), ((sin (7*pi/4)*(list_da[i]/2-1),l_absatz+f[i]+list_lf[i]/2,cos (7*pi/4)*(list_da[i]/2-1)),), ((sin (pi/4)*(list_di[i]/2),lw,cos (pi/4)*(list_di[i]/2)),), # nach Erstem Radius Kante ((sin (3*pi/4)*(list_di[i]/2),lw,cos (3*pi/4)*(list_di[i]/2)),), ((sin (5*pi/4)*(list_di[i]/2),lw,cos (5*pi/4)*(list_di[i]/2)),), ((sin (7*pi/4)*(list_di[i]/2),lw,cos (7*pi/4)*(list_di[i]/2)),), ((sin (pi/4)*(list_da[i]/2),lw,cos (pi/4)*(list_da[i]/2)),), ((sin (3*pi/4)*(list_da[i]/2),lw,cos (3*pi/4)*(list_da[i]/2)),), ((sin (5*pi/4)*(list_da[i]/2),lw,cos (5*pi/4)*(list_da[i]/2)),), ((sin (7*pi/4)*(list_da[i]/2),lw,cos (7*pi/4)*(list_da[i]/2)),), ((sin (pi/4)*(list_da[i]/2-1),lw,cos (pi/4)*(list_da[i]/2-1)),), ((sin (3*pi/4)*(list_da[i]/2-1),lw,cos (3*pi/4)*(list_da[i]/2-1)),), ((sin (5*pi/4)*(list_da[i]/2-1),lw,cos (5*pi/4)*(list_da[i]/2-1)),), ((sin (7*pi/4)*(list_da[i]/2-1),lw,cos (7*pi/4)*(list_da[i]/2-1)),), ((sin (pi/4)*(list_da[i]/2+h_absatz[i]),l_absatz-20,cos (pi/4)*(list_da[i]/2+h_absatz[i])),), ((sin (3*pi/4)*(list_da[i]/2+h_absatz[i]),l_absatz-20,cos (3*pi/4)*(list_da[i]/2+h_absatz[i])),), ((sin (5*pi/4)*(list_da[i]/2+h_absatz[i]),l_absatz-20,cos (5*pi/4)*(list_da[i]/2+h_absatz[i])),), ((sin (7*pi/4)*(list_da[i]/2+h_absatz[i]),l_absatz-20,cos (7*pi/4)*(list_da[i]/2+h_absatz[i])),), ((sin (pi/4)*(list_da[i]/2+radius[i]),l_absatz-20,cos (pi/4)*(list_da[i]/2+radius[i])),), ((sin (3*pi/4)*(list_da[i]/2+radius[i]),l_absatz-20,cos (3*pi/4)*(list_da[i]/2+radius[i])),), ((sin (5*pi/4)*(list_da[i]/2+radius[i]),l_absatz-20,cos (5*pi/4)*(list_da[i]/2+radius[i])),), ((sin (7*pi/4)*(list_da[i]/2+radius[i]),l_absatz-20,cos (7*pi/4)*(list_da[i]/2+radius[i])),), ((sin (pi/4)*(list_da[i]/2-1),l_absatz-20,cos (pi/4)*(list_da[i]/2-1)),), ((sin (3*pi/4)*(list_da[i]/2-1),l_absatz-20,cos (3*pi/4)*(list_da[i]/2-1)),), ((sin (5*pi/4)*(list_da[i]/2-1),l_absatz-20,cos (5*pi/4)*(list_da[i]/2-1)),), ((sin (7*pi/4)*(list_da[i]/2-1),l_absatz-20,cos (7*pi/4)*(list_da[i]/2-1)),), ((sin (pi/4)*(list_di[i]/2),l_absatz-20,cos (pi/4)*(list_di[i]/2)),), ((sin (3*pi/4)*(list_di[i]/2),l_absatz-20,cos (3*pi/4)*(list_di[i]/2)),), ((sin (5*pi/4)*(list_di[i]/2),l_absatz-20,cos (5*pi/4)*(list_di[i]/2)),), ((sin (7*pi/4)*(list_di[i]/2),l_absatz-20,cos (7*pi/4)*(list_di[i]/2)),), ((sin (pi/4)*(list_di[i]/2),l_absatz+f[i]+4,cos (pi/4)*(list_di[i]/2)),), # nach Erstem Radius Kante ((sin (3*pi/4)*(list_di[i]/2),l_absatz+f[i]+4,cos (3*pi/4)*(list_di[i]/2)),), ((sin (5*pi/4)*(list_di[i]/2),l_absatz+f[i]+4,cos (5*pi/4)*(list_di[i]/2)),), ((sin (7*pi/4)*(list_di[i]/2),l_absatz+f[i]+4,cos (7*pi/4)*(list_di[i]/2)),), ((sin (pi/4)*(list_da[i]/2),l_absatz+f[i]+4,cos (pi/4)*(list_da[i]/2)),), ((sin (3*pi/4)*(list_da[i]/2),l_absatz+f[i]+4,cos (3*pi/4)*(list_da[i]/2)),), ((sin (5*pi/4)*(list_da[i]/2),l_absatz+f[i]+4,cos (5*pi/4)*(list_da[i]/2)),), ((sin (7*pi/4)*(list_da[i]/2),l_absatz+f[i]+4,cos (7*pi/4)*(list_da[i]/2)),), ((sin (pi/4)*(list_da[i]/2-1),l_absatz+f[i]+4,cos (pi/4)*(list_da[i]/2-1)),), ((sin (3*pi/4)*(list_da[i]/2-1),l_absatz+f[i]+4,cos (3*pi/4)*(list_da[i]/2-1)),), ((sin (5*pi/4)*(list_da[i]/2-1),l_absatz+f[i]+4,cos (5*pi/4)*(list_da[i]/2-1)),), ((sin (7*pi/4)*(list_da[i]/2-1),l_absatz+f[i]+4,cos (7*pi/4)*(list_da[i]/2-1)),),) else: Tangentialedge_Welle= e_w.findAt(((sin (pi/4)*(list_da[i]/2+h_absatz[i]),l_absatz-2,cos (pi/4)*(list_da[i]/2+h_absatz[i])),), # Kante vor Absatz ((sin (3*pi/4)*(list_da[i]/2+h_absatz[i]),l_absatz-2,cos (3*pi/4)*(list_da[i]/2+h_absatz[i])),), ((sin (5*pi/4)*(list_da[i]/2+h_absatz[i]),l_absatz-2,cos (5*pi/4)*(list_da[i]/2+h_absatz[i])),), ((sin (7*pi/4)*(list_da[i]/2+h_absatz[i]),l_absatz-2,cos (7*pi/4)*(list_da[i]/2+h_absatz[i])),), ((sin (pi/4)*(list_da[i]/2+radius[i]),l_absatz-2,cos (pi/4)*(list_da[i]/2+radius[i])),), ((sin (3*pi/4)*(list_da[i]/2+radius[i]),l_absatz-2,cos (3*pi/4)*(list_da[i]/2+radius[i])),), ((sin (5*pi/4)*(list_da[i]/2+radius[i]),l_absatz-2,cos (5*pi/4)*(list_da[i]/2+radius[i])),), ((sin (7*pi/4)*(list_da[i]/2+radius[i]),l_absatz-2,cos (7*pi/4)*(list_da[i]/2+radius[i])),), ((sin (pi/4)*(list_da[i]/2-1),l_absatz-2,cos (pi/4)*(list_da[i]/2-1)),), ((sin (3*pi/4)*(list_da[i]/2-1),l_absatz-2,cos (3*pi/4)*(list_da[i]/2-1)),), ((sin (5*pi/4)*(list_da[i]/2-1),l_absatz-2,cos (5*pi/4)*(list_da[i]/2-1)),), ((sin (7*pi/4)*(list_da[i]/2-1),l_absatz-2,cos (7*pi/4)*(list_da[i]/2-1)),), ((sin (pi/4)*(list_da[i]/2+h_absatz[i]),l_absatz,cos (pi/4)*(list_da[i]/2+h_absatz[i])),), # Absatzkante ((sin (3*pi/4)*(list_da[i]/2+h_absatz[i]),l_absatz,cos (3*pi/4)*(list_da[i]/2+h_absatz[i])),), ((sin (5*pi/4)*(list_da[i]/2+h_absatz[i]),l_absatz,cos (5*pi/4)*(list_da[i]/2+h_absatz[i])),), ((sin (7*pi/4)*(list_da[i]/2+h_absatz[i]),l_absatz,cos (7*pi/4)*(list_da[i]/2+h_absatz[i])),), ((sin (pi/4)*(list_da[i]/2+radius[i]),l_absatz,cos (pi/4)*(list_da[i]/2+radius[i])),), ((sin (3*pi/4)*(list_da[i]/2+radius[i]),l_absatz,cos (3*pi/4)*(list_da[i]/2+radius[i])),), ((sin (5*pi/4)*(list_da[i]/2+radius[i]),l_absatz,cos (5*pi/4)*(list_da[i]/2+radius[i])),), ((sin (7*pi/4)*(list_da[i]/2+radius[i]),l_absatz,cos (7*pi/4)*(list_da[i]/2+radius[i])),), ((sin (pi/4)*(list_da[i]/2-t1[i]),l_absatz+radius[i],cos (pi/4)*(list_da[i]/2-t1[i])),), ((sin (3*pi/4)*(list_da[i]/2-t1[i]),l_absatz+radius[i],cos (3*pi/4)*(list_da[i]/2-t1[i])),), ((sin (5*pi/4)*(list_da[i]/2-t1[i]),l_absatz+radius[i],cos (5*pi/4)*(list_da[i]/2-t1[i])),), ((sin (7*pi/4)*(list_da[i]/2-t1[i]),l_absatz+radius[i],cos (7*pi/4)*(list_da[i]/2-t1[i])),), ((sin (pi/4)*(list_da[i]/2-1),l_absatz+radius[i],cos (pi/4)*(list_da[i]/2-1)),), ((sin (3*pi/4)*(list_da[i]/2-1),l_absatz+radius[i],cos (3*pi/4)*(list_da[i]/2-1)),), ((sin (5*pi/4)*(list_da[i]/2-1),l_absatz+radius[i],cos (5*pi/4)*(list_da[i]/2-1)),), ((sin (7*pi/4)*(list_da[i]/2-1),l_absatz+radius[i],cos (7*pi/4)*(list_da[i]/2-1)),), ((sin (pi/4)*(list_da[i]/2-t1[i]),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos (pi/4)*(list_da[i]/2-t1[i])),), ((sin (3*pi/4)*(list_da[i]/2-t1[i]),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos (3*pi/4)*(list_da[i]/2-t1[i])),), ((sin (5*pi/4)*(list_da[i]/2-t1[i]),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos (5*pi/4)*(list_da[i]/2-t1[i])),), ((sin (7*pi/4)*(list_da[i]/2-t1[i]),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos (7*pi/4)*(list_da[i]/2-t1[i])),), ((sin (pi/4)*(list_da[i]/2-1),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos (pi/4)*(list_da[i]/2-1)),), ((sin (3*pi/4)*(list_da[i]/2-1),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos (3*pi/4)*(list_da[i]/2-1)),), ((sin (5*pi/4)*(list_da[i]/2-1),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos (5*pi/4)*(list_da[i]/2-1)),), ((sin (7*pi/4)*(list_da[i]/2-1),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos (7*pi/4)*(list_da[i]/2-1)),), ((sin (pi/4)*(list_da[i]/2-(radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos (pi/4)*(list_da[i]/2-(radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i])))),), ((sin (3*pi/4)*(list_da[i]/2-(radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos (3*pi/4)*(list_da[i]/2-(radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i])))),), ((sin (5*pi/4)*(list_da[i]/2-(radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos (5*pi/4)*(list_da[i]/2-(radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i])))),), ((sin (7*pi/4)*(list_da[i]/2-(radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos (7*pi/4)*(list_da[i]/2-(radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i])))),), ((sin (pi/4)*(list_da[i]/2-1),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos (pi/4)*(list_da[i]/2-1)),), ((sin (3*pi/4)*(list_da[i]/2-1),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos (3*pi/4)*(list_da[i]/2-1)),), ((sin (5*pi/4)*(list_da[i]/2-1),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos (5*pi/4)*(list_da[i]/2-1)),), ((sin (7*pi/4)*(list_da[i]/2-1),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos (7*pi/4)*(list_da[i]/2-1)),), ((sin (pi/4)*(list_da[i]/2),l_absatz+f[i],cos (pi/4)*(list_da[i]/2)),), ((sin (3*pi/4)*(list_da[i]/2),l_absatz+f[i],cos (3*pi/4)*(list_da[i]/2)),), ((sin (5*pi/4)*(list_da[i]/2),l_absatz+f[i],cos (5*pi/4)*(list_da[i]/2)),), ((sin (7*pi/4)*(list_da[i]/2),l_absatz+f[i],cos (7*pi/4)*(list_da[i]/2)),), ((sin (pi/4)*(list_da[i]/2-1),l_absatz+f[i],cos (pi/4)*(list_da[i]/2-1)),), ((sin (3*pi/4)*(list_da[i]/2-1),l_absatz+f[i],cos (3*pi/4)*(list_da[i]/2-1)),), ((sin (5*pi/4)*(list_da[i]/2-1),l_absatz+f[i],cos (5*pi/4)*(list_da[i]/2-1)),), ((sin (7*pi/4)*(list_da[i]/2-1),l_absatz+f[i],cos (7*pi/4)*(list_da[i]/2-1)),), ((sin (pi/4)*(list_da[i]/2),l_absatz+f[i]+list_lf[i],cos (pi/4)*(list_da[i]/2)),), ((sin (3*pi/4)*(list_da[i]/2),l_absatz+f[i]+list_lf[i],cos (3*pi/4)*(list_da[i]/2)),), ((sin (5*pi/4)*(list_da[i]/2),l_absatz+f[i]+list_lf[i],cos (5*pi/4)*(list_da[i]/2)),), ((sin (7*pi/4)*(list_da[i]/2),l_absatz+f[i]+list_lf[i],cos (7*pi/4)*(list_da[i]/2)),), ((sin (pi/4)*(list_da[i]/2-1),l_absatz+f[i]+list_lf[i],cos (pi/4)*(list_da[i]/2-1)),), ((sin (3*pi/4)*(list_da[i]/2-1),l_absatz+f[i]+list_lf[i],cos (3*pi/4)*(list_da[i]/2-1)),), ((sin (5*pi/4)*(list_da[i]/2-1),l_absatz+f[i]+list_lf[i],cos (5*pi/4)*(list_da[i]/2-1)),), ((sin (7*pi/4)*(list_da[i]/2-1),l_absatz+f[i]+list_lf[i],cos (7*pi/4)*(list_da[i]/2-1)),), ((sin (pi/4)*(list_da[i]/2),l_absatz+f[i]+list_lf[i]/2,cos (pi/4)*(list_da[i]/2)),), ((sin (3*pi/4)*(list_da[i]/2),l_absatz+f[i]+list_lf[i]/2,cos (3*pi/4)*(list_da[i]/2)),), ((sin (5*pi/4)*(list_da[i]/2),l_absatz+f[i]+list_lf[i]/2,cos (5*pi/4)*(list_da[i]/2)),), ((sin (7*pi/4)*(list_da[i]/2),l_absatz+f[i]+list_lf[i]/2,cos (7*pi/4)*(list_da[i]/2)),), ((sin (pi/4)*(list_da[i]/2-1),l_absatz+f[i]+list_lf[i]/2,cos (pi/4)*(list_da[i]/2-1)),), ((sin (3*pi/4)*(list_da[i]/2-1),l_absatz+f[i]+list_lf[i]/2,cos (3*pi/4)*(list_da[i]/2-1)),), ((sin (5*pi/4)*(list_da[i]/2-1),l_absatz+f[i]+list_lf[i]/2,cos (5*pi/4)*(list_da[i]/2-1)),), ((sin (7*pi/4)*(list_da[i]/2-1),l_absatz+f[i]+list_lf[i]/2,cos (7*pi/4)*(list_da[i]/2-1)),), ((sin (pi/4)*(list_da[i]/2),lw,cos (pi/4)*(list_da[i]/2)),), ((sin (3*pi/4)*(list_da[i]/2),lw,cos (3*pi/4)*(list_da[i]/2)),), ((sin (5*pi/4)*(list_da[i]/2),lw,cos (5*pi/4)*(list_da[i]/2)),), ((sin (7*pi/4)*(list_da[i]/2),lw,cos (7*pi/4)*(list_da[i]/2)),), ((sin (pi/4)*(list_da[i]/2-1),lw,cos (pi/4)*(list_da[i]/2-1)),), ((sin (3*pi/4)*(list_da[i]/2-1),lw,cos (3*pi/4)*(list_da[i]/2-1)),), ((sin (5*pi/4)*(list_da[i]/2-1),lw,cos (5*pi/4)*(list_da[i]/2-1)),), ((sin (7*pi/4)*(list_da[i]/2-1),lw,cos (7*pi/4)*(list_da[i]/2-1)),), ((sin (pi/4)*(list_da[i]/2+h_absatz[i]),l_absatz-20,cos (pi/4)*(list_da[i]/2+h_absatz[i])),), ((sin (3*pi/4)*(list_da[i]/2+h_absatz[i]),l_absatz-20,cos (3*pi/4)*(list_da[i]/2+h_absatz[i])),), ((sin (5*pi/4)*(list_da[i]/2+h_absatz[i]),l_absatz-20,cos (5*pi/4)*(list_da[i]/2+h_absatz[i])),), ((sin (7*pi/4)*(list_da[i]/2+h_absatz[i]),l_absatz-20,cos (7*pi/4)*(list_da[i]/2+h_absatz[i])),), ((sin (pi/4)*(list_da[i]/2+radius[i]),l_absatz-20,cos (pi/4)*(list_da[i]/2+radius[i])),), ((sin (3*pi/4)*(list_da[i]/2+radius[i]),l_absatz-20,cos (3*pi/4)*(list_da[i]/2+radius[i])),), ((sin (5*pi/4)*(list_da[i]/2+radius[i]),l_absatz-20,cos (5*pi/4)*(list_da[i]/2+radius[i])),), ((sin (7*pi/4)*(list_da[i]/2+radius[i]),l_absatz-20,cos (7*pi/4)*(list_da[i]/2+radius[i])),), ((sin (pi/4)*(list_da[i]/2-1),l_absatz-20,cos (pi/4)*(list_da[i]/2-1)),), ((sin (3*pi/4)*(list_da[i]/2-1),l_absatz-20,cos (3*pi/4)*(list_da[i]/2-1)),), ((sin (5*pi/4)*(list_da[i]/2-1),l_absatz-20,cos (5*pi/4)*(list_da[i]/2-1)),), ((sin (7*pi/4)*(list_da[i]/2-1),l_absatz-20,cos (7*pi/4)*(list_da[i]/2-1)),), ((sin (pi/4)*(list_da[i]/2),l_absatz+f[i]+4,cos (pi/4)*(list_da[i]/2)),), ((sin (3*pi/4)*(list_da[i]/2),l_absatz+f[i]+4,cos (3*pi/4)*(list_da[i]/2)),), ((sin (5*pi/4)*(list_da[i]/2),l_absatz+f[i]+4,cos (5*pi/4)*(list_da[i]/2)),), ((sin (7*pi/4)*(list_da[i]/2),l_absatz+f[i]+4,cos (7*pi/4)*(list_da[i]/2)),), ((sin (pi/4)*(list_da[i]/2-1),l_absatz+f[i]+4,cos (pi/4)*(list_da[i]/2-1)),), ((sin (3*pi/4)*(list_da[i]/2-1),l_absatz+f[i]+4,cos (3*pi/4)*(list_da[i]/2-1)),), ((sin (5*pi/4)*(list_da[i]/2-1),l_absatz+f[i]+4,cos (5*pi/4)*(list_da[i]/2-1)),), ((sin (7*pi/4)*(list_da[i]/2-1),l_absatz+f[i]+4,cos (7*pi/4)*(list_da[i]/2-1)),),) welle.Set(edges=Tangentialedge_Welle, name='Tangential_02') Axialedge_Fr_horiz = e_w.findAt(((sin (pi/2)*(list_da[i]/2-t1[i]), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))-0.01, cos (pi/2)*(list_da[i]/2-t1[i])),), ((sin (pi/2)*(list_da[i]/2-1), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))-0.01, cos (pi/2)*(list_da[i]/2-1)),), ((sin (pi/2)*(list_di[i]/2), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))-0.01, cos (pi/2)*(list_di[i]/2)),), ((sin (pi)*(list_da[i]/2-t1[i]), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))-0.01, cos (pi)*(list_da[i]/2-t1[i])),), ((sin (pi)*(list_da[i]/2-1), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))-0.01, cos (pi)*(list_da[i]/2-1)),), ((sin (pi)*(list_di[i]/2), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))-0.01, cos (pi)*(list_di[i]/2)),), ((sin (3*pi/2)*(list_da[i]/2-t1[i]), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))-0.01, cos (3*pi/2)*(list_da[i]/2-t1[i])),), ((sin (3*pi/2)*(list_da[i]/2-1), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))-0.01, cos (3*pi/2)*(list_da[i]/2-1)),), ((sin (3*pi/2)*(list_di[i]/2), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))-0.01, cos (3*pi/2)*(list_di[i]/2)),), ((sin (2*pi)*(list_da[i]/2-t1[i]), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))-0.01, cos (2*pi)*(list_da[i]/2-t1[i])),), ((sin (2*pi)*(list_da[i]/2-1), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))-0.01, cos (2*pi)*(list_da[i]/2-1)),), ((sin (2*pi)*(list_di[i]/2), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))-0.01, cos (2*pi)*(list_di[i]/2)),),) welle.Set(edges=Axialedge_Fr_horiz, name='Axial_01_FR_Horizontale') if list_da[i] < 120: Axialedge_Fr_2R = e_w.findAt(((sin (pi/2)*(list_da[i]/2-(radius[i]*sin((90-7.5)*pi/180)-(radius[i]-t1[i]))), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))+radius[i]*sin(7.5*pi/180), cos (pi/2)*(list_da[i]/2-(radius[i]*sin((90-7.5)*pi/180)-(radius[i]-t1[i])))),), ((sin (pi/2)*(list_da[i]/2-1), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))+0.01, cos (pi/2)*(list_da[i]/2-1)),), ((sin (pi/2)*(list_di[i]/2), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))+0.01, cos (pi/2)*(list_di[i]/2)),), ((sin (pi)*(list_da[i]/2-(radius[i]*sin((90-7.5)*pi/180)-(radius[i]-t1[i]))), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))+radius[i]*sin(7.5*pi/180), cos (pi)*(list_da[i]/2-(radius[i]*sin((90-7.5)*pi/180)-(radius[i]-t1[i])))),), ((sin (pi)*(list_da[i]/2-1), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))+0.01, cos (pi)*(list_da[i]/2-1)),), ((sin (pi)*(list_di[i]/2), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))+0.01, cos (pi)*(list_di[i]/2)),), ((sin (3*pi/2)*(list_da[i]/2-(radius[i]*sin((90-7.5)*pi/180)-(radius[i]-t1[i]))), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))+radius[i]*sin(7.5*pi/180), cos (3*pi/2)*(list_da[i]/2-(radius[i]*sin((90-7.5)*pi/180)-(radius[i]-t1[i])))),), ((sin (3*pi/2)*(list_da[i]/2-1), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))+0.01, cos (3*pi/2)*(list_da[i]/2-1)),), ((sin (3*pi/2)*(list_di[i]/2), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))+0.01, cos (3*pi/2)*(list_di[i]/2)),), ((sin (2*pi)*(list_da[i]/2-(radius[i]*sin((90-7.5)*pi/180)-(radius[i]-t1[i]))), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))+radius[i]*sin(7.5*pi/180), cos (2*pi)*(list_da[i]/2-(radius[i]*sin((90-7.5)*pi/180)-(radius[i]-t1[i])))),), ((sin (2*pi)*(list_da[i]/2-1), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))+0.01, cos (2*pi)*(list_da[i]/2-1)),), ((sin (2*pi)*(list_di[i]/2), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))+0.01, cos (2*pi)*(list_di[i]/2)),),) else: Axialedge_Fr_2R_1 = e_w.findAt(((sin (pi/2)*(list_da[i]/2-1), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))+0.01, cos (pi/2)*(list_da[i]/2-1)),), ((sin (pi/2)*(list_di[i]/2), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))+0.01, cos (pi/2)*(list_di[i]/2)),), ((sin (pi)*(list_da[i]/2-1), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))+0.01, cos (pi)*(list_da[i]/2-1)),), ((sin (pi)*(list_di[i]/2), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))+0.01, cos (pi)*(list_di[i]/2)),), ((sin (3*pi/2)*(list_da[i]/2-1), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))+0.01, cos (3*pi/2)*(list_da[i]/2-1)),), ((sin (3*pi/2)*(list_di[i]/2), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))+0.01, cos (3*pi/2)*(list_di[i]/2)),), ((sin (2*pi)*(list_da[i]/2-1), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))+0.01, cos (2*pi)*(list_da[i]/2-1)),), ((sin (2*pi)*(list_di[i]/2), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))+0.01, cos (2*pi)*(list_di[i]/2)),),) Axialedge_Fr_2R_2 = e_w.getByBoundingBox(sin (pi/2)*(list_da[i]/2-t1[i]-0.1), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))-0.1, -0.1,sin (pi/2)*(list_da[i]/2+0.1), (l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+0.1, 0.1) Axialedge_Fr_2R_3 = e_w.getByBoundingBox(-0.1, ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))-0.1, cos (pi)*(list_da[i]/2-t1[i]+0.2),0.1, (l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+0.1, cos (pi)*(list_da[i]/2-0.6)) Axialedge_Fr_2R_4 = e_w.getByBoundingBox(sin (3*pi/2)*(list_da[i]/2-t1[i]+0.2), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))-0.3, -0.1,sin (3*pi/2)*(list_da[i]/2-0.6), (l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+0.3, 0.1) Axialedge_Fr_2R_5 = e_w.getByBoundingBox(-0.1, ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))-0.1, cos (2*pi)*(list_da[i]/2-t1[i]-0.1),0.1, (l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+0.1, cos (2*pi)*(list_da[i]/2+0.1)) Axialedge_Fr_2R = [Axialedge_Fr_2R_1, Axialedge_Fr_2R_2, Axialedge_Fr_2R_5, Axialedge_Fr_2R_3, Axialedge_Fr_2R_4] welle.Set(edges=Axialedge_Fr_2R, name='Axial_02_FR_2Radius') Axialedge_Fr_schraeg = e_w.findAt(((sin(pi/2)*(list_da[i]/2-0.1*tan(15*pi/180)),l_absatz+f[i]-0.1,cos (pi/2)*(list_da[i]/2-0.1*tan(15*pi/180))),), ((sin(pi/2)*(list_da[i]/2-1),l_absatz+f[i]-0.1,cos(pi/2)*(list_da[i]/2-1)),), ((sin(pi/2)*(list_di[i]/2),l_absatz+f[i]-0.1,cos(pi/2)*(list_di[i]/2)),), ((sin(pi)*(list_da[i]/2-0.1*tan(15*pi/180)),l_absatz+f[i]-0.1,cos (pi)*(list_da[i]/2-0.1*tan(15*pi/180))),), ((sin(pi)*(list_da[i]/2-1),l_absatz+f[i]-0.1,cos(pi)*(list_da[i]/2-1)),), ((sin(pi)*(list_di[i]/2),l_absatz+f[i]-0.1,cos(pi)*(list_di[i]/2)),), ((sin(3*pi/2)*(list_da[i]/2-0.1*tan(15*pi/180)),l_absatz+f[i]-0.1,cos (3*pi/2)*(list_da[i]/2-0.1*tan(15*pi/180))),), ((sin(3*pi/2)*(list_da[i]/2-1),l_absatz+f[i]-0.1,cos(3*pi/2)*(list_da[i]/2-1)),), ((sin(3*pi/2)*(list_di[i]/2),l_absatz+f[i]-0.1,cos(3*pi/2)*(list_di[i]/2)),), ((sin(2*pi)*(list_da[i]/2-0.1*tan(15*pi/180)),l_absatz+f[i]-0.1,cos (2*pi)*(list_da[i]/2-0.1*tan(15*pi/180))),), ((sin(2*pi)*(list_da[i]/2-1),l_absatz+f[i]-0.1,cos(2*pi)*(list_da[i]/2-1)),), ((sin(2*pi)*(list_di[i]/2),l_absatz+f[i]-0.1,cos(2*pi)*(list_di[i]/2)),),) welle.Set(edges=Axialedge_Fr_schraeg, name='Axial_03_FR_schraeg') Axialedge_Fr_1R = e_w.findAt(((sin (pi/2)*(list_da[i]/2-1),l_absatz,cos (pi/2)*(list_da[i]/2-1)),), ((sin (pi/2)*(list_di[i]/2),l_absatz,cos (pi/2)*(list_di[i]/2)),), ((sin (pi)*(list_da[i]/2-1),l_absatz,cos (pi)*(list_da[i]/2-1)),), ((sin (pi)*(list_di[i]/2),l_absatz,cos (pi)*(list_di[i]/2)),), ((sin (3*pi/2)*(list_da[i]/2-1),l_absatz,cos (3*pi/2)*(list_da[i]/2-1)),), ((sin (3*pi/2)*(list_di[i]/2),l_absatz,cos (3*pi/2)*(list_di[i]/2)),), ((sin (2*pi)*(list_da[i]/2-1),l_absatz,cos (2*pi)*(list_da[i]/2-1)),), ((sin (2*pi)*(list_di[i]/2),l_absatz,cos (2*pi)*(list_di[i]/2)),),) welle.Set(edges=Axialedge_Fr_1R, name='Axial_04_FR_1Radius') Axialedge_Fugenflaeche_u = e_w.findAt(((sin (pi/2)*(list_da[i]/2),l_absatz+f[i]+6,cos (pi/2)*(list_da[i]/2)),), ((sin (pi/2)*(list_da[i]/2-1),l_absatz+f[i]+6,cos (pi/2)*(list_da[i]/2-1)),), ((sin (pi/2)*(list_di[i]/2),l_absatz+f[i]+6,cos (pi/2)*(list_di[i]/2)),), ((sin (pi)*(list_da[i]/2),l_absatz+f[i]+6,cos (pi)*(list_da[i]/2)),), ((sin (pi)*(list_da[i]/2-1),l_absatz+f[i]+6,cos (pi)*(list_da[i]/2-1)),), ((sin (pi)*(list_di[i]/2),l_absatz+f[i]+6,cos (pi)*(list_di[i]/2)),), ((sin (3*pi/2)*(list_da[i]/2),l_absatz+f[i]+6,cos (3*pi/2)*(list_da[i]/2)),), ((sin (3*pi/2)*(list_da[i]/2-1),l_absatz+f[i]+6,cos (3*pi/2)*(list_da[i]/2-1)),), ((sin (3*pi/2)*(list_di[i]/2),l_absatz+f[i]+6,cos (3*pi/2)*(list_di[i]/2)),), ((sin (2*pi)*(list_da[i]/2),l_absatz+f[i]+6,cos (2*pi)*(list_da[i]/2)),), ((sin (2*pi)*(list_da[i]/2-1),l_absatz+f[i]+6,cos (2*pi)*(list_da[i]/2-1)),), ((sin (2*pi)*(list_di[i]/2),l_absatz+f[i]+6,cos (2*pi)*(list_di[i]/2)),),) welle.Set(edges=Axialedge_Fugenflaeche_u, name='Axial_05_Fugenflaeche_unten') Axialedge_Fugenflaeche_o = e_w.findAt(((sin (pi/2)*(list_da[i]/2),l_absatz+f[i]+list_lf[i]-2,cos (pi/2)*(list_da[i]/2)),), ((sin (pi/2)*(list_da[i]/2-1),l_absatz+f[i]+list_lf[i]-2,cos (pi/2)*(list_da[i]/2-1)),), ((sin (pi/2)*(list_di[i]/2),l_absatz+f[i]+list_lf[i]-2,cos (pi/2)*(list_di[i]/2)),), ((sin (pi)*(list_da[i]/2),l_absatz+f[i]+list_lf[i]-2,cos (pi)*(list_da[i]/2)),), ((sin (pi)*(list_da[i]/2-1),l_absatz+f[i]+list_lf[i]-2,cos (pi)*(list_da[i]/2-1)),), ((sin (pi)*(list_di[i]/2),l_absatz+f[i]+list_lf[i]-2,cos (pi)*(list_di[i]/2)),), ((sin (3*pi/2)*(list_da[i]/2),l_absatz+f[i]+list_lf[i]-2,cos (3*pi/2)*(list_da[i]/2)),), ((sin (3*pi/2)*(list_da[i]/2-1),l_absatz+f[i]+list_lf[i]-2,cos (3*pi/2)*(list_da[i]/2-1)),), ((sin (3*pi/2)*(list_di[i]/2),l_absatz+f[i]+list_lf[i]-2,cos (3*pi/2)*(list_di[i]/2)),), ((sin (2*pi)*(list_da[i]/2),l_absatz+f[i]+list_lf[i]-2,cos (2*pi)*(list_da[i]/2)),), ((sin (2*pi)*(list_da[i]/2-1),l_absatz+f[i]+list_lf[i]-2,cos (2*pi)*(list_da[i]/2-1)),), ((sin (2*pi)*(list_di[i]/2),l_absatz+f[i]+list_lf[i]-2,cos (2*pi)*(list_di[i]/2)),),) welle.Set(edges=Axialedge_Fugenflaeche_o, name='Axial_06_Fugenflaeche_oben') Axialedge_Kante_vor_Absatz = e_w.findAt(((sin (pi/2)*(list_da[i]/2+radius[i]),l_absatz-4,cos (pi/2)*(list_da[i]/2+radius[i])),), ((sin (pi/2)*(list_da[i]/2-1),l_absatz-4,cos (pi/2)*(list_da[i]/2-1)),), ((sin (pi/2)*(list_di[i]/2),l_absatz-4,cos (pi/2)*(list_di[i]/2)),), ((sin (pi)*(list_da[i]/2+radius[i]),l_absatz-4,cos (pi)*(list_da[i]/2+radius[i])),), ((sin (pi)*(list_da[i]/2-1),l_absatz-4,cos (pi)*(list_da[i]/2-1)),), ((sin (pi)*(list_di[i]/2),l_absatz-4,cos (pi)*(list_di[i]/2)),), ((sin (3*pi/2)*(list_da[i]/2+radius[i]),l_absatz-4,cos (3*pi/2)*(list_da[i]/2+radius[i])),), ((sin (3*pi/2)*(list_da[i]/2-1),l_absatz-4,cos (3*pi/2)*(list_da[i]/2-1)),), ((sin (3*pi/2)*(list_di[i]/2),l_absatz-4,cos (3*pi/2)*(list_di[i]/2)),), ((sin (2*pi)*(list_da[i]/2+radius[i]),l_absatz-4,cos (2*pi)*(list_da[i]/2+radius[i])),), ((sin (2*pi)*(list_da[i]/2-1),l_absatz-4,cos (2*pi)*(list_da[i]/2-1)),), ((sin (2*pi)*(list_di[i]/2),l_absatz-4,cos (2*pi)*(list_di[i]/2)),), ((sin(pi/2)*(list_da[i]/2+h_absatz[i]),l_absatz-4,cos(pi/2)*(list_da[i]/2+h_absatz[i])),), ((sin(pi)*(list_da[i]/2+h_absatz[i]),l_absatz-4,cos(pi)*(list_da[i]/2+h_absatz[i])),), ((sin(3*pi/2)*(list_da[i]/2+h_absatz[i]),l_absatz-4,cos(3*pi/2)*(list_da[i]/2+h_absatz[i])),), ((sin(2*pi)*(list_da[i]/2+h_absatz[i]),l_absatz-4,cos(2*pi)*(list_da[i]/2+h_absatz[i])),),) welle.Set(edges=Axialedge_Kante_vor_Absatz, name='Axial_07_Kante_vor_Absatz') Axialedge_Auslauf = e_w.findAt(((sin (pi/2)*(list_da[i]/2),lw-2,cos (pi/2)*(list_da[i]/2)),), ((sin (pi/2)*(list_da[i]/2-1),lw-2,cos (pi/2)*(list_da[i]/2-1)),), ((sin (pi/2)*(list_di[i]/2),lw-2,cos (pi/2)*(list_di[i]/2)),), ((sin (pi)*(list_da[i]/2),lw-2,cos (pi)*(list_da[i]/2)),), ((sin (pi)*(list_da[i]/2-1),lw-2,cos (pi)*(list_da[i]/2-1)),), ((sin (pi)*(list_di[i]/2),lw-2,cos (pi)*(list_di[i]/2)),), ((sin (3*pi/2)*(list_da[i]/2),lw-2,cos (3*pi/2)*(list_da[i]/2)),), ((sin (3*pi/2)*(list_da[i]/2-1),lw-2,cos (3*pi/2)*(list_da[i]/2-1)),), ((sin (3*pi/2)*(list_di[i]/2),lw-2,cos (3*pi/2)*(list_di[i]/2)),), ((sin (2*pi)*(list_da[i]/2),lw-2,cos (2*pi)*(list_da[i]/2)),), ((sin (2*pi)*(list_da[i]/2-1),lw-2,cos (2*pi)*(list_da[i]/2-1)),), ((sin (2*pi)*(list_di[i]/2),lw-2,cos (2*pi)*(list_di[i]/2)),),) welle.Set(edges=Axialedge_Auslauf, name='Axial_08_Auslauf') Axialedge_Absatz = e_w.findAt(((sin (pi/2)*(list_da[i]/2+radius[i]),4,cos (pi/2)*(list_da[i]/2+radius[i])),), ((sin (pi/2)*(list_da[i]/2-1),4,cos (pi/2)*(list_da[i]/2-1)),), ((sin (pi/2)*(list_di[i]/2),4,cos (pi/2)*(list_di[i]/2)),), ((sin (pi/2)*(list_da[i]/2+h_absatz[i]),4,cos (pi/2)*(list_da[i]/2+h_absatz[i])),), ((sin (pi)*(list_da[i]/2+radius[i]),4,cos (pi)*(list_da[i]/2+radius[i])),), ((sin (pi)*(list_da[i]/2-1),4,cos (pi)*(list_da[i]/2-1)),), ((sin (pi)*(list_di[i]/2),4,cos (pi)*(list_di[i]/2)),), ((sin (pi)*(list_da[i]/2+h_absatz[i]),4,cos (pi)*(list_da[i]/2+h_absatz[i])),), ((sin (3*pi/2)*(list_da[i]/2+radius[i]),4,cos (3*pi/2)*(list_da[i]/2+radius[i])),), ((sin (3*pi/2)*(list_da[i]/2-1),4,cos (3*pi/2)*(list_da[i]/2-1)),), ((sin (3*pi/2)*(list_di[i]/2),4,cos (3*pi/2)*(list_di[i]/2)),), ((sin (3*pi/2)*(list_da[i]/2+h_absatz[i]),4,cos (3*pi/2)*(list_da[i]/2+h_absatz[i])),), ((sin (2*pi)*(list_da[i]/2+radius[i]),4,cos (2*pi)*(list_da[i]/2+radius[i])),), ((sin (2*pi)*(list_da[i]/2-1),4,cos (2*pi)*(list_da[i]/2-1)),), ((sin (2*pi)*(list_di[i]/2),4,cos (2*pi)*(list_di[i]/2)),), ((sin (2*pi)*(list_da[i]/2+h_absatz[i]),4,cos (2*pi)*(list_da[i]/2+h_absatz[i])),),) welle.Set(edges=Axialedge_Absatz, name='Axial_09_Absatz') Axialedge_Absatzkante = e_w.findAt(((sin(pi/2)*(list_da[i]/2+h_absatz[i]),l_absatz-0.1,cos(pi/2)*(list_da[i]/2+h_absatz[i])),), ((sin(pi)*(list_da[i]/2+h_absatz[i]),l_absatz-0.1,cos(pi)*(list_da[i]/2+h_absatz[i])),), ((sin(3*pi/2)*(list_da[i]/2+h_absatz[i]),l_absatz-0.1,cos(3*pi/2)*(list_da[i]/2+h_absatz[i])),), ((sin(2*pi)*(list_da[i]/2+h_absatz[i]),l_absatz-0.1,cos(2*pi)*(list_da[i]/2+h_absatz[i])),), ((sin(pi/2)*(list_da[i]/2+radius[i]),l_absatz-0.1,cos(pi/2)*(list_da[i]/2+radius[i])),), ((sin(pi)*(list_da[i]/2+radius[i]),l_absatz-0.1,cos(pi)*(list_da[i]/2+radius[i])),), ((sin(3*pi/2)*(list_da[i]/2+radius[i]),l_absatz-0.1,cos(3*pi/2)*(list_da[i]/2+radius[i])),), ((sin(2*pi)*(list_da[i]/2+radius[i]),l_absatz-0.1,cos(2*pi)*(list_da[i]/2+radius[i])),),) welle.Set(edges=Axialedge_Absatzkante, name='Axial_10_Absatzkante') Axialedge_Anfang_Fuge = e_w.findAt(((sin (pi/2)*(list_da[i]/2),l_absatz+f[i]+3,cos (pi/2)*(list_da[i]/2)),), ((sin (pi/2)*(list_da[i]/2-1),l_absatz+f[i]+3,cos (pi/2)*(list_da[i]/2-1)),), ((sin (pi/2)*(list_di[i]/2),l_absatz+f[i]+3,cos (pi/2)*(list_di[i]/2)),), ((sin (pi)*(list_da[i]/2),l_absatz+f[i]+3,cos (pi)*(list_da[i]/2)),), ((sin (pi)*(list_da[i]/2-1),l_absatz+f[i]+3,cos (pi)*(list_da[i]/2-1)),), ((sin (pi)*(list_di[i]/2),l_absatz+f[i]+3,cos (pi)*(list_di[i]/2)),), ((sin (3*pi/2)*(list_da[i]/2),l_absatz+f[i]+3,cos (3*pi/2)*(list_da[i]/2)),), ((sin (3*pi/2)*(list_da[i]/2-1),l_absatz+f[i]+3,cos (3*pi/2)*(list_da[i]/2-1)),), ((sin (3*pi/2)*(list_di[i]/2),l_absatz+f[i]+3,cos (3*pi/2)*(list_di[i]/2)),), ((sin (2*pi)*(list_da[i]/2),l_absatz+f[i]+3,cos (2*pi)*(list_da[i]/2)),), ((sin (2*pi)*(list_da[i]/2-1),l_absatz+f[i]+3,cos (2*pi)*(list_da[i]/2-1)),), ((sin (2*pi)*(list_di[i]/2),l_absatz+f[i]+3,cos (2*pi)*(list_di[i]/2)),),) welle.Set(edges=Axialedge_Anfang_Fuge, name='Axial_11_Anfang_Fugenflaeche') # radial Radialedge_FR_Fugenflaeche_a = e_w.findAt(((sin(pi/2)*(list_da[i]/2-1+0.1),l_absatz+radius[i],cos(pi/2)*(list_da[i]/2-1+0.1)),), ((sin(pi/2)*(list_da[i]/2-1+0.1),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos(pi/2)*(list_da[i]/2-1+0.1)),), ((sin(pi/2)*(list_da[i]/2-1+0.1),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos(pi/2)*(list_da[i]/2-1+0.1)),), ((sin(pi/2)*(list_da[i]/2-1+0.1),l_absatz+f[i],cos(pi/2)*(list_da[i]/2-1+0.1)),), ((sin(pi/2)*(list_da[i]/2-1+0.1),l_absatz+f[i]+4,cos(pi/2)*(list_da[i]/2-1+0.1)),), ((sin(pi/2)*(list_da[i]/2-1+0.1),(l_absatz+f[i]+list_lf[i]),cos(pi/2)*(list_da[i]/2-1+0.1)),), ((sin(pi/2)*(list_da[i]/2-1+0.1),(l_absatz+f[i]+list_lf[i]/2),cos(pi/2)*(list_da[i]/2-1+0.1)),), ((sin(pi)*(list_da[i]/2-1+0.1),l_absatz+radius[i],cos(pi)*(list_da[i]/2-1+0.1)),), ((sin(pi)*(list_da[i]/2-1+0.1),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos(pi)*(list_da[i]/2-1+0.1)),), ((sin(pi)*(list_da[i]/2-1+0.1),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos(pi)*(list_da[i]/2-1+0.1)),), ((sin(pi)*(list_da[i]/2-1+0.1),l_absatz+f[i],cos(pi)*(list_da[i]/2-1+0.1)),), ((sin(pi)*(list_da[i]/2-1+0.1),l_absatz+f[i]+4,cos(pi)*(list_da[i]/2-1+0.1)),), ((sin(pi)*(list_da[i]/2-1+0.1),(l_absatz+f[i]+list_lf[i]),cos(pi)*(list_da[i]/2-1+0.1)),), ((sin(pi)*(list_da[i]/2-1+0.1),(l_absatz+f[i]+list_lf[i]/2),cos(pi)*(list_da[i]/2-1+0.1)),), ((sin(3*pi/2)*(list_da[i]/2-1+0.1),l_absatz+radius[i],cos(3*pi/2)*(list_da[i]/2-1+0.1)),), ((sin(3*pi/2)*(list_da[i]/2-1+0.1),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos(3*pi/2)*(list_da[i]/2-1+0.1)),), ((sin(3*pi/2)*(list_da[i]/2-1+0.1),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos(3*pi/2)*(list_da[i]/2-1+0.1)),), ((sin(3*pi/2)*(list_da[i]/2-1+0.1),l_absatz+f[i],cos(3*pi/2)*(list_da[i]/2-1+0.1)),), ((sin(3*pi/2)*(list_da[i]/2-1+0.1),l_absatz+f[i]+4,cos(3*pi/2)*(list_da[i]/2-1+0.1)),), ((sin(3*pi/2)*(list_da[i]/2-1+0.1),(l_absatz+f[i]+list_lf[i]),cos(3*pi/2)*(list_da[i]/2-1+0.1)),), ((sin(3*pi/2)*(list_da[i]/2-1+0.1),(l_absatz+f[i]+list_lf[i]/2),cos(3*pi/2)*(list_da[i]/2-1+0.1)),), ((sin(2*pi)*(list_da[i]/2-1+0.1),l_absatz+radius[i],cos(2*pi)*(list_da[i]/2-1+0.1)),), ((sin(2*pi)*(list_da[i]/2-1+0.1),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos(2*pi)*(list_da[i]/2-1+0.1)),), ((sin(2*pi)*(list_da[i]/2-1+0.1),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos(2*pi)*(list_da[i]/2-1+0.1)),), ((sin(2*pi)*(list_da[i]/2-1+0.1),l_absatz+f[i],cos(2*pi)*(list_da[i]/2-1+0.1)),), ((sin(2*pi)*(list_da[i]/2-1+0.1),l_absatz+f[i]+4,cos(2*pi)*(list_da[i]/2-1+0.1)),), ((sin(2*pi)*(list_da[i]/2-1+0.1),(l_absatz+f[i]+list_lf[i]),cos(2*pi)*(list_da[i]/2-1+0.1)),), ((sin(2*pi)*(list_da[i]/2-1+0.1),(l_absatz+f[i]+list_lf[i]/2),cos(2*pi)*(list_da[i]/2-1+0.1)),),) welle.Set(edges=Radialedge_FR_Fugenflaeche_a, name='Radial_00_Freistich+Fugenflaeche_aussen') Radialedge_Rundung = e_w.findAt(((sin(pi/2)*(list_da[i]/2),l_absatz-2,cos(pi/2)*(list_da[i]/2)),), ((sin(pi/2)*(list_da[i]/2),l_absatz-20,cos(pi/2)*(list_da[i]/2)),), ((sin(pi)*(list_da[i]/2),l_absatz-2,cos(pi)*(list_da[i]/2)),), ((sin(pi)*(list_da[i]/2),l_absatz-20,cos(pi)*(list_da[i]/2)),), ((sin(3*pi/2)*(list_da[i]/2),l_absatz-2,cos(3*pi/2)*(list_da[i]/2)),), ((sin(3*pi/2)*(list_da[i]/2),l_absatz-20,cos(3*pi/2)*(list_da[i]/2)),), ((sin(2*pi)*(list_da[i]/2),l_absatz-2,cos(2*pi)*(list_da[i]/2)),), ((sin(2*pi)*(list_da[i]/2),l_absatz-20,cos(2*pi)*(list_da[i]/2)),),) welle.Set(edges=Radialedge_Rundung, name='Radial_01_Rundung') Radialedge_FR_Fugenflaeche_i = e_w.findAt(((sin(pi/2)*(list_da[i]/2-1-0.1),l_absatz+radius[i],cos(pi/2)*(list_da[i]/2-1-0.1)),), ((sin(pi/2)*(list_da[i]/2-1-0.1),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos(pi/2)*(list_da[i]/2-1-0.1)),), ((sin(pi/2)*(list_da[i]/2-1-0.1),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos(pi/2)*(list_da[i]/2-1-0.1)),), ((sin(pi/2)*(list_da[i]/2-1-0.1),l_absatz+f[i],cos(pi/2)*(list_da[i]/2-1-0.1)),), ((sin(pi/2)*(list_da[i]/2-1-0.1),l_absatz+f[i]+4,cos(pi/2)*(list_da[i]/2-1-0.1)),), ((sin(pi/2)*(list_da[i]/2-1-0.1),(l_absatz+f[i]+list_lf[i]),cos(pi/2)*(list_da[i]/2-1-0.1)),), ((sin(pi/2)*(list_da[i]/2-1-0.1),(l_absatz+f[i]+list_lf[i]/2),cos(pi/2)*(list_da[i]/2-1-0.1)),), ((sin(pi/2)*(list_da[i]/2-1-0.1),l_absatz-2,cos(pi/2)*(list_da[i]/2-1-0.1)),), ((sin(pi/2)*(list_da[i]/2-1-0.1),l_absatz-20,cos(pi/2)*(list_da[i]/2-1-0.1)),), ((sin(pi/2)*(list_da[i]/2-1-0.1),lw,cos(pi/2)*(list_da[i]/2-1-0.1)),), ((sin(pi)*(list_da[i]/2-1-0.1),l_absatz+radius[i],cos(pi)*(list_da[i]/2-1-0.1)),), ((sin(pi)*(list_da[i]/2-1-0.1),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos(pi)*(list_da[i]/2-1-0.1)),), ((sin(pi)*(list_da[i]/2-1-0.1),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos(pi)*(list_da[i]/2-1-0.1)),), ((sin(pi)*(list_da[i]/2-1-0.1),l_absatz+f[i],cos(pi)*(list_da[i]/2-1-0.1)),), ((sin(pi)*(list_da[i]/2-1-0.1),l_absatz+f[i]+4,cos(pi)*(list_da[i]/2-1-0.1)),), ((sin(pi)*(list_da[i]/2-1-0.1),(l_absatz+f[i]+list_lf[i]),cos(pi)*(list_da[i]/2-1-0.1)),), ((sin(pi)*(list_da[i]/2-1-0.1),(l_absatz+f[i]+list_lf[i]/2),cos(pi)*(list_da[i]/2-1-0.1)),), ((sin(pi)*(list_da[i]/2-1-0.1),l_absatz-2,cos(pi)*(list_da[i]/2-1-0.1)),), ((sin(pi)*(list_da[i]/2-1-0.1),l_absatz-20,cos(pi)*(list_da[i]/2-1-0.1)),), ((sin(pi)*(list_da[i]/2-1-0.1),lw,cos(pi)*(list_da[i]/2-1-0.1)),), ((sin(3*pi/2)*(list_da[i]/2-1-0.1),l_absatz+radius[i],cos(3*pi/2)*(list_da[i]/2-1-0.1)),), ((sin(3*pi/2)*(list_da[i]/2-1-0.1),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos(3*pi/2)*(list_da[i]/2-1-0.1)),), ((sin(3*pi/2)*(list_da[i]/2-1-0.1),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos(3*pi/2)*(list_da[i]/2-1-0.1)),), ((sin(3*pi/2)*(list_da[i]/2-1-0.1),l_absatz+f[i],cos(3*pi/2)*(list_da[i]/2-1-0.1)),), ((sin(3*pi/2)*(list_da[i]/2-1-0.1),l_absatz+f[i]+4,cos(3*pi/2)*(list_da[i]/2-1-0.1)),), ((sin(3*pi/2)*(list_da[i]/2-1-0.1),(l_absatz+f[i]+list_lf[i]),cos(3*pi/2)*(list_da[i]/2-1-0.1)),), ((sin(3*pi/2)*(list_da[i]/2-1-0.1),(l_absatz+f[i]+list_lf[i]/2),cos(3*pi/2)*(list_da[i]/2-1-0.1)),), ((sin(3*pi/2)*(list_da[i]/2-1-0.1),l_absatz-2,cos(3*pi/2)*(list_da[i]/2-1-0.1)),), ((sin(3*pi/2)*(list_da[i]/2-1-0.1),l_absatz-20,cos(3*pi/2)*(list_da[i]/2-1-0.1)),), ((sin(3*pi/2)*(list_da[i]/2-1-0.1),lw,cos(3*pi/2)*(list_da[i]/2-1-0.1)),), ((sin(2*pi)*(list_da[i]/2-1-0.1),l_absatz+radius[i],cos(2*pi)*(list_da[i]/2-1-0.1)),), ((sin(2*pi)*(list_da[i]/2-1-0.1),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos(2*pi)*(list_da[i]/2-1-0.1)),), ((sin(2*pi)*(list_da[i]/2-1-0.1),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos(2*pi)*(list_da[i]/2-1-0.1)),), ((sin(2*pi)*(list_da[i]/2-1-0.1),l_absatz+f[i],cos(2*pi)*(list_da[i]/2-1-0.1)),), ((sin(2*pi)*(list_da[i]/2-1-0.1),l_absatz+f[i]+4,cos(2*pi)*(list_da[i]/2-1-0.1)),), ((sin(2*pi)*(list_da[i]/2-1-0.1),(l_absatz+f[i]+list_lf[i]),cos(2*pi)*(list_da[i]/2-1-0.1)),), ((sin(2*pi)*(list_da[i]/2-1-0.1),(l_absatz+f[i]+list_lf[i]/2),cos(2*pi)*(list_da[i]/2-1-0.1)),), ((sin(2*pi)*(list_da[i]/2-1-0.1),l_absatz-2,cos(2*pi)*(list_da[i]/2-1-0.1)),), ((sin(2*pi)*(list_da[i]/2-1-0.1),l_absatz-20,cos(2*pi)*(list_da[i]/2-1-0.1)),), ((sin(2*pi)*(list_da[i]/2-1-0.1),lw,cos(2*pi)*(list_da[i]/2-1-0.1)),),) welle.Set(edges=Radialedge_FR_Fugenflaeche_i, name='Radial_02_Freistich+Fugenflaeche_innen') Radialedge_Anfangsteil = e_w.findAt(((sin(pi/2)*(list_da[i]/2+h_absatz[i]-0.01),0,cos(pi/2)*(list_da[i]/2+h_absatz[i]-0.01)),), ((sin(pi/2)*(list_da[i]/2+radius[i]-0.01),0,cos(pi/2)*(list_da[i]/2+radius[i]-0.01)),), ((sin(pi/2)*(list_da[i]/2-1-0.01),0,cos(pi/2)*(list_da[i]/2-1-0.01)),), ((sin(pi/2)*(list_di[i]/2+0.01),0,cos(pi/2)*(list_di[i]/2+0.01)),), ((sin(pi)*(list_da[i]/2+h_absatz[i]-0.01),0,cos(pi)*(list_da[i]/2+h_absatz[i]-0.01)),), ((sin(pi)*(list_da[i]/2+radius[i]-0.01),0,cos(pi)*(list_da[i]/2+radius[i]-0.01)),), ((sin(pi)*(list_da[i]/2-1-0.01),0,cos(pi)*(list_da[i]/2-1-0.01)),), ((sin(pi)*(list_di[i]/2+0.01),0,cos(pi)*(list_di[i]/2+0.01)),), ((sin(3*pi/2)*(list_da[i]/2+h_absatz[i]-0.01),0,cos(3*pi/2)*(list_da[i]/2+h_absatz[i]-0.01)),), ((sin(3*pi/2)*(list_da[i]/2+radius[i]-0.01),0,cos(3*pi/2)*(list_da[i]/2+radius[i]-0.01)),), ((sin(3*pi/2)*(list_da[i]/2-1-0.01),0,cos(3*pi/2)*(list_da[i]/2-1-0.01)),), ((sin(3*pi/2)*(list_di[i]/2+0.01),0,cos(3*pi/2)*(list_di[i]/2+0.01)),), ((sin(2*pi)*(list_da[i]/2+h_absatz[i]-0.01),0,cos(2*pi)*(list_da[i]/2+h_absatz[i]-0.01)),), ((sin(2*pi)*(list_da[i]/2+radius[i]-0.01),0,cos(2*pi)*(list_da[i]/2+radius[i]-0.01)),), ((sin(2*pi)*(list_da[i]/2-1-0.01),0,cos(2*pi)*(list_da[i]/2-1-0.01)),), ((sin(2*pi)*(list_di[i]/2+0.01),0,cos(2*pi)*(list_di[i]/2+0.01)),),) welle.Set(edges=Radialedge_Anfangsteil, name='Radial_03_Anfangsteil') Radialedge_Absatzkante = e_w.findAt(((sin(pi/2)*(list_da[i]/2+radius[i]+0.01),l_absatz-2,cos(pi/2)*(list_da[i]/2+radius[i]+0.01)),), ((sin(pi/2)*(list_da[i]/2+radius[i]+0.01),l_absatz-20,cos(pi/2)*(list_da[i]/2+radius[i]+0.01)),), ((sin(pi/2)*(list_da[i]/2+radius[i]+0.01),l_absatz,cos(pi/2)*(list_da[i]/2+radius[i]+0.01)),), ((sin(pi)*(list_da[i]/2+radius[i]+0.01),l_absatz-2,cos(pi)*(list_da[i]/2+radius[i]+0.01)),), ((sin(pi)*(list_da[i]/2+radius[i]+0.01),l_absatz-20,cos(pi)*(list_da[i]/2+radius[i]+0.01)),), ((sin(pi)*(list_da[i]/2+radius[i]+0.01),l_absatz,cos(pi)*(list_da[i]/2+radius[i]+0.01)),), ((sin(3*pi/2)*(list_da[i]/2+radius[i]+0.01),l_absatz-2,cos(3*pi/2)*(list_da[i]/2+radius[i]+0.01)),), ((sin(3*pi/2)*(list_da[i]/2+radius[i]+0.01),l_absatz-20,cos(3*pi/2)*(list_da[i]/2+radius[i]+0.01)),), ((sin(3*pi/2)*(list_da[i]/2+radius[i]+0.01),l_absatz,cos(3*pi/2)*(list_da[i]/2+radius[i]+0.01)),), ((sin(2*pi)*(list_da[i]/2+radius[i]+0.01),l_absatz-2,cos(2*pi)*(list_da[i]/2+radius[i]+0.01)),), ((sin(2*pi)*(list_da[i]/2+radius[i]+0.01),l_absatz-20,cos(2*pi)*(list_da[i]/2+radius[i]+0.01)),), ((sin(2*pi)*(list_da[i]/2+radius[i]+0.01),l_absatz,cos(2*pi)*(list_da[i]/2+radius[i]+0.01)),),) welle.Set(edges=Radialedge_Absatzkante, name='Radial_04_Absatzkante') Radialedge_Auslauf = e_w.findAt(((sin(pi/2)*(list_da[i]/2-0.01),lw,cos(pi/2)*(list_da[i]/2-0.01)),), ((sin(pi)*(list_da[i]/2-0.01),lw,cos(pi)*(list_da[i]/2-0.01)),), ((sin(3*pi/2)*(list_da[i]/2-0.01),lw,cos(3*pi/2)*(list_da[i]/2-0.01)),), ((sin(2*pi)*(list_da[i]/2-0.01),lw,cos(2*pi)*(list_da[i]/2-0.01)),),) welle.Set(edges=Radialedge_Auslauf, name='Radial_05_Auslauf') # Sets erstellen # Nabe #tangential nabe = mdb.models[Modell_Name[i]].parts[Naben_Name[i]] session.viewports['Viewport: 1'].setValues(displayedObject=nabe) e_n = nabe.edges Tangentialedge_Nabe = e_n.findAt(((sin (pi/4)*list_Da[i]/2, l_absatz+1.5,cos (pi/4)*list_Da[i]/2),), ((sin (pi/4)*(list_Di[i]/2+1), l_absatz+1.5,cos (pi/4)*(list_Di[i]/2+1)),), ((sin (pi/4)*list_Di[i]/2, l_absatz+1.5,cos (pi/4)*list_Di[i]/2),), ((sin (pi/4)*list_Da[i]/2, l_absatz+f[i],cos (pi/4)*list_Da[i]/2),), ((sin (pi/4)*(list_Di[i]/2+1), l_absatz+f[i],cos (pi/4)*(list_Di[i]/2+1)),), ((sin (pi/4)*list_Di[i]/2, l_absatz+f[i],cos (pi/4)*list_Di[i]/2),), ((sin (pi/4)*list_Da[i]/2, l_absatz+f[i]+4,cos (pi/4)*list_Da[i]/2),), ((sin (pi/4)*(list_Di[i]/2+1), l_absatz+f[i]+4,cos (pi/4)*(list_Di[i]/2+1)),), ((sin (pi/4)*list_Di[i]/2, l_absatz+f[i]+4,cos (pi/4)*list_Di[i]/2),), ((sin (pi/4)*list_Da[i]/2, l_absatz+f[i]+list_lf[i],cos (pi/4)*list_Da[i]/2),), ((sin (pi/4)*(list_Di[i]/2+1), l_absatz+f[i]+list_lf[i],cos (pi/4)*(list_Di[i]/2+1)),), ((sin (pi/4)*list_Di[i]/2, l_absatz+f[i]+list_lf[i],cos (pi/4)*list_Di[i]/2),), ((sin (pi/4)*list_Da[i]/2, l_absatz+f[i]+list_lf[i]/2,cos (pi/4)*list_Da[i]/2),), ((sin (pi/4)*(list_Di[i]/2+1), l_absatz+f[i]+list_lf[i]/2,cos (pi/4)*(list_Di[i]/2+1)),), ((sin (pi/4)*list_Di[i]/2, l_absatz+f[i]+list_lf[i]/2,cos (pi/4)*list_Di[i]/2),), ((sin (3*pi/4)*list_Da[i]/2, l_absatz+1.5,cos (3*pi/4)*list_Da[i]/2),), ((sin (3*pi/4)*(list_Di[i]/2+1), l_absatz+1.5,cos (3*pi/4)*(list_Di[i]/2+1)),), ((sin (3*pi/4)*list_Di[i]/2, l_absatz+1.5,cos (3*pi/4)*list_Di[i]/2),), ((sin (3*pi/4)*list_Da[i]/2, l_absatz+f[i],cos (3*pi/4)*list_Da[i]/2),), ((sin (3*pi/4)*(list_Di[i]/2+1), l_absatz+f[i],cos (3*pi/4)*(list_Di[i]/2+1)),), ((sin (3*pi/4)*list_Di[i]/2, l_absatz+f[i],cos (3*pi/4)*list_Di[i]/2),), ((sin (3*pi/4)*list_Da[i]/2, l_absatz+f[i]+4,cos (3*pi/4)*list_Da[i]/2),), ((sin (3*pi/4)*(list_Di[i]/2+1), l_absatz+f[i]+4,cos (3*pi/4)*(list_Di[i]/2+1)),), ((sin (3*pi/4)*list_Di[i]/2, l_absatz+f[i]+4,cos (3*pi/4)*list_Di[i]/2),), ((sin (3*pi/4)*list_Da[i]/2, l_absatz+f[i]+list_lf[i],cos (3*pi/4)*list_Da[i]/2),), ((sin (3*pi/4)*(list_Di[i]/2+1), l_absatz+f[i]+list_lf[i],cos (3*pi/4)*(list_Di[i]/2+1)),), ((sin (3*pi/4)*list_Di[i]/2, l_absatz+f[i]+list_lf[i],cos (3*pi/4)*list_Di[i]/2),), ((sin (3*pi/4)*list_Da[i]/2, l_absatz+f[i]+list_lf[i]/2,cos (3*pi/4)*list_Da[i]/2),), ((sin (3*pi/4)*(list_Di[i]/2+1), l_absatz+f[i]+list_lf[i]/2,cos (3*pi/4)*(list_Di[i]/2+1)),), ((sin (3*pi/4)*list_Di[i]/2, l_absatz+f[i]+list_lf[i]/2,cos (3*pi/4)*list_Di[i]/2),), ((sin (5*pi/4)*list_Da[i]/2, l_absatz+1.5,cos (5*pi/4)*list_Da[i]/2),), ((sin (5*pi/4)*(list_Di[i]/2+1), l_absatz+1.5,cos (5*pi/4)*(list_Di[i]/2+1)),), ((sin (5*pi/4)*list_Di[i]/2, l_absatz+1.5,cos (5*pi/4)*list_Di[i]/2),), ((sin (5*pi/4)*list_Da[i]/2, l_absatz+f[i],cos (5*pi/4)*list_Da[i]/2),), ((sin (5*pi/4)*(list_Di[i]/2+1), l_absatz+f[i],cos (5*pi/4)*(list_Di[i]/2+1)),), ((sin (5*pi/4)*list_Di[i]/2, l_absatz+f[i],cos (5*pi/4)*list_Di[i]/2),), ((sin (5*pi/4)*list_Da[i]/2, l_absatz+f[i]+4,cos (5*pi/4)*list_Da[i]/2),), ((sin (5*pi/4)*(list_Di[i]/2+1), l_absatz+f[i]+4,cos (5*pi/4)*(list_Di[i]/2+1)),), ((sin (5*pi/4)*list_Di[i]/2, l_absatz+f[i]+4,cos (5*pi/4)*list_Di[i]/2),), ((sin (5*pi/4)*list_Da[i]/2, l_absatz+f[i]+list_lf[i],cos (5*pi/4)*list_Da[i]/2),), ((sin (5*pi/4)*(list_Di[i]/2+1), l_absatz+f[i]+list_lf[i],cos (5*pi/4)*(list_Di[i]/2+1)),), ((sin (5*pi/4)*list_Di[i]/2, l_absatz+f[i]+list_lf[i],cos (5*pi/4)*list_Di[i]/2),), ((sin (5*pi/4)*list_Da[i]/2, l_absatz+f[i]+list_lf[i]/2,cos (5*pi/4)*list_Da[i]/2),), ((sin (5*pi/4)*(list_Di[i]/2+1), l_absatz+f[i]+list_lf[i]/2,cos (5*pi/4)*(list_Di[i]/2+1)),), ((sin (5*pi/4)*list_Di[i]/2, l_absatz+f[i]+list_lf[i]/2,cos (5*pi/4)*list_Di[i]/2),), ((sin (7*pi/4)*list_Da[i]/2, l_absatz+1.5,cos (7*pi/4)*list_Da[i]/2),), ((sin (7*pi/4)*(list_Di[i]/2+1), l_absatz+1.5,cos (7*pi/4)*(list_Di[i]/2+1)),), ((sin (7*pi/4)*list_Di[i]/2, l_absatz+1.5,cos (7*pi/4)*list_Di[i]/2),), ((sin (7*pi/4)*list_Da[i]/2, l_absatz+f[i],cos (7*pi/4)*list_Da[i]/2),), ((sin (7*pi/4)*(list_Di[i]/2+1), l_absatz+f[i],cos (7*pi/4)*(list_Di[i]/2+1)),), ((sin (7*pi/4)*list_Di[i]/2, l_absatz+f[i],cos (7*pi/4)*list_Di[i]/2),), ((sin (7*pi/4)*list_Da[i]/2, l_absatz+f[i]+4,cos (7*pi/4)*list_Da[i]/2),), ((sin (7*pi/4)*(list_Di[i]/2+1), l_absatz+f[i]+4,cos (7*pi/4)*(list_Di[i]/2+1)),), ((sin (7*pi/4)*list_Di[i]/2, l_absatz+f[i]+4,cos (7*pi/4)*list_Di[i]/2),), ((sin (7*pi/4)*list_Da[i]/2, l_absatz+f[i]+list_lf[i],cos (7*pi/4)*list_Da[i]/2),), ((sin (7*pi/4)*(list_Di[i]/2+1), l_absatz+f[i]+list_lf[i],cos (7*pi/4)*(list_Di[i]/2+1)),), ((sin (7*pi/4)*list_Di[i]/2, l_absatz+f[i]+list_lf[i],cos (7*pi/4)*list_Di[i]/2),), ((sin (7*pi/4)*list_Da[i]/2, l_absatz+f[i]+list_lf[i]/2,cos (7*pi/4)*list_Da[i]/2),), ((sin (7*pi/4)*(list_Di[i]/2+1), l_absatz+f[i]+list_lf[i]/2,cos (7*pi/4)*(list_Di[i]/2+1)),), ((sin (7*pi/4)*list_Di[i]/2, l_absatz+f[i]+list_lf[i]/2,cos (7*pi/4)*list_Di[i]/2),),) nabe.Set(edges=Tangentialedge_Nabe, name='Tangential_01') #axial Axialedge_Kontaktfrei = e_n.findAt(((sin (pi/2)*list_Di[i]/2, l_absatz+1.5+0.01, cos (pi/2)*list_Di[i]/2),), ((sin (pi/2)*(list_Di[i]/2+1), l_absatz+1.5+0.01, cos (pi/2)*(list_Di[i]/2+1)),), ((sin (pi/2)*list_Da[i]/2, l_absatz+1.5+0.01, cos (pi/2)*list_Da[i]/2),), ((sin (pi)*list_Di[i]/2, l_absatz+1.5+0.01, cos (pi)*list_Di[i]/2),), ((sin (pi)*(list_Di[i]/2+1), l_absatz+1.5+0.01, cos (pi)*(list_Di[i]/2+1)),), ((sin (pi)*list_Da[i]/2, l_absatz+1.5+0.01, cos (pi)*list_Da[i]/2),), ((sin (3*pi/2)*list_Di[i]/2, l_absatz+1.5+0.01, cos (3*pi/2)*list_Di[i]/2),), ((sin (3*pi/2)*(list_Di[i]/2+1), l_absatz+1.5+0.01, cos (3*pi/2)*(list_Di[i]/2+1)),), ((sin (3*pi/2)*list_Da[i]/2, l_absatz+1.5+0.01, cos (3*pi/2)*list_Da[i]/2),), ((sin (2*pi)*list_Di[i]/2, l_absatz+1.5+0.01, cos (2*pi)*list_Di[i]/2),), ((sin (2*pi)*(list_Di[i]/2+1), l_absatz+1.5+0.01, cos (2*pi)*(list_Di[i]/2+1)),), ((sin (2*pi)*list_Da[i]/2, l_absatz+1.5+0.01, cos (2*pi)*list_Da[i]/2),),) nabe.Set(edges=Axialedge_Kontaktfrei, name='Axial_00_Kontaktfrei') Axialedge_Kontaktflaeche_u = e_n.findAt(((sin (pi/2)*list_Di[i]/2, l_absatz+f[i]+6, cos (pi/2)*list_Di[i]/2),), ((sin (pi/2)*(list_Di[i]/2+1), l_absatz+f[i]+6, cos (pi/2)*(list_Di[i]/2+1)),), ((sin (pi/2)*list_Da[i]/2, l_absatz+f[i]+6, cos (pi/2)*list_Da[i]/2),), ((sin (pi)*list_Di[i]/2, l_absatz+f[i]+6, cos (pi)*list_Di[i]/2),), ((sin (pi)*(list_Di[i]/2+1), l_absatz+f[i]+6, cos (pi)*(list_Di[i]/2+1)),), ((sin (pi)*list_Da[i]/2, l_absatz+f[i]+6, cos (pi)*list_Da[i]/2),), ((sin (3*pi/2)*list_Di[i]/2, l_absatz+f[i]+6, cos (3*pi/2)*list_Di[i]/2),), ((sin (3*pi/2)*(list_Di[i]/2+1), l_absatz+f[i]+6, cos (3*pi/2)*(list_Di[i]/2+1)),), ((sin (3*pi/2)*list_Da[i]/2, l_absatz+f[i]+6, cos (3*pi/2)*list_Da[i]/2),), ((sin (2*pi)*list_Di[i]/2, l_absatz+f[i]+6, cos (2*pi)*list_Di[i]/2),), ((sin (2*pi)*(list_Di[i]/2+1), l_absatz+f[i]+6, cos (2*pi)*(list_Di[i]/2+1)),), ((sin (2*pi)*list_Da[i]/2, l_absatz+f[i]+6, cos (2*pi)*list_Da[i]/2),),) nabe.Set(edges=Axialedge_Kontaktflaeche_u, name='Axial_01_Kontaktflaeche_unten') Axialedge_Kontaktflaeche_o = e_n.findAt(((sin (pi/2)*list_Di[i]/2, l_absatz+f[i]+list_lf[i]-2, cos (pi/2)*list_Di[i]/2),), ((sin (pi/2)*(list_Di[i]/2+1), l_absatz+f[i]+list_lf[i]-2, cos (pi/2)*(list_Di[i]/2+1)),), ((sin (pi/2)*list_Da[i]/2, l_absatz+f[i]+list_lf[i]-2, cos (pi/2)*list_Da[i]/2),), ((sin (pi)*list_Di[i]/2, l_absatz+f[i]+list_lf[i]-2, cos (pi)*list_Di[i]/2),), ((sin (pi)*(list_Di[i]/2+1), l_absatz+f[i]+list_lf[i]-2, cos (pi)*(list_Di[i]/2+1)),), ((sin (pi)*list_Da[i]/2, l_absatz+f[i]+list_lf[i]-2, cos (pi)*list_Da[i]/2),), ((sin (3*pi/2)*list_Di[i]/2, l_absatz+f[i]+list_lf[i]-2, cos (3*pi/2)*list_Di[i]/2),), ((sin (3*pi/2)*(list_Di[i]/2+1), l_absatz+f[i]+list_lf[i]-2, cos (3*pi/2)*(list_Di[i]/2+1)),), ((sin (3*pi/2)*list_Da[i]/2, l_absatz+f[i]+list_lf[i]-2, cos (3*pi/2)*list_Da[i]/2),), ((sin (2*pi)*list_Di[i]/2, l_absatz+f[i]+list_lf[i]-2, cos (2*pi)*list_Di[i]/2),), ((sin (2*pi)*(list_Di[i]/2+1), l_absatz+f[i]+list_lf[i]-2, cos (2*pi)*(list_Di[i]/2+1)),), ((sin (2*pi)*list_Da[i]/2, l_absatz+f[i]+list_lf[i]-2, cos (2*pi)*list_Da[i]/2),),) nabe.Set(edges=Axialedge_Kontaktflaeche_o, name='Axial_02_Kontaktflaeche_oben') Axialedge_Kontaktflaeche_Anfang = e_n.findAt(((sin (pi/2)*list_Di[i]/2, l_absatz+f[i]+0.01, cos (pi/2)*list_Di[i]/2),), ((sin (pi/2)*(list_Di[i]/2+1), l_absatz+f[i]+0.01, cos (pi/2)*(list_Di[i]/2+1)),), ((sin (pi/2)*list_Da[i]/2, l_absatz+f[i]+0.01, cos (pi/2)*list_Da[i]/2),), ((sin (pi)*list_Di[i]/2, l_absatz+f[i]+0.01, cos (pi)*list_Di[i]/2),), ((sin (pi)*(list_Di[i]/2+1), l_absatz+f[i]+0.01, cos (pi)*(list_Di[i]/2+1)),), ((sin (pi)*list_Da[i]/2, l_absatz+f[i]+0.01, cos (pi)*list_Da[i]/2),), ((sin (3*pi/2)*list_Di[i]/2, l_absatz+f[i]+0.01, cos (3*pi/2)*list_Di[i]/2),), ((sin (3*pi/2)*(list_Di[i]/2+1), l_absatz+f[i]+0.01, cos (3*pi/2)*(list_Di[i]/2+1)),), ((sin (3*pi/2)*list_Da[i]/2, l_absatz+f[i]+0.01, cos (3*pi/2)*list_Da[i]/2),), ((sin (2*pi)*list_Di[i]/2, l_absatz+f[i]+0.01, cos (2*pi)*list_Di[i]/2),), ((sin (2*pi)*(list_Di[i]/2+1), l_absatz+f[i]+0.01, cos (2*pi)*(list_Di[i]/2+1)),), ((sin (2*pi)*list_Da[i]/2, l_absatz+f[i]+0.01, cos (2*pi)*list_Da[i]/2),),) nabe.Set(edges=Axialedge_Kontaktflaeche_Anfang, name='Axial_03_Kontaktflaeche_Anfang') #radial Radialedge_innen = e_n.findAt(((sin (pi/2)*(list_Di[i]/2+0.01), l_absatz+1.5, cos (pi/2)*(list_Di[i]/2+0.01)),), ((sin (pi/2)*(list_Di[i]/2+0.01), l_absatz+f[i], cos (pi/2)*(list_Di[i]/2+0.01)),), ((sin (pi/2)*(list_Di[i]/2+0.01), l_absatz+f[i]+4, cos (pi/2)*(list_Di[i]/2+0.01)),), ((sin (pi/2)*(list_Di[i]/2+0.01), l_absatz+f[i]+list_lf[i], cos (pi/2)*(list_Di[i]/2+0.01)),), ((sin (pi/2)*(list_Di[i]/2+0.01), l_absatz+f[i]+list_lf[i]/2, cos (pi/2)*(list_Di[i]/2+0.01)),), ((sin (pi)*(list_Di[i]/2+0.01), l_absatz+1.5, cos (pi)*(list_Di[i]/2+0.01)),), ((sin (pi)*(list_Di[i]/2+0.01), l_absatz+f[i], cos (pi)*(list_Di[i]/2+0.01)),), ((sin (pi)*(list_Di[i]/2+0.01), l_absatz+f[i]+4, cos (pi)*(list_Di[i]/2+0.01)),), ((sin (pi)*(list_Di[i]/2+0.01), l_absatz+f[i]+list_lf[i], cos (pi)*(list_Di[i]/2+0.01)),), ((sin (pi)*(list_Di[i]/2+0.01), l_absatz+f[i]+list_lf[i]/2, cos (pi)*(list_Di[i]/2+0.01)),), ((sin (3*pi/2)*(list_Di[i]/2+0.01), l_absatz+1.5, cos (3*pi/2)*(list_Di[i]/2+0.01)),), ((sin (3*pi/2)*(list_Di[i]/2+0.01), l_absatz+f[i], cos (3*pi/2)*(list_Di[i]/2+0.01)),), ((sin (3*pi/2)*(list_Di[i]/2+0.01), l_absatz+f[i]+4, cos (3*pi/2)*(list_Di[i]/2+0.01)),), ((sin (3*pi/2)*(list_Di[i]/2+0.01), l_absatz+f[i]+list_lf[i], cos (3*pi/2)*(list_Di[i]/2+0.01)),), ((sin (3*pi/2)*(list_Di[i]/2+0.01), l_absatz+f[i]+list_lf[i]/2, cos (3*pi/2)*(list_Di[i]/2+0.01)),), ((sin (2*pi)*(list_Di[i]/2+0.01), l_absatz+1.5, cos (2*pi)*(list_Di[i]/2+0.01)),), ((sin (2*pi)*(list_Di[i]/2+0.01), l_absatz+f[i], cos (2*pi)*(list_Di[i]/2+0.01)),), ((sin (2*pi)*(list_Di[i]/2+0.01), l_absatz+f[i]+4, cos (2*pi)*(list_Di[i]/2+0.01)),), ((sin (2*pi)*(list_Di[i]/2+0.01), l_absatz+f[i]+list_lf[i], cos (2*pi)*(list_Di[i]/2+0.01)),), ((sin (2*pi)*(list_Di[i]/2+0.01), l_absatz+f[i]+list_lf[i]/2, cos (2*pi)*(list_Di[i]/2+0.01)),),) nabe.Set(edges=Radialedge_innen, name='Radial_00_innen') Radialedge_aussen = e_n.findAt(((sin (pi/2)*(list_Da[i]/2-0.01), l_absatz+1.5, cos (pi/2)*(list_Da[i]/2-0.01)),), ((sin (pi/2)*(list_Da[i]/2-0.01), l_absatz+f[i], cos (pi/2)*(list_Da[i]/2-0.01)),), ((sin (pi/2)*(list_Da[i]/2-0.01), l_absatz+f[i]+4, cos (pi/2)*(list_Da[i]/2-0.01)),), ((sin (pi/2)*(list_Da[i]/2-0.01), l_absatz+f[i]+list_lf[i], cos (pi/2)*(list_Da[i]/2-0.01)),), ((sin (pi/2)*(list_Da[i]/2-0.01), l_absatz+f[i]+list_lf[i]/2, cos (pi/2)*(list_Da[i]/2-0.01)),), ((sin (pi)*(list_Da[i]/2-0.01), l_absatz+1.5, cos (pi)*(list_Da[i]/2-0.01)),), ((sin (pi)*(list_Da[i]/2-0.01), l_absatz+f[i], cos (pi)*(list_Da[i]/2-0.01)),), ((sin (pi)*(list_Da[i]/2-0.01), l_absatz+f[i]+4, cos (pi)*(list_Da[i]/2-0.01)),), ((sin (pi)*(list_Da[i]/2-0.01), l_absatz+f[i]+list_lf[i], cos (pi)*(list_Da[i]/2-0.01)),), ((sin (pi)*(list_Da[i]/2-0.01), l_absatz+f[i]+list_lf[i]/2, cos (pi)*(list_Da[i]/2-0.01)),), ((sin (3*pi/2)*(list_Da[i]/2-0.01), l_absatz+1.5, cos (3*pi/2)*(list_Da[i]/2-0.01)),), ((sin (3*pi/2)*(list_Da[i]/2-0.01), l_absatz+f[i], cos (3*pi/2)*(list_Da[i]/2-0.01)),), ((sin (3*pi/2)*(list_Da[i]/2-0.01), l_absatz+f[i]+4, cos (3*pi/2)*(list_Da[i]/2-0.01)),), ((sin (3*pi/2)*(list_Da[i]/2-0.01), l_absatz+f[i]+list_lf[i], cos (3*pi/2)*(list_Da[i]/2-0.01)),), ((sin (3*pi/2)*(list_Da[i]/2-0.01), l_absatz+f[i]+list_lf[i]/2, cos (3*pi/2)*(list_Da[i]/2-0.01)),), ((sin (2*pi)*(list_Da[i]/2-0.01), l_absatz+1.5, cos (2*pi)*(list_Da[i]/2-0.01)),), ((sin (2*pi)*(list_Da[i]/2-0.01), l_absatz+f[i], cos (2*pi)*(list_Da[i]/2-0.01)),), ((sin (2*pi)*(list_Da[i]/2-0.01), l_absatz+f[i]+4, cos (2*pi)*(list_Da[i]/2-0.01)),), ((sin (2*pi)*(list_Da[i]/2-0.01), l_absatz+f[i]+list_lf[i], cos (2*pi)*(list_Da[i]/2-0.01)),), ((sin (2*pi)*(list_Da[i]/2-0.01), l_absatz+f[i]+list_lf[i]/2, cos (2*pi)*(list_Da[i]/2-0.01)),),) nabe.Set(edges=Radialedge_aussen, name='Radial_01_aussen') # Sets erstellen # Welle innen if i > 8 and i < 13: #tangential welle_innen = mdb.models[Modell_Name[i]].parts[Welle_innen_Name[i]] session.viewports['Viewport: 1'].setValues(displayedObject=welle_innen) e_wi = welle_innen.edges Tangentialedge_Wi = e_wi.findAt(((sin (pi/4)*da_i/2, lw,cos (pi/4)*da_i/2),), ((sin (pi/4)*da_i/2, (l_absatz+f[i]+list_lf[i]),cos (pi/4)*da_i/2),), ((sin (pi/4)*da_i/2, (l_absatz+f[i]+list_lf[i]/2),cos (pi/4)*da_i/2),), ((sin (pi/4)*da_i/2, (l_absatz+f[i]),cos (pi/4)*da_i/2),), ((sin (pi/4)*da_i/2, (l_absatz+f[i]+4),cos (pi/4)*da_i/2),), ((sin (pi/4)*da_i/2, (l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos (pi/4)*da_i/2),), ((sin (pi/4)*da_i/2, (l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos (pi/4)*da_i/2),), ((sin (pi/4)*da_i/2, l_absatz+radius[i],cos (pi/4)*da_i/2),), ((sin (pi/4)*da_i/2, l_absatz-2,cos (pi/4)*da_i/2),), ((sin (pi/4)*da_i/2, l_absatz-20,cos (pi/4)*da_i/2),), ((sin (pi/4)*da_i/2, 0,cos (pi/4)*da_i/2),), ((sin (pi/4)*(da_i/2-1), lw,cos (pi/4)*(da_i/2-1)),), ((sin (pi/4)*(da_i/2-1), (l_absatz+f[i]+list_lf[i]),cos (pi/4)*(da_i/2-1)),), ((sin (pi/4)*(da_i/2-1), (l_absatz+f[i]+list_lf[i]/2),cos (pi/4)*(da_i/2-1)),), ((sin (pi/4)*(da_i/2-1), (l_absatz+f[i]),cos (pi/4)*(da_i/2-1)),), ((sin (pi/4)*(da_i/2-1), (l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos (pi/4)*(da_i/2-1)),), ((sin (pi/4)*(da_i/2-1), (l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos (pi/4)*(da_i/2-1)),), ((sin (pi/4)*(da_i/2-1), l_absatz+radius[i],cos (pi/4)*(da_i/2-1)),), ((sin (pi/4)*(da_i/2-1), l_absatz-2,cos (pi/4)*(da_i/2-1)),), ((sin (pi/4)*(da_i/2-1), l_absatz-20,cos (pi/4)*(da_i/2-1)),), ((sin (pi/4)*(da_i/2-1), 0,cos (pi/4)*(da_i/2-1)),), ((sin (3*pi/4)*da_i/2, lw,cos (3*pi/4)*da_i/2),), ((sin (3*pi/4)*da_i/2, (l_absatz+f[i]+list_lf[i]),cos (3*pi/4)*da_i/2),), ((sin (3*pi/4)*da_i/2, (l_absatz+f[i]+list_lf[i]/2),cos (3*pi/4)*da_i/2),), ((sin (3*pi/4)*da_i/2, (l_absatz+f[i]),cos (3*pi/4)*da_i/2),), ((sin (3*pi/4)*da_i/2, (l_absatz+f[i]+4),cos (3*pi/4)*da_i/2),), ((sin (3*pi/4)*da_i/2, (l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos (3*pi/4)*da_i/2),), ((sin (3*pi/4)*da_i/2, (l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos (3*pi/4)*da_i/2),), ((sin (3*pi/4)*da_i/2, l_absatz+radius[i],cos (3*pi/4)*da_i/2),), ((sin (3*pi/4)*da_i/2, l_absatz-2,cos (3*pi/4)*da_i/2),), ((sin (3*pi/4)*da_i/2, l_absatz-20,cos (3*pi/4)*da_i/2),), ((sin (3*pi/4)*da_i/2, 0,cos (3*pi/4)*da_i/2),), ((sin (3*pi/4)*(da_i/2-1), lw,cos (3*pi/4)*(da_i/2-1)),), ((sin (3*pi/4)*(da_i/2-1), (l_absatz+f[i]+list_lf[i]),cos (3*pi/4)*(da_i/2-1)),), ((sin (3*pi/4)*(da_i/2-1), (l_absatz+f[i]+list_lf[i]/2),cos (3*pi/4)*(da_i/2-1)),), ((sin (3*pi/4)*(da_i/2-1), (l_absatz+f[i]),cos (3*pi/4)*(da_i/2-1)),), ((sin (3*pi/4)*(da_i/2-1), (l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos (3*pi/4)*(da_i/2-1)),), ((sin (3*pi/4)*(da_i/2-1), (l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos (3*pi/4)*(da_i/2-1)),), ((sin (3*pi/4)*(da_i/2-1), l_absatz+radius[i],cos (3*pi/4)*(da_i/2-1)),), ((sin (3*pi/4)*(da_i/2-1), l_absatz-2,cos (3*pi/4)*(da_i/2-1)),), ((sin (3*pi/4)*(da_i/2-1), l_absatz-20,cos (3*pi/4)*(da_i/2-1)),), ((sin (3*pi/4)*(da_i/2-1), 0,cos (3*pi/4)*(da_i/2-1)),), ((sin (5*pi/4)*da_i/2, lw,cos (5*pi/4)*da_i/2),), ((sin (5*pi/4)*da_i/2, (l_absatz+f[i]+list_lf[i]),cos (5*pi/4)*da_i/2),), ((sin (5*pi/4)*da_i/2, (l_absatz+f[i]+list_lf[i]/2),cos (5*pi/4)*da_i/2),), ((sin (5*pi/4)*da_i/2, (l_absatz+f[i]),cos (5*pi/4)*da_i/2),), ((sin (5*pi/4)*da_i/2, (l_absatz+f[i]+4),cos (5*pi/4)*da_i/2),), ((sin (5*pi/4)*da_i/2, (l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos (5*pi/4)*da_i/2),), ((sin (5*pi/4)*da_i/2, (l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos (5*pi/4)*da_i/2),), ((sin (5*pi/4)*da_i/2, l_absatz+radius[i],cos (5*pi/4)*da_i/2),), ((sin (5*pi/4)*da_i/2, l_absatz-2,cos (5*pi/4)*da_i/2),), ((sin (5*pi/4)*da_i/2, l_absatz-20,cos (5*pi/4)*da_i/2),), ((sin (5*pi/4)*da_i/2, 0,cos (5*pi/4)*da_i/2),), ((sin (5*pi/4)*(da_i/2-1), lw,cos (5*pi/4)*(da_i/2-1)),), ((sin (5*pi/4)*(da_i/2-1), (l_absatz+f[i]+list_lf[i]),cos (5*pi/4)*(da_i/2-1)),), ((sin (5*pi/4)*(da_i/2-1), (l_absatz+f[i]+list_lf[i]/2),cos (5*pi/4)*(da_i/2-1)),), ((sin (5*pi/4)*(da_i/2-1), (l_absatz+f[i]),cos (5*pi/4)*(da_i/2-1)),), ((sin (5*pi/4)*(da_i/2-1), (l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos (5*pi/4)*(da_i/2-1)),), ((sin (5*pi/4)*(da_i/2-1), (l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos (5*pi/4)*(da_i/2-1)),), ((sin (5*pi/4)*(da_i/2-1), l_absatz+radius[i],cos (5*pi/4)*(da_i/2-1)),), ((sin (5*pi/4)*(da_i/2-1), l_absatz-2,cos (5*pi/4)*(da_i/2-1)),), ((sin (5*pi/4)*(da_i/2-1), l_absatz-20,cos (5*pi/4)*(da_i/2-1)),), ((sin (5*pi/4)*(da_i/2-1), 0,cos (5*pi/4)*(da_i/2-1)),), ((sin (7*pi/4)*da_i/2, lw,cos (7*pi/4)*da_i/2),), ((sin (7*pi/4)*da_i/2, (l_absatz+f[i]+list_lf[i]),cos (7*pi/4)*da_i/2),), ((sin (7*pi/4)*da_i/2, (l_absatz+f[i]+list_lf[i]/2),cos (7*pi/4)*da_i/2),), ((sin (7*pi/4)*da_i/2, (l_absatz+f[i]),cos (7*pi/4)*da_i/2),), ((sin (7*pi/4)*da_i/2, (l_absatz+f[i]+4),cos (7*pi/4)*da_i/2),), ((sin (7*pi/4)*da_i/2, (l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos (7*pi/4)*da_i/2),), ((sin (7*pi/4)*da_i/2, (l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos (7*pi/4)*da_i/2),), ((sin (7*pi/4)*da_i/2, l_absatz+radius[i],cos (7*pi/4)*da_i/2),), ((sin (7*pi/4)*da_i/2, l_absatz-2,cos (7*pi/4)*da_i/2),), ((sin (7*pi/4)*da_i/2, l_absatz-20,cos (7*pi/4)*da_i/2),), ((sin (7*pi/4)*da_i/2, 0,cos (7*pi/4)*da_i/2),), ((sin (7*pi/4)*(da_i/2-1), lw,cos (7*pi/4)*(da_i/2-1)),), ((sin (7*pi/4)*(da_i/2-1), (l_absatz+f[i]+list_lf[i]),cos (7*pi/4)*(da_i/2-1)),), ((sin (7*pi/4)*(da_i/2-1), (l_absatz+f[i]+list_lf[i]/2),cos (7*pi/4)*(da_i/2-1)),), ((sin (7*pi/4)*(da_i/2-1), (l_absatz+f[i]),cos (7*pi/4)*(da_i/2-1)),), ((sin (7*pi/4)*(da_i/2-1), (l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos (7*pi/4)*(da_i/2-1)),), ((sin (7*pi/4)*(da_i/2-1), (l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos (7*pi/4)*(da_i/2-1)),), ((sin (7*pi/4)*(da_i/2-1), l_absatz+radius[i],cos (7*pi/4)*(da_i/2-1)),), ((sin (7*pi/4)*(da_i/2-1), l_absatz-2,cos (7*pi/4)*(da_i/2-1)),), ((sin (7*pi/4)*(da_i/2-1), l_absatz-20,cos (7*pi/4)*(da_i/2-1)),), ((sin (7*pi/4)*(da_i/2-1), 0,cos (7*pi/4)*(da_i/2-1)),),) welle_innen.Set(edges=Tangentialedge_Wi, name='Tangential_01') # axial Axialedge_Fr_horiz_Wi = e_wi.findAt(((sin (pi/2)*da_i/2, l_absatz+radius[i]+0.01, cos (pi/2)*da_i/2),), ((sin (pi/2)*(da_i/2-1), l_absatz+radius[i]+0.01, cos (pi/2)*(da_i/2-1)),), ((sin (pi)*da_i/2, l_absatz+radius[i]+0.01, cos (pi)*da_i/2),), ((sin (pi)*(da_i/2-1), l_absatz+radius[i]+0.01, cos (pi)*(da_i/2-1)),), ((sin (3*pi/2)*da_i/2, l_absatz+radius[i]+0.01, cos (3*pi/2)*da_i/2),), ((sin (3*pi/2)*(da_i/2-1), l_absatz+radius[i]+0.01, cos (3*pi/2)*(da_i/2-1)),), ((sin (2*pi)*da_i/2, l_absatz+radius[i]+0.01, cos (2*pi)*da_i/2),), ((sin (2*pi)*(da_i/2-1), l_absatz+radius[i]+0.01, cos (2*pi)*(da_i/2-1)),), ((0, l_absatz+radius[i]+0.01,0),),) welle_innen.Set(edges=Axialedge_Fr_horiz_Wi, name='Axial_00_FR_Horizontale') Axialedge_Fr_2R_Wi = e_wi.findAt(((sin (pi/2)*(da_i/2), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))+0.01, cos (pi/2)*(da_i/2)),), ((sin (pi/2)*(da_i/2-1), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))+0.01, cos (pi/2)*(da_i/2-1)),), ((sin (pi)*(da_i/2), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))+0.01, cos (pi)*(da_i/2)),), ((sin (pi)*(da_i/2-1), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))+0.01, cos (pi)*(da_i/2-1)),), ((sin (3*pi/2)*(da_i/2), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))+0.01, cos (3*pi/2)*(da_i/2)),), ((sin (3*pi/2)*(da_i/2-1), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))+0.01, cos (3*pi/2)*(da_i/2-1)),), ((sin (2*pi)*(da_i/2), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))+0.01, cos (2*pi)*(da_i/2)),), ((sin (2*pi)*(da_i/2-1), ((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))+0.01, cos (2*pi)*(da_i/2-1)),), ((0,((l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)))+0.01,0),),) welle_innen.Set(edges=Axialedge_Fr_2R_Wi, name='Axial_01_FR_Radius2') Axialedge_Fr_schraeg_Wi = e_wi.findAt(((sin(pi/2)*(da_i/2),l_absatz+f[i]-0.1,cos(pi/2)*(da_i/2)),), ((sin(pi/2)*(da_i/2-1),l_absatz+f[i]-0.1,cos(pi/2)*(da_i/2-1)),), ((sin(pi)*(da_i/2),l_absatz+f[i]-0.1,cos(pi)*(da_i/2)),), ((sin(pi)*(da_i/2-1),l_absatz+f[i]-0.1,cos(pi)*(da_i/2-1)),), ((sin(3*pi/2)*(da_i/2),l_absatz+f[i]-0.1,cos(3*pi/2)*(da_i/2)),), ((sin(3*pi/2)*(da_i/2-1),l_absatz+f[i]-0.1,cos(3*pi/2)*(da_i/2-1)),), ((sin(2*pi)*(da_i/2),l_absatz+f[i]-0.1,cos(2*pi)*(da_i/2)),), ((sin(2*pi)*(da_i/2-1),l_absatz+f[i]-0.1,cos(2*pi)*(da_i/2-1)),), ((0,l_absatz+f[i]-0.1,0),),) welle_innen.Set(edges=Axialedge_Fr_schraeg_Wi, name='Axial_02_FR_Schraege') Axialedge_Fr_1R_Wi = e_wi.findAt(((sin (pi/2)*(da_i/2),l_absatz,cos (pi/2)*(da_i/2)),), ((sin (pi)*(da_i/2),l_absatz,cos (pi)*(da_i/2)),), ((sin (3*pi/2)*(da_i/2),l_absatz,cos (3*pi/2)*(da_i/2)),), ((sin (2*pi)*(da_i/2),l_absatz,cos (2*pi)*(da_i/2)),), ((sin (pi/2)*(da_i/2-1),l_absatz,cos (pi/2)*(da_i/2-1)),), ((sin (pi)*(da_i/2-1),l_absatz,cos (pi)*(da_i/2-1)),), ((sin (3*pi/2)*(da_i/2-1),l_absatz,cos (3*pi/2)*(da_i/2-1)),), ((sin (2*pi)*(da_i/2-1),l_absatz,cos (2*pi)*(da_i/2-1)),), ((0,l_absatz,0),),) welle_innen.Set(edges=Axialedge_Fr_1R_Wi, name='Axial_03_FR_Radius1') Axialedge_Fugenflaeche_u_Wi = e_wi.findAt(((sin (pi/2)*(da_i/2),l_absatz+f[i]+6,cos (pi/2)*(da_i/2)),), ((sin (pi/2)*(da_i/2-1),l_absatz+f[i]+6,cos (pi/2)*(da_i/2-1)),), ((sin (pi)*(da_i/2),l_absatz+f[i]+6,cos (pi)*(da_i/2)),), ((sin (pi)*(da_i/2-1),l_absatz+f[i]+6,cos (pi)*(da_i/2-1)),), ((sin (3*pi/2)*(da_i/2),l_absatz+f[i]+6,cos (3*pi/2)*(da_i/2)),), ((sin (3*pi/2)*(da_i/2-1),l_absatz+f[i]+6,cos (3*pi/2)*(da_i/2-1)),), ((sin (2*pi)*(da_i/2),l_absatz+f[i]+6,cos (2*pi)*(da_i/2)),), ((sin (2*pi)*(da_i/2-1),l_absatz+f[i]+6,cos (2*pi)*(da_i/2-1)),), ((0,l_absatz+f[i]+6,0),),) welle_innen.Set(edges=Axialedge_Fugenflaeche_u_Wi, name='Axial_05_Fugenflaeche_unten') Axialedge_Fugenflaeche_o_Wi = e_wi.findAt(((sin (pi/2)*(da_i/2),l_absatz+f[i]+list_lf[i]/2+1,cos (pi/2)*(da_i/2)),), ((sin (pi/2)*(da_i/2-1),l_absatz+f[i]+list_lf[i]/2+1,cos (pi/2)*(da_i/2-1)),), ((sin (pi)*(da_i/2),l_absatz+f[i]+list_lf[i]/2+1,cos (pi)*(da_i/2)),), ((sin (pi)*(da_i/2-1),l_absatz+f[i]+list_lf[i]/2+1,cos (pi)*(da_i/2-1)),), ((sin (3*pi/2)*(da_i/2),l_absatz+f[i]+list_lf[i]/2+1,cos (3*pi/2)*(da_i/2)),), ((sin (3*pi/2)*(da_i/2-1),l_absatz+f[i]+list_lf[i]/2+1,cos (3*pi/2)*(da_i/2-1)),), ((sin (2*pi)*(da_i/2),l_absatz+f[i]+list_lf[i]/2+1,cos (2*pi)*(da_i/2)),), ((sin (2*pi)*(da_i/2-1),l_absatz+f[i]+list_lf[i]/2+1,cos (2*pi)*(da_i/2-1)),), ((0,l_absatz+f[i]+list_lf[i]/2+1,0),),) welle_innen.Set(edges=Axialedge_Fugenflaeche_o_Wi, name='Axial_06_Fugenflaeche_oben') Axialedge_Kante_vor_Absatz_Wi = e_wi.findAt(((sin (pi/2)*(da_i/2),l_absatz-4,cos (pi/2)*(da_i/2)),), ((sin (pi/2)*(da_i/2-1),l_absatz-4,cos (pi/2)*(da_i/2-1)),), ((sin (pi)*(da_i/2),l_absatz-4,cos (pi)*(da_i/2)),), ((sin (pi)*(da_i/2-1),l_absatz-4,cos (pi)*(da_i/2-1)),), ((sin (3*pi/2)*(da_i/2),l_absatz-4,cos (3*pi/2)*(da_i/2)),), ((sin (3*pi/2)*(da_i/2-1),l_absatz-4,cos (3*pi/2)*(da_i/2-1)),), ((sin (2*pi)*(da_i/2),l_absatz-4,cos (2*pi)*(da_i/2)),), ((sin (2*pi)*(da_i/2-1),l_absatz-4,cos (2*pi)*(da_i/2-1)),), ((0,l_absatz-4,0),),) welle_innen.Set(edges=Axialedge_Kante_vor_Absatz_Wi, name='Axial_06_Kante_vor_Absatz') Axialedge_Auslauf_Wi = e_wi.findAt(((sin (pi/2)*(da_i/2),lw-2,cos (pi/2)*(da_i/2)),), ((sin (pi/2)*(da_i/2-1),lw-2,cos (pi/2)*(da_i/2-1)),), ((sin (pi)*(da_i/2),lw-2,cos (pi)*(da_i/2)),), ((sin (pi)*(da_i/2-1),lw-2,cos (pi)*(da_i/2-1)),), ((sin (3*pi/2)*(da_i/2),lw-2,cos (3*pi/2)*(da_i/2)),), ((sin (3*pi/2)*(da_i/2-1),lw-2,cos (3*pi/2)*(da_i/2-1)),), ((sin (2*pi)*(da_i/2),lw-2,cos (2*pi)*(da_i/2)),), ((sin (2*pi)*(da_i/2-1),lw-2,cos (2*pi)*(da_i/2-1)),), ((0,lw-2,0),),) welle_innen.Set(edges=Axialedge_Auslauf_Wi, name='Axial_07_Auslauf') Axialedge_Absatz_Wi = e_wi.findAt(((sin (pi/2)*(da_i/2),4,cos (pi/2)*(da_i/2)),), ((sin (pi/2)*(da_i/2-1),4,cos (pi/2)*(da_i/2-1)),), ((sin (pi)*(da_i/2),4,cos (pi)*(da_i/2)),), ((sin (pi)*(da_i/2-1),4,cos (pi)*(da_i/2-1)),), ((sin (3*pi/2)*(da_i/2),4,cos (3*pi/2)*(da_i/2)),), ((sin (3*pi/2)*(da_i/2-1),4,cos (3*pi/2)*(da_i/2-1)),), ((sin (2*pi)*(da_i/2),4,cos (2*pi)*(da_i/2)),), ((sin (2*pi)*(da_i/2-1),4,cos (2*pi)*(da_i/2-1)),), ((0,4,0),),) welle_innen.Set(edges=Axialedge_Absatz_Wi, name='Axial_08_Absatz') Axialedge_Fugenflaeche_Anfang_Wi = e_wi.findAt(((sin (pi/2)*(da_i/2),l_absatz+f[i]+3,cos (pi/2)*(da_i/2)),), ((sin (pi/2)*(da_i/2-1),l_absatz+f[i]+3,cos (pi/2)*(da_i/2-1)),), ((sin (pi)*(da_i/2),l_absatz+f[i]+3,cos (pi)*(da_i/2)),), ((sin (pi)*(da_i/2-1),l_absatz+f[i]+3,cos (pi)*(da_i/2-1)),), ((sin (3*pi/2)*(da_i/2),l_absatz+f[i]+3,cos (3*pi/2)*(da_i/2)),), ((sin (3*pi/2)*(da_i/2-1),l_absatz+f[i]+3,cos (3*pi/2)*(da_i/2-1)),), ((sin (2*pi)*(da_i/2),l_absatz+f[i]+3,cos (2*pi)*(da_i/2)),), ((sin (2*pi)*(da_i/2-1),l_absatz+f[i]+3,cos (2*pi)*(da_i/2-1)),), ((0,l_absatz+f[i]+3,0),),) welle_innen.Set(edges=Axialedge_Fugenflaeche_Anfang_Wi, name='Axial_09_Fugenflaeche_Anfang') #radial Radialedge_FR_Fugenflaeche_a_Wi = e_wi.findAt(((sin(pi/2)*(da_i/2-0.5),l_absatz+radius[i],cos(pi/2)*(da_i/2-0.5)),), ((sin(pi/2)*(da_i/2-0.5),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos(pi/2)*(da_i/2-0.5)),), ((sin(pi/2)*(da_i/2-0.5),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos(pi/2)*(da_i/2-0.5)),), ((sin(pi/2)*(da_i/2-0.5),l_absatz+f[i],cos(pi/2)*(da_i/2-0.5)),), ((sin(pi/2)*(da_i/2-0.5),l_absatz+f[i]+4,cos(pi/2)*(da_i/2-0.5)),), ((sin(pi/2)*(da_i/2-0.5),(l_absatz+f[i]+list_lf[i]),cos(pi/2)*(da_i/2-0.5)),), ((sin(pi/2)*(da_i/2-0.5),(l_absatz+f[i]+list_lf[i]/2),cos(pi/2)*(da_i/2-0.5)),), ((sin(pi)*(da_i/2-0.5),l_absatz+radius[i],cos(pi)*(da_i/2-0.5)),), ((sin(pi)*(da_i/2-0.5),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos(pi)*(da_i/2-0.5)),), ((sin(pi)*(da_i/2-0.5),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos(pi)*(da_i/2-0.5)),), ((sin(pi)*(da_i/2-0.5),l_absatz+f[i],cos(pi)*(da_i/2-0.5)),), ((sin(pi)*(da_i/2-0.5),l_absatz+f[i]+4,cos(pi)*(da_i/2-0.5)),), ((sin(pi)*(da_i/2-0.5),(l_absatz+f[i]+list_lf[i]),cos(pi)*(da_i/2-0.5)),), ((sin(pi)*(da_i/2-0.5),(l_absatz+f[i]+list_lf[i]/2),cos(pi)*(da_i/2-0.5)),), ((sin(3*pi/2)*(da_i/2-0.5),l_absatz+radius[i],cos(3*pi/2)*(da_i/2-0.5)),), ((sin(3*pi/2)*(da_i/2-0.5),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos(3*pi/2)*(da_i/2-0.5)),), ((sin(3*pi/2)*(da_i/2-0.5),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos(3*pi/2)*(da_i/2-0.5)),), ((sin(3*pi/2)*(da_i/2-0.5),l_absatz+f[i],cos(3*pi/2)*(da_i/2-0.5)),), ((sin(3*pi/2)*(da_i/2-0.5),l_absatz+f[i]+4,cos(3*pi/2)*(da_i/2-0.5)),), ((sin(3*pi/2)*(da_i/2-0.5),(l_absatz+f[i]+list_lf[i]),cos(3*pi/2)*(da_i/2-0.5)),), ((sin(3*pi/2)*(da_i/2-0.5),(l_absatz+f[i]+list_lf[i]/2),cos(3*pi/2)*(da_i/2-0.5)),), ((sin(2*pi)*(da_i/2-0.5),l_absatz+radius[i],cos(2*pi)*(da_i/2-0.5)),), ((sin(2*pi)*(da_i/2-0.5),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos(2*pi)*(da_i/2-0.5)),), ((sin(2*pi)*(da_i/2-0.5),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos(2*pi)*(da_i/2-0.5)),), ((sin(2*pi)*(da_i/2-0.5),l_absatz+f[i],cos(2*pi)*(da_i/2-0.5)),), ((sin(2*pi)*(da_i/2-0.5),l_absatz+f[i]+4,cos(2*pi)*(da_i/2-0.5)),), ((sin(2*pi)*(da_i/2-0.5),(l_absatz+f[i]+list_lf[i]),cos(2*pi)*(da_i/2-0.5)),), ((sin(2*pi)*(da_i/2-0.5),(l_absatz+f[i]+list_lf[i]),cos(2*pi)*(da_i/2-0.5)),),) welle_innen.Set(edges=Radialedge_FR_Fugenflaeche_a_Wi, name='Radial_00_Freistich+Fugenflaeche_aussen') Radialedge_FR_Fugenflaeche_i_Wi = e_wi.findAt(((sin(pi/2)*(da_i/2-3),l_absatz+radius[i],cos(pi/2)*(da_i/2-3)),), ((sin(pi/2)*(da_i/2-3),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos(pi/2)*(da_i/2-3)),), ((sin(pi/2)*(da_i/2-3),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos(pi/2)*(da_i/2-3)),), ((sin(pi/2)*(da_i/2-3),l_absatz+f[i],cos(pi/2)*(da_i/2-3)),), ((sin(pi/2)*(da_i/2-3),l_absatz+f[i]+4,cos(pi/2)*(da_i/2-3)),), ((sin(pi/2)*(da_i/2-3),(l_absatz+f[i]+list_lf[i]),cos(pi/2)*(da_i/2-3)),), ((sin(pi/2)*(da_i/2-3),(l_absatz+f[i]+list_lf[i]/2),cos(pi/2)*(da_i/2-3)),), ((sin(pi/2)*(da_i/2-3),lw,cos(pi/2)*(da_i/2-3)),), ((sin(pi)*(da_i/2-3),l_absatz+radius[i],cos(pi)*(da_i/2-3)),), ((sin(pi)*(da_i/2-3),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos(pi)*(da_i/2-3)),), ((sin(pi)*(da_i/2-3),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos(pi)*(da_i/2-3)),), ((sin(pi)*(da_i/2-3),l_absatz+f[i],cos(pi)*(da_i/2-3)),), ((sin(pi)*(da_i/2-3),l_absatz+f[i]+4,cos(pi)*(da_i/2-3)),), ((sin(pi)*(da_i/2-3),(l_absatz+f[i]+list_lf[i]),cos(pi)*(da_i/2-3)),), ((sin(pi)*(da_i/2-3),(l_absatz+f[i]+list_lf[i]/2),cos(pi)*(da_i/2-3)),), ((sin(pi)*(da_i/2-3),lw,cos(pi)*(da_i/2-3)),), ((sin(3*pi/2)*(da_i/2-3),l_absatz+radius[i],cos(3*pi/2)*(da_i/2-3)),), ((sin(3*pi/2)*(da_i/2-3),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos(3*pi/2)*(da_i/2-3)),), ((sin(3*pi/2)*(da_i/2-3),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos(3*pi/2)*(da_i/2-3)),), ((sin(3*pi/2)*(da_i/2-3),l_absatz+f[i],cos(3*pi/2)*(da_i/2-3)),), ((sin(3*pi/2)*(da_i/2-3),l_absatz+f[i]+4,cos(3*pi/2)*(da_i/2-3)),), ((sin(3*pi/2)*(da_i/2-3),(l_absatz+f[i]+list_lf[i]),cos(3*pi/2)*(da_i/2-3)),), ((sin(3*pi/2)*(da_i/2-3),(l_absatz+f[i]+list_lf[i]/2),cos(3*pi/2)*(da_i/2-3)),), ((sin(3*pi/2)*(da_i/2-3),lw,cos(3*pi/2)*(da_i/2-3)),), ((sin(2*pi)*(da_i/2-3),l_absatz+radius[i],cos(2*pi)*(da_i/2-3)),), ((sin(2*pi)*(da_i/2-3),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos(2*pi)*(da_i/2-3)),), ((sin(2*pi)*(da_i/2-3),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos(2*pi)*(da_i/2-3)),), ((sin(2*pi)*(da_i/2-3),l_absatz+f[i],cos(2*pi)*(da_i/2-3)),), ((sin(2*pi)*(da_i/2-3),l_absatz+f[i]+4,cos(2*pi)*(da_i/2-3)),), ((sin(2*pi)*(da_i/2-3),(l_absatz+f[i]+list_lf[i]),cos(2*pi)*(da_i/2-3)),), ((sin(2*pi)*(da_i/2-3),(l_absatz+f[i]+list_lf[i]/2),cos(2*pi)*(da_i/2-3)),), ((sin(2*pi)*(da_i/2-3),lw,cos(2*pi)*(da_i/2-3)),),) welle_innen.Set(edges=Radialedge_FR_Fugenflaeche_i_Wi, name='Radial_01_Freistich+Fugenflaeche_innen') Radialedge_Anfangsteil_Wi = e_wi.findAt(((sin(pi/2)*(da_i/2-0.5),0,cos(pi/2)*(da_i/2-0.5)),), ((sin(pi/2)*(da_i/2-3),0,cos(pi/2)*(da_i/2-3)),), ((sin(pi)*(da_i/2-0.5),0,cos(pi)*(da_i/2-0.5)),), ((sin(pi)*(da_i/2-3),0,cos(pi)*(da_i/2-3)),), ((sin(3*pi/2)*(da_i/2-0.5),0,cos(3*pi/2)*(da_i/2-0.5)),), ((sin(3*pi/2)*(da_i/2-3),0,cos(3*pi/2)*(da_i/2-3)),), ((sin(2*pi)*(da_i/2-0.5),0,cos(2*pi)*(da_i/2-0.5)),), ((sin(2*pi)*(da_i/2-3),0,cos(2*pi)*(da_i/2-3)),),) welle_innen.Set(edges=Radialedge_Anfangsteil_Wi, name='Radial_02_Anfangsteil') Radialedge_Rundung_Wi = e_wi.findAt(((sin(pi/2)*(da_i/2-0.01),l_absatz-2,cos(pi/2)*(da_i/2-0.01)),), ((sin(pi/2)*(da_i/2-0.01),l_absatz-20,cos(pi/2)*(da_i/2-0.01)),), ((sin(pi)*(da_i/2-0.01),l_absatz-2,cos(pi)*(da_i/2-0.01)),), ((sin(pi)*(da_i/2-0.01),l_absatz-20,cos(pi)*(da_i/2-0.01)),), ((sin(3*pi/2)*(da_i/2-0.01),l_absatz-2,cos(3*pi/2)*(da_i/2-0.01)),), ((sin(3*pi/2)*(da_i/2-0.01),l_absatz-20,cos(3*pi/2)*(da_i/2-0.01)),), ((sin(2*pi)*(da_i/2-0.01),l_absatz-2,cos(2*pi)*(da_i/2-0.01)),), ((sin(2*pi)*(da_i/2-0.01),l_absatz-20,cos(2*pi)*(da_i/2-0.01)),),) welle_innen.Set(edges=Radialedge_Rundung_Wi, name='Radial_03_Rundung') Radialedge_Auslauf_Wi = e_wi.findAt(((sin(pi/2)*(da_i/2-0.01),lw,cos(pi/2)*(da_i/2-0.01)),), ((sin(pi)*(da_i/2-0.01),lw,cos(pi)*(da_i/2-0.01)),), ((sin(3*pi/2)*(da_i/2-0.01),lw,cos(3*pi/2)*(da_i/2-0.01)),), ((sin(2*pi)*(da_i/2-0.01),lw,cos(2*pi)*(da_i/2-0.01)),),) welle_innen.Set(edges=Radialedge_Auslauf_Wi, name='Radial_04_Auslauf') # Interaction #RP erzeugen # Referenzpunkt definieren wnv = mdb.models[Modell_Name[i]].rootAssembly wnv.ReferencePoint(point=(0.0, -1000.0+l_absatz+1.5, 0.0)) # Erzeugt Referenzpunkt mdb.models[Modell_Name[i]].rootAssembly.features.changeKey( fromName='RP-1', toName='RP-Welle') # Umbenennung wnv.ReferencePoint(point=(0.0, l_absatz+f[i]+list_lf[i], 0.0)) mdb.models[Modell_Name[i]].rootAssembly.features.changeKey( fromName='RP-1', toName='RP-Nabe') if i > 8 and i < 13: wnv.ReferencePoint(point=(0.0, 0.0, 0.0)) mdb.models[Modell_Name[i]].rootAssembly.features.changeKey( fromName='RP-1', toName='RP-Welle_innen') #Set RP-Welle wnv = mdb.models[Modell_Name[i]].rootAssembly r1_w = wnv.referencePoints refPoints1=(wnv.referencePoints.findAt((0.0, -1000.0+l_absatz+1.5, 0.0)), ) wnv.Set(referencePoints=refPoints1, name='RP-Welle') #Set RP-Nabe wnv = mdb.models[Modell_Name[i]].rootAssembly r1_n = wnv.referencePoints refPoints2=(wnv.referencePoints.findAt((0.0, l_absatz+f[i]+list_lf[i], 0.0)), ) wnv.Set(referencePoints=refPoints2, name='RP-Nabe') if i > 8 and i < 13: #Set RP-Welle innen wnv = mdb.models[Modell_Name[i]].rootAssembly r1_wi = wnv.referencePoints refPoints3=(wnv.referencePoints.findAt((0.0, 0.0, 0.0)), ) wnv.Set(referencePoints=refPoints3, name='RP-Welle_innen') #Zwangsbedingungen definieren # Welle wnv = mdb.models[Modell_Name[i]].rootAssembly region1=wnv.sets['RP-Welle'] wnv = mdb.models[Modell_Name[i]].rootAssembly region2=wnv.instances[Wellen_Name[i]].sets['Einspannung Welle'] mdb.models[Modell_Name[i]].Coupling(name='Welleneinspannung', controlPoint=region1, surface=region2, influenceRadius=WHOLE_SURFACE, couplingType=KINEMATIC, localCsys=None, u1=OFF, u2=ON, u3=OFF, ur1=ON, ur2=ON, ur3=ON) #Zwangsbedingungen definieren # Nabe wnv = mdb.models[Modell_Name[i]].rootAssembly region1=wnv.sets['RP-Nabe'] wnv = mdb.models[Modell_Name[i]].rootAssembly region2=wnv.instances[Naben_Name[i]].sets['Einspannung Nabe'] mdb.models[Modell_Name[i]].Coupling(name='Nabeneinspannung', controlPoint=region1, surface=region2, influenceRadius=WHOLE_SURFACE, couplingType=KINEMATIC, localCsys=None, u1=OFF, u2=ON, u3=OFF, ur1=OFF, ur2=ON, ur3=OFF) if i > 8 and i < 13: #Zwangsbedingungen definieren # innenliegende Welle wnv = mdb.models[Modell_Name[i]].rootAssembly region1=wnv.sets['RP-Welle_innen'] wnv = mdb.models[Modell_Name[i]].rootAssembly region2=wnv.instances[Welle_innen_Name[i]].sets['Einspannung Welle innen'] mdb.models[Modell_Name[i]].Coupling(name='Wellen_innen_Einspannung', controlPoint=region1, surface=region2, influenceRadius=WHOLE_SURFACE, couplingType=KINEMATIC, localCsys=None, u1=OFF, u2=ON, u3=OFF, ur1=ON, ur2=ON, ur3=ON) #Wechselwirkung definieren # Normales Verhalten mdb.models[Modell_Name[i]].ContactProperty('Normal') mdb.models[Modell_Name[i]].interactionProperties['Normal'].NormalBehavior( pressureOverclosure=HARD, allowSeparation=ON, contactStiffness=10000000.0, contactStiffnessScaleFactor=1.0, clearanceAtZeroContactPressure=0.0, stiffnessBehavior=LINEAR, constraintEnforcementMethod=PENALTY) # Wechselwirkung definieren # Normales & Tangentiales Verhalten mdb.models[Modell_Name[i]].ContactProperty( 'Normal_Tangential (my 0,2)') mdb.models[Modell_Name[i]].interactionProperties['Normal_Tangential (my 0,2)'].NormalBehavior( pressureOverclosure=HARD, allowSeparation=ON, contactStiffness=10000000.0, contactStiffnessScaleFactor=1.0, clearanceAtZeroContactPressure=0.0, stiffnessBehavior=LINEAR, constraintEnforcementMethod=PENALTY) mdb.models[Modell_Name[i]].interactionProperties['Normal_Tangential (my 0,2)'].TangentialBehavior( formulation=PENALTY, directionality=ISOTROPIC, slipRateDependency=OFF, pressureDependency=OFF, temperatureDependency=OFF, dependencies=0, table=((reibwert, ), ), shearStressLimit=None, maximumElasticSlip=ABSOLUTE_DISTANCE, absoluteDistance=0.0001, elasticSlipStiffness=None) # Wechselwirkung der Kontaktflächen definieren wnv = mdb.models[Modell_Name[i]].rootAssembly session.viewports['Viewport: 1'].setValues(displayedObject=wnv) session.viewports['Viewport: 1'].assemblyDisplay.setValues( step='Querpressverbund') wnv = mdb.models[Modell_Name[i]].rootAssembly region1=wnv.instances[Naben_Name[i]].surfaces['Kontaktflaeche PV'] wnv = mdb.models[Modell_Name[i]].rootAssembly region2=wnv.instances[Wellen_Name[i]].surfaces['Kontaktflaeche PV'] mdb.models[Modell_Name[i]].SurfaceToSurfaceContactStd(name='PV', createStepName='Querpressverbund', master=region1, slave=region2, sliding=FINITE, thickness=ON, interactionProperty='Normal', interferenceType=UNIFORM, overclosure=-uebermass[i]/2, interferenceDirectionType=COMPUTED, amplitude=None, adjustMethod=NONE, initialClearance=OMIT, datumAxis=None, clearanceRegion=None) # Wechselwirkung der Kontaktflächen modifizieren wnv = mdb.models[Modell_Name[i]].rootAssembly session.viewports['Viewport: 1'].setValues(displayedObject=wnv) session.viewports['Viewport: 1'].assemblyDisplay.setValues(step='Reibung') mdb.models[Modell_Name[i]].TabularAmplitude(name='Const', timeSpan=STEP, smooth=SOLVER_DEFAULT, data=((0.0, 1.0), (1.0, 1.0))) mdb.models[Modell_Name[i]].interactions['PV'].setValuesInStep( stepName='Reibung', interactionProperty='Normal_Tangential (my 0,2)') mdb.models[Modell_Name[i]].interactions['PV'].setValuesInStep( stepName='Reibung', interferenceType=UNIFORM, overclosure=-uebermass[i]/2, interferenceDirectionType=COMPUTED, amplitude='Const') if i > 8 and i < 13: # Wechselwirkung der Kontaktflächen innen definieren wnv = mdb.models[Modell_Name[i]].rootAssembly session.viewports['Viewport: 1'].setValues(displayedObject=wnv) session.viewports['Viewport: 1'].assemblyDisplay.setValues( step='Querpressverbund') wnv = mdb.models[Modell_Name[i]].rootAssembly region1=wnv.instances[Wellen_Name[i]].surfaces['Kontaktflaeche Innen'] wnv = mdb.models[Modell_Name[i]].rootAssembly region2=wnv.instances[Welle_innen_Name[i]].surfaces['Kontaktflaeche'] mdb.models[Modell_Name[i]].SurfaceToSurfaceContactStd(name='PV innen', createStepName='Querpressverbund', master=region1, slave=region2, sliding=FINITE, thickness=ON, interactionProperty='Normal', interferenceType=UNIFORM, overclosure=-uebermass_wi[i]/2, interferenceDirectionType=COMPUTED, amplitude=None, adjustMethod=NONE, initialClearance=OMIT, datumAxis=None, clearanceRegion=None) # Wechselwirkung der Kontaktflächen innen modifizieren wnv = mdb.models[Modell_Name[i]].rootAssembly session.viewports['Viewport: 1'].setValues(displayedObject=wnv) session.viewports['Viewport: 1'].assemblyDisplay.setValues(step='Reibung') mdb.models[Modell_Name[i]].interactions['PV innen'].setValuesInStep( stepName='Reibung', interactionProperty='Normal_Tangential (my 0,2)') mdb.models[Modell_Name[i]].interactions['PV innen'].setValuesInStep( stepName='Reibung', interferenceType=UNIFORM, overclosure=-uebermass_wi[i]/2, interferenceDirectionType=COMPUTED, amplitude='Const') # Randbedingung definieren # Feste Einspannung # Welle session.viewports['Viewport: 1'].assemblyDisplay.setValues(step='Initial') wnv = mdb.models[Modell_Name[i]].rootAssembly region = wnv.sets['RP-Welle'] mdb.models[Modell_Name[i]].DisplacementBC(name='Einspannung Welle', createStepName='Initial', region=region, u1=SET, u2=SET, u3=SET, ur1=SET, ur2=SET, ur3=SET, amplitude=UNSET, distributionType=UNIFORM, fieldName='', localCsys=None) # Randbedingung definieren # Feste Einspannung # Nabe wnv = mdb.models[Modell_Name[i]].rootAssembly region = wnv.sets['RP-Nabe'] mdb.models[Modell_Name[i]].EncastreBC(name='Einspannung Nabe', createStepName='Initial', region=region, localCsys=None) if i > 8 and i < 13: # Randbedingung definieren # Feste Einspannung # Welle session.viewports['Viewport: 1'].assemblyDisplay.setValues(step='Initial') wnv = mdb.models[Modell_Name[i]].rootAssembly region = wnv.sets['RP-Welle_innen'] mdb.models[Modell_Name[i]].DisplacementBC(name='Einspannung Welle innen', createStepName='Initial', region=region, u1=SET, u2=SET, u3=SET, ur1=SET, ur2=SET, ur3=SET, amplitude=UNSET, distributionType=UNIFORM, fieldName='', localCsys=None) # Netzerstellung # Welle allgemein welle = mdb.models[Modell_Name[i]].parts[Wellen_Name[i]] session.viewports['Viewport: 1'].setValues(displayedObject=welle) session.viewports['Viewport: 1'].partDisplay.setValues(mesh=ON) session.viewports['Viewport: 1'].partDisplay.meshOptions.setValues( meshTechnique=ON) session.viewports['Viewport: 1'].partDisplay.geometryOptions.setValues( referenceRepresentation=OFF) #tangential welle = mdb.models[Modell_Name[i]].parts[Wellen_Name[i]] e_w = welle.edges welle.seedEdgeByNumber(edges=Tangentialedge_Welle_Absatz, number=FU, constraint=FIXED) welle.seedEdgeByNumber(edges=Tangentialedge_Welle, number=FU, constraint=FIXED) #axial welle.seedEdgeByNumber(edges=Axialedge_Fr_horiz, number=DR, constraint=FIXED) if list_da[i] > 120: for t in range (5): welle.seedEdgeByNumber(edges=Axialedge_Fr_2R[t], number=DR, constraint=FIXED) else: welle.seedEdgeByNumber(edges=Axialedge_Fr_2R, number=DR, constraint=FIXED) welle.seedEdgeByNumber(edges=Axialedge_Fr_schraeg, number=VI, constraint=FIXED) Edge1 = e_w.findAt(((sin (pi/2)*(list_da[i]/2-1),l_absatz,cos (pi/2)*(list_da[i]/2-1)),), ((sin (pi)*(list_da[i]/2-1),l_absatz,cos (pi)*(list_da[i]/2-1)),), ((sin (pi)*(list_di[i]/2),l_absatz,cos (pi)*(list_di[i]/2)),),) Edge2 = e_w.findAt(((sin (pi/2)*(list_di[i]/2),l_absatz,cos (pi/2)*(list_di[i]/2)),), ((sin (3*pi/2)*(list_da[i]/2-1),l_absatz,cos (3*pi/2)*(list_da[i]/2-1)),), ((sin (3*pi/2)*(list_di[i]/2),l_absatz,cos (3*pi/2)*(list_di[i]/2)),), ((sin (2*pi)*(list_da[i]/2-1),l_absatz,cos (2*pi)*(list_da[i]/2-1)),), ((sin (2*pi)*(list_di[i]/2),l_absatz,cos (2*pi)*(list_di[i]/2)),),) welle.seedEdgeByBias(biasMethod=SINGLE, end1Edges=Edge1, end2Edges=Edge2, ratio=9.0, number=FU, constraint=FINER) welle.seedEdgeByNumber(edges=Axialedge_Fugenflaeche_o, number=VI, constraint=FIXED) welle.seedEdgeByNumber(edges=Axialedge_Fugenflaeche_u, number=VI, constraint=FIXED) Edge1 = e_w.findAt(((sin (pi/2)*(list_da[i]/2-1),l_absatz-4,cos (pi/2)*(list_da[i]/2-1)),), ((sin (pi)*(list_da[i]/2-1),l_absatz-4,cos (pi)*(list_da[i]/2-1)),), ((sin (pi)*(list_di[i]/2),l_absatz-4,cos (pi)*(list_di[i]/2)),), ((sin (3*pi/2)*(list_da[i]/2+radius[i]),l_absatz-4,cos (3*pi/2)*(list_da[i]/2+radius[i])),), ((sin (2*pi)*(list_da[i]/2+radius[i]),l_absatz-4,cos (2*pi)*(list_da[i]/2+radius[i])),), ((sin(pi/2)*(list_da[i]/2+h_absatz[i]),l_absatz-4,cos(pi/2)*(list_da[i]/2+h_absatz[i])),), ((sin(3*pi/2)*(list_da[i]/2+h_absatz[i]),l_absatz-4,cos(3*pi/2)*(list_da[i]/2+h_absatz[i])),), ((sin(2*pi)*(list_da[i]/2+h_absatz[i]),l_absatz-4,cos(2*pi)*(list_da[i]/2+h_absatz[i])),),) Edge2 = e_w.findAt(((sin (pi/2)*(list_da[i]/2+radius[i]),l_absatz-4,cos (pi/2)*(list_da[i]/2+radius[i])),), ((sin (pi/2)*(list_di[i]/2),l_absatz-4,cos (pi/2)*(list_di[i]/2)),), ((sin(pi)*(list_da[i]/2+h_absatz[i]),l_absatz-4,cos(pi)*(list_da[i]/2+h_absatz[i])),), ((sin (pi)*(list_da[i]/2+radius[i]),l_absatz-4,cos (pi)*(list_da[i]/2+radius[i])),), ((sin (3*pi/2)*(list_da[i]/2-1),l_absatz-4,cos (3*pi/2)*(list_da[i]/2-1)),), ((sin (3*pi/2)*(list_di[i]/2),l_absatz-4,cos (3*pi/2)*(list_di[i]/2)),), ((sin (2*pi)*(list_da[i]/2-1),l_absatz-4,cos (2*pi)*(list_da[i]/2-1)),), ((sin (2*pi)*(list_di[i]/2),l_absatz-4,cos (2*pi)*(list_di[i]/2)),),) welle.seedEdgeByBias(biasMethod=SINGLE, end1Edges=Edge1, end2Edges=Edge2, ratio=4.0, number=FU, constraint=FINER) welle.seedEdgeByNumber(edges=Axialedge_Auslauf, number=EI, constraint=FIXED) welle.seedEdgeByNumber(edges=Axialedge_Absatz, number=EI, constraint=FIXED) Edge1 = e_w.findAt(((sin(pi/2)*(list_da[i]/2+h_absatz[i]),l_absatz-0.1,cos(pi/2)*(list_da[i]/2+h_absatz[i])),), ((sin(3*pi/2)*(list_da[i]/2+h_absatz[i]),l_absatz-0.1,cos(3*pi/2)*(list_da[i]/2+h_absatz[i])),), ((sin(2*pi)*(list_da[i]/2+h_absatz[i]),l_absatz-0.1,cos(2*pi)*(list_da[i]/2+h_absatz[i])),), ((sin(3*pi/2)*(list_da[i]/2+radius[i]),l_absatz-0.1,cos(3*pi/2)*(list_da[i]/2+radius[i])),), ((sin(2*pi)*(list_da[i]/2+radius[i]),l_absatz-0.1,cos(2*pi)*(list_da[i]/2+radius[i])),),) Edge2 = e_w.findAt(((sin(pi/2)*(list_da[i]/2+radius[i]),l_absatz-0.1,cos(pi/2)*(list_da[i]/2+radius[i])),), ((sin(pi)*(list_da[i]/2+h_absatz[i]),l_absatz-0.1,cos(pi)*(list_da[i]/2+h_absatz[i])),), ((sin(pi)*(list_da[i]/2+radius[i]),l_absatz-0.1,cos(pi)*(list_da[i]/2+radius[i])),),) welle.seedEdgeByBias(biasMethod=SINGLE, end1Edges=Edge1, end2Edges=Edge2, ratio=4.0, number=ZW, constraint=FINER) Edge1 = e_w.findAt(((sin (pi/2)*(list_di[i]/2),l_absatz+f[i]+3,cos (pi/2)*(list_di[i]/2)),), ((sin (pi)*(list_da[i]/2),l_absatz+f[i]+3,cos (pi)*(list_da[i]/2)),), ((sin (3*pi/2)*(list_da[i]/2-1),l_absatz+f[i]+3,cos (3*pi/2)*(list_da[i]/2-1)),), ((sin (3*pi/2)*(list_di[i]/2),l_absatz+f[i]+3,cos (3*pi/2)*(list_di[i]/2)),), ((sin (2*pi)*(list_da[i]/2-1),l_absatz+f[i]+3,cos (2*pi)*(list_da[i]/2-1)),), ((sin (2*pi)*(list_di[i]/2),l_absatz+f[i]+3,cos (2*pi)*(list_di[i]/2)),),) Edge2 = e_w.findAt(((sin (pi/2)*(list_da[i]/2),l_absatz+f[i]+3,cos (pi/2)*(list_da[i]/2)),), ((sin (pi/2)*(list_da[i]/2-1),l_absatz+f[i]+3,cos (pi/2)*(list_da[i]/2-1)),), ((sin (pi)*(list_da[i]/2-1),l_absatz+f[i]+3,cos (pi)*(list_da[i]/2-1)),), ((sin (pi)*(list_di[i]/2),l_absatz+f[i]+3,cos (pi)*(list_di[i]/2)),), ((sin (3*pi/2)*(list_da[i]/2),l_absatz+f[i]+3,cos (3*pi/2)*(list_da[i]/2)),), ((sin (2*pi)*(list_da[i]/2),l_absatz+f[i]+3,cos (2*pi)*(list_da[i]/2)),),) welle.seedEdgeByBias(biasMethod=SINGLE, end1Edges=Edge1, end2Edges=Edge2, ratio=4.0, number=SE, constraint=FINER) #radial Edge1 = e_w.findAt(((sin(3*pi/2)*(list_da[i]/2-1+0.1),l_absatz+radius[i],cos(3*pi/2)*(list_da[i]/2-1+0.1)),), ((sin(3*pi/2)*(list_da[i]/2-1+0.1),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos(3*pi/2)*(list_da[i]/2-1+0.1)),), ((sin(3*pi/2)*(list_da[i]/2-1+0.1),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos(3*pi/2)*(list_da[i]/2-1+0.1)),), ((sin(3*pi/2)*(list_da[i]/2-1+0.1),l_absatz+f[i],cos(3*pi/2)*(list_da[i]/2-1+0.1)),), ((sin(3*pi/2)*(list_da[i]/2-1+0.1),l_absatz+f[i]+4,cos(3*pi/2)*(list_da[i]/2-1+0.1)),), ((sin(3*pi/2)*(list_da[i]/2-1+0.1),(l_absatz+f[i]+list_lf[i]),cos(3*pi/2)*(list_da[i]/2-1+0.1)),), ((sin(3*pi/2)*(list_da[i]/2-1+0.1),(l_absatz+f[i]+list_lf[i]/2),cos(3*pi/2)*(list_da[i]/2-1+0.1)),), ((sin(2*pi)*(list_da[i]/2-1+0.1),l_absatz+radius[i],cos(2*pi)*(list_da[i]/2-1+0.1)),), ((sin(2*pi)*(list_da[i]/2-1+0.1),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos(2*pi)*(list_da[i]/2-1+0.1)),), ((sin(2*pi)*(list_da[i]/2-1+0.1),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos(2*pi)*(list_da[i]/2-1+0.1)),), ((sin(2*pi)*(list_da[i]/2-1+0.1),l_absatz+f[i],cos(2*pi)*(list_da[i]/2-1+0.1)),), ((sin(2*pi)*(list_da[i]/2-1+0.1),l_absatz+f[i]+4,cos(2*pi)*(list_da[i]/2-1+0.1)),), ((sin(2*pi)*(list_da[i]/2-1+0.1),(l_absatz+f[i]+list_lf[i]),cos(2*pi)*(list_da[i]/2-1+0.1)),), ((sin(2*pi)*(list_da[i]/2-1+0.1),(l_absatz+f[i]+list_lf[i]/2),cos(2*pi)*(list_da[i]/2-1+0.1)),),) Edge2 =e_w.findAt(((sin(pi/2)*(list_da[i]/2-1+0.1),l_absatz+radius[i],cos(pi/2)*(list_da[i]/2-1+0.1)),), ((sin(pi/2)*(list_da[i]/2-1+0.1),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos(pi/2)*(list_da[i]/2-1+0.1)),), ((sin(pi/2)*(list_da[i]/2-1+0.1),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos(pi/2)*(list_da[i]/2-1+0.1)),), ((sin(pi/2)*(list_da[i]/2-1+0.1),l_absatz+f[i],cos(pi/2)*(list_da[i]/2-1+0.1)),), ((sin(pi/2)*(list_da[i]/2-1+0.1),l_absatz+f[i]+4,cos(pi/2)*(list_da[i]/2-1+0.1)),), ((sin(pi/2)*(list_da[i]/2-1+0.1),(l_absatz+f[i]+list_lf[i]),cos(pi/2)*(list_da[i]/2-1+0.1)),), ((sin(pi/2)*(list_da[i]/2-1+0.1),(l_absatz+f[i]+list_lf[i]/2),cos(pi/2)*(list_da[i]/2-1+0.1)),), ((sin(pi)*(list_da[i]/2-1+0.1),l_absatz+radius[i],cos(pi)*(list_da[i]/2-1+0.1)),), ((sin(pi)*(list_da[i]/2-1+0.1),(l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin((winkel/2)*pi/180)),cos(pi)*(list_da[i]/2-1+0.1)),), ((sin(pi)*(list_da[i]/2-1+0.1),(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),cos(pi)*(list_da[i]/2-1+0.1)),), ((sin(pi)*(list_da[i]/2-1+0.1),l_absatz+f[i],cos(pi)*(list_da[i]/2-1+0.1)),), ((sin(pi)*(list_da[i]/2-1+0.1),l_absatz+f[i]+4,cos(pi)*(list_da[i]/2-1+0.1)),), ((sin(pi)*(list_da[i]/2-1+0.1),(l_absatz+f[i]+list_lf[i]),cos(pi)*(list_da[i]/2-1+0.1)),), ((sin(pi)*(list_da[i]/2-1+0.1),(l_absatz+f[i]+list_lf[i]/2),cos(pi)*(list_da[i]/2-1+0.1)),),) welle.seedEdgeByBias(biasMethod=SINGLE, end1Edges=Edge1, end2Edges=Edge2, ratio=4.0, number=FU, constraint=FINER) Edge1 = e_w.findAt(((sin(3*pi/2)*(list_da[i]/2),l_absatz-2,cos(3*pi/2)*(list_da[i]/2)),), ((sin(3*pi/2)*(list_da[i]/2),l_absatz-20,cos(3*pi/2)*(list_da[i]/2)),), ((sin(2*pi)*(list_da[i]/2),l_absatz-2,cos(2*pi)*(list_da[i]/2)),), ((sin(2*pi)*(list_da[i]/2),l_absatz-20,cos(2*pi)*(list_da[i]/2)),),) Edge2 = e_w.findAt(((sin(pi/2)*(list_da[i]/2),l_absatz-2,cos(pi/2)*(list_da[i]/2)),), ((sin(pi/2)*(list_da[i]/2),l_absatz-20,cos(pi/2)*(list_da[i]/2)),), ((sin(pi)*(list_da[i]/2),l_absatz-2,cos(pi)*(list_da[i]/2)),), ((sin(pi)*(list_da[i]/2),l_absatz-20,cos(pi)*(list_da[i]/2)),),) welle.seedEdgeByBias(biasMethod=SINGLE, end1Edges=Edge1, end2Edges=Edge2, ratio=4.0, number=FU, constraint=FINER) welle.seedEdgeByNumber(edges=Radialedge_FR_Fugenflaeche_i, number=EI, constraint=FIXED) welle.seedEdgeByNumber(edges=Radialedge_Anfangsteil, number=EI, constraint=FIXED) Edge1 = e_w.findAt(((sin(pi/2)*(list_da[i]/2+radius[i]+0.01),l_absatz-2,cos(pi/2)*(list_da[i]/2+radius[i]+0.01)),), ((sin(pi/2)*(list_da[i]/2+radius[i]+0.01),l_absatz-20,cos(pi/2)*(list_da[i]/2+radius[i]+0.01)),), ((sin(pi/2)*(list_da[i]/2+radius[i]+0.01),l_absatz,cos(pi/2)*(list_da[i]/2+radius[i]+0.01)),), ((sin(pi)*(list_da[i]/2+radius[i]+0.01),l_absatz-2,cos(pi)*(list_da[i]/2+radius[i]+0.01)),), ((sin(pi)*(list_da[i]/2+radius[i]+0.01),l_absatz-20,cos(pi)*(list_da[i]/2+radius[i]+0.01)),), ((sin(3*pi/2)*(list_da[i]/2+radius[i]+0.01),l_absatz,cos(3*pi/2)*(list_da[i]/2+radius[i]+0.01)),), ((sin(2*pi)*(list_da[i]/2+radius[i]+0.01),l_absatz,cos(2*pi)*(list_da[i]/2+radius[i]+0.01)),),) Edge2 = e_w.findAt(((sin(pi)*(list_da[i]/2+radius[i]+0.01),l_absatz,cos(pi)*(list_da[i]/2+radius[i]+0.01)),), ((sin(3*pi/2)*(list_da[i]/2+radius[i]+0.01),l_absatz-2,cos(3*pi/2)*(list_da[i]/2+radius[i]+0.01)),), ((sin(3*pi/2)*(list_da[i]/2+radius[i]+0.01),l_absatz-20,cos(3*pi/2)*(list_da[i]/2+radius[i]+0.01)),), ((sin(2*pi)*(list_da[i]/2+radius[i]+0.01),l_absatz-2,cos(2*pi)*(list_da[i]/2+radius[i]+0.01)),), ((sin(2*pi)*(list_da[i]/2+radius[i]+0.01),l_absatz-20,cos(2*pi)*(list_da[i]/2+radius[i]+0.01)),),) welle.seedEdgeByBias(biasMethod=SINGLE, end1Edges=Edge1, end2Edges=Edge2, ratio=4.0, number=EI, constraint=FINER) welle.seedEdgeByNumber(edges=Radialedge_Auslauf, number=EI, constraint=FIXED) welle = mdb.models[Modell_Name[i]].parts[Wellen_Name[i]] welle.generateMesh() welle.generateMesh(seedConstraintOverride=ON) # Nabe #tangential nabe = mdb.models[Modell_Name[i]].parts[Naben_Name[i]] e_n = nabe.edges nabe.seedEdgeByNumber(edges=Tangentialedge_Nabe, number=FU, constraint=FIXED) #axial nabe.seedEdgeByNumber(edges=Axialedge_Kontaktfrei, number=ZW, constraint=FIXED) nabe.seedEdgeByNumber(edges=Axialedge_Kontaktflaeche_o, number=VI, constraint=FIXED) nabe.seedEdgeByNumber(edges=Axialedge_Kontaktflaeche_u, number=VI, constraint=FIXED) Edge1 = e_n.findAt(((sin (pi/2)*list_Di[i]/2, l_absatz+f[i]+0.01, cos (pi/2)*list_Di[i]/2),), ((sin (pi/2)*(list_Di[i]/2+1), l_absatz+f[i]+0.01, cos (pi/2)*(list_Di[i]/2+1)),), ((sin (pi)*(list_Di[i]/2+1), l_absatz+f[i]+0.01, cos (pi)*(list_Di[i]/2+1)),), ((sin (pi)*list_Da[i]/2, l_absatz+f[i]+0.01, cos (pi)*list_Da[i]/2),), ((sin (3*pi/2)*list_Di[i]/2, l_absatz+f[i]+0.01, cos (3*pi/2)*list_Di[i]/2),), ((sin (2*pi)*list_Di[i]/2, l_absatz+f[i]+0.01, cos (2*pi)*list_Di[i]/2),),) Edge2 = e_n.findAt(((sin (pi/2)*list_Da[i]/2, l_absatz+f[i]+0.01, cos (pi/2)*list_Da[i]/2),), ((sin (pi)*list_Di[i]/2, l_absatz+f[i]+0.01, cos (pi)*list_Di[i]/2),), ((sin (3*pi/2)*(list_Di[i]/2+1), l_absatz+f[i]+0.01, cos (3*pi/2)*(list_Di[i]/2+1)),), ((sin (3*pi/2)*list_Da[i]/2, l_absatz+f[i]+0.01, cos (3*pi/2)*list_Da[i]/2),), ((sin (2*pi)*(list_Di[i]/2+1), l_absatz+f[i]+0.01, cos (2*pi)*(list_Di[i]/2+1)),), ((sin (2*pi)*list_Da[i]/2, l_absatz+f[i]+0.01, cos (2*pi)*list_Da[i]/2),),) nabe.seedEdgeByBias(biasMethod=SINGLE, end1Edges=Edge1, end2Edges=Edge2, ratio=4.0, number=SE, constraint=FINER) #radial Edge1 = e_n.findAt(((sin (pi/2)*(list_Di[i]/2+0.01), l_absatz+f[i], cos (pi/2)*(list_Di[i]/2+0.01)),), ((sin (pi/2)*(list_Di[i]/2+0.01), l_absatz+f[i]+list_lf[i], cos (pi/2)*(list_Di[i]/2+0.01)),), ((sin (pi/2)*(list_Di[i]/2+0.01), l_absatz+f[i]+list_lf[i]/2, cos (pi/2)*(list_Di[i]/2+0.01)),), ((sin (pi)*(list_Di[i]/2+0.01), l_absatz+1.5, cos (pi)*(list_Di[i]/2+0.01)),), ((sin (pi)*(list_Di[i]/2+0.01), l_absatz+f[i], cos (pi)*(list_Di[i]/2+0.01)),), ((sin (3*pi/2)*(list_Di[i]/2+0.01), l_absatz+f[i]+list_lf[i], cos (3*pi/2)*(list_Di[i]/2+0.01)),), ((sin (pi)*(list_Di[i]/2+0.01), l_absatz+f[i]+list_lf[i]/2, cos (pi)*(list_Di[i]/2+0.01)),), ((sin (2*pi)*(list_Di[i]/2+0.01), l_absatz+f[i]+list_lf[i], cos (2*pi)*(list_Di[i]/2+0.01)),), ((sin (pi/2)*(list_Di[i]/2+0.01), l_absatz+f[i]+4, cos (pi/2)*(list_Di[i]/2+0.01)),), ((sin (pi)*(list_Di[i]/2+0.01), l_absatz+f[i]+4, cos (pi)*(list_Di[i]/2+0.01)),),) Edge2 = e_n.findAt(((sin (pi/2)*(list_Di[i]/2+0.01), l_absatz+1.5, cos (pi/2)*(list_Di[i]/2+0.01)),), ((sin (pi)*(list_Di[i]/2+0.01), l_absatz+f[i]+list_lf[i], cos (pi)*(list_Di[i]/2+0.01)),), ((sin (3*pi/2)*(list_Di[i]/2+0.01), l_absatz+1.5, cos (3*pi/2)*(list_Di[i]/2+0.01)),), ((sin (3*pi/2)*(list_Di[i]/2+0.01), l_absatz+f[i], cos (3*pi/2)*(list_Di[i]/2+0.01)),), ((sin (2*pi)*(list_Di[i]/2+0.01), l_absatz+1.5, cos (2*pi)*(list_Di[i]/2+0.01)),), ((sin (2*pi)*(list_Di[i]/2+0.01), l_absatz+f[i], cos (2*pi)*(list_Di[i]/2+0.01)),), ((sin (3*pi/2)*(list_Di[i]/2+0.01), l_absatz+f[i]+4, cos (3*pi/2)*(list_Di[i]/2+0.01)),), ((sin (2*pi)*(list_Di[i]/2+0.01), l_absatz+f[i]+4, cos (2*pi)*(list_Di[i]/2+0.01)),), ((sin (3*pi/2)*(list_Di[i]/2+0.01), l_absatz+f[i]+list_lf[i]/2, cos (3*pi/2)*(list_Di[i]/2+0.01)),), ((sin (2*pi)*(list_Di[i]/2+0.01), l_absatz+f[i]+list_lf[i]/2, cos (2*pi)*(list_Di[i]/2+0.01)),),) nabe.seedEdgeByBias(biasMethod=SINGLE, end1Edges=Edge1, end2Edges=Edge2, ratio=9.0, number=EI, constraint=FINER) nabe.seedEdgeByNumber(edges=Radialedge_aussen, number=EI, constraint=FIXED) nabe = mdb.models[Modell_Name[i]].parts[Naben_Name[i]] nabe.generateMesh() # Netzerstellung # Welle innen if i > 8 and i < 13: #tangential welle_innen = mdb.models[Modell_Name[i]].parts[Welle_innen_Name[i]] e_wi = welle_innen.edges welle_innen.seedEdgeByNumber(edges=Tangentialedge_Wi, number=FU, constraint=FIXED) #axial welle_innen.seedEdgeByNumber(edges=Axialedge_Fr_horiz_Wi, number=DR, constraint=FIXED) welle_innen.seedEdgeByNumber(edges=Axialedge_Fr_2R_Wi, number=DR, constraint=FIXED) welle_innen.seedEdgeByNumber(edges=Axialedge_Fr_schraeg_Wi, number=VI, constraint=FIXED) Edge1 = e_wi.findAt(((sin (pi/2)*(da_i/2),l_absatz,cos (pi/2)*(da_i/2)),), ((sin (3*pi/2)*(da_i/2),l_absatz,cos (3*pi/2)*(da_i/2)),), ((sin (2*pi)*(da_i/2),l_absatz,cos (2*pi)*(da_i/2)),), ((sin (3*pi/2)*(da_i/2-1),l_absatz,cos (3*pi/2)*(da_i/2-1)),), ((sin (2*pi)*(da_i/2-1),l_absatz,cos (2*pi)*(da_i/2-1)),),) Edge2 = e_wi.findAt(((sin (pi/2)*(da_i/2-1),l_absatz,cos (pi/2)*(da_i/2-1)),), ((sin (pi)*(da_i/2),l_absatz,cos (pi)*(da_i/2)),), ((sin (pi)*(da_i/2-1),l_absatz,cos (pi)*(da_i/2-1)),), ((0,l_absatz,0),),) welle_innen.seedEdgeByBias(biasMethod=SINGLE, end1Edges=Edge1, end2Edges=Edge2, ratio=9, number=FU, constraint=FINER) welle_innen.seedEdgeByNumber(edges=Axialedge_Fugenflaeche_o_Wi, number=VI, constraint=FIXED) welle_innen.seedEdgeByNumber(edges=Axialedge_Fugenflaeche_u_Wi, number=VI, constraint=FIXED) Edge1 = e_wi.findAt(((sin (pi/2)*(da_i/2),l_absatz-5,cos (pi/2)*(da_i/2)),), ((sin (3*pi/2)*(da_i/2),l_absatz-5,cos (3*pi/2)*(da_i/2)),), ((sin (3*pi/2)*(da_i/2-1),l_absatz-5,cos (3*pi/2)*(da_i/2-1)),), ((sin (2*pi)*(da_i/2),l_absatz-5,cos (2*pi)*(da_i/2)),), ((sin (2*pi)*(da_i/2-1),l_absatz-5,cos (2*pi)*(da_i/2-1)),),) Edge2 = e_wi.findAt(((sin (pi/2)*(da_i/2-1),l_absatz-5,cos (pi/2)*(da_i/2-1)),), ((sin (pi)*(da_i/2),l_absatz-5,cos (pi)*(da_i/2)),), ((sin (pi)*(da_i/2-1),l_absatz-5,cos (pi)*(da_i/2-1)),), ((0,l_absatz-5,0),),) welle_innen.seedEdgeByBias(biasMethod=SINGLE, end1Edges=Edge1, ratio=4, number=FU, constraint=FIXED) welle_innen.seedEdgeByNumber(edges=Axialedge_Auslauf_Wi, number=EI, constraint=FIXED) welle_innen.seedEdgeByNumber(edges=Axialedge_Absatz_Wi, number=EI, constraint=FIXED) Edge2 = e_wi.findAt(((sin (pi/2)*(da_i/2),l_absatz+f[i]+3,cos (pi/2)*(da_i/2)),), ((sin (3*pi/2)*(da_i/2),l_absatz+f[i]+3,cos (3*pi/2)*(da_i/2)),), ((sin (3*pi/2)*(da_i/2-1),l_absatz+f[i]+3,cos (3*pi/2)*(da_i/2-1)),), ((sin (2*pi)*(da_i/2),l_absatz+f[i]+3,cos (2*pi)*(da_i/2)),), ((sin (2*pi)*(da_i/2-1),l_absatz+f[i]+3,cos (2*pi)*(da_i/2-1)),),) Edge1 = e_wi.findAt(((sin (pi/2)*(da_i/2-1),l_absatz+f[i]+3,cos (pi/2)*(da_i/2-1)),), ((sin (pi)*(da_i/2),l_absatz+f[i]+3,cos (pi)*(da_i/2)),), ((sin (pi)*(da_i/2-1),l_absatz+f[i]+3,cos (pi)*(da_i/2-1)),), ((0,l_absatz+f[i]+3,0),),) welle_innen.seedEdgeByBias(biasMethod=SINGLE, end1Edges=Edge1, end2Edges=Edge2, ratio=4, number=SE, constraint=FIXED) #radial welle_innen.seedEdgeByNumber(edges=Radialedge_FR_Fugenflaeche_a_Wi, number=EI, constraint=FIXED) welle_innen.seedEdgeByNumber(edges=Radialedge_Rundung_Wi, number=EI, constraint=FIXED) welle_innen.seedEdgeByNumber(edges=Radialedge_FR_Fugenflaeche_i_Wi, number=EI, constraint=FIXED) welle_innen.seedEdgeByNumber(edges=Radialedge_Anfangsteil_Wi, number=EI, constraint=FIXED) welle_innen.seedEdgeByNumber(edges=Radialedge_Auslauf_Wi, number=EI, constraint=FIXED) welle_innen = mdb.models[Modell_Name[i]].parts[Welle_innen_Name[i]] welle_innen.generateMesh() welle_innen.generateMesh(seedConstraintOverride=ON) welle_innen = mdb.models[Modell_Name[i]].parts[Welle_innen_Name[i]] c_wi = welle_innen.cells pickedRegions = c_wi.getByBoundingBox(-da_i/2,l_absatz-22,-da_i/2,da_i/2,l_absatz-2,da_i/2) welle_innen.deleteMesh(regions=pickedRegions) welle_innen = mdb.models[Modell_Name[i]].parts[Welle_innen_Name[i]] e_wi = welle_innen.edges Edge1 = e_wi.findAt(((sin (pi/2)*(da_i/2),l_absatz-5,cos (pi/2)*(da_i/2)),), ((sin (3*pi/2)*(da_i/2),l_absatz-5,cos (3*pi/2)*(da_i/2)),), ((sin (3*pi/2)*(da_i/2-1),l_absatz-5,cos (3*pi/2)*(da_i/2-1)),), ((sin (2*pi)*(da_i/2),l_absatz-5,cos (2*pi)*(da_i/2)),), ((sin (2*pi)*(da_i/2-1),l_absatz-5,cos (2*pi)*(da_i/2-1)),),) Edge2 = e_wi.findAt(((sin (pi/2)*(da_i/2-1),l_absatz-5,cos (pi/2)*(da_i/2-1)),), ((sin (pi)*(da_i/2),l_absatz-5,cos (pi)*(da_i/2)),), ((sin (pi)*(da_i/2-1),l_absatz-5,cos (pi)*(da_i/2-1)),), ((0,l_absatz-5,0),),) welle_innen.seedEdgeByBias(biasMethod=SINGLE, end1Edges=Edge1, ratio=4, number=FU, constraint=FIXED) welle_innen = mdb.models[Modell_Name[i]].parts[Welle_innen_Name[i]] welle_innen.generateMesh() # Field Output mdb.models[Modell_Name[i]].fieldOutputRequests['F-Output-1'].setValues( variables=('S', 'MISES', 'MISESMAX', 'E', 'EE', 'U', 'RF', 'CF', 'CSTRESS', 'CDISP'), frequency=LAST_INCREMENT) mdb.models[Modell_Name[i]].historyOutputRequests['H-Output-1'].setValues( variables=('CSTRESS', ), frequency=LAST_INCREMENT) # Node Sets # Wellenoberfläche session.viewports['Viewport: 1'].partDisplay.setValues(mesh=ON) welle = mdb.models[Modell_Name[i]].parts[Wellen_Name[i]] n_w = welle.nodes nodes_axial_1 = n_w.getByBoundingBox(-0.1, list_da[i]+22-0.01,list_da[i]/2+radius[i],0.1,list_da[i]+22+0.01,list_da[i]/2+h_absatz[i]) nodes_axial_2 = n_w.getByBoundingSphere((0, l_absatz+radius[i], list_da[i]/2+(radius[i]-t1[i])), radius[i]) nodes_axial_3 = n_w.getByBoundingBox(-0.1, list_da[i]+22+radius[i],list_da[i]/2-t1[i]-0.01,0.1,(l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))-2*radius[i]*sin(winkel/2*pi/180),list_da[i]/2-t1[i]+0.01) nodes_axial_4 = n_w.getByBoundingSphere((0, (l_absatz+f[i])-(((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180))+2*radius[i]*sin(winkel/2*pi/180)), list_da[i]/2+(radius[i]-t1[i])), radius[i]+0.01) nodes_axial_5 = n_w.getByBoundingBox(-0.1, (l_absatz+f[i])-((radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i]))/tan(winkel*pi/180)),list_da[i]/2-(radius[i]*sin((90-winkel)*pi/180)-(radius[i]-t1[i])),0.1,l_absatz+f[i]+0.01,list_da[i]/2) nodes_axial_6 = n_w.getByBoundingBox(-0.1, l_absatz+f[i], list_da[i]/2-0.01, 0.1, l_absatz+f[i]+4.01, list_da[i]/2+0.01) nodes_axial = (nodes_axial_1, nodes_axial_2, nodes_axial_3, nodes_axial_4, nodes_axial_5, nodes_axial_6) welle.Set(nodes=nodes_axial, name='Wellenfreistich_axial') session.viewports['Viewport: 1'].partDisplay.setValues(mesh=OFF) # Node Sets # Wellenoberfläche_Startpunkt session.viewports['Viewport: 1'].partDisplay.setValues(mesh=ON) welle = mdb.models[Modell_Name[i]].parts[Wellen_Name[i]] n_w = welle.nodes nodes_axial_start = n_w.getByBoundingBox(-0.1, list_da[i]+22-0.01,list_da[i]/2+h_absatz[i]-0.01,0.1,list_da[i]+22+0.01,list_da[i]/2+h_absatz[i]+0.01) welle.Set(nodes=nodes_axial_start, name='Wellenfreistich_axial_Startpunkt') session.viewports['Viewport: 1'].partDisplay.setValues(mesh=OFF) # Node Sets # Nabenkante radial session.viewports['Viewport: 1'].partDisplay.setValues(mesh=ON) welle = mdb.models[Modell_Name[i]].parts[Wellen_Name[i]] n_w = welle.nodes nodes_radial = n_w.getByBoundingBox(-0.1, list_da[i]+22+f[i]-0.01,-list_da[i]/2,0.1,list_da[i]+22+f[i]+0.01,list_da[i]/2+0.01) welle.Set(nodes=nodes_radial, name='Nabenkante_radial') session.viewports['Viewport: 1'].partDisplay.setValues(mesh=OFF) # Node Sets # Nabenkante radial Startpunkt session.viewports['Viewport: 1'].partDisplay.setValues(mesh=ON) welle = mdb.models[Modell_Name[i]].parts[Wellen_Name[i]] n_w = welle.nodes nodes_radial_start = n_w.getByBoundingBox(-0.1, list_da[i]+22+f[i]-0.01,list_da[i]/2-0.01,0.1,list_da[i]+22+f[i]+0.01,list_da[i]/2+0.01) welle.Set(nodes=nodes_radial_start, name='Nabenkante_radial_Startpunkt') session.viewports['Viewport: 1'].partDisplay.setValues(mesh=OFF) # Job erzeugen # wnv = mdb.models[Modell_Name[i]].rootAssembly wnv.regenerate() wnv = mdb.models[Modell_Name[i]].rootAssembly session.viewports['Viewport: 1'].setValues(displayedObject=wnv) session.viewports['Viewport: 1'].assemblyDisplay.setValues(mesh=OFF) session.viewports['Viewport: 1'].assemblyDisplay.meshOptions.setValues( meshTechnique=OFF) mdb.Job(name=Modell_Name[i], model=Modell_Name[i], description='', type=ANALYSIS, atTime=None, waitMinutes=0, waitHours=0, queue=None, memory=90, memoryUnits=PERCENTAGE, getMemoryFromAnalysis=True, explicitPrecision=SINGLE, nodalOutputPrecision=SINGLE, echoPrint=OFF, modelPrint=OFF, contactPrint=OFF, historyPrint=OFF, userSubroutine='', scratch='', resultsFormat=ODB, multiprocessingMode=DEFAULT, numCpus=1, numGPUs=0) #Arbeitsverzeichnis wechseln os.chdir('C:/temp/V-B_B1') # Wechsel os.getcwd() # Anzeige welches Arbeitsverzeichnis mdb.jobs[Modell_Name[i]].writeInput(consistencyChecking=OFF) #mdb.jobs[Modell_Name[i]].submit(consistencyChecking=OFF) mdb.jobs[Modell_Name[i]].waitForCompletion() o3 = session.openOdb(name=odb_Pfad[i]) session.viewports['Viewport: 1'].setValues(displayedObject=o3) session.viewports['Viewport: 1'].makeCurrent() session.mdbData.summary() session.viewports['Viewport: 1'].odbDisplay.display.setValues(plotState=( CONTOURS_ON_UNDEF, )) #Paths erzeugen Wellenfreistich axial session.viewports['Viewport: 1'].partDisplay.setValues(mesh=ON) set_nodes = 'Nabenkante_radial' set_startnode = 'Nabenkante_radial_Startpunkt' instancename = Wellen_Name[i] pathname = 'Nabenkante_radial' dict_nodes = {} list_nodes = [] set_nodes = set_nodes.upper() set_startnode = set_startnode.upper() instancename = instancename.upper() vps = session.viewports.values()[0] odbName = vps.displayedObject.name odb = session.odbs[odbName] coord_point = odb.rootAssembly.instances[instancename].nodeSets[set_startnode].nodes[0].coordinates for node in odb.rootAssembly.instances[instancename].nodeSets[set_nodes].nodes: dist = sqrt((pow(node.coordinates[0]-coord_point[0],2))+ \ (pow(node.coordinates[1]-coord_point[1],2))+ \ (pow(node.coordinates[2]-coord_point[2],2))) dict_nodes[dist] = node.label x = dict_nodes.keys() x.sort() for y in x: list_nodes.append(dict_nodes[y]) list_nodes = tuple(list_nodes) session.Path(name=pathname, type=NODE_LIST, expression=((instancename, (list_nodes)), )) #Paths erzeugen Wellenfreistich axial session.viewports['Viewport: 1'].partDisplay.setValues(mesh=ON) set_nodes = 'Wellenfreistich_axial' set_startnode = 'Wellenfreistich_axial_Startpunkt' instancename = Wellen_Name[i] pathname = 'Wellenfreistich_axial' dict_nodes = {} list_nodes = [] set_nodes = set_nodes.upper() set_startnode = set_startnode.upper() instancename = instancename.upper() vps = session.viewports.values()[0] odbName = vps.displayedObject.name odb = session.odbs[odbName] coord_point = odb.rootAssembly.instances[instancename].nodeSets[set_startnode].nodes[0].coordinates for node in odb.rootAssembly.instances[instancename].nodeSets[set_nodes].nodes: dist = sqrt((pow(node.coordinates[0]-coord_point[0],2))+ \ (pow(node.coordinates[1]-coord_point[1],2))+ \ (pow(node.coordinates[2]-coord_point[2],2))) dict_nodes[dist] = node.label x = dict_nodes.keys() x.sort() for y in x: list_nodes.append(dict_nodes[y]) list_nodes = tuple(list_nodes) session.Path(name=pathname, type=NODE_LIST, expression=((instancename, (list_nodes)), )) # XY DATA # Wellenoberflaeche_axial session.viewports['Viewport: 1'].odbDisplay.setFrame(step=1, frame=1) session.viewports['Viewport: 1'].odbDisplay.setFrame(step=0, frame=1) pth = session.paths['Wellenfreistich_axial'] session.XYDataFromPath(name='Wellenfreistich_axial', path=pth, includeIntersections=False, projectOntoMesh=False, pathStyle=PATH_POINTS, numIntervals=10, projectionTolerance=0, shape=UNDEFORMED, labelType=TRUE_DISTANCE_Y, removeDuplicateXYPairs=True, includeAllElements=False) xyp = session.XYPlot('Wellenfreistich_axial') chartName = xyp.charts.keys()[0] chart = xyp.charts[chartName] pth = session.paths['Wellenfreistich_axial'] xy1 = xyPlot.XYDataFromPath(path=pth, includeIntersections=False, projectOntoMesh=False, pathStyle=PATH_POINTS, numIntervals=10, projectionTolerance=0, shape=UNDEFORMED, labelType=TRUE_DISTANCE_Y, removeDuplicateXYPairs=True, includeAllElements=False) c1 = session.Curve(xyData=xy1) chart.setValues(curvesToPlot=(c1, ), ) session.charts[chartName].autoColor(lines=True, symbols=True) session.viewports['Viewport: 1'].setValues(displayedObject=xyp) # XY DATA # Wellenoberflaeche_axial session.viewports['Viewport: 1'].odbDisplay.setFrame(step=1, frame=1) session.viewports['Viewport: 1'].odbDisplay.setFrame(step=0, frame=1) pth = session.paths['Nabenkante_radial'] session.XYDataFromPath(name='Nabenkante_radial', path=pth, includeIntersections=False, projectOntoMesh=False, pathStyle=PATH_POINTS, numIntervals=10, projectionTolerance=0, shape=UNDEFORMED, labelType=TRUE_DISTANCE_Y, removeDuplicateXYPairs=True, includeAllElements=False) xyp = session.XYPlot('Nabenkante_radial') chartName = xyp.charts.keys()[0] chart = xyp.charts[chartName] pth = session.paths['Nabenkante_radial'] xy1 = xyPlot.XYDataFromPath(path=pth, includeIntersections=False, projectOntoMesh=False, pathStyle=PATH_POINTS, numIntervals=10, projectionTolerance=0, shape=UNDEFORMED, labelType=TRUE_DISTANCE_Y, removeDuplicateXYPairs=True, includeAllElements=False) c1 = session.Curve(xyData=xy1) chart.setValues(curvesToPlot=(c1, ), ) session.charts[chartName].autoColor(lines=True, symbols=True) session.viewports['Viewport: 1'].setValues(displayedObject=xyp)