def Run(self):
"""
Runs a simulation for a given couple (sy, n) and returns the (disp, force) couple.
"""
#MODEL DEFINITION
n = self.n
sy = self.sy
E = self.settings['E']
nu = self.settings['nu']
inner_radius = self.settings['inner_radius']
outer_radius = self.settings['outer_radius']
disp = self.settings['displacement'] / 2.
nFrames = self.settings['nFrames']
Nr = self.settings['Nr']
Nt = self.settings['Nt']
Na = self.settings['Na']
Ne = self.settings['Ne']
thickness = self.settings['thickness'] / 2.
unloading = settings['unloading']
export_fields = settings['export_fields']
export_fields = False
print E, nu, sy, n
material = Hollomon(labels="SAMPLE_MAT", E=E, nu=nu, sy=sy, n=n)
m = RingCompression(material=material,
inner_radius=inner_radius,
outer_radius=outer_radius,
disp=disp,
thickness=thickness,
nFrames=nFrames,
Nr=Nr,
Nt=Nt,
Na=Na,
workdir=workdir,
label=label,
elType=elType,
abqlauncher=abqlauncher,
cpus=cpus,
is_3D=is_3D,
unloading=unloading,
export_fields=export_fields)
# SIMULATION
m.MakeMesh()
m.MakeInp()
m.Run()
m.PostProc()
outputs = m.outputs
force = -4. * outputs['history']['force']
disp = -2. * outputs['history']['disp']
self.disp = disp
self.force = force
def Run(self):
"""
Runs a simulation for a given couple (sy, n) and returns the (disp, force) couple.
"""
#MODEL DEFINITION
sy = self.sy
n = self.n
E = self.settings['E']
nu = self.settings['nu']
inner_radius = self.settings['inner_radius']
outer_radius = self.settings['outer_radius']
disp = self.settings['displacement']/2.
nFrames = self.settings['nFrames']
Nr = self.settings['Nr']
Nt = self.settings['Nt']
Na = self.settings['Na']
Ne = self.settings['Ne']
thickness = self.settings['thickness']
print E, nu, sy, n
material = Hollomon(
labels = "SAMPLE_MAT",
E = E, nu = nu,
sy = sy, n = n)
m = RingCompression( material = material ,
inner_radius = inner_radius,
outer_radius = outer_radius,
disp = disp,
thickness = thickness,
nFrames = nFrames,
Nr = Nr,
Nt = Nt,
Na = Na,
workdir = workdir,
label = label,
elType = elType,
abqlauncher = abqlauncher,
cpus = cpus,
<<<<<<< HEAD
is_3D = True)
#TASKS run_sim = True plot = True #MODEL DEFINITION m = RingCompression( material = material , inner_radius = inner_radius, outer_radius = outer_radius, disp = disp/2, thickness = thickness/2, nFrames = nFrames, Nr = Nr, Nt = Nt, Na = Na, unloading = unloading, export_fields = export_fields, workdir = workdir, label = label, elType = elType, abqlauncher = abqlauncher, cpus = cpus, is_3D = is_3D, compart = compart) # SIMULATION m.MakeMesh() if run_sim: m.MakeInp() m.Run() m.PostProc()
#TASKS run_sim = True plot = True #MODEL DEFINITION m = RingCompression( material = material,inner_radius = inner_radius, outer_radius = outer_radius, disp = disp/2, thickness = thickness, nFrames = nFrames, Nr = Nr, Nt = Nt, Na = Na, workdir = "D:\donnees_pyth/workdir/", label = label, elType = elType, abqlauncher = abqlauncher, cpus = 1, compart = True, is_3D = True) # SIMULATION m.MakeMesh() if run_sim: m.MakeInp() m.Run() m.PostProc()
return -np.array(disp_exp), -np.array(force_exp)
disp_exp, force_exp = read_file(filename)
#MODEL DEFINITION
disp = displacement / 2
material = Hollomon(labels="SAMPLE_MAT", E=E, nu=nu, sy=sy, n=n)
m = RingCompression(material=material,
inner_radius=inner_radius,
outer_radius=outer_radius,
disp=disp,
thickness=thickness,
nFrames=nFrames,
Nr=Nr,
Nt=Nt,
Na=Na,
workdir=workdir,
label=label,
elType=elType,
abqlauncher=abqlauncher,
cpus=cpus,
is_3D=is_3D)
# SIMULATION
m.MakeMesh()
if run_sim:
m.MakeInp()
m.Run()
m.PostProc()
if node == 'epua-pd45':
abqlauncher = 'C:\SIMULIA/Abaqus/Commands/abaqus'
#TASKS
Run_simu = True
#MODEL DEFINITION
disp = displacement / 2.
model = RingCompression(inner_radius=inner_radius,
outer_radius=outer_radius,
disp=disp,
thickness=thickness,
nFrames=nFrames,
Nr=Nr,
Nt=Nt,
Na=Na,
workdir=workdir,
label=label,
elType=elType,
abqlauncher=abqlauncher,
cpus=cpus,
is_3D=is_3D,
compart=True)
# SIMULATION
model.MakeMesh()
mesh = model.mesh
centroids = mesh.centroids()
nodes = mesh.nodes
nodes_postion = np.array([nodes.x, nodes.y, nodes.z]).transpose()
conn = np.array(mesh.connectivity)
def Run(self):
"""
Runs a simulation for a given couple (sy, n) and returns the (disp, force) couple.
"""
#MODEL DEFINITION
sy_mean = self.sy_mean * np.ones(settings['Ne'])
n = self.n * np.ones(settings['Ne'])
Ssat = self.Ssat * np.ones(settings['Ne'])
E = self.settings['E']
nu = self.settings['nu']
inner_radius = self.settings['inner_radius']
outer_radius = self.settings['outer_radius']
disp = self.settings['displacement'] / 2.
nFrames = self.settings['nFrames']
Nr = self.settings['Nr']
Nt = self.settings['Nt']
Na = self.settings['Na']
Ne = self.settings['Ne']
thickness = self.settings['thickness'] / 2
#TASKS
run_sim = True
plot = True
print E[0], nu[0], Ssat[0], n[0], sy_mean[0]
#print E[0], nu[0], n[0], sy_mean[0]
ray_param = sy_mean / 1.253314
sy = np.random.rayleigh(ray_param, Ne)
labels = ['mat_{0}'.format(i + 1) for i in xrange(len(sy))]
material = [
materials.Bilinear(labels=labels[i],
E=E[i],
nu=nu[i],
Ssat=Ssat[i],
n=n[i],
sy=sy[i]) for i in xrange(Ne)
]
m = RingCompression(material=material,
inner_radius=inner_radius,
outer_radius=outer_radius,
disp=disp,
thickness=thickness,
nFrames=nFrames,
Nr=Nr,
Nt=Nt,
Na=Na,
unloading=unloading,
export_fields=export_fields,
workdir=workdir,
label=label,
elType=elType,
abqlauncher=abqlauncher,
cpus=cpus,
is_3D=is_3D,
compart=compart)
# SIMULATION
m.MakeMesh()
if run_sim:
m.MakeInp()
m.Run()
m.PostProc()
outputs = m.outputs
force = -4. * outputs['history']['force']
disp = -2 * outputs['history']['disp']
self.disp = disp
self.force = force
#TASKS Run_simu = True #MODEL DEFINITION disp = displacement/2. model = RingCompression( inner_radius = inner_radius, outer_radius = outer_radius, disp = disp, thickness = thickness, nFrames = nFrames, Nr = Nr, Nt = Nt, Na = Na, workdir = workdir, label = label, elType = elType, abqlauncher = abqlauncher, cpus = cpus, is_3D = is_3D, compart = True) # SIMULATION model.MakeMesh() mesh = model.mesh centroids = mesh.centroids() nodes = mesh.nodes nodes_postion = np.array([nodes.x, nodes.y, nodes.z]).transpose() conn = np.array(mesh.connectivity)