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_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)