def test_panair_io_02(self):
model = PanairGrid(log=None, debug=False)
msg = (
'$circular sections - nacelle with composite panels\n'
# =kn\n'
'1.\n'
# =kt\n'
'1.\n'
# =nopt netname\n'
'0. cowlu\n'
# =nm\n'
'20.\n'
# =xs(1) ri(1) xs(2) ri(2) xs(*) ri(*)\n'
' 2.0000 2.3000 1.5756 2.3000 1.1486 2.3000\n'
' 0.7460 2.3030 0.4069 2.3286 0.1624 2.3790\n'
' 0.0214 2.4542 -0.0200 2.5485 0.0388 2.6522\n'
' 0.2056 2.7554 0.4869 2.8522 0.8883 2.9413\n'
' 1.4250 3.0178 2.1188 3.0656 2.9586 3.0658\n'
' 3.8551 3.0175 4.6715 2.9439 5.3492 2.8700\n'
' 6.0000 2.7842 6.4687 2.7442\n'
# =nn\n'
'5.\n'
# =th(1) th(2) th(3) th(4) th(5)\n'
'-90. -45. 0. 45. 90.\n')
section = msg.split('\n')
model._read_circular_section(section)
model.write_panair('junk_circ.inp')
os.remove('junk_circ.inp')
def test_panair_io_03(self):
"""tests the SWB model"""
log = get_logger(level='warning')
in_filename = os.path.join(TEST_PATH, 'SWB', 'SWB.inp')
#out_filename = os.path.join(TEST_PATH, 'M100', 'M100_out.inp')
model = PanairGrid(log=log, debug=False)
model.read_panair(in_filename)
(points, elements, regions, kt, cp_nrom) = model.get_points_elements_regions()
model.write_panair('junk_swb.inp')
os.remove('junk_swb.inp')
def test_panair_io_01(self):
in_filename = os.path.join(test_path, 'M100', 'M100.inp')
out_filename = os.path.join(test_path, 'M100', 'M100_out.inp')
#with open(infile_name, 'w') as f:
#f.write(lines)
model = PanairGrid(log=None, debug=False)
model.read_panair(in_filename)
#assert len(cart3d.points) == 7, 'npoints=%s' % len(cart3d.points)
#assert len(cart3d.elements) == 6, 'nelements=%s' % len(cart3d.elements)
#assert len(cart3d.regions) == 6, 'nregions=%s' % len(cart3d.regions)
#assert len(cart3d.loads) == 0, 'nloads=%s' % len(cart3d.loads)
model.write_panair(out_filename)
(points, elements, regions, kt, cp_nrom) = model.get_points_elements_regions()
os.remove(out_filename)
def test_panair_io_01(self):
log = get_logger(level='warning')
in_filename = os.path.join(test_path, 'M100', 'M100.inp')
out_filename = os.path.join(test_path, 'M100', 'M100_out.inp')
#with open(infile_name, 'w') as f:
#f.write(lines)
model = PanairGrid(log=log, debug=False)
model.read_panair(in_filename)
#assert len(cart3d.points) == 7, 'npoints=%s' % len(cart3d.points)
#assert len(cart3d.elements) == 6, 'nelements=%s' % len(cart3d.elements)
#assert len(cart3d.regions) == 6, 'nregions=%s' % len(cart3d.regions)
#assert len(cart3d.loads) == 0, 'nloads=%s' % len(cart3d.loads)
model.write_panair(out_filename)
(points, elements, regions, kt,
cp_nrom) = model.get_points_elements_regions()
os.remove(out_filename)
def test_panair_io_01(self):
"""test the M100 model"""
log = get_logger(level='warning')
in_filename = os.path.join(TEST_PATH, 'M100', 'M100.inp')
out_filename = os.path.join(TEST_PATH, 'M100', 'M100_out.inp')
#with open(infile_name, 'w') as f:
#f.write(lines)
model = PanairGrid(log=log, debug=False)
model.read_panair(in_filename)
model.write_panair(out_filename)
(points, elements, regions, kt, cp_nrom) = model.get_points_elements_regions()
model.write_panair('junk_m100.inp')
os.remove('junk_m100.inp')
model.print_options()
model.print_abutments()
model.print_grid_summary()
model.print_out_header()
os.remove(out_filename)
def write_panair(self, panair_filename, panair_case_filename):
pan = PanairGrid()
pan.mach = 0.5
pan.is_end = True
pan.ncases = 2
pan.alphas = [0., 5.]
i = 0
pan.nNetworks = 1
kt = 1
cp_norm = 1
for name, comps in sorted(self.components.items()):
#panair_file.write('$ name = %r\n' % name)
for comp in comps:
str(comp)
namei = name + str(i)
# lifting_surface_xyz, lifting_surface_nx, lifting_surface_ny
x = deepcopy(comp.xyz[:, 0])
y = deepcopy(comp.xyz[:, 1])
z = deepcopy(comp.xyz[:, 2])
x = x.reshape((comp.nx, comp.ny))
y = y.reshape((comp.nx, comp.ny))
z = z.reshape((comp.nx, comp.ny))
xyz = np.dstack([x, y, z])
assert xyz.shape[2] == 3
patch = PanairPatch(pan.nNetworks, namei, kt, cp_norm, xyz, self.log)
pan.patches[i] = patch
pan.nNetworks += 1
i += 1
if 'wing' in name.lower(): # make a wing cap
namei = 'cap%i' % i
#assert comp.lifting_surface_nx == 6, comp.lifting_surface_nx
assert comp.ny == 33, comp.ny
#print(x.shape)
xend = deepcopy(x[-1, :])
#print(xend)
yend = deepcopy(y[-1, :])
zend = deepcopy(z[-1, :])
imid = comp.ny // 2
x = np.zeros((imid+1, 2), dtype='float32')
y = np.zeros((imid+1, 2), dtype='float32')
z = np.zeros((imid+1, 2), dtype='float32')
#print(imid, xend[imid], xend.min())
xflip = list(xend[0:imid+1])
yflip = list(yend[0:imid+1])
zflip = list(zend[0:imid+1])
x[:, 0] = xflip[::-1]
y[:, 0] = yflip[::-1]
z[:, 0] = zflip[::-1]
x[:, 1] = xend[imid:]
y[:, 1] = yend[imid:]
z[:, 1] = zend[imid:]
#print(x)
#x = xend[0:imid:-1].extend(x[imid:])
#y = yend[0:imid:-1].extend(y[imid:])
#z = zend[0:imid:-1].extend(z[imid:])
#print(x)
x = x.reshape((2, imid+1))
y = y.reshape((2, imid+1))
z = z.reshape((2, imid+1))
#print(xend)
xyz = np.dstack([x, y, z])
assert xyz.shape[2] == 3
patch = PanairPatch(pan.nNetworks, namei, kt, cp_norm,
xyz, self.log)
pan.patches[i] = patch
pan.nNetworks += 1
i += 1
#i += 1
pan.write_panair(panair_filename)
def write_panair(self, panair_filename, panair_case_filename):
#panair_file = open(panair_filename, 'wb')
pan = PanairGrid()
pan.mach = 0.5
pan.isEnd = True
pan.ncases = 2
pan.alphas = [0., 5.]
i = 0
pan.nNetworks = 1
kt = 1
cpNorm = 1
for name, comps in sorted(iteritems(self.components)):
#panair_file.write('$ name = %r\n' % name)
for comp in comps:
namei = name + str(i)
x = deepcopy(comp.lifting_surface_xyz[:, 0])
y = deepcopy(comp.lifting_surface_xyz[:, 1])
z = deepcopy(comp.lifting_surface_xyz[:, 2])
x = x.reshape((comp.lifting_surface_nx, comp.lifting_surface_ny))
y = y.reshape((comp.lifting_surface_nx, comp.lifting_surface_ny))
z = z.reshape((comp.lifting_surface_nx, comp.lifting_surface_ny))
patch = PanairPatch(pan.nNetworks, namei, kt, cpNorm, x, y, z, self.log)
pan.patches[i] = patch
pan.nNetworks += 1
i += 1
if 'wing' in name.lower(): # make a wing cap
namei = 'cap%i' % i
#assert comp.lifting_surface_nx == 6, comp.lifting_surface_nx
assert comp.lifting_surface_ny == 33, comp.lifting_surface_ny
#print(x.shape)
xend = deepcopy(x[-1, :])
print(xend)
yend = deepcopy(y[-1, :])
zend = deepcopy(z[-1, :])
imid = comp.lifting_surface_ny // 2
x = np.zeros((imid+1, 2), dtype='float32')
y = np.zeros((imid+1, 2), dtype='float32')
z = np.zeros((imid+1, 2), dtype='float32')
print(imid, xend[imid], xend.min())
xflip = list(xend[0:imid+1])
yflip = list(yend[0:imid+1])
zflip = list(zend[0:imid+1])
x[:, 0] = xflip[::-1]
y[:, 0] = yflip[::-1]
z[:, 0] = zflip[::-1]
x[:, 1] = xend[imid:]
y[:, 1] = yend[imid:]
z[:, 1] = zend[imid:]
print (x)
#x = xend[0:imid:-1].extend(x[imid:])
#y = yend[0:imid:-1].extend(y[imid:])
#z = zend[0:imid:-1].extend(z[imid:])
#print(x)
x = x.reshape((2, imid+1))
y = y.reshape((2, imid+1))
z = z.reshape((2, imid+1))
#print(xend)
patch = PanairPatch(pan.nNetworks, namei, kt, cpNorm, x.T, y.T, z.T, self.log)
pan.patches[i] = patch
pan.nNetworks += 1
i += 1
#i += 1
pan.write_panair(panair_filename)