def setUp(self):
self.model1 = parse.parse_xml("test/examples.xml", "test_perzyna")
mu = 40000.0
K = 84000.0
elastic = elasticity.IsotropicLinearElasticModel(
mu, "shear", K, "bulk")
sy = 100.0
r = 100.0
d = 1000.0
KK = 1000.0
n = 5.0
eta = 500.0
surface = surfaces.IsoKinJ2()
iso = hardening.VoceIsotropicHardeningRule(sy, r, d)
kin = hardening.LinearKinematicHardeningRule(KK)
hmodel = hardening.CombinedHardeningRule(iso, kin)
gmodel = visco_flow.GPowerLaw(n, eta)
vmodel = visco_flow.PerzynaFlowRule(surface, hmodel, gmodel)
flow = general_flow.TVPFlowRule(elastic, vmodel)
self.model2 = models.GeneralIntegrator(elastic, flow)
self.T = 550.0 + 273.15
self.tmax = 10.0
self.nsteps = 100.0
self.emax = np.array([0.1, 0, 0, 0, 0, 0])
def setUp(self):
self.hist0 = np.zeros((7, ))
self.hist0 = np.zeros((7, ))
self.E = 92000.0
self.nu = 0.3
self.mu = self.E / (2 * (1 + self.nu))
self.K = self.E / (3 * (1 - 2 * self.nu))
self.s0 = 180.0
self.R = 150.0
self.d = 10.0
self.n = 2.0
self.eta = 200.0
self.elastic = elasticity.IsotropicLinearElasticModel(
self.mu, "shear", self.K, "bulk")
surface = surfaces.IsoJ2()
hrule = hardening.VoceIsotropicHardeningRule(self.s0, self.R, self.d)
g = visco_flow.GPowerLaw(self.n, self.eta)
vmodel = visco_flow.PerzynaFlowRule(surface, hrule, g)
flow = general_flow.TVPFlowRule(self.elastic, vmodel)
self.model = models.GeneralIntegrator(self.elastic, flow)
self.efinal = np.array([0.05, 0, 0, 0.02, 0, -0.01])
self.tfinal = 10.0
self.T = 300.0
self.nsteps = 100
def setUp(self):
self.s0 = 100.0
self.Kp = 1500.0
self.n = 2.0
self.eta = 100.0
self.surface = surfaces.IsoJ2()
self.hrule = hardening.LinearIsotropicHardeningRule(self.s0, self.Kp)
self.g = visco_flow.GPowerLaw(self.n, self.eta)
self.model = visco_flow.PerzynaFlowRule(self.surface, self.hrule,
self.g)
self.hist0 = np.zeros((1, ))
self.T = 300.0
def setUp(self):
self.s0 = 180.0
self.R = 150.0
self.d = 10.0
self.n = 5.0
self.eta = 20.0
self.surface = surfaces.IsoJ2()
self.hrule = hardening.VoceIsotropicHardeningRule(
self.s0, self.R, self.d)
self.g = visco_flow.GPowerLaw(self.n, self.eta)
self.model = visco_flow.PerzynaFlowRule(self.surface, self.hrule,
self.g)
self.hist0 = np.zeros((1, ))
self.T = 300.0
def test_rd(self):
s0 = 10.0
R = 150.0
d = 10.0
n = 3.0
eta = 200.0
surface = surfaces.IsoJ2()
hrule = hardening.VoceIsotropicHardeningRule(s0, R, d)
g = visco_flow.GPowerLaw(n, eta)
vmodel = visco_flow.PerzynaFlowRule(surface, hrule, g)
flow = general_flow.TVPFlowRule(self.elastic1, vmodel)
model = models.GeneralIntegrator(self.elastic1, flow)
self.very_close(model, self.emodel1)
model.set_elastic_model(self.elastic2)
self.very_close(model, self.emodel2)
H = 0.0001
elastic = elasticity.IsotropicLinearElasticModel(E, "youngs", nu,
"poissons")
surface = surfaces.IsoJ2()
iso = hardening.LinearIsotropicHardeningRule(sY, H)
flow = ri_flow.RateIndependentAssociativeFlow(surface, iso)
model1 = models.SmallStrainRateIndependentPlasticity(elastic, flow)
model2 = models.SmallStrainPerfectPlasticity(elastic, surface, sY)
n = 45.0
eta = 200.0
iso2 = hardening.LinearIsotropicHardeningRule(0, H)
gpower = visco_flow.GPowerLaw(n, eta)
vflow = visco_flow.PerzynaFlowRule(surface, iso2, gpower)
gflow = general_flow.TVPFlowRule(elastic, vflow)
model3 = models.GeneralIntegrator(elastic, gflow)
# T is in hours, strain in percent, stress in MPa
A = 1.85e-9
n = 2.5
m = 0.3
smodel = creep.NortonBaileyCreep(A, n, m)
cmodel = creep.J2CreepModel(smodel)
model4 = models.SmallStrainCreepPlasticity(elastic, model1, cmodel)
erate = 1.0e-4
emax = 0.025
def setUp(self):
self.n = 11.0
self.eta = 200.0
self.model = visco_flow.GPowerLaw(self.n, self.eta)
self.T = 300.0
