#/usr/bin/ipython import pylab from Network import * from Folder import * import os #------------------Liste des donnees d'entre--------------- SimNum = 1 repertoir = "/data/H.Le/simulation" + str(SimNum) + "/S/network/" NumFile = len(os.listdir(repertoir)) NumPlot = 2 Zoom = 2 Border = 0 #---------------------------------------------------------- Fold = Folder(Folder=repertoir, ParameterFile="Parameters.log", ResFile="Positions", Ext=".res") N = Network(filename=Fold.OpenDataFile(NumPlot), time=NumPlot * Fold.dt, sizex=Fold.Lx, sizey=Fold.Ly, sizez=Fold.Lz, distance=Fold.eps) N.plot2D_fiber(cracking=Fold.PCrack) N.AdjustWindows(zoom=Zoom, border=Border) N.show()
print(Fold2.OpenStressFile()) rm = "cd ./cracking && rm hole" + str(S) mkdir = "cd ./cracking && mkdir hole" + str(S) os.system(rm) os.system(mkdir) NumFile = len(os.listdir(Fichier)) if entier == 1: Nimage = NumFile for num in range(0, Nimage): n = NumFile - num - 2 print(n) time = n * Fold.dt try: N = Network(filename=Fold.OpenDataFile(n), time=time, sizex=Fold.Lx, sizey=Fold.Ly, sizez=Fold.Lz, distance=Fold.eps) except: print("fin de la boucle pour " + str(n)) break N.plot2D_fiber(cracking=Fold.PCrack) N.AdjustWindows(zoom=1.8, border=0) if num == numshow: N.show() #for i in range(0,Fold.pas): numero = str(n) #+i) if len(numero) == 1:
#/usr/bin/ipython import pylab from Folder import * from Network import * import os time=100 for SimNum=6 repertoir="/data/H.Le/simulation"+str(SimNum)+"/S/network/" k=1.1 Fold=Folder(Folder=repertoir,ParameterFile="Parameters.log",ResFile="Positions",Ext=".res") numero=time/(Fold.dt*Fold.pas) N=Network(filename=Fold.OpenDataFile(numero),time=numero*Fold.dt,sizex=Fold.Lx,sizey=Fold.Ly,sizez=Fold.Lz,distance=Fold.eps,k=k) N.tabl() N.StressProfile(direction=1,hole="false",mean="false",large=0,Nbande=10)