def test_knownframe(self):
xs1 = xs.generalSection(A=1, Ix=1)
s1 = sc.Structure(xs1, ma.Steel(E=144))
s1.addNode(0, 0, fixity='fixed')
s1.addNode(24, 0, fixity='roller')
s1.addNode(42, 0, fixity='roller')
s1.addMember(0, 1)
s1.addMember(2, 1)
def test_bucklingCapacity(self):
xs1 = xs.AISC('W21X44')
s1 = sc.Structure(xs1, ma.Steel(Fy=60))
s1.addNode(0, 0)
s1.addNode(0, 15)
s1.addMember(0, 1)
member = s1.members[0]
# almost equal because rounding.
self.assertAlmostEqual(member.bucklingCapacity(1), 143.906124819)
## test another example
xs1 = xs.AISC('W16X67')
s1 = sc.Structure(xs1, ma.Steel(Fy=60))
s1.addNode(0, 0)
s1.addNode(0, 15)
s1.addMember(0, 1)
member = s1.members[0]
self.assertAlmostEqual(member.bucklingCapacity(1), 661.66119741611089)
def __init__(self, SN, EN, material=ma.Steel(), cross=xs.generalSection()):
self.cross = cross
self.material = material
# assign start and end node properties
self.SN = SN
self.EN = EN
#the axial force
self.axial = 0
def test_Frame(self):
# global stiffness matrix must be symmetric
xs1 = xs.generalSection(A=30, Ix=700)
s1 = sc.Structure(xs1, ma.Steel(E=29000))
s1.addNode(0, 0)
s1.addNode(0, 144)
s1.addNode(144, 144)
s1.addMember(0, 1)
s1.addMember(1, 2)
m1 = s1.members[0]
stiff1 = m1.T.T * m1.kframe * m1.T
v = np.all(stiff1 == stiff1.T)
self.assertEqual(v, True)
import cross_sections as xs import structural_classes as sc import materials as ma xs1 = xs.generalSection(A=30, Ix=700) s1 = sc.Structure(xs1, ma.Steel(E=29000)) s1.addNode(0, 0) s1.addNode(0, 144) s1.addNode(144, 144) s1.addMember(0, 1) s1.addMember(1, 2) m1 = s1.members[0] stiff1 = m1.T.T * m1.kframe * m1.T s1.plot()