Elemento Bar

In [39]:

    

"""
2.41 (Beer)
"""

from nusa import *
from math import pi

# Input data 
Es = 200e9
As = (pi/4.)*(40e-3)**2
Eb = 105e9
Ab = (pi/4.)*(30e-3)**2
Ps = 60e3
Pb = 40e3

# Model
m1 = BarModel("Bar Model")
# Nodes
n1 = Node((0,0))
n2 = Node((0.180,0))
n3 = Node((0.180+0.120,0))
n4 = Node((0.180+0.120+0.100,0))
n5 = Node((0.180+0.120+0.200,0))
# Elements
e1 = Bar((n1,n2),Es,As)
e2 = Bar((n2,n3),Es,As)
e3 = Bar((n3,n4),Eb,Ab)
e4 = Bar((n4,n5),Eb,Ab)

# Add elements 
for nd in (n1,n2,n3,n4,n5):
    m1.add_node(nd)
for el in (e1,e2,e3,e4):
    m1.add_element(el)

m1.add_constraint(n1,ux=0) # fixed 
m1.add_constraint(n5,ux=0) # fixed
m1.add_force(n2,(Ps,))
m1.add_force(n4,(Pb,))

m1.solve() # Solve model

# Reacciones en A y E
print("RA = {0} N".format(n1.fx))
print("RE = {0} N".format(n5.fx))

# Deflexión en C
print("Desplazamiento en C = {0} m".format(n3.ux))


    

    




RA = -62808.87926367081 N
RE = -37191.12073632918 N
Desplazamiento en C = 4.6324687280023236e-05 m

In [ ]:

    




    

In [ ]: