def setUp(self): self.L = crystallography.CubicLattice(1.0) self.L.add_slip_system([1, 1, 0], [1, 1, 1]) self.Q = rotations.Orientation(35.0, 17.0, 14.0, angle_type="degrees") self.S = tensors.Symmetric( np.array([[100.0, -25.0, 10.0], [-25.0, -17.0, 15.0], [10.0, 15.0, 35.0]])) self.nslip = self.L.ntotal self.static = 20.0 self.H = history.History() self.T = 300.0 self.s0 = [self.static] * self.nslip self.model = slipharden.FixedStrengthHardening(self.s0) self.g0 = 1.0 self.n = 3.0 self.sliprule = sliprules.PowerLawSlipRule(self.model, self.g0, self.n) self.fixed = history.History()
def setUp(self): self.L = crystallography.CubicLattice(1.0) self.L.add_slip_system([1, 1, 0], [1, 1, 1]) self.Q = rotations.Orientation(35.0, 17.0, 14.0, angle_type="degrees") self.S = tensors.Symmetric( np.array([[100.0, -25.0, 10.0], [-25.0, -17.0, 15.0], [10.0, 15.0, 35.0]])) self.nslip = self.L.ntotal self.current = -35.0 self.H = history.History() for i in range(self.nslip): self.H.add_scalar("strength" + str(i)) self.H.set_scalar("strength" + str(i), self.current) self.T = 300.0 self.s = np.array([self.current] * self.nslip) self.k = 1000.0 self.sat = 40.0 self.model = slipharden.FASlipHardening([self.k] * self.nslip, [self.sat] * self.nslip) self.g0 = 1.0 self.n = 3.0 self.sliprule = sliprules.PowerLawSlipRule(self.model, self.g0, self.n) self.fixed = history.History()
def setUp(self): self.strength = 35.0 self.H = history.History() self.H.add_scalar("strength") self.H.set_scalar("strength", self.strength) self.H.add_scalar("whatever") self.H.set_scalar("whatever", 0.5) self.tau0 = 10.0 self.tau_sat = 50.0 self.b = 2.5 self.strengthmodel = slipharden.VoceSlipHardening( self.tau_sat, self.b, self.tau0) self.g0 = 1.0 self.n = 3.0 self.slipmodel = sliprules.PowerLawSlipRule(self.strengthmodel, self.g0, self.n) self.imodel = inelasticity.AsaroInelasticity(self.slipmodel) self.L = crystallography.CubicLattice(1.0) self.L.add_slip_system([1, 1, 0], [1, 1, 1]) self.Q = rotations.Orientation(35.0, 17.0, 14.0, angle_type="degrees") self.S = tensors.Symmetric( np.array([[100.0, -25.0, 10.0], [-25.0, -17.0, 15.0], [10.0, 15.0, 35.0]])) self.T = 300.0 self.mu = 29000.0 self.E = 120000.0 self.nu = 0.3 self.emodel = elasticity.CubicLinearElasticModel( self.E, self.nu, self.mu, "moduli") self.dn = np.array([[4.1, 2.8, -1.2], [3.1, 7.1, 0.2], [4, 2, 3]]) self.dn = 0.5 * (self.dn + self.dn.T) self.d = tensors.Symmetric(self.dn) self.wn = np.array([[-9.36416517, 2.95527444, 8.70983194], [-1.54693052, 8.7905658, -5.10895168], [-8.52740468, -0.7741642, 2.89544992]]) self.wn = 0.5 * (self.wn - self.wn.T) self.w = tensors.Skew(self.wn) self.dmodel = crystaldamage.NilDamageModel() self.model = kinematics.DamagedStandardKinematicModel( self.emodel, self.imodel, self.dmodel) self.fspin = self.model.spin(self.S, self.d, self.w, self.Q, self.H, self.L, self.T, history.History()) self.fixed = self.model.decouple(self.S, self.d, self.w, self.Q, self.H, self.L, self.T, history.History())
def setUp(self): self.L = crystallography.CubicLattice(1.0) self.L.add_slip_system([1, 1, 0], [1, 1, 1]) self.Q = rotations.Orientation(35.0, 17.0, 14.0, angle_type="degrees") self.S = tensors.Symmetric( np.array([[100.0, -25.0, 10.0], [-25.0, -17.0, 15.0], [10.0, 15.0, 35.0]])) self.strength = 35.0 self.H = history.History() self.H.add_scalar("strength") self.H.set_scalar("strength", self.strength) self.T = 300.0 self.tau0 = 10.0 self.tau_sat = 50.0 self.b = 2.5 self.strengthmodel = slipharden.VoceSlipHardening( self.tau_sat, self.b, self.tau0) self.strengths = [self.strengthmodel] self.static = self.tau0 self.strength_values = [self.strength + self.static] self.g0 = 1.0 self.n = 3.0 self.model = sliprules.PowerLawSlipRule(self.strengthmodel, self.g0, self.n) self.tau = 33.0 self.fixed = history.History()
def setUp(self): self.L = crystallography.CubicLattice(1.0) self.L.add_slip_system([1, 1, 0], [1, 1, 1]) self.Q = rotations.Orientation(35.0, 17.0, 14.0, angle_type="degrees") self.S = tensors.Symmetric( np.array([[100.0, -25.0, 10.0], [-25.0, -17.0, 15.0], [10.0, 15.0, 35.0]])) self.strength0 = 15.0 self.strength1 = 20.0 self.strength2 = 20.0 self.H = history.History() self.H.add_scalar("strength_#0") self.H.add_scalar("strength_#1") self.H.add_scalar("strength_#2") self.H.set_scalar("strength_#0", self.strength0) self.H.set_scalar("strength_#1", self.strength1) self.H.set_scalar("strength_#2", self.strength2) self.T = 300.0 self.tau0_1 = 10.0 self.tau_sat_1 = 50.0 self.b_1 = 2.5 self.tau0_2 = 15.0 self.tau_sat_2 = 55.0 self.b_2 = 3.0 self.tau0_3 = 10.0 self.tau_sat_3 = 45.0 self.b_3 = 5.0 self.backstrength = slipharden.VoceSlipHardening( self.tau_sat_1, self.b_1, self.tau0_1) self.isostrength = slipharden.VoceSlipHardening( self.tau_sat_2, self.b_2, self.tau0_2) self.flowresistance = slipharden.VoceSlipHardening( self.tau_sat_3, self.b_3, self.tau0_3) self.strengths = [ self.backstrength, self.isostrength, self.flowresistance ] self.strength_values = [ self.strength0 + self.tau0_1, self.strength1 + self.tau0_2, self.strength2 + self.tau0_3 ] self.g0 = 1.0 self.n = 3.0 self.model = sliprules.KinematicPowerLawSlipRule( self.backstrength, self.isostrength, self.flowresistance, self.g0, self.n) self.tau = 80.0 self.fixed = history.History()
def setUp(self): self.a = 1.3 self.lattice = crystallography.CubicLattice(self.a) self.S = np.array([[20.0,-15.0,12.0],[-15.0,-40.0,5.0],[12.0,5.0,60.0]]) self.ST = tensors.Symmetric(self.S) self.Q = rotations.Orientation(30.0,43.0,10.0, angle_type = "degrees") self.QM = self.Q.to_matrix() self.lattice.add_slip_system([1,1,0],[1,1,1])
def setUp(self): self.L = crystallography.CubicLattice(1.0) self.L.add_slip_system([1, 1, 0], [1, 1, 1]) self.nslip = self.L.ntotal self.Q = rotations.Orientation(35.0, 17.0, 14.0, angle_type="degrees") self.S = tensors.Symmetric( np.array([[100.0, -25.0, 10.0], [-25.0, -17.0, 15.0], [10.0, 15.0, 35.0]])) * 2 self.static = 20.0 self.s0 = [self.static] * self.nslip self.k = 1000.0 self.sat = 40.0 self.m = 1.5 self.hmodel = slipharden.VocePerSystemHardening( self.s0, [self.k] * self.nslip, [self.sat] * self.nslip, [self.m] * self.nslip) self.g0 = 1.0 self.n = 3.0 self.sliprule = sliprules.PowerLawSlipRule(self.hmodel, self.g0, self.n) self.T = 300.0 self.c = 10.0 self.beta = 2.0 self.dmodel = crystaldamage.WorkPlaneDamage() self.nfunc = crystaldamage.SigmoidTransformation(self.c, self.beta) self.sfunc = crystaldamage.SigmoidTransformation(self.c, self.beta) self.model = crystaldamage.PlanarDamageModel(self.dmodel, self.sfunc, self.nfunc, self.L) self.huse = history.History() self.hmodel.populate_history(self.huse) self.model.populate_history(self.huse) for i in range(12): self.huse.set_scalar("strength" + str(i), 2.0) for j in range(4): self.huse.set_scalar("slip_damage_" + str(j), self.c * 0.4) self.hbase = self.huse.subset(["strength" + str(i) for i in range(12)]) self.hdmg = self.huse.subset( ["slip_damage_" + str(i) for i in range(4)]) self.fixed = history.History()
def setUp(self): self.L = crystallography.CubicLattice(1.0) self.L.add_slip_system([1, 1, 0], [1, 1, 1]) self.Q = rotations.Orientation(35.0, 17.0, 14.0, angle_type="degrees") self.S = tensors.Symmetric( np.array([[100.0, -25.0, 10.0], [-25.0, -17.0, 15.0], [10.0, 15.0, 35.0]])) * 3 self.H = history.History() self.T = 300.0 E = 160000.0 nu = 0.31 self.g0 = 1.0e-4 self.n = 10.0 K = E / 50.0 s0 = 75.0 M = matrix.SquareMatrix(self.L.ntotal, type="dense", data=[K] * (self.L.ntotal * self.L.ntotal)) self.isostrength = slipharden.GeneralLinearHardening(M, [s0 / 2] * self.L.ntotal, absval=False) self.flowresistance = slipharden.GeneralLinearHardening(M, [s0 / 2] * self.L.ntotal, absval=False) self.backstrength = slipharden.GeneralLinearHardening(M, [0] * self.L.ntotal, absval=False) self.strengths = [ self.backstrength, self.isostrength, self.flowresistance ] self.model = sliprules.KinematicPowerLawSlipRule( self.backstrength, self.isostrength, self.flowresistance, self.g0, self.n) self.model.populate_history(self.H) self.model.init_history(self.H) self.strength_values = np.linspace(0, 10, 36) + 5.0 self.H.copy_data(self.strength_values) self.fixed = history.History()
def setUp(self): self.model = inelasticity.NoInelasticity() self.L = crystallography.CubicLattice(1.0) self.L.add_slip_system([1, 1, 0], [1, 1, 1]) self.Q = rotations.Orientation(35.0, 17.0, 14.0, angle_type="degrees") self.S = tensors.Symmetric( np.array([[100.0, -25.0, 10.0], [-25.0, -17.0, 15.0], [10.0, 15.0, 35.0]])) self.T = 300.0 self.H = history.History() self.fixed = history.History()
def setUp(self): self.L = crystallography.CubicLattice(1.0) self.L.add_slip_system([1, 1, 0], [1, 1, 1]) self.Q = rotations.Orientation(35.0, 17.0, 14.0, angle_type="degrees") self.S = tensors.Symmetric( np.array([[100.0, -25.0, 10.0], [-25.0, -17.0, 15.0], [10.0, 15.0, 35.0]])) self.strength_0 = 35.0 self.strength_1 = -5.0 self.H = history.History() self.H.add_scalar("strength0") self.H.set_scalar("strength0", self.strength_0) self.H.add_scalar("strength1") self.H.set_scalar("strength1", self.strength_1) self.T = 300.0 self.tau0_1 = 10.0 self.tau_sat_1 = 50.0 self.b_1 = 2.5 self.tau0_2 = 5.0 self.tau_sat_2 = -25.0 self.b_2 = 1.5 self.static = self.tau0_1 + self.tau0_2 self.strength = self.strength_0 + self.strength_1 self.model1 = slipharden.VoceSlipHardening(self.tau_sat_1, self.b_1, self.tau0_1) self.model2 = slipharden.VoceSlipHardening(self.tau_sat_2, self.b_2, self.tau0_2) self.model = slipharden.SumSlipSingleStrengthHardening( [self.model1, self.model2]) self.g0 = 1.0 self.n = 3.0 self.sliprule = sliprules.PowerLawSlipRule(self.model, self.g0, self.n) self.fixed = history.History() self.nye = tensors.RankTwo([[1.1, 1.2, 1.3], [2.1, 2.2, 2.3], [3.1, 3.2, 3.3]]) self.nye_part = 0.0
def setUp(self): self.L = crystallography.CubicLattice(1.0) self.L.add_slip_system([1, 1, 0], [1, 1, 1]) self.nslip = self.L.ntotal self.Q = rotations.Orientation(35.0, 17.0, 14.0, angle_type="degrees") self.S = tensors.Symmetric( np.array([[100.0, -25.0, 10.0], [-25.0, -17.0, 15.0], [10.0, 15.0, 35.0]])) self.static = 20.0 self.s0 = [self.static] * self.nslip self.k = 1000.0 self.sat = 40.0 self.m = 1.5 self.hmodel = slipharden.VocePerSystemHardening( self.s0, [self.k] * self.nslip, [self.sat] * self.nslip, [self.m] * self.nslip) self.g0 = 1.0 self.n = 3.0 self.sliprule = sliprules.PowerLawSlipRule(self.hmodel, self.g0, self.n) self.T = 300.0 self.model = crystaldamage.NilDamageModel() self.huse = history.History() self.hmodel.populate_history(self.huse) self.model.populate_history(self.huse) for i in range(12): self.huse.set_scalar("strength" + str(i), 25.0) self.huse.set_scalar("whatever", 0.5) self.hbase = self.huse.subset(["strength" + str(i) for i in range(12)]) self.hdmg = self.huse.subset(["whatever"]) self.fixed = history.History()
def setUp(self): self.strength_0 = 35.0 self.H = history.History() self.H.add_scalar("strength0") self.H.set_scalar("strength0", self.strength_0) self.strength_1 = 25.0 self.H.add_scalar("strength1") self.H.set_scalar("strength1", self.strength_1) self.tau0_0 = 10.0 self.tau_sat_0 = 50.0 self.b_0 = 2.5 self.tau0_1 = 5.0 self.tau_sat_1 = 25.0 self.b_1 = 1.0 self.strengthmodel = slipharden.SumSlipSingleStrengthHardening([ slipharden.VoceSlipHardening(self.tau_sat_0, self.b_0, self.tau0_0), slipharden.VoceSlipHardening(self.tau_sat_1, self.b_1, self.tau0_1) ]) self.g0 = 1.0 self.n = 3.0 self.slipmodel = sliprules.PowerLawSlipRule(self.strengthmodel, self.g0, self.n) self.model = inelasticity.AsaroInelasticity(self.slipmodel) self.L = crystallography.CubicLattice(1.0) self.L.add_slip_system([1, 1, 0], [1, 1, 1]) self.Q = rotations.Orientation(35.0, 17.0, 14.0, angle_type="degrees") self.S = tensors.Symmetric( np.array([[100.0, -25.0, 10.0], [-25.0, -17.0, 15.0], [10.0, 15.0, 35.0]])) self.T = 300.0 self.fixed = history.History()
def setUp(self): self.vname = "strength" self.L = crystallography.CubicLattice(1.0) self.L.add_slip_system([1, 1, 0], [1, 1, 1]) self.Q = rotations.Orientation(35.0, 17.0, 14.0, angle_type="degrees") self.S = tensors.Symmetric( np.array([[100.0, -25.0, 10.0], [-25.0, -17.0, 15.0], [10.0, 15.0, 35.0]])) self.strength = 35.0 self.H = history.History() self.H.add_scalar("strength") self.H.set_scalar("strength", self.strength) self.T = 300.0 self.tau0 = 10.0 self.tau_sat = 50.0 self.b = 2.5 self.static = self.tau0 self.k = 1.2 self.model = slipharden.VoceSlipHardening(self.tau_sat, self.b, self.tau0, k=self.k) self.g0 = 1.0 self.n = 3.0 self.sliprule = sliprules.PowerLawSlipRule(self.model, self.g0, self.n) self.nye = tensors.RankTwo([[1.1, 1.2, 1.3], [2.1, 2.2, 2.3], [3.1, 3.2, 3.3]]) self.nye_part = self.k * np.sqrt( la.norm(self.nye.data.reshape((3, 3)), ord='fro')) self.fixed = history.History() self.fixed.add_ranktwo("nye") self.fixed.set_ranktwo("nye", self.nye)
def setUp(self): self.A = 1.0e-5 self.n = 3.1 self.model1 = inelasticity.PowerLawInelasticity(self.A, self.n) self.L = crystallography.CubicLattice(1.0) self.L.add_slip_system([1, 1, 0], [1, 1, 1]) self.Q = rotations.Orientation(35.0, 17.0, 14.0, angle_type="degrees") self.S = tensors.Symmetric( np.array([[100.0, -25.0, 10.0], [-25.0, -17.0, 15.0], [10.0, 15.0, 35.0]])) self.T = 300.0 self.strength = 35.0 self.H = history.History() self.H.add_scalar("strength") self.H.set_scalar("strength", self.strength) self.tau0 = 10.0 self.tau_sat = 50.0 self.b = 2.5 self.strengthmodel = slipharden.VoceSlipHardening( self.tau_sat, self.b, self.tau0) self.g0 = 1.0 self.n = 3.0 self.slipmodel = sliprules.PowerLawSlipRule(self.strengthmodel, self.g0, self.n) self.model2 = inelasticity.AsaroInelasticity(self.slipmodel) self.model = inelasticity.CombinedInelasticity( [self.model1, self.model2]) self.fixed = history.History()
def setUp(self): self.N = 10 self.tau0 = 10.0 self.tau_sat = 50.0 self.b = 2.5 self.strengthmodel = slipharden.VoceSlipHardening(self.tau_sat, self.b, self.tau0) self.g0 = 1.0 self.n = 3.0 self.slipmodel = sliprules.PowerLawSlipRule(self.strengthmodel, self.g0, self.n) self.imodel = inelasticity.AsaroInelasticity(self.slipmodel) self.L = crystallography.CubicLattice(1.0) self.L.add_slip_system([1,1,0],[1,1,1]) self.Q = rotations.Orientation(35.0,17.0,14.0, angle_type = "degrees") self.mu = 29000.0 self.E = 120000.0 self.nu = 0.3 self.emodel = elasticity.CubicLinearElasticModel(self.E, self.nu, self.mu, "moduli") self.kmodel = kinematics.StandardKinematicModel(self.emodel, self.imodel) self.model = singlecrystal.SingleCrystalModel(self.kmodel, self.L, miter = 120) self.orientations = rotations.random_orientations(self.N) self.D = np.array([0.01,-0.002,-0.003,0.012,-0.04,0.01]) self.W = np.array([0.02,-0.02,0.03]) self.T = 300.0 self.dt = 2.0
def setUp(self): self.vname = "wee" self.L = crystallography.CubicLattice(1.0) self.L.add_slip_system([1, 1, 0], [1, 1, 1]) self.Q = rotations.Orientation(35.0, 17.0, 14.0, angle_type="degrees") self.S = tensors.Symmetric( np.array([[100.0, -25.0, 10.0], [-25.0, -17.0, 15.0], [10.0, 15.0, 35.0]])) self.strength = 35.0 self.H = history.History() self.H.add_scalar(self.vname) self.H.set_scalar(self.vname, self.strength) self.T = 300.0 self.tau0 = 10.0 self.tau_sat = 50.0 self.b = 2.5 self.static = self.tau0 self.model = slipharden.VoceSlipHardening(self.tau_sat, self.b, self.tau0) self.model.set_varnames([self.vname]) self.g0 = 1.0 self.n = 3.0 self.sliprule = sliprules.PowerLawSlipRule(self.model, self.g0, self.n) self.fixed = history.History() self.nye = tensors.RankTwo([[1.1, 1.2, 1.3], [2.1, 2.2, 2.3], [3.1, 3.2, 3.3]]) self.nye_part = 0.0
def setUp(self): self.L = crystallography.CubicLattice(1.0) self.L.add_slip_system([1, 1, 0], [1, 1, 1]) self.Q = rotations.Orientation(35.0, 17.0, 14.0, angle_type="degrees") self.S = tensors.Symmetric( np.array([[100.0, -25.0, 10.0], [-25.0, -17.0, 15.0], [10.0, 15.0, 35.0]])) self.nslip = self.L.ntotal self.static = 20.0 self.current = 25.0 self.H = history.History() for i in range(self.nslip): self.H.add_scalar("strength" + str(i)) self.H.set_scalar("strength" + str(i), self.current) self.T = 300.0 self.M = matrix.SquareMatrix(self.nslip, type="block", data=[0.1, 0.2, 0.3, 0.4], blocks=[6, 6]) self.s0 = [self.static] * self.nslip self.model = slipharden.GeneralLinearHardening(self.M, self.s0, absval=False) self.g0 = 1.0 self.n = 3.0 self.sliprule = sliprules.PowerLawSlipRule(self.model, self.g0, self.n) self.fixed = history.History()
def setUp(self): self.tau0 = 10.0 self.tau_sat = 50.0 self.b = 2.5 self.strengthmodel = slipharden.VoceSlipHardening( self.tau_sat, self.b, self.tau0) self.g0 = 1.0 self.n = 3.0 self.slipmodel = sliprules.PowerLawSlipRule(self.strengthmodel, self.g0, self.n) self.imodel = inelasticity.AsaroInelasticity(self.slipmodel) self.L = crystallography.CubicLattice(1.0) self.L.add_slip_system([1, 1, 0], [1, 1, 1]) self.Q = rotations.Orientation(35.0, 17.0, 14.0, angle_type="degrees") self.mu = 29000.0 self.E = 120000.0 self.nu = 0.3 self.emodel = elasticity.CubicLinearElasticModel( self.E, self.nu, self.mu, "moduli") self.kmodel = kinematics.StandardKinematicModel( self.emodel, self.imodel) self.model = singlecrystal.SingleCrystalModel(self.kmodel, self.L, initial_rotation=self.Q) self.model_no_rot = singlecrystal.SingleCrystalModel( self.kmodel, self.L, initial_rotation=self.Q, update_rotation=False, verbose=False) self.T = 300.0 self.stress_n = np.array([120.0, -60.0, 170.0, 35.0, 80.0, -90.0]) self.stress_np1 = np.array([15.0, -40.0, 120.0, 70.0, -10.0, -50.0]) self.d = np.array([0.1, -0.2, 0.25, 0.11, -0.05, 0.075]) self.w = np.array([0.1, 0.2, -0.2]) self.strength_n = 25.0 self.strength_np1 = 30.0 self.S_np1 = tensors.Symmetric(common.usym(self.stress_np1)) self.S_n = tensors.Symmetric(common.usym(self.stress_n)) self.D = tensors.Symmetric(common.usym(self.d)) self.W = tensors.Skew(common.uskew(self.w)) self.strength_n = 25.0 self.strength_np1 = 30.0 self.H_n = history.History() self.H_n.add_scalar("strength") self.H_n.set_scalar("strength", self.strength_n) self.H_np1 = history.History() self.H_np1.add_scalar("strength") self.H_np1.set_scalar("strength", self.strength_np1) self.dt = 2.0 self.fixed = self.kmodel.decouple(self.S_n, self.D, self.W, self.Q, self.H_n, self.L, self.T, history.History()) self.ts = singlecrystal.SCTrialState(self.D, self.W, self.S_n, self.H_n, self.Q, self.L, self.T, self.dt, self.fixed) self.x = np.zeros((self.model.nparams, )) self.x[:6] = self.stress_np1 self.x[6] = self.strength_np1 self.Ddir = np.array([0.01, -0.005, -0.003, 0.01, 0.02, -0.003]) * 2 self.Wdir = np.array([0.02, -0.03, 0.01]) * 2 self.nsteps = 10
from neml import elasticity, drivers import matplotlib.pyplot as plt nthreads = 4 if __name__ == "__main__": N = 50 emax = 0.122 nsteps = 300 erate = 1.0e-4 nthreads = 4 orientations = rotations.random_orientations(N) lattice = crystallography.CubicLattice(1.0) lattice.add_slip_system([1, 1, 0], [1, 1, 1]) nslip = lattice.ntotal t0 = 100.0 b = 10.0 sat = 50 g0 = 1.0e-4 n = 6.0 mu = 29000.0 E = 120000.0 nu = 0.3
def setUp(self): self.L = crystallography.CubicLattice(1.0) self.L.add_slip_system([1, 1, 0], [1, 1, 1]) self.nslip = self.L.ntotal self.strength = 35.0 self.c = 10.0 self.beta = 2.0 self.H = history.History() for i in range(12): self.H.add_scalar("strength" + str(i)) self.H.set_scalar("strength" + str(i), self.strength) for j in range(4): self.H.add_scalar("slip_damage_" + str(j)) self.H.set_scalar("slip_damage_" + str(j), self.c * 0.4) self.static = 20.0 self.s0 = [self.static] * self.nslip self.k = 1000.0 self.sat = 40.0 self.m = 1.5 self.strengthmodel = slipharden.VocePerSystemHardening( self.s0, [self.k] * self.nslip, [self.sat] * self.nslip, [self.m] * self.nslip) self.g0 = 1.0 self.n = 3.0 self.slipmodel = sliprules.PowerLawSlipRule(self.strengthmodel, self.g0, self.n) self.imodel = inelasticity.AsaroInelasticity(self.slipmodel) self.Q = rotations.Orientation(35.0, 17.0, 14.0, angle_type="degrees") self.S = tensors.Symmetric( np.array([[100.0, -25.0, 10.0], [-25.0, -17.0, 15.0], [10.0, 15.0, 35.0]])) self.T = 300.0 self.mu = 29000.0 self.E = 120000.0 self.nu = 0.3 self.emodel = elasticity.CubicLinearElasticModel( self.E, self.nu, self.mu, "moduli") self.dn = np.array([[4.1, 2.8, -1.2], [3.1, 7.1, 0.2], [4, 2, 3]]) self.dn = 0.5 * (self.dn + self.dn.T) self.d = tensors.Symmetric(self.dn) self.wn = np.array([[-9.36416517, 2.95527444, 8.70983194], [-1.54693052, 8.7905658, -5.10895168], [-8.52740468, -0.7741642, 2.89544992]]) self.wn = 0.5 * (self.wn - self.wn.T) self.w = tensors.Skew(self.wn) self.dmodel = crystaldamage.WorkPlaneDamage() self.nfunc = crystaldamage.SigmoidTransformation(self.c, self.beta) self.sfunc = crystaldamage.SigmoidTransformation(self.c, self.beta) self.dmodel = crystaldamage.PlanarDamageModel(self.dmodel, self.nfunc, self.sfunc, self.L) self.model = kinematics.DamagedStandardKinematicModel( self.emodel, self.imodel, self.dmodel) self.fspin = self.model.spin(self.S, self.d, self.w, self.Q, self.H, self.L, self.T, history.History()) self.fixed = self.model.decouple(self.S, self.d, self.w, self.Q, self.H, self.L, self.T, history.History())