template=GetTemplate(TemplateName="EngData")
system=template.CreateSystem(Position="Default")
container=system.GetContainer(ComponentName="Engineering Data")
mats=container.CreateMaterial(Name="LSR2030")
strain=[]
unia=[]
biax=[]
shea=[]
SetUserPathRoot (DirectoryPath="U:/praktikanten/Cornelia/FEM/ANSYS/Experimentelle Daten")
fileName=AbsUserPathName("testdataset1.txt")
dataFile=open(fileName,"r")
for line in dataFile:
	(tstrain,tunia,tbiax,tshea)=line.split(",")
	tstrain=Quantity(float(tstrain),"m m^-1")
	tunia=Quantity(float(tunia),"MPa")
	tbiax=Quantity(float(tbiax),"MPa")
	tshea=Quantity(float(tshea),"MPa")
	strain.append(tstrain)
	unia.append(tunia)
	biax.append(tbiax)
	shea.append(tshea)
C10=[1]
C01=[1]
D1=[0]
miss=mats.CreateProperty(Name="Mooney-Rivlin")
miss.SetData(
 SheetName="Mooney-Rivlin",
 Variables=["Material Constant C10", "Material Constant C01", "Incompressibility Parameter D1"], 
 Values=[C10,C01,D1])

up=mats.CreateProperty(Name="Uniaxial Test Data")
up.SetData(
Variables=["strain","stress"],
Values=[strain,unia])
bp=mats.CreateProperty(Name="Biaxial Test Data")
bp.SetData(
Variables=["strain","stress"],
Values=[strain, biax])
sp=mats.CreateProperty(Name="Shear Test Data")
sp.SetData(
Variables=["strain","stress"],
Values=[strain, shea])

mrd=miss.GetPropertyData(Name="Mooney-Rivlin")
curveFit=mrd.CreateCurveFitting(
Type="Mooney-Rivlin",
Definition="3 Parameter")
uniax=LSR2030.GetProperty(Name="Uniaxial Test Data")
curveFit.AddTestData(TestData=uniax)
ebiax=LSR2030.GetProperty(Name="Biaxial Test Data")
curveFit.AddTestData(TestData=ebiax)
shear=LSR2030.GetProperty(Name="Shear Test Data")
curveFit.AddTestData(TestData=shear)
curveFit.Solve()
curveFit.CopyCoefficients(Destination=mrd)