"""
Parameter study: Young's modulus of the holder

This script performs a parameter study varying the Young's modulus of the
holder component while all other parameters remain constant.

The script calls the Abaqus model generation script HSK_63_axi.py and
generates a list of job names for subsequent evaluation.

Author: F. Birk
Context: Dissertation – numerical investigation of HSK shaft–hub connections
"""

import numpy as np

joblist = []

for g_Emod_H in range(10000, 400000, 10000):

    g_Emod_A = 210000.0
    g_q = 0
    g_F_sp = 18000.0
    g_Umin = 25000.0
    g_reib = 0.15

    string = str(g_Emod_H).replace('.', '_')
    g_Jobname = 'Berechnung-' + string + '-mit-er'
    joblist.append(g_Jobname)

    execfile('HSK_63_axi.py', __main__.__dict__)
    Mdb()

    g_q = 1
    g_Jobname = 'Berechnung-' + string + '-mit-fr'
    joblist.append(g_Jobname)

    execfile('HSK_63_axi.py', __main__.__dict__)

paramsFile = open('joblist_E.txt', 'w')
for job in joblist:
    paramsFile.write(str(job) + '\n')
paramsFile.close()