def __init__(self):
OP2_F06_Common.__init__(self)
self.card_count = {}
self.additional_matrices = {}
self.matrices = {}
self.subcase_key = defaultdict(list)
self.end_options = {}
self._results = ResultSet(self.get_all_results())
def __init__(self):
OP2_F06_Common.__init__(self)
self.card_count = {}
self.additional_matrices = {}
self.matrices = {}
self._results = ResultSet(self.get_all_results())
def __init__(self):
OP2_F06_Common.__init__(self)
self.card_count = {}
self.additional_matrices = {}
self.matrices = {}
self.subcase_key = defaultdict(list)
self._results = ResultSet(self.get_all_results())
class F06Writer(OP2_F06_Common):
def __init__(self):
OP2_F06_Common.__init__(self)
self.card_count = {}
self.additional_matrices = {}
self.matrices = {}
self._results = ResultSet(self.get_all_results())
def get_all_results(self):
all_results = ['stress', 'strain', 'element_forces', 'constraint_forces'] + self.get_table_types()
return all_results
def _clear_results(self):
self._results.clear()
def add_results(self, results):
if isinstance(results, string_types):
results = [results]
all_results = self.get_all_results()
for result in results:
result = str(result)
if result not in all_results:
raise RuntimeError('%r is not a valid result to remove; all_results=%s' % (result, all_results))
if 'stress' == result:
stress_results = []
for result in all_results:
if 'stress' in result.lower():
stress_results.append(result)
#stress_results = [result if 'stress' in result.lower() for result in all_results]
self._results.update(stress_results)
elif 'strain' == result:
strain_results = []
for result in all_results:
if 'strain' in result.lower():
strain_results.append(result)
#strain_results = [result if 'strain' in result.lower() for result in all_results]
self._results.update(strain_results)
elif 'stress' in result.lower():
self._results.add('stress')
elif 'strain' in result.lower():
self._results.add('strain')
elif 'spc_forces' == result or 'mpc_forces' == result or 'constraint_forces' == result:
self._results.add('constraint_forces')
elif 'force' in result.lower(): # could use more validation...
self._results.add('element_forces')
# thermalLoad_VU_3D, thermalLoad_1D, thermalLoad_CONV, thermalLoad_2D_3D
self._results.add(result)
def set_results(self, results):
if isinstance(results, string_types):
results = [results]
self._clear_results()
self.add_results(results)
def remove_results(self, results):
self._results.remove(results)
def make_f06_header(self):
"""If this class is inherited, the F06 Header may be overwritten"""
return make_f06_header()
def make_stamp(self, Title, today):
"""If this class is inherited, the PAGE stamp may be overwritten"""
return make_stamp(Title, today)
def make_grid_point_singularity_table(self, failed):
msg = ''
if failed:
msg += '0 G R I D P O I N T S I N G U L A R I T Y T A B L E\n'
msg += '0 POINT TYPE FAILED STIFFNESS OLD USET NEW USET\n'
msg += ' ID DIRECTION RATIO EXCLUSIVE UNION EXCLUSIVE UNION\n'
for (nid, dof) in failed:
msg += ' %8s G %s 0.00E+00 B F SB SB *\n' % (nid, dof)
else:
msg += 'No constraints have been applied...\n'
page_stamp = self.make_stamp(self.Title, self.date)
msg += page_stamp % self.page_num
self.page_num += 1
return msg
def write_summary(self, f06, card_count=None):
summary_header = ' M O D E L S U M M A R Y\n\n'
summary = ''
self.cards_to_read = set([
# rigid elements
'RBAR', 'RBAR1', 'RBE1', 'RBE2', 'RBE3',
# spc/mpc constraints
'SPC', 'SPCADD', 'SPC1', 'SPCD', 'SPCAX',
'MPC', 'MPCADD',
'SUPORT', 'SUPORT1',
# aero cards
'CAERO1', 'CAERO2', 'CAERO3', 'CAERO4', 'CAERO5',
# temperature cards
'CHBDYE', 'CHBDYG', 'CHBDYP',
'CONV',
])
blocks = [
['POINTS', ['GRID', 'GRDSET', ]],
['ENTRIES', ['SPOINT']],
['ELEMENTS',
[
# these are sorted
# elements
'CONM1', 'CONM2', 'CMASS1', 'CMASS2', 'CMASS3', 'CMASS4',
# springs
'CELAS1', 'CELAS2', 'CELAS3', 'CELAS4', 'CELAS5',
# bushings
'CBUSH', 'CBUSH1D', 'CBUSH2D',
# dampers
'CDAMP1', 'CDAMP2', 'CDAMP3', 'CDAMP4', 'CDAMP5',
# bar flags
'BAROR', 'CBARAO',
# bars
'CBAR', 'CROD', 'CTUBE', 'CBEAM', 'CBEAM3', 'CONROD', 'CBEND',
# shells
'CTRIA3', 'CTRIA6', 'CTRIAR', 'CTRIAX', 'CTRIAX6',
'CQUAD4', 'CQUAD8', 'CQUADR', 'CQUADX', 'CQUAD',
# solids
'CTETRA', 'CPENTA', 'CHEXA',
# other
'CSHEAR', 'CVISC', 'CRAC2D', 'CRAC3D',
'CGAP', 'CFAST', 'RBE2', 'RBE3',
# thermal
'CHBDYP', 'CHBDYG', 'CONV',
]],
]
#print("self.card_count", self.card_count)
if card_count is None:
card_count = self.card_count
for block in blocks:
block_name, keys = block
key_count = 0
for key in sorted(keys):
try:
value = card_count[key]
summary += ' NUMBER OF %-8s %-8s = %8s\n' % (key, block_name, value)
key_count += 1
except KeyError:
pass
if key_count:
summary += ' \n'
if summary:
f06.write(summary_header)
f06.write(summary)
page_stamp = self.make_stamp(self.Title, self.date)
f06.write(page_stamp % self.page_num)
self.page_num += 1
def write_f06(self, f06_outname, is_mag_phase=False,
delete_objects=True, end_flag=False):
"""
Writes an F06 file based on the data we have stored in the object
:param self: the F06 object
:param f06_outname: the name of the F06 file to write
:param is_mag_phase: should complex data be written using Magnitude/Phase
instead of Real/Imaginary (default=False; Real/Imag)
Real objects don't use this parameter.
:param delete_objects: should objects be deleted after they're written
to reduce memory (default=True)
:param end_flag: should a dummy Nastran "END" table be made
(default=False)
"""
print("F06:")
if isinstance(f06_outname, str):
if PY2:
f06 = open(f06_outname, 'wb')
else:
f06 = open(f06_outname, 'w')
self.write_summary(f06)
else:
assert isinstance(f06_outname, file), 'type(f06_outname)= %s' % f06_outname
f06 = f06_outname
f06_outname = f06.name
print('f06_outname =', f06_outname)
page_stamp = self.make_stamp(self.Title, self.date)
if self.grid_point_weight.reference_point is not None:
print("grid_point_weight")
self.page_num = self.grid_point_weight.write_f06(f06, page_stamp, self.page_num)
assert isinstance(self.page_num, int), self.grid_point_weight.__class__.__name__
if self.oload_resultant is not None:
self.page_num = self.oload_resultant.write_f06(f06, page_stamp, self.page_num)
assert isinstance(self.page_num, int), self.oload_resultant.__class__.__name__
#is_mag_phase = False
header = [' DEFAULT \n',
'\n', '']
# eigenvalues are written first
f06.write(page_stamp % self.page_num)
self.page_num += 1
for ikey, result in sorted(iteritems(self.eigenvalues)):
print('%-18s case=%r' % (result.__class__.__name__, ikey))
self.page_num = result.write_f06(f06, header, page_stamp,
page_num=self.page_num)
assert isinstance(self.page_num, int), 'pageNum=%r' % str(self.page_num)
if delete_objects:
del result
self.page_num += 1
# then eigenvectors
# has a special header
iSubcases = sorted(self.iSubcaseNameMap.keys())
res_keys = []
for isubcase in iSubcases:
for subtitle in self.subtitles[isubcase]:
res_keys.append((isubcase, subtitle))
#print('res_keys=%s' % res_keys)
for res_key in res_keys:
isubcase = res_key[0]
subtitle = res_key[1]
label = self.labels[(isubcase, subtitle)]
is_compressed = len(self.subtitles[isubcase]) == 1
if is_compressed:
res_key = isubcase
if res_key not in self.eigenvectors:
continue
result = self.eigenvectors[res_key]
header[0] = ' %s\n' % subtitle
header[1] = '0 SUBCASE %i\n' % (isubcase)
#header[2] = complex/nonlinear
star = ' '
if hasattr(self, 'data'):
star = '*'
res_length = 18
res_format = '%%-%is SUBCASE=%%i' % res_length
res_format_vectorized = '%%-%is SUBCASE=%%i SUBTITLE=%%s' % res_length
if is_compressed:
print(star + res_format % (result.__class__.__name__, isubcase))
else:
print(star + res_format_vectorized % (result.__class__.__name__, isubcase, subtitle))
self.page_num = result.write_f06(header, page_stamp,
self.page_num, f=f06, is_mag_phase=is_mag_phase)
assert isinstance(self.page_num, int), 'pageNum=%r' % str(self.page_num)
if delete_objects:
del result
self.page_num += 1
# finally, we writte all the other tables
# nastran puts the tables in order of the Case Control deck,
# but we're lazy so we just hardcode the order
# subcase name, subcase ID, transient word & value
header_old = [' DEFAULT \n',
'\n', ' \n']
header = copy.deepcopy(header_old)
res_types = [
self.accelerations,
self.displacements, self.displacementsPSD, self.displacementsATO, self.displacementsRMS,
self.scaledDisplacements, # ???
self.force_vectors,
self.load_vectors,
self.temperatures,
self.velocities, #self.eigenvectors,
self.mpc_forces,
self.spc_forces,
self.thermal_load_vectors,
self.strain_energy,
#------------------------------------------
# OEF - forces
# alphabetical order...
# bars
self.cbar_force,
self.cbar100_force,
# beam
self.cbend_force,
self.cbeam_force,
# alphabetical
self.celas1_force,
self.celas2_force,
self.celas3_force,
self.celas4_force,
self.conrod_force,
self.cquad4_force,
self.cquad8_force,
self.cquadr_force,
self.conrod_force,
self.crod_force,
self.cshear_force,
self.ctria3_force,
self.ctria6_force,
self.ctriar_force,
self.ctube_force,
# springs
self.celas1_force,
self.celas2_force,
self.celas3_force,
self.celas4_force,
# dampers
self.cdamp1_force,
self.cdamp2_force,
self.cdamp3_force,
self.cdamp4_force,
# other
self.cbush_force,
self.cgap_force,
self.cvisc_force,
self.chexa_pressure_force,
self.cpenta_pressure_force,
self.ctetra_pressure_force,
#------------------------------------------
# OES - strain
# 1. cbar
# 2. cbeam
# 3. crod/ctube/conrod
# springs
self.celas1_strain,
self.celas2_strain,
self.celas3_strain,
self.celas4_strain,
self.nonlinear_celas1_stress,
self.nonlinear_celas3_stress,
# bars/beams
self.cbar_strain,
self.cbeam_strain,
# plates
self.cquad4_composite_strain,
self.cquad8_composite_strain,
self.cquadr_composite_strain,
self.ctria3_composite_strain,
self.ctria6_composite_strain,
self.ctriar_composite_strain,
self.nonlinear_ctria3_strain,
self.nonlinear_cquad4_strain,
self.ctriax_strain,
# rods
self.nonlinear_crod_strain,
self.nonlinear_ctube_strain,
self.nonlinear_conrod_strain,
self.chexa_strain,
self.conrod_strain,
self.cpenta_strain,
self.cquad4_strain,
self.cquad8_strain,
self.cquadr_strain,
self.crod_strain,
self.cshear_strain,
self.ctetra_strain,
self.ctria3_strain,
self.ctria6_strain,
self.ctriar_strain,
self.ctube_strain,
# bush
self.cbush_strain,
self.nonlinear_cbush_stress,
self.cbush1d_stress_strain,
#------------------------------------------
# cbars/cbeams
self.cbar_stress,
self.nonlinear_cbeam_stress,
self.cbeam_stress,
# bush
self.cbush_stress,
# rods
self.nonlinear_crod_stress,
self.nonlinear_ctube_stress,
self.nonlinear_conrod_stress,
# shear
# OES - stress
self.celas1_stress,
self.celas2_stress,
self.celas3_stress,
self.celas4_stress,
self.chexa_stress,
self.conrod_stress,
self.cpenta_stress,
self.cquad4_stress,
self.cquad8_stress,
self.cquadr_stress,
self.crod_stress,
self.cshear_stress,
self.ctetra_stress,
self.ctria3_stress,
self.ctria6_stress,
self.ctriar_stress,
self.ctube_stress,
self.cquad4_composite_stress,
self.cquad8_composite_stress,
self.cquadr_composite_stress,
self.ctria3_composite_stress,
self.ctria6_composite_stress,
self.ctriar_composite_stress,
self.nonlinear_ctria3_stress,
self.nonlinear_cquad4_stress,
self.ctriax_stress,
self.hyperelastic_cquad4_strain,
#------------------------------------------
self.grid_point_stresses, self.grid_point_volume_stresses, self.grid_point_forces,
]
for res_key in res_keys:
#title = self.Title
isubcase = res_key[0]
subtitle = res_key[1]
label = self.labels[(isubcase, subtitle)]
is_compressed = len(self.subtitles[isubcase]) == 1
#print(self.subtitles[isubcase])
if is_compressed:
res_key = isubcase
#header[0] = ' %-127s\n' % subtitle
#header[1] = '0 %-72s SUBCASE %-15i\n' % (label, isubcase)
#header[1] = '0 %-72s SUBCASE %-15i\n' % ('',isubcase)
res_length = self._get_result_length(res_types, res_key)
if res_length == 0:
# skipped subcase; no saved results
continue
res_format = '%%-%is SUBCASE=%%i%%s' % res_length
res_format_vectorized = '%%-%is SUBCASE=%%i SUBTITLE=%%s %%s' % res_length
for res_type in res_types:
#print("res_type ", res_type)
header = ['', '']
#header[0] = ' %s\n' % subtitle
header[0] = ' %-126s\n' % subtitle
header[1] = '0 %-32s SUBCASE %-15i\n \n' % (label, isubcase)
#print("res_type = %s" % res_type)
if res_key in res_type:
#header = copy.deepcopy(headerOld) # fixes bug in case
result = res_type[res_key]
if result.nonlinear_factor is not None:
header.append('')
try:
element_name = ''
if hasattr(result, 'element_name'):
element_name = ' - ' + result.element_name
star = '*'
if hasattr(result, 'data'):
star = ' '
if is_compressed:
print(star + res_format % (result.__class__.__name__, isubcase, element_name))
else:
print(star + res_format_vectorized % (result.__class__.__name__, isubcase, subtitle, element_name))
self.page_num = result.write_f06(header, page_stamp, page_num=self.page_num, f=f06, is_mag_phase=False)
assert isinstance(self.page_num, int), 'pageNum=%r' % str(self.page_num)
except:
#print("result name = %r" % result.name())
raise
if delete_objects:
del result
self.page_num += 1
f06.write(make_end(end_flag))
f06.close()
class F06Writer(OP2_F06_Common):
def __init__(self):
OP2_F06_Common.__init__(self)
self.card_count = {}
self.additional_matrices = {}
self.matrices = {}
self.subcase_key = defaultdict(list)
self.end_options = {}
self._results = ResultSet(self.get_all_results())
def get_all_results(self):
all_results = [
'stress', 'strain', 'element_forces', 'constraint_forces'
] + self.get_table_types()
return all_results
def clear_results(self):
self._results.clear()
def add_results(self, results):
if isinstance(results, string_types):
results = [results]
all_results = self.get_all_results()
for result in results:
result = str(result)
if result not in all_results:
raise RuntimeError(
'%r is not a valid result to remove; all_results=%s' %
(result, all_results))
if 'stress' == result:
stress_results = []
for result in all_results:
if 'stress' in result.lower():
stress_results.append(result)
#stress_results = [result if 'stress' in result.lower() for result in all_results]
self._results.update(stress_results)
elif 'strain' == result:
strain_results = []
for result in all_results:
if 'strain' in result.lower():
strain_results.append(result)
#strain_results = [result if 'strain' in result.lower() for result in all_results]
self._results.update(strain_results)
elif 'stress' in result.lower():
self._results.add('stress')
elif 'strain' in result.lower():
self._results.add('strain')
elif result in ('spc_forces', 'mpc_forces', 'constraint_forces'):
self._results.add('constraint_forces')
elif 'force' in result.lower(): # could use more validation...
self._results.add('element_forces')
# thermalLoad_VU_3D, thermalLoad_1D, thermalLoad_CONV, thermalLoad_2D_3D
self._results.add(result)
def set_results(self, results):
if isinstance(results, string_types):
results = [results]
self.clear_results()
self.add_results(results)
def remove_results(self, results):
self._results.remove(results)
def make_f06_header(self):
"""If this class is inherited, the F06 Header may be overwritten"""
return make_f06_header()
def make_stamp(self, Title, today):
"""If this class is inherited, the PAGE stamp may be overwritten"""
return make_stamp(Title, today)
def make_grid_point_singularity_table(self, failed):
msg = ''
if failed:
msg += '0 G R I D P O I N T S I N G U L A R I T Y T A B L E\n'
msg += '0 POINT TYPE FAILED STIFFNESS OLD USET NEW USET\n'
msg += ' ID DIRECTION RATIO EXCLUSIVE UNION EXCLUSIVE UNION\n'
for (nid, dof) in failed:
msg += ' %8s G %s 0.00E+00 B F SB SB *\n' % (
nid, dof)
else:
msg += 'No constraints have been applied...\n'
page_stamp = self.make_stamp(self.title, self.date)
msg += page_stamp % self.page_num
self.page_num += 1
return msg
def _write_summary(self, f06, card_count=None):
summary_header = ' M O D E L S U M M A R Y\n\n'
summary = ''
self.cards_to_read = set([
# rigid elements
'RBAR',
'RBAR1',
'RBE1',
'RBE2',
'RBE3',
# spc/mpc constraints
'SPC',
'SPCADD',
'SPC1',
'SPCD',
'SPCAX',
'MPC',
'MPCADD',
'SUPORT',
'SUPORT1',
# aero cards
'CAERO1',
'CAERO2',
'CAERO3',
'CAERO4',
'CAERO5',
# temperature cards
'CHBDYE',
'CHBDYG',
'CHBDYP',
'CONV',
])
blocks = [
['POINTS', [
'GRID',
'GRDSET',
]],
['ENTRIES', ['SPOINT']],
[
'ELEMENTS',
[
# these are sorted
# elements
'CONM1',
'CONM2',
'CMASS1',
'CMASS2',
'CMASS3',
'CMASS4',
# springs
'CELAS1',
'CELAS2',
'CELAS3',
'CELAS4',
'CELAS5',
# bushings
'CBUSH',
'CBUSH1D',
'CBUSH2D',
# dampers
'CDAMP1',
'CDAMP2',
'CDAMP3',
'CDAMP4',
'CDAMP5',
# bar flags
'BAROR',
'CBARAO',
# bars
'CBAR',
'CROD',
'CTUBE',
'CBEAM',
'CBEAM3',
'CONROD',
'CBEND',
# shells
'CTRIA3',
'CTRIA6',
'CTRIAR',
'CTRIAX',
'CTRIAX6',
'CQUAD4',
'CQUAD8',
'CQUADR',
'CQUADX',
'CQUAD',
# solids
'CTETRA',
'CPENTA',
'CHEXA',
# other
'CSHEAR',
'CVISC',
'CRAC2D',
'CRAC3D',
'CGAP',
'CFAST',
'RBE2',
'RBE3',
# thermal
'CHBDYP',
'CHBDYG',
'CONV',
]
],
]
#print("self.card_count", self.card_count)
if card_count is None:
card_count = self.card_count
for block in blocks:
block_name, keys = block
key_count = 0
for key in sorted(keys):
try:
value = card_count[key]
summary += ' NUMBER OF %-8s %-8s = %8s\n' % (
key, block_name, value)
key_count += 1
except KeyError:
pass
if key_count:
summary += ' \n'
if summary:
f06.write(summary_header)
f06.write(summary)
page_stamp = self.make_stamp(self.title, self.date)
f06.write(page_stamp % self.page_num)
self.page_num += 1
def write_f06(self,
f06_outname,
is_mag_phase=False,
is_sort1=True,
delete_objects=True,
end_flag=False,
quiet=True,
repr_check=False,
close=True):
"""
Writes an F06 file based on the data we have stored in the object
Parameters
----------
f06_outname : str
the name of the F06 file to write
is_mag_phase : bool; default=False
should complex data be written using Magnitude/Phase
instead of Real/Imaginary
Real objects don't use this parameter
is_sort1 : bool; default=True
writes output in SORT1 format if the output is transient;
ignored for static analyses
delete_objects : bool; default=True
should objects be deleted after they're written to reduce memory
end_flag : bool; default=False
should a dummy Nastran "END" table be made
quiet : bool; default=False
suppress print messages
repr_check: bool; defualt=False
calls the object repr as a validation test (prints nothing)
close : bool; default=True
close the f06 file
"""
if not quiet:
print("F06:")
if isinstance(f06_outname, str):
if PY2:
f06 = open(f06_outname, 'wb')
else:
f06 = open(f06_outname, 'w')
self._write_summary(f06)
elif hasattr(f06_outname, 'read') and hasattr(f06_outname, 'write'):
#f06 = f06_outname
#else:
#print('type(f06_outname) =', type(f06_outname))
#assert isinstance(f06_outname, file), 'type(f06_outname)= %s' % f06_outname
f06 = f06_outname
f06_outname = f06.name
if not quiet:
print('f06_outname =', f06_outname)
page_stamp = self.make_stamp(self.title, self.date)
if self.grid_point_weight.reference_point is not None:
if not quiet:
print(" grid_point_weight")
self.page_num = self.grid_point_weight.write_f06(
f06, page_stamp, self.page_num)
if repr_check:
str(self.grid_point_weight)
assert isinstance(self.page_num,
int), self.grid_point_weight.__class__.__name__
if self.oload_resultant is not None:
self.page_num = self.oload_resultant.write_f06(
f06, page_stamp, self.page_num)
if repr_check:
str(self.oload_resultant)
assert isinstance(self.page_num,
int), self.oload_resultant.__class__.__name__
# writes all results for
self._write_f06_subcase_based(f06,
page_stamp,
delete_objects=delete_objects,
is_mag_phase=is_mag_phase,
is_sort1=is_sort1,
quiet=quiet,
repr_check=repr_check)
#self._write_f06_time_based(f06, page_stamp)
self.write_matrices(f06, page_stamp, self.page_num)
f06.write(make_end(end_flag, self.end_options))
if close:
f06.close()
def write_matrices(self, f06, page_stamp, page_num):
if hasattr(self, 'monitor1'):
page_num = self.monitor1.write(f06,
page_stamp=page_stamp,
page_num=page_num)
print('MONPNT1 from [PMRF, PERF, PFRF, AGRF]')
for name, matrix in iteritems(self.matrices):
if name == 'MP3F':
page_num = self.monitor3.write(f06,
page_stamp=page_stamp,
page_num=page_num)
print('MONPNT3 from MP3F')
elif name in ['PMRF', 'PERF', 'PFRF', 'AGRF']:
pass
else:
print(matrix)
matrix.write(f06)
def _write_f06_subcase_based(self,
f06,
page_stamp,
delete_objects=True,
is_mag_phase=False,
is_sort1=True,
quiet=False,
repr_check=False):
"""
Helper function for ``write_f06`` that does the real work
Parameters
----------
f06 : file
the opened file object
page_stamp : ???
???
is_mag_phase : bool; default=False
should complex data be written using Magnitude/Phase
instead of Real/Imaginary
Real objects don't use this parameter
is_sort1 : bool; default=True
writes output in SORT1 format if the output is transient;
ignored for static analyses
delete_objects : bool; default=True
should objects be deleted after they're written to reduce memory
quiet : bool; default=False
suppress print messages
repr_check: bool; defualt=False
calls the object repr as a validation test (prints nothing)
"""
is_failed = False
header = [
' DEFAULT \n',
'\n', ''
]
# eigenvalues are written first
f06.write(page_stamp % self.page_num)
self.page_num += 1
for ikey, result in sorted(iteritems(self.eigenvalues)):
if not quiet:
print('%-18s case=%r' % (result.__class__.__name__, ikey))
self.page_num = result.write_f06(f06,
header,
page_stamp,
page_num=self.page_num)
if repr_check:
str(result)
assert isinstance(self.page_num,
int), 'pageNum=%r' % str(self.page_num)
if delete_objects:
del result
self.page_num += 1
# then eigenvectors
# has a special header
# isubcases = sorted(self.iSubcaseNameMap.keys())
# TODO: superelement version...need the nominal...
res_keys_subcase = self.subcase_key
if len(res_keys_subcase) == 0:
self.log.warning('no cases to write...self.subcase_key=%r' %
self.subcase_key)
return
for isubcase, res_keys in sorted(iteritems(res_keys_subcase)):
for res_key in res_keys:
if isinstance(res_key, tuple):
is_compressed = False
else:
# int
is_compressed = True
isubcase = res_key
if res_key not in self.eigenvectors:
continue
result = self.eigenvectors[res_key]
if repr_check:
str(result)
subtitle = result.subtitle
header[0] = ' %s\n' % subtitle
header[
1] = '0 SUBCASE %i\n' % isubcase
#header[2] = complex/nonlinear
res_length = 18
res_format = '*%%-%is SUBCASE=%%i' % res_length
res_format_vectorized = ' %%-%is SUBCASE=%%i SUBTITLE=%%s' % res_length
class_name = result.__class__.__name__
if hasattr(result, 'data'):
if not quiet:
print(res_format_vectorized %
(class_name, isubcase, subtitle))
else:
print(res_format % (class_name, isubcase))
self.page_num = result.write_f06(f06,
header,
page_stamp,
self.page_num,
is_mag_phase=is_mag_phase,
is_sort1=True)
assert isinstance(self.page_num,
int), 'pageNum=%r' % str(self.page_num)
if delete_objects:
del result
self.page_num += 1
# finally, we writte all the other tables
# nastran puts the tables in order of the Case Control deck,
# but we're lazy so we just hardcode the order
# subcase name, subcase ID, transient word & value
header_old = [
' DEFAULT \n',
'\n', ' \n'
]
header = copy.deepcopy(header_old)
res_types = [
self.displacements,
self.displacements_ROUGV1,
self.displacements_PSD,
self.displacements_ATO,
self.displacements_RMS,
self.displacements_CRM,
self.displacements_NO,
self.displacements_scaled, # ???
self.accelerations,
self.accelerations_ROUGV1,
self.accelerations_PSD,
self.accelerations_ATO,
self.accelerations_RMS,
self.accelerations_CRM,
self.accelerations_NO,
self.velocities,
self.velocities_ROUGV1,
self.velocities_PSD,
self.velocities_ATO,
self.velocities_RMS,
self.velocities_CRM,
self.velocities_NO,
self.force_vectors,
self.load_vectors,
self.temperatures,
#self.eigenvectors,
self.eigenvectors_RADCONS,
self.eigenvectors_RADEFFM,
self.eigenvectors_RADEATC,
self.eigenvectors_ROUGV1,
self.mpc_forces,
self.mpc_forces_PSD,
self.mpc_forces_ATO,
self.mpc_forces_RMS,
self.mpc_forces_RAQCONS,
#self.mpc_forces_RAQEATC,
self.spc_forces,
self.spc_forces_PSD,
self.spc_forces_ATO,
self.spc_forces_RMS,
self.thermal_load_vectors,
#self.strain_energy,
self.cquad4_strain_energy,
self.cquad8_strain_energy,
self.cquadr_strain_energy,
self.cquadx_strain_energy,
self.ctria3_strain_energy,
self.ctria6_strain_energy,
self.ctriar_strain_energy,
self.ctriax_strain_energy,
self.ctriax6_strain_energy,
self.ctetra_strain_energy,
self.cpenta_strain_energy,
self.chexa_strain_energy,
self.cpyram_strain_energy,
self.crod_strain_energy,
self.ctube_strain_energy,
self.conrod_strain_energy,
self.cbar_strain_energy,
self.cbeam_strain_energy,
self.cgap_strain_energy,
self.cbush_strain_energy,
self.celas1_strain_energy,
self.celas2_strain_energy,
self.celas3_strain_energy,
self.celas4_strain_energy,
self.cdum8_strain_energy,
self.dmig_strain_energy,
self.cbend_strain_energy,
self.genel_strain_energy,
self.cshear_strain_energy,
#------------------------------------------
# OEF - forces
# alphabetical order...
# bars
self.cbar_force,
self.cbar_force_ATO,
self.cbar_force_CRM,
self.cbar_force_PSD,
self.cbar_force_RMS,
self.cbar_force_NO,
self.cbar_force_10nodes,
# beam
self.cbend_force,
self.cbeam_force,
self.cbeam_force_ATO,
self.cbeam_force_CRM,
self.cbeam_force_PSD,
self.cbeam_force_RMS,
self.cbeam_force_NO,
# alphabetical
self.celas1_force,
self.celas2_force,
self.celas3_force,
self.celas4_force,
self.cquad4_force,
self.cquad4_force_ATO,
self.cquad4_force_CRM,
self.cquad4_force_PSD,
self.cquad4_force_RMS,
self.cquad4_force_NO,
self.cquad8_force,
self.cquadr_force,
self.conrod_force,
self.crod_force,
self.cshear_force,
self.ctria3_force,
self.ctria6_force,
self.ctriar_force,
self.ctube_force,
# springs
self.celas1_force,
self.celas2_force,
self.celas3_force,
self.celas4_force,
# dampers
self.cdamp1_force,
self.cdamp2_force,
self.cdamp3_force,
self.cdamp4_force,
# other
self.cbush_force,
self.cbush_force_ATO,
self.cbush_force_PSD,
self.cbush_force_CRM,
self.cbush_force_RMS,
self.cbush_force_NO,
self.cgap_force,
self.cvisc_force,
self.chexa_pressure_force,
self.cpenta_pressure_force,
self.ctetra_pressure_force,
self.coneax_force,
#------------------------------------------
# OES - strain
# 1. cbar
# 2. cbeam
# 3. crod/ctube/conrod
# springs
self.celas1_strain,
self.celas2_strain,
self.celas3_strain,
self.celas4_strain,
self.nonlinear_celas1_stress,
self.nonlinear_celas3_stress,
# bars/beams
self.cbar_strain,
self.cbar_strain_10nodes,
self.cbeam_strain,
# plates
self.cquad4_composite_strain,
self.cquad8_composite_strain,
self.cquadr_composite_strain,
self.ctria3_composite_strain,
self.ctria6_composite_strain,
self.ctriar_composite_strain,
self.nonlinear_ctria3_strain,
self.nonlinear_cquad4_strain,
self.ctriax_strain,
# rods
self.nonlinear_crod_strain,
self.nonlinear_ctube_strain,
self.nonlinear_conrod_strain,
self.chexa_strain,
self.conrod_strain,
self.cpenta_strain,
self.cquad4_strain,
self.cquad8_strain,
self.cquadr_strain,
self.crod_strain,
self.cshear_strain,
self.ctetra_strain,
self.ctria3_strain,
self.ctria6_strain,
self.ctriar_strain,
self.ctube_strain,
# bush
self.cbush_strain,
self.nonlinear_cbush_stress,
self.cbush1d_stress_strain,
#------------------------------------------
# cbars/cbeams
self.cbar_stress,
self.cbar_stress_10nodes,
self.nonlinear_cbeam_stress,
self.cbeam_stress,
# bush
self.cbush_stress,
# rods
self.nonlinear_crod_stress,
self.nonlinear_ctube_stress,
self.nonlinear_conrod_stress,
# shear
# OES - stress
self.celas1_stress,
self.celas2_stress,
self.celas3_stress,
self.celas4_stress,
self.chexa_stress,
self.conrod_stress,
self.cpenta_stress,
self.cquad4_stress,
self.cquad8_stress,
self.cquadr_stress,
self.crod_stress,
self.cshear_stress,
self.ctetra_stress,
self.ctria3_stress,
self.ctria6_stress,
self.ctriar_stress,
self.ctube_stress,
self.cquad4_composite_stress,
self.cquad8_composite_stress,
self.cquadr_composite_stress,
self.ctria3_composite_stress,
self.ctria6_composite_stress,
self.ctriar_composite_stress,
self.nonlinear_ctria3_stress,
self.nonlinear_cquad4_stress,
self.ctriax_stress,
self.hyperelastic_cquad4_strain,
#------------------------------------------
#OEF - Fluxes - tCode=4 thermal=1
self.thermalLoad_CONV,
#self.thermalLoad_CHBDY,
self.chbdye_thermal_load,
self.chbdyg_thermal_load,
self.chbdyp_thermal_load,
#self.thermalLoad_1D,
self.crod_thermal_load,
self.cbeam_thermal_load,
self.ctube_thermal_load,
self.conrod_thermal_load,
self.cbar_thermal_load,
self.cbend_thermal_load,
#self.thermalLoad_2D_3D,
self.cquad4_thermal_load,
self.ctriax6_thermal_load,
self.cquad8_thermal_load,
self.ctria3_thermal_load,
self.ctria6_thermal_load,
self.ctetra_thermal_load,
self.cthexa_thermal_load,
self.cpenta_thermal_load,
self.thermalLoad_VU,
self.thermalLoad_VU_3D,
self.thermalLoad_VUBeam,
#------------------------------------------
self.grid_point_stresses,
self.grid_point_volume_stresses,
self.grid_point_forces,
]
for isubcase, res_keys in sorted(iteritems(res_keys_subcase)):
# print(res_keys)
for res_key in res_keys:
if isinstance(res_key, tuple):
is_compressed = False
else:
is_compressed = True
res_length = self._get_result_length(res_types, res_key)
if res_length == 0:
# skipped subcase; no saved results
continue
res_format = '*%%-%is SUBCASE=%%i%%s' % res_length
res_format_vectorized = ' %%-%is SUBCASE=%%i SUBTITLE=%%s %%s' % res_length
for res_type in res_types:
if res_key not in res_type:
continue
result = res_type[res_key]
if repr_check:
str(result)
subtitle = result.subtitle
label = result.label
header = ['', '']
header[0] = ' %-126s\n' % subtitle
header[
1] = '0 %-32s SUBCASE %-15i\n \n' % (
label, isubcase)
if result.nonlinear_factor is not None:
header.append('')
try:
element_name = ''
if hasattr(result, 'element_name'):
element_name = ' - ' + result.element_name
class_name = result.__class__.__name__
if hasattr(result, 'data'):
if not quiet:
print(res_format_vectorized %
(class_name, isubcase, subtitle,
element_name))
else:
print(res_format %
(class_name, isubcase, element_name))
try:
self.page_num = result.write_f06(
f06,
header,
page_stamp,
page_num=self.page_num,
is_mag_phase=is_mag_phase,
is_sort1=is_sort1)
except Exception as e:
print_exc(file=sys.stdout)
is_failed = True
assert isinstance(
self.page_num,
int), 'pageNum=%r' % str(self.page_num)
except:
#print("result name = %r" % result.name())
raise
if delete_objects:
del result
self.page_num += 1
if is_failed:
if PY2:
raise
raise RuntimeError(
'Python 3 cannot reraise Exceptions like this...see previous traceback'
)
class F06Writer(OP2_F06_Common):
def __init__(self):
OP2_F06_Common.__init__(self)
self.card_count = {}
self.additional_matrices = {}
self.matrices = {}
self.subcase_key = defaultdict(list)
self.end_options = {}
self._results = ResultSet(self.get_all_results())
def get_all_results(self):
all_results = ['stress', 'strain', 'element_forces', 'constraint_forces'] + self.get_table_types()
return all_results
def clear_results(self):
self._results.clear()
def add_results(self, results):
if isinstance(results, string_types):
results = [results]
all_results = self.get_all_results()
for result in results:
result = str(result)
if result not in all_results:
raise RuntimeError('%r is not a valid result to remove; all_results=%s' % (result, all_results))
if 'stress' == result:
stress_results = []
for result in all_results:
if 'stress' in result.lower():
stress_results.append(result)
#stress_results = [result if 'stress' in result.lower() for result in all_results]
self._results.update(stress_results)
elif 'strain' == result:
strain_results = []
for result in all_results:
if 'strain' in result.lower():
strain_results.append(result)
#strain_results = [result if 'strain' in result.lower() for result in all_results]
self._results.update(strain_results)
elif 'stress' in result.lower():
self._results.add('stress')
elif 'strain' in result.lower():
self._results.add('strain')
elif result in ('spc_forces', 'mpc_forces', 'constraint_forces'):
self._results.add('constraint_forces')
elif 'force' in result.lower(): # could use more validation...
self._results.add('element_forces')
# thermalLoad_VU_3D, thermalLoad_1D, thermalLoad_CONV, thermalLoad_2D_3D
self._results.add(result)
def set_results(self, results):
if isinstance(results, string_types):
results = [results]
self.clear_results()
self.add_results(results)
def remove_results(self, results):
self._results.remove(results)
def make_f06_header(self):
"""If this class is inherited, the F06 Header may be overwritten"""
return make_f06_header()
def make_stamp(self, Title, today):
"""If this class is inherited, the PAGE stamp may be overwritten"""
return make_stamp(Title, today)
def make_grid_point_singularity_table(self, failed):
msg = ''
if failed:
msg += '0 G R I D P O I N T S I N G U L A R I T Y T A B L E\n'
msg += '0 POINT TYPE FAILED STIFFNESS OLD USET NEW USET\n'
msg += ' ID DIRECTION RATIO EXCLUSIVE UNION EXCLUSIVE UNION\n'
for (nid, dof) in failed:
msg += ' %8s G %s 0.00E+00 B F SB SB *\n' % (nid, dof)
else:
msg += 'No constraints have been applied...\n'
page_stamp = self.make_stamp(self.title, self.date)
msg += page_stamp % self.page_num
self.page_num += 1
return msg
def _write_summary(self, f06, card_count=None):
summary_header = ' M O D E L S U M M A R Y\n\n'
summary = ''
self.cards_to_read = set([
# rigid elements
'RBAR', 'RBAR1', 'RBE1', 'RBE2', 'RBE3',
# spc/mpc constraints
'SPC', 'SPCADD', 'SPC1', 'SPCD', 'SPCAX',
'MPC', 'MPCADD',
'SUPORT', 'SUPORT1',
# aero cards
'CAERO1', 'CAERO2', 'CAERO3', 'CAERO4', 'CAERO5',
# temperature cards
'CHBDYE', 'CHBDYG', 'CHBDYP',
'CONV',
])
blocks = [
['POINTS', ['GRID', 'GRDSET', ]],
['ENTRIES', ['SPOINT']],
['ELEMENTS',
[
# these are sorted
# elements
'CONM1', 'CONM2', 'CMASS1', 'CMASS2', 'CMASS3', 'CMASS4',
# springs
'CELAS1', 'CELAS2', 'CELAS3', 'CELAS4', 'CELAS5',
# bushings
'CBUSH', 'CBUSH1D', 'CBUSH2D',
# dampers
'CDAMP1', 'CDAMP2', 'CDAMP3', 'CDAMP4', 'CDAMP5',
# bar flags
'BAROR', 'CBARAO',
# bars
'CBAR', 'CROD', 'CTUBE', 'CBEAM', 'CBEAM3', 'CONROD', 'CBEND',
# shells
'CTRIA3', 'CTRIA6', 'CTRIAR', 'CTRIAX', 'CTRIAX6',
'CQUAD4', 'CQUAD8', 'CQUADR', 'CQUADX', 'CQUAD',
# solids
'CTETRA', 'CPENTA', 'CHEXA',
# other
'CSHEAR', 'CVISC', 'CRAC2D', 'CRAC3D',
'CGAP', 'CFAST', 'RBE2', 'RBE3',
# thermal
'CHBDYP', 'CHBDYG', 'CONV',
]],
]
#print("self.card_count", self.card_count)
if card_count is None:
card_count = self.card_count
for block in blocks:
block_name, keys = block
key_count = 0
for key in sorted(keys):
try:
value = card_count[key]
summary += ' NUMBER OF %-8s %-8s = %8s\n' % (key, block_name, value)
key_count += 1
except KeyError:
pass
if key_count:
summary += ' \n'
if summary:
f06.write(summary_header)
f06.write(summary)
page_stamp = self.make_stamp(self.title, self.date)
f06.write(page_stamp % self.page_num)
self.page_num += 1
def write_f06(self, f06_outname, is_mag_phase=False, is_sort1=True,
delete_objects=True, end_flag=False, quiet=True, repr_check=False,
close=True):
"""
Writes an F06 file based on the data we have stored in the object
Parameters
----------
f06_outname : str
the name of the F06 file to write
is_mag_phase : bool; default=False
should complex data be written using Magnitude/Phase
instead of Real/Imaginary
Real objects don't use this parameter
is_sort1 : bool; default=True
writes output in SORT1 format if the output is transient;
ignored for static analyses
delete_objects : bool; default=True
should objects be deleted after they're written to reduce memory
end_flag : bool; default=False
should a dummy Nastran "END" table be made
quiet : bool; default=False
suppress print messages
repr_check: bool; defualt=False
calls the object repr as a validation test (prints nothing)
close : bool; default=True
close the f06 file
"""
if not quiet:
print("F06:")
if isinstance(f06_outname, str):
if PY2:
f06 = open(f06_outname, 'wb')
else:
f06 = open(f06_outname, 'w')
self._write_summary(f06)
elif hasattr(f06_outname, 'read') and hasattr(f06_outname, 'write'):
#f06 = f06_outname
#else:
#print('type(f06_outname) =', type(f06_outname))
#assert isinstance(f06_outname, file), 'type(f06_outname)= %s' % f06_outname
f06 = f06_outname
f06_outname = f06.name
if not quiet:
print('f06_outname =', f06_outname)
page_stamp = self.make_stamp(self.title, self.date)
if self.grid_point_weight.reference_point is not None:
if not quiet:
print(" grid_point_weight")
self.page_num = self.grid_point_weight.write_f06(f06, page_stamp, self.page_num)
if repr_check:
str(self.grid_point_weight)
assert isinstance(self.page_num, int), self.grid_point_weight.__class__.__name__
if self.oload_resultant is not None:
self.page_num = self.oload_resultant.write_f06(f06, page_stamp, self.page_num)
if repr_check:
str(self.oload_resultant)
assert isinstance(self.page_num, int), self.oload_resultant.__class__.__name__
# writes all results for
self._write_f06_subcase_based(f06, page_stamp, delete_objects=delete_objects,
is_mag_phase=is_mag_phase, is_sort1=is_sort1,
quiet=quiet, repr_check=repr_check)
#self._write_f06_time_based(f06, page_stamp)
f06.write(make_end(end_flag, self.end_options))
if close:
f06.close()
def _write_f06_subcase_based(self, f06, page_stamp, delete_objects=True,
is_mag_phase=False, is_sort1=True, quiet=False,
repr_check=False):
"""
Helper function for ``write_f06`` that does the real work
Parameters
----------
f06 : file
the opened file object
page_stamp : ???
???
is_mag_phase : bool; default=False
should complex data be written using Magnitude/Phase
instead of Real/Imaginary
Real objects don't use this parameter
is_sort1 : bool; default=True
writes output in SORT1 format if the output is transient;
ignored for static analyses
delete_objects : bool; default=True
should objects be deleted after they're written to reduce memory
quiet : bool; default=False
suppress print messages
repr_check: bool; defualt=False
calls the object repr as a validation test (prints nothing)
"""
is_failed = False
header = [' DEFAULT \n',
'\n', '']
# eigenvalues are written first
f06.write(page_stamp % self.page_num)
self.page_num += 1
for ikey, result in sorted(iteritems(self.eigenvalues)):
if not quiet:
print('%-18s case=%r' % (result.__class__.__name__, ikey))
self.page_num = result.write_f06(f06, header, page_stamp,
page_num=self.page_num)
if repr_check:
str(result)
assert isinstance(self.page_num, int), 'pageNum=%r' % str(self.page_num)
if delete_objects:
del result
self.page_num += 1
# then eigenvectors
# has a special header
# isubcases = sorted(self.iSubcaseNameMap.keys())
# TODO: superelement version...need the nominal...
res_keys_subcase = self.subcase_key
if len(res_keys_subcase) == 0:
self.log.warning('no cases to write...self.subcase_key=%r' % self.subcase_key)
return
for isubcase, res_keys in sorted(iteritems(res_keys_subcase)):
for res_key in res_keys:
if isinstance(res_key, tuple):
is_compressed = False
else:
# int
is_compressed = True
isubcase = res_key
if res_key not in self.eigenvectors:
continue
result = self.eigenvectors[res_key]
if repr_check:
str(result)
subtitle = result.subtitle
header[0] = ' %s\n' % subtitle
header[1] = '0 SUBCASE %i\n' % isubcase
#header[2] = complex/nonlinear
res_length = 18
res_format = '*%%-%is SUBCASE=%%i' % res_length
res_format_vectorized = ' %%-%is SUBCASE=%%i SUBTITLE=%%s' % res_length
class_name = result.__class__.__name__
if hasattr(result, 'data'):
if not quiet:
print(res_format_vectorized % (class_name, isubcase, subtitle))
else:
print(res_format % (class_name, isubcase))
self.page_num = result.write_f06(f06, header, page_stamp,
self.page_num, is_mag_phase=is_mag_phase, is_sort1=True)
assert isinstance(self.page_num, int), 'pageNum=%r' % str(self.page_num)
if delete_objects:
del result
self.page_num += 1
# finally, we writte all the other tables
# nastran puts the tables in order of the Case Control deck,
# but we're lazy so we just hardcode the order
# subcase name, subcase ID, transient word & value
header_old = [' DEFAULT \n',
'\n', ' \n']
header = copy.deepcopy(header_old)
res_types = [
self.displacements, self.displacements_ROUGV1,
self.displacements_PSD, self.displacements_ATO, self.displacements_RMS,
self.displacements_CRM, self.displacements_NO,
self.displacements_scaled, # ???
self.accelerations, self.accelerations_ROUGV1,
self.accelerations_PSD, self.accelerations_ATO, self.accelerations_RMS,
self.accelerations_CRM, self.accelerations_NO,
self.velocities, self.velocities_ROUGV1,
self.velocities_PSD, self.velocities_ATO, self.velocities_RMS,
self.velocities_CRM, self.velocities_NO,
self.force_vectors,
self.load_vectors,
self.temperatures,
#self.eigenvectors,
self.eigenvectors_RADCONS,
self.eigenvectors_RADEFFM,
self.eigenvectors_RADEATC,
self.eigenvectors_ROUGV1,
self.mpc_forces, self.mpc_forces_PSD, self.mpc_forces_ATO, self.mpc_forces_RMS,
self.mpc_forces_RAQCONS,
#self.mpc_forces_RAQEATC,
self.spc_forces, self.spc_forces_PSD, self.spc_forces_ATO, self.spc_forces_RMS,
self.thermal_load_vectors,
#self.strain_energy,
self.cquad4_strain_energy, self.cquad8_strain_energy,
self.cquadr_strain_energy, self.cquadx_strain_energy,
self.ctria3_strain_energy, self.ctria6_strain_energy,
self.ctriar_strain_energy, self.ctriax_strain_energy,
self.ctriax6_strain_energy,
self.ctetra_strain_energy, self.cpenta_strain_energy,
self.chexa_strain_energy, self.cpyram_strain_energy,
self.crod_strain_energy, self.ctube_strain_energy,
self.conrod_strain_energy,
self.cbar_strain_energy, self.cbeam_strain_energy,
self.cgap_strain_energy, self.cbush_strain_energy,
self.celas1_strain_energy, self.celas2_strain_energy,
self.celas3_strain_energy, self.celas4_strain_energy,
self.cdum8_strain_energy, self.dmig_strain_energy,
self.cbend_strain_energy,
self.genel_strain_energy, self.cshear_strain_energy,
#------------------------------------------
# OEF - forces
# alphabetical order...
# bars
self.cbar_force, self.cbar_force_ATO, self.cbar_force_CRM, self.cbar_force_PSD, self.cbar_force_RMS, self.cbar_force_NO,
self.cbar_force_10nodes,
# beam
self.cbend_force,
self.cbeam_force, self.cbeam_force_ATO, self.cbeam_force_CRM, self.cbeam_force_PSD, self.cbeam_force_RMS, self.cbeam_force_NO,
# alphabetical
self.celas1_force,
self.celas2_force,
self.celas3_force,
self.celas4_force,
self.cquad4_force, self.cquad4_force_ATO, self.cquad4_force_CRM, self.cquad4_force_PSD, self.cquad4_force_RMS, self.cquad4_force_NO,
self.cquad8_force,
self.cquadr_force,
self.conrod_force,
self.crod_force,
self.cshear_force,
self.ctria3_force,
self.ctria6_force,
self.ctriar_force,
self.ctube_force,
# springs
self.celas1_force,
self.celas2_force,
self.celas3_force,
self.celas4_force,
# dampers
self.cdamp1_force,
self.cdamp2_force,
self.cdamp3_force,
self.cdamp4_force,
# other
self.cbush_force, self.cbush_force_ATO, self.cbush_force_PSD, self.cbush_force_CRM, self.cbush_force_RMS, self.cbush_force_NO,
self.cgap_force,
self.cvisc_force,
self.chexa_pressure_force,
self.cpenta_pressure_force,
self.ctetra_pressure_force,
self.coneax_force,
#------------------------------------------
# OES - strain
# 1. cbar
# 2. cbeam
# 3. crod/ctube/conrod
# springs
self.celas1_strain,
self.celas2_strain,
self.celas3_strain,
self.celas4_strain,
self.nonlinear_celas1_stress,
self.nonlinear_celas3_stress,
# bars/beams
self.cbar_strain,
self.cbar_strain_10nodes,
self.cbeam_strain,
# plates
self.cquad4_composite_strain,
self.cquad8_composite_strain,
self.cquadr_composite_strain,
self.ctria3_composite_strain,
self.ctria6_composite_strain,
self.ctriar_composite_strain,
self.nonlinear_ctria3_strain,
self.nonlinear_cquad4_strain,
self.ctriax_strain,
# rods
self.nonlinear_crod_strain,
self.nonlinear_ctube_strain,
self.nonlinear_conrod_strain,
self.chexa_strain,
self.conrod_strain,
self.cpenta_strain,
self.cquad4_strain,
self.cquad8_strain,
self.cquadr_strain,
self.crod_strain,
self.cshear_strain,
self.ctetra_strain,
self.ctria3_strain,
self.ctria6_strain,
self.ctriar_strain,
self.ctube_strain,
# bush
self.cbush_strain,
self.nonlinear_cbush_stress,
self.cbush1d_stress_strain,
#------------------------------------------
# cbars/cbeams
self.cbar_stress,
self.cbar_stress_10nodes,
self.nonlinear_cbeam_stress,
self.cbeam_stress,
# bush
self.cbush_stress,
# rods
self.nonlinear_crod_stress,
self.nonlinear_ctube_stress,
self.nonlinear_conrod_stress,
# shear
# OES - stress
self.celas1_stress,
self.celas2_stress,
self.celas3_stress,
self.celas4_stress,
self.chexa_stress,
self.conrod_stress,
self.cpenta_stress,
self.cquad4_stress,
self.cquad8_stress,
self.cquadr_stress,
self.crod_stress,
self.cshear_stress,
self.ctetra_stress,
self.ctria3_stress,
self.ctria6_stress,
self.ctriar_stress,
self.ctube_stress,
self.cquad4_composite_stress,
self.cquad8_composite_stress,
self.cquadr_composite_stress,
self.ctria3_composite_stress,
self.ctria6_composite_stress,
self.ctriar_composite_stress,
self.nonlinear_ctria3_stress,
self.nonlinear_cquad4_stress,
self.ctriax_stress,
self.hyperelastic_cquad4_strain,
#------------------------------------------
#OEF - Fluxes - tCode=4 thermal=1
self.thermalLoad_CONV,
#self.thermalLoad_CHBDY,
self.chbdye_thermal_load,
self.chbdyg_thermal_load,
self.chbdyp_thermal_load,
#self.thermalLoad_1D,
self.crod_thermal_load,
self.cbeam_thermal_load,
self.ctube_thermal_load,
self.conrod_thermal_load,
self.cbar_thermal_load,
self.cbend_thermal_load,
#self.thermalLoad_2D_3D,
self.cquad4_thermal_load,
self.ctriax6_thermal_load,
self.cquad8_thermal_load,
self.ctria3_thermal_load,
self.ctria6_thermal_load,
self.ctetra_thermal_load,
self.cthexa_thermal_load,
self.cpenta_thermal_load,
self.thermalLoad_VU,
self.thermalLoad_VU_3D,
self.thermalLoad_VUBeam,
#------------------------------------------
self.grid_point_stresses, self.grid_point_volume_stresses, self.grid_point_forces,
]
for isubcase, res_keys in sorted(iteritems(res_keys_subcase)):
# print(res_keys)
for res_key in res_keys:
if isinstance(res_key, tuple):
is_compressed = False
else:
is_compressed = True
res_length = self._get_result_length(res_types, res_key)
if res_length == 0:
# skipped subcase; no saved results
continue
res_format = '*%%-%is SUBCASE=%%i%%s' % res_length
res_format_vectorized = ' %%-%is SUBCASE=%%i SUBTITLE=%%s %%s' % res_length
for res_type in res_types:
if res_key not in res_type:
continue
result = res_type[res_key]
if repr_check:
str(result)
subtitle = result.subtitle
label = result.label
header = ['', '']
header[0] = ' %-126s\n' % subtitle
header[1] = '0 %-32s SUBCASE %-15i\n \n' % (label, isubcase)
if result.nonlinear_factor is not None:
header.append('')
try:
element_name = ''
if hasattr(result, 'element_name'):
element_name = ' - ' + result.element_name
class_name = result.__class__.__name__
if hasattr(result, 'data'):
if not quiet:
print(res_format_vectorized % (class_name, isubcase, subtitle, element_name))
else:
print(res_format % (class_name, isubcase, element_name))
try:
self.page_num = result.write_f06(f06, header, page_stamp, page_num=self.page_num,
is_mag_phase=is_mag_phase, is_sort1=is_sort1)
except Exception as e:
print_exc(file=sys.stdout)
is_failed = True
assert isinstance(self.page_num, int), 'pageNum=%r' % str(self.page_num)
except:
#print("result name = %r" % result.name())
raise
if delete_objects:
del result
self.page_num += 1
if is_failed:
if PY2:
raise
raise RuntimeError('Python 3 cannot reraise Exceptions like this...see previous traceback')
