repsample_script.py

import os
os.chdir('x:/YuxiangWang/AbaqusFolder/YoshiModel/')

from setfiber import *

repsampleArr = np.genfromtxt('./csvs/repsample.csv', delimiter=',')
runFiber = True
runForce = True
doAnalysis = True
skipWait = True

stimBlock = setHoldDisplArray(np.r_[.025:.15:6j])
stimBlockForce = setHoldForceArray(np.r_[1.3e-3:7.8e-3:6j])

def runSimulation(baseModelName, materialBlock):
    # Run displacement controlled analysis
    fiber = Fiber(baseModelName=baseModelName, materialBlock=materialBlock,
                  stimBlock=stimBlock,
                  runFiber=runFiber, doAnalysis=doAnalysis, skipWait=skipWait)
    # Now to force controlled
    # Note: need to refactor the code to tune force levels yet 
    # if cylinder radius is not 0.5
    if runForce:
        fiberForce = FiberForce(fiber, runFiber=runFiber, doAnalysis=doAnalysis,
                                stimBlock=stimBlockForce,
                                skipWait=skipWait)
    return


for sample_id, (tau1, tau2, g1, g2, ginf, mu, alpha, thickness) in enumerate(
        repsampleArr):
    materialBlock = copy.deepcopy(materialBlockDefault)
    materialBlock['skin_tau_array'] = [np.inf, tau2, tau1]
    materialBlock['skin_g_array'] = [ginf, g2, g1]
    materialBlock['skin_mu'] = mu * alpha / 2
    materialBlock['skin_alpha'] = alpha
    materialBlock['thicknessAll'][0] = thickness
    materialBlock['cylinderRadius'] = .50
    runSimulation('RepSample%d' % sample_id, materialBlock)


# Get distribution csvs
from readodb import getSurfaceDistribution
for fname in os.listdir('./odbs'):
    if fname.endswith('.odb') and fname.startswith('RepSample'):
        getSurfaceDistribution(fname[:-4])