Struct PyFemAssembly 

pub struct PyFemAssembly {
    pub n_dofs: usize,
    /* private fields */
}

FEM assembly for a truss structure.

Assembles a global stiffness matrix from bar elements, applies boundary conditions, and solves with a direct (dense) solver suitable for demonstration/testing with small meshes.

Fields

n_dofs: usize

Number of degrees of freedom (nodes × 3 for 3-D).

Implementations

impl PyFemAssembly

pub fn new(n_nodes: usize) -> Self

Create a new FEM assembly with n_nodes 3-D nodes.

pub fn add_bar_element( &mut self, node_a: usize, node_b: usize, ea: f64, length: f64, )

Add a bar element between two nodes with given EA and undeformed length.

pub fn apply_force(&mut self, dof: usize, force: f64)

Apply an external force to a DOF (node_index * 3 + direction).

pub fn fix_dof(&mut self, dof: usize)

Fix a DOF (zero displacement Dirichlet BC).

pub fn assemble(&mut self)

Assemble the global stiffness matrix from all bar elements.

For each bar element the 2×2 (1-D axial) stiffness matrix contribution k_e = EA/L * [[1, -1\], [-1, 1]] is added to the corresponding DOF rows/columns along the element’s axis direction.

pub fn solve(&mut self) -> bool

Solve the linear system K·u = f using Gaussian elimination.

The system is solved in-place. Returns true on success, false if the matrix is singular or not yet assembled.

pub fn displacement(&self, dof: usize) -> f64

Get the displacement at DOF dof, or 0.0 if out of bounds.

pub fn displacements(&self) -> &[f64]

Get all displacements as a slice.

pub fn element_force(&self, elem_idx: usize) -> Option<f64>

Compute the axial force in element elem_idx.

Returns None if the index is out of bounds or the system is unsolved.

Trait Implementations

impl Clone for PyFemAssembly

fn clone(&self) -> PyFemAssembly

Returns a duplicate of the value.

fn clone_from(&mut self, source: &Self)

Performs copy-assignment from source.

impl Debug for PyFemAssembly

fn fmt(&self, f: &mut Formatter<'_>) -> Result

Formats the value using the given formatter.