class TestCyclicResultReader(object):
# avoids errors in collection
result_file = os.path.join(data_path, 'cyclic_%s.rst' % rver)
try:
# test if raw results are being read properly by using the normal result reader
result = ResultFile(result_file, ignore_cyclic=True)
if rver == 'v182':
ansys = pyansys.ANSYS('/usr/ansys_inc/v182/ansys/bin/ansys182',
override=True, jobname=rver,
loglevel='DEBUG',
interactive_plotting=False,
prefer_pexpect=True)
else:
ansys = pyansys.ANSYS('/usr/ansys_inc/v150/ansys/bin/ansys150',
override=True, jobname=rver,
loglevel='DEBUG',
interactive_plotting=False,
prefer_pexpect=True)
# copy result file to ansys's temporary path
copyfile(result_file, os.path.join(ansys.path, '%s.rst' % rver))
# setup ansys for output without line breaks
ansys.Post1()
ansys.Set(1, 1)
ansys.Header('OFF', 'OFF', 'OFF', 'OFF', 'OFF', 'OFF')
nsigfig = 10
ansys.Format('', 'E', nsigfig + 9, nsigfig)
ansys.Page(1E9, '', -1, 240)
except: # for travis and appveyor
pass
def test_prnsol_u(self):
# verify cyclic displacements
table = self.ansys.Prnsol('u').splitlines()
if self.ansys.using_corba:
array = np.genfromtxt(table[7:])
else:
array = np.genfromtxt(table[9:])
ansys_nnum = array[:, 0].astype(np.int)
ansys_disp = array[:, 1:-1]
nnum, disp = self.result.nodal_solution(0)
# cyclic model will only output the master sector
ansys_nnum = ansys_nnum[:nnum.size]
ansys_disp = ansys_disp[:nnum.size]
assert np.allclose(ansys_nnum, nnum)
assert np.allclose(ansys_disp, disp)
def test_presol_s(self):
# verify element stress
element_stress, _, enode = self.result.element_stress(0)
element_stress = np.vstack(element_stress)
enode = np.hstack(enode)
# parse ansys result
table = self.ansys.Presol('S').splitlines()
ansys_element_stress = []
line_length = len(table[15])
for line in table:
if len(line) == line_length:
ansys_element_stress.append(line)
ansys_element_stress = np.genfromtxt(ansys_element_stress)
ansys_enode = ansys_element_stress[:, 0].astype(np.int)
ansys_element_stress = ansys_element_stress[:, 1:]
assert np.allclose(element_stress, ansys_element_stress)
assert np.allclose(enode, ansys_enode)
def test_prnsol_s(self):
# verify cyclic displacements
table = self.ansys.Prnsol('s').splitlines()
if self.ansys.using_corba:
array = np.genfromtxt(table[7:])
else:
array = np.genfromtxt(table[10:])
ansys_nnum = array[:, 0].astype(np.int)
ansys_stress = array[:, 1:]
nnum, stress = self.result.nodal_stress(0)
# v150 includes nodes in the geometry that aren't in the result
mask = np.in1d(nnum, ansys_nnum)
nnum = nnum[mask]
stress = stress[mask]
assert np.allclose(ansys_nnum, nnum)
assert np.allclose(ansys_stress, stress)
def test_prnsol_prin(self):
# verify principal stress
table = self.ansys.Prnsol('prin').splitlines()
if self.ansys.using_corba:
array = np.genfromtxt(table[7:])
else:
array = np.genfromtxt(table[10:])
ansys_nnum = array[:, 0].astype(np.int)
ansys_stress = array[:, 1:]
nnum, stress = self.result.principal_nodal_stress(0)
# v150 includes nodes in the geometry that aren't in the result
mask = np.in1d(nnum, ansys_nnum)
nnum = nnum[mask]
stress = stress[mask]
assert np.allclose(ansys_nnum, nnum)
assert np.allclose(ansys_stress, stress, atol=1E-2)
@pytest.mark.skipif(not system_supports_plotting(), reason="Requires active X Server")
def test_plot(self):
filename = '/tmp/temp.png'
self.result.plot_nodal_solution(0, screenshot=filename,
off_screen=True)
# self.result.plot_nodal_stress(0, 'Sx', screenshot=filename,
# off_screen=True)
self.result.plot_principal_nodal_stress(0, 'EQV', screenshot=filename,
off_screen=True)
def test_exit(self):
self.ansys.Exit()
def read_binary(filename, **kwargs):
"""
Reads ANSYS-written binary files:
- Jobname.RST: result file from structural analysis
- Jobname.EMAT: Stores data related to element matrices
- Jobname.FULL Stores the full stiffness-mass matrix
Parameters
----------
filename : str
Filename to read.
**kwargs : keyword arguments
See the individual classes for additional keyword arguments.
Examples
--------
>>> import pyansys
>>> result = pyansys.read_binary('file.rst')
>>> full_file = pyansys.read_binary('file.full')
>>> emat_file = pyansys.read_binary('file.emat')
Notes
-----
The following file types are unsupported
- Jobname.DSUB file, storing displacements related to substructure
matrices
- Jobname.SUB file, storing data related to substructure matrices
- Jobname.RFRQ file, storing data related to a mode-superposition
harmonic analysis
- The Jobname.RDSP file, storing data related to a
mode-superposition transient analysis.
- Jobname.MODE file, storing data related to a modal analysis
- Jobname.RMG A magnetic analysis
- Jobname.RFL A FLOTRAN analysis (a legacy results file)
- Jobname.RTH A thermal analysis
"""
if not os.path.isfile(filename):
raise FileNotFoundError('%s is not a file or cannot be found' %
str(filename))
# return BinaryFile(filename)
file_format = read_standard_header(filename)['file format']
if file_format == 2:
from pyansys.emat import EmatFile
return EmatFile(filename, **kwargs)
elif file_format == 4:
from pyansys.full import FullFile
return FullFile(filename, **kwargs)
elif file_format == 12:
from pyansys.rst import ResultFile
result = ResultFile(filename, **kwargs)
# check if it's a cyclic result file
ignore_cyclic = kwargs.pop('ignore_cyclic', False)
if result.header['nSector'] != 1 and not ignore_cyclic:
from pyansys.cyclic_reader import CyclicResult
return CyclicResult(filename)
return result
elif file_format == 16:
from pyansys.db import Database
return Database(filename, debug=kwargs.pop('debug', False))
else:
if file_format in ANSYS_BINARY_FILE_TYPES:
file_type = ANSYS_BINARY_FILE_TYPES[file_format]
else:
file_type = str(file_format)
raise RuntimeError('ANSYS binary "%s" not supported' % file_type)