Tek-Tips is the largest IT community on the Internet today!

Members share and learn making Tek-Tips Forums the best source of peer-reviewed technical information on the Internet!

  • 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
Joined
Mar 8, 2010
Messages
4
Location
US
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...
 
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

Back
Top