#to study the effect of prestrain on fld curve predicted by m-k: biaxial prestrain
test=ParStudy(par='beta',name='prestrain')
test.define(DISCRETE,par='beta',domain=(-0.7,-0.6,-0.5,-0.4,-0.3,-0.2,-0.1,0.1,0.2,0.3,0.4,0.5,0.6,0.7))
test.sample(NUMBER,par='beta',number=14)
test.combine(MESH)
test.generate(template='damage_prestrain_biaxial')
test.execute(ALL)
test.output(file=ODB)
test.gather(results='lep1',step=2,request=HISTORY,variable='LEP1',element=1);
test.gather(results='lep2',step=2,request=HISTORY,variable='LEP2',element=1);
test.report(XYPLOT,file='output_biaxial.dat',par='beta',results=('lep1','lep2'))