###################################################################
# Mesh convergence study for plate pentration verification        #
# problem in ABAQUS/Explicit                                      #
#                                                                 #
# Study parameters:                                               #
#       fstrain: Shear failure strain                             #
###################################################################

# Create the study
matl = ParStudy(par='fstrain')

# Define the parameter
matl.define(DISCRETE, par='fstrain', domain=('0.50','0.17'))

# CONSIDER SHEAR FAILURE
# Sample the parameter
matl.sample(INTERVAL, par='fstrain', interval=1)

# Combine the samples into design points
matl.combine(MESH, name='matl_2_study')

# Generate analysis job data
matl.generate(template='pp_mat_2_study')

# Execute all analysis jobs sequentially
matl.execute(ALL)

# Study output at end of step 1
matl.output(step=1)

# Gather the results for the velocity of the projectile
matl.gather(results='vel', variable='V', node=99999)

# Report the gathered results to an output file
matl.report(FILE, results=('vel.2'), par=('fstrain'),
	    file='pp_mat_2_study.psr')

# read same set of results from ODB file

# Study output at end of step 1
matl.output(step=1, file=ODB)

# Gather the results for the velocity of the projectile
matl.gather(results='o vel', variable='V2', node=99999, request=HISTORY)

# Report the gathered results to an output file
matl.report(FILE, results=('o vel'), par=('fstrain'),
	    file='pp_mat_2_study.psr')