238 lines
6.3 KiB
Python
Executable File
238 lines
6.3 KiB
Python
Executable File
#!/usr/bin/env python3
|
|
|
|
import os
|
|
import sys
|
|
from optparse import OptionParser
|
|
|
|
import numpy as np
|
|
|
|
import damask
|
|
|
|
scriptName = os.path.splitext(os.path.basename(__file__))[0]
|
|
scriptID = ' '.join([scriptName,damask.version])
|
|
sys.path.append(damask.solver.Marc().libraryPath())
|
|
|
|
#-------------------------------------------------------------------------------------------------
|
|
def outMentat(cmd,locals):
|
|
if cmd[0:3] == '(!)':
|
|
exec(cmd[3:])
|
|
elif cmd[0:3] == '(?)':
|
|
cmd = eval(cmd[3:])
|
|
py_mentat.py_send(cmd)
|
|
else:
|
|
py_mentat.py_send(cmd)
|
|
return
|
|
|
|
#-------------------------------------------------------------------------------------------------
|
|
def outFile(cmd,locals,dest):
|
|
if cmd[0:3] == '(!)':
|
|
exec(cmd[3:])
|
|
elif cmd[0:3] == '(?)':
|
|
cmd = eval(cmd[3:])
|
|
dest.write(cmd+'\n')
|
|
else:
|
|
dest.write(cmd+'\n')
|
|
return
|
|
|
|
#-------------------------------------------------------------------------------------------------
|
|
def output(cmds,locals,dest):
|
|
for cmd in cmds:
|
|
if isinstance(cmd,list):
|
|
output(cmd,locals,dest)
|
|
else:
|
|
if dest == 'Mentat':
|
|
outMentat(str(cmd),locals)
|
|
else:
|
|
outFile(str(cmd),locals,dest)
|
|
return
|
|
|
|
|
|
|
|
#-------------------------------------------------------------------------------------------------
|
|
def init():
|
|
return [
|
|
"|"+' '.join([scriptID] + sys.argv[1:]),
|
|
"*draw_manual", # prevent redrawing in Mentat, should be much faster
|
|
"*new_model yes",
|
|
"*reset",
|
|
"*select_clear",
|
|
"*set_element_class hex8",
|
|
"*set_nodes off",
|
|
"*elements_solid",
|
|
"*show_view 4",
|
|
"*reset_view",
|
|
"*view_perspective",
|
|
"*redraw",
|
|
]
|
|
|
|
|
|
#-------------------------------------------------------------------------------------------------
|
|
def mesh(r,d):
|
|
return [
|
|
"*add_nodes",
|
|
"%f %f %f"%(0.0,0.0,0.0),
|
|
"%f %f %f"%(0.0,0.0,d[2]),
|
|
"%f %f %f"%(0.0,d[1],d[2]),
|
|
"%f %f %f"%(0.0,d[1],0.0),
|
|
"%f %f %f"%(-d[0],0.0,0.0),
|
|
"%f %f %f"%(-d[0],0.0,d[2]),
|
|
"%f %f %f"%(-d[0],d[1],d[2]),
|
|
"%f %f %f"%(-d[0],d[1],0.0),
|
|
"*add_elements",
|
|
"1",
|
|
"2",
|
|
"3",
|
|
"4",
|
|
"5",
|
|
"6",
|
|
"7",
|
|
"8",
|
|
"*sub_divisions",
|
|
"%i %i %i"%(r[2],r[1],r[0]),
|
|
"*subdivide_elements",
|
|
"all_existing",
|
|
"*set_sweep_tolerance",
|
|
"%f"%(float(min(d))/max(r)/2.0),
|
|
"*sweep_all",
|
|
"*renumber_all",
|
|
"*set_move_scale_factor x -1",
|
|
"*move_elements",
|
|
"all_existing",
|
|
"*flip_elements",
|
|
"all_existing",
|
|
"*fill_view",
|
|
]
|
|
|
|
|
|
#-------------------------------------------------------------------------------------------------
|
|
def material():
|
|
cmds = [\
|
|
"*new_mater standard",
|
|
"*mater_option general:state:solid",
|
|
"*mater_option structural:type:hypo_elast",
|
|
"*mater_name",
|
|
"hypela2",
|
|
"*add_mater_elements",
|
|
"all_existing",
|
|
"*geometry_type mech_three_solid",
|
|
# "*geometry_option red_integ_capacity:on", reduced integration with one IP gave trouble being always OUTDATED...
|
|
"*add_geometry_elements",
|
|
"all_existing",
|
|
]
|
|
|
|
return cmds
|
|
|
|
|
|
#-------------------------------------------------------------------------------------------------
|
|
def geometry():
|
|
cmds = [\
|
|
"*geometry_type mech_three_solid",
|
|
# "*geometry_option red_integ_capacity:on",
|
|
"*add_geometry_elements",
|
|
"all_existing",
|
|
# we are NOT using reduced integration (type 117) but opt for /elementhomogeneous/ in the respective phase description (material.config)
|
|
"*element_type 7",
|
|
"all_existing",
|
|
]
|
|
|
|
return cmds
|
|
|
|
|
|
#-------------------------------------------------------------------------------------------------
|
|
def initial_conditions(homogenization,microstructures):
|
|
elements = []
|
|
element = 0
|
|
for id in microstructures:
|
|
element += 1
|
|
if len(elements) < id:
|
|
for i in range(id-len(elements)):
|
|
elements.append([])
|
|
elements[id-1].append(element)
|
|
|
|
cmds = [\
|
|
"*new_icond",
|
|
"*icond_name _temperature",
|
|
"*icond_type state_variable",
|
|
"*icond_param_value state_var_id 1",
|
|
"*icond_dof_value var 300",
|
|
"*add_icond_elements",
|
|
"all_existing",
|
|
"*new_icond",
|
|
"*icond_name _homogenization",
|
|
"*icond_type state_variable",
|
|
"*icond_param_value state_var_id 2",
|
|
"*icond_dof_value var %i"%homogenization,
|
|
"*add_icond_elements",
|
|
"all_existing",
|
|
]
|
|
|
|
for grain,elementList in enumerate(elements):
|
|
cmds.append([\
|
|
"*new_icond",
|
|
"*icond_name microstructure_%i"%(grain+1),
|
|
"*icond_type state_variable",
|
|
"*icond_param_value state_var_id 3",
|
|
"*icond_dof_value var %i"%(grain+1),
|
|
"*add_icond_elements",
|
|
elementList,
|
|
"#",
|
|
])
|
|
return cmds
|
|
|
|
|
|
#--------------------------------------------------------------------------------------------------
|
|
# MAIN
|
|
#--------------------------------------------------------------------------------------------------
|
|
|
|
parser = OptionParser(option_class=damask.extendableOption, usage='%prog options [file[s]]', description = """
|
|
Generate MSC.Marc FE hexahedral mesh from geom file.
|
|
|
|
""", version = scriptID)
|
|
|
|
parser.add_option('-p', '--port',
|
|
dest = 'port',
|
|
type = 'int', metavar = 'int',
|
|
help = 'Mentat connection port [%default]')
|
|
|
|
parser.set_defaults(port = None,
|
|
)
|
|
|
|
(options, filenames) = parser.parse_args()
|
|
|
|
if options.port is not None:
|
|
try:
|
|
import py_mentat
|
|
except ImportError:
|
|
parser.error('no valid Mentat release found.')
|
|
|
|
# --- loop over input files ------------------------------------------------------------------------
|
|
|
|
if filenames == []: filenames = [None]
|
|
|
|
for name in filenames:
|
|
damask.util.report(scriptName,name)
|
|
|
|
geom = damask.Geom.from_file(StringIO(''.join(sys.stdin.read())) if name is None else name)
|
|
microstructure = geom.get_microstructure().flatten(order='F')
|
|
|
|
cmds = [\
|
|
init(),
|
|
mesh(geom.grid,geom.size),
|
|
material(),
|
|
geometry(),
|
|
initial_conditions(geom.homogenization,microstructure),
|
|
'*identify_sets',
|
|
'*show_model',
|
|
'*redraw',
|
|
'*draw_automatic',
|
|
]
|
|
|
|
outputLocals = {}
|
|
if options.port:
|
|
py_mentat.py_connect('',options.port)
|
|
output(cmds,outputLocals,'Mentat')
|
|
py_mentat.py_disconnect()
|
|
else:
|
|
with sys.stdout if name is None else open(os.path.splitext(name)[0]+'.proc','w') as f:
|
|
output(cmds,outputLocals,f)
|