Skip to content
Snippets Groups Projects
Commit 924e1997 authored by Noels Ludovic's avatar Noels Ludovic
Browse files

Merge branch 'dev_mm' into 'main'

[Feature] - USF Cell added in Lattice_Cell_script

See merge request !5
parents 815bc5d7 aca42065
No related branches found
No related tags found
1 merge request!5[Feature] - USF Cell added in Lattice_Cell_script
......@@ -525,7 +525,33 @@ def Geo_generator(name, rho, vfr, length, width, height, sb, Radius, vfraction,
[[sb/2,sb/2,0],[sb/4,3*sb/4,sb/2],[r,0,0],[ri,0,0]], #9
[[sb/2,sb/2,sb],[3*sb/4,sb/4,sb/2],[r,0,0],[ri,0,0]], #14
[[sb/2,sb/2,sb],[3*sb/4,3*sb/4,sb/2],[r,0,0],[ri,0,0]]] #15
elif name== "10":
x_cell=x_cell+sb
y_cell=y_cell+sb
#z=0.5x and x=y
l=sb/2
z_cell=l
if vfr=="b":
r=Radius
A=[[[0,0,0],[0,sb,0],[r,0,0]], #2
[[sb,0,0],[sb,sb,0 ],[r,0,0]], #3
[[sb,sb,l],[0,0,0],[r,0,0]], #4
[[0,sb,0],[sb,0,l],[r,0,0]], #5
[[sb,sb/2,0],[sb/2,0,l/2],[r,0,0]], #6
[[sb/2,sb,l/2],[sb,sb/2,0],[r,0,0]], #7
[[sb/2,0,l/2],[sb/2,sb,l/2],[r,0,0]], #8
[[sb/2,0,l/2],[0,sb/2,l],[r,0,0]], #9
[[0,sb/2,l],[sb/2,sb,l/2],[r,0,0]], #10
[[0,0,l],[0,sb,l],[r,0,0]], #11
[[sb,0,l],[sb,sb,l],[r,0,0]]] #12
r0=r
......
......@@ -24,7 +24,7 @@ import matplotlib.pyplot as plt
from Auto_geometry_parametric import *
#############Can be changed##########################
cellType = "3" # from "2", "3", "4", "5", "6", "7", "8", "9", "91" (or "92" if typeRand ="b")
cellType = "3" # from "2", "3", "4", "5", "6", "7", "8", "9", "10", "91" (or "92" if typeRand ="b")
typeRand = "a" # "a" for random volume fraction and "b" for random radius
nx =2 #number of cell along x
ny =2 #number of cell along y
......@@ -47,8 +47,8 @@ bow_angle_max=85 #maximum bow angle of cell 8, should be lower than 90
#############Cannot be changed##############################
vf_tpye=["a", "b"]
GeoID = []
GeoID.append(["2", "3", "4", "5", "6", "7", "8", "9", "91"])
GeoID.append(["2", "3", "4", "5", "6", "7", "8", "9", "91", "92"])
GeoID.append(["2", "3", "4", "5", "6", "7", "8", "9", "91", "10"])
GeoID.append(["2", "3", "4", "5", "6", "7", "8", "9", "91", "92", "10"])
cell_num = np.array([nx,ny,nz])
####################################################
......
......@@ -50,6 +50,7 @@ Auxetic 8
Octahedron v1 9
Octahedron v2 91
Octahedron v3 (hollow) 92
Cirp USF 10
"""
print (initial_data)
......
0% Loading or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment