from odbAccess import *
from odbSection import *
import csv
import time
import sys
from math import *
import math
import os


DB = 'abaqus.odb'
odb = openOdb(DB)

myassembly=odb.rootAssembly
Instance=myassembly.instances['PART-1-1']

liste_surfaces=[]
for surface in myassembly.surfaces.keys():
  liste_surfaces.append(surface)
 
nb_surfaces=len(liste_surfaces)

i=0
while i < nb_surfaces:
  print liste_surfaces[i] 
  i+=1 

  
liste_faces=[]
nb_faces=0
i=0
while i < nb_surfaces:
  liste_faces.append(myassembly.surfaces[liste_surfaces[i]].faces)
  #print liste_faces
 #Ici, je souhaite recuperer les faces d'elements associees a la surface i
  i+=1

##############################################

# i=0
# M_Element_GLOBAL=[]
## On recupere tous les elts et les met dans la liste M_Element_GLOBAL
# for node in Instance.elements :
	# M_Element_GLOBAL.append(Instance.elements[i].label)
	# i+=1
# i=0	
# file_element=open('elements.csv','w')
# for element in M_Element_GLOBAL:
  # file_element.write(str(M_Element_GLOBAL[i]) + ';')
  # i+=1
# file_element.write('\n')
  
###################
# Pour avoir les valeurs de RADFLA par increment (frame) pour chaque face d'element k
# Comment connaitre la valeur max de k ????????????????????
# Et savoir quelle face d'element est k ??????????
# Et a quelle surface geometrique k est associe ?????????
##############
# M_RADFLA_GLOBAL=[]
# i=0  
# for Step in odb.steps.keys() :
    # step_temp = odb.steps.keys()[i]
    # j=0
    # for incr in odb.steps[step_temp].frames:
	   # for k in ???????????????????????
        # M_RADFLA_GLOBAL.append(odb.steps[step_temp].frames[j].fieldOutputs['RADFLA'].values[k].data)
		# k+=1
    # j+=1	
# i+=1 		


M_RADFLA_GLOBAL=[]
i=0  
for Step in odb.steps.keys() :
    step_temp = odb.steps.keys()[i]
    j=0
    for incr in odb.steps[step_temp].frames:
     k=0
     while k < 34:
          M_RADFLA_GLOBAL.append(odb.steps[step_temp].frames[j].fieldOutputs['RADFLA'].values[k].data)
		  #Radfla_global est mal renseigne: il faudrait qu'il soit a 2 dim (instant,face d'element)
          k+=1
     j+=1	
i+=1