#/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)
