LinearElasticity2

A cantilever beam undergoing bending deformation in dynamic.

Note that this simulation is also performed in examples/Dynamic/Dynamic1.py.

Static Scene

Interactive Scene

0.000 s
0.010 s
0.020 s
0.030 s
0.040 s
0.050 s
0.060 s
0.070 s
0.080 s
0.090 s
0.100 s
0.110 s
0.120 s
0.130 s
0.140 s
0.150 s
0.160 s
0.170 s
0.180 s
0.190 s
0.200 s
0.210 s
0.220 s
0.230 s
0.240 s
0.250 s
0.260 s
0.270 s
0.280 s
0.290 s
0.300 s
0.310 s
0.320 s
0.330 s
0.340 s
0.350 s
0.360 s
0.370 s
0.380 s
0.390 s
0.400 s
0.410 s
0.420 s
0.430 s
0.440 s
0.450 s
0.460 s
0.470 s
0.480 s
0.490 s
0.500 s
Generate movie 0/51
Generate movie 1/51 (1.96 %) 10.36 s
Generate movie 2/51 (3.92 %) 10.03 s
Generate movie 3/51 (5.88 %) 10.03 s
Generate movie 4/51 (7.84 %) 9.68 s
Generate movie 5/51 (9.80 %) 9.29 s
Generate movie 6/51 (11.76 %) 8.93 s
Generate movie 7/51 (13.73 %) 9.19 s
Generate movie 8/51 (15.69 %) 8.93 s
Generate movie 9/51 (17.65 %) 8.76 s
Generate movie 10/51 (19.61 %) 8.56 s
Generate movie 11/51 (21.57 %) 8.29 s
Generate movie 12/51 (23.53 %) 8.07 s
Generate movie 13/51 (25.49 %) 7.92 s
Generate movie 14/51 (27.45 %) 7.81 s
Generate movie 15/51 (29.41 %) 7.56 s
Generate movie 16/51 (31.37 %) 7.28 s
Generate movie 17/51 (33.33 %) 7.10 s
Generate movie 18/51 (35.29 %) 6.89 s
Generate movie 19/51 (37.25 %) 6.66 s
Generate movie 20/51 (39.22 %) 6.45 s
Generate movie 21/51 (41.18 %) 6.04 s
Generate movie 22/51 (43.14 %) 6.02 s
Generate movie 23/51 (45.10 %) 5.66 s
Generate movie 24/51 (47.06 %) 5.61 s
Generate movie 25/51 (49.02 %) 5.18 s
Generate movie 26/51 (50.98 %) 5.09 s
Generate movie 27/51 (52.94 %) 4.94 s
Generate movie 28/51 (54.90 %) 4.79 s
Generate movie 29/51 (56.86 %) 4.55 s
Generate movie 30/51 (58.82 %) 4.28 s
Generate movie 31/51 (60.78 %) 4.08 s
Generate movie 32/51 (62.75 %) 4.05 s
Generate movie 33/51 (64.71 %) 3.73 s
Generate movie 34/51 (66.67 %) 3.53 s
Generate movie 35/51 (68.63 %) 3.32 s
Generate movie 36/51 (70.59 %) 3.12 s
Generate movie 37/51 (72.55 %) 2.93 s
Generate movie 38/51 (74.51 %) 2.72 s
Generate movie 39/51 (76.47 %) 2.49 s
Generate movie 40/51 (78.43 %) 2.24 s
Generate movie 41/51 (80.39 %) 2.03 s
Generate movie 42/51 (82.35 %) 1.83 s
Generate movie 43/51 (84.31 %) 1.66 s
Generate movie 44/51 (86.27 %) 1.43 s
Generate movie 45/51 (88.24 %) 1.23 s
Generate movie 46/51 (90.20 %) 1.02 s
Generate movie 47/51 (92.16 %) 827.89 ms
Generate movie 48/51 (94.12 %) 609.39 ms
Generate movie 49/51 (96.08 %) 405.91 ms
Generate movie 50/51 (98.04 %) 199.98 ms
Generate movie 51/51 (100.00 %) 0.00 µs

from EasyFEA import Folder, Display, ElemType, Models, Simulations, PyVista, np
from EasyFEA.fem import Field, BiLinearForm, FeArray, Sym_Grad, Trace
from EasyFEA.Geoms import Domain

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

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

    dim = 3

    # outputs
    makeMovie = True
    folder = Folder.Join(Folder.RESULTS_DIR, "WeakForms", "LinearElasticity2")

    # geom
    L = 120  # mm
    h = 13

    # model
    elastic = Models.ElasIsot(dim, 210000, 0.3, planeStress=True, thickness=h)
    lmbda = elastic.get_lambda()
    mu = elastic.get_mu()
    rho = 8100 * 1e-9

    # load
    Tmax = 0.5
    N = 50
    dt = Tmax / N
    time = -dt

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

    contour = Domain((0, 0), (L, h), h / 3)

    if dim == 2:
        mesh = contour.Mesh_2D([], ElemType.QUAD4, isOrganised=True)
    else:
        mesh = contour.Mesh_Extrude(
            [], [0, 0, h], [3], ElemType.HEXA8, isOrganised=True
        )

    nodesX0 = mesh.Nodes_Conditions(lambda x, y, z: x == 0)
    nodesXL = mesh.Nodes_Conditions(lambda x, y, z: x == L)

    # ----------------------------------------------
    # Formulations
    # ----------------------------------------------

    field = Field(mesh.groupElem, dim)

    def S(u: Field) -> FeArray:
        Eps = Sym_Grad(u)
        return 2 * mu * Eps + lmbda * Trace(Eps) * np.eye(dim)

    @BiLinearForm
    def computeK(u: Field, v: Field):
        Sig = S(u)
        Eps = Sym_Grad(v)
        return Sig.ddot(Eps)

    @BiLinearForm
    def computeM(u: Field, v: Field):
        return rho * u.dot(v)

    @BiLinearForm
    def computeC(u: Field, v: Field):
        K = computeK(u, v)
        M = computeM(u, v)
        C = K * 1e-3 + M * 1e-3
        return C

    weakForms = Models.WeakForms(field, computeK, computeC, computeM)

    # ----------------------------------------------
    # Simulations
    # ----------------------------------------------

    simu = Simulations.WeakFormSimu(mesh, weakForms)

    # static simulation
    simu.add_dirichlet(nodesX0, [0] * dim, simu.Get_unknowns())
    simu.add_dirichlet(nodesXL, [-10], ["y"])

    simu.Solve()

    # dynamic simulation
    simu.Solver_Set_Hyperbolic_Algorithm(dt)
    simu.Save_Iter()

    simu.Bc_Init()
    simu.add_dirichlet(nodesX0, [0] * dim, simu.Get_unknowns())

    while time <= Tmax:

        time += dt

        simu.Solve()
        simu.Save_Iter()

        print(f"{time:.3f} s", end="\r")

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

    if makeMovie:
        PyVista.Movie_simu(
            simu,
            "uy",
            folder,
            "uy.gif",
            deformFactor=1,
            plotMesh=True,
            N=400,
            nodeValues=True,
        )

    simu.Set_Iter(0)
    PyVista.Plot(simu, "uy", 1, plotMesh=True).show()

    Display.plt.show()