}, 'stiffness matrix': { 'evaluation': 'closed form' } } } #--------------------------------------------------------------- # Mesh info print("Meshing...\n") L = 10000 n = 40 mesh = mesh_engine.Mesh() mesh.Elements = n mesh.Nodes = n + 1 mesh.NodesElement = 2 mesh.dofsNode = 1 mesh.d = 1 mesh.points = np.zeros(shape=(mesh.Nodes, 3)) mesh.points[:, 0] = np.linspace(0, L, mesh.Nodes) mesh.elements = element_table(mesh) mesh.elementMaterialTag = np.ones(n, dtype=int) mesh.elementType = np.array([["diff_bar"] for i in range(n)]) # Load info print("Applying BCs...\n")
""" Created on Sat Mar 21 17:38:10 2020 @author: Umberto Bettinardi [email protected] """ import numpy as np from FEM import FEM_engine as engine from FEM import mesh_engine as mesher from FEM import material_data as material from FEM import boundaries #------------------------------------------------------ mesh_info = mesher.Mesh() material_info = material.Material() load_info = boundaries.BoundaryConditions() #------------------------------------------------------ # Mesh info mesh_info.Elements = 4 mesh_info.Nodes = 4 mesh_info.NodesElement = 2 mesh_info.dofsNode = 1 mesh_info.elements = np.array([[0, 1], [1, 2], [1, 2], [2, 3]]) #------------------------------------------------------ # Material info material_info.k_spring = np.array([20, 30, 2, 10])