<\/p>\n
from abaqusConstants import *\r\nfrom math import *\r\nimport odbAccess\r\nimport os\r\nimport multiprocessing as mp\r\n\r\n################# CONFIGS #####################\r\nodbFiles = ['MICRO_COMPR_XX', 'MICRO_COMPR_YY', 'MICRO_SHEAR_XY'] # seperate new odb files with comma\r\nmaxPercent = 0.0001 # Percent of elements with maximum values which is taken to determinate average maximum value \r\ncpus = 4 # How many CPUs do you want to use - how many odbs will be accessed in the same time. \r\n###############################################\r\n\r\ndef avgMaxFromPercentList(l, percent, revSort):\r\n l.sort(reverse=revSort) \r\n maxElements = int(ceil(len(l) * maxPercent))\r\n lmax = l[:maxElements]\r\n avg = sum(lmax) \/ float(len(lmax))\r\n return avg\r\n\r\n\r\ndef Analise(resultFile):\r\n print resultFile + \"\\t[START]\"\r\n\r\n try:\r\n os.remove(resultFile + '_MAX.txt')\r\n except OSError:\r\n pass\r\n \r\n odb = odbAccess.openOdb(resultFile+'.odb', readOnly=True)\r\n step=odb.steps['Step-1']\r\n\r\n frames = step.frames\r\n\r\n maxStrain = 0\r\n maxStress = 0\r\n # Stress\/Strain Curve\r\n for frame in frames:\r\n # Max stress \/ strain \r\n stress=frame.fieldOutputs['S'].getSubset(position=INTEGRATION_POINT).values\r\n strain=frame.fieldOutputs['PEEQ'].getSubset(position=INTEGRATION_POINT).values \r\n \r\n count = len(stress)\r\n stressList = [None] * count\r\n strainList = [None] * count\r\n for i in range(0, count):\r\n stressList[i] = stress[i].mises \r\n strainList[i] = strain[i].data \r\n \r\n\r\n newStress = avgMaxFromPercentList(stressList, maxPercent, True)\r\n if newStress > maxStress:\r\n maxStress = newStress\r\n \r\n newStrain = avgMaxFromPercentList(strainList, maxPercent, True)\r\n if newStrain > maxStrain:\r\n maxStrain = newStrain\r\n \r\n \r\n file = open(resultFile+'_MAX.txt', 'w')\r\n file.write('{0}\\t{1}\\n'.format(maxStrain, maxStress))\r\n file.close()\r\n odb.close()\r\n \r\n print resultFile + \"\\t[DONE]\"\r\n\r\n \r\nif __name__ == '__main__': \r\n pool = mp.Pool(processes=cpus)\r\n pool.map(Analise, odbFiles)\r\n pool.close() # no more tasks\r\n pool.join() # wrap up current tasks\r\n \r\n print \"Finished!\"\r\n\r\n<\/pre>\nSkrypt pobiera warto\u015b\u0107 maksymaln\u0105 napr\u0119\u017ce\u0144 Misesa i ekwiwalentnych odkszta\u0142ce\u0144 z n% element\u00f3w z najwi\u0119kszymi warto\u015bciami, w\u015br\u00f3d wszystkich krok\u00f3w czasowych symulacji.<\/p>\n","protected":false},"excerpt":{"rendered":"
from abaqusConstants import * from math import * import odbAccess import os import multiprocessing as mp ################# CONFIGS ##################### odbFiles = [’MICRO_COMPR_XX’, 'MICRO_COMPR_YY’, 'MICRO_SHEAR_XY’] # seperate new odb files with comma maxPercent = 0.0001 # Percent of elements with maximum values which is taken to determinate average maximum value… Continue reading