Elas2

Bending bracket component.

• 
• 
• 

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

Element type: TETRA10
Ne = 2088, Nn = 3537

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

Isotropic:
E = 2.10e+05, v = 0.3

solver : scipy

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

Unspecified.

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


W def = 4.04

Svm max = 10.53

Evm max = 0.01 %

Ux max = 2.12e-03
Ux min = -2.88e-03

Uy max = 5.15e-05
Uy min = -1.01e-02

Uz max = 2.47e-04
Uz min = -2.41e-04

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

Mesh : 108.415 ms
Boundary Conditions : 31.948 µs
Matrix : 474.501 ms
Solver : 625.634 ms
PostProcessing : 6.346 ms

import matplotlib.pyplot as plt
import numpy as np

from EasyFEA import Display, Models, ElemType, Simulations
from EasyFEA.Geoms import Point, Points

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

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

    # geom
    dim = 3
    L = 120  # mm
    h = L * 0.3

    # model
    E = 210000  # MPa (Young's modulus)
    v = 0.3  # Poisson's ratio
    coef = 1

    # load
    load = 800

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

    # Define points and crack geometry for the mesh
    pt1 = Point(isOpen=True, r=-10)
    pt2 = Point(x=L)
    pt3 = Point(x=L, y=h)
    pt4 = Point(x=h, y=h, r=10)
    pt5 = Point(x=h, y=L)
    pt6 = Point(y=L)
    pt7 = Point(x=h, y=h)

    contour = Points([pt1, pt2, pt3, pt4, pt5, pt6], h / 3)

    if dim == 2:
        mesh = contour.Mesh_2D([], ElemType.TRI3)
    elif dim == 3:
        mesh = contour.Mesh_Extrude([], [0, 0, -h], [4], elemType=ElemType.TETRA10)

    nodes_x0 = mesh.Nodes_Conditions(lambda x, y, z: x == 0)
    nodes_xL = mesh.Nodes_Conditions(lambda x, y, z: x == L)

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

    material = Models.Elastic.Isotropic(dim, E, v, planeStress=True, thickness=h)
    simu = Simulations.Elastic(mesh, material)

    simu.add_dirichlet(nodes_x0, [0] * dim, simu.Get_unknowns())
    simu.add_surfLoad(nodes_xL, [-800 / (h * h)], ["y"])

    sol = simu.Solve()
    simu.Save_Iter()

    # ----------------------------------------------
    # Results
    # ----------------------------------------------
    print(simu)

    Display.Plot_Mesh(simu, h / 2 / np.abs(sol).max())
    Display.Plot_BoundaryConditions(simu)
    Display.Plot_Result(simu, "Svm", nodeValues=True, coef=1 / coef, ncolors=20)

    plt.show()