🔧 Model

Models are stored in the directory fibermat.model.

Degrees of Freedom

Mechanical model

stiffness

constraint

Example

from fibermat import *

# Generate a set of fibers
mat = Mat(100)
# Build the fiber network
net = Net(mat)
# Stack fibers
stack = Stack(net)
# Create the fiber mesh
mesh = Mesh(stack)
# Instantiate the model
model = Timoshenko(mesh)

# Permutation of indices
perm = sp.sparse.csgraph.reverse_cuthill_mckee(model.P, symmetric_mode=True)
# Visualize the system
fig, ax = plt.subplots(1, 2, figsize=(2 * 6.4, 4.8))
plot_system(model.stiffness(), model.constraint(), perm=None, ax=ax[0])
plot_system(model.stiffness(), model.constraint(), perm=perm, ax=ax[1])
plt.show()