#####################################################################
# Study of sensitivity of shell elements with respect to skewness,  # 
# element type, and mesh density                                    # 
# Parameters used in the parametric study:                          #
#    level: mesh density level (level 1.=4X4, level 2.=8X8)         #
#    delta: the skew angle (in degrees)                             #
#    elemType: the shell element type                               #
#####################################################################

# create the study
skewXpl = parStudy(par=('level', 'elemType', 'delta'),
		   directory=ON, verbose=ON)

# define the parameters
skewXpl.define(DISCRETE, par='level', domain=(1.,2.))
skewXpl.define(DISCRETE, par='delta', domain=(90., 80., 60., 40., 30.))
skewXpl.define(DISCRETE, par='elemType', domain=('s4r', 's4rs'))

# sample each parameter
skewXpl.sample(INTERVAL, par='level', interval=1)
skewXpl.sample(NUMBER, par='delta', number=3)
skewXpl.sample(INTERVAL, par='elemType', interval=1)

# combine the samples into design points
skewXpl.combine(MESH)

# generate analysis job data
skewXpl.generate(template='skew_discr') 

# execute all runs sequentially
skewXpl.execute(ALL)

# study output at end of step 1
skewXpl.output(step=1, inc=LAST)

# gather displacement results at node 405
skewXpl.gather(results='N405_U', variable='U' ,node=405)

# screen report of U3 displacement results at node 405 
skewXpl.report(PRINT, results=('N405_U.3'))

# for QA purposes
skewXpl.report(FILE, results=('N405_U.3'), par=(),
	       file='skew_discr.psr')

# read same set of results from ODB file

# study output at end of step 1
skewXpl.output(step=1, inc=LAST, file=ODB)

# gather displacement results at node 405
skewXpl.gather(results='o N405_U', variable='U' ,node=405)

# screen report of U3 displacement results at node 405 
skewXpl.report(PRINT, results=('o N405_U.3'))

# for QA purposes
skewXpl.report(FILE, results=('o N405_U.3'), par=(),
	       file='skew_discr.psr')