Tutorial

import truss
import numpy as np
import matplotlib.pyplot as plt
from matplotlib import cm
from matplotlib import patches

# CONSTANTS
modulus = 210.e9     #Pa
rho = 2700.          #kg/m**3
surface = .0001       #m**2
yield_stress = 400.e6   #Pa
from scipy import optimize

# MODEL
m = truss.core.Model()

# NODES
A = m.add_node((0.,0.), label = "A")
B = m.add_node((0.,1.), label = "B")
C = m.add_node((1.,0.), label = "C")
D = m.add_node((1.,1.), label = "D")

# BOUNDARY CONDITIONS
A.block[0] = True
A.block[1] = True
B.block[0] = True

#EXTERNAL FORCES
C.force = np.array([1., 1.])*1.e3

#BARS
AB = m.add_bar(A, B, modulus = modulus, density = rho, section = surface)
AC = m.add_bar(A, C, modulus = modulus, density = rho, section = surface)
BC = m.add_bar(B, C, modulus = modulus, density = rho, section = surface)
CD = m.add_bar(C, D, modulus = modulus, density = rho, section = surface)
BD = m.add_bar(B, D, modulus = modulus, density = rho, section = surface)

#SOLVING
m.solve()

#PLOTTING
xlim, ylim = m.bbox(deformed = True)
fig = plt.figure(0)
plt.clf()
ax = fig.add_subplot(1,1,1)
ax.set_aspect("equal")
#ax.axis("off")
m.draw(ax, deformed = True, field = "stress", label = True, force_scale = 1.e-4)
plt.xlim(xlim)
plt.ylim(ylim)
plt.grid()
plt.show()

(Source code, png, hires.png, pdf)