def test_solid_03(self):
"""checks linear static solid material"""
mid = 2
pid = 4
rho = 0.1
tableID = 42
cards = [
#$ Solid Nodes
['GRID', 11, 0, 0., 0., 0., 0],
['GRID', 12, 0, 1., 0., 0., 0],
['GRID', 13, 0, 1., 1., 0., 0],
['GRID', 14, 0, 0., 1., 0., 0],
['GRID', 15, 0, 0., 0., 2., 0],
['GRID', 16, 0, 1., 0., 2., 0],
['GRID', 17, 0, 1., 1., 2., 0],
['GRID', 18, 0, 0., 1., 2., 0],
# Solids
['CHEXA', 7, pid, 11, 12, 13, 14, 15, 16, 17, 18],
['CTETRA', 8, pid, 11, 12, 13, 15],
# Solid Nodes
['GRID', 21, 0, 0., 0., 0., 0,],
['GRID', 22, 0, 1., 0., 0., 0,],
['GRID', 23, 0, 1., 1., 0., 0,],
['GRID', 24, 0, 0., 0., 2., 0,],
['GRID', 25, 0, 1., 0., 2., 0,],
['GRID', 26, 0, 1., 1., 2., 0,],
['CPENTA', 9, pid, 21, 22, 23, 24, 25, 26],
# static
['PSOLID', pid, mid, 0],
['MAT1', mid, 1.0, 2.0, 3.0, rho],
['MATS1', mid, tableID, 'PLASTIC', 0.0, 1, 1, 100000., ],
#['TABLEST'],
['TABLES1', tableID, 1, None, None, None, None, None, None,
1.0, 10.0, 2.0, 10.0, 'ENDT'],
]
card_count = {
'GRID' : 14,
'CPENTA6': 1,
'CTETRA4': 1,
'PSOLID': 1,
'CHEXA8' : 1,
'MAT1': 1,
'MATS1': 1,
'TABLES1': 1,
}
model = BDF(debug=False)
model.allocate(card_count)
for fields in cards:
model.add_card(fields, fields[0], is_list=True)
model.cross_reference()
mat = model.materials[mid]
print('----MAT----', type(mat))
print(mat)
print('E = %s' % mat.get_E_by_material_index())
print('E = %s' % mat.get_E_by_material_id())
def test_solid_02(self):
mid = 2
pid = 4
rho = 0.1
cards = [
#$ Solid Nodes
['GRID', 11, 0, 0., 0., 0., 0],
['GRID', 12, 0, 1., 0., 0., 0],
['GRID', 13, 0, 1., 1., 0., 0],
['GRID', 14, 0, 0., 1., 0., 0],
['GRID', 15, 0, 0., 0., 2., 0],
['GRID', 16, 0, 1., 0., 2., 0],
['GRID', 17, 0, 1., 1., 2., 0],
['GRID', 18, 0, 0., 1., 2., 0],
# Solids
['CHEXA', 7, pid, 11, 12, 13, 14, 15, 16, 17, 18],
['CTETRA', 8, pid, 11, 12, 13, 15],
# Solid Nodes
['GRID', 21, 0, 0., 0., 0., 0,],
['GRID', 22, 0, 1., 0., 0., 0,],
['GRID', 23, 0, 1., 1., 0., 0,],
['GRID', 24, 0, 0., 0., 2., 0,],
['GRID', 25, 0, 1., 0., 2., 0,],
['GRID', 26, 0, 1., 1., 2., 0,],
['CPENTA', 9, pid, 21, 22, 23, 24, 25, 26],
# hyperelastic
['PLSOLID', pid, mid, 'GRID'],
['MATHP', mid, None, None, None, rho],
]
card_count = {
'GRID' : 14,
'CTETRA4': 1,
'CPENTA6': 1,
'CHEXA8' : 1,
'PLSOLID': 1,
'MATHP': 1,
}
model = BDF(debug=True)
model.allocate(card_count)
for fields in cards:
model.add_card(fields, fields[0], is_list=True)
model.cross_reference()
# CTETRA
eid = 8
nsm = 0.
V = 1. / 3.
self.check_solid(model, eid, 'CTETRA4', pid, 'PLSOLID', mid, 'MATHP', nsm, rho, V)
eid = 9
V = 1.0
self.check_solid(model, eid, 'CPENTA6', pid, 'PLSOLID', mid, 'MATHP', nsm, rho, V)
eid = 7
V = 2.0
self.check_solid(model, eid, 'CHEXA8', pid, 'PLSOLID', mid, 'MATHP', nsm, rho, V)
def test_cord2r_1(self):
grid = ['GRID 20143 7 -9.31-4 .11841 .028296']
coord = [
'CORD2R 7 1.135 .089237 -.0676 .135 .089237 -.0676',
' 1.135 .089237 .9324'
]
model = BDF(debug=False)
card = model.process_card(grid)
model.add_card(card, card[0])
card = model.process_card(coord)
model.add_card(card, card[0])
model.cross_reference()
g = model.Node(20143)
#print(g.Position(debug=False))
# by running it through Patran...
#GRID 20143 1.1067 .207647 -.068531
diff = g.Position() - array([1.106704, .207647, -0.068531])
msg = 'diff=%s' % diff
assert allclose(diff, 0.), msg
coord = model.Coord(7)
coord.T()
self.assertTrue(array_equal(coord.T(), coord.beta_n(2)))
def getNodes(self, grids, grids_expected, coords, debug=False):
model = BDF(debug=False)
for grid in grids:
(nid, cid, x, y, z) = grid
model.add_card(['GRID', nid, cid, x, y, z], 'GRID')
#gridObj = model.Node(nid)
#if debug:
#print(gridObj)
for coord in coords:
#print coord
(cid, rid, x, y, z) = coord
model.add_card(['CORD2R', cid, rid] + x + y + z, 'CORD2R')
coordObj = model.Coord(cid)
if debug:
print(coordObj)
model.cross_reference()
for (i, grid) in enumerate(grids_expected):
(nid, cid, x, y, z) = grid
nodes = model.grid
pos = nodes.get_positions([nid])
n = array([x, y, z])
msg = 'i=%s expected=%s actual=%s\n' % (i, n, pos)
print(msg)
assert allclose(n, pos), msg
def test_cord2_bad_01(self):
model = BDF(debug=False)
cards = [
['CORD2R', 1, 0, 0., 0., 0., 0., 0., 0., 0., 0.,
0.], # fails on self.k
['CORD2R', 2, 0, 0., 0., 0., 1., 0., 0., 0., 0.,
0.], # fails on normalize self.j
['CORD2R', 3, 0, 0., 0., 0., 1., 0., 0., 1., 1., 0.], # passes
['CORD2R', 4, 0, 0., 1., 0., 1., 0., 0., 1., 1., 0.], # passes
['CORD2R', 5, 4, 0., 1., 0., 1., 0., 0., 1., 1., 0.], # passes
]
for card in cards:
cid = card[1]
if cid in [1, 2]:
with self.assertRaises(RuntimeError):
model.add_card(card, card[0], is_list=True)
else:
model.add_card(card, card[0], is_list=True)
# this runs because it's got rid=0
cord4 = model.Coord(4)
cord4.transformToGlobal([0., 0., 0.])
# this doesn't run because rid != 0
cord5 = model.Coord(5)
with self.assertRaises(RuntimeError):
cord5.transformToGlobal([0., 0., 0.])
model.cross_reference()
def test_cord1_01(self):
model = BDF(debug=False)
card_count = {
'CCORD1R' : 1,
'GRID' : 3,
}
model.allocate(card_count)
cards = [
['CORD1R', 1, 1, 2, 3], # fails on self.k
['GRID', 1, 0, 0., 0., 0.],
['GRID', 2, 0, 1., 0., 0.],
['GRID', 3, 0, 1., 1., 0.],
]
for card in cards:
model.add_card(card, card[0], comment='comment', is_list=True)
c1 = model.Coord(1)
self.assertEquals(c1.G1(), 1)
self.assertEquals(c1.G2(), 2)
self.assertEquals(c1.G3(), 3)
model.cross_reference()
self.assertEquals(c1.G1(), 1)
self.assertEquals(c1.G2(), 2)
self.assertEquals(c1.G3(), 3)
self.assertEquals(c1.NodeIDs(), [1, 2, 3])
def getNodes(self, grids, grids_expected, coords):
model = BDF(debug=False)
card_count = {
'GRID' : len(grids),
'CORD2R' : len(coords),
}
model.allocate(card_count)
for grid in grids:
(nid, cid, x, y, z) = grid
model.add_card(['GRID', nid, cid, x, y, z], 'GRID')
for coord in coords:
(cid, rid, x, y, z) = coord
model.add_card(['CORD2R', cid, rid] + x + y + z, 'CORD2R')
#coordObj = model.Coord(cid)
model.cross_reference()
for (i, grid) in enumerate(grids_expected):
(nid, cid, x, y, z) = grid
nodes = model.grid
pos = nodes.get_position_by_node_id([nid])[0]
n = array([x, y, z])
msg = 'i=%s expected=%s actual=%s\n' % (i, n, pos)
print(msg)
assert allclose(n, pos), msg
def test_cord2r_02(self):
grid = ['GRID 20143 7 -9.31-4 .11841 .028296']
coord = ['CORD2R 7 1.135 .089237 -.0676 .135 .089237 -.0676',
' 1.135 .089237 .9324']
model = BDF(debug=False)
card_count = {
'GRID' : 1,
'CORD2R' : 1,
}
model.allocate(card_count)
card = model.process_card(grid)
model.add_card(card, card[0])
card = model.process_card(coord)
model.add_card(card, card[0])
model.cross_reference()
g = model.Node(20143)
#xyz = g.Position()
xyz = model.coords.get_global_position_by_node_id(20143, g.cp[0])[0]
# by running it through Patran...
#GRID 20143 1.1067 .207647 -.068531
expected = array([1.106704, .207647, -0.068531])
diff = xyz - expected
msg = '\nexpected=%s \nactual =%s \ndiff =%s' % (expected, xyz, diff)
assert allclose(diff, 0.), msg
coord = model.Coord(7)
coord.T()
self.assertTrue(array_equal(coord.T(), coord.beta_n(2)))
def add_card(self, card_lines, card_name, comment='', is_list=True):
card_name = card_name.upper()
self._increase_card_count(card_name)
if card_name in ['DEQATN']:
card_obj = card_lines
card = card_lines
else:
if is_list:
fields = card_lines
else:
fields = to_fields(card_lines, card_name)
# apply OPENMDAO syntax
if self._is_dynamic_syntax:
fields = [self._parse_dynamic_syntax(field) if '%' in
field[0:1] else field for field in fields]
card = wipe_empty_fields([interpret_value(field, fields)
if field is not None
else None for field in fields])
else: # leave everything as strings
card = wipe_empty_fields(fields)
card_obj = BDFCard(card)
if card_name == 'HYPER':
hyper = HYPER(card_obj, comment)
self.hyper[hyper.pid] = hyper
return
elif card_name == 'FLOW':
flow = FLOW(card_obj, comment)
self.flow[flow.flow_id] = flow
return
BDF.add_card(self, card, card_name, comment=comment, is_list=True)
def test_cord1c_01(self):
lines = ['cord1c,2,1,4,3']
grids = [
['GRID', 4, 0, 0.0, 0., 0.],
['GRID', 3, 0, 0.0, 0., 1.],
['GRID', 1, 0, 0.0, 1., 0.],
]
card_count = {
'CORD1C' : 1,
'GRID' : 3,
}
model = BDF(debug=False)
model.allocate(card_count)
model.add_card(lines, 'CORD1C', is_list=False)
for grid in grids:
model.add_card(grid, grid[0], is_list=True)
model.build()
size = 8
bdf_file = StringIO()
card = model.coords.slice_by_coord_id(2)
self.assertEquals(card.get_cid_by_coord_id(), 2)
self.assertEquals(card.get_rid_by_coord_id(), 0)
card.write_card(bdf_file, size=8, is_double=False)
def test_cord1c_01(self):
lines = ['cord1c,2,1,4,3']
grids = [
['GRID', 4, 0, 0.0, 0., 0.],
['GRID', 3, 0, 0.0, 0., 1.],
['GRID', 1, 0, 0.0, 1., 0.],
]
card_count = {
'CORD1C': 1,
'GRID': 3,
}
model = BDF(debug=False)
model.allocate(card_count)
model.add_card(lines, 'CORD1C', is_list=False)
for grid in grids:
model.add_card(grid, grid[0], is_list=True)
model.build()
size = 8
bdf_file = StringIO()
card = model.coords.slice_by_coord_id(2)
self.assertEquals(card.get_cid_by_coord_id(), 2)
self.assertEquals(card.get_rid_by_coord_id(), 0)
card.write_card(bdf_file, size=8, is_double=False)
def test_cord2_rcs_03(self):
"""
all points are located at <30,40,50>
"""
model = BDF(debug=False)
card_count = {
'GRID' : 3,
'CORD2R' : 1,
'CORD2C' : 1,
'CORD2S' : 2,
}
model.allocate(card_count)
cards = [
[
'CORD2S* 2 0 0. 0.',
'* 0. 0. 0. 1.* ',
'* 1. 0. 1.',],
[
#'$ Femap with NX Nastran Coordinate System 30 : rectangular in spherical',
'CORD2R* 30 2 14. 30.',
'* 70. 13.431863852 32.1458443949 75.2107442927* ',
'* 14.4583462334 33.4569982885 68.2297989286',],
[
#'$ Femap with NX Nastran Coordinate System 31 : cylindrical in spherical',
'CORD2C* 31 2 3. 42.',
'* -173. 2.86526881213 45.5425615252 159.180363517* ',
'* 3.65222385965 29.2536614627 -178.631312271',],
[
#'$ Femap with NX Nastran Coordinate System 32 : spherical in spherical',
'CORD2S* 32 2 22. 14.',
'* 85. 22.1243073983 11.9537753718 77.9978191005* ',
'* 21.0997242967 13.1806120497 88.4824763008',],
[
'GRID* 30 30 40.7437952957 -23.6254877994',
'* -33.09784854 0',],
[
'GRID* 31 31 62.9378078196 15.9774797923',
'* 31.0484428362 0',],
[
'GRID* 32 32 53.8270847449 95.8215692632',
'* 159.097767463 0',],
]
for lines in cards:
card = model.process_card(lines)
model.add_card(card, card[0])
model.build()
for nid in model.nodes:
a = array([30., 40., 50.])
b = model.Node(nid).get_position()
self.assertTrue(allclose(array([30., 40., 50.]), model.Node(nid).get_position()), str(a - b))
def test_cord2_rcs_02(self):
"""
all points are located at <30,40,50>
"""
model = BDF(debug=False)
card_count = {
'GRID' : 3,
'CORD2R' : 1,
'CORD2C' : 2,
'CORD2S' : 1,
}
model.allocate(card_count)
cards = [
[
'CORD2C* 1 0 0. 0.',
'* 0. 0. 0. 1.* ',
'* 1. 0. 1.',],
[
#'$ Femap with NX Nastran Coordinate System 20 : rectangular defined in cylindrical',
'CORD2R* 20 1 7. 20.',
'* -6. 7.07106781187 28.1301023542 -6.* ',
'* 7.70710678119 20. -5.29289321881',],
[
#'$ Femap with NX Nastran Coordinate System 21 : cylindrical defined in cylindrical',
'CORD2C* 21 1 15. -30.',
'* 12. 14.6565766735 -30.3177805524 12.9355733712* ',
'* 14.6234241583 -26.4257323272 11.9304419665',],
[
#'$ Femap with NX Nastran Coordinate System 22 : spherical defined in cylindrical',
'CORD2S* 22 1 5. -75.',
'* 20. 5.66032384035 -82.9319986389 19.8502545865* ',
'* 4.88876051026 -73.8006653677 19.0116094889',],
[
'GRID* 20 20 64.2559135157 -14.9400459772',
'* 27.3271005317 0',],
[
'GRID* 21 21 52.8328862418 -28.8729017195',
'* 34.615939507 0',],
[
'GRID* 22 22 61.1042111232 158.773483595',
'* -167.4951724 0',],
]
for lines in cards:
card = model.process_card(lines)
model.add_card(card, card[0])
model.build()
for nid in model.nodes:
a = array([30., 40., 50.])
b = model.Node(nid).get_position()
self.assertTrue(allclose(array([30., 40., 50.]), model.Node(nid).get_position()), str(a - b))
def test_cord2_rcs_03(self):
"""
all points are located at <30,40,50>
"""
model = BDF(debug=False)
cards = [
[
'CORD2S* 2 0 0. 0.',
'* 0. 0. 0. 1.* ',
'* 1. 0. 1.',
],
[ #'$ Femap with NX Nastran Coordinate System 30 : rectangular in spherical',
'CORD2R* 30 2 14. 30.',
'* 70. 13.431863852 32.1458443949 75.2107442927* ',
'* 14.4583462334 33.4569982885 68.2297989286',
],
[ #'$ Femap with NX Nastran Coordinate System 31 : cylindrical in spherical',
'CORD2C* 31 2 3. 42.',
'* -173. 2.86526881213 45.5425615252 159.180363517* ',
'* 3.65222385965 29.2536614627 -178.631312271',
],
[ #'$ Femap with NX Nastran Coordinate System 32 : spherical in spherical',
'CORD2S* 32 2 22. 14.',
'* 85. 22.1243073983 11.9537753718 77.9978191005* ',
'* 21.0997242967 13.1806120497 88.4824763008',
],
[
'GRID* 30 30 40.7437952957 -23.6254877994',
'* -33.09784854 0',
],
[
'GRID* 31 31 62.9378078196 15.9774797923',
'* 31.0484428362 0',
],
[
'GRID* 32 32 53.8270847449 95.8215692632',
'* 159.097767463 0',
],
]
for lines in cards:
card = model.process_card(lines)
model.add_card(card, card[0])
model.cross_reference()
for nid in model.nodes:
a = array([30., 40., 50.])
b = model.Node(nid).Position()
self.assertTrue(
allclose(array([30., 40., 50.]),
model.Node(nid).Position()), str(a - b))
def test_cord2_rcs_02(self):
"""
all points are located at <30,40,50>
"""
model = BDF(debug=False)
cards = [
[
'CORD2C* 1 0 0. 0.',
'* 0. 0. 0. 1.* ',
'* 1. 0. 1.',
],
[ #'$ Femap with NX Nastran Coordinate System 20 : rectangular defined in cylindrical',
'CORD2R* 20 1 7. 20.',
'* -6. 7.07106781187 28.1301023542 -6.* ',
'* 7.70710678119 20. -5.29289321881',
],
[ #'$ Femap with NX Nastran Coordinate System 21 : cylindrical defined in cylindrical',
'CORD2C* 21 1 15. -30.',
'* 12. 14.6565766735 -30.3177805524 12.9355733712* ',
'* 14.6234241583 -26.4257323272 11.9304419665',
],
[ #'$ Femap with NX Nastran Coordinate System 22 : spherical defined in cylindrical',
'CORD2S* 22 1 5. -75.',
'* 20. 5.66032384035 -82.9319986389 19.8502545865* ',
'* 4.88876051026 -73.8006653677 19.0116094889',
],
[
'GRID* 20 20 64.2559135157 -14.9400459772',
'* 27.3271005317 0',
],
[
'GRID* 21 21 52.8328862418 -28.8729017195',
'* 34.615939507 0',
],
[
'GRID* 22 22 61.1042111232 158.773483595',
'* -167.4951724 0',
],
]
for lines in cards:
card = model.process_card(lines)
model.add_card(card, card[0])
model.cross_reference()
for nid in model.nodes:
a = array([30., 40., 50.])
b = model.Node(nid).Position()
self.assertTrue(
allclose(array([30., 40., 50.]),
model.Node(nid).Position()), str(a - b))
def test_solid_03(self):
"""checks linear static solid material"""
mid = 2
pid = 4
rho = 0.1
cards = [
#$ Solid Nodes
['GRID', 11, 0, 0., 0., 0., 0],
['GRID', 12, 0, 1., 0., 0., 0],
['GRID', 13, 0, 1., 1., 0., 0],
['GRID', 14, 0, 0., 1., 0., 0],
['GRID', 15, 0, 0., 0., 2., 0],
['GRID', 16, 0, 1., 0., 2., 0],
['GRID', 17, 0, 1., 1., 2., 0],
['GRID', 18, 0, 0., 1., 2., 0],
# Solids
['CHEXA', 7, pid, 11, 12, 13, 14, 15, 16, 17, 18],
['CTETRA', 8, pid, 11, 12, 13, 15],
# Solid Nodes
['GRID', 21, 0, 0., 0., 0., 0,],
['GRID', 22, 0, 1., 0., 0., 0,],
['GRID', 23, 0, 1., 1., 0., 0,],
['GRID', 24, 0, 0., 0., 2., 0,],
['GRID', 25, 0, 1., 0., 2., 0,],
['GRID', 26, 0, 1., 1., 2., 0,],
['CPENTA', 9, pid, 21, 22, 23, 24, 25, 26],
# static
['PSOLID', pid, mid, 0],
['MAT1', mid, 1.0, 2.0, 3.0, rho],
['MATS1', mid, None, 'PLASTIC', 0.0, 1, 1, 100000., ],
]
card_count = {
'GRID' : 14,
'CTETRA4': 1,
'CPENTA6': 1,
'CHEXA8' : 1,
'PSOLID': 1,
'MAT1': 1,
'MATS1': 1,
}
model = BDF(debug=True)
model.allocate(card_count)
for fields in cards:
model.add_card(fields, fields[0], is_list=True)
model.cross_reference()
def test_crod_01(self):
model = BDF(debug=debug)
model.allocate({'CROD': 1})
lines = ['CROD 10 100 10 2']
model.add_card(lines, 'CROD', is_list=False)
size = 8
f = StringIO()
#card = CROD(card)
card = model.crod.slice_by_element_id([10])
card.write_card(f, size)
#card.rawFields()
self.assertEquals(card.get_element_id_by_element_index(), 10)
self.assertEquals(card.get_property_id_by_element_index(), 100)
def test_crod_01(self):
model = BDF(debug=debug)
model.allocate({"CROD": 1})
lines = ["CROD 10 100 10 2"]
model.add_card(lines, "CROD", is_list=False)
size = 8
f = StringIO()
# card = CROD(card)
card = model.crod.slice_by_element_id([10])
card.write_card(f, size)
# card.rawFields()
self.assertEquals(card.get_element_id_by_element_index(), 10)
self.assertEquals(card.get_property_id_by_element_index(), 100)
def test_cord2_rcs_01(self):
"""
all points are located at <30,40,50>
"""
model = BDF(debug=False)
card_count = {
'GRID' : 3,
'CORD2R' : 1,
'CORD2C' : 1,
'CORD2S' : 1,
}
model.allocate(card_count)
cards = [
[
#'$ Femap with NX Nastran Coordinate System 10 : rectangular defined in a rectangular',
'CORD2R* 10 0 10. 5.',
'* 3. 10.3420201433 4.53015368961 3.81379768136* ',
'* 10.7198463104 5.68767171433 3.09449287122',],
[
#'$ Femap with NX Nastran Coordinate System 11 : cylindrical defined in rectangular',
'CORD2C* 11 0 7. 3.',
'* 9. 7.64278760969 2.73799736977 9.71984631039* ',
'* 7.75440650673 3.37968226211 8.46454486422',],
[
#'$ Femap with NX Nastran Coordinate System 12 : spherical defined in rectangular',
'CORD2S* 12 0 12. 8.',
'* 5. 12.6427876097 7.86697777844 5.75440650673* ',
'* 12.6634139482 8.58906867688 4.53861076379',],
[
'GRID* 10 10 42.9066011565 34.2422137135',
'* 28.6442730262 0',],
[
'GRID* 11 11 48.8014631871 78.8394787869',
'* 34.6037164304 0',],
[
'GRID* 12 12 58.0775343829 44.7276544324',
'* 75.7955331161 0',],
]
for lines in cards:
card = model.process_card(lines)
model.add_card(card, card[0])
model.build()
for nid in model.nodes:
a = array([30., 40., 50.])
b = model.Node(nid).get_position()
self.assertTrue(allclose(array([30., 40., 50.]), model.Node(nid).get_position()), str(a - b))
def test_cord1r_02(self):
model = BDF(debug=False)
card_count = {
'CORD1R': 1,
'GRID': 3,
}
#model.allocate(card_count)
cards = [
['CORD1R', 1, 1, 2, 3], # fails on self.k
['GRID', 1, 0, 0., 0., 0.],
['GRID', 2, 0, 1., 0., 0.],
['GRID', 3, 0, 1., 1., 0.],
]
for card in cards:
model.add_card(card, card[0], comment='comment\n', is_list=True)
model.build()
c1 = model.coords.slice_by_coord_id(1)
def test_cord1r_02(self):
model = BDF(debug=False)
card_count = {
'CORD1R' : 1,
'GRID' : 3,
}
#model.allocate(card_count)
cards = [
['CORD1R', 1, 1, 2, 3], # fails on self.k
['GRID', 1, 0, 0., 0., 0.],
['GRID', 2, 0, 1., 0., 0.],
['GRID', 3, 0, 1., 1., 0.],
]
for card in cards:
model.add_card(card, card[0], comment='comment\n', is_list=True)
model.build()
c1 = model.coords.slice_by_coord_id(1)
def test_conrod_01(self):
model = BDF(debug=debug)
model.allocate({'CONROD': 1})
eid = 10
nid1 = 2
nid2 = 3
mid = 5
A = 27.0
lines = ['conrod,%i, %i, %i, %i, %f' % (eid, nid1, nid2, mid, A)]
model.add_card(lines, 'CONROD', is_list=False)
size = 8
card = model.conrod.slice_by_element_id([eid])
f = StringIO()
card.write_card(f, size)
#card.rawFields()
self.assertEquals(card.get_element_id_by_element_index(), eid)
self.assertEquals(card.get_material_id_by_element_index(), mid)
def test_conrod_01(self):
model = BDF(debug=debug)
model.allocate({"CONROD": 1})
eid = 10
nid1 = 2
nid2 = 3
mid = 5
A = 27.0
lines = ["conrod,%i, %i, %i, %i, %f" % (eid, nid1, nid2, mid, A)]
model.add_card(lines, "CONROD", is_list=False)
size = 8
card = model.conrod.slice_by_element_id([eid])
f = StringIO()
card.write_card(f, size)
# card.rawFields()
self.assertEquals(card.get_element_id_by_element_index(), eid)
self.assertEquals(card.get_material_id_by_element_index(), mid)
def test_cord2_rcs_01(self):
"""
all points are located at <30,40,50>
"""
model = BDF(debug=False)
cards = [
[ #'$ Femap with NX Nastran Coordinate System 10 : rectangular defined in a rectangular',
'CORD2R* 10 0 10. 5.',
'* 3. 10.3420201433 4.53015368961 3.81379768136* ',
'* 10.7198463104 5.68767171433 3.09449287122',
],
[ #'$ Femap with NX Nastran Coordinate System 11 : cylindrical defined in rectangular',
'CORD2C* 11 0 7. 3.',
'* 9. 7.64278760969 2.73799736977 9.71984631039* ',
'* 7.75440650673 3.37968226211 8.46454486422',
],
[ #'$ Femap with NX Nastran Coordinate System 12 : spherical defined in rectangular',
'CORD2S* 12 0 12. 8.',
'* 5. 12.6427876097 7.86697777844 5.75440650673* ',
'* 12.6634139482 8.58906867688 4.53861076379',
],
[
'GRID* 10 10 42.9066011565 34.2422137135',
'* 28.6442730262 0',
],
[
'GRID* 11 11 48.8014631871 78.8394787869',
'* 34.6037164304 0',
],
[
'GRID* 12 12 58.0775343829 44.7276544324',
'* 75.7955331161 0',
],
]
for lines in cards:
card = model.process_card(lines)
model.add_card(card, card[0])
model.cross_reference()
for nid in model.nodes:
a = array([30., 40., 50.])
b = model.Node(nid).Position()
self.assertTrue(
allclose(array([30., 40., 50.]),
model.Node(nid).Position()), str(a - b))
def test_grid_01(self):
model = BDF(debug=False)
cards = [
#['CORD1R', 1, 1, 2, 3], # fails on self.k
['GRID', 1, 0, 0., 0., 0.],
['GRID', 2, 0, 1., 0., 0.],
['GRID', 4, 0, 1., 2., 3.],
]
card_count = {
'GRID': 3,
}
model.allocate(card_count)
for card in cards:
model.add_card(card, card[0], comment='comment', is_list=True)
#+------+-----+----+----+----+----+----+----+------+
#| 0 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 |
#+======+=====+====+====+====+====+====+====+======+
#| GRID | NID | CP | X1 | X2 | X3 | CD | PS | SEID |
#+------+-----+----+----+----+----+----+----+------+
node = model.Node(4)
def test_grid_01(self):
model = BDF(debug=False)
card_lines = [
#['CORD1R', 1, 1, 2, 3], # fails on self.k
['GRID', 1, 0, 0., 0., 0.],
['GRID', 2, 0, 1., 0., 0.],
['GRID', 4, 0, 1., 2., 3.],
]
#card_count = {
#'GRID': 3,
#}
cards, card_count = model.add_cards_lines(card_lines)
model.allocate(card_count, cards)
for card in cards:
model.add_card(card, card[0], comment='comment', is_list=True)
model.build()
#+------+-----+----+----+----+----+----+----+------+
#| 0 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 |
#+======+=====+====+====+====+====+====+====+======+
#| GRID | NID | CP | X1 | X2 | X3 | CD | PS | SEID |
#+------+-----+----+----+----+----+----+----+------+
node = model.grid.slice_by_node_id(4)
def test_grid_01(self):
model = BDF(debug=False)
card_lines = [
#['CORD1R', 1, 1, 2, 3], # fails on self.k
['GRID', 1, 0, 0., 0., 0.],
['GRID', 2, 0, 1., 0., 0.],
['GRID', 4, 0, 1., 2., 3.],
]
#card_count = {
#'GRID': 3,
#}
cards, card_count = model.add_cards_lines(card_lines)
model.allocate(card_count, cards)
for card in cards:
model.add_card(card, card[0], comment='comment', is_list=True)
model.build()
#+------+-----+----+----+----+----+----+----+------+
#| 0 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 |
#+======+=====+====+====+====+====+====+====+======+
#| GRID | NID | CP | X1 | X2 | X3 | CD | PS | SEID |
#+------+-----+----+----+----+----+----+----+------+
node = model.grid.slice_by_node_id(4)
def test_grid_01(self):
model = BDF(debug=False)
cards = [
#['CORD1R', 1, 1, 2, 3], # fails on self.k
['GRID', 1, 0, 0., 0., 0.],
['GRID', 2, 0, 1., 0., 0.],
['GRID', 4, 0, 1., 2., 3.],
]
for card in cards:
model.add_card(card, card[0], comment='comment', is_list=True)
#+------+-----+----+----+----+----+----+----+------+
#| 0 | 1 | 2 | 3 | 4 | 5 | 6 | 7 | 8 |
#+======+=====+====+====+====+====+====+====+======+
#| GRID | NID | CP | X1 | X2 | X3 | CD | PS | SEID |
#+------+-----+----+----+----+----+----+----+------+
node = model.Node(4)
self.assertEqual(node.get_field(1), 4)
self.assertEqual(node.get_field(2), 0)
self.assertEqual(node.get_field(3), 1.)
self.assertEqual(node.get_field(4), 2.)
self.assertEqual(node.get_field(5), 3.)
node.update_field(1, 5)
node.update_field(2, 6)
node.update_field(3, 7.)
node.update_field(4, 8.)
node.update_field(5, 9.)
with self.assertRaises(KeyError):
node.update_field(9, 'dummy')
self.assertEqual(node.get_field(1), 5)
self.assertEqual(node.get_field(2), 6)
self.assertEqual(node.get_field(3), 7.)
self.assertEqual(node.get_field(4), 8.)
self.assertEqual(node.get_field(5), 9.)
with self.assertRaises(KeyError):
node.get_field(9)
def test_cord2_bad_01(self):
model = BDF(debug=False)
cards = [
['CORD2R', 1, 0, 0., 0., 0.,
0., 0., 0.,
0., 0., 0.], # fails on self.k
['CORD2R', 2, 0, 0., 0., 0.,
1., 0., 0.,
0., 0., 0.], # fails on normalize self.j
['CORD2R', 3, 0, 0., 0., 0.,
1., 0., 0.,
1., 1., 0.], # passes
['CORD2R', 4, 0, 0., 1., 0.,
1., 0., 0.,
1., 1., 0.], # passes
['CORD2R', 5, 4, 0., 1., 0.,
1., 0., 0.,
1., 1., 0.], # passes
]
for card in cards:
cid = card[1]
if cid in [1, 2]:
with self.assertRaises(RuntimeError):
model.add_card(card, card[0], is_list=True)
else:
model.add_card(card, card[0], is_list=True)
# this runs because it's got rid=0
cord4 = model.Coord(4)
cord4.transformToGlobal([0., 0., 0.])
# this doesn't run because rid != 0
cord5 = model.Coord(5)
with self.assertRaises(RuntimeError):
cord5.transformToGlobal([0., 0., 0.])
model.cross_reference()
def add_card(self, card_lines, card_name, comment='', is_list=True):
card_name = card_name.upper()
self._increase_card_count(card_name)
if card_name in ['DEQATN']:
card_obj = card_lines
card = card_lines
else:
if is_list:
fields = card_lines
else:
fields = to_fields(card_lines, card_name)
# apply OPENMDAO syntax
if self._is_dynamic_syntax:
fields = [
self._parse_dynamic_syntax(field)
if '%' in field[0:1] else field for field in fields
]
card = wipe_empty_fields([
interpret_value(field, fields)
if field is not None else None for field in fields
])
else: # leave everything as strings
card = wipe_empty_fields(fields)
card_obj = BDFCard(card)
if card_name == 'HYPER':
hyper = HYPER(card_obj, comment)
self.hyper[hyper.pid] = hyper
return
elif card_name == 'FLOW':
flow = FLOW(card_obj, comment)
self.flow[flow.flow_id] = flow
return
BDF.add_card(self, card, card_name, comment=comment, is_list=True)
def test_cord2c_01(self):
lines_a = [
'CORD2C* 3 0 0. 0.',
'* 0. 0. 0. 1.*',
'* 1. 0. 1.'
]
lines_b = [
'CORD2R 4 3 10. 0. 5. 10. 90. 5.',
' 10. 0. 6.'
]
card_count = {
'CORD2C': 1,
'CORD2R': 1,
}
model = BDF(debug=False)
cards = {
'CORD2C': ['', lines_a],
'CORD2R': ['', lines_b],
}
model.add_cards(cards, card_count)
#model.allocate(cards, card_count)
card = model.add_card(lines_a, 'CORD2C', is_list=False)
card = model.add_card(lines_b, 'CORD2R', is_list=False)
model.build()
cord2r = model.coords.slice_by_coord_id(3)
#print(type(cord2r))
self.assertEquals(cord2r.get_cid_by_coord_id(), 3)
self.assertEquals(cord2r.get_rid_by_coord_id(), 0)
cord2r = model.coords.slice_by_coord_id(4)
#print(type(cord2r))
self.assertEquals(cord2r.get_cid_by_coord_id(), 4)
self.assertEquals(cord2r.get_rid_by_coord_id(), 3)
i = model.coords.get_coord_index_by_coord_id(4)
T = model.coords.T[i, :, :].reshape(3, 3)
Ti = T[0, :]
Tj = T[1, :]
Tk = T[2, :]
msg = 'i=%s expected=(0., 0., 1.)' % Ti
self.assertTrue(allclose(Ti, array([0., 0., 1.])), msg)
delta = Tj - array([1., 1., 0.]) / 2**0.5
self.assertTrue(allclose(Tj, array([1., 1., 0.]) / 2**0.5), str(delta))
delta = Tk - array([-1., 1., 0.]) / 2**0.5
self.assertTrue(allclose(Tk,
array([-1., 1., 0.]) / 2**0.5), str(delta))
def test_cord2c_01(self):
lines_a = [
'CORD2C* 3 0 0. 0.',
'* 0. 0. 0. 1.*',
'* 1. 0. 1.'
]
lines_b = [
'CORD2R 4 3 10. 0. 5. 10. 90. 5.',
' 10. 0. 6.'
]
card_count = {
'CORD2C' : 1,
'CORD2R' : 1,
}
model = BDF(debug=False)
cards = {
'CORD2C' : ['', lines_a],
'CORD2R' : ['', lines_b],
}
model.add_cards(cards, card_count)
#model.allocate(cards, card_count)
card = model.add_card(lines_a, 'CORD2C', is_list=False)
card = model.add_card(lines_b, 'CORD2R', is_list=False)
model.build()
cord2r = model.coords.slice_by_coord_id(3)
#print(type(cord2r))
self.assertEquals(cord2r.get_cid_by_coord_id(), 3)
self.assertEquals(cord2r.get_rid_by_coord_id(), 0)
cord2r = model.coords.slice_by_coord_id(4)
#print(type(cord2r))
self.assertEquals(cord2r.get_cid_by_coord_id(), 4)
self.assertEquals(cord2r.get_rid_by_coord_id(), 3)
i = model.coords.get_coord_index_by_coord_id(4)
T = model.coords.T[i, :, :].reshape(3, 3)
Ti = T[0, :]
Tj = T[1, :]
Tk = T[2, :]
msg = 'i=%s expected=(0., 0., 1.)' % Ti
self.assertTrue(allclose(Ti, array([0., 0., 1.])), msg)
delta = Tj - array([1., 1., 0.]) / 2**0.5
self.assertTrue(allclose(Tj, array([1., 1., 0.]) / 2**0.5), str(delta))
delta = Tk - array([-1., 1., 0.]) / 2**0.5
self.assertTrue(allclose(Tk, array([-1., 1., 0.]) / 2**0.5), str(delta))
def test_solid_03(self):
"""checks linear static solid material"""
mid = 2
pid = 4
rho = 0.1
cards = [
#$ Solid Nodes
['GRID', 11, 0, 0., 0., 0., 0],
['GRID', 12, 0, 1., 0., 0., 0],
['GRID', 13, 0, 1., 1., 0., 0],
['GRID', 14, 0, 0., 1., 0., 0],
['GRID', 15, 0, 0., 0., 2., 0],
['GRID', 16, 0, 1., 0., 2., 0],
['GRID', 17, 0, 1., 1., 2., 0],
['GRID', 18, 0, 0., 1., 2., 0],
# Solids
['CHEXA', 7, pid, 11, 12, 13, 14, 15, 16, 17, 18],
['CTETRA', 8, pid, 11, 12, 13, 15],
# Solid Nodes
[
'GRID',
21,
0,
0.,
0.,
0.,
0,
],
[
'GRID',
22,
0,
1.,
0.,
0.,
0,
],
[
'GRID',
23,
0,
1.,
1.,
0.,
0,
],
[
'GRID',
24,
0,
0.,
0.,
2.,
0,
],
[
'GRID',
25,
0,
1.,
0.,
2.,
0,
],
[
'GRID',
26,
0,
1.,
1.,
2.,
0,
],
['CPENTA', 9, pid, 21, 22, 23, 24, 25, 26],
# static
['PSOLID', pid, mid, 0],
['MAT1', mid, 1.0, 2.0, 3.0, rho],
[
'MATS1',
mid,
None,
'PLASTIC',
0.0,
1,
1,
100000.,
],
]
card_count = {
'GRID': 14,
'CTETRA4': 1,
'CPENTA6': 1,
'CHEXA8': 1,
'PSOLID': 1,
'MAT1': 1,
'MATS1': 1,
}
model = BDF(debug=True)
model.allocate(card_count)
for fields in cards:
model.add_card(fields, fields[0], is_list=True)
model.cross_reference()
def test_pcomp_02(self):
"""
symmetrical, nsm=0.0 and nsm=1.0
"""
model = BDF()
pid = 1
z0 = 0.
nsm = 0.
sb = 0.
ft = 'HOFF'
TRef = 0.
ge = 0.
lam = 'SYM' # isSymmetrical SYM
Mid = [1, 2, 3]
Theta = [0., 10., 20.]
T = [.1, .2, .3]
Sout = ['YES', 'YES', 'NO'] # 0-NO, 1-YES
pcomp = ['PCOMP', pid, z0, nsm, sb, ft, TRef, ge, lam,
Mid[0], T[0], Theta[0], Sout[0],
Mid[1], T[1], Theta[1], Sout[1],
Mid[2], T[2], Theta[2], Sout[2]]
#----
mid = 1
E = 3.0e7
G = None
nu = None
rho = 1.0
a = None
St = None
Sc = None
Ss = None
Mcsid = None
mat1_a = ['MAT1', mid, E, G, nu, rho, a, TRef, ge, St, Sc, Ss, Mcsid]
mat1_b = ['MAT1', mid + 1, E, G, nu, rho, a, TRef, ge, St, Sc, Ss, Mcsid]
mat1_c = ['MAT1', mid + 2, E, G, nu, rho, a, TRef, ge, St, Sc, Ss, Mcsid]
card_count = {
'PCOMP' : 1,
'MAT1' : 3,
}
model.allocate(card_count)
model.add_card(pcomp, pcomp[0], is_list=True)
model.add_card(mat1_a, 'MAT1', is_list=True)
model.add_card(mat1_b, 'MAT1', is_list=True)
model.add_card(mat1_c, 'MAT1', is_list=True)
model.build()
p = model.properties.properties_shell.pcomp[0]
self.assertTrue(p.is_symmetrical_by_property_id())
self.assertTrue(p.is_symmetrical_by_property_index())
self.assertEqual(p.get_nplies_by_property_id(), 6)
self.assertAlmostEqual(p.get_thickness_by_property_id()[0], 1.2)
self.assertAlmostEqual(p.get_thickness_by_property_id_ply(pid, 0), 0.1)
self.assertAlmostEqual(p.get_thickness_by_property_id_ply(pid, 1), 0.2)
self.assertAlmostEqual(p.get_thickness_by_property_id_ply(pid, 2), 0.3)
self.assertAlmostEqual(p.get_thickness_by_property_id_ply(pid, 3), 0.1)
self.assertAlmostEqual(p.get_thickness_by_property_id_ply(pid, 4), 0.2)
self.assertAlmostEqual(p.get_thickness_by_property_id_ply(pid, 5), 0.3)
with self.assertRaises(IndexError):
p.get_thickness_by_property_id_ply(pid, 6)
self.assertAlmostEqual(p.get_theta_by_property_id_ply(pid, 0), 0.)
self.assertAlmostEqual(p.get_theta_by_property_id_ply(pid, 1), 10.)
self.assertAlmostEqual(p.get_theta_by_property_id_ply(pid, 2), 20.)
self.assertAlmostEqual(p.get_theta_by_property_id_ply(pid, 3), 0.)
self.assertAlmostEqual(p.get_theta_by_property_id_ply(pid, 4), 10.)
self.assertAlmostEqual(p.get_theta_by_property_id_ply(pid, 5), 20.)
with self.assertRaises(IndexError):
p.get_theta_by_property_id_ply(pid, 6)
self.assertEqual(p.get_material_id_by_property_id_ply(pid, 0), 1)
self.assertEqual(p.get_material_id_by_property_id_ply(pid, 1), 2)
self.assertEqual(p.get_material_id_by_property_id_ply(pid, 2), 3)
self.assertEqual(p.get_material_id_by_property_id_ply(pid, 3), 1)
self.assertEqual(p.get_material_id_by_property_id_ply(pid, 4), 2)
self.assertEqual(p.get_material_id_by_property_id_ply(pid, 5), 3)
with self.assertRaises(IndexError):
p.get_material_id_by_property_id_ply(pid, 6)
self.assertTrue(allclose(
p.get_material_ids_by_property_id(pid),
[1, 2, 3, 1, 2, 3])
)
self.assertEqual(p.get_sout_by_property_id_ply(pid, 0), 'YES')
self.assertEqual(p.get_sout_by_property_id_ply(pid, 1), 'YES')
self.assertEqual(p.get_sout_by_property_id_ply(pid, 2), 'NO')
self.assertEqual(p.get_sout_by_property_id_ply(pid, 3), 'YES')
self.assertEqual(p.get_sout_by_property_id_ply(pid, 4), 'YES')
self.assertEqual(p.get_sout_by_property_id_ply(pid, 5), 'NO')
with self.assertRaises(IndexError):
p.get_sout_by_property_id_ply(pid, 6)
#---------------
#Rho
self.assertAlmostEqual(p.get_density_by_property_id_ply(pid, 0), 1.0)
self.assertAlmostEqual(p.get_density_by_property_id_ply(pid, 1), 1.0)
self.assertAlmostEqual(p.get_density_by_property_id_ply(pid, 2), 1.0)
self.assertAlmostEqual(p.get_density_by_property_id_ply(pid, 3), 1.0)
self.assertAlmostEqual(p.get_density_by_property_id_ply(pid, 4), 1.0)
self.assertAlmostEqual(p.get_density_by_property_id_ply(pid, 5), 1.0)
with self.assertRaises(IndexError):
p.get_density_by_property_id_ply(pid, 6)
# MassPerArea
#self.assertAlmostEqual(p.get_mass_per_area_by_property_id(), 1.2)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id(pid), 1.2)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 0), 0.1)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 1), 0.2)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 2), 0.3)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 3), 0.1)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 4), 0.2)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 5), 0.3)
with self.assertRaises(IndexError):
p.get_mass_per_area_by_property_id_ply(pid, 6)
self.assertEqual(p.get_nonstructural_mass_by_property_id(pid), 0.0)
#----------------------
# change the nsm to 1.0
p.set_nonstructural_mass_by_property_id(pid, 1.0)
self.assertEqual(p.get_nonstructural_mass_by_property_id(pid), 1.0)
# MassPerArea
self.assertAlmostEqual(p.get_mass_per_area_by_property_id(pid), 2.2)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 0, method='nplies'), 0.1+1/6.)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 1, method='nplies'), 0.2+1/6.)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 2, method='nplies'), 0.3+1/6.)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 3, method='nplies'), 0.1+1/6.)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 4, method='nplies'), 0.2+1/6.)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 5, method='nplies'), 0.3+1/6.)
with self.assertRaises(IndexError):
p.get_mass_per_area_by_property_id_ply(pid, 6)
def test_pcomp_01(self):
"""
asymmetrical, nsm=0.0 and nsm=1.0
"""
#self.pid = data[0]
#self.z0 = data[1]
#self.nsm = data[2]
#self.sb = data[3]
#self.ft = data[4]
#self.TRef = data[5]
#self.ge = data[6]
#self.lam = data[7]
#Mid = data[8]
#T = data[9]
#Theta = data[10]
#Sout = data[11]
pid = 1
z0 = 0.
nsm = 0.
sb = 0.
ft = 'HILL'
TRef = 0.
ge = 0.
#lam = 'NO' # isSymmetrical YES/NO
lam = 'BEND' # isSymmetrical YES/NO
Mid = [1, 2, 3]
Theta = [0., 10., 20.]
T = [.1, .2, .3]
Sout = ['YES', 'YES', 'NO'] # 0-NO, 1-YES
card_lines = [
'PCOMP', pid, z0, nsm, sb, ft, TRef, ge, lam,
Mid[0], T[0], Theta[0], Sout[0],
Mid[1], T[1], Theta[1], Sout[1],
Mid[2], T[2], Theta[2], Sout[2],
]
model = BDF()
card_count = {
'PCOMP' : 1,
'MAT1' : 3,
}
model.allocate(card_count)
model.add_card(card_lines, 'PCOMP', comment='', is_list=True)
# material...
mid = 1
E = 1e7
G = None
nu = None
rho = 1.0
a = None
St = None
Sc = None
Ss = None
Mcsid = None
mat1_a = ['MAT1', mid, E, G, nu, rho, a, TRef, ge, St, Sc, Ss, Mcsid]
mat1_b = ['MAT1', mid + 1, E, G, nu, rho, a, TRef, ge, St, Sc, Ss, Mcsid]
mat1_c = ['MAT1', mid + 2, E, G, nu, rho, a, TRef, ge, St, Sc, Ss, Mcsid]
model.add_card(mat1_a, 'MAT1', comment='', is_list=True)
model.add_card(mat1_b, 'MAT1', comment='', is_list=True)
model.add_card(mat1_c, 'MAT1', comment='', is_list=True)
#card = BDFCard(mat1)
#m = MAT1(model)
#m.allocate(1)
#m.add(card)
#m.build()
model.build()
#card = BDFCard(data)
#p = PCOMP(model)
#p = model.properties.pcomp
#p.add(card)
#p.build()
p = model.properties_shell.pcomp
m = model.materials.mat1
#self.assertFalse(p.is_symmetrical())
self.assertEqual(p.get_nplies_by_property_id(), 3)
self.assertAlmostEqual(p.get_thickness_by_property_id(pid), 0.6)
self.assertAlmostEqual(p.get_thickness_by_property_id_ply(pid, 0), 0.1)
self.assertAlmostEqual(p.get_thickness_by_property_id_ply(pid, 1), 0.2)
self.assertAlmostEqual(p.get_thickness_by_property_id_ply(pid, 2), 0.3)
with self.assertRaises(IndexError):
p.get_thickness_by_property_id_ply(pid, 3)
self.assertAlmostEqual(p.get_theta_by_property_id_ply(pid, 0), 0.)
self.assertAlmostEqual(p.get_theta_by_property_id_ply(pid, 1), 10.)
self.assertAlmostEqual(p.get_theta_by_property_id_ply(pid, 2), 20.)
with self.assertRaises(IndexError):
p.get_theta_by_property_id_ply(pid, 3)
self.assertEqual(p.get_material_id_by_property_id_ply(pid, 0), 1)
self.assertEqual(p.get_material_id_by_property_id_ply(pid, 1), 2)
self.assertEqual(p.get_material_id_by_property_id_ply(pid, 2), 3)
with self.assertRaises(IndexError):
p.get_material_id_by_property_id_ply(pid, 3)
#print('get_material_id_by_property_id = ', p.get_material_id_by_property_id(pid))
self.assertEqual(p.get_material_ids_by_property_id(pid)[0], 1)
self.assertEqual(p.get_material_ids_by_property_id(pid)[0], 1)
self.assertEqual(p.get_material_ids_by_property_id(pid)[1], 2)
self.assertEqual(p.get_material_ids_by_property_id(pid)[2], 3)
self.assertEqual(p.get_sout_by_property_id_ply(pid, 0), 'YES')
self.assertEqual(p.get_sout_by_property_id_ply(pid, 1), 'YES')
self.assertEqual(p.get_sout_by_property_id_ply(pid, 2), 'NO')
with self.assertRaises(IndexError):
p.get_sout_by_property_id_ply(pid, 3)
#for iply in range(len(p.plies)):
#mid = p.plies[iply][0]
#p.plies[iply][0] = m # MAT1
##p.mids = [m, m, m]
from StringIO import StringIO
f = StringIO()
#print(m.write_card(f, size=8, material_id=None))
p.write_card(f)
m.write_card(f)
print(f.getvalue())
#Mid
self.assertAlmostEqual(p.get_material_id_by_property_id_ply(pid, 0), 1)
self.assertAlmostEqual(p.get_material_id_by_property_id_ply(pid, 1), 2)
self.assertAlmostEqual(p.get_material_id_by_property_id_ply(pid, 2), 3)
with self.assertRaises(IndexError):
p.get_density_by_property_id_ply(pid, 3)
#Rho
self.assertAlmostEqual(p.get_density_by_property_id_ply(pid, 0), 1.0)
self.assertAlmostEqual(p.get_density_by_property_id_ply(pid, 1), 1.0)
self.assertAlmostEqual(p.get_density_by_property_id_ply(pid, 2), 1.0)
with self.assertRaises(IndexError):
p.get_density_by_property_id_ply(pid, 3)
# MassPerArea
self.assertAlmostEqual(p.get_mass_per_area_by_property_id(), 0.6)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id(pid), 0.6)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id([pid]), 0.6)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 0), 0.1)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 1), 0.2)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 2), 0.3)
#with self.assertRaises(IndexError):
#p.MassPerArea(3)
#----------------------
# change the nsm to 1.0
p.set_nonstructural_mass_by_property_id(pid, 1.0)
self.assertEqual(p.get_nonstructural_mass_by_property_id(), 1.0)
# MassPerArea
self.assertAlmostEqual(p.get_mass_per_area_by_property_id(), 1.6)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 0, method='nplies'), 0.1+1/3.)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 1, method='nplies'), 0.2+1/3.)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 2, method='nplies'), 0.3+1/3.)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 0, method='t'), 0.1+1/6.)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 1, method='t'), 0.2+2/6.)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 2, method='t'), 0.3+3/6.)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 0, method='rho*t'), 0.1+1/6.)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 1, method='rho*t'), 0.2+2/6.)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 2, method='rho*t'), 0.3+3/6.)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 0, method='t'), 0.1+1/6.)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 1, method='t'), 0.2+2/6.)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 2, method='t'), 0.3+3/6.)
with self.assertRaises(IndexError):
p.get_mass_per_area_by_property_id_ply(pid, 3, method='nplies')
z0 = p.get_z0_by_property_id(pid)
z = p.get_z_locations_by_property_id(pid)
z_expected = array([0., T[0], T[0]+T[1], T[0]+T[1]+T[2]])
for za, ze in zip(z, z_expected):
self.assertAlmostEqual(za, ze)
#z0 =
p.z0[0] = 1.0
z_expected = 1.0 + z_expected
z = p.get_z_locations_by_property_id(pid)
for za, ze in zip(z, z_expected):
self.assertAlmostEqual(za, ze)
def get_mass(self, nid1, nid2, xyz1, xyz2, eid, pid, mid, A, J, c, nsm, E,
G, nu, rho, L):
"""tests a CROD and a CONROD"""
card_count = {
'CONROD': 1,
'CTUBE': 1,
'PTUBE': 1,
'CROD': 1,
'PROD': 1,
'GRID': 2,
'MAT1': 1,
}
model = BDF(debug=debug)
model.allocate(card_count)
lines = [
'conrod,%i, %i, %i, %i, %f, %f, %f, %f' %
(eid, nid1, nid2, mid, A, J, c, nsm)
]
model.add_card(lines, 'conrod', is_list=False)
lines = ['crod,%i, %i, %i, %i' % (eid + 1, pid, nid1, nid2)]
model.add_card(lines, 'crod', is_list=False)
#lines = ['ctube,%i, %i, %i, %i' % (eid+2, pid+1, nid1, nid2)]
#model.add_card(lines, 'ctube', is_list=False)
lines = ['prod,%i, %i, %f, %f, %f, %f' % (pid, mid, A, J, c, nsm)]
model.add_card(lines, 'prod', is_list=False)
OD1 = sqrt(4 * A / pi)
t = 0.
OD2 = OD1
#lines = ['ptube,%i, %i, %f, %f, %f, %f' % (pid+1, mid, OD1, t, nsm, OD2)]
#model.add_card(lines, 'ptube', is_list=False)
lines = ['mat1,%i, %.2e, %.2e, %f, %f' % (mid, E, G, nu, rho)]
model.add_card(lines, 'mat1', is_list=False)
lines = [
'grid,%i, %i, %f, %f, %f' % (nid1, 0, xyz1[0], xyz1[1], xyz1[2])
]
model.add_card(lines, 'grid', is_list=False)
lines = [
'grid,%i, %i, %f, %f, %f' % (nid2, 0, xyz2[0], xyz2[1], xyz2[2])
]
model.add_card(lines, 'grid', is_list=False)
model.cross_reference()
mass = L * (rho * A + nsm)
# conrod
conrod = model.conrod[eid]
self.assertEquals(conrod.get_element_id_by_element_index(), eid)
#self.assertEquals(conrod.get_property_id_by_element_index(), None)
self.assertEquals(conrod.get_material_id_by_element_index(), mid)
self.assertEquals(conrod.get_length_by_element_index(), L)
#self.assertEquals(conrod.Nsm(), nsm)
self.assertEquals(conrod.get_mass_by_element_index(), mass)
#self.assertEquals(conrod.E(), E)
#self.assertEquals(conrod.G(), G)
#self.assertEquals(conrod.Area(), A)
#self.assertEquals(conrod.J(), J)
#self.assertEquals(conrod.C(), c)
#self.assertEquals(conrod.Rho(), rho)
# crod
conrod = model.conrod[eid + 1]
self.assertEquals(crod.get_element_id_by_element_index(), eid + 1)
self.assertEquals(crod.get_property_id_by_element_index(), pid)
self.assertEquals(crod.get_material_id_by_element_index(), mid)
self.assertEquals(crod.get_length_by_element_index(), L)
#self.assertEquals(crod.Nsm(), nsm)
self.assertEquals(crod.get_mass_by_element_index(), mass)
#self.assertEquals(crod.E(), E)
#self.assertEquals(crod.G(), G)
#self.assertEquals(crod.Area(), A)
#self.assertEquals(crod.J(), J)
#self.assertEquals(crod.C(), c)
#self.assertEquals(crod.Rho(), rho)
#self.assertEquals(crod.Nu(), nu)
# prod
self.assertEquals(prod.Pid(), pid)
self.assertEquals(prod.Mid(), mid)
self.assertEquals(prod.Nsm(), nsm)
self.assertEquals(prod.E(), E)
self.assertEquals(prod.G(), G)
self.assertEquals(prod.Area(), A)
self.assertEquals(prod.J(), J)
self.assertEquals(prod.C(), c)
self.assertEquals(prod.Rho(), rho)
# ctube
if 0:
self.assertEquals(ctube.Eid(), eid + 2)
self.assertEquals(ctube.Pid(), pid + 1)
self.assertEquals(ctube.Mid(), mid)
self.assertEquals(ctube.Length(), L)
self.assertEquals(ctube.Nsm(), nsm)
self.assertAlmostEquals(ctube.Mass(), mass, 5)
self.assertEquals(ctube.E(), E)
self.assertEquals(ctube.G(), G)
self.assertAlmostEquals(ctube.Area(), A, 5)
ctube.J()
self.assertEquals(ctube.Rho(), rho)
# ptube
self.assertEquals(ptube.Pid(), pid + 1)
self.assertEquals(ptube.Mid(), mid)
self.assertEquals(ptube.Nsm(), nsm)
self.assertEquals(ptube.E(), E)
self.assertEquals(ptube.G(), G)
self.assertAlmostEquals(ptube.Area(), A, 5)
ptube.J()
self.assertEquals(ptube.Rho(), rho)
def test_solid_01(self):
"""checks linear static solid material"""
mid = 2
pid = 4
rho = 0.1
cards = [
#$ Solid Nodes
['GRID', 11, 0, 0., 0., 0., 0],
['GRID', 12, 0, 1., 0., 0., 0],
['GRID', 13, 0, 1., 1., 0., 0],
['GRID', 14, 0, 0., 1., 0., 0],
['GRID', 15, 0, 0., 0., 2., 0],
['GRID', 16, 0, 1., 0., 2., 0],
['GRID', 17, 0, 1., 1., 2., 0],
['GRID', 18, 0, 0., 1., 2., 0],
# Solids
['CHEXA', 7, pid, 11, 12, 13, 14, 15, 16, 17, 18],
['CTETRA', 8, pid, 11, 12, 13, 15],
# Solid Nodes
['GRID', 21, 0, 0., 0., 0., 0,],
['GRID', 22, 0, 1., 0., 0., 0,],
['GRID', 23, 0, 1., 1., 0., 0,],
['GRID', 24, 0, 0., 0., 2., 0,],
['GRID', 25, 0, 1., 0., 2., 0,],
['GRID', 26, 0, 1., 1., 2., 0,],
['CPENTA', 9, pid, 21, 22, 23, 24, 25, 26],
# static
['PSOLID', pid, mid, 0],
['MAT1', mid, 1.0, 2.0, 3.0, rho]
]
card_count = {
'GRID' : 14,
'CTETRA4': 1,
'CPENTA6': 1,
'CHEXA8' : 1,
'PSOLID': 1,
'MAT1': 1,
}
model = BDF(debug=True)
model.allocate(card_count)
for fields in cards:
model.add_card(fields, fields[0], is_list=True)
model.cross_reference()
# CTETRA
eid = 8
mid = 2
pid = 4
nsm = 0.
V = 1. / 3.
rho = 0.1
self.check_solid(model, eid, 'CTETRA4', pid, 'PSOLID', mid, 'MAT1', nsm, rho, V)
eid = 9
V = 1.0
self.check_solid(model, eid, 'CPENTA6', pid, 'PSOLID', mid, 'MAT1', nsm, rho, V)
eid = 7
V = 2.0
self.check_solid(model, eid, 'CHEXA8', pid, 'PSOLID', mid, 'MAT1', nsm, rho, V)
def test_solid_02(self):
mid = 2
pid = 4
rho = 0.1
cards = [
#$ Solid Nodes
['GRID', 11, 0, 0., 0., 0., 0],
['GRID', 12, 0, 1., 0., 0., 0],
['GRID', 13, 0, 1., 1., 0., 0],
['GRID', 14, 0, 0., 1., 0., 0],
['GRID', 15, 0, 0., 0., 2., 0],
['GRID', 16, 0, 1., 0., 2., 0],
['GRID', 17, 0, 1., 1., 2., 0],
['GRID', 18, 0, 0., 1., 2., 0],
# Solids
['CHEXA', 7, pid, 11, 12, 13, 14, 15, 16, 17, 18],
['CTETRA', 8, pid, 11, 12, 13, 15],
# Solid Nodes
[
'GRID',
21,
0,
0.,
0.,
0.,
0,
],
[
'GRID',
22,
0,
1.,
0.,
0.,
0,
],
[
'GRID',
23,
0,
1.,
1.,
0.,
0,
],
[
'GRID',
24,
0,
0.,
0.,
2.,
0,
],
[
'GRID',
25,
0,
1.,
0.,
2.,
0,
],
[
'GRID',
26,
0,
1.,
1.,
2.,
0,
],
['CPENTA', 9, pid, 21, 22, 23, 24, 25, 26],
# hyperelastic
['PLSOLID', pid, mid, 'GRID'],
['MATHP', mid, None, None, None, rho],
]
card_count = {
'GRID': 14,
'CTETRA4': 1,
'CPENTA6': 1,
'CHEXA8': 1,
'PLSOLID': 1,
'MATHP': 1,
}
model = BDF(debug=True)
model.allocate(card_count)
for fields in cards:
model.add_card(fields, fields[0], is_list=True)
model.cross_reference()
# CTETRA
eid = 8
nsm = 0.
V = 1. / 3.
self.check_solid(model, eid, 'CTETRA4', pid, 'PLSOLID', mid, 'MATHP',
nsm, rho, V)
eid = 9
V = 1.0
self.check_solid(model, eid, 'CPENTA6', pid, 'PLSOLID', mid, 'MATHP',
nsm, rho, V)
eid = 7
V = 2.0
self.check_solid(model, eid, 'CHEXA8', pid, 'PLSOLID', mid, 'MATHP',
nsm, rho, V)
def test_cord2r_bad_01(self):
model = BDF(debug=True)
card_count = {
'GRID' : 4,
'CORD2R' : 3,
}
model.allocate(card_count)
grids = [
['GRID', 1, 0],
['GRID', 20, 0],
['GRID', 30, 0],
['GRID', 11, 5],
]
for grid in grids:
model.add_card(grid, 'GRID', is_list=True)
cards = [
['CORD2R', 1, 0, 0., 0., 0.,
0., 0., 0.,
0., 0., 0.], # fails on self.k
['CORD2R', 2, 0, 0., 0., 0.,
1., 0., 0.,
0., 0., 0.], # fails on normalize self.j
['CORD2R', 3, 0, 0., 0., 0.,
1., 0., 0.,
1., 1., 0.], # passes
['CORD2R', 4, 0, 0., 1., 0.,
1., 0., 0.,
1., 1., 0.], # passes
['CORD2R', 5, 4, 0., 1., 0.,
1., 0., 0.,
1., 1., 0.], # passes
]
for card in cards:
cid = card[1]
if cid in [1, 2]:
with self.assertRaises(RuntimeError):
model.add_card(card, card[0], is_list=True)
else:
model.add_card(card, card[0], is_list=True)
model.build()
# this runs because it's got rid=0
coords = model.coords #.slice_by_coord_id(4)
coords.transform_node_id_to_global_xyz(30) # nid
coords.transform_node_id_to_global_xyz([30, 1, 11]) # [nid, nid, nid]
coords.transform_node_id_to_local_by_coord_id([30, 1, 11], 4) # [nid, ...], cp_goal
coords.transform_node_id_to_local_by_coord_id([30, 1, 11], 0) # [nid, ...], cp_goal
xyz = [0., 0., 0.]
coords.transform_xyz_to_global_by_coord_id(xyz, 4) # [xyz], cp_initial
xyz = [
[0., 0., 0.],
[1., 1., 1.],
]
coords.transform_xyz_to_global_by_coord_id(xyz, 4) # [xyz], cp_initial
# global from global
coords.transform_xyz_to_global_by_coord_id(xyz, 0) # [xyz], cp_initial
# this doesn't run because rid != 0
with self.assertRaises(RuntimeError):
# cp=0 doesn't exist
coords.transform_xyz_to_global_by_coord_id(xyz, 2)
def get_mass(self, nid1, nid2, xyz1, xyz2, eid, pid, mid, A, J, c, nsm, E,
G, nu, rho, L):
"""tests a CROD and a CONROD"""
card_count = {
'CONROD': 1,
'CTUBE': 1,
'PTUBE': 1,
'CROD': 1,
'PROD': 1,
'GRID': 2,
'MAT1': 1,
}
model = BDF(debug=debug)
model.allocate(card_count)
lines = [
'conrod,%i, %i, %i, %i, %f, %f, %f, %f' %
(eid, nid1, nid2, mid, A, J, c, nsm)
]
model.add_card(lines, 'conrod', is_list=False)
lines = ['crod,%i, %i, %i, %i' % (eid + 1, pid, nid1, nid2)]
model.add_card(lines, 'crod', is_list=False)
lines = ['ctube,%i, %i, %i, %i' % (eid + 2, pid + 1, nid1, nid2)]
model.add_card(lines, 'ctube', is_list=False)
lines = ['prod,%i, %i, %f, %f, %f, %f' % (pid, mid, A, J, c, nsm)]
model.add_card(lines, 'prod', is_list=False)
OD1 = sqrt(4 * A / pi)
t = 0.
OD2 = OD1
lines = [
'ptube,%i, %i, %f, %f, %f, %f' % (pid + 1, mid, OD1, t, nsm, OD2)
]
model.add_card(lines, 'ptube', is_list=False)
lines = ['mat1,%i, %.2e, %.2e, %f, %f' % (mid, E, G, nu, rho)]
model.add_card(lines, 'mat1', is_list=False)
lines = [
'grid,%i, %i, %f, %f, %f' % (nid1, 0, xyz1[0], xyz1[1], xyz1[2])
]
model.add_card(lines, 'grid', is_list=False)
lines = [
'grid,%i, %i, %f, %f, %f' % (nid2, 0, xyz2[0], xyz2[1], xyz2[2])
]
model.add_card(lines, 'grid', is_list=False)
model.build()
mass = L * (rho * A + nsm)
f = StringIO()
model.write_bdf(out_filename=f,
interspersed=True,
size=8,
precision='single',
enddata=None)
print(f.getvalue())
#positions = model.get_positions()
grid_cid0 = None
# conrod
conrod = model.conrod.slice_by_element_id(eid)
self.assertEquals(conrod.get_element_id_by_element_index(), eid)
#self.assertEquals(conrod.get_property_id_by_element_id(), None)
self.assertEquals(conrod.get_material_id_by_element_id(eid), mid)
self.assertEquals(
conrod.get_length_by_element_index(i=None, grid_cid0=grid_cid0), L)
#self.assertEquals(conrod.Nsm(), nsm)
rhoi = conrod.get_density_by_element_id(eid)
Ai = conrod.get_area_by_element_id(eid)
Li = conrod.get_length_by_element_id(eid, grid_cid0=grid_cid0)
nsmi = conrod.get_non_structural_mass_by_element_id(eid)
massa = conrod.get_mass_by_element_index()
mass_msg_conrod = 'mass = L * (rho * A + nsm)\n'
mass_msg_conrod += 'L=%s expected=%s\n' % (Li, L)
mass_msg_conrod += 'rho=%s expected=%s\n' % (rhoi, rho)
mass_msg_conrod += 'A=%s expected=%s\n' % (Ai, A)
mass_msg_conrod += 'nsm=%s expected=%s\n' % (nsmi, nsm)
mass_msg_conrod += 'mass=%s actual=%s expected=%s\n' % (
Li * (rhoi * Ai + nsmi), massa, mass)
#mass_msg_conrod += 'mass=%s expected=%s\n' % (Li * (rhoi*Ai + nsmi), mass)
self.assertEquals(massa, mass, mass_msg_conrod)
#self.assertEquals(conrod.E(), E)
#self.assertEquals(conrod.G(), G)
#self.assertEquals(conrod.area(), A)
#self.assertEquals(conrod.J(), J)
#self.assertEquals(conrod.C(), c)
#self.assertEquals(conrod.Rho(), rho)
# crod
crod_eid = eid + 1
crod = model.crod.slice_by_element_id([crod_eid])
self.assertEquals(crod.get_element_id_by_element_index(), crod_eid)
self.assertEquals(crod.get_property_id_by_element_id(crod_eid), pid)
self.assertEquals(crod.get_material_id_by_element_id(crod_eid), mid)
rhoi = crod.get_density_by_element_id(crod_eid)
Ai = crod.get_area_by_element_id(crod_eid)
Li = crod.get_length_by_element_id(crod_eid, grid_cid0=grid_cid0)
nsmi = crod.get_non_structural_mass_by_element_id(crod_eid)
self.assertEquals(Li, L)
#self.assertEquals(crod.Nsm(), nsm)
massa = crod.get_mass_by_element_id(crod_eid)
mass_msg_crod = 'mass = L * (rho * A + nsm)\n'
mass_msg_crod += 'L=%s expected=%s\n' % (Li, L)
mass_msg_crod += 'rho=%s expected=%s\n' % (rhoi, rho)
mass_msg_crod += 'A=%s expected=%s\n' % (Ai, A)
mass_msg_crod += 'nsm=%s expected=%s\n' % (nsmi, nsm)
mass_msg_crod += 'mass=%s actual=%s expected=%s\n' % (
Li * (rhoi * Ai + nsmi), massa, mass)
self.assertEquals(massa, mass, mass_msg_crod)
#self.assertEquals(crod.E(), E)
#self.assertEquals(crod.G(), G)
#self.assertEquals(crod.area(), A)
#self.assertEquals(crod.J(), J)
#self.assertEquals(crod.C(), c)
#self.assertEquals(crod.Rho(), rho)
#self.assertEquals(crod.Nu(), nu)
# prod
prod = model.prod.slice_by_property_id([pid])
self.assertEquals(prod.property_id[0], pid)
self.assertEquals(prod.get_material_id_by_property_id(pid), mid)
self.assertEquals(prod.get_non_structural_mass_by_property_id(pid),
nsm)
self.assertEquals(prod.get_E_by_property_id(pid), E)
self.assertEquals(prod.get_G_by_property_id(pid), G)
self.assertEquals(prod.get_area_by_property_id(pid), A)
self.assertEquals(prod.get_J_by_property_id(pid), J)
self.assertEquals(prod.get_c_by_property_id(pid), c)
self.assertEquals(prod.get_density_by_property_id(pid), rho)
# ctube
if 1:
ctube_eid = eid + 2
ptube_pid = pid + 1
assert ctube_eid == 12, ctube_eid
assert ptube_pid == 68, ptube_pid
ctube = model.ctube.slice_by_element_id(ctube_eid)
self.assertEquals(ctube.get_element_id_by_element_index(),
ctube_eid)
self.assertEquals(ctube.get_property_id_by_element_id(ctube_eid),
ptube_pid)
self.assertEquals(ctube.get_material_id_by_element_id(ctube_eid),
mid)
self.assertEquals(
ctube.get_length_by_element_id(ctube_eid, grid_cid0), L)
self.assertEquals(
ctube.get_non_structural_mass_by_element_id(ctube_eid), nsm)
self.assertAlmostEquals(ctube.get_mass_by_element_id(ctube_eid),
mass, 5)
self.assertEquals(ctube.get_E_by_element_id(ctube_eid), E)
self.assertEquals(ctube.get_G_by_element_id(ctube_eid), G)
self.assertAlmostEquals(ctube.get_area_by_element_id(ctube_eid), A,
5)
ctube.get_J_by_element_id(ctube_eid)
self.assertEquals(ctube.get_density_by_element_id(), rho)
# ptube
ptube = model.ptube.slice_by_property_id(pid + 1)
self.assertEquals(ptube.get_property_id_by_property_index(),
pid + 1)
self.assertEquals(ptube.get_material_id_by_property_id(), mid)
self.assertEquals(ptube.get_non_structural_mass_by_property_id(),
nsm)
self.assertEquals(ptube.get_E_by_property_id(), E)
self.assertEquals(ptube.get_G_by_property_id(), G)
self.assertAlmostEquals(ptube.get_area_by_property_id(), A, 5)
ptube.get_J_by_property_id()
self.assertEquals(ptube.get_density_by_property_id(), rho)
def test_solid_01(self):
"""checks linear static solid material"""
mid = 2
pid = 4
rho = 0.1
cards = [
#$ Solid Nodes
['GRID', 11, 0, 0., 0., 0., 0],
['GRID', 12, 0, 1., 0., 0., 0],
['GRID', 13, 0, 1., 1., 0., 0],
['GRID', 14, 0, 0., 1., 0., 0],
['GRID', 15, 0, 0., 0., 2., 0],
['GRID', 16, 0, 1., 0., 2., 0],
['GRID', 17, 0, 1., 1., 2., 0],
['GRID', 18, 0, 0., 1., 2., 0],
# Solids
['CHEXA', 7, pid, 11, 12, 13, 14, 15, 16, 17, 18],
['CTETRA', 8, pid, 11, 12, 13, 15],
# Solid Nodes
[
'GRID',
21,
0,
0.,
0.,
0.,
0,
],
[
'GRID',
22,
0,
1.,
0.,
0.,
0,
],
[
'GRID',
23,
0,
1.,
1.,
0.,
0,
],
[
'GRID',
24,
0,
0.,
0.,
2.,
0,
],
[
'GRID',
25,
0,
1.,
0.,
2.,
0,
],
[
'GRID',
26,
0,
1.,
1.,
2.,
0,
],
['CPENTA', 9, pid, 21, 22, 23, 24, 25, 26],
# static
['PSOLID', pid, mid, 0],
['MAT1', mid, 1.0, 2.0, 3.0, rho]
]
card_count = {
'GRID': 14,
'CTETRA4': 1,
'CPENTA6': 1,
'CHEXA8': 1,
'PSOLID': 1,
'MAT1': 1,
}
model = BDF(debug=True)
model.allocate(card_count)
for fields in cards:
model.add_card(fields, fields[0], is_list=True)
model.cross_reference()
# CTETRA
eid = 8
mid = 2
pid = 4
nsm = 0.
V = 1. / 3.
rho = 0.1
self.check_solid(model, eid, 'CTETRA4', pid, 'PSOLID', mid, 'MAT1',
nsm, rho, V)
eid = 9
V = 1.0
self.check_solid(model, eid, 'CPENTA6', pid, 'PSOLID', mid, 'MAT1',
nsm, rho, V)
eid = 7
V = 2.0
self.check_solid(model, eid, 'CHEXA8', pid, 'PSOLID', mid, 'MAT1', nsm,
rho, V)
def test_solid_03(self):
"""checks linear static solid material"""
mid = 2
pid = 4
rho = 0.1
tableID = 42
cards = [
#$ Solid Nodes
['GRID', 11, 0, 0., 0., 0., 0],
['GRID', 12, 0, 1., 0., 0., 0],
['GRID', 13, 0, 1., 1., 0., 0],
['GRID', 14, 0, 0., 1., 0., 0],
['GRID', 15, 0, 0., 0., 2., 0],
['GRID', 16, 0, 1., 0., 2., 0],
['GRID', 17, 0, 1., 1., 2., 0],
['GRID', 18, 0, 0., 1., 2., 0],
# Solids
['CHEXA', 7, pid, 11, 12, 13, 14, 15, 16, 17, 18],
['CTETRA', 8, pid, 11, 12, 13, 15],
# Solid Nodes
[
'GRID',
21,
0,
0.,
0.,
0.,
0,
],
[
'GRID',
22,
0,
1.,
0.,
0.,
0,
],
[
'GRID',
23,
0,
1.,
1.,
0.,
0,
],
[
'GRID',
24,
0,
0.,
0.,
2.,
0,
],
[
'GRID',
25,
0,
1.,
0.,
2.,
0,
],
[
'GRID',
26,
0,
1.,
1.,
2.,
0,
],
['CPENTA', 9, pid, 21, 22, 23, 24, 25, 26],
# static
['PSOLID', pid, mid, 0],
['MAT1', mid, 1.0, 2.0, 3.0, rho],
[
'MATS1',
mid,
tableID,
'PLASTIC',
0.0,
1,
1,
100000.,
],
#['TABLEST'],
[
'TABLES1', tableID, 1, None, None, None, None, None, None, 1.0,
10.0, 2.0, 10.0, 'ENDT'
],
]
card_count = {
'GRID': 14,
'CPENTA6': 1,
'CTETRA4': 1,
'PSOLID': 1,
'CHEXA8': 1,
'MAT1': 1,
'MATS1': 1,
'TABLES1': 1,
}
model = BDF(debug=False)
model.allocate(card_count)
for fields in cards:
model.add_card(fields, fields[0], is_list=True)
model.cross_reference()
mat = model.materials[mid]
print('----MAT----', type(mat))
print(mat)
print('E = %s' % mat.get_E_by_material_index())
print('E = %s' % mat.get_E_by_material_id())
def get_mass(self, nid1, nid2, xyz1, xyz2, eid, pid, mid, A, J, c, nsm, E, G, nu, rho, L):
"""tests a CROD and a CONROD"""
card_count = {"CONROD": 1, "CTUBE": 1, "PTUBE": 1, "CROD": 1, "PROD": 1, "GRID": 2, "MAT1": 1}
model = BDF(debug=debug)
model.allocate(card_count)
lines = ["conrod,%i, %i, %i, %i, %f, %f, %f, %f" % (eid, nid1, nid2, mid, A, J, c, nsm)]
model.add_card(lines, "conrod", is_list=False)
lines = ["crod,%i, %i, %i, %i" % (eid + 1, pid, nid1, nid2)]
model.add_card(lines, "crod", is_list=False)
lines = ["ctube,%i, %i, %i, %i" % (eid + 2, pid + 1, nid1, nid2)]
model.add_card(lines, "ctube", is_list=False)
lines = ["prod,%i, %i, %f, %f, %f, %f" % (pid, mid, A, J, c, nsm)]
model.add_card(lines, "prod", is_list=False)
OD1 = sqrt(4 * A / pi)
t = 0.0
OD2 = OD1
lines = ["ptube,%i, %i, %f, %f, %f, %f" % (pid + 1, mid, OD1, t, nsm, OD2)]
model.add_card(lines, "ptube", is_list=False)
lines = ["mat1,%i, %.2e, %.2e, %f, %f" % (mid, E, G, nu, rho)]
model.add_card(lines, "mat1", is_list=False)
lines = ["grid,%i, %i, %f, %f, %f" % (nid1, 0, xyz1[0], xyz1[1], xyz1[2])]
model.add_card(lines, "grid", is_list=False)
lines = ["grid,%i, %i, %f, %f, %f" % (nid2, 0, xyz2[0], xyz2[1], xyz2[2])]
model.add_card(lines, "grid", is_list=False)
model.build()
mass = L * (rho * A + nsm)
f = StringIO()
model.write_bdf(out_filename=f, interspersed=True, size=8, precision="single", enddata=None)
print(f.getvalue())
# positions = model.get_positions()
grid_cid0 = None
# conrod
conrod = model.conrod.slice_by_element_id(eid)
self.assertEquals(conrod.get_element_id_by_element_index(), eid)
# self.assertEquals(conrod.get_property_id_by_element_id(), None)
self.assertEquals(conrod.get_material_id_by_element_id(eid), mid)
self.assertEquals(conrod.get_length_by_element_index(i=None, grid_cid0=grid_cid0), L)
# self.assertEquals(conrod.Nsm(), nsm)
rhoi = conrod.get_density_by_element_id(eid)
Ai = conrod.get_area_by_element_id(eid)
Li = conrod.get_length_by_element_id(eid, grid_cid0=grid_cid0)
nsmi = conrod.get_non_structural_mass_by_element_id(eid)
massa = conrod.get_mass_by_element_index()
mass_msg_conrod = "mass = L * (rho * A + nsm)\n"
mass_msg_conrod += "L=%s expected=%s\n" % (Li, L)
mass_msg_conrod += "rho=%s expected=%s\n" % (rhoi, rho)
mass_msg_conrod += "A=%s expected=%s\n" % (Ai, A)
mass_msg_conrod += "nsm=%s expected=%s\n" % (nsmi, nsm)
mass_msg_conrod += "mass=%s actual=%s expected=%s\n" % (Li * (rhoi * Ai + nsmi), massa, mass)
# mass_msg_conrod += 'mass=%s expected=%s\n' % (Li * (rhoi*Ai + nsmi), mass)
self.assertEquals(massa, mass, mass_msg_conrod)
# self.assertEquals(conrod.E(), E)
# self.assertEquals(conrod.G(), G)
# self.assertEquals(conrod.area(), A)
# self.assertEquals(conrod.J(), J)
# self.assertEquals(conrod.C(), c)
# self.assertEquals(conrod.Rho(), rho)
# crod
crod_eid = eid + 1
crod = model.crod.slice_by_element_id([crod_eid])
self.assertEquals(crod.get_element_id_by_element_index(), crod_eid)
self.assertEquals(crod.get_property_id_by_element_id(crod_eid), pid)
self.assertEquals(crod.get_material_id_by_element_id(crod_eid), mid)
rhoi = crod.get_density_by_element_id(crod_eid)
Ai = crod.get_area_by_element_id(crod_eid)
Li = crod.get_length_by_element_id(crod_eid, grid_cid0=grid_cid0)
nsmi = crod.get_non_structural_mass_by_element_id(crod_eid)
self.assertEquals(Li, L)
# self.assertEquals(crod.Nsm(), nsm)
massa = crod.get_mass_by_element_id(crod_eid)
mass_msg_crod = "mass = L * (rho * A + nsm)\n"
mass_msg_crod += "L=%s expected=%s\n" % (Li, L)
mass_msg_crod += "rho=%s expected=%s\n" % (rhoi, rho)
mass_msg_crod += "A=%s expected=%s\n" % (Ai, A)
mass_msg_crod += "nsm=%s expected=%s\n" % (nsmi, nsm)
mass_msg_crod += "mass=%s actual=%s expected=%s\n" % (Li * (rhoi * Ai + nsmi), massa, mass)
self.assertEquals(massa, mass, mass_msg_crod)
# self.assertEquals(crod.E(), E)
# self.assertEquals(crod.G(), G)
# self.assertEquals(crod.area(), A)
# self.assertEquals(crod.J(), J)
# self.assertEquals(crod.C(), c)
# self.assertEquals(crod.Rho(), rho)
# self.assertEquals(crod.Nu(), nu)
# prod
prod = model.prod.slice_by_property_id([pid])
self.assertEquals(prod.property_id[0], pid)
self.assertEquals(prod.get_material_id_by_property_id(pid), mid)
self.assertEquals(prod.get_non_structural_mass_by_property_id(pid), nsm)
self.assertEquals(prod.get_E_by_property_id(pid), E)
self.assertEquals(prod.get_G_by_property_id(pid), G)
self.assertEquals(prod.get_area_by_property_id(pid), A)
self.assertEquals(prod.get_J_by_property_id(pid), J)
self.assertEquals(prod.get_c_by_property_id(pid), c)
self.assertEquals(prod.get_density_by_property_id(pid), rho)
# ctube
if 1:
ctube_eid = eid + 2
ptube_pid = pid + 1
assert ctube_eid == 12, ctube_eid
assert ptube_pid == 68, ptube_pid
ctube = model.ctube.slice_by_element_id(ctube_eid)
self.assertEquals(ctube.get_element_id_by_element_index(), ctube_eid)
self.assertEquals(ctube.get_property_id_by_element_id(ctube_eid), ptube_pid)
self.assertEquals(ctube.get_material_id_by_element_id(ctube_eid), mid)
self.assertEquals(ctube.get_length_by_element_id(ctube_eid, grid_cid0), L)
self.assertEquals(ctube.get_non_structural_mass_by_element_id(ctube_eid), nsm)
self.assertAlmostEquals(ctube.get_mass_by_element_id(ctube_eid), mass, 5)
self.assertEquals(ctube.get_E_by_element_id(ctube_eid), E)
self.assertEquals(ctube.get_G_by_element_id(ctube_eid), G)
self.assertAlmostEquals(ctube.get_area_by_element_id(ctube_eid), A, 5)
ctube.get_J_by_element_id(ctube_eid)
self.assertEquals(ctube.get_density_by_element_id(), rho)
# ptube
ptube = model.ptube.slice_by_property_id(pid + 1)
self.assertEquals(ptube.get_property_id_by_property_index(), pid + 1)
self.assertEquals(ptube.get_material_id_by_property_id(), mid)
self.assertEquals(ptube.get_non_structural_mass_by_property_id(), nsm)
self.assertEquals(ptube.get_E_by_property_id(), E)
self.assertEquals(ptube.get_G_by_property_id(), G)
self.assertAlmostEquals(ptube.get_area_by_property_id(), A, 5)
ptube.get_J_by_property_id()
self.assertEquals(ptube.get_density_by_property_id(), rho)
def test_pcomp_01(self):
"""
asymmetrical, nsm=0.0 and nsm=1.0
"""
#self.pid = data[0]
#self.z0 = data[1]
#self.nsm = data[2]
#self.sb = data[3]
#self.ft = data[4]
#self.TRef = data[5]
#self.ge = data[6]
#self.lam = data[7]
#Mid = data[8]
#T = data[9]
#Theta = data[10]
#Sout = data[11]
pid = 1
z0 = 0.
nsm = 0.
sb = 0.
ft = 'HILL'
TRef = 0.
ge = 0.
#lam = 'NO' # isSymmetrical YES/NO
lam = 'BEND' # isSymmetrical YES/NO
Mid = [1, 2, 3]
Theta = [0., 10., 20.]
T = [.1, .2, .3]
Sout = ['YES', 'YES', 'NO'] # 0-NO, 1-YES
card_lines = [
'PCOMP',
pid,
z0,
nsm,
sb,
ft,
TRef,
ge,
lam,
Mid[0],
T[0],
Theta[0],
Sout[0],
Mid[1],
T[1],
Theta[1],
Sout[1],
Mid[2],
T[2],
Theta[2],
Sout[2],
]
model = BDF(debug=False)
card_count = {
'PCOMP': 1,
'MAT1': 3,
}
model.allocate(card_count)
model.add_card(card_lines, 'PCOMP', comment='', is_list=True)
# material...
mid = 1
E = 1e7
G = None
nu = None
rho = 1.0
a = None
St = None
Sc = None
Ss = None
Mcsid = None
mat1_a = ['MAT1', mid, E, G, nu, rho, a, TRef, ge, St, Sc, Ss, Mcsid]
mat1_b = [
'MAT1', mid + 1, E, G, nu, rho, a, TRef, ge, St, Sc, Ss, Mcsid
]
mat1_c = [
'MAT1', mid + 2, E, G, nu, rho, a, TRef, ge, St, Sc, Ss, Mcsid
]
model.add_card(mat1_a, 'MAT1', comment='', is_list=True)
model.add_card(mat1_b, 'MAT1', comment='', is_list=True)
model.add_card(mat1_c, 'MAT1', comment='', is_list=True)
#card = BDFCard(mat1)
#m = MAT1(model)
#m.allocate(1)
#m.add(card)
#m.build()
model.build()
#card = BDFCard(data)
#p = PCOMP(model)
#p = model.properties.pcomp
#p.add(card)
#p.build()
p = model.properties_shell.pcomp
m = model.materials.mat1
#self.assertFalse(p.is_symmetrical())
self.assertEqual(p.get_nplies_by_property_id(), 3)
self.assertAlmostEqual(p.get_thickness_by_property_id(pid), 0.6)
self.assertAlmostEqual(p.get_thickness_by_property_id_ply(pid, 0), 0.1)
self.assertAlmostEqual(p.get_thickness_by_property_id_ply(pid, 1), 0.2)
self.assertAlmostEqual(p.get_thickness_by_property_id_ply(pid, 2), 0.3)
with self.assertRaises(IndexError):
p.get_thickness_by_property_id_ply(pid, 3)
self.assertAlmostEqual(p.get_theta_by_property_id_ply(pid, 0), 0.)
self.assertAlmostEqual(p.get_theta_by_property_id_ply(pid, 1), 10.)
self.assertAlmostEqual(p.get_theta_by_property_id_ply(pid, 2), 20.)
with self.assertRaises(IndexError):
p.get_theta_by_property_id_ply(pid, 3)
self.assertEqual(p.get_material_id_by_property_id_ply(pid, 0), 1)
self.assertEqual(p.get_material_id_by_property_id_ply(pid, 1), 2)
self.assertEqual(p.get_material_id_by_property_id_ply(pid, 2), 3)
with self.assertRaises(IndexError):
p.get_material_id_by_property_id_ply(pid, 3)
#print('get_material_id_by_property_id = ', p.get_material_id_by_property_id(pid))
self.assertEqual(p.get_material_ids_by_property_id(pid)[0], 1)
self.assertEqual(p.get_material_ids_by_property_id(pid)[0], 1)
self.assertEqual(p.get_material_ids_by_property_id(pid)[1], 2)
self.assertEqual(p.get_material_ids_by_property_id(pid)[2], 3)
self.assertEqual(p.get_sout_by_property_id_ply(pid, 0), 'YES')
self.assertEqual(p.get_sout_by_property_id_ply(pid, 1), 'YES')
self.assertEqual(p.get_sout_by_property_id_ply(pid, 2), 'NO')
with self.assertRaises(IndexError):
p.get_sout_by_property_id_ply(pid, 3)
#for iply in range(len(p.plies)):
#mid = p.plies[iply][0]
#p.plies[iply][0] = m # MAT1
##p.mids = [m, m, m]
from StringIO import StringIO
f = StringIO()
#print(m.write_card(f, size=8, material_id=None))
p.write_card(f)
m.write_card(f)
print(f.getvalue())
#Mid
self.assertAlmostEqual(p.get_material_id_by_property_id_ply(pid, 0), 1)
self.assertAlmostEqual(p.get_material_id_by_property_id_ply(pid, 1), 2)
self.assertAlmostEqual(p.get_material_id_by_property_id_ply(pid, 2), 3)
with self.assertRaises(IndexError):
p.get_density_by_property_id_ply(pid, 3)
#Rho
self.assertAlmostEqual(p.get_density_by_property_id_ply(pid, 0), 1.0)
self.assertAlmostEqual(p.get_density_by_property_id_ply(pid, 1), 1.0)
self.assertAlmostEqual(p.get_density_by_property_id_ply(pid, 2), 1.0)
with self.assertRaises(IndexError):
p.get_density_by_property_id_ply(pid, 3)
# MassPerArea
self.assertAlmostEqual(p.get_mass_per_area_by_property_id(), 0.6)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id(pid), 0.6)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id([pid]), 0.6)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 0),
0.1)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 1),
0.2)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 2),
0.3)
#with self.assertRaises(IndexError):
#p.MassPerArea(3)
#----------------------
# change the nsm to 1.0
p.set_nonstructural_mass_by_property_id(pid, 1.0)
self.assertEqual(p.get_nonstructural_mass_by_property_id(), 1.0)
# MassPerArea
self.assertAlmostEqual(p.get_mass_per_area_by_property_id(), 1.6)
self.assertAlmostEqual(
p.get_mass_per_area_by_property_id_ply(pid, 0, method='nplies'),
0.1 + 1 / 3.)
self.assertAlmostEqual(
p.get_mass_per_area_by_property_id_ply(pid, 1, method='nplies'),
0.2 + 1 / 3.)
self.assertAlmostEqual(
p.get_mass_per_area_by_property_id_ply(pid, 2, method='nplies'),
0.3 + 1 / 3.)
self.assertAlmostEqual(
p.get_mass_per_area_by_property_id_ply(pid, 0, method='t'),
0.1 + 1 / 6.)
self.assertAlmostEqual(
p.get_mass_per_area_by_property_id_ply(pid, 1, method='t'),
0.2 + 2 / 6.)
self.assertAlmostEqual(
p.get_mass_per_area_by_property_id_ply(pid, 2, method='t'),
0.3 + 3 / 6.)
self.assertAlmostEqual(
p.get_mass_per_area_by_property_id_ply(pid, 0, method='rho*t'),
0.1 + 1 / 6.)
self.assertAlmostEqual(
p.get_mass_per_area_by_property_id_ply(pid, 1, method='rho*t'),
0.2 + 2 / 6.)
self.assertAlmostEqual(
p.get_mass_per_area_by_property_id_ply(pid, 2, method='rho*t'),
0.3 + 3 / 6.)
self.assertAlmostEqual(
p.get_mass_per_area_by_property_id_ply(pid, 0, method='t'),
0.1 + 1 / 6.)
self.assertAlmostEqual(
p.get_mass_per_area_by_property_id_ply(pid, 1, method='t'),
0.2 + 2 / 6.)
self.assertAlmostEqual(
p.get_mass_per_area_by_property_id_ply(pid, 2, method='t'),
0.3 + 3 / 6.)
with self.assertRaises(IndexError):
p.get_mass_per_area_by_property_id_ply(pid, 3, method='nplies')
z0 = p.get_z0_by_property_id(pid)
z = p.get_z_locations_by_property_id(pid)
z_expected = array([0., T[0], T[0] + T[1], T[0] + T[1] + T[2]])
for za, ze in zip(z, z_expected):
self.assertAlmostEqual(za, ze)
#z0 =
p.z0[0] = 1.0
z_expected = 1.0 + z_expected
z = p.get_z_locations_by_property_id(pid)
for za, ze in zip(z, z_expected):
self.assertAlmostEqual(za, ze)
def get_mass(self,nid1, nid2, xyz1, xyz2, eid, pid, mid, A, J, c, nsm, E, G, nu, rho, L):
"""tests a CROD and a CONROD"""
card_count = {
'CONROD' : 1,
'CTUBE' : 1,
'PTUBE' : 1,
'CROD' : 1,
'PROD' : 1,
'GRID' : 2,
'MAT1' : 1,
}
model = BDF(debug=debug)
model.allocate(card_count)
lines = ['conrod,%i, %i, %i, %i, %f, %f, %f, %f' % (eid, nid1, nid2, mid, A, J, c, nsm)]
model.add_card(lines, 'conrod', is_list=False)
lines = ['crod,%i, %i, %i, %i' % (eid+1, pid, nid1, nid2)]
model.add_card(lines, 'crod', is_list=False)
#lines = ['ctube,%i, %i, %i, %i' % (eid+2, pid+1, nid1, nid2)]
#model.add_card(lines, 'ctube', is_list=False)
lines = ['prod,%i, %i, %f, %f, %f, %f' % (pid, mid, A, J, c, nsm)]
model.add_card(lines, 'prod', is_list=False)
OD1 = sqrt(4*A/pi)
t = 0.
OD2 = OD1
#lines = ['ptube,%i, %i, %f, %f, %f, %f' % (pid+1, mid, OD1, t, nsm, OD2)]
#model.add_card(lines, 'ptube', is_list=False)
lines = ['mat1,%i, %.2e, %.2e, %f, %f' % (mid, E, G, nu, rho)]
model.add_card(lines, 'mat1', is_list=False)
lines = ['grid,%i, %i, %f, %f, %f' % (nid1, 0, xyz1[0], xyz1[1], xyz1[2])]
model.add_card(lines, 'grid', is_list=False)
lines = ['grid,%i, %i, %f, %f, %f' % (nid2, 0, xyz2[0], xyz2[1], xyz2[2])]
model.add_card(lines, 'grid', is_list=False)
model.cross_reference()
mass = L * (rho * A + nsm)
# conrod
conrod = model.conrod[eid]
self.assertEquals(conrod.get_element_id_by_element_index(), eid)
#self.assertEquals(conrod.get_property_id_by_element_index(), None)
self.assertEquals(conrod.get_material_id_by_element_index(), mid)
self.assertEquals(conrod.get_length_by_element_index(), L)
#self.assertEquals(conrod.Nsm(), nsm)
self.assertEquals(conrod.get_mass_by_element_index(), mass)
#self.assertEquals(conrod.E(), E)
#self.assertEquals(conrod.G(), G)
#self.assertEquals(conrod.Area(), A)
#self.assertEquals(conrod.J(), J)
#self.assertEquals(conrod.C(), c)
#self.assertEquals(conrod.Rho(), rho)
# crod
conrod = model.conrod[eid+1]
self.assertEquals(crod.get_element_id_by_element_index(), eid+1)
self.assertEquals(crod.get_property_id_by_element_index(), pid)
self.assertEquals(crod.get_material_id_by_element_index(), mid)
self.assertEquals(crod.get_length_by_element_index(), L)
#self.assertEquals(crod.Nsm(), nsm)
self.assertEquals(crod.get_mass_by_element_index(), mass)
#self.assertEquals(crod.E(), E)
#self.assertEquals(crod.G(), G)
#self.assertEquals(crod.Area(), A)
#self.assertEquals(crod.J(), J)
#self.assertEquals(crod.C(), c)
#self.assertEquals(crod.Rho(), rho)
#self.assertEquals(crod.Nu(), nu)
# prod
self.assertEquals(prod.Pid(), pid)
self.assertEquals(prod.Mid(), mid)
self.assertEquals(prod.Nsm(), nsm)
self.assertEquals(prod.E(), E)
self.assertEquals(prod.G(), G)
self.assertEquals(prod.Area(), A)
self.assertEquals(prod.J(), J)
self.assertEquals(prod.C(), c)
self.assertEquals(prod.Rho(), rho)
# ctube
if 0:
self.assertEquals(ctube.Eid(), eid+2)
self.assertEquals(ctube.Pid(), pid+1)
self.assertEquals(ctube.Mid(), mid)
self.assertEquals(ctube.Length(), L)
self.assertEquals(ctube.Nsm(), nsm)
self.assertAlmostEquals(ctube.Mass(), mass, 5)
self.assertEquals(ctube.E(), E)
self.assertEquals(ctube.G(), G)
self.assertAlmostEquals(ctube.Area(), A, 5)
ctube.J()
self.assertEquals(ctube.Rho(), rho)
# ptube
self.assertEquals(ptube.Pid(), pid+1)
self.assertEquals(ptube.Mid(), mid)
self.assertEquals(ptube.Nsm(), nsm)
self.assertEquals(ptube.E(), E)
self.assertEquals(ptube.G(), G)
self.assertAlmostEquals(ptube.Area(), A, 5)
ptube.J()
self.assertEquals(ptube.Rho(), rho)
def test_pcomp_02(self):
"""
symmetrical, nsm=0.0 and nsm=1.0
"""
model = BDF(debug=False)
pid = 1
z0 = 0.
nsm = 0.
sb = 0.
ft = 'HOFF'
TRef = 0.
ge = 0.
lam = 'SYM' # isSymmetrical SYM
Mid = [1, 2, 3]
Theta = [0., 10., 20.]
T = [.1, .2, .3]
Sout = ['YES', 'YES', 'NO'] # 0-NO, 1-YES
pcomp = [
'PCOMP', pid, z0, nsm, sb, ft, TRef, ge, lam, Mid[0], T[0],
Theta[0], Sout[0], Mid[1], T[1], Theta[1], Sout[1], Mid[2], T[2],
Theta[2], Sout[2]
]
#----
mid = 1
E = 3.0e7
G = None
nu = None
rho = 1.0
a = None
St = None
Sc = None
Ss = None
Mcsid = None
mat1_a = ['MAT1', mid, E, G, nu, rho, a, TRef, ge, St, Sc, Ss, Mcsid]
mat1_b = [
'MAT1', mid + 1, E, G, nu, rho, a, TRef, ge, St, Sc, Ss, Mcsid
]
mat1_c = [
'MAT1', mid + 2, E, G, nu, rho, a, TRef, ge, St, Sc, Ss, Mcsid
]
card_count = {
'PCOMP': 1,
'MAT1': 3,
}
model.allocate(card_count)
model.add_card(pcomp, pcomp[0], is_list=True)
model.add_card(mat1_a, 'MAT1', is_list=True)
model.add_card(mat1_b, 'MAT1', is_list=True)
model.add_card(mat1_c, 'MAT1', is_list=True)
model.build()
p = model.properties.properties_shell.pcomp[0]
self.assertTrue(p.is_symmetrical_by_property_id())
self.assertTrue(p.is_symmetrical_by_property_index())
self.assertEqual(p.get_nplies_by_property_id(), 6)
self.assertAlmostEqual(p.get_thickness_by_property_id()[0], 1.2)
self.assertAlmostEqual(p.get_thickness_by_property_id_ply(pid, 0), 0.1)
self.assertAlmostEqual(p.get_thickness_by_property_id_ply(pid, 1), 0.2)
self.assertAlmostEqual(p.get_thickness_by_property_id_ply(pid, 2), 0.3)
self.assertAlmostEqual(p.get_thickness_by_property_id_ply(pid, 3), 0.1)
self.assertAlmostEqual(p.get_thickness_by_property_id_ply(pid, 4), 0.2)
self.assertAlmostEqual(p.get_thickness_by_property_id_ply(pid, 5), 0.3)
with self.assertRaises(IndexError):
p.get_thickness_by_property_id_ply(pid, 6)
self.assertAlmostEqual(p.get_theta_by_property_id_ply(pid, 0), 0.)
self.assertAlmostEqual(p.get_theta_by_property_id_ply(pid, 1), 10.)
self.assertAlmostEqual(p.get_theta_by_property_id_ply(pid, 2), 20.)
self.assertAlmostEqual(p.get_theta_by_property_id_ply(pid, 3), 0.)
self.assertAlmostEqual(p.get_theta_by_property_id_ply(pid, 4), 10.)
self.assertAlmostEqual(p.get_theta_by_property_id_ply(pid, 5), 20.)
with self.assertRaises(IndexError):
p.get_theta_by_property_id_ply(pid, 6)
self.assertEqual(p.get_material_id_by_property_id_ply(pid, 0), 1)
self.assertEqual(p.get_material_id_by_property_id_ply(pid, 1), 2)
self.assertEqual(p.get_material_id_by_property_id_ply(pid, 2), 3)
self.assertEqual(p.get_material_id_by_property_id_ply(pid, 3), 1)
self.assertEqual(p.get_material_id_by_property_id_ply(pid, 4), 2)
self.assertEqual(p.get_material_id_by_property_id_ply(pid, 5), 3)
with self.assertRaises(IndexError):
p.get_material_id_by_property_id_ply(pid, 6)
self.assertTrue(
allclose(p.get_material_ids_by_property_id(pid),
[1, 2, 3, 1, 2, 3]))
self.assertEqual(p.get_sout_by_property_id_ply(pid, 0), 'YES')
self.assertEqual(p.get_sout_by_property_id_ply(pid, 1), 'YES')
self.assertEqual(p.get_sout_by_property_id_ply(pid, 2), 'NO')
self.assertEqual(p.get_sout_by_property_id_ply(pid, 3), 'YES')
self.assertEqual(p.get_sout_by_property_id_ply(pid, 4), 'YES')
self.assertEqual(p.get_sout_by_property_id_ply(pid, 5), 'NO')
with self.assertRaises(IndexError):
p.get_sout_by_property_id_ply(pid, 6)
#---------------
#Rho
self.assertAlmostEqual(p.get_density_by_property_id_ply(pid, 0), 1.0)
self.assertAlmostEqual(p.get_density_by_property_id_ply(pid, 1), 1.0)
self.assertAlmostEqual(p.get_density_by_property_id_ply(pid, 2), 1.0)
self.assertAlmostEqual(p.get_density_by_property_id_ply(pid, 3), 1.0)
self.assertAlmostEqual(p.get_density_by_property_id_ply(pid, 4), 1.0)
self.assertAlmostEqual(p.get_density_by_property_id_ply(pid, 5), 1.0)
with self.assertRaises(IndexError):
p.get_density_by_property_id_ply(pid, 6)
# MassPerArea
#self.assertAlmostEqual(p.get_mass_per_area_by_property_id(), 1.2)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id(pid), 1.2)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 0),
0.1)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 1),
0.2)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 2),
0.3)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 3),
0.1)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 4),
0.2)
self.assertAlmostEqual(p.get_mass_per_area_by_property_id_ply(pid, 5),
0.3)
with self.assertRaises(IndexError):
p.get_mass_per_area_by_property_id_ply(pid, 6)
self.assertEqual(p.get_nonstructural_mass_by_property_id(pid), 0.0)
#----------------------
# change the nsm to 1.0
p.set_nonstructural_mass_by_property_id(pid, 1.0)
self.assertEqual(p.get_nonstructural_mass_by_property_id(pid), 1.0)
# MassPerArea
self.assertAlmostEqual(p.get_mass_per_area_by_property_id(pid), 2.2)
self.assertAlmostEqual(
p.get_mass_per_area_by_property_id_ply(pid, 0, method='nplies'),
0.1 + 1 / 6.)
self.assertAlmostEqual(
p.get_mass_per_area_by_property_id_ply(pid, 1, method='nplies'),
0.2 + 1 / 6.)
self.assertAlmostEqual(
p.get_mass_per_area_by_property_id_ply(pid, 2, method='nplies'),
0.3 + 1 / 6.)
self.assertAlmostEqual(
p.get_mass_per_area_by_property_id_ply(pid, 3, method='nplies'),
0.1 + 1 / 6.)
self.assertAlmostEqual(
p.get_mass_per_area_by_property_id_ply(pid, 4, method='nplies'),
0.2 + 1 / 6.)
self.assertAlmostEqual(
p.get_mass_per_area_by_property_id_ply(pid, 5, method='nplies'),
0.3 + 1 / 6.)
with self.assertRaises(IndexError):
p.get_mass_per_area_by_property_id_ply(pid, 6)
