"""This module is used to control the output to the Abaqus result
(`.fil`) file
.. note:: Uses direct editing of input and should be called after all
cae options have been set.
.. codeauthor:: Knut Andreas Meyer
"""
from __future__ import print_function
from rollover.utils import naming_mod as names
from rollover.utils import inp_file_edit as inp_edit
[docs]def add(the_model, num_cycles):
""" Add .fil output to input file for each rolling step.
For the first, add node coordinates and displacements. For the
remaining, add only displacements.
If rail reference point is used, add node coord and displacements
If rail substructure is used, the instance set names change...
:param the_model: The model to which the output is to be added
:type the_model: Model object (Abaqus)
:param num_cycles: Number of rollover cycles to simulate
:type num_cycles: int
"""
assy = the_model.rootAssembly
if assy.isOutOfDate:
assy.regenerate()
kwb = the_model.keywordBlock
use_substr = names.rail_substructure in the_model.parts.keys()
rail_rp = names.rail_rp_set if names.rail_rp_set in assy.sets.keys() else None
# Setup output after first rollover
add_to_step(kwb, 'COORD, U', names.get_step_rolling(1), rail_rp, use_substr)
for cycle_nr in range(2,num_cycles+1):
add_to_step(kwb, '', names.get_step_return(cycle_nr), rail_rp, use_substr)
add_to_step(kwb, 'U', names.get_step_rolling(cycle_nr), rail_rp, use_substr)
[docs]def add_to_step(kwb, varstr, step_name, rail_rp=None, use_substr=False):
""" Add output specified to given step.
:param kwb: The keyword block that can be used to edit input file
directly.
:type kwb: KeywordBlock object (Abaqus)
:param varstr: The string specifying which variables to add to
output
:type varstr: str
:param step_name: The name of the step to which output should be
added/modified
:type step_name: str
:param rail_rp: The name of the rail reference point set
:type rail_rp: str
:param use_substr: Is a rail substructure used?
:type use_substr: bool
"""
sep = '_' if use_substr else '.'
wheel_cn_set = names.wheel_inst + sep + names.wheel_contact_nodes
wheel_rp_set = names.wheel_inst + sep + names.wheel_rp_set
sets = [wheel_cn_set, wheel_rp_set]
if rail_rp is not None:
sets.append(rail_rp)
for set in sets:
add_str = get_node_file_output_str(set, varstr)
inp_edit.add_at_end_of_cat(kwb, add_str, category='Step', name=step_name)
[docs]def get_node_file_output_str(nset, varstr, frequency=99999999):
""" Get the string to add to the input file.
:param nset: Name of node set from which output should be output
:type nset: str
:param varstr: comma separated list of variables to output
:type varstr: str
:param frequency: How often to write output (increments)
:type frequency: int
:returns: The string to add
:rtype: str
"""
output_str = ('*NODE FILE, NSET=' + nset +
', FREQUENCY=%0.0f \n'
+ varstr) % (frequency)
return output_str