import matplotlib.pyplot as plt import matplotlib from mpl_toolkits.mplot3d import Axes3D import numpy as np import cuboct import pfea import frame3dd import subprocess import json from math import * import DM23dd import pickle mat_matrix = DM23dd.mat_matrix_from_json('data/bridge.json') subdiv = 1 #Subdivide the material matrix dims = np.shape(mat_matrix) new_mat_matrix = np.zeros( ((dims[0] - 2) * subdiv + 2, (dims[1] - 2) * subdiv + 2, (dims[2] - 2) * subdiv + 2)) for i in range(1, dims[0] - 1): for j in range(1, dims[1] - 1): for k in range(1, dims[2] - 1): if mat_matrix[i][j][k] == 1: dex = ((i - 1) * subdiv + 1, (j - 1) * subdiv + 1, (k - 1) * subdiv + 1) for l in range(0, subdiv): for m in range(0, subdiv): for n in range(0, subdiv):
import matplotlib.pyplot as plt import matplotlib from mpl_toolkits.mplot3d import Axes3D import numpy as np import cuboct import pfea import frame3dd import subprocess import json from math import * import DM23dd import pickle mat_matrix = DM23dd.mat_matrix_from_json('data/hut.json') subdiv = 1 #Subdivide the material matrix dims = np.shape(mat_matrix) new_mat_matrix = np.zeros(((dims[0]-2)*subdiv+2,(dims[1]-2)*subdiv+2,(dims[2]-2)*subdiv+2)) for i in range(1,dims[0]-1): for j in range(1,dims[1]-1): for k in range(1,dims[2]-1): if mat_matrix[i][j][k] == 1: dex = ((i-1)*subdiv+1,(j-1)*subdiv+1,(k-1)*subdiv+1) for l in range(0,subdiv): for m in range(0,subdiv): for n in range(0,subdiv): new_mat_matrix[l+dex[0]][m+dex[1]][n+dex[2]] = 1 #print(new_mat_matrix)
import matplotlib.pyplot as plt import matplotlib from mpl_toolkits.mplot3d import Axes3D import numpy as np import cuboct import pfea import frame3dd import subprocess import json from math import * import DM23dd import pickle mat_matrix = DM23dd.mat_matrix_from_json('data/bridge.json') subdiv = 1 #Subdivide the material matrix dims = np.shape(mat_matrix) new_mat_matrix = np.zeros(((dims[0]-2)*subdiv+2,(dims[1]-2)*subdiv+2,(dims[2]-2)*subdiv+2)) for i in range(1,dims[0]-1): for j in range(1,dims[1]-1): for k in range(1,dims[2]-1): if mat_matrix[i][j][k] == 1: dex = ((i-1)*subdiv+1,(j-1)*subdiv+1,(k-1)*subdiv+1) for l in range(0,subdiv): for m in range(0,subdiv): for n in range(0,subdiv): new_mat_matrix[l+dex[0]][m+dex[1]][n+dex[2]] = 1 #print(new_mat_matrix)