Beam6

A cantilever beam undergoing bending deformation in dynamic.

• 
• 
• 
• 

Generate movie 01/20 (5.00 %) 4.02 s
Generate movie 02/20 (10.00 %) 1.74 s
Generate movie 03/20 (15.00 %) 1.63 s
Generate movie 04/20 (20.00 %) 1.57 s
Generate movie 05/20 (25.00 %) 1.43 s
Generate movie 06/20 (30.00 %) 1.33 s
Generate movie 07/20 (35.00 %) 1.25 s
Generate movie 08/20 (40.00 %) 1.14 s
Generate movie 09/20 (45.00 %) 1.04 s
Generate movie 10/20 (50.00 %) 947.30 ms
Generate movie 11/20 (55.00 %) 862.89 ms
Generate movie 12/20 (60.00 %) 769.08 ms
Generate movie 13/20 (65.00 %) 665.48 ms
Generate movie 14/20 (70.00 %) 585.85 ms
Generate movie 15/20 (75.00 %) 480.77 ms
Generate movie 16/20 (80.00 %) 387.35 ms
Generate movie 17/20 (85.00 %) 283.17 ms
Generate movie 18/20 (90.00 %) 189.88 ms
Generate movie 19/20 (95.00 %) 95.04 ms
Generate movie 20/20 (100.00 %) 0.00 µs


==================== Mesh ====================

Element type: SEG3
Ne = 10, Nn = 21

==================== Model ====================

<EasyFEA.Models.Beam._beam.BeamStructure object at 0x7422441281d0>

solver:scipy

============= Boundary Conditions =============

Unspecified.

=================== Results ===================



Ux max = 0.00e+00
Ux min = 0.00e+00

Uy max = 0.00e+00
Uy min = -7.69e-01

=================== TicTac ===================

Mesh: 18.599 ms
Boundary Conditions: 20.504 µs
Matrix: 11.366 ms
Solver: 99.710 ms
Display: 14.532 ms
PyVista_Interface: 3.169 s

import matplotlib.pyplot as plt

from EasyFEA import (
    Folder,
    Display,
    Models,
    Mesher,
    ElemType,
    Simulations,
    Paraview,
    PyVista,
)
from EasyFEA.Geoms import Domain, Line

if __name__ == "__main__":
    Display.Clear()

    # ----------------------------------------------
    # Configuration
    # ----------------------------------------------

    # outputs
    folder = Folder.Results_Dir()
    makeParaview = False
    makeMovie = True
    result = "uy"

    # geom
    L = 120
    nL = 10
    h = 13
    b = 13
    e = 3

    # model
    E = 210000
    uy = 0.3

    # load
    load = 800

    # time
    Tmax = 2.0
    N = 50
    dt = Tmax / N

    # ----------------------------------------------
    # Mesh
    # ----------------------------------------------

    elemType = ElemType.SEG3
    beamDim = 2  # must be >= 2

    # Create a section object for the beam mesh
    domain1 = Domain((0, 0), (b, h), h / 10)
    domain2 = Domain((e, e), (b - e, h - e), h / 10)
    section = domain1.Mesh_2D([domain2])

    p1 = (0, 0)
    p2 = (L, 0)
    line = Line(p1, p2, L / nL)
    beam = Models.Beam.Isotropic(beamDim, line, section, E, uy)

    mesh = Mesher().Mesh_Beams([beam], elemType=elemType)

    # ----------------------------------------------
    # Simulation
    # ----------------------------------------------

    # Initialize the beam structure with the defined beam segments
    beamStructure = Models.Beam.BeamStructure([beam])

    # Create the beam simulation
    simu = Simulations.Beam(mesh, beamStructure)
    dof_n = simu.Get_dof_n()

    # Apply boundary conditions
    simu.add_dirichlet(mesh.Nodes_Point(p1), [0] * dof_n, simu.Get_unknowns())
    simu.add_neumann(mesh.Nodes_Point(p2), [-load], ["y"])

    # Solve the beam problem and get displacement results
    sol = simu.Solve()
    simu.Save_Iter()

    simu.Solver_Set_Hyperbolic_Algorithm(dt)
    simu.Bc_Init()
    simu.add_dirichlet(mesh.Nodes_Point(p1), [0] * dof_n, simu.Get_unknowns())

    for _ in range(N):
        simu.Solve()
        simu.Save_Iter()

    # ----------------------------------------------
    # Results
    # ----------------------------------------------

    Display.Plot_BoundaryConditions(simu)

    if makeParaview:
        Paraview.Save_simu(simu, folder)

    deform = 10
    if makeMovie:
        PyVista.Movie_simu(
            simu,
            f"{result}",
            folder,
            f"{result}.gif",
            N=20,
            deformFactor=deform,
            plotMesh=True,
        )

    print(simu)

    Display.Plot_Result(simu, result, deformFactor=deform)
    ax = Display.Plot_Mesh(section)
    ax.set_title("Section")
    Display.Plot_Mesh(simu, deformFactor=deform)

    plt.show()