def test_cquad4_pshell_mat1(self):
"""Tests a CQUAD4/PSHELL/MAT1"""
model = BDF(debug=True, log=None, mode='msc')
model.bdf_filename = 'cquad4_pshell_mat1.bdf'
model.add_grid(1, [0., 0., 0., ])
model.add_grid(2, [1., 0., 0., ])
model.add_grid(3, [1., 0., 2., ])
model.add_grid(4, [0., 0., 2., ])
area = 2.0
nsm_area = 0.0
thickness = 0.3
rho = 0.1
mass = area * (thickness * rho + nsm_area)
mid = 3
nids = [1, 2, 3, 4]
model.add_cquad4(1, 1, nids, theta_mcid=0.0, zoffset=0., tflag=0,
T1=None, T2=None, T3=None, T4=None, comment='')
model.add_cquad4(2, 2, nids, theta_mcid=0.0, zoffset=0., tflag=0,
T1=None, T2=None, T3=None, T4=None, comment='')
model.add_cquad4(3, 3, nids, theta_mcid=0.0, zoffset=0., tflag=0,
T1=None, T2=None, T3=None, T4=None, comment='')
model.add_cquad4(4, 4, nids, theta_mcid=0.0, zoffset=0., tflag=0,
T1=None, T2=None, T3=None, T4=None, comment='')
model.add_cquad4(5, 5, nids, theta_mcid=0.0, zoffset=0., tflag=0,
T1=None, T2=None, T3=None, T4=None, comment='')
model.add_pshell(1, mid1=mid, t=thickness, mid2=None, twelveIt3=1.0,
mid3=None, tst=0.833333, nsm=0.0, z1=None, z2=None,
mid4=None, comment='')
model.add_pshell(2, mid1=None, t=thickness, mid2=mid, twelveIt3=1.0,
mid3=None, tst=0.833333, nsm=0.0, z1=None, z2=None,
mid4=None, comment='')
model.add_pshell(3, mid1=None, t=thickness, mid2=None, twelveIt3=1.0,
mid3=mid, tst=0.833333, nsm=0.0, z1=None, z2=None,
mid4=None, comment='')
model.add_pshell(4, mid1=None, t=thickness, mid2=None, twelveIt3=1.0,
mid3=None, tst=0.833333, nsm=0.0, z1=None, z2=None,
mid4=mid, comment='')
model.add_pcomp(5, [mid], [thickness], thetas=None,
souts=None, nsm=0., sb=0., ft=None, tref=0., ge=0., lam=None, z0=None, comment='')
E = 1.0E7
G = None
nu = 0.3
model.add_mat1(mid, E, G, nu, rho=1.0, a=0.0, tref=0.0, ge=0.0,
St=0.0, Sc=0.0, Ss=0.0, mcsid=0, comment='')
spc_id = 3
load_id = 2
components = '123456'
nodes = [1, 2]
model.add_spc1(spc_id, components, nodes, comment='')
model.add_force(load_id, 3, 1.0, [0., 0., 1.], cid=0, comment='')
model.add_force(load_id, 4, 1.0, [0., 0., 1.], cid=0, comment='')
setup_case_control(model)
solver = Solver(model)
#with self.assertRaises(RuntimeError):
solver.run()
def test_cquad8_pshell_mat1(self):
"""Tests a CQUAD8/PSHELL/MAT1"""
# 4--7--3
# | |
# 8 6
# | |
# 1--5--2
model = BDF(debug=True, log=None, mode='msc')
model.bdf_filename = 'cquad8_pshell_mat1.bdf'
model.add_grid(1, [0., 0., 0.])
model.add_grid(2, [1., 0., 0.])
model.add_grid(3, [1., 0., 2.])
model.add_grid(4, [0., 0., 2.])
model.add_grid(5, [0.5, 0., 0.])
model.add_grid(6, [1., 0., 1.])
model.add_grid(7, [0.5, 0., 2.])
model.add_grid(8, [0., 0., 1.])
mid = 3
nids = [1, 2, 3, 4, 5, 6, 7, 8]
model.add_cquad8(1,
1,
nids,
theta_mcid=0.0,
zoffset=0.,
tflag=0,
T1=None,
T2=None,
T3=None,
T4=None,
comment='')
model.add_cquad8(2,
2,
nids,
theta_mcid=0.0,
zoffset=0.,
tflag=0,
T1=None,
T2=None,
T3=None,
T4=None,
comment='')
model.add_cquad8(3,
3,
nids,
theta_mcid=0.0,
zoffset=0.,
tflag=0,
T1=None,
T2=None,
T3=None,
T4=None,
comment='')
model.add_cquad8(4,
4,
nids,
theta_mcid=0.0,
zoffset=0.,
tflag=0,
T1=None,
T2=None,
T3=None,
T4=None,
comment='')
model.add_cquad8(5,
5,
nids,
theta_mcid=0.0,
zoffset=0.,
tflag=0,
T1=None,
T2=None,
T3=None,
T4=None,
comment='')
model.add_pshell(1,
mid1=mid,
t=0.3,
mid2=None,
twelveIt3=1.0,
mid3=None,
tst=0.833333,
nsm=0.0,
z1=None,
z2=None,
mid4=None,
comment='')
model.add_pshell(2,
mid1=None,
t=0.3,
mid2=mid,
twelveIt3=1.0,
mid3=None,
tst=0.833333,
nsm=0.0,
z1=None,
z2=None,
mid4=None,
comment='')
model.add_pshell(3,
mid1=None,
t=0.3,
mid2=None,
twelveIt3=1.0,
mid3=mid,
tst=0.833333,
nsm=0.0,
z1=None,
z2=None,
mid4=None,
comment='')
model.add_pshell(4,
mid1=None,
t=0.3,
mid2=None,
twelveIt3=1.0,
mid3=None,
tst=0.833333,
nsm=0.0,
z1=None,
z2=None,
mid4=mid,
comment='')
model.add_pcomp(5, [mid], [0.1],
thetas=None,
souts=None,
nsm=0.,
sb=0.,
ft=None,
tref=0.,
ge=0.,
lam=None,
z0=None,
comment='')
E = 1.0E7
G = None
nu = 0.3
model.add_mat1(mid,
E,
G,
nu,
rho=0.0,
a=0.0,
tref=0.0,
ge=0.0,
St=0.0,
Sc=0.0,
Ss=0.0,
mcsid=0,
comment='')
spc_id = 3
load_id = 2
components = '123456'
nodes = [1, 2, 5]
model.add_spc1(spc_id, components, nodes, comment='')
model.add_force(load_id, 3, 1.0, [0., 0., 1.], cid=0, comment='')
model.add_force(load_id, 4, 1.0, [0., 0., 1.], cid=0, comment='')
setup_case_control(model)
solver = Solver(model)
with self.assertRaises(RuntimeError):
solver.run()
