Beam7

Frame with two beams in dynamic.

• 
• 
• 
• 

The beam's vertical axis has been selected incorrectly (collinear with the beam x-axis).
Axis [-1.  0.  0.] has been assigned for beam12.

Generate movie 01/20 (5.00 %) 2.38 s
Generate movie 02/20 (10.00 %) 1.77 s
Generate movie 03/20 (15.00 %) 1.67 s
Generate movie 04/20 (20.00 %) 1.58 s
Generate movie 05/20 (25.00 %) 1.44 s
Generate movie 06/20 (30.00 %) 1.36 s
Generate movie 07/20 (35.00 %) 1.27 s
Generate movie 08/20 (40.00 %) 1.20 s
Generate movie 09/20 (45.00 %) 1.08 s
Generate movie 10/20 (50.00 %) 995.32 ms
Generate movie 11/20 (55.00 %) 880.84 ms
Generate movie 12/20 (60.00 %) 778.02 ms
Generate movie 13/20 (65.00 %) 678.31 ms
Generate movie 14/20 (70.00 %) 586.86 ms
Generate movie 15/20 (75.00 %) 486.20 ms
Generate movie 16/20 (80.00 %) 390.55 ms
Generate movie 17/20 (85.00 %) 298.77 ms
Generate movie 18/20 (90.00 %) 197.34 ms
Generate movie 19/20 (95.00 %) 97.91 ms
Generate movie 20/20 (100.00 %) 0.00 µs


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

Element type: SEG3
Ne = 14, Nn = 30

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

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

solver:scipy

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

Unspecified.

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



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

Uy max = 1.34e-02
Uy min = -3.88e-01

Uz max = 0.00e+00
Uz min = -1.95e+00

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

Mesh: 15.069 ms
Boundary Conditions: 161.886 µs
Matrix: 24.778 ms
Solver: 157.148 ms
Display: 16.139 ms
PyVista_Interface: 3.206 s

import matplotlib.pyplot as plt

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

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

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

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

    # geom
    L = 120
    h = 13
    b = 13

    # model
    E = 210000
    v = 0.3

    # load
    load = 800

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

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

    elemType = ElemType.SEG3

    # Create a section object for the beam mesh
    mesher = Mesher()
    section = mesher.Mesh_2D(Circle((0, 0), h))

    p1 = (0, 0)
    p2 = (0, L)
    p3 = (L / 2, L)
    line1 = Line(p1, p2, L / 9)
    line2 = Line(p2, p3, L / 9)
    beam1 = Models.Beam.Isotropic(3, line1, section, E, v)
    beam2 = Models.Beam.Isotropic(3, line2, section, E, v)
    beams = [beam1, beam2]

    mesh = mesher.Mesh_Beams(beams=beams, elemType=elemType)

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

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

    # 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(p3), [-load, load], ["y", "z"])
    if beamStructure.nBeam > 1:
        simu.add_connection_fixed(mesh.Nodes_Point(p2))

    # 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())
    if beamStructure.nBeam > 1:
        simu.add_connection_fixed(mesh.Nodes_Point(p2))

    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()