def test_celas4_cd(self):
"""Tests a CELAS4"""
model = BDF(debug=True, log=None, mode='msc')
model.bdf_filename = 'celas4.bdf'
model.add_spoint([1, 2])
#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_celas4(eid, k, nids, comment='')
load_id = 2
spc_id = 3
# model.add_sload(load_id, 2, 20.)
#fxyz = [1., 0., 0.]
mag = 20.
model.add_sload(load_id, 2, mag)
components = 0
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_celas3(self):
"""Tests a CELAS3/PELAS"""
model = BDF(debug=True, log=None, mode='msc')
model.bdf_filename = 'celas3.bdf'
#model.add_grid(1, [0., 0., 0.])
#model.add_grid(2, [0., 0., 0.])
model.add_spoint([1, 2])
nids = [1, 2]
eid = 1
pid = 2
model.add_celas3(eid, pid, nids, 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='')
nids = 2
mags = mag
model.add_sload(load_id, nids, mags, comment='')
components = 0
nodes = 1
model.add_spc1(spc_id, components, nodes, comment='')
setup_case_control(model)
solver = Solver(model)
#model.sol = 103
#solver.run()
model.sol = 101
solver.run()
# F = k * d
# d = F / k
# 20=1000.*d
d = mag / k
assert np.allclose(solver.xa_[0], d), solver.xa_
os.remove(solver.f06_filename)
os.remove(solver.op2_filename)
