def setUp(self):
self.n = 12.0
self.K = 150.0
self.k = 6.0
self.C = 24800.0
self.g0 = 300.0
self.Q = 86 - self.k
self.gs = 300.0
self.b = 10.0
self.beta = 0.0
self.surface = surfaces.IsoKinJ2()
self.iso = hardening.VoceIsotropicHardeningRule(self.k, self.Q, self.b)
cs = [self.C]
gs = [hardening.SatGamma(self.gs, self.g0, self.beta)]
As = [0.0]
ns = [1.0]
self.hardening = hardening.Chaboche(self.iso, cs, gs, As, ns)
self.fluidity = visco_flow.ConstantFluidity(self.K)
self.model = visco_flow.ChabocheFlowRule(self.surface, self.hardening,
self.fluidity, self.n)
self.hist0 = np.zeros((7, ))
self.T = 300.0
def setUp(self):
self.s0 = 200.0
self.K = 1000.0
self.n = 2
self.cs = list(np.array(range(self.n)) / (2 * self.n) + 10.0)
self.g0s = np.array(range(1, self.n + 1)) * 10.0
self.gss = 2.0 * self.g0s
self.betas = np.array(range(self.n)) / self.n
self.gammas = [
hardening.SatGamma(a, b, c)
for a, b, c in zip(self.g0s, self.gss, self.betas)
]
self.As = [0.0] * self.n
self.ns = [1.0] * self.n
self.iso = hardening.LinearIsotropicHardeningRule(self.s0, self.K)
self.model = hardening.Chaboche(self.iso, self.cs, self.gammas,
self.As, self.ns)
self.hist0 = np.zeros((1 + self.n * 6, ))
self.conform = 7
self.T = 300.0
def setUp(self):
self.g0 = 100.0
self.gs = 200.0
self.beta = 2.5
self.T = 300.0
self.model = hardening.SatGamma(self.gs, self.g0, self.beta)
n = 12.0
eta = 150.0
k = 6.0
C = 24800.0
g0 = 300.0
Q = 86 - k
gs = 300.0
b = 10.0
beta = 0.0
elastic = elasticity.IsotropicLinearElasticModel(mu, "shear", K, "bulk")
surface = surfaces.IsoKinJ2()
iso = hardening.VoceIsotropicHardeningRule(k, Q, b)
cs = [C]
gs = [hardening.SatGamma(gs, g0, beta)]
As = [0.0]
ns = [1.0]
hmodel = hardening.Chaboche(iso, cs, gs, As, ns)
fluidity = visco_flow.ConstantFluidity(eta)
vmodel = visco_flow.ChabocheFlowRule(
surface, hmodel, fluidity, n)
flow = general_flow.TVPFlowRule(elastic, vmodel)
model = models.GeneralIntegrator(elastic, flow, verbose = False)
# Uniaxial stress/strain curves at decades of strain rates
erates = np.logspace(-6,2,9)