---

Hallo,

ich habe ein eingescanntes Bauteil in *.stl vorliegen. Hiervon ausgehend möchte ich ein FE-Netz erstellen. Mich verwundert es, dass eine scheinbar einfache Konvertierung von Dreieicken (*.stl) in andere Dreiecke doch nicht so einfach möglich ist. Mir stehen folgende Tools zur Verfügung:

CATIA V5 R17
Abaqus CAE
HyperMesh 8.0/9.0 (leider ohne stl-Einlese-Lizenz)

Wie geht ihr bei sowas vor? Die Rekonstruktion des Körpers (Schnitte, Hilfslinien, Flächen erzeugen usw.) erschein mir zu umständlich.

------------------
Gruß

Christian

---

Poste mal einen aussagekräftigen Ausschnitt aus der Datei oder ein Beispiel. Man kann bestimmt mit überschaubarem Aufwand ein kleines Script schreiben, dass die Dreiecke in der stl in Dreieckselemente im Format für die .inp erzeugt. Falls man so ein Oberflächennetz eines Volumens bekäme, könnte man in CAE dieses leicht mit Tets füllen.

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Hallo Mustaine!

Hier mal ein kleiner Auszug:

facet normal -0.032117 0.909769 -0.413871
outer loop
    vertex -161.052353 -163.713348 -18.519823
    vertex -162.983444 -163.613007 -18.149408
    vertex -163.374664 -162.984009 -16.736383
endloop
endfacet
facet normal -0.287154 0.440071 -0.850811
outer loop
    vertex -18.913666 121.900505 -86.940163
    vertex -14.050548 124.272934 -87.354378
    vertex -16.267523 121.486763 -88.047249
endloop
endfacet
facet normal -0.286058 0.019311 -0.958018
outer loop
    vertex 392.822968 40.932137 -255.861374
    vertex 396.009735 46.636719 -256.697937
    vertex 398.822845 39.056053 -257.690735
endloop
endfacet
facet normal -0.917578 -0.022571 -0.396915
outer loop
    vertex 53.229713 77.057732 -116.101082
    vertex 52.353962 81.681900 -114.339493
    vertex 53.454327 78.442978 -116.699112
endloop
endfacet
facet normal -0.182382 -0.852094 -0.490584
outer loop
    vertex 734.503479 -272.289734 -296.838684
    vertex 743.632263 -272.738068 -299.453766
    vertex 741.842285 -274.068970 -296.476685
endloop
endfacet
facet normal 0.082182 0.952793 -0.292288
outer loop
    vertex -167.880722 91.992455 -45.120731
    vertex -164.858871 91.867889 -44.677132
    vertex -165.484756 91.219521 -46.966652
endloop
endfacet
facet normal 0.657533 -0.540560 0.524829
outer loop
    vertex 950.443726 -179.879227 -268.017731
    vertex 947.048767 -184.457870 -268.480286
    vertex 951.100952 -180.981766 -269.976746
endloop
endfacet
facet normal -0.014108 0.981940 -0.188666
outer loop
    vertex 822.884094 146.388870 -334.973816
    vertex 816.727417 146.340652 -334.764343
    vertex 814.954529 147.191040 -330.205811
endloop
endfacet
facet normal -0.264995 -0.281028 -0.922389
outer loop
    vertex 414.065033 13.750456 -256.447784
    vertex 420.721466 13.982141 -258.430695
    vertex 420.061462 9.874003 -256.989441
endloop
endfacet
facet normal -0.291821 -0.075686 0.953474
outer loop
    vertex 716.659485 -56.020805 -225.348892
    vertex 721.169617 -59.645687 -224.256256
    vertex 723.767883 -51.490807 -222.813690
endloop
endfacet
facet normal 0.108355 -0.172453 0.979040
outer loop
    vertex 620.516602 -176.611618 -267.941162
    vertex 617.304443 -184.785065 -269.025360
    vertex 620.543579 -186.320297 -269.654266
endloop
endfacet
facet normal 0.402324 -0.246686 0.881636
outer loop
    vertex -165.842026 91.890747 -43.660610
    vertex -166.116318 91.998306 -43.505344
    vertex -166.350555 91.894669 -43.427452
endloop
endfacet
facet normal -0.587756 0.118395 -0.800329
outer loop
    vertex -166.236359 92.293907 -42.957569
    vertex -167.399719 92.520813 -42.069641
    vertex -166.848328 95.540359 -42.027889
endloop
endfacet
facet normal 0.803106 -0.184290 0.566620
outer loop
    vertex 981.770081 -144.612259 -300.010284
    vertex 984.204834 -135.378769 -300.458069
    vertex 980.826355 -136.487488 -296.030182
endloop
endfacet
facet normal -0.898856 0.127124 0.419400
outer loop
    vertex -160.877411 -175.230286 -8.962713
    vertex -161.116257 -174.223099 -9.779893
    vertex -161.233231 -175.230286 -9.725312
endloop
endfacet
facet normal 0.128169 -0.195843 0.972223
outer loop
    vertex -165.773376 92.059715 -43.538185
    vertex -165.682220 92.477730 -43.465996
    vertex -166.116318 91.998306 -43.505344
endloop
endfacet

So ein Script wäre echt hilfreich. Daran habe ich noch gar nicht gedacht.

------------------
Gruß

Christian

---

Hallo Christian,

ich würde Dir nicht empfehlen, die Stl-Dreiecke direkt als Oberflächengitter zu verwenden, da die i.d.R nur irgendwie die Geometrie beschreiben und für FE-Anwendungen viel zu verzerrt sind.
Wenn Du willst, kann ich mir das Teil mal ansehen und vernetzen (für Abaqus) (kostenlos). Ich entwickle gerade einen Vernetzer und da wäre es interessanter Test.
Gruß
Ulrich

------------------
Dr.-Ing. Ulrich Heck
ulrich_heck@dhcae.de
http://www.dhcae.de

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Hallo Ulrich,

vielen Dank für das Angebot. Die stl-Dreiecke werde ich nicht direkt zur Berechnung verwenden sondern als Basis, um davon ausgehend ein besseres Netz zu erzeugen. Ich werde dein Angebot im Hinterkopf behalten.

Die Script-Variante finde ich im Moment doch reizvoller.

Gruß

Christian

------------------
Gruß

Christian

---

Die multiplen Knoten an einer Coordinate in CAE über Edit mesh->Nodes->Merge entfernen.

Nicht besonders elegant, aber es sollte reichen...

Code:

---

"""
Aufruf:
abaqus python /scriptname/ /dateiname/
"""

import sys

input=file(sys.argv[1],'r')
zeilen=input.readlines()
input.close()

filename1 = sys.argv[1].split('.')[0]+'.inp'
output1=file(filename1, 'w')
filename2 = sys.argv[1].split('.')[0]+'.inc'
output2=file(filename2, 'w')

elemnr = 1
nodenr = 1

print >> output1, '*Node'
print >> output2, '*Element, Type=S3'

for i in range(len(zeilen)):
    if zeilen[i].split()[0]=='outer':
        print >> output1, str(nodenr)+', '+zeilen[i+1].split()[1]+', '+zeilen[i+1].split()[2]+', '+zeilen[i+1].split()[3]
        n1 = nodenr
        nodenr = nodenr+1
        print >> output1, str(nodenr)+', '+zeilen[i+2].split()[1]+', '+zeilen[i+2].split()[2]+', '+zeilen[i+2].split()[3]
        n2 = nodenr
        nodenr = nodenr+1
        print >> output1, str(nodenr)+', '+zeilen[i+3].split()[1]+', '+zeilen[i+3].split()[2]+', '+zeilen[i+3].split()[3]
        n3 = nodenr
        nodenr = nodenr+1

        print >> output2, str(elemnr)+', '+str(n1)+', '+str(n2)+', '+str(n3)
        elemnr = elemnr+1
       

print >> output1, '**\n*Include, Input='+filename2

print '\nDone\n'

output1.close()
output2.close()

---

Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP

---

Hallo Mustaine,

vielen Dank. Hat geklappt!

Gruß

Christian

------------------
Gruß

Christian

[Diese Nachricht wurde von c-dot am 12. Jan. 2009 editiert.]