#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'))