def test_dispplot(self): plot.dispPlot(self.model)
model.v = 1. # Initial slider velocity, generally is vlp(t=0) model.vref = 1. # Reference velocity, generally vlp(t=0) state1 = staterelations.DieterichState() state1.b = 0.01 # Empirical coefficient for the evolution effect state1.Dc = 10. # Critical slip distance model.state_relations = [state1] # Which state relation we want to use # We want to solve for 40 seconds at 100Hz model.time = np.arange(0,40.01,0.01) # We want to slide at 1 um/s for 10 s, then at 10 um/s for 31 lp_velocity = np.ones_like(model.time) lp_velocity[10*100:] = 10. # Velocity after 10 seconds is 10 um/s # Set the model load point velocity, must be same shape as model.model_time model.loadpoint_velocity = lp_velocity # Run the model! model.solve() # Make the phase plot plot.phasePlot(model) # Make a plot in displacement plot.dispPlot(model) # Make a plot in time plot.timePlot(model)
model.v = 1. # Initial slider velocity, generally is vlp(t=0) model.vref = 1. # Reference velocity, generally vlp(t=0) state1 = staterelations.DieterichState() state1.b = 0.01 # Empirical coefficient for the evolution effect state1.Dc = 10. # Critical slip distance model.state_relations = [state1] # Which state relation we want to use # We want to solve for 40 seconds at 100Hz model.time = np.arange(0, 40.01, 0.01) # We want to slide at 1 um/s for 10 s, then at 10 um/s for 31 lp_velocity = np.ones_like(model.time) lp_velocity[10 * 100:] = 10. # Velocity after 10 seconds is 10 um/s # Set the model load point velocity, must be same shape as model.model_time model.loadpoint_velocity = lp_velocity # Run the model! model.solve() # Make the phase plot plot.phasePlot(model) # Make a plot in displacement plot.dispPlot(model) # Make a plot in time plot.timePlot(model)