def setUp(self):
self.G1 = 36000.0
self.G2 = 25000.0
self.K1 = 10000.0
self.K2 = 12000.0
self.elastic1 = elasticity.IsotropicLinearElasticModel(
self.G1, "shear", self.K1, "bulk")
self.elastic2 = elasticity.IsotropicLinearElasticModel(
self.G2, "shear", self.K2, "bulk")
def setUp(self):
self.E = 100000.0
self.nu = 0.3
self.mu = self.E / (2 * (1 + self.nu))
self.K = self.E / (3 * (1 - 2.0 * self.nu))
self.model_Ev = elasticity.IsotropicLinearElasticModel(
self.E, "youngs", self.nu, "poissons")
self.model_GK = elasticity.IsotropicLinearElasticModel(
self.mu, "shear", self.K, "bulk")
self.T = 300.0
def setUp(self):
self.hist0 = np.zeros((13, ))
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.Kp = 1000.0
self.H = 1000.0
self.elastic = elasticity.IsotropicLinearElasticModel(
self.mu, "shear", self.K, "bulk")
surface = surfaces.IsoKinJ2()
iso = hardening.LinearIsotropicHardeningRule(self.s0, self.Kp)
kin = hardening.LinearKinematicHardeningRule(self.H)
hrule = hardening.CombinedHardeningRule(iso, kin)
flow = ri_flow.RateIndependentAssociativeFlow(surface, hrule)
self.model = models.SmallStrainRateIndependentPlasticity(
self.elastic, flow)
self.efinal = np.array([0.1, -0.05, 0.02, -0.03, 0.1, -0.15])
self.tfinal = 10.0
self.T = 300.0
self.nsteps = 10
def setUp(self):
self.E = 92000.0
self.nu = 0.3
self.s0 = 180.0
self.Kp = 1000.0
self.H = 1000.0
self.elastic = elasticity.IsotropicLinearElasticModel(self.E, "youngs",
self.nu, "poissons")
surface = surfaces.IsoKinJ2()
iso = hardening.LinearIsotropicHardeningRule(self.s0, self.Kp)
kin = hardening.LinearKinematicHardeningRule(self.H)
hrule = hardening.CombinedHardeningRule(iso, kin)
flow = ri_flow.RateIndependentAssociativeFlow(surface, hrule)
self.bmodel = models.SmallStrainRateIndependentPlasticity(self.elastic,
flow)
self.fn = interpolate.PolynomialInterpolate([-6.653e-9,2.952e-4,-6.197e-1])
self.C = 32.06
self.lmr = larsonmiller.LarsonMillerRelation(self.fn, self.C)
self.effective = damage.VonMisesEffectiveStress()
self.model = damage.LarsonMillerCreepDamageModel_sd(
self.elastic, self.lmr, self.effective, self.bmodel)
self.stress = np.array([100,-50.0,300.0,-99,50.0,125.0])
self.T = 100.0
self.s_np1 = self.stress
self.s_n = np.array([-25,150,250,-25,-100,25])
self.d_np1 = 0.5
self.d_n = 0.4
self.e_np1 = np.array([0.1,-0.01,0.15,-0.05,-0.1,0.15])
self.e_n = np.array([-0.05,0.025,-0.1,0.2,0.11,0.13])
self.T_np1 = self.T
self.T_n = 90.0
self.t_np1 = 1.0
self.t_n = 0.0
self.dt = self.t_np1 - self.t_n
self.u_n = 0.0
self.p_n = 0.0
# This is a rather boring baseline history state to probe, but I can't
# think of a better way to get a "generic" history from a generic model
self.hist_n = np.array([self.d_n] + list(self.bmodel.init_store()))
self.x_trial = np.array([50,-25,150,-150,190,100.0] + [0.41])
self.nsteps = 10
self.etarget = np.array([0.1,-0.025,0.02,0.015,-0.02,-0.05])
self.ttarget = 10.0
def setUp(self):
self.model1 = parse.parse_xml("test/examples.xml", "test_powerdamage")
E = 92000.0
nu = 0.3
s0 = 180.0
Kp = 1000.0
elastic = elasticity.IsotropicLinearElasticModel(
E, "youngs", nu, "poissons")
surface = surfaces.IsoJ2()
hrule = hardening.LinearIsotropicHardeningRule(s0, Kp)
flow = ri_flow.RateIndependentAssociativeFlow(surface, hrule)
bmodel = models.SmallStrainRateIndependentPlasticity(elastic, flow)
a = 2.2
A = 2e-5
self.model2 = damage.NEMLPowerLawDamagedModel_sd(elastic, A, a, bmodel)
self.T = 300.0
self.tmax = 10.0
self.nsteps = 100.0
self.emax = np.array([0.05, 0, 0, 0, 0, 0])
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.model1 = parse.parse_xml("test/examples.xml", "test_j2comb")
mu = 40000.0
K = 84000.0
ys = 100.0
r = 100.0
d = 1000.0
KH = 1000.0
elastic = elasticity.IsotropicLinearElasticModel(
mu, "shear", K, "bulk")
surface = surfaces.IsoKinJ2()
iso = hardening.VoceIsotropicHardeningRule(ys, r, d)
kin = hardening.LinearKinematicHardeningRule(KH)
hrule = hardening.CombinedHardeningRule(iso, kin)
flow = ri_flow.RateIndependentAssociativeFlow(surface, hrule)
self.model2 = models.SmallStrainRateIndependentPlasticity(
elastic, flow)
self.T = 300.0
self.tmax = 10.0
self.nsteps = 100.0
self.emax = np.array([0.1, 0, 0, 0, 0, 0])
def setUp(self):
self.model1 = parse.parse_xml("test/examples.xml", "test_j2isocomb")
E = 150000.0
nu = 0.3
ys = 100.0
H = 100.0
r = 100.0
d = 1000.0
elastic = elasticity.IsotropicLinearElasticModel(
E, "youngs", nu, "poissons")
surface = surfaces.IsoJ2()
hrule1 = hardening.LinearIsotropicHardeningRule(ys, H)
hrule2 = hardening.VoceIsotropicHardeningRule(0.0, r, d)
hrule = hardening.CombinedIsotropicHardeningRule([hrule1, hrule2])
flow = ri_flow.RateIndependentAssociativeFlow(surface, hrule)
self.model2 = models.SmallStrainRateIndependentPlasticity(
elastic, flow)
self.T = 300.0
self.tmax = 10.0
self.nsteps = 100.0
self.emax = np.array([0.1, 0, 0, 0, 0, 0])
def shear():
"""
Recreate some of the results of Bazant 2013
"""
F = lambda t: np.array([[1,t,0],[0,1.0,0],[0,0,1]])
tmax = 10.0
nsteps = 250
G = 100.0
K = 10.0
elastic = elasticity.IsotropicLinearElasticModel(G, "shear", K, "bulk")
model = models.SmallStrainElasticity(elastic)
res = drivers.def_grad_driver(model, F, tmax, nsteps)
theta = [Fi[0,1] / 2.0 for Fi in res['F']]
nstress = [Si[5] / G / np.sqrt(2) for Si in res['stress']]
plt.plot(theta, nstress, label = "Truesdell")
model = models.SmallStrainElasticity(elastic, truesdell = False)
res = drivers.def_grad_driver(model, F, tmax, nsteps)
theta = [Fi[0,1] / 2.0 for Fi in res['F']]
nstress = [Si[5] / G / np.sqrt(2) for Si in res['stress']]
plt.plot(theta, nstress, label = "Jaumann")
plt.xlabel(r"$\tan\left(\theta\right / 2)$")
plt.ylabel(r"$\tau_{12}/G$")
plt.legend(loc = 'best')
plt.show()
def rotation():
"""
Rotate a fixed stress tensor
"""
def F(t):
cut = 1.0
bF = lambda tt: np.array([[1+tt,0,0],[0,1.0-tt/4.0,0],[0,0,1-tt/4.0]])
R = lambda tt: np.array([
[np.cos(tt*(np.pi/2)), -np.sin(tt*(np.pi/2)),0],
[np.sin(tt*(np.pi/2)), np.cos(tt*(np.pi/2)),0],
[0, 0, 1]])
if t <= cut:
F = bF(t)
else:
F = np.dot(R(t-cut), bF(cut))
return F
tmax = 2.0
nsteps = 1000
E = 100.0
nu = 0.25
elastic = elasticity.IsotropicLinearElasticModel(E, "youngs", nu, "poissons")
model = models.SmallStrainElasticity(elastic)
res = drivers.def_grad_driver(model, F, tmax, 2*nsteps)
print(res['stress'][nsteps])
print(res['stress'][2*nsteps])
def simple_ex():
E = 92000.0
nu = 0.3
s0 = 180.0
Kp = 1000.0
H = 1000.0
elastic = elasticity.IsotropicLinearElasticModel(E, "youngs", nu,
"poissons")
surface = surfaces.IsoKinJ2()
iso = hardening.LinearIsotropicHardeningRule(s0, Kp)
kin = hardening.LinearKinematicHardeningRule(H)
hrule = hardening.CombinedHardeningRule(iso, kin)
flow = ri_flow.RateIndependentAssociativeFlow(surface, hrule)
bmodel = models.SmallStrainRateIndependentPlasticity(elastic, flow)
A = 0.0e-6
a = 2.2
model_off = damage.NEMLPowerLawDamagedModel_sd(elastic, A, a, bmodel)
A = 2e-5
model_on = damage.NEMLPowerLawDamagedModel_sd(elastic, A, a, bmodel)
res_off = drivers.uniaxial_test(model_off, 1.0e-2, emax=0.13)
res_on = drivers.uniaxial_test(model_on, 1.0e-2, emax=0.13)
plt.plot(res_off['strain'], res_off['stress'], 'k-')
plt.plot(res_on['strain'], res_on['stress'], 'r-')
plt.show()
def example4():
# T is in hours, strain in percent, stress in MPa
A = 1.85e-10
n = 2.5
m = 0.3
smodel = creep.PowerLawCreep(A, n)
cmodel = creep.J2CreepModel(smodel)
E = 150000.0
nu = 0.3
sY = 200.0
H = E / 25.0
elastic = elasticity.IsotropicLinearElasticModel(E, "youngs", nu,
"poissons")
surface = surfaces.IsoJ2()
iso = hardening.LinearIsotropicHardeningRule(sY, H)
flow = ri_flow.RateIndependentAssociativeFlow(surface, iso)
pmodel = models.SmallStrainRateIndependentPlasticity(elastic, flow)
model = models.SmallStrainCreepPlasticity(elastic, pmodel, cmodel)
res = drivers.creep(model, 205.0, 3600.0, 100.0, verbose = False,
nsteps_up = 500)
plt.plot(res['strain'], res['stress'])
plt.show()
def setUp(self):
string_rep = '<test_j2iso type="SmallStrainRateIndependentPlasticity"><elastic type="IsotropicLinearElasticModel"><m1>103561.64383561644</m1><m1_type>youngs</m1_type><m2>0.2945205479452055</m2><m2_type>poissons</m2_type></elastic><flow type="RateIndependentAssociativeFlow"><surface type="IsoJ2"/><hardening type="LinearIsotropicHardeningRule"><s0>100.0</s0><K>1000.0</K></hardening></flow></test_j2iso>'
self.model1 = parse.parse_string(string_rep)
mu = 40000.0
K = 84000.0
ys = 100.0
H = 1000.0
elastic = elasticity.IsotropicLinearElasticModel(
mu, "shear", K, "bulk")
surface = surfaces.IsoJ2()
hrule = hardening.LinearIsotropicHardeningRule(ys, H)
flow = ri_flow.RateIndependentAssociativeFlow(surface, hrule)
self.model2 = models.SmallStrainRateIndependentPlasticity(
elastic, flow)
self.T = 300.0
self.tmax = 10.0
self.nsteps = 100
self.emax = np.array([0.1, 0, 0, 0, 0, 0])
def setUp(self):
self.model1 = parse.parse_xml(localize("examples.xml"), "test_perfect")
E = [-100, 100000]
nu = 0.3
youngs = interpolate.PolynomialInterpolate(E)
poissons = interpolate.ConstantInterpolate(nu)
elastic = elasticity.IsotropicLinearElasticModel(
youngs, "youngs", poissons, "poissons")
surface = surfaces.IsoJ2()
Ts = [100.0, 300.0, 500.0, 700.0]
Sys = [1000.0, 120.0, 60.0, 30.0]
yields = interpolate.PiecewiseLinearInterpolate(Ts, Sys)
self.model2 = models.SmallStrainPerfectPlasticity(
elastic, surface, yields)
self.T = 550.0
self.tmax = 10.0
self.nsteps = 50
self.emax = np.array([0.1, 0.05, 0, -0.025, 0, 0])
def setUp(self):
self.ro = 5
self.ri = 4.5
self.h = 2.5
self.nr = 10
self.nt = 20
self.nz = 5
self.E = 150000.0
self.nu = 0.35
self.tube = receiver.Tube(self.ro, self.ro - self.ri, self.h, self.nr,
self.nt, self.nz)
self.times = np.array([0, 1])
self.tube.set_times(self.times)
self.tube.set_pressure_bc(
receiver.PressureBC(self.times, self.times * 0))
self.emodel = elasticity.IsotropicLinearElasticModel(
self.E, "youngs", self.nu, "poissons")
self.mat = models.SmallStrainElasticity(self.emodel)
self.d = 0.25
self.force_exact = (np.pi * (self.ro**2.0 - self.ri**2.0) * self.E *
self.d / self.h)
self.stiffness_exact = (np.pi * (self.ro**2.0 - self.ri**2.0) *
self.E / self.h)
self.solver = structural.PythonTubeSolver(verbose=False)
def plastic_model(E, Y, a, nu=0.3):
elmodel = elasticity.IsotropicLinearElasticModel(E, "youngs", nu,
"poissons")
surf = surfaces.IsoJ2()
return models.SmallStrainPerfectPlasticity(elmodel, surf, Y, alpha=a)
def verify_Q():
E = 30000.0
nu = 0.3
sy = 100.0
Q = 50.0
b = 100.0
C = 0.0
g = 0.0
mu = E / (2 * (1.0 + nu))
K = E / (3 * (1 - 2 * nu))
elastic = elasticity.IsotropicLinearElasticModel(mu, "shear", K, "bulk")
surface = surfaces.IsoKinJ2()
iso = hardening.VoceIsotropicHardeningRule(sy, Q, b)
gmodels = [hardening.ConstantGamma(g)]
hrule = hardening.Chaboche(iso, [C], gmodels, [0.0], [1.0])
flow = ri_flow.RateIndependentNonAssociativeHardening(surface, hrule)
model = models.SmallStrainRateIndependentPlasticity(elastic,
flow,
verbose=False)
res = drivers.uniaxial_test(model, 1.0e-2, emax=0.2)
stress = res['stress']
print("Q: %f / %f" % (Q, stress[-1] - sy))
def verify_Cg():
E = 30000.0
nu = 0.3
sy = 100.0
Q = 0.0
b = 0.0
C = 1000.0
g = 10.0
mu = E / (2 * (1.0 + nu))
K = E / (3 * (1 - 2 * nu))
elastic = elasticity.IsotropicLinearElasticModel(mu, "shear", K, "bulk")
surface = surfaces.IsoKinJ2()
iso = hardening.VoceIsotropicHardeningRule(sy, Q, b)
hrule = hardening.Chaboche(iso, [C], [hardening.ConstantGamma(g)], [0.0],
[1.0])
flow = ri_flow.RateIndependentNonAssociativeHardening(surface, hrule)
model = models.SmallStrainRateIndependentPlasticity(elastic,
flow,
verbose=False)
res = drivers.strain_cyclic(model, 0.4, -1.0, 1.0e-4, 1)
strain = res['strain']
stress = res['stress']
mv = np.max(np.abs(stress))
hu = mv - sy
print("C/y: %f / %f" % ((C / g), hu))
def setUp(self):
self.model1 = parse.parse_xml("test/examples.xml",
"test_creep_plasticity")
A = 1.85e-10
n = 2.5
smodel = creep.PowerLawCreep(A, n)
cmodel = creep.J2CreepModel(smodel)
E = 150000.0
nu = 0.3
sY = 200.0
H = E / 50.0
elastic = elasticity.IsotropicLinearElasticModel(
E, "youngs", nu, "poissons")
surface = surfaces.IsoJ2()
iso = hardening.LinearIsotropicHardeningRule(sY, H)
flow = ri_flow.RateIndependentAssociativeFlow(surface, iso)
pmodel = models.SmallStrainRateIndependentPlasticity(elastic, flow)
self.model2 = models.SmallStrainCreepPlasticity(
elastic, pmodel, cmodel)
self.T = 300.0
self.tmax = 10.0
self.nsteps = 100.0
self.emax = np.array([0.1, 0, 0, 0, 0, 0])
def verify_warp3d():
E = 30000.0
nu = 0.3
sy = 100.0
Q = 50.0
b = 100.0
C = 1000.0
g = 10.0
mu = E / (2 * (1.0 + nu))
K = E / (3 * (1 - 2 * nu))
elastic = elasticity.IsotropicLinearElasticModel(mu, "shear", K, "bulk")
surface = surfaces.IsoKinJ2()
iso = hardening.VoceIsotropicHardeningRule(sy, Q, b)
hrule = hardening.Chaboche(iso, [C], [hardening.ConstantGamma(g)], [0.0],
[1.0])
flow = ri_flow.RateIndependentNonAssociativeHardening(surface, hrule)
model = models.SmallStrainRateIndependentPlasticity(elastic,
flow,
verbose=False)
res = drivers.strain_cyclic(model, 0.0075, -1.0, 1.0e-4, 50)
strain = res['strain']
stress = res['stress']
data_warp = np.load('data_fa_warp.npy')
plt.plot(strain, stress, 'k-')
plt.plot(data_warp[0], data_warp[1], 'r-')
plt.show()
def setUp(self):
self.model1 = parse.parse_xml("test/examples.xml", "test_nonassri")
mu = 40000.0
K = 84000.0
ys = 100.0
r = 100.0
d = 1000.0
cs = [5.0, 10.0]
gs = [1000.0, 1000.0]
As = [0.0, 0.0]
ns = [1.0, 1.0]
elastic = elasticity.IsotropicLinearElasticModel(
mu, "shear", K, "bulk")
surface = surfaces.IsoKinJ2()
iso = hardening.VoceIsotropicHardeningRule(ys, r, d)
gammas = [hardening.ConstantGamma(g) for g in gs]
hmodel = hardening.Chaboche(iso, cs, gammas, As, ns)
flow = ri_flow.RateIndependentNonAssociativeHardening(surface, hmodel)
self.model2 = models.SmallStrainRateIndependentPlasticity(
elastic, flow)
self.T = 300.0
self.tmax = 10.0
self.nsteps = 100
self.emax = np.array([0.05, 0, 0, 0, 0, 0])
def setUp(self):
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.elastic = elasticity.IsotropicLinearElasticModel(
self.mu, "shear", self.K, "bulk")
surface = surfaces.IsoJ2()
hrule = hardening.VoceIsotropicHardeningRule(self.s0, self.R, self.d)
flow = ri_flow.RateIndependentAssociativeFlow(surface, hrule)
self.model = models.SmallStrainRateIndependentPlasticity(
self.elastic, flow, check_kt=False)
self.efinal = np.array([0.1, -0.05, 0.02, -0.03, 0.1, -0.15])
self.tfinal = 10.0
self.T = 300.0
self.nsteps = 10
def setUp(self):
self.model1 = parse.parse_xml("test/examples.xml", "test_pcreep")
E = [-100, 100000]
nu = 0.3
youngs = interpolate.PolynomialInterpolate(E)
poissons = interpolate.ConstantInterpolate(nu)
elastic = elasticity.IsotropicLinearElasticModel(
youngs, "youngs", poissons, "poissons")
surface = surfaces.IsoJ2()
Ts = [100.0, 300.0, 500.0, 700.0]
Sys = [1000.0, 120.0, 60.0, 30.0]
yields = interpolate.PiecewiseLinearInterpolate(Ts, Sys)
pmodel = models.SmallStrainPerfectPlasticity(elastic, surface, yields)
self.T = 550.0
self.tmax = 10.0
self.nsteps = 50.0
self.emax = np.array([0.1, 0.05, 0, -0.025, 0, 0])
A = 1.85e-10
n = 2.5
smodel = creep.PowerLawCreep(A, n)
cmodel = creep.J2CreepModel(smodel)
self.model2 = models.SmallStrainCreepPlasticity(
elastic, pmodel, cmodel)
def setUp(self):
self.hist0 = np.zeros((13, ))
self.A = 1.85e-10
self.n = 2.5
self.smodel = creep.PowerLawCreep(self.A, self.n)
self.cmodel = creep.J2CreepModel(self.smodel)
self.E = 150000.0
self.nu = 0.3
self.sY = 200.0
self.H = self.E / 50.0
self.elastic = elasticity.IsotropicLinearElasticModel(
self.E, "youngs", self.nu, "poissons")
self.surface = surfaces.IsoJ2()
self.iso = hardening.LinearIsotropicHardeningRule(self.sY, self.H)
self.flow = ri_flow.RateIndependentAssociativeFlow(
self.surface, self.iso)
self.pmodel = models.SmallStrainRateIndependentPlasticity(
self.elastic, self.flow)
self.model = models.SmallStrainCreepPlasticity(self.elastic,
self.pmodel,
self.cmodel)
self.efinal = np.array([0.1, -0.05, 0.02, -0.03, 0.1, -0.15])
self.tfinal = 10.0
self.T = 300.0
self.nsteps = 10
def complete(self):
self.E = 92000.0
self.nu = 0.3
self.s0 = 180.0
self.Kp = 1000.0
self.H = 1000.0
self.elastic = elasticity.IsotropicLinearElasticModel(self.E, "youngs",
self.nu, "poissons")
surface = surfaces.IsoKinJ2()
iso = hardening.LinearIsotropicHardeningRule(self.s0, self.Kp)
kin = hardening.LinearKinematicHardeningRule(self.H)
hrule = hardening.CombinedHardeningRule(iso, kin)
flow = ri_flow.RateIndependentAssociativeFlow(surface, hrule)
self.bmodel = models.SmallStrainRateIndependentPlasticity(self.elastic,
flow)
self.xi = 0.478
self.phi = 1.914
self.A = 10000000.0
self.model = damage.ModularCreepDamageModel_sd(
self.elastic,
self.A, self.xi, self.phi,
self.effective_model(),
self.bmodel)
self.stress = np.array([100,-50.0,300.0,-99,50.0,125.0])
self.T = 100.0
self.s_np1 = self.stress
self.s_n = np.array([-25,150,250,-25,-100,25])
self.d_np1 = 0.5
self.d_n = 0.4
self.e_np1 = np.array([0.1,-0.01,0.15,-0.05,-0.1,0.15])
self.e_n = np.array([-0.05,0.025,-0.1,0.2,0.11,0.13])
self.T_np1 = self.T
self.T_n = 90.0
self.t_np1 = 1.0
self.t_n = 0.0
self.u_n = 0.0
self.p_n = 0.0
# This is a rather boring baseline history state to probe, but I can't
# think of a better way to get a "generic" history from a generic model
self.hist_n = np.array([self.d_n] + list(self.bmodel.init_store()))
self.x_trial = np.array([50,-25,150,-150,190,100.0] + [0.41])
self.nsteps = 10
self.etarget = np.array([0.1,-0.025,0.02,0.015,-0.02,-0.05])
self.ttarget = 10.0
def setUp(self):
self.hist0 = np.zeros((6, ))
self.A = 1.85e-10
self.n = 2.5
self.smodel = creep.PowerLawCreep(self.A, self.n)
self.cmodel = creep.J2CreepModel(self.smodel)
self.E = 150000.0
self.nu = 0.3
self.sY = 200.0
self.elastic = elasticity.IsotropicLinearElasticModel(
self.E, "youngs", self.nu, "poissons")
self.surface = surfaces.IsoJ2()
self.pmodel = models.SmallStrainPerfectPlasticity(
self.elastic, self.surface, self.sY)
self.model = models.SmallStrainCreepPlasticity(self.elastic,
self.pmodel,
self.cmodel)
self.efinal = np.array([0.1, -0.05, 0.02, -0.05, 0.1, -0.15])
self.tfinal = 10.0
self.T = 300.0
self.nsteps = 10
def setUp(self):
self.eps0 = 1.0e2
self.D = 100.0
self.n = 5.2
self.s0 = 150.0
self.K = 10000.0
self.vmodel = walker.TestFlowRule(self.eps0, self.D, self.n, self.s0,
self.K)
E = 92000.0
nu = 0.3
self.emodel = elasticity.IsotropicLinearElasticModel(
E, "youngs", nu, "poissons")
self.lv = 1
self.model = walker.WalkerKremplSwitchRule(self.emodel, self.vmodel,
self.lv, self.eps0)
self.T_n = 300.0
self.e_n = np.zeros((6, ))
self.t_n = 0.0
self.h_n = np.zeros((2, ))
self.eps = 1.0e-6
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):
E = 200000.0
nu = 0.27
mu = E / (2 * (1.0 + nu))
K = E / (3 * (1 - 2 * nu))
s0 = 300.0
Kp = 0.0
c = [30000.0]
r = [60.0]
A = [0.0]
n = [1.0]
elastic = elasticity.IsotropicLinearElasticModel(
mu, "shear", K, "bulk")
surface = surfaces.IsoKinJ2()
iso = hardening.LinearIsotropicHardeningRule(s0, Kp)
gmodels = [hardening.ConstantGamma(g) for g in r]
hrule = hardening.Chaboche(iso, c, gmodels, A, n)
flow = ri_flow.RateIndependentNonAssociativeHardening(surface, hrule)
self.model = models.SmallStrainRateIndependentPlasticity(
elastic, flow, verbose=False, check_kt=False)
self.umodel = uniaxial.UniaxialModel(self.model)
self.de = 0.001
self.dt = 1.0
self.dT = 0.0
self.T0 = 0.0
self.nsteps = 100
def setUp(self):
self.mu = 29000.0
self.K = 64000.0
self.T = 325.0
self.model = elasticity.IsotropicLinearElasticModel(self.mu,
"shear", self.K, "bulk")
