###################################################################
# Study of sensitivity of shell elements with respect to skewness #
# Parameters used in the parametric study:                        #
#    delta: the skew angle (in degrees)                           #
#    elemType: the shell element type                             #
###################################################################

# create the study
skewStudy = ParStudy(par=('elemType', 'delta'),
		     directory=ON, verbose=ON)

# define the parameters
skewStudy.define(DISCRETE, par='delta', domain=(90., 80., 60., 40., 30.))
skewStudy.define(DISCRETE, par='elemType', domain=('s8r', 's4r', 's4'))

# sample the parameters
skewStudy.sample(INTERVAL, par=('delta','elemType'), interval=1)

# combine the samples to give the designs
skewStudy.combine(MESH)

# generate analysis data
skewStudy.generate(template='skewshell_parametric') 

# execute all analysis jobs sequentially
skewStudy.execute(ALL)

# study output at end of step 1
skewStudy.output(step=1)

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

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

# read same set of results from ODB file

# parametric study output at end of step 1
skewStudy.output(step=1, file=ODB)

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

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

skewStudy.report(FILE, results=('o N405_U.3'), par=(), 
		 file='skewshell_parametric.psr')