############################################################################
# Mesh convergence study for plate pentration verification                 #
# problem in ABAQUS/Explicit                                               # 
#                                                                          #
# Study parameters:                                                        #
#       elno: Number of elements in the horizontal direction of the plate  #
#       rows: Number of elements through the thickness of the plate        #
############################################################################

# Create the study
mesh = ParStudy(par=('elno','rows'))

# Define the parameters
mesh.define(DISCRETE, par='elno', domain=(50,70,90))
mesh.define(DISCRETE, par='rows', domain=(5,7,9))

# Sample the parameters
mesh.sample(INTERVAL, par='elno', interval=1)
mesh.sample(INTERVAL, par='rows', interval=1)

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

# Generate analysis job data
mesh.generate(template='pp_mesh_study')

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

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

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

# Report the gathered results to an output file
mesh.report(FILE, results=('vel.2'), par=('elno','rows'),
	    file='pp_mesh_study.psr')

# read same set of results from ODB file

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

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

# Report the gathered results to an output file
mesh.report(FILE, results=('o vel'), par=('elno','rows'),
	    file='pp_mesh_study.psr')