##to study the effect of prestrain on fld curve predicted by msfld
test=ParStudy(par='beta_step2',name='failure')
test.define(DISCRETE,par='beta_step2',domain=(-0.9,-0.8,-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_step2',number=16)
test.combine(MESH)
test.generate(template='damage_msfld_m0p6_higher')
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.gather(results='peeq',step=2,request=HISTORY,variable='PEEQ',element=1);
test.report(XYPLOT,file='output_m0p6_higher.dat',par='beta_step2',results=('peeq','lep1','lep2'))