Output Mesh Data to ASCII format

pichon wrote:
Hi,

Have you already found the solution?

Thank You

Hello,

No I do not have the solution.  Unfortunately, I am still waiting for somebody to help.

JMB

Grigory Zhivotovsky wrote:

JMB wrote:
Hello all,

I would like to output a Salome generated mesh using a TUI to an ASCII file.  I have already been able to export the nodes using a snippet of code I saw somewhere here.  Could someone please provide the snippet for writing the elements too please?

Thanks
JMB

Hello,

I suggest you a small example how to write mesh data in ASCII file. Sure, you can rearrange data as well as to add data for polygons, pyramids etc.

Best Regards,

Grigory

Grigory,

Thank you very much for your help.  I will try it!

JMB

Grigory Zhivotovsky wrote:

JMB wrote:
Hello all,

I would like to output a Salome generated mesh using a TUI to an ASCII file.  I have already been able to export the nodes using a snippet of code I saw somewhere here.  Could someone please provide the snippet for writing the elements too please?
Thanks
JMB

Hello,

I suggest you a small example how to write mesh data in ASCII file. Sure, you can rearrange data as well as to add data for polygons, pyramids etc.

Best Regards,

Grigory

Hello Grigory,

Thank you for your assistance and the script that you emailed.  I tried your script, but it would not work.  I got an error:
>>> execfile("/home/.../write_mesh.py"
Traceback (most recent call last):
  File "", line 1, in ?
  File "/home/.../write_mesh.py", line 27, in ?
    M_nodes = Mesh_1.NbNodes()
AttributeError: Mesh instance has no attribute 'NbNodes'
>>>

I found the error and fixed it.  The script needs something like :
   mesh=Mesh_1.GetMesh()
and all references to Mesh_1.* changed now to mesh.* (such as "nodes = mesh.GetNodesId()" ).  Then your script works fine.  Does this seem logical to you?

Next I would like to use this script after I have created a geometric model of my own manually; meshed it like I want (manually); added groups, etc.  ( I want to be able to write the groups as well.  So later, I would like some help with that too please ).  However when I tried it I get an error:

>>> execfile("/home/../write_mesh-1.py"
Traceback (most recent call last):
  File "", line 1, in ?
  File "/home/.../write_mesh-1.py", line 18, in ?
    mesh=Mesh_1.GetMesh()
NameError: name 'Mesh_1' is not defined

The script "write_mesh-1.py" is a slightly modified version with all the geometry creation & meshing lines removed as shown below.  What do I need to do to make "write_mesh-1.py" a generic type of script that can be used after ANY model has been meshed?

+++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
import salome
salome.salome_init()
import geompy
import smesh
import SMESH
import StdMeshers
import os
import sys

dir=os.getenv("HOME"
File = dir + "/mesh_info.txt"
if os.access(File, os.F_OK):
            os.remove(File)
sys.stdout = open (File, "w"

# Get mesh entities and write info into ascii file
error=0
mesh=Mesh_1.GetMesh()
M_nodes = mesh.NbNodes()
if M_nodes == 0:
    print "ERROR!!! The number of NODES is equal 0!!!"; error=error+1
else:
    nodes = mesh.GetNodesId()
    print "Node,  X, Y, Z"
    for ind in nodes:      
        coord=mesh.GetNodeXYZ(ind)
        print ind, coord[0], coord[1],  coord[2]
M_edges= mesh.NbEdges()
if M_edges != 0:
    edges = mesh.GetElementsByType(SMESH.EDGE)
    print "Edge,  N1, N2"
    for ind in edges:
        nnode1=mesh.GetElemNode(ind,0)
        nnode2=mesh.GetElemNode(ind,1)
        print ind,nnode1,nnode2
M_faces= mesh.NbTriangles()+mesh.NbQuadrangles()
if M_faces != 0:
    fac = mesh.GetElementsByType(SMESH.FACE)
    print "Faces,  N1, N2, N3, N4"
    for ind in fac:
        nb_nodes=mesh.GetElemNbNodes(ind)
        nnode1=mesh.GetElemNode(ind,0)
        nnode2=mesh.GetElemNode(ind,1)
        nnode3=mesh.GetElemNode(ind,2)
        if nb_nodes ==3:
            print ind,nnode1,nnode2,nnode3
        else:
            nnode4=mesh.GetElemNode(ind,3)
            print ind,nnode1,nnode2,nnode3,nnode4
M_volum = mesh.NbTetras()+ mesh.NbPrisms()
if M_volum != 0:
    vol = mesh.GetElementsByType(SMESH.VOLUME)
    print "Volumes,  N1, N2, N3, N4, N5, N6, N7, N8"
    for ind in vol:
        nb_nodes=mesh.GetElemNbNodes(ind)
        nnode1=mesh.GetElemNode(ind,0)
        nnode2=mesh.GetElemNode(ind,1)
        nnode3=mesh.GetElemNode(ind,2)
        nnode4=mesh.GetElemNode(ind,3)
        if nb_nodes ==4:
            print ind,nnode1,nnode2,nnode3,nnode4
        else:
            nnode5=mesh.GetElemNode(ind,4)
            nnode6=mesh.GetElemNode(ind,5)
            nnode7=mesh.GetElemNode(ind,6)
            nnode8=mesh.GetElemNode(ind,7)          
            print ind,nnode1,nnode2,nnode3,nnode4,nnode5,nnode6,nnode7,nnode8