def test_cquad4_bad_normal(self):
"""test that the code crashes with a bad normal"""
model = BDF(debug=False, log=None, mode='msc')
model.bdf_filename = 'cquad4_bad_normal.bdf'
mid = 3
model.add_grid(1, [0., 0., 0.])
model.add_grid(3, [1., 0., 0.])
model.add_grid(2, [1., 0., 2.])
model.add_grid(4, [0., 0., 2.])
nids = [1, 2, 3, 4]
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_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]
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):
# invalid normal vector
solver.run()
def test_celas1(self):
"""Tests a CELAS1/PELAS"""
model = BDF(debug=True, log=None, mode='msc')
model.add_grid(1, [0., 0., 0.])
model.add_grid(2, [0., 0., 0.])
nids = [1, 2]
eid = 1
pid = 2
model.add_celas1(eid, pid, nids, c1=1, c2=1, comment='')
k = 1000.
load_id = 2
spc_id = 3
model.add_pelas(pid, k, ge=0., s=0., comment='')
# model.add_sload(load_id, 2, 20.)
fxyz = [1., 0., 0.]
mag = 20.
model.add_force(load_id, 2, mag, fxyz, cid=0, comment='')
components = 123456
nodes = 1
model.add_spc1(spc_id, components, nodes, comment='')
setup_case_control(model)
solver = Solver(model)
solver.run()
# F = k * d
# d = F / k
# 20=1000.*d
d = mag / k
assert np.allclose(solver.xa_[0], d)
def test_celas2_cd(self):
"""Tests a CELAS2"""
model = BDF(debug=True, log=None, mode='msc')
model.bdf_filename = 'celas2.bdf'
model.add_grid(1, [0., 0., 0.])
cd = 100
model.add_grid(2, [0., 0., 0.], cd=cd)
origin = [0., 0., 0.]
zaxis = [0., 0., 1.]
xzplane = [0., 1., 0.]
model.add_cord2r(cd, origin, zaxis, xzplane, rid=0, setup=True, comment='')
nids = [1, 2]
eid = 1
k = 1000.
model.add_celas2(eid, k, nids, c1=1, c2=2, ge=0., s=0., comment='')
load_id = 2
spc_id = 3
# model.add_sload(load_id, 2, 20.)
fxyz = [1., 0., 0.]
mag = 20.
model.add_force(load_id, 2, mag, fxyz, cid=0, comment='')
components = 123456
nodes = 1
model.add_spc1(spc_id, components, nodes, comment='')
setup_case_control(model)
solver = Solver(model)
solver.run()
# F = k * d
d = mag / k
assert np.allclose(solver.xa_[0], d)
def test_ctube(self):
"""Tests a CTUBE/PTUBE"""
model = BDF(debug=True, log=None, mode='msc')
model.add_grid(1, [0., 0., 0.])
model.add_grid(2, [1., 0., 0.])
nids = [1, 2]
eid = 1
pid = 2
mid = 3
E = 3.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)
model.add_ctube(eid, pid, nids)
OD1 = 1.0
model.add_ptube(pid, mid, OD1, t=0.1, nsm=0., OD2=None, comment='')
load_id = 2
spc_id = 3
nid = 2
mag = 1.
fxyz = [1., 0., 0.]
model.add_force(load_id, nid, mag, fxyz, cid=0)
components = 123456
nodes = 1
model.add_spc1(spc_id, components, nodes, comment='')
setup_case_control(model)
solver = Solver(model)
solver.run()
def test_crod(self):
"""Tests a CROD/PROD"""
model = BDF(debug=True, log=None, mode='msc')
model.bdf_filename = 'crod.bdf'
model.add_grid(1, [0., 0., 0.])
model.add_grid(2, [1., 0., 0.])
nids = [1, 2]
eid = 1
pid = 2
mid = 3
E = 3.0e7
G = None
nu = 0.3
model.add_mat1(mid, E, G, nu, rho=0.1, a=0.0, tref=0.0, ge=0.0, St=0.0,
Sc=0.0, Ss=0.0, mcsid=0)
model.add_crod(eid, pid, nids)
model.add_prod(pid, mid, A=1.0, j=2., c=0., nsm=0.)
load_id = 2
spc_id = 3
nid = 2
mag = 1.
fxyz = [1., 0., 0.]
model.add_force(load_id, nid, mag, fxyz, cid=0)
components = 123456
nodes = 1
model.add_spc1(spc_id, components, nodes, comment='')
setup_case_control(model)
solver = Solver(model)
solver.run()
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_cbeam(self):
"""Tests a CBEAM/PBEAM"""
model = BDF(debug=True, log=None, mode='msc')
model.bdf_filename = 'cbeam.bdf'
model.add_grid(1, [0., 0., 0.])
model.add_grid(2, [1., 0., 0.])
L = 1.0
nids = [1, 2]
eid = 1
pid = 2
mid = 3
E = 3.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)
x = [0., 1., 0.]
g0 = None
model.add_cbeam(eid, pid, nids, x, g0, offt='GGG', bit=None,
pa=0, pb=0, wa=None, wb=None, sa=0, sb=0, comment='')
#beam_type = 'BAR'
xxb = [0.]
#dims = [[1., 2.]]
#model.add_pbeaml(pid, mid, beam_type, xxb, dims,
#so=None, nsm=None, group='MSCBML0', comment='')
A = 1.0
so = ['YES']
area = [1.0]
i1 = [1.]
i2 = [2.]
i12 = [0.]
j = [2.]
model.add_pbeam(pid, mid, xxb, so, area, i1, i2, i12, j, nsm=None,
c1=None, c2=None, d1=None, d2=None, e1=None, e2=None, f1=None, f2=None,
k1=1., k2=1., s1=0., s2=0., nsia=0., nsib=None, cwa=0., cwb=None,
m1a=0., m2a=0., m1b=None, m2b=None, n1a=0., n2a=0., n1b=None, n2b=None,
comment='')
load_id = 2
spc_id = 3
nid = 2
mag = 1.
fxyz = [1., 0., 0.]
model.add_force(load_id, nid, mag, fxyz, cid=0)
components = 123456
nodes = 1
model.add_spc1(spc_id, components, nodes, comment='')
setup_case_control(model)
solver = Solver(model)
solver.run()
# F = kx
k_axial = A * E / L
fmag = 1.0
dx = fmag / k_axial
assert dx == solver.xg[6], f'dx={dx} xg={xg}'
def test_crod_spcd(self):
"""Tests a CROD/PROD with an SPCD and no free DOFs"""
model = BDF(debug=True, log=None, mode='msc')
model.bdf_filename = 'crod_spcd.bdf'
model.add_grid(1, [0., 0., 0.])
model.add_grid(2, [1., 0., 0.])
nids = [1, 2]
eid = 1
pid = 2
mid = 3
E = 42.
G = None
nu = 0.3
model.add_mat1(mid, E, G, nu, rho=0.1, a=0.0, tref=0.0, ge=0.0, St=0.0,
Sc=0.0, Ss=0.0, mcsid=0)
model.add_crod(eid, pid, nids)
model.add_prod(pid, mid, A=1.0, j=0., c=0., nsm=0.)
load_id = 2
spc_id = 3
nid = 2
mag = 1.
fxyz = [0., 0., 0.]
model.add_force(load_id, nid, mag, fxyz, cid=0)
nodes = 1
components = 123456
#model.add_spc1(spc_id, components, nodes, comment='')
model.add_spc(spc_id, nodes, components, 0., comment='')
#components = 123456
#nodes = 1
#enforced = 0.1
#model.add_spc(spc_id, nodes, components, enforced, comment='')
nodes = 2
components = 1
enforced = 0.1
model.add_spcd(load_id, nodes, components, enforced, comment='')
components = '1'
model.add_spcd(9999999, nodes, components, enforced, comment='')
components = 123456
#model.add_spc1(spc_id, components, nodes, comment='')
model.add_spc(spc_id, nodes, components, 0., comment='')
setup_case_control(model)
solver = Solver(model)
solver.run()
# F = kx
dx = enforced
k = E
F = k * dx
assert len(solver.xa_) == 0, solver.xa_
assert len(solver.Fa_) == 0, solver.Fa_
assert np.allclose(solver.xg[6], dx), solver.xa_
assert np.allclose(solver.Fg[6], F), solver.xa_
def test_conm2(self):
"""Tests a CMASS1/PMASS"""
model = BDF(debug=True, log=None, mode='msc')
model.bdf_filename = 'celas1.bdf'
model.add_grid(1, [0., 0., 0.], cd=1)
model.add_grid(2, [1., 0., 0.])
model.add_grid(3, [0.5, 1., 0.], cd=3)
model.add_grid(4, [0.5, 0., 1.])
origin = [0., 0., 0.]
zaxis = [0., 0., 1]
xzplane = [1., 1., 0.]
model.add_cord2r(1,
origin,
zaxis,
xzplane,
rid=0,
setup=True,
comment='')
origin = [0.5, 1., 0.]
zaxis = [0., 0., 1]
xzplane = [0.5, 2., 0.]
model.add_cord2r(3,
origin,
zaxis,
xzplane,
rid=0,
setup=True,
comment='')
model.add_conm2(1, 1, mass=2.0, cid=0, X=None, I=None, comment='')
model.add_conm2(2, 2, mass=3.0, cid=0, X=None, I=None, comment='')
model.add_conm2(3, 3, mass=3.0, cid=0, X=None, I=None, comment='')
model.add_conm2(4, 4, mass=5.0, cid=0, X=None, I=None, comment='')
nids = [1, 2]
eid = 1
k = 1000.
model.add_celas2(eid, k, nids, c1=1, c2=2, ge=0., s=0., comment='')
load_id = 2
spc_id = 3
# model.add_sload(load_id, 2, 20.)
fxyz = [1., 0., 0.]
mag = 20.
model.add_force(load_id, 2, mag, fxyz, cid=0, comment='')
components = 123456
nodes = 1
model.add_spc1(spc_id, components, nodes, comment='')
setup_case_control(model)
solver = Solver(model)
model.sol = 101
solver.run()
def test_cbar2(self):
"""Tests a CBAR/PBAR"""
model = BDF(debug=True, log=None, mode='msc')
model.bdf_filename = 'cbar.bdf'
model.add_grid(1, [0., 0., 0.])
model.add_grid(2, [0.5, 0., 0.])
model.add_grid(3, [1., 0., 0.])
L = 1.0
pid = 2
mid = 3
E = 3.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)
eid = 1
x = [0., 1., 0.]
g0 = None
nids = [1, 2]
model.add_cbar(eid, pid, nids, x, g0,
offt='GGG', pa=0, pb=0, wa=None, wb=None, comment='')
eid = 2
nids = [2, 3]
model.add_cbar(eid, pid, nids, x, g0,
offt='GGG', pa=0, pb=0, wa=None, wb=None, comment='')
A = 1.
k_axial = A * E / L
model.add_pbar(pid, mid,
A=A, i1=1., i2=1., i12=1., j=1.,
nsm=0.,
c1=0., c2=0., d1=0., d2=0.,
e1=0., e2=0., f1=0., f2=0.,
k1=1.e8, k2=1.e8, comment='')
load_id = 2
spc_id = 3
nid = 3
mag = 1.
fxyz = [1., 0., 0.]
model.add_force(load_id, nid, mag, fxyz, cid=0)
components = 123456
nodes = 1
model.add_spc1(spc_id, components, nodes, comment='')
setup_case_control(model)
solver = Solver(model)
solver.run()
# F = kx
fmag = 1.0
dx = fmag / k_axial
assert dx == solver.xg[6*2], f'dx={dx} xg={xg}'
def test_ctube(self):
"""Tests a CTUBE/PTUBE"""
model = BDF(debug=True, log=None, mode='msc')
model.bdf_filename = 'ctube.bdf'
model.add_grid(1, [0., 0., 0.])
model.add_grid(2, [1., 0., 0.])
nids = [1, 2]
eid = 1
pid = 2
mid = 3
E = 3.0e7
G = None
nu = 0.3
model.add_mat1(mid,
E,
G,
nu,
rho=0.1,
a=0.0,
tref=0.0,
ge=0.0,
St=0.0,
Sc=0.0,
Ss=0.0,
mcsid=0)
model.add_ctube(eid, pid, nids)
OD1 = 1.0
t = 0.1
model.add_ptube(pid, mid, OD1, t=t, nsm=0., OD2=None, comment='')
load_id = 2
spc_id = 3
nid = 2
mag = 1.
fxyz = [1., 0., 0.]
model.add_force(load_id, nid, mag, fxyz, cid=0)
components = 123456
nodes = 1
model.add_spc1(spc_id, components, nodes, comment='')
setup_case_control(model)
solver = Solver(model)
solver.run()
ID = OD1 - 2 * t
# F = k * x
# x = F / k
L = 1.0
A = np.pi * (OD1**2 - ID**2) / 4
kaxial = A * E / L
F = 1.0
dx = F / kaxial
assert np.allclose(solver.xg[6], dx), f'dx={dx} xg[6]={solver.xg[6]}'
assert np.allclose(solver.Fg[6], F), solver.Fg
def test_conrod(self):
"""Tests a CONROD"""
model = BDF(debug=True, log=None, mode='msc')
L = 1.
model.add_grid(1, [0., 0., 0.])
model.add_grid(2, [L, 0., 0.])
nids = [1, 2]
eid = 1
mid = 3
E = 200.
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)
A = 1.0
J = 2.0
model.add_conrod(eid, mid, nids, A=A, j=J, c=0.0, nsm=0.0)
load_id = 2
spc_id = 3
nid = 2
mag_axial = 10000.
fxyz = [1., 0., 0.]
model.add_force(load_id, nid, mag_axial, fxyz, cid=0)
mag_torsion = 20000.
fxyz = [1., 0., 0.]
model.add_moment(load_id, nid, mag_torsion, fxyz, cid=0)
components = 123456
nodes = 1
model.add_spc1(spc_id, components, nodes, comment='')
#nodes = 2
#model.add_spc1(spc_id, 23456, nodes, comment='')
setup_case_control(model)
solver = Solver(model)
solver.run()
# TODO: why is this a torsional result???
G = E / (2 * (1 + nu))
kaxial = A * E / L
ktorsion = G * J / L
# F = k * d
daxial = mag_axial / kaxial
dtorsion = mag_torsion / ktorsion
assert np.allclose(solver.xa_[0], daxial), f'daxial={daxial} kaxial={kaxial} xa_={solver.xa_}'
assert np.allclose(solver.xa_[3], dtorsion), f'dtorsion={dtorsion} ktorsion={ktorsion} xa_={solver.xa_}'
def test_crod_mpc(self):
"""Tests a CROD/PROD"""
model = BDF(debug=True, log=None, mode='msc')
model.add_grid(1, [0., 0., 0.])
model.add_grid(2, [1., 0., 0.])
model.add_grid(3, [1., 0., 0.])
nids = [1, 2]
eid = 1
pid = 2
mid = 3
E = 3.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)
model.add_crod(eid, pid, nids)
model.add_prod(pid, mid, A=1.0, j=2., c=0., nsm=0.)
mpc_id = 10
nodes = [2, 3]
components = [1, 3]
coefficients = [1., -1.]
model.add_mpc(mpc_id, nodes, components, coefficients)
load_id = 2
spc_id = 3
nid = 3
mag = 1.
fxyz = [0., 0., 1.]
model.add_force(load_id, nid, mag, fxyz, cid=0)
components = 123456
nodes = 1
model.add_spc1(spc_id, components, nodes, comment='')
setup_case_control(model, extra_case_lines=['MPC=10'])
solver = Solver(model)
solver.run()
def test_crod_aset(self):
"""
Tests a CROD/PROD using an ASET
same answer as ``test_crod``
"""
model = BDF(debug=True, log=None, mode='msc')
model.add_grid(1, [0., 0., 0.])
model.add_grid(2, [1., 0., 0.])
nids = [1, 2]
eid = 1
pid = 2
mid = 3
E = 3.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)
model.add_crod(eid, pid, nids)
model.add_prod(pid, mid, A=1.0, j=2., c=0., nsm=0.)
load_id = 2
spc_id = 3
nid = 2
mag = 1.
fxyz = [1., 0., 0.]
model.add_force(load_id, nid, mag, fxyz, cid=0)
components = 123456
nodes = 1
model.add_spc1(spc_id, components, nodes, comment='')
aset_ids = 2
aset_components = '456'
model.add_aset(aset_ids, aset_components)
aset_ids = [2, 2]
aset_components = ['456', '123']
model.add_aset1(aset_ids, aset_components, comment='')
setup_case_control(model)
solver = Solver(model)
solver.run()
def test_cbar(self):
"""Tests a CBAR/PBAR"""
model = BDF(debug=True, log=None, mode='msc')
model.add_grid(1, [0., 0., 0.])
model.add_grid(2, [1., 0., 0.])
nids = [1, 2]
eid = 1
pid = 2
mid = 3
E = 3.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)
x = [0., 1., 0.]
g0 = None
model.add_cbar(eid, pid, nids, x, g0,
offt='GGG', pa=0, pb=0, wa=None, wb=None, comment='')
model.add_pbar(pid, mid,
A=1., i1=1., i2=1., i12=1., j=1.,
nsm=0.,
c1=0., c2=0., d1=0., d2=0.,
e1=0., e2=0., f1=0., f2=0.,
k1=1.e8, k2=1.e8, comment='')
load_id = 2
spc_id = 3
nid = 2
mag = 1.
fxyz = [1., 0., 0.]
model.add_force(load_id, nid, mag, fxyz, cid=0)
components = 123456
nodes = 1
model.add_spc1(spc_id, components, nodes, comment='')
setup_case_control(model)
solver = Solver(model)
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()
def test_cquad4_bad_jacobian(self):
"""
Tests that the code crashes with a really terrible CQUAD4
The Jacobian is defined between [-1, 1]
"""
model = BDF(debug=None, log=None, mode='msc')
model.bdf_filename = 'cquad4_bad_jacobian.bdf'
mid = 3
model.add_grid(1, [0., 0., 0.])
model.add_grid(2, [0.5, 100., 0.])
model.add_grid(3, [1., 0., 0.])
model.add_grid(4, [0.5, 1., 0.])
nids = [1, 2, 3, 4]
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_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]
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()
