Continue to Site

Eng-Tips is the largest engineering community on the Internet

Intelligent Work Forums for Engineering Professionals

  • Congratulations MintJulep on being selected by the Eng-Tips community for having the most helpful posts in the forums last week. Way to Go!

Save XYZ data

Status
Not open for further replies.

ColdCalifornian

Mechanical
Mar 8, 2010
4
I have an edge that I have partitioned and created a series of points on. I would like to output the xyz location of each point. Does anyone know how to do this? Right now I query each point and write it on a piece of paper...
 
Replies continue below

Recommended for you

Put the points in a set and run a data check and write the coordinates of that set to the dat file via *Node print, nset=....<CR>coords

Tara
 
If you don't want to rerun your simulation you can also use Probe to get the values and then save the results to a file.
 
Hi,

I'm struggling to obtain the external triangles of a tetrahedral mesh...

I did meshing operation with a Python script:

##### CREATE MODEL ######
myModel = mdb.Model(name=modelName)

##### IMPORT IGES MODEL ######
mdb.openIges('model.igs'
, msbo=True, scaleFromFile=OFF, trimCurve=DEFAULT)
myModel.PartFromGeometryFile(combine=False, convertToAnalytical=1
, dimensionality=THREE_D, geometryFile=mdb.acis, name='model',
stitchAfterCombine=False, stitchEdges=1, stitchTolerance=1.0, type=
DEFORMABLE_BODY)

[..assign materials, assembly, section, loads., etc etc....]


##### CREATE TETRAHEDRAL MESH #######
myModel.parts['part'].seedPart(size=8.0)
pickedRegions = model.parts[part'].cells
myModel.parts['part'].setMeshControls(elemShape=TET, regions=pickedRegions,
sizeGrowth=MAXIMUM, technique=FREE)
elemType = mesh.ElemType(elemCode=C3D4, elemLibrary=STANDARD,
distortionControl=DEFAULT)
pickedRegions = (myModel.parts['part'].cells,)
myModel.parts['part'].setElementType(regions=pickedRegions, elemTypes=(elemType,))
myModel.parts['part'].generateMesh()
myModel.rootAssembly.regenerate()
#########################################

I need to export the 3d faces that compose the external surface. I tried with:

for ef1 in myModel.parts['part'].elementFaces:
print ef1

#RESULT:
({'face': FACE1, 'label': 1})
({'face': FACE2, 'label': 1})
({'face': FACE3, 'label': 1})
({'face': FACE4, 'label': 1})
({'face': FACE1, 'label': 2})
({'face': FACE2, 'label': 2})
.... etc etc ......

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

for face in myModel.parts['part'].faces:
print face

#RESULT:
({'index': 0, 'instanceName': None, 'pointOn': ((16.626792, 29.168272, -50.806711),)})
({'index': 1, 'instanceName': None, 'pointOn': ((15.252681, -118.266439, -14.106842),)})
({'index': 2, 'instanceName': None, 'pointOn': ((-8.57784, -116.351098, -12.873274),)})
({'index': 3, 'instanceName': None, 'pointOn': ((-27.789027, -117.360474, 8.417591),)})
({'index': 4, 'instanceName': None, 'pointOn': ((-5.268457, -119.022915, -34.536412),)})
({'index': 5, 'instanceName': None, 'pointOn': ((-8.440513, -115.967772, -39.54762),)})
.... etc etc .....

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

but these results are not useful for my purposes.. I'm very beginning with Abaqus scripting, and I didn't find useful info in the documentation

please can you provide any suggestions?

thank you in advance, best regards
 
Status
Not open for further replies.

Part and Inventory Search

Sponsor