def setUp(self):
self.model = visco_flow.YaguchiGr91FlowRule()
self.pTrange = np.linspace(473.0, 873.0)
self.hist0 = np.zeros((14, ))
self.T = 500.0
def setUp(self):
self.vmodel = visco_flow.YaguchiGr91FlowRule()
E = 92000.0
nu = 0.3
mu = E / (2 * (1 + nu))
K = E / (3 * (1 - 2 * nu))
self.emodel = elasticity.IsotropicLinearElasticModel(
mu, "shear", K, "bulk")
self.model = general_flow.TVPFlowRule(self.emodel, self.vmodel)
self.T_n = 300.0
self.e_n = np.zeros((6, ))
self.t_n = 0.0
self.h_n = np.zeros((2 + 12, ))
def setUp(self):
self.model1 = parse.parse_xml("test/examples.xml", "test_yaguchi")
mu_poly = [-0.11834615, 115.5, 48807.69]
K_poly = [-0.256417, 250.25, 105750.0]
shear = interpolate.PolynomialInterpolate(mu_poly)
bulk = interpolate.PolynomialInterpolate(K_poly)
elastic = elasticity.IsotropicLinearElasticModel(
shear, "shear", bulk, "bulk")
vmodel = visco_flow.YaguchiGr91FlowRule()
flow = general_flow.TVPFlowRule(elastic, vmodel)
self.model2 = models.GeneralIntegrator(elastic, flow)
self.T = 500.0
self.tmax = 10.0
self.nsteps = 100.0
self.emax = np.array([0.1, 0, 0, 0, 0, 0])
def setUp(self):
vmodel = visco_flow.YaguchiGr91FlowRule()
E_poly = np.array([-3.077e-1, 3.003e2, 1.269e5])
nu = 0.3
mu_poly = list(E_poly / (2 * (1 + nu)))
K_poly = list(E_poly / (3 * (1 - 2 * nu)))
shear = interpolate.PolynomialInterpolate(mu_poly)
bulk = interpolate.PolynomialInterpolate(K_poly)
self.elastic = elasticity.IsotropicLinearElasticModel(
shear, "shear", bulk, "bulk")
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 = 550.0
self.nsteps = 100
plt.tight_layout()
plt.show()
# Isochronous curve
time = 50000.0 * 3600.0
res = drivers.isochronous_curve(model, time)
plt.figure()
plt.plot(res['strain'], res['stress'], 'k-')
plt.xlabel("Strain (mm/mm)")
plt.ylabel("Stress (MPa)")
plt.title("Isochronous stress-strain curve")
plt.tight_layout()
plt.show()
# Need a temperature sensitive model for this one!
vmodel = visco_flow.YaguchiGr91FlowRule()
E_poly = np.array([-3.077e-1, 3.003e2, 1.269e5])
nu = 0.3
mu_poly = list(E_poly / (2 * (1 + nu)))
K_poly = list(E_poly / (3 * (1 - 2 * nu)))
shear = interpolate.PolynomialInterpolate(mu_poly)
bulk = interpolate.PolynomialInterpolate(K_poly)
elastic = elasticity.IsotropicLinearElasticModel(shear, "shear", bulk,
"bulk")
flow = general_flow.TVPFlowRule(elastic, vmodel)
model = models.GeneralIntegrator(elastic, flow)