Weź druki z banku/banków Kadry – C2 – pokój 3 (generalnie opis jest mylący, bo wchodzi się przez 3 ale idzie się do 4) Jeśli na zaświadczeniach są wplywy z umów o dzieło – to C1 – pokój 21 (obsługa stron do 13:00) Podpis Dyrektora ds. Pracowniczych – C2 –… Continue reading
Bytescout Bitmap Visualizer for Visual Studio
Założenia: zdefiniowane sety z węzłów na każdej krawędzi jako LEFT, RIGHT, TOP, BOTTOM i stworzone dwa punkty referencyjne w złożeniu oraz dwa sety z nimi o nazwach REF-PT-1 i REF-PT-2
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 |
# -*- coding: mbcs -*- # Do not delete the following import lines from abaqus import * from abaqusConstants import * import collections import __main__ import section import regionToolset import displayGroupMdbToolset as dgm import part import material import assembly import optimization import step import interaction import load import mesh import job import sketch import visualization import xyPlot import displayGroupOdbToolset as dgo import connectorBehavior modelName = 'COMPRESSION_XY' periodicityLeftRight = True periodicityBottomTop = True #modelName = getInput('Enter model name:') root = mdb.models[modelName].rootAssembly # Clear all sets and equations connected to previous constaints names = [set[0] for set in root.sets.items()] for name in names: if "PERIOD-" in name: del root.sets[name] names = [c[0] for c in mdb.models[modelName].constraints.items()] for name in names: if "PERIOD-" in name: del mdb.models[modelName].constraints[name] # Get nodes from edges. Sets need to be prepared before this step leftSet = root.sets['LEFT'].nodes rightSet = root.sets['RIGHT'].nodes bottomSet = root.sets['BOTTOM'].nodes topSet = root.sets['TOP'].nodes # Connect Left to Right node if periodicityLeftRight: idsLR = collections.OrderedDict() for nodeL in leftSet: for nodeR in rightSet: if nodeL.coordinates[1] == nodeR.coordinates[1]: idsLR[nodeL.label] = nodeR.label sortedLeft = tuple(idsLR.keys()) sortedRight = tuple(idsLR.values()) # Connect Bottom to Top node if periodicityBottomTop: idsBT = collections.OrderedDict() for nodeB in bottomSet: for nodeT in topSet: if nodeB.coordinates[0] == nodeT.coordinates[0]: idsBT[nodeB.label] = nodeT.label sortedBottom = tuple(idsBT.keys()) sortedTop = tuple(idsBT.values()) i = 1 for n1, n2 in zip(sortedLeft, sortedRight): leftSetName = "PERIOD-LEFT"+str(i) rightSetName = "PERIOD-RIGHT"+str(i) root.SetFromNodeLabels(name=leftSetName, nodeLabels=(('PART-1-1', (n1,)), )) root.SetFromNodeLabels(name=rightSetName, nodeLabels=(('PART-1-1', (n2,)), )) eqName = "PERIOD-1-"+str(i) mdb.models[modelName].Equation(name=eqName, terms=((1.0, leftSetName, 1), (-1.0, rightSetName, 1), (1.0, 'REF-PT-1', 1), (0, 'REF-PT-1', 1))) # eqName = "PERIOD-2-"+str(i) # mdb.models[modelName].Equation(name=eqName, terms=((1.0, leftSetName, 2), (-1.0, rightSetName, 2), (1.0, 'REF-PT-2', 2), (0, 'REF-PT-2', 2))) i += 1 for n1, n2 in zip(sortedBottom, sortedTop): bottomSetName = "PERIOD-BOTTOM"+str(i) topSetName = "PERIOD-TOP"+str(i) root.SetFromNodeLabels(name=bottomSetName, nodeLabels=(('PART-1-1', (n1,)), )) root.SetFromNodeLabels(name=topSetName, nodeLabels=(('PART-1-1', (n2,)), )) eqName = "PERIOD-2-"+str(i) mdb.models[modelName].Equation(name=eqName, terms=((1.0, bottomSetName, 2), (-1.0, topSetName, 2), (1.0, 'REF-PT-2', 2), (0, 'REF-PT-2', 2))) # eqName = "PERIOD-1-"+str(i) # mdb.models[modelName].Equation(name=eqName, terms=((1.0, bottomSetName, 1), (-1.0, topSetName, 1), (1.0, 'REF-PT-1', 1), (0, 'REF-PT-1', 1))) i += 1 |
Format wyjściowy biblioteki deal.ii
|
1 2 3 4 5 6 |
(0,0) 6.73636e-09 (0,1) 2.76187e-09 (0,2) -3.27721e-09 (0,3) -3.6236e-11 (0,4) 1.13877e-08 (0,5) -5.36419e-10 |
Matlab wymaga numeracji od 1 i linii ze spacjami zaraz przecinków:
|
1 2 3 4 5 6 |
1 1 6.73636e-09 1 2 2.76187e-09 1 3 -3.27721e-09 1 4 -3.6236e-11 1 5 1.13877e-08 1 6 -5.36419e-10 |
Skrypt dokonujący konwersji:
|
1 2 3 4 5 6 7 8 9 |
#!/usr/bin/python import sys import re with open(sys.argv[1]+".new", 'a+') as w: with open(sys.argv[1]) as file: for line in file: m = re.match(r'\((.*?),(.*?)\)(.*?)$', line) w.write('{0} {1} {2}\n'.format(int(m.group(1))+1, int(m.group(2))+1, m.group(3))) |
Uruchomienie:
|
1 |
./matrix_convert.py A.dat |
Disable Windows 10 Popup
Skrypt dla Greasemonkey podmieniający Youtubową playerkę na stronie naekranie.pl na oryginalną YT. Działa również w chromowym Temperedmonkey – ale nie pozwala na full screena (Work in progress)
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 |
// ==UserScript== // @name NaEkranie.pl // @description Youtube without commercials for naekranie.pl // @namespace https://kszyh.kisim.eu.org // @include http://naekranie.pl/* // @version 1.0 // @grant none // @require http://ajax.googleapis.com/ajax/libs/jquery/2.1.3/jquery.js // ==/UserScript== $(document).ready(function() { $("iframe").each(function() { var src= $(this).attr('src'); if(src.indexOf("youtube") > -1){ var index = src.indexOf('?'); src = src.substring(0, index); src = src.replace('embed', 'v'); $(this).closest('.media-browser').replaceWith("<embed allowfullscreen='true' width='100%' height='400' src='"+src+"&vq=hd1080'>"); } }); }); |
Download: ver. 1.0
Iron Python:
|
1 |
ipy -X:TabCompletion |
Python: Do .pythonrc
|
1 2 |
import rlcompleter, readline readline.parse_and_bind('tab:complete') |
export PYTHONSTARTUP=~/.pythonrc lub: echo „export PYTHONSTARTUP=~/.pythonrc” >> ~/.bashrc
|
1 2 |
*EL FILE,FREQ=1 S |
|
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 |
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 cpus = 4 # How many CPUs do you want to use - how many odbs will be accessed in the same time. ############################################### def avgMaxFromPercentList(l, percent, revSort): l.sort(reverse=revSort) maxElements = int(ceil(len(l) * maxPercent)) lmax = l[:maxElements] avg = sum(lmax) / float(len(lmax)) return avg def Analise(resultFile): print resultFile + "\t[START]" try: os.remove(resultFile + '_MAX.txt') except OSError: pass odb = odbAccess.openOdb(resultFile+'.odb', readOnly=True) step=odb.steps['Step-1'] frames = step.frames maxStrain = 0 maxStress = 0 # Stress/Strain Curve for frame in frames: # Max stress / strain stress=frame.fieldOutputs['S'].getSubset(position=INTEGRATION_POINT).values strain=frame.fieldOutputs['PEEQ'].getSubset(position=INTEGRATION_POINT).values count = len(stress) stressList = [None] * count strainList = [None] * count for i in range(0, count): stressList[i] = stress[i].mises strainList[i] = strain[i].data newStress = avgMaxFromPercentList(stressList, maxPercent, True) if newStress > maxStress: maxStress = newStress newStrain = avgMaxFromPercentList(strainList, maxPercent, True) if newStrain > maxStrain: maxStrain = newStrain file = open(resultFile+'_MAX.txt', 'w') file.write('{0}\t{1}\n'.format(maxStrain, maxStress)) file.close() odb.close() print resultFile + "\t[DONE]" if __name__ == '__main__': pool = mp.Pool(processes=cpus) pool.map(Analise, odbFiles) pool.close() # no more tasks pool.join() # wrap up current tasks print "Finished!" |
Skrypt pobiera wartość maksymalną naprężeń Misesa i ekwiwalentnych odkształceń z n% elementów z największymi wartościami, wśród wszystkich kroków czasowych symulacji.
Styl do tworzenia bibliografii wg. formatu Computer Methods in Materials Science (Informatyka w technologii materiałów) dla Endnote.
