Esempio n. 1
0
import sys
sys.path.append('/Users/dsoto/current/swDataFlat/code')
import roxanne as rx

dataFile = open('../data/separated/039-converted.data','r')
rx.getHeader(dataFile)
data = rx.readDataFileArray(dataFile)

position = data['positionLateralMicron']
force = data['forceLateralMicroNewton']

position = map(float,position)
force = map(float,force)

import matplotlib.pyplot as plt

plt.plot(position,force)

plt.title('039 - SWS 12 Shear Bending')
plt.xlabel('Position (microns)')
plt.ylabel('Force (microNewtons)')
plt.grid()
plt.show()

plt.savefig('039-shearBending.pdf')
Esempio n. 2
0
def main():
    import sys
    sys.path.append('/Users/dsoto/current/swDataFlat/roxanne')
    import roxanne
    import glob
    import os.path
    import numpy

    parseFileIn = open('../037-sps06-ls/data/parsed.dat',
                       'r')
    parseDict = roxanne.readDataFileArray(parseFileIn)
    #print parseDict.keys()

    fOut = open('analyzed.data','w')
    fileNameList = glob.glob('../037-sps06-ls/data/separated/*.data')
    outputList = ['fileName',
                  'pulloffAngle',
                  'forceMaxAdhesion',
                  'forceMaxShear',
                  'forcePreload']

    # removed angle pulloff from output list
    sep = '\t'
    headerString = sep.join(outputList)
    fOut.write(headerString)
    fOut.write('\n')

    for dataFileName in fileNameList:

        print 'processing file : ' + dataFileName

        # read in data file
        fileIn = open(dataFileName)
        headerDict = roxanne.getHeader(fileIn)
        dataDict = roxanne.readDataFileArray(fileIn)
        # this dataDict needs to be converted to floats
        for key in dataDict.keys():
            dataDict[key] = map(float, dataDict[key])

        # angle value should be in headerDict
        pulloffAngle = float(headerDict['pulloffAngle'])

        # get indices from the parsed.data file
        # get location of filename in parseDict
        # and pull data from same index in the other arrays
        parsedFileNames = parseDict['dataFileName']
        dataFileNameNoPath = os.path.split(dataFileName)[1]
        #fileNameNoPath = fileNames
        indexFileName = parsedFileNames.index(dataFileNameNoPath)

        indexContact     = parseDict['indexContact'][indexFileName]
        indexPreload     = parseDict['indexMaxPreload'][indexFileName]
        indexMaxAdhesion = parseDict['indexMaxAdhesion'][indexFileName]

        indexContact     = int(indexContact)
        indexPreload     = int(indexPreload)
        indexMaxAdhesion = int(indexMaxAdhesion)

        # pull force values at those indices
        # data will be in dataDict
        normalForceContactMuN = dataDict['forceNormalMicroNewton'][indexContact]
        normalForcePreloadMuN = dataDict['forceNormalMicroNewton'][indexPreload]
        normalForcePulloffMuN = dataDict['forceNormalMicroNewton'][indexMaxAdhesion]

        shearForceContactMuN  = dataDict['forceLateralMicroNewton'][indexContact]
        shearForcePreloadMuN  = dataDict['forceLateralMicroNewton'][indexPreload]
        shearForcePulloffMuN  = dataDict['forceLateralMicroNewton'][indexMaxAdhesion]

        # correct for force at contact
        # adhesion force = maxAdhesion - force at contact
        maxAdhesionMuN = (normalForcePulloffMuN -
                          normalForceContactMuN)
        # shear force = maxShear - force at contact
        maxShearMuN = (shearForcePulloffMuN -
                       shearForceContactMuN)

        # force at preload - force at contact = force of preload
        forcePreload = normalForcePreloadMuN - normalForceContactMuN

        fOut.write(dataFileNameNoPath + '\t')
        fOut.write('% 5.1f\t' % pulloffAngle )
        fOut.write('% 5.3f\t' % maxAdhesionMuN)
        fOut.write('% 5.3f\t' % maxShearMuN)
        fOut.write('% 5.3f\t' % forcePreload)
        fOut.write('\n')

    fOut.close()