# -*- coding: mbcs -*-
from part import *																								#How to read script
from material import *																							#1.Define Materials
from section import *																								#2.Create Sketch Biegeblech
from assembly import *																							#3.Create sketch Niederhalter
from step import *																								#4.Create sketch Biegeform_2
from interaction import *																							#5.Create sketch Biegekern
from load import *																								#6.Create sketch final position
from mesh import *																								#7.Create part Endposition
from optimization import *																						#8.Create part Biegeform_2
from job import *																									#9.Create part Biegeblech (incl. instances, mesh, seed, partition, elementtype, sets, surfaces, section assignments)
from sketch import *																								#10.Create part Niederhalter
from visualization import *																						#11.Create part Biegekern
from connectorBehavior import *																					#12.Create instances 
from abaqusConstants import * 																					#13.Create position constraints
import regionToolset 																							#14.Create Reference Points
import sys, math																									#15.Create inertias
import os																										#16.Create surfaces
																												#17.Create constraints
																												#18.Create step
# insert parameters and always use data type float!!!! example=xx.x												#19.Create contact property
																												#20.Create interaction properties
for rk in [2.0]:						#Biegekernradius [mm]														#21.Create boundary conditions
	for alpha in [90.0]:				#Stempelgeschwindigkeit [mm/s]												#22.Create load
																												#23.Define Field Output
		t = 1.0;						#Blechdicke [mm]															#24.Submit and run job
		b = 30.0;					#Blechbreite [mm]															#25.Import Explicit to Standard
		#s = 2;						#Stempelform []																#26.New properties for Springback Analysis
		V = 10.0;					#Stempelgeschwindigkeit														#27.Submit and run job Springback Analysis
		Zeit=(12.0/V);				#Zeitangabe 12/Stempelgeschiwndigkeit
		alpha = 30.0;				#Biegewinkel [°] 
		materialname='AL_d0';		#Materialangabe mit Namen wie definiert => Walzrichtung beachten
		mass_scale = 5e-08;			#Massenskalierung
		
		
		Simu_name =  '3D-Simulation_t'+str(int(t))+'_b'+str(int(b))+'_'+materialname+'_rk'+str(int(rk))+'_v'+str(int(V))+'_alpha'+str(int(alpha));
		#Simu_name = Dicke+Breite+Material+Biegekernradius+Stempelgeschwindigkeit+Biegewinkel

		Simu_3D = mdb.Model(modelType=STANDARD_EXPLICIT , name=Simu_name)	
		############################################### 1.Define Materials ########################################################################
		# Define Material
		#Al_d0
		AL_d0_EModul = 68000.00;

		Al_d0 = Simu_3D.Material(name='AL_d0')
		Al_d0.Density(table=((2.7e-09, ), ))
		Al_d0.Elastic(table=((AL_d0_EModul, 0.34), ))
		Al_d0.Plastic(table=((137.1987687, 0.0), (156.7938262, 0.01), (174.3143622, 0.02), (189.9800054, 0.03), 
			(203.9871327, 0.04), (216.5113306, 0.05), (227.7095963, 0.06), (237.7223056, 0.07), (246.6749731, 0.08), (254.679825, 0.09), (261.8372061, 0.1), 
			(268.2368379, 0.11), (273.958943, 0.12), (279.075251, 0.13), (283.6498973, 0.14), (287.7402274, 0.15), (291.3975159, 0.16), (294.6676085, 0.17), 
			(297.5914976, 0.18), (300.2058357, 0.19), (302.5433947, 0.2), (304.6334772, 0.21), (306.5022835, 0.22), (308.1732399, 0.23)))
		Al_d0.plastic.Potential(table=((1.0, 1.016956108, 0.929516003, 1.186884852, 1.0, 1.0), ))		

		#Al_d90
		AL_d90_EModul = 68000.00;
	
		Al_d90 = Simu_3D.Material(name='AL_d90')
		Al_d90.Density(table=((2.7e-09, ), ))
		Al_d90.Elastic(table=((AL_d90_EModul, 0.34), ))
		Al_d90.Plastic(table=((137.1987687, 0.0), (156.7938262, 0.01), (174.3143622, 0.02), (189.9800054, 0.03), 
			(203.9871327, 0.04), (216.5113306, 0.05), (227.7095963, 0.06), (237.7223056, 0.07), (246.6749731, 0.08), (254.679825, 0.09), (261.8372061, 0.1), 
			(268.2368379, 0.11), (273.958943, 0.12), (279.075251, 0.13), (283.6498973, 0.14), (287.7402274, 0.15), (291.3975159, 0.16), (294.6676085, 0.17), 
			(297.5914976, 0.18), (300.2058357, 0.19), (302.5433947, 0.2), (304.6334772, 0.21), (306.5022835, 0.22), (308.1732399, 0.23)))
		Al_d90.plastic.Potential(table=((1.016956108, 1, 0.929516003, 1.186884852, 1.0, 1.0), ))		
				
		#DC04_d0
		DC04_d0_EModul = 210000.00;

		DC04_d0 = Simu_3D.Material(name='DC04_d0')
		DC04_d0.Density(table=((7.86e-09, ), ))
		DC04_d0.Elastic(table=((DC04_d0_EModul, 0.3), ))
		DC04_d0.Plastic(table=((132.3685966, 0.0), (196.8397638, 0.01), (218.1851987, 0.02), (233.8219806, 0.03), 
			(246.6211611, 0.04), (257.6548438, 0.05), (267.4595395, 0.06), (276.3480311, 0.07), (284.5210199, 0.08), (292.1161215, 0.09), (299.2324151, 0.1), 
			(305.9439592, 0.11), (312.307781, 0.12), (318.3688675, 0.13), (324.1634267, 0.14), (329.7210971, 0.15), (335.0664935, 0.16), (340.2203152, 0.17), 
			(345.20016, 0.18), (350.0211335, 0.19), (354.6963132, 0.2), (359.2371076, 0.21), (363.6535378, 0.22), (367.9544609, 0.23), (372.1477492, 0.24), 
			(376.2404364, 0.25), (380.2388361, 0.26), (384.1486404, 0.27), (387.9750013, 0.28), (391.7225999, 0.29), (395.3957033, 0.3), (398.9982142, 0.31), 
			(402.5337126, 0.32), (406.0054915, 0.33), (409.4165878, 0.34), (412.7698096, 0.35), (416.0677587, 0.36), (419.3128517, 0.37), (422.5073374, 0.38), 
			(425.6533126, 0.39), (428.752736, 0.4), (431.8074407, 0.41), (434.8191448, 0.42), (437.7894613, 0.43), (440.7199069, 0.44), (443.6119095, 0.45), 
			(446.4668157, 0.46), (449.2858966, 0.47), (452.0703539, 0.48), (454.821325, 0.49), (457.5398873, 0.5), (460.2270631, 0.51), (462.8838231, 0.52), 
			(465.5110897, 0.53), (468.1097409, 0.54), (470.6806126, 0.55), (473.2245017, 0.56), (475.7421685, 0.57), (478.2343388, 0.58), (480.7017064, 0.59),
			(483.1449348, 0.6), (485.564659, 0.61), (487.9614871, 0.62), (490.336002, 0.63), (492.688763, 0.64), (495.0203065, 0.65), (497.331148, 0.66), 
			(499.6217827, 0.67), (501.8926867, 0.68), (504.1443184, 0.69), (506.3771188, 0.7), (508.5915126, 0.71), (510.7879093, 0.72), (512.9667038, 0.73), 
			(515.1282768, 0.74), (517.272996, 0.75), (519.4012162, 0.76), (521.5132805, 0.77), (523.6095203, 0.78), (525.690256, 0.79), (527.7557978, 0.8)))
		DC04_d0.plastic.RateDependent(table=((0.995, 0.0), (1.0, 0.004), (1.005003134, 0.008), (1.007941375, 0.012), 
			(1.010031299, 0.016), (1.011655353, 0.02), (1.012984241, 0.024), (1.014109161, 0.028), (1.015084621, 0.032), (1.015945816, 0.036), (1.0167168, 0.04), 
			(1.017414745, 0.044), (1.018052336, 0.048), (1.018639217, 0.052), (1.019182885, 0.056), (1.019689288, 0.06), (1.020163225, 0.064), (1.02060862, 0.068), 
			(1.021028729, 0.072), (1.021426276, 0.076), (1.021803571, 0.08), (1.022162582, 0.084), (1.022505007, 0.088), (1.022832315, 0.092), (1.023145789, 0.096), 
			(1.023446554, 0.1), (1.023735605, 0.104), (1.024013823, 0.108), (1.024281993, 0.112), (1.024540818, 0.116), (1.02479093, 0.12), (1.025032898, 0.124), 
			(1.025267238, 0.128), (1.025494417, 0.132), (1.025714862, 0.136), (1.025928961, 0.14), (1.026137072, 0.144), (1.026339521, 0.148), (1.026536609, 0.152), 
			(1.026728613, 0.156), (1.026915791, 0.16), (1.027098379, 0.164), (1.027276598, 0.168), (1.027450654, 0.172), (1.027620736, 0.176), (1.027787023, 0.18), 
			(1.027949682, 0.184), (1.028108867, 0.188), (1.028264724, 0.192), (1.02841739, 0.196), (1.028566994, 0.2)), type=YIELD_RATIO)
		DC04_d0.plastic.Potential(table=((1.0, 1, 1.187434209, 1.072127101, 1, 1.0), ))	
			
		#DC04_d90

		#has to be defined

		
		############################################### 2.Create Sketch Biegeblech ###########################################################################
		Blech_Sketch = Simu_3D.ConstrainedSketch(name='Biegeblech', sheetSize=200.0)
			# Construction lines
		Blech_Sketch.ConstructionLine(point1=(0.0, 0.0), point2=(0.0, -13.75))
		Blech_Sketch.VerticalConstraint(addUndoState=False, entity=Blech_Sketch.geometry[2])
		Blech_Sketch.ConstructionLine(point1=(0.0, 0.0), point2=(16.25, 0.0))
		Blech_Sketch.HorizontalConstraint(addUndoState=False, entity=Blech_Sketch.geometry[3])
		Blech_Sketch.FixedConstraint(entity=Blech_Sketch.geometry[3])
		Blech_Sketch.FixedConstraint(entity=Blech_Sketch.geometry[2])
			# Biegeblech
		Blech_Sketch.rectangle(point1=(-20, 8.75), point2=(40.0, 0.0))
		Blech_Sketch.CoincidentConstraint(addUndoState=False, entity1=Blech_Sketch.vertices[2], entity2=Blech_Sketch.geometry[3])
		Blech_Sketch.DistanceDimension(entity1=Blech_Sketch.geometry[7], entity2=Blech_Sketch.geometry[3], textPoint=(-30.0, 4.5), value=t)
		Blech_Sketch.DistanceDimension(entity1=Blech_Sketch.geometry[4], entity2= Blech_Sketch.geometry[6], textPoint=(-5.5, 7.5), value=40.0)		#Biegeblechlaenge				
		Blech_Sketch.Spot(point=(3.8, 0.0))
		Blech_Sketch.CoincidentConstraint(addUndoState=False, entity1=Blech_Sketch.vertices[4], entity2=Blech_Sketch.geometry[3])
		Blech_Sketch.DistanceDimension(entity1=Blech_Sketch.vertices[4], entity2=Blech_Sketch.geometry[6], textPoint=(4.5, -1.8), value=15)											
		Blech_Sketch.Spot(point=(5.7, 0.0))
		Blech_Sketch.CoincidentConstraint(addUndoState=False, entity1=Blech_Sketch.vertices[5], entity2=Blech_Sketch.geometry[3])
		Blech_Sketch.HorizontalDimension(textPoint=(1.0, 2.5), value=7.0, vertex1=Blech_Sketch.vertices[5], vertex2=Blech_Sketch.vertices[4])
		Blech_Sketch.CoincidentConstraint(entity1=Blech_Sketch.vertices[4], entity2=Blech_Sketch.geometry[2])

		############################################### 3.Create sketch Niederhalter ##########################################################################
		Niederhalter_Sketch = Simu_3D.ConstrainedSketch(name='Niederhalter', sheetSize=200.0)
		Niederhalter_Sketch.sketchOptions.setValues( gridOrigin=(0, 0))
		Niederhalter_Sketch.retrieveSketch(sketch=Blech_Sketch)
		Niederhalter_Sketch.setAsConstruction(objectList=(Niederhalter_Sketch.geometry[6], Niederhalter_Sketch.geometry[7], Niederhalter_Sketch.geometry[8], Niederhalter_Sketch.geometry[9]))
		Niederhalter_Sketch.rectangle(point1=(-5.5, 6.0), point2=(12.0, 3.0))
		Niederhalter_Sketch.DistanceDimension(entity1= Niederhalter_Sketch.geometry[10], entity2=Niederhalter_Sketch.geometry[12],
			textPoint=(2.5, 9.0), value=19.5)					#19.5 [mm] = Length Niederhalter
		Niederhalter_Sketch.DistanceDimension(entity1=Niederhalter_Sketch.geometry[13], entity2=Niederhalter_Sketch.geometry[11], 
			textPoint=(17.0, 4.5), value=2.5)					#2.5 [mm] = height Niederhalter
		Niederhalter_Sketch.Spot(point=(7.5,3.0))
		Niederhalter_Sketch.CoincidentConstraint(addUndoState=False, entity1=Niederhalter_Sketch.vertices[10], entity2=Niederhalter_Sketch.geometry[11])
		Niederhalter_Sketch.Spot(point=(1.0, 2.85))
		Niederhalter_Sketch.CoincidentConstraint(addUndoState=False, entity1=Niederhalter_Sketch.vertices[11], entity2=Niederhalter_Sketch.geometry[11])
		Niederhalter_Sketch.HorizontalDimension(textPoint=(4.0, 4.0), value=7.0, vertex1=Niederhalter_Sketch.vertices[11], 
			vertex2=Niederhalter_Sketch.vertices[10])	# distance Sucherlöcher Niederhalter
		Niederhalter_Sketch.CoincidentConstraint(entity1=Niederhalter_Sketch.vertices[11], entity2=Niederhalter_Sketch.geometry[4])
		Niederhalter_Sketch.DistanceDimension(entity1=Niederhalter_Sketch.geometry[11], entity2=Niederhalter_Sketch.geometry[5], 
			textPoint=(-8.0, 1.0), value=t)				# distance Sucherloch zu Biegeebene (Blechdicke t)
		Niederhalter_Sketch.DistanceDimension(entity1=Niederhalter_Sketch.vertices[11], entity2=Niederhalter_Sketch.geometry[10], 
			textPoint=(-3.5, 5.5), value=4.6)			# distance Sucherloch zu CenterLine

		################################################ 4.Create sketch Biegekern #############################################################################
		Biegekern_Sketch = Simu_3D.ConstrainedSketch(name='Biegekern', sheetSize= 200.0)
		Biegekern_Sketch.sketchOptions.setValues(gridOrigin=(0, 0))
		Biegekern_Sketch.retrieveSketch(sketch=Blech_Sketch)
		Biegekern_Sketch.setAsConstruction(objectList=(Biegekern_Sketch.geometry[6], Biegekern_Sketch.geometry[7], Biegekern_Sketch.geometry[8], Biegekern_Sketch.geometry[9]))
		Biegekern_Sketch.Line(point1=( 29.0, 0.0), point2=(1.9, 0.0))
		Biegekern_Sketch.HorizontalConstraint(addUndoState=False, entity=Biegekern_Sketch.geometry[10])
		Biegekern_Sketch.ParallelConstraint(addUndoState=False, entity1=Biegekern_Sketch.geometry[5], entity2=Biegekern_Sketch.geometry[10])
		Biegekern_Sketch.CoincidentConstraint(addUndoState=False, entity1=Biegekern_Sketch.vertices[6], entity2=Biegekern_Sketch.geometry[5])
		Biegekern_Sketch.CoincidentConstraint(addUndoState=False, entity1=Biegekern_Sketch.vertices[7], entity2=Biegekern_Sketch.geometry[5])
		Biegekern_Sketch.Line(point1=(30.0, 0.0), point2=(30.0, -27.0))
		Biegekern_Sketch.VerticalConstraint(addUndoState=False, entity=Biegekern_Sketch.geometry[11])
		Biegekern_Sketch.PerpendicularConstraint(addUndoState=False, entity1=Biegekern_Sketch.geometry[10], entity2=Biegekern_Sketch.geometry[11])
		Biegekern_Sketch.Line(point1=( 30.0, -27.0), point2=(4.5, -27.0))
		Biegekern_Sketch.HorizontalConstraint(addUndoState=False, entity=Biegekern_Sketch.geometry[12])
		Biegekern_Sketch.PerpendicularConstraint(addUndoState=False, entity1=Biegekern_Sketch.geometry[11], entity2=Biegekern_Sketch.geometry[12])
		Biegekern_Sketch.Line(point1=(4.5, -27.0), point2=(1.9, 0.0))
		Biegekern_Sketch.FilletByRadius(curve1=Biegekern_Sketch.geometry[13], curve2=Biegekern_Sketch.geometry[10], nearPoint1=(2.2, -3.3), 
			nearPoint2=(3.1, 0.15), radius=rk)			#Biegekernradius rk
		Biegekern_Sketch.Spot(point=(23.0, 0.0))
		Biegekern_Sketch.CoincidentConstraint(addUndoState=False, entity1=Biegekern_Sketch.vertices[15], entity2=Biegekern_Sketch.geometry[5])
		Biegekern_Sketch.Spot(point=(18.5, 0.0))
		Biegekern_Sketch.CoincidentConstraint(addUndoState=False, entity1=Biegekern_Sketch.vertices[16], entity2=Biegekern_Sketch.geometry[5])
		Biegekern_Sketch.CoincidentConstraint(entity1=Biegekern_Sketch.vertices[8], entity2=Biegekern_Sketch.vertices[6])
		Biegekern_Sketch.DistanceDimension(entity1=Biegekern_Sketch.geometry[5], entity2=Biegekern_Sketch.geometry[12], textPoint=(37.0, -14.1), value=25.0)
		Biegekern_Sketch.AngularDimension(line1=Biegekern_Sketch.geometry[13], line2=Biegekern_Sketch.geometry[5], textPoint=(11.1, -9.5), value=80.0)
		Biegekern_Sketch.DistanceDimension(entity1=Biegekern_Sketch.vertices[13], entity2=Biegekern_Sketch.geometry[11], textPoint=(10.7224779129028, 4.23918437957764), value=25.0-rk)
		Biegekern_Sketch.DistanceDimension(entity1=Biegekern_Sketch.geometry[4], entity2=Biegekern_Sketch.geometry[11], textPoint=(9.67, -12.39), value=19.9)
		Biegekern_Sketch.CoincidentConstraint(entity1=Biegekern_Sketch.geometry[4], entity2=Biegekern_Sketch.vertices[16])

		################################################## 5.Create sketch Biegeform_2 ###########################################################################
		Biegeform2_Sketch = Simu_3D.ConstrainedSketch(name='Biegeform_2', sheetSize=300.0)
		Biegeform2_Sketch.sketchOptions.setValues(gridOrigin=(0, 0))
		Biegeform2_Sketch.retrieveSketch(sketch=Biegekern_Sketch)
		Biegeform2_Sketch.setAsConstruction(objectList=(Biegeform2_Sketch.geometry[12], Biegeform2_Sketch.geometry[13], Biegeform2_Sketch.geometry[14], 
			Biegeform2_Sketch.geometry[15], Biegeform2_Sketch.geometry[16]))
		Biegeform2_Sketch.ConstructionLine(point1=(-7.24, 23.53), point2=(-7.24, 28.0549999999674))
		Biegeform2_Sketch.VerticalConstraint(addUndoState=False, entity=Biegeform2_Sketch.geometry[17])
		Biegeform2_Sketch.Line(point1=(-7.24, 12.67), point2=(-7.24, 34.6614646911621))
		Biegeform2_Sketch.VerticalConstraint(addUndoState=False, entity=Biegeform2_Sketch.geometry[18])
		Biegeform2_Sketch.ParallelConstraint(addUndoState=False, entity1=Biegeform2_Sketch.geometry[17], entity2=Biegeform2_Sketch.geometry[18])
		Biegeform2_Sketch.CoincidentConstraint(addUndoState=False, entity1=Biegeform2_Sketch.vertices[17], entity2=Biegeform2_Sketch.geometry[17])
		Biegeform2_Sketch.Line(point1=(-7.24, 34.6614646911621), point2=(-19.0049999999441, 34.6614646911621))
		Biegeform2_Sketch.HorizontalConstraint(addUndoState=False, entity=Biegeform2_Sketch.geometry[19])
		Biegeform2_Sketch.PerpendicularConstraint(addUndoState=False, entity1=Biegeform2_Sketch.geometry[18], entity2=Biegeform2_Sketch.geometry[19])
		Biegeform2_Sketch.Line(point1=(-19.0049999999441, 34.6614646911621), point2=(-19.0049999999441, 3.62))
		Biegeform2_Sketch.VerticalConstraint(addUndoState=False, entity=Biegeform2_Sketch.geometry[20])
		Biegeform2_Sketch.PerpendicularConstraint(addUndoState=False, entity1=Biegeform2_Sketch.geometry[19], entity2=Biegeform2_Sketch.geometry[20])
		Biegeform2_Sketch.Line(point1=(-19.0049999999441, 3.62), point2=(-10.86, 3.62))
		Biegeform2_Sketch.HorizontalConstraint(addUndoState=False, entity=Biegeform2_Sketch.geometry[21])
		Biegeform2_Sketch.PerpendicularConstraint(addUndoState=False, entity1=Biegeform2_Sketch.geometry[20], entity2=Biegeform2_Sketch.geometry[21])
			
		#Eingabe der Teilkreise für den Analytical Rigid Body Biegeform 2 (Spline nicht möglich!!)

		Biegeform2_Sketch.ArcByCenterEnds(center=(-14.2828578948975, 8.12394332885742), direction=COUNTERCLOCKWISE, point1=(-10.86, 3.62), point2=(-10.1229019165039, 4.24623107910156))
		Biegeform2_Sketch.VerticalDimension(textPoint=(-9.01970958709717, 3.90205526351929), value=0.07, vertex1=Biegeform2_Sketch.vertices[22], vertex2= Biegeform2_Sketch.vertices[21])
		Biegeform2_Sketch.HorizontalDimension(textPoint=(-10.1520080566406, 4.53326463699341), value=0.3, vertex1= Biegeform2_Sketch.vertices[22], vertex2=Biegeform2_Sketch.vertices[21])
		Biegeform2_Sketch.ArcByCenterEnds(center=(-11.4033117294312, 6.78747701644897), direction=COUNTERCLOCKWISE, point1=(-9.91122834294817, 4.26669308601393), point2=(-9.37093639373779, 4.64011573791504))
		Biegeform2_Sketch.VerticalDimension(textPoint=(-7.7111644744873, 4.42030715942383), value=0.2, vertex1=Biegeform2_Sketch.vertices[22], vertex2=Biegeform2_Sketch.vertices[24])
		Biegeform2_Sketch.HorizontalDimension(textPoint=(-9.7004337310791, 3.44227695465088), value=0.28, vertex1=Biegeform2_Sketch.vertices[22], vertex2=Biegeform2_Sketch.vertices[24])
		Biegeform2_Sketch.ArcByCenterEnds(center=(-10.8104095458984, 6.55671501159668), direction=COUNTERCLOCKWISE, point1=(-9.63122834294817, 4.46669308601393), point2=(-9.38496780395508, 4.78370189666748))
		Biegeform2_Sketch.VerticalDimension(textPoint=(-7.60900735855103, 4.58540439605713), value=0.29, vertex1=Biegeform2_Sketch.vertices[24], vertex2=Biegeform2_Sketch.vertices[26])
		Biegeform2_Sketch.HorizontalDimension(textPoint=(-9.36159896850586, 3.55892252922058), value=0.22, vertex1=Biegeform2_Sketch.vertices[24], vertex2=Biegeform2_Sketch.vertices[26])
		Biegeform2_Sketch.ArcByCenterEnds(center=(-12.0839595794678, 6.58004426956177), direction=COUNTERCLOCKWISE, point1=(-9.41122834294817, 4.75669308601393), point2=(-9.17465591430664, 5.1569676399231))
		Biegeform2_Sketch.VerticalDimension(textPoint=( -7.16501712799072, 4.86535358428955), value=0.32, vertex1=Biegeform2_Sketch.vertices[26], vertex2=Biegeform2_Sketch.vertices[28])
		Biegeform2_Sketch.HorizontalDimension(textPoint=(-9.33823204040527, 3.13899827003479), value=0.17, vertex1=Biegeform2_Sketch.vertices[26], vertex2=Biegeform2_Sketch.vertices[28])
		Biegeform2_Sketch.ArcByCenterEnds(center=(-12.1540632247925, 6.97663927078247), direction=COUNTERCLOCKWISE, point1=(-9.24122834294817, 5.07669308601393), point2=(-9.05781650543213, 5.55356311798096))
		Biegeform2_Sketch.VerticalDimension(textPoint=(-7.85436964035034, 5.30860710144043), value=0.34, vertex1=Biegeform2_Sketch.vertices[30], vertex2=Biegeform2_Sketch.vertices[28])
		Biegeform2_Sketch.HorizontalDimension(textPoint=(-9.11623573303223, 3.48893523216248), value=0.14, vertex1=Biegeform2_Sketch.vertices[28], vertex2=Biegeform2_Sketch.vertices[30])
		Biegeform2_Sketch.ArcByCenterEnds(center=(-11.4763946533203, 7.13994312286377), direction=COUNTERCLOCKWISE, point1=(-9.10122834294816, 5.41669308601393), point2=(-8.97602939605713, 5.81018304824829))
		Biegeform2_Sketch.VerticalDimension(textPoint=(-8.14646911621094, 5.7168664932251), value=0.35, vertex1=Biegeform2_Sketch.vertices[30], vertex2=Biegeform2_Sketch.vertices[32])
		Biegeform2_Sketch.HorizontalDimension(textPoint=(-8.95266151428223, 6.32342386245728), value=0.11, vertex1=Biegeform2_Sketch.vertices[30], vertex2=Biegeform2_Sketch.vertices[32])
		Biegeform2_Sketch.ArcByCenterEnds(center=(-12.0839595794678, 6.8949875831604), direction=COUNTERCLOCKWISE, point1=(-8.99122834294817, 5.76669308601393), point2=(-8.8358211517334, 6.17178440093994))
		Biegeform2_Sketch.VerticalDimension(textPoint=(-8.01794528961182, 6.06680345535278), value=0.71, vertex1=Biegeform2_Sketch.vertices[34], vertex2=Biegeform2_Sketch.vertices[32])
		Biegeform2_Sketch.HorizontalDimension(textPoint=(-8.8825569152832, 6.8949875831604), value=0.15, vertex1=Biegeform2_Sketch.vertices[34], vertex2=Biegeform2_Sketch.vertices[32])
		Biegeform2_Sketch.ArcByCenterEnds(center=(-12.1190118789673, 7.17493677139282), direction=COUNTERCLOCKWISE, point1=(-8.84122834294817, 6.47669308601393), point2=(-8.7774019241333, 7.23325967788696))
		Biegeform2_Sketch.VerticalDimension(textPoint=(-8.01794528961182, 6.85999393463135), value=0.81, vertex1=Biegeform2_Sketch.vertices[36], vertex2=Biegeform2_Sketch.vertices[34])
		Biegeform2_Sketch.HorizontalDimension(textPoint=(-8.82413768768311, 7.74650049209595), value=0.06, vertex1=Biegeform2_Sketch.vertices[34], vertex2=Biegeform2_Sketch.vertices[36])
			
		#Teilkreis Radiusangabe

		Biegeform2_Sketch.RadialDimension(curve=Biegeform2_Sketch.geometry[22], radius=0.7, textPoint=(-10.3518867492676, 5.31571626663208))
		Biegeform2_Sketch.RadialDimension(curve=Biegeform2_Sketch.geometry[23], radius=1.0, textPoint=(-10.038444519043, 4.83036851882935))
		Biegeform2_Sketch.RadialDimension(curve=Biegeform2_Sketch.geometry[24], radius=1.8, textPoint=(-10.0256509780884, 5.05388402938843))
		Biegeform2_Sketch.RadialDimension(curve=Biegeform2_Sketch.geometry[25], radius=2.9, textPoint=(-9.76338386535645, 5.20076560974121))
		Biegeform2_Sketch.RadialDimension(curve=Biegeform2_Sketch.geometry[26], radius=3.5, textPoint=(-9.64184474945068, 5.44982528686523))
		Biegeform2_Sketch.RadialDimension(curve=Biegeform2_Sketch.geometry[27], radius=5.1, textPoint=(-9.60986137390137, 5.80106401443481))
		Biegeform2_Sketch.RadialDimension(curve=Biegeform2_Sketch.geometry[28], radius=6.3, textPoint=(-9.6162576675415, 6.33750057220459))
		Biegeform2_Sketch.RadialDimension(curve=Biegeform2_Sketch.geometry[29], radius=5.4, textPoint=(-9.66103553771973, 6.93779945373535))
		Biegeform2_Sketch.CoincidentConstraint(entity1=Biegeform2_Sketch.vertices[17], entity2=Biegeform2_Sketch.vertices[36])
		Biegeform2_Sketch.DistanceDimension(entity1=Biegeform2_Sketch.geometry[20], entity2=Biegeform2_Sketch.geometry[17], textPoint=(-14.0524339675903, 38.599365234375), 
			value=15.0)						#Länge Biegeform in Vorschubsrichtung
		Biegeform2_Sketch.DistanceDimension(entity1=Biegeform2_Sketch.geometry[19], entity2=Biegeform2_Sketch.geometry[21], textPoint=(-28.0429954528809, 17.1948585510254), 
			value=40.0)						#Höhe Biegeform
		Biegeform2_Sketch.DistanceDimension(entity1=Biegeform2_Sketch.geometry[21], entity2=Biegeform2_Sketch.geometry[7], textPoint=(-34.0164451599121, 1.98553276062012), 
			value=t)						#Abstand zu Unterseite Blech (Blechdicke)

		#'Positionierung' Biegeform mit Niederhalterspalt	

		Biegeform2_Sketch.ConstructionLine(point1=(-5, 21), point2=(-5, 16))
		Biegeform2_Sketch.VerticalConstraint(addUndoState=False, entity=Biegeform2_Sketch.geometry[30])
		Biegeform2_Sketch.DistanceDimension(entity1=Biegeform2_Sketch.geometry[30], entity2=Biegeform2_Sketch.vertices[12], textPoint=(-5, -4), value=rk)	
		Biegeform2_Sketch.DistanceDimension(entity1=Biegeform2_Sketch.geometry[17], entity2=Biegeform2_Sketch.geometry[30], textPoint=(-4.8831033706665, 47.3591232299805), 
			value=1.13)						#--------------Niederhalterspalt-----------------#

		############################################## 6.Create Sketch final position as a function of angle #########################################################
		Positionierung_Sketch = Simu_3D.ConstrainedSketch(name='Positionierung', sheetSize=200.0)
		Positionierung_Sketch.sketchOptions.setValues(gridOrigin=(0, 0))
		Positionierung_Sketch.retrieveSketch(sketch=Biegeform2_Sketch)
		Positionierung_Sketch.setAsConstruction(objectList=(Positionierung_Sketch.geometry[20], Positionierung_Sketch.geometry[21], Positionierung_Sketch.geometry[22], 
			Positionierung_Sketch.geometry[23], Positionierung_Sketch.geometry[24], Positionierung_Sketch.geometry[25], Positionierung_Sketch.geometry[26], 
			Positionierung_Sketch.geometry[27], Positionierung_Sketch.geometry[28], Positionierung_Sketch.geometry[29], Positionierung_Sketch.geometry[30], 
			Positionierung_Sketch.geometry[31]))
		Positionierung_Sketch.ConstructionLine(point1=(-4.1, -1.0), point2=(-5.00623700492539, -0.577229978707188))
		Positionierung_Sketch.Line( point1=(-5.65106912632541, 0.262615698307409), point2=(-8.71110765567573, -3.4965101343405))
		Positionierung_Sketch.Line(point1=(-8.71110765567573, -3.4965101343405), point2=(-7.81655766250333, -4.22470002935734))
		Positionierung_Sketch.Line(point1=(-7.81655766250333, -4.22470002935734), point2=(-4.87553322185471, -0.36869324271029))
		#Positionierung_Sketch.AngularDimension(line1=Positionierung_Sketch.geometry[13], line2=Positionierung_Sketch.geometry[36], textPoint=(-8.937, -1.3126), 
			#value=alpha)					#Biegewinkel
		Positionierung_Sketch.Line(point1=(-13.5025, -1.42153031604876), point2=(-13.5025, 1.0))
		Positionierung_Sketch.CoincidentConstraint(addUndoState=False, entity1=Positionierung_Sketch.vertices[43], entity2=Positionierung_Sketch.geometry[23])
		Positionierung_Sketch.dragEntity(entity=Positionierung_Sketch.geometry[34]
			, points=((-10.3563370453384, -0.415478326172973), (-10.3805198669434, -0.389714181423187), (-8.7900390625, -1.42640423774719), (-8.2511157989502, -1.67573475837708), (-8.17224884033203, -1.7151026725769)))
		Positionierung_Sketch.ParallelConstraint(entity1=Positionierung_Sketch.geometry[36], entity2=Positionierung_Sketch.geometry[34])
		Positionierung_Sketch.DistanceDimension(entity1=Positionierung_Sketch.geometry[36], entity2=Positionierung_Sketch.geometry[34], textPoint=(-10.1439189910889, -3.05361366271973), value=t)
		Positionierung_Sketch.CoincidentConstraint(entity1=Positionierung_Sketch.vertices[38], entity2=Positionierung_Sketch.geometry[33])
		Positionierung_Sketch.CoincidentConstraint(entity1=Positionierung_Sketch.vertices[41], entity2=Positionierung_Sketch.geometry[33])
		Positionierung_Sketch.CoincidentConstraint(entity1=Positionierung_Sketch.geometry[33], entity2=Positionierung_Sketch.vertices[12])
		Positionierung_Sketch.AngularDimension(line1=Positionierung_Sketch.geometry[8], line2=Positionierung_Sketch.geometry[33], textPoint=(-1.60257291793823, -1.70021998882294), value=alpha)
		Positionierung_Sketch.CoincidentConstraint(entity1=Positionierung_Sketch.vertices[38], entity2=Positionierung_Sketch.geometry[18])
		Positionierung_Sketch.ParallelConstraint( entity1= Positionierung_Sketch.geometry[35], entity2=Positionierung_Sketch.geometry[33])
		Positionierung_Sketch.DistanceDimension( entity1=Positionierung_Sketch.geometry[35], entity2=Positionierung_Sketch.geometry[33], textPoint=(-6.56813097000122, -3.60456895828247), value=7.0)
		Positionierung_Sketch.CoincidentConstraint(entity1=Positionierung_Sketch.vertices[42], entity2=Positionierung_Sketch.geometry[13])
		Positionierung_Sketch.DistanceDimension(entity1=Positionierung_Sketch.geometry[37], entity2=Positionierung_Sketch.geometry[8], textPoint=(-6.61416435241699, 4.37695074081421), value=15.0)
		Positionierung_Sketch.PerpendicularConstraint(entity1=Positionierung_Sketch.geometry[34], entity2=Positionierung_Sketch.geometry[33])
		Positionierung_Sketch.setAsConstruction(objectList=(Positionierung_Sketch.geometry[36], Positionierung_Sketch.geometry[34], Positionierung_Sketch.geometry[35]))
		#Delete dimension for position in sketch Biegeform
		Positionierung_Sketch.delete(objectList=(Positionierung_Sketch.dimensions[35], ))
		
		#Abfrage welcher Tangentialbereich zutrifft
		if 0.0<=alpha<=25.6:
			Positionierung_Sketch.TangentConstraint(entity1=Positionierung_Sketch.geometry[24], entity2=Positionierung_Sketch.geometry[36])
		elif 25.6<alpha<=45.3:
			Positionierung_Sketch.TangentConstraint(entity1=Positionierung_Sketch.geometry[25], entity2=Positionierung_Sketch.geometry[36])
		elif 45.0<alpha<=58.6:
			Positionierung_Sketch.TangentConstraint(entity1=Positionierung_Sketch.geometry[26], entity2=Positionierung_Sketch.geometry[36])			
		elif 58.6<alpha<=65.6:
			Positionierung_Sketch.TangentConstraint(entity1=Positionierung_Sketch.geometry[27], entity2=Positionierung_Sketch.geometry[36])		
		elif 65.6<alpha<=70.6:
			Positionierung_Sketch.TangentConstraint(entity1=Positionierung_Sketch.geometry[28], entity2=Positionierung_Sketch.geometry[36])
		elif 70.6<alpha<=74.6:
			Positionierung_Sketch.TangentConstraint(entity1=Positionierung_Sketch.geometry[29], entity2=Positionierung_Sketch.geometry[36])		
		elif 74.6<alpha<=81.4:
			Positionierung_Sketch.TangentConstraint(entity1=Positionierung_Sketch.geometry[30], entity2=Positionierung_Sketch.geometry[36])		
		elif 81.4<alpha<90:
			Positionierung_Sketch.TangentConstraint(entity1=Positionierung_Sketch.geometry[31], entity2=Positionierung_Sketch.geometry[36])		
		elif alpha==90.0:
			Weg = 12
			#Positionierung_Sketch.delete(objectList=(Positionierung_Sketch.geometry[39], ))
		else:
			print('Falsche Eingabe')
		
		################################################ 7.Create part Analytical Rigid Body endposition #######################################################
		if alpha <90.0:
				#3D Shell Sketch
			Pos3D_Sketch = Simu_3D.ConstrainedSketch(name='Positionierung3D', sheetSize=200.0)
			Pos3D_Sketch.sketchOptions.setValues(gridOrigin=(0, 0))
			Pos3D_Sketch.retrieveSketch(sketch=Positionierung_Sketch)
				#3D Shell part
			Pos3D_part = Simu_3D.Part(dimensionality=THREE_D, name='Positionierung', type=ANALYTIC_RIGID_SURFACE)
			Pos3D_part.AnalyticRigidSurfExtrude(depth=0.01, sketch=Pos3D_Sketch)
			del Pos3D_Sketch			
				#Instance and Surface, Set
			Simu_3D.rootAssembly.DatumCsysByDefault(CARTESIAN)
			Endpos_instance = Simu_3D.rootAssembly.Instance(dependent=ON, name='Positionierung-1', part=Pos3D_part)
			Pos3D_part.Surface(name='Positionierung', side2Faces=Pos3D_part.faces.getSequenceFromMask(('[#1 ]', ), ))	
				#Reference point
			#Endpos_RP = Pos3D_part.ReferencePoint(point=Pos3D_part.InterestingPoint(Pos3D_part.edges[2], MIDDLE))
			Inertia_RP = Simu_3D.parts['Positionierung'].ReferencePoint(point=Simu_3D.parts['Positionierung'].InterestingPoint(Simu_3D.parts['Positionierung'].edges[2], MIDDLE))
			Inerita_RP_id = Inertia_RP.id
			Inertia_RP_set = Simu_3D.parts['Positionierung'].Set(name='Inertia_RP_set', referencePoints=(Simu_3D.parts['Positionierung'].referencePoints[Inerita_RP_id], ))
			Endpos_RP = Simu_3D.rootAssembly.ReferencePoint(point=Simu_3D.rootAssembly.instances['Positionierung-1'].vertices[2])
			Endpos_RP_set = Simu_3D.rootAssembly.Set(name='Endpos_RP_set', referencePoints=(Simu_3D.rootAssembly.referencePoints[4], ))
			#vector = ('%g,%g,%g' % (Endpos_RP.coordinates[0], Endpos_RP.coordinates[1], Endpos_RP.coordinates[2])) 
					
			
			RP_dis_ID = Endpos_RP.id
			RP_dis_ID_str = str(RP_dis_ID)
			#testxyz = Simu_3D.rootAssembly.features['RP-3']
			#test1 = Simu_3D.rootAssembly.referencePoints.values(RP_dis_ID_str)
			#test2 = Simu_3D.rootAssembly.features.keys(RP_dis_ID_str)
			#test3 = Simu_3D.rootAssembly.referencePoints.keys()[0]
			#test4 = Simu_3D.rootAssembly.referencePoints[RP_dis_ID]
			#test5 = Simu_3D.rootAssembly.referencePoints.getCoordinates(RP_dis_ID_str)
			
			#test = (test1,test2,test3,test4)
			
		#################################################### 8.Create part biegeform_2 ##########################################################################
			#Sketch
		Biegeform2_3D_Sketch = Simu_3D.ConstrainedSketch(name='Biegeform_2_3D', sheetSize=200.0)
		Biegeform2_3D_Sketch.sketchOptions.setValues(gridOrigin=(0, 0))
		Biegeform2_3D_Sketch.retrieveSketch(sketch=Biegeform2_Sketch)
			#3D Shell
		Biegeform2_3D_part = Simu_3D.Part(dimensionality=THREE_D, name='Biegeform_2', type=ANALYTIC_RIGID_SURFACE)
		Biegeform2_3D_part.AnalyticRigidSurfExtrude(depth=b, sketch=Biegeform2_3D_Sketch)	
		del Biegeform2_3D_Sketch	
			
		#################################################### 9.Create part Biegeblech ############################################################################
		Blech_part_sketch = Simu_3D.ConstrainedSketch(name='Biegeblech_3D', sheetSize=100.0)
		Blech_part_sketch.sketchOptions.setValues(gridOrigin=(0, 0))
		Blech_part_sketch.retrieveSketch(sketch=Blech_Sketch)
		Blech_part = Simu_3D.Part(dimensionality=THREE_D, name='Biegeblech', type=DEFORMABLE_BODY)
		Blech_part.BaseSolidExtrude(depth=b, sketch=Blech_part_sketch)
		Blech_part_sketch2 = Simu_3D.ConstrainedSketch(gridSpacing=2.5, name='Biegeblech_3D', sheetSize=100.0, transform=Blech_part.MakeSketchTransform
			(sketchPlane=Blech_part.faces[3], sketchPlaneSide=SIDE1, sketchUpEdge=Blech_part.edges[3], sketchOrientation=RIGHT, origin=(-5.0, 1.0, 15.0)))
		Blech_part.projectReferencesOntoSketch(filter=COPLANAR_EDGES, sketch=Blech_part_sketch2)
		Blech_part_sketch2.CircleByCenterPerimeter(center=(12.5, 0.0), point1=(13.75, 1.25))
		Blech_part_sketch2.CircleByCenterPerimeter(center=(3.75, 0.0), point1=(5.0, 1.25))
		Blech_part_sketch2.RadialDimension(curve=Blech_part_sketch2.geometry[7], radius=1.5, textPoint=(0.489141941070557, -1.13224792480469))
		Blech_part_sketch2.RadialDimension(curve=Blech_part_sketch2.geometry[6], radius=1.5, textPoint=(9.48728609085083, -2.21023178100586))
		Blech_part_sketch2.HorizontalDimension(textPoint=(8.34752321243286, 4.01812839508057), value=7.0, vertex1=Blech_part_sketch2.vertices[4], vertex2=Blech_part_sketch2.vertices[6])
		Blech_part_sketch2.DistanceDimension(entity1=Blech_part_sketch2.geometry[2], entity2=Blech_part_sketch2.vertices[4], textPoint=(14.2862968444824, 4.01812839508057), value=8.0)
		Blech_part_sketch2.DistanceDimension(entity1=Blech_part_sketch2.vertices[4], entity2=Blech_part_sketch2.geometry[3], textPoint=(23.1044807434082, -9.09736442565918), value=b/2)
		Blech_part_sketch2.DistanceDimension(entity1=Blech_part_sketch2.vertices[6], entity2=Blech_part_sketch2.geometry[3], textPoint=(-3.88995361328125, -9.45669174194336), value=b/2)
		Blech_part.CutExtrude(flipExtrudeDirection=OFF, sketch=Blech_part_sketch2, sketchOrientation=RIGHT, sketchPlane=Blech_part.faces[3], sketchPlaneSide=SIDE1, sketchUpEdge=Blech_part.edges[3])
		del Blech_part_sketch2
				
			#create sketch partition Biegeblech
		Part_sketch = Simu_3D.ConstrainedSketch(gridSpacing=2.5, name='Biegebl_Part', sheetSize=100.0, transform=Blech_part.
			MakeSketchTransform(sketchPlane=Blech_part.faces[3], sketchPlaneSide=SIDE1, sketchUpEdge=Blech_part.edges[9], sketchOrientation=RIGHT, origin=(-5, 1.0, 15.0)))
		Blech_part.projectReferencesOntoSketch(filter=COPLANAR_EDGES, sketch=Part_sketch)
		Part_sketch.Line(point1=(1.875, 15.0), point2=(1.875, -15.0))
		Part_sketch.VerticalConstraint(addUndoState=False, entity=Part_sketch.geometry[8])
		Part_sketch.PerpendicularConstraint(addUndoState=False, entity1=Part_sketch.geometry[3], entity2=Part_sketch.geometry[8])
		Part_sketch.CoincidentConstraint(addUndoState=False, entity1=Part_sketch.vertices[8], entity2=Part_sketch.geometry[3])
		Part_sketch.CoincidentConstraint(addUndoState=False, entity1=Part_sketch.vertices[9], entity2=Part_sketch.geometry[2])
		Part_sketch.Line(point1=(-3.125, 15.0), point2=(-3.125, -15.0))
		Part_sketch.VerticalConstraint(addUndoState=False, entity=Part_sketch.geometry[9])
		Part_sketch.PerpendicularConstraint(addUndoState=False, entity1=Part_sketch.geometry[3], entity2=Part_sketch.geometry[9])
		Part_sketch.CoincidentConstraint(addUndoState=False, entity1=Part_sketch.vertices[10], entity2=Part_sketch.geometry[3])
		Part_sketch.CoincidentConstraint(addUndoState=False, entity1=Part_sketch.vertices[11], entity2=Part_sketch.geometry[2])
		Part_sketch.Line(point1=(-8.125, 15.0), point2=(-8.125, -15.0))
		Part_sketch.VerticalConstraint(addUndoState=False, entity=Part_sketch.geometry[10])
		Part_sketch.PerpendicularConstraint(addUndoState=False, entity1=Part_sketch.geometry[3], entity2=Part_sketch.geometry[10])
		Part_sketch.CoincidentConstraint(addUndoState=False, entity1=Part_sketch.vertices[12], entity2=Part_sketch.geometry[3])
		Part_sketch.CoincidentConstraint(addUndoState=False, entity1=Part_sketch.vertices[13], entity2=Part_sketch.geometry[2])
		Part_sketch.DistanceDimension(entity1=Part_sketch.vertices[6], entity2=Part_sketch.geometry[8], textPoint=(3.30954675982666, 4.71877098083496), 
			value=2.1)				#distance center inner Sucherloch to transition mesh
		Part_sketch.DistanceDimension(entity1=Part_sketch.geometry[8], entity2=Part_sketch.geometry[9], textPoint=(0.35892420123291, 4.81228637695313), 
			value=1.0)				#length transition mesh
		Part_sketch.DistanceDimension(entity1=Part_sketch.geometry[10], entity2=Part_sketch.geometry[9], textPoint=(-2.02967328717041, 4.76552963256836), 
			value=8.0)				#length bending area
		Blech_part.PartitionFaceBySketchThruAll(faces=Blech_part.faces.getSequenceFromMask(('[#ff ]', ), ), 
			sketch=Part_sketch, sketchPlane=Blech_part.faces[3], sketchPlaneSide=SIDE1, sketchUpEdge=Blech_part.edges[9])
		del Part_sketch		

			#create partition Biegeblech
		Blech_part.PartitionCellByPlanePointNormal(cells=Blech_part.cells.getSequenceFromMask(('[#1 ]', ), ), normal=Blech_part.edges[1], point=Blech_part.vertices[1])
		Blech_part.PartitionCellByPlanePointNormal(cells=Blech_part.cells.getSequenceFromMask(('[#1 ]', ), ), normal=Blech_part.edges[12], point=Blech_part.vertices[11])
		Blech_part.PartitionCellByPlanePointNormal(cells=Blech_part.cells.getSequenceFromMask(('[#1 ]', ), ), normal=Blech_part.edges[39], point=Blech_part.vertices[10])		

			#create sets and surfaces Biegeblech
				#create Sets for seed
		Blech_part.Set(edges=Blech_part.edges.getSequenceFromMask(('[#aa0 ]', ), ), name='Dicke_Seed')
		Seed_Breite = Blech_part.Set(edges=Blech_part.edges.getSequenceFromMask(('[#1b00555 #1e ]', ), ), name='Breite_Seed')
		Blech_part.Set(edges=Blech_part.edges.getSequenceFromMask(('[#c0000 #60 ]', ), ), name='Laenge_Biegebereich_Seed')
		Blech_part.Set(edges=Blech_part.edges.getSequenceFromMask(('[#82005000 ]', ), ), name='Laenge_Frei_Seed')
		Blech_part.Set(edges=Blech_part.edges.getSequenceFromMask(('[#10428000 ]', ), ), name='Laenge_Klemm_Seed')
		Blech_part.Set( edges=Blech_part.edges.getSequenceFromMask(('[#c000000 #81 ]', ), ), name='Laenge_uebergang_Seed')	
				#create Surfaces
		Blech_Surface_Oberseite = Blech_part.Surface(name='Blech_Oberseite', side1Faces=Blech_part.faces.getSequenceFromMask(('[#47000 ]', ), ))
		Blech_Surface_Unterseite = Blech_part.Surface(name='Blech_Unterseite', side1Faces=Blech_part.faces.getSequenceFromMask(('[#201c0 ]', ), ))
		Blech_Surface_Sucher = Blech_part.Surface(name='Sucherloecher', side1Faces=Blech_part.faces.getSequenceFromMask(('[#18000 ]', ), ))	
				#create Sets
		Sucherloecher_set = Blech_part.Set(faces=Blech_part.faces.getSequenceFromMask(('[#18000 ]', ), ), name='Sucherloecher')
		
			#create global orientation Biegeblech
		Blech_part.MaterialOrientation(additionalRotationType=ROTATION_NONE, axis=AXIS_1, fieldName='', localCsys=None, orientationType=GLOBAL,
			region=Region(cells=Blech_part.cells.getSequenceFromMask(mask=('[#f ]', ), )), stackDirection=STACK_3)
			
			#create Section Assigment Biegeblech
		Simu_3D.HomogeneousSolidSection(material=materialname, name='Biegeblech', thickness=None)
		Blech_part.SectionAssignment(offset=0.0, offsetField='', offsetType=MIDDLE_SURFACE, region=Region(cells=Blech_part.
			cells.getSequenceFromMask(mask=('[#f ]', ), )), sectionName='Biegeblech', thicknessAssignment=FROM_SECTION)	
		
			#create Controls (Mesh) Biegeblech
		Blech_part.setMeshControls(algorithm=ADVANCING_FRONT, regions=Blech_part.cells.getSequenceFromMask(('[#8 ]', ), ), technique=SWEEP)
		Blech_part.setMeshControls(algorithm=ADVANCING_FRONT, regions=Blech_part.cells.getSequenceFromMask(('[#2 ]', ), ), technique=SWEEP)
		Blech_part.setMeshControls(regions=Blech_part.cells.getSequenceFromMask(('[#4 ]', ), ), technique=STRUCTURED)
		Blech_part.setMeshControls(algorithm=ADVANCING_FRONT, regions=Blech_part.cells.getSequenceFromMask(('[#2 ]', ), ), technique=SWEEP)
			
			#create Element Type Biegeblech
		Blech_part.setElementType(elemTypes=(ElemType(elemCode=C3D8R, elemLibrary=EXPLICIT, secondOrderAccuracy=OFF, kinematicSplit=AVERAGE_STRAIN, hourglassControl=DEFAULT,
			distortionControl=DEFAULT), ElemType(elemCode=C3D6, elemLibrary=EXPLICIT), ElemType(elemCode=C3D4, elemLibrary=EXPLICIT)), regions=(Blech_part.cells.getSequenceFromMask(('[#8 ]', ), ), ))
		Blech_part.setElementType(elemTypes=(ElemType(elemCode=C3D8R, elemLibrary=EXPLICIT, secondOrderAccuracy=OFF, kinematicSplit=AVERAGE_STRAIN, hourglassControl=DEFAULT, 
			distortionControl=DEFAULT), ElemType(elemCode=C3D6, elemLibrary=EXPLICIT, secondOrderAccuracy=OFF, distortionControl=DEFAULT), ElemType(elemCode=C3D4, elemLibrary=EXPLICIT, secondOrderAccuracy=OFF, distortionControl=DEFAULT)), 
			regions=(Blech_part.cells.getSequenceFromMask(('[#4 ]', ), ), ))
		Blech_part.setElementType(elemTypes=(ElemType(elemCode=C3D8R, elemLibrary=EXPLICIT, secondOrderAccuracy=OFF, kinematicSplit=AVERAGE_STRAIN, hourglassControl=DEFAULT, 
			distortionControl=DEFAULT), ElemType(elemCode=C3D6, elemLibrary=EXPLICIT), ElemType(elemCode=C3D4, elemLibrary=EXPLICIT)), regions=(Blech_part.cells.getSequenceFromMask(('[#2 ]', ), ), ))
		Blech_part.setElementType(elemTypes=(ElemType(elemCode=C3D8R, elemLibrary=EXPLICIT, secondOrderAccuracy=OFF, kinematicSplit=AVERAGE_STRAIN, hourglassControl=DEFAULT, 
			distortionControl=DEFAULT), ElemType(elemCode=C3D6, elemLibrary=EXPLICIT),ElemType(elemCode=C3D4, elemLibrary=EXPLICIT)), regions=(Blech_part.cells.getSequenceFromMask(('[#1 ]', ), ), ))		
			
			#create Seed Biegeblech (Einteilung in Dickenrichtung funktioniert noch nicht!!!!!!!!)
		Blech_part.seedEdgeByNumber(constraint=FINER, edges=Blech_part.edges.getSequenceFromMask(('[#1b00555 #1e ]', ), ), number=20)						#Seed Breite Blech
		Blech_part.seedEdgeByNumber(constraint=FINER, edges=Blech_part.edges.getSequenceFromMask(('[#aa0 ]', ), ), number=10*int(t))						#Seed Dicke gesamt
		Blech_part.seedEdgeBySize(constraint=FINER, deviationFactor=0.1, edges=Blech_part.edges.getSequenceFromMask(('[#c000000 #81 ]', ), ), size=0.2)	#Seed Längsrichtung Übergangsbereich
		Blech_part.seedEdgeBySize(constraint=FINER, deviationFactor=0.1, edges=Blech_part.edges.getSequenceFromMask(('[#c0000 #60 ]', ), ), size=0.1)		#Seed Längsrichtung Biegezon
		number = Blech_part.seedEdgeBySize(constraint=FINER, deviationFactor=0.1, edges=Blech_part.edges.getSequenceFromMask(('[#82005000 ]', ), ), size=0.5)		#Seed Längsrichtung freie Biegerand
		Blech_part.seedEdgeBySize(constraint=FINER, deviationFactor=0.1, edges=Blech_part.edges.getSequenceFromMask(('[#10428000 ]', ), ), size=0.6)		#Seed Längsrichtung Klemmbereich
			
			#create Mesh Biegeblech
		Blech_part.generateMesh( )
		
		###################################################### 10.Create part Niederhalter ########################################################################
		Niederh3D_sketch = Simu_3D.ConstrainedSketch(name='Niederhalter3D', sheetSize= 200.0)
		Niederh3D_sketch.sketchOptions.setValues(gridOrigin=(0, 0))
		Niederh3D_sketch.retrieveSketch(sketch=Niederhalter_Sketch)
		Niederh3D_part = Simu_3D.Part(dimensionality=THREE_D, name='Niederhalter', type=ANALYTIC_RIGID_SURFACE)
		Niederh3D_part.AnalyticRigidSurfExtrude(depth=30.0, sketch=Niederh3D_sketch)
		del Niederh3D_sketch

		###################################################### 11.Create part Biegekern ############################################################################
		Biegekern3D_sketch = Simu_3D.ConstrainedSketch(name='Biegekern3D', sheetSize= 200.0)
		Biegekern3D_sketch.sketchOptions.setValues(gridOrigin=(0, 0))
		Biegekern3D_sketch.retrieveSketch(sketch=Biegekern_Sketch)
		Biegekern3D_part = Simu_3D.Part(dimensionality=THREE_D, name='Biegekern', type=ANALYTIC_RIGID_SURFACE)
		Biegekern3D_part.AnalyticRigidSurfExtrude(depth=b, sketch=Biegekern3D_sketch)
		del Biegekern3D_sketch

		###################################################### 12.Create instances #################################################################################
		Biegekern_Instance = Simu_3D.rootAssembly.Instance(dependent=ON, name='Biegekern-1', part=Biegekern3D_part)
		Niederh_Instance = Simu_3D.rootAssembly.Instance(dependent=ON, name='Niederhalter-1', part=Niederh3D_part)
		Biegeform2_Instance = Simu_3D.rootAssembly.Instance(dependent=ON, name='Biegeform_2', part=Biegeform2_3D_part)
		Blech_Instance = Simu_3D.rootAssembly.Instance(dependent=ON, name='Biegeblech-1', part=Blech_part)
		
		################################################# 13.Create position constraints ###########################################################################
		Simu_3D.rootAssembly.EdgeToEdge(fixedAxis=Niederh_Instance.edges[4], flip=ON, movableAxis=Blech_Instance.edges[28])
		

		################################################# 14.Create Reference Points ###############################################################################
		Biegeform2_3D_part.ReferencePoint(point=(-14.445, 41.0, 0.0))
		Niederh3D_part.ReferencePoint(point=(5.15, 1.0, 0.0))
		Biegekern3D_part.ReferencePoint(point=(7.9, 0.0, 0.0))
			#create sets of reference points
		Biegeform2_RP_set = Biegeform2_3D_part.Set(name='Biegeform_2_RP', referencePoints=(Biegeform2_3D_part.referencePoints[2], ))
		Niederh_RP_set = Niederh3D_part.Set(name='Niederhalter_RP', referencePoints=(Niederh3D_part.referencePoints[2], ))
		Biegekern_RP_set = Biegekern3D_part.Set(name='Biegekern', referencePoints=(Biegekern3D_part.referencePoints[2], ))
		
		#################################################### 15.Create inertias ####################################################################################
		Biegeform2_3D_part.engineeringFeatures.PointMassInertia(alpha=0.0, composite=0.0, mass=0.00032, name='Biegeform_2', region=Biegeform2_3D_part.sets['Biegeform_2_RP'])
		Niederh3D_part.engineeringFeatures.PointMassInertia(alpha=0.0, composite=0.0, mass=0.00052, name='Inertia_Niederhalter', region=Niederh3D_part.sets['Niederhalter_RP'])
		#Pos3D_part.engineeringFeatures.PointMassInertia(alpha=0.0, composite=0.0, mass=0.00052, name='Inertia_Endpos', region=Pos3D_part.sets['Inertia_RP_set'])
		
		####################################################### 16.Create surfaces #################################################################################
		Biegekern_Surface  = Biegekern3D_part.Surface(name='Biegekern', side1Faces=Biegekern3D_part.faces.getSequenceFromMask(( '[#1f ]', ), ))
		Biegeform_Surface = Biegeform2_3D_part.Surface(name='Biegeform_2', side1Faces=Biegeform2_3D_part.faces.getSequenceFromMask(('[#fff ]', ), ))
		Niederh_Surface = Niederh3D_part.Surface(name='Niederhalter', side1Faces=Niederh3D_part.faces.getSequenceFromMask(('[#f ]', ), ))

		########################################################17.Create constraints ##############################################################################
		mdb.models[Simu_name].RigidBody(name='Biegeform_2', refPointAtCOM=ON, refPointRegion=Region(referencePoints=(mdb.models[Simu_name].rootAssembly.instances['Biegeform_2'].
			referencePoints[2], )), surfaceRegion=mdb.models[Simu_name].rootAssembly.instances['Biegeform_2'].surfaces['Biegeform_2'])
		mdb.models[Simu_name].RigidBody(name='Niederhalter', refPointAtCOM=ON, refPointRegion=Region(referencePoints=(mdb.models[Simu_name].rootAssembly.instances['Niederhalter-1'].
			referencePoints[2], )), surfaceRegion=mdb.models[Simu_name].rootAssembly.instances['Niederhalter-1'].surfaces['Niederhalter'])
		mdb.models[Simu_name].RigidBody(name='Biegekern', refPointAtCOM=ON, refPointRegion=Region(referencePoints=(mdb.models[Simu_name].rootAssembly.instances['Biegekern-1'].
			referencePoints[2], )), surfaceRegion=mdb.models[Simu_name].rootAssembly.instances['Biegekern-1'].surfaces['Biegekern'])
		#mdb.models[Simu_name].RigidBody(name='Endposition', refPointAtCOM=ON, refPointRegion=Region(referencePoints=(mdb.models[Simu_name].rootAssembly.instances['Positionierung-1'].
			#referencePoints[3], )), surfaceRegion=mdb.models[Simu_name].rootAssembly.instances['Positionierung-1'].surfaces['Positionierung'])
			
		####################################################### 18.Create step ######################################################################################
		Simu_3D.ExplicitDynamicsStep(massScaling=((SEMI_AUTOMATIC, MODEL, AT_BEGINNING, 0.0, mass_scale, BELOW_MIN, 0, 0, 0.0, 0.0, 0, None), ), name='Biegevorgang', previous='Initial', timePeriod=Zeit)	
			
		######################################################## 19.Create contact property #########################################################################
		Simu_3D.ContactProperty('Friction')
		Simu_3D.interactionProperties['Friction'].TangentialBehavior(dependencies=0, directionality=ISOTROPIC, elasticSlipStiffness=None, formulation=PENALTY, fraction=0.005, 
			maximumElasticSlip=FRACTION, pressureDependency=OFF, shearStressLimit=None, slipRateDependency=OFF, table=((0.12, ), ), temperatureDependency=OFF)
		Simu_3D.interactionProperties['Friction'].NormalBehavior(allowSeparation=ON, constraintEnforcementMethod=DEFAULT, pressureOverclosure=HARD)

		######################################################## 20.Create interaction properties ###################################################################
		mdb.models[Simu_name].SurfaceToSurfaceContactExp(clearanceRegion=None, createStepName='Initial', datumAxis=None, initialClearance=OMIT, interactionProperty='Friction', 
			master=mdb.models[Simu_name].rootAssembly.instances['Biegeblech-1'].surfaces['Blech_Oberseite'], mechanicalConstraint=PENALTY, name='Blech_Biegeform', 
			slave=mdb.models[Simu_name].rootAssembly.instances['Biegeform_2'].surfaces['Biegeform_2'], sliding=FINITE)
		mdb.models[Simu_name].SurfaceToSurfaceContactExp(clearanceRegion=None, createStepName='Initial', datumAxis=None, initialClearance=OMIT, interactionProperty='Friction', 
			master=mdb.models[Simu_name].rootAssembly.instances['Biegeblech-1'].surfaces['Blech_Unterseite'], mechanicalConstraint=PENALTY, name='Blech_Biegekern', 
			slave=mdb.models[Simu_name].rootAssembly.instances['Biegekern-1'].surfaces['Biegekern'], sliding=FINITE)
		mdb.models[Simu_name].SurfaceToSurfaceContactExp(clearanceRegion=None, createStepName='Initial', datumAxis=None, initialClearance=OMIT, interactionProperty='Friction', 
			master=mdb.models[Simu_name].rootAssembly.instances['Biegeblech-1'].surfaces['Blech_Oberseite'], mechanicalConstraint=PENALTY, name='Blech_Niederhalter', 
			slave=mdb.models[Simu_name].rootAssembly.instances['Niederhalter-1'].surfaces['Niederhalter'], sliding=FINITE)		 
		

		######################################################### 21.Create boundary conditions ######################################################################
		mdb.models[Simu_name].EncastreBC(createStepName='Initial', localCsys=None, name='Biegekern', region=mdb.models[Simu_name].rootAssembly.instances['Biegekern-1'].sets['Biegekern'])
		mdb.models[Simu_name].VelocityBC(amplitude=UNSET, createStepName='Initial', distributionType=UNIFORM, fieldName='', localCsys=None, name='Biegeform', region=mdb.models[Simu_name].
			rootAssembly.instances['Biegeform_2'].sets['Biegeform_2_RP'], v1=0.0, v2=0, v3=0.0, vr1=0.0, vr2=0.0, vr3=0.0)
		mdb.models[Simu_name].DisplacementBC(amplitude=UNSET, createStepName='Initial', distributionType=UNIFORM, fieldName='', localCsys=None, name='Sucherloecher', 
			region=mdb.models[Simu_name].rootAssembly.instances['Biegeblech-1'].sets['Sucherloecher'], u1=SET, u2=UNSET, u3=UNSET, ur1=UNSET, ur2=UNSET, ur3=UNSET)	
		mdb.models[Simu_name].VelocityBC(amplitude=UNSET, createStepName='Initial', distributionType=UNIFORM, fieldName='', localCsys=None, name='Niederhalter', region=mdb.models[Simu_name].
			rootAssembly.instances['Niederhalter-1'].sets['Niederhalter_RP'], v1=0.0, v2=0.0, v3=0.0, vr1=0.0, vr2=0.0, vr3=0.0)
		mdb.models[Simu_name].boundaryConditions['Biegeform'].setValuesInStep(stepName='Biegevorgang', v2=-V)
		mdb.models[Simu_name].boundaryConditions['Niederhalter'].setValuesInStep(stepName='Biegevorgang', v2=FREED)
		#mdb.models[Simu_name].DisplacementBC(amplitude=UNSET, createStepName='Initial', distributionType=UNIFORM, fieldName='', localCsys=None, name='Endposition', region=mdb.models[Simu_name].
			#rootAssembly.instances['Positionierung-1'].sets['Inertia_RP_set'], u1=SET, u2=SET, u3=SET, ur1=SET, ur2=SET, ur3=SET)
		
		######################################################### 22.Create loads ####################################################################################
		Simu_3D.ConcentratedForce(cf2=-1678.0, createStepName='Biegevorgang', distributionType=UNIFORM, field='', localCsys=None, name='Niederhalterkraft', region=Simu_3D.
			rootAssembly.instances['Niederhalter-1'].sets['Niederhalter_RP'])

		######################################################### 23.Define Field Output #############################################################################
		Simu_3D.fieldOutputRequests['F-Output-1'].setValues(numIntervals=500, variables=('S', 'SVAVG', 'PE', 'PEVAVG', 'PEEQ', 'PEEQVAVG', 'LE', 'U', 'V', 'A', 'RF', 'CSTRESS', 'EVF', 'COORD'))
		
		######################################################## 24.Create list of nodes #############################################################################
		
			# #create tuple with nodes for path
		# edge_nodes = Blech_part.edges.findAt((-18.0,0.0,30.0),)
		# nodes_path1 = edge_nodes.getNodes()
		# list = []
		
		# for i in range (0,29):					#for i in range (0,imax) with imax = (length.edge/seed)+2 = (13.9/0.5)+2 = 29
			# x = nodes_path1[i]
			# label = x.label
			# list.append(label)
		
		# list =  str(list).strip('[]')				#convert list into string
		# list = tuple(list)						#convert list into tuple
		# #path1 = session.Path(name='Path1',type=NODE_LIST,expression=list)
		
		######################################################### 25.Create and run job ##############################################################################
		mdb.Job(activateLoadBalancing=False, atTime=None, contactPrint=OFF, description='', echoPrint=OFF, explicitPrecision=SINGLE, historyPrint=OFF, 
			memory=90, memoryUnits=PERCENTAGE, model=Simu_name, modelPrint=OFF, multiprocessingMode=DEFAULT, name=Simu_name, nodalOutputPrecision=
			SINGLE, numCpus=1, numDomains=1, parallelizationMethodExplicit=DOMAIN, queue=None, scratch='', type=ANALYSIS, userSubroutine='', waitHours=0, waitMinutes=0)

		#run job
		mdb.jobs[Simu_name].submit()                    
		#mdb.jobs[Simu_name].waitForCompletion()
		
		######################################################### 25.Import Explicit to Standard #####################################################################
			#define name of simu
		SB_name = 'SB_'+Simu_name
			#copy old simu and delete not required objects and features
		SB_Simu = mdb.Model(name=SB_name, objectToCopy=Simu_3D)
		del SB_Simu.rootAssembly.features['Edge to Edge-1']
		del SB_Simu.rootAssembly.features['Biegekern-1']
		del SB_Simu.rootAssembly.features['Niederhalter-1']
		del SB_Simu.rootAssembly.features['Biegeform_2']
		del SB_Simu.parts['Biegeform_2']
		del SB_Simu.parts['Biegekern']
		del SB_Simu.parts['Niederhalter']
		del SB_Simu.interactions['Blech_Biegeform']
		del SB_Simu.interactions['Blech_Biegekern']
		del SB_Simu.interactions['Blech_Niederhalter']
		del SB_Simu.interactionProperties['Friction']
		del SB_Simu.constraints['Biegeform_2']
		del SB_Simu.constraints['Biegekern']
		del SB_Simu.constraints['Niederhalter']
		del SB_Simu.loads['Niederhalterkraft']
		del SB_Simu.boundaryConditions['Biegeform']
		del SB_Simu.boundaryConditions['Biegekern']
		del SB_Simu.boundaryConditions['Niederhalter']
		del SB_Simu.sketches['Biegeform_2']
		del SB_Simu.sketches['Biegeform_2_3D']
		del SB_Simu.sketches['Biegekern']
		del SB_Simu.sketches['Biegekern3D']
		del SB_Simu.sketches['Niederhalter']
		del SB_Simu.sketches['Niederhalter3D']
		del SB_Simu.sketches['Positionierung']
		
		########################################################## 26.New properties ##################################################################################
		
			# change to static step
		del SB_Simu.steps['Biegevorgang']
		SB_Simu.StaticStep(name='Springback', nlgeom=ON, previous='Initial')
			#define initial state
		SB_Simu.InitialState(createStepName='Initial', endIncrement=STEP_END, endStep=LAST_STEP, fileName=Simu_name, instances=(
			SB_Simu.rootAssembly.instances['Biegeblech-1'], ), name='import_stress', updateReferenceConfiguration=ON)
		mdb.Job(atTime=None, contactPrint=OFF, description='', echoPrint=OFF, explicitPrecision=SINGLE, getMemoryFromAnalysis=True, historyPrint=OFF, 
			memory=90, memoryUnits=PERCENTAGE, model=SB_name, modelPrint=OFF, multiprocessingMode=DEFAULT, name=SB_name, nodalOutputPrecision=SINGLE
			, numCpus=1, numGPUs=0, queue=None, scratch='', type=ANALYSIS, userSubroutine='', waitHours=0, waitMinutes=0)
		
		########################################################### 27.Submit and run job Springback Analysis ########################################################
		mdb.jobs[SB_name].submit
		#mdb.jobs['Springback'].waitForCompletion()