2017-12-25-172238.ipynb

Kernel: Python 3 (Anaconda)

# Imports and settings
%matplotlib inline
import csv, math, os.path
import numpy as np
import pylab
from matplotlib import pyplot as plt
plt.rcParams['figure.figsize'] = (16.0, 10.0)

# Function to read in the data
def read_data(Filename):
    # Read in data
    a = np.loadtxt(Filename, delimiter=";", skiprows=2)
    # Sort on the 5th column, return the result
    return (a[a[:,4].argsort()])





# Some constants that can be used in the remainder of the code
area = 3.5**2*math.pi # cross section area of the sample
cutoff_peak = 0.005   # we assume the load peak is at a strain (x-value) lower than this value
cutoff_min = 0.05     # we assume the minimum after the peak is at a strain lower than this value



#Function to calculate the strain and load from the loaded data
def calc_strain_load(data):
    
    strain = data[:,4]
    print("Här är kolumnen för strain-värdena: " + str(strain))
    
    load = data [:,2]
    spanning = load/area
    
    print("Här är kolumnen för spänning-värdena: " + str(spanning))    
    
    return strain,spanning


# Function to compute the peak strain and return the result
def find_peak(strain, spanning):
    
    #Vi skapar en ny matris med endast värden under 0,005 (cutoff_peak)
    newlength = 0 
    count1 = 0
    while strain[count1] < cutoff_peak: #vi loopar över raderna i strain kolumnen för att se om strain är mindre än cut off peak.
        newlength = newlength + 1
        count1 = count1 + 1
    
    topp = max(spanning[:newlength]) #Nu beräknar vi max värdet på spänningen genom att kolla i den nya matrisen

    rad = 0
    while spanning[rad] != topp: #Nu hittar vi raden som är lika med topp spänningen för att kunna plocka ut både spänning och                                                                                                                               strain
        rad = rad + 1
        
    peak_x = strain[rad] #Här ger vi de värdena nya variabelnamn
    peak_y = spanning[rad]
    
    print ("Max-värdet för detta dragtest har x-värdet: " + str(peak_x) + " och y-värdet: " + str(peak_y) + " och Rel är: " + str(topp) + " MPa")
    return peak_x, peak_y, rad


# Function to compute E-modul and return the result
def find_gradient(peak_x, peak_y):
    #Vi hittar lutningen på grafen då vi beräknar utgångspunktens X,Y minus slutpunktens X,Y. Utgångspunkten är noll, därmed blir det slutpunktens X,Y värden delat på varandra. delta y / delta x 
    e_modul = peak_y/peak_x
    print("E-modulen för detta dragtest är: " + str(e_modul) + " KPa")
    return e_modul

# Function to compute R_el and return the result
def find_strain_limit(strain, spanning, rad):
    #Vi skapar en ny matris igen med endast värden under 0,05 (cutoff_min)
    newlength1 = 0 
    count1 = 0
    while strain[count1] < cutoff_min: #vi loopar över raderna i strain kolumnen för att se om strain är mindre än cut off min.
        newlength1 = newlength1 + 1
        count1 = count1 + 1
    
    Botten = min(spanning[rad:newlength1]) #Nu beräknar vi min värdet på spänningen genom att kolla i den nya matrisen. Den nya matrisen går utifrån topp värdet vi beräknat tidigare och till värdet 0.05. 

    rad2 = 0
    while spanning[rad2] != Botten : #Nu hittar vi raden som är lika med topp spänningen för att kunna plocka ut både spänning och                                                                                                                               strain
        rad2 = rad2 + 1
    
    limit_x = strain [rad2]
    limit_y = spanning [rad2]
    
    print("Min-värdet för detta dragtest har x-värdet: " + str(limit_x) + " Newton" + " och y-värdet: " + str(limit_y) + " Newton")
    
    return limit_x, limit_y

# Function to process a file using the other functions defined above
def process_file(Filename):
    data = read_data(Filename)
    strain,spanning = calc_strain_load(data)
    peak_x,peak_y,rad = find_peak(strain,spanning)
    gradient = find_gradient(peak_x,peak_y)
    limit_x,limit_y = find_strain_limit(strain,spanning,rad)
    plt.plot(strain,spanning)
    plt.ylabel("Spänning")
    plt.xlabel("Töjning")
    pylab.plot(peak_x,peak_y,'go', label='Max-värde')
    pylab.plot(limit_x,limit_y,'ro', label='Min-värde')
    pylab.legend() 
    plt.show()


# Call function process_file for each of the filenames...

for j in range (1,29):
    Filename = ("Specimen_RawData_" + str(j) +".csv") #Vi skapar ett namn för filerna för att hitta dem i datorn.
    
    if os.path.exists(Filename): # if the file exists 
        print("Här är namnet på din inlästa fil " + (Filename)) 
        process_file (Filename)
        

Här är namnet på din inlästa fil Specimen_RawData_1.csv
Här är kolumnen för strain-värdena: [  1.00000000e-05   2.00000000e-05   2.00000000e-05 ...,   1.92480000e-01
   1.92540000e-01   1.96490000e-01]
Här är kolumnen för spänning-värdena: [  2.68350149e+00   2.13714557e-01   2.21033242e+00 ...,   5.90276293e+02
   5.83956587e+02   9.26199206e+01]
Max-värdet för detta dragtest har x-värdet: 0.00269 och y-värdet: 544.288637077 och Rel är: 544.288637077 MPa
E-modulen för detta dragtest är: 202337.7833 KPa
Min-värdet för detta dragtest har x-värdet: 0.0139 Newton och y-värdet: 486.984666476 Newton

Här är namnet på din inlästa fil Specimen_RawData_3.csv
Här är kolumnen för strain-värdena: [  1.00000000e-05   2.00000000e-05   2.00000000e-05 ...,   2.10870000e-01
   2.11020000e-01   2.11190000e-01]
Här är kolumnen för spänning-värdena: [  7.90816071  10.40604622   8.80214741 ...,   0.69187135   0.44613534
   0.57343435]
Max-värdet för detta dragtest har x-värdet: 0.00228 och y-värdet: 480.146694524 och Rel är: 480.146694524 MPa
E-modulen för detta dragtest är: 210590.655493 KPa
Min-värdet för detta dragtest har x-värdet: 0.00229 Newton och y-värdet: 451.362867997 Newton

Här är namnet på din inlästa fil Specimen_RawData_4.csv
Här är kolumnen för strain-värdena: [ -1.00000000e-05   0.00000000e+00   1.00000000e-05 ...,   1.93960000e-01
   1.94010000e-01   1.95940000e-01]
Här är kolumnen för spänning-värdena: [ -2.80645408e+00  -1.14716648e+00  -2.21463389e-01 ...,   5.89596953e+02
   5.82688510e+02   3.52175330e+00]
Max-värdet för detta dragtest har x-värdet: 0.00246 och y-värdet: 519.772555936 och Rel är: 519.772555936 MPa
E-modulen för detta dragtest är: 211289.656884 KPa
Min-värdet för detta dragtest har x-värdet: 0.01416 Newton och y-värdet: 480.124391785 Newton

Här är namnet på din inlästa fil Specimen_RawData_6.csv
Här är kolumnen för strain-värdena: [  0.00000000e+00   1.00000000e-05   1.00000000e-05 ...,   2.01470000e-01
   2.01630000e-01   2.01790000e-01]
Här är kolumnen för spänning-värdena: [ 3.27991626  4.00349673  4.77833497 ...,  1.31627296  0.02057971
  1.39737183]
Max-värdet för detta dragtest har x-värdet: 0.00244 och y-värdet: 511.896155276 och Rel är: 511.896155276 MPa
E-modulen för detta dragtest är: 209793.506261 KPa
Min-värdet för detta dragtest har x-värdet: 0.01226 Newton och y-värdet: 452.704125038 Newton

Här är namnet på din inlästa fil Specimen_RawData_7.csv
Här är kolumnen för strain-värdena: [  1.00000000e-05   2.00000000e-05   2.00000000e-05 ...,   2.00680000e-01
   2.00800000e-01   2.00960000e-01]
Här är kolumnen för spänning-värdena: [ 1.30768691  4.8936741   4.49513765 ...,  0.39349806  1.4105826
  0.60414645]
Max-värdet för detta dragtest har x-värdet: 0.00248 och y-värdet: 510.60984164 och Rel är: 510.60984164 MPa
E-modulen för detta dragtest är: 205891.065177 KPa
Min-värdet för detta dragtest har x-värdet: 0.01214 Newton och y-värdet: 449.286668509 Newton

Här är namnet på din inlästa fil Specimen_RawData_8.csv
Här är kolumnen för strain-värdena: [  0.00000000e+00   1.00000000e-05   1.00000000e-05 ...,   1.94660000e-01
   1.94820000e-01   1.94920000e-01]
Här är kolumnen för spänning-värdena: [   4.10123761    4.78609342    7.73900564 ...,  593.46514731  591.4847409
  586.64254725]
Max-värdet för detta dragtest har x-värdet: 0.00256 och y-värdet: 529.051778146 och Rel är: 529.051778146 MPa
E-modulen för detta dragtest är: 206660.850838 KPa
Min-värdet för detta dragtest har x-värdet: 0.01166 Newton och y-värdet: 485.303642345 Newton

Här är namnet på din inlästa fil Specimen_RawData_9.csv
Här är kolumnen för strain-värdena: [  0.00000000e+00   0.00000000e+00   1.00000000e-05 ...,   1.99330000e-01
   1.99480000e-01   1.99610000e-01]
Här är kolumnen för spänning-värdena: [ 3.541607    3.91770065  6.67643501 ...,  1.55812481  0.26404546
  1.73262801]
Max-värdet för detta dragtest har x-värdet: 0.00248 och y-värdet: 492.844573747 och Rel är: 492.844573747 MPa
E-modulen för detta dragtest är: 198727.650704 KPa
Min-värdet för detta dragtest har x-värdet: 0.0083 Newton och y-värdet: 484.69562369 Newton

Här är namnet på din inlästa fil Specimen_RawData_10.csv
Här är kolumnen för strain-värdena: [ 0.       0.       0.      ...,  0.19361  0.19363  0.19463]
Här är kolumnen för spänning-värdena: [   2.82890675    2.50961049    4.36161692 ...,  587.19852601  574.88370272
 -135.02519271]
Max-värdet för detta dragtest har x-värdet: 0.00262 och y-värdet: 541.221905294 och Rel är: 541.221905294 MPa
E-modulen för detta dragtest är: 206573.246295 KPa
Min-värdet för detta dragtest har x-värdet: 0.01233 Newton och y-värdet: 485.037210734 Newton

Här är namnet på din inlästa fil Specimen_RawData_11.csv
Här är kolumnen för strain-värdena: [ -1.00000000e-05   0.00000000e+00   0.00000000e+00 ...,   1.90430000e-01
   1.90590000e-01   1.90610000e-01]
Här är kolumnen för spänning-värdena: [   2.98835531    5.71564195    5.64213939 ...,  584.56855489  576.86201841
  579.06281842]
Max-värdet för detta dragtest har x-värdet: 0.00237 och y-värdet: 498.090654572 och Rel är: 498.090654572 MPa
E-modulen för detta dragtest är: 210164.833153 KPa
Min-värdet för detta dragtest har x-värdet: 0.00242 Newton och y-värdet: 459.688574365 Newton

Här är namnet på din inlästa fil Specimen_RawData_13.csv
Här är kolumnen för strain-värdena: [  0.00000000e+00   0.00000000e+00   1.00000000e-05 ...,   1.77460000e-01
   1.77490000e-01   1.80300000e-01]
Här är kolumnen för spänning-värdena: [ -6.65788651e-01  -1.26412211e+00   3.26941931e-02 ...,   5.72901167e+02
   5.63488748e+02  -1.13045391e+01]
Max-värdet för detta dragtest har x-värdet: 0.00249 och y-värdet: 509.788863797 och Rel är: 509.788863797 MPa
E-modulen för detta dragtest är: 204734.483453 KPa
Min-värdet för detta dragtest har x-värdet: 0.01312 Newton och y-värdet: 455.788009697 Newton

Här är namnet på din inlästa fil Specimen_RawData_16.csv
Här är kolumnen för strain-värdena: [  0.00000000e+00   0.00000000e+00   2.00000000e-05 ...,   2.00830000e-01
   2.00980000e-01   2.01140000e-01]
Här är kolumnen för spänning-värdena: [  2.94340398e+00   2.64043533e+00   6.37999055e+00 ...,   3.91066432e-03
   1.22596858e+00   7.71422060e-02]
Max-värdet för detta dragtest har x-värdet: 0.00253 och y-värdet: 533.868339406 och Rel är: 533.868339406 MPa
E-modulen för detta dragtest är: 211015.153915 KPa
Min-värdet för detta dragtest har x-värdet: 0.00929 Newton och y-värdet: 482.575864598 Newton

Här är namnet på din inlästa fil Specimen_RawData_18.csv
Här är kolumnen för strain-värdena: [ -4.00000000e-05  -2.00000000e-05  -2.00000000e-05 ...,   1.97050000e-01
   1.97080000e-01   1.98510000e-01]
Här är kolumnen för spänning-värdena: [ -14.27591828  -10.87739327   -9.74484357 ...,  579.22598407  571.77217214
  -32.1345323 ]
Max-värdet för detta dragtest har x-värdet: 0.00242 och y-värdet: 506.655069188 och Rel är: 506.655069188 MPa
E-modulen för detta dragtest är: 209361.598838 KPa
Min-värdet för detta dragtest har x-värdet: 0.01229 Newton och y-värdet: 444.601290418 Newton

Här är namnet på din inlästa fil Specimen_RawData_20.csv
Här är kolumnen för strain-värdena: [ -3.00000000e-05  -2.00000000e-05  -2.00000000e-05 ...,   2.04850000e-01
   2.04980000e-01   2.09360000e-01]
Här är kolumnen för spänning-värdena: [   5.65997956    6.79035617    5.4201584  ...,  586.0246163   579.57500501
  108.50960657]
Max-värdet för detta dragtest har x-värdet: 0.00239 och y-värdet: 527.018054708 och Rel är: 527.018054708 MPa
E-modulen för detta dragtest är: 220509.646321 KPa
Min-värdet för detta dragtest har x-värdet: 0.00333 Newton och y-värdet: 481.638455755 Newton

Här är namnet på din inlästa fil Specimen_RawData_21.csv
Här är kolumnen för strain-värdena: [ -2.00000000e-05  -1.00000000e-05   1.00000000e-05 ...,   1.97170000e-01
   1.97200000e-01   1.97420000e-01]
Här är kolumnen för spänning-värdena: [   7.79027796   10.12439784   10.90040746 ...,  575.75269339  562.30488725
 -130.18587398]
Max-värdet för detta dragtest har x-värdet: 0.00248 och y-värdet: 514.104079711 och Rel är: 514.104079711 MPa
E-modulen för detta dragtest är: 207300.032141 KPa
Min-värdet för detta dragtest har x-värdet: 0.01504 Newton och y-värdet: 453.268351268 Newton

Här är namnet på din inlästa fil Specimen_RawData_22.csv
Här är kolumnen för strain-värdena: [ -1.00000000e-05   0.00000000e+00   1.00000000e-05 ...,   2.01420000e-01
   2.01560000e-01   2.01710000e-01]
Här är kolumnen för spänning-värdena: [  7.58163817   6.24479044  11.11746154 ...,   1.1473042   -0.32014335
   1.08690613]
Max-värdet för detta dragtest har x-värdet: 0.00235 och y-värdet: 497.684715403 och Rel är: 497.684715403 MPa
E-modulen för detta dragtest är: 211780.729959 KPa
Min-värdet för detta dragtest har x-värdet: 0.00246 Newton och y-värdet: 448.211726929 Newton

Här är namnet på din inlästa fil Specimen_RawData_23.csv
Här är kolumnen för strain-värdena: [ -1.00000000e-05   3.00000000e-05   5.00000000e-05 ...,   1.94600000e-01
   1.94650000e-01   1.94710000e-01]
Här är kolumnen för spänning-värdena: [  -4.27436883   -2.33856063    5.87941746 ...,  581.59238889  578.79394739
  570.30945064]
Max-värdet för detta dragtest har x-värdet: 0.00251 och y-värdet: 516.379763875 och Rel är: 516.379763875 MPa
E-modulen för detta dragtest är: 205728.989592 KPa
Min-värdet för detta dragtest har x-värdet: 0.00739 Newton och y-värdet: 452.557918681 Newton

Här är namnet på din inlästa fil Specimen_RawData_25.csv
Här är kolumnen för strain-värdena: [ -1.00000000e-05   0.00000000e+00   2.00000000e-05 ...,   2.04500000e-01
   2.04620000e-01   2.04800000e-01]
Här är kolumnen för spänning-värdena: [-3.83677121  0.25366362 -1.30463062 ...,  0.18449059 -0.92169317
  0.13741944]
Max-värdet för detta dragtest har x-värdet: 0.00233 och y-värdet: 490.759197839 och Rel är: 490.759197839 MPa
E-modulen för detta dragtest är: 210626.265167 KPa
Min-värdet för detta dragtest har x-värdet: 0.0066 Newton och y-värdet: 454.48237478 Newton

Här är namnet på din inlästa fil Specimen_RawData_27.csv
Här är kolumnen för strain-värdena: [ -1.00000000e-05   0.00000000e+00   0.00000000e+00 ...,   1.95010000e-01
   1.95170000e-01   1.97580000e-01]
Här är kolumnen för spänning-värdena: [  49.79921791   50.86272294   53.99989579 ...,  585.38043348  583.38504261
  151.5245136 ]
Max-värdet för detta dragtest har x-värdet: 0.00234 och y-värdet: 539.853892813 och Rel är: 539.853892813 MPa
E-modulen för detta dragtest är: 230706.7918 KPa
Min-värdet för detta dragtest har x-värdet: 0.01068 Newton och y-värdet: 490.098326491 Newton

Här är namnet på din inlästa fil Specimen_RawData_28.csv
Här är kolumnen för strain-värdena: [ 0.       0.       0.      ...,  0.20714  0.2073   0.20745]
Här är kolumnen för spänning-värdena: [ 1.65498742  2.31109607  4.31857259 ..., -0.15568056  0.72555763
 -0.08624119]
Max-värdet för detta dragtest har x-värdet: 0.00261 och y-värdet: 541.912784299 och Rel är: 541.912784299 MPa
E-modulen för detta dragtest är: 207629.419272 KPa
Min-värdet för detta dragtest har x-värdet: 0.0061 Newton och y-värdet: 487.078446025 Newton