import sys
sys.path.append("C:\\franc3d\\")

import os
import shutil

import PyF3D

f3d = PyF3D.F3DApp()

# existing files:  cube.inp, cube_RETAINED_ELEMS.txt, cube_02.crk

f3d.Submodeler(
    model_type="ABAQUS",
    orig_file_name='cube.inp',
    submodel_file_name='cube_LOCAL.inp',
    global_file_name='cube_GLOBAL.inp',
    elem_file_name='cube_RETAINED_ELEMS.txt')

f3d.OpenMeshModel(
    model_type="ABAQUS",
    file_name='cube_LOCAL.inp',
    global_name='cube_GLOBAL.inp')

f3d.InsertFileFlaw(
    file_name='cube_02.crk')

f3d.RunAnalysis(
    model_type="ABAQUS",
    file_name='cube_02crk.fdb',
    merge_tol=0.0001,
    connection_type="MERGE",
    command='"abaqus.bat" job=cube_02crk_full -interactive -analysis ',
    global_model='cube_GLOBAL.inp',
    merge_surf_labels=["AUTO_CUT_SURF"],
    global_surf_labels=["GLOBAL_CONNECT_SURF"],
    del_file_list=["cid","com","mdl","stt","res","prt","sim","txt","script","sta","msg"])

f3d.ComputeSif()

f3d.SetUnits(
    model_length="INCH",
    model_stress="KSI",
    temperature="F")

f3d.SetGrowthParams(
    growth_type="QUASI_STATIC",
    fem_units="INCH|KSI|F|UNKNOWN|PASSES",
    load_step_map=["VERSION: 1","SUB_RANGE: 1","0 0 8","LABELS: 1","1 Load Step 1"],
    kink_angle_strategy="MTS",
    quasi_static_n=2,
    quasi_static_loads=["ONE_STEP: 1 NONE 1 1 0"],
    dist_step=0.0225,
    scale_node=0.5,
    cycles_step=1000,
    front_mult=[1])

f3d.GrowCrack(
    dist_step=0.025,
    scale_node=0.5,
    cycles_step=1000,
    front_mult=[1],
    temp_radius=0.015,
    discard=[[1,1]],
    extrapolate=[[2,2]],
    check_fit=False)

f3d.AutoGrowth(
    model_type="ABAQUS",
    cur_step=1,
    file_name='cube_02crk',
    const_median_step=0.025,
    num_steps=3,
    step_type="SCONST",
    const_step_size=0.025,
    temp_radius_type="ABSOLUTE",
    temp_radius=0.015,
    extrapolate=[[3,3]],
    merge_tol=0.0001,
    connection_type="MERGE",
    command='"abaqus.bat" job=cube_02crk_STEP_001_full -interactive -analysis ',
    global_model='cube_GLOBAL.inp',
    merge_surf_labels=["AUTO_CUT_SURF"],
    global_surf_labels=["GLOBAL_CONNECT_SURF"],
    del_file_list=["cid","com","mdl","stt","res","prt","sim","txt","script","sta","msg"])

f3d.ComputeSif()

f3d.WriteSifPath(
    file_name='cube_cube_02crk_sif_hist.sif',
    load_step=1,
    flags=["TAB","A","KI"])