Hallo zusammen,
ich habe ein einfaches Bauteil, welches aus Stahl und Alu besteht. Dieses soll von 450° auf Raumtemp abgekühlt werden, dabei entstehen Eigenspannungen (Druck im Stahl, Zug im Alu). Diese möchte ich mit Ansys berechnen.Dafür gibt es nun mehrere Wege: transiente Abkühleun per Konvektion und dann Wechsel auf Strukturdynamik oder gleich solid227 Elemente. ich hab nun diesbezüglich das Forum durchsucht und auch einige Einträge gefunden, die in die ähnliche Richtung gehen nur fällt es mir etwas schwer die Programme nachzuvollziehen, deswegen diese neue Anfrage.
Das Material habe ich für Alu elastisch-plastisch modelliert:
TB,PLASTIC,1,9,11,MISO
TBTEMP,20.0
TBPT,DEFI,0,70
TBPT,DEFI,0.002,71
TBPT,DEFI,0.005,73
TBPT,DEFI,0.01,78
TBPT,DEFI,0.02,81
TBPT,DEFI,0.05,85.000
TBPT,DEFI,0.1,100.000 (sind eigentlich 11 Datenpaare)
etc.
wird für Temperaturen bis 500° wiederholt.(Daten aus Warmzugproben)
Den Stahl habe ich elastisch modelliert.
Für beide Werkstoffe habe ich auch temp.-abhängige E-Moduli und Wäremausdehnungskoeffizienten angegeben
Beispiel Stahl:
MPTEMP,1,20,50,100,200,250,300
MPTEMP,7,350,400,500
MPDATA,Ex,2,1,210000,205000,190000,186000,182000,179000
MPDATA,EX,2,7,174000,170000,166000
MPDATA,ALPX,2,1,13.1e-6,13.8e-6,14.2e-6,14.7e-6,15.1e-6,18.5e-6,
MPDATA,ALPX,2,7,1.58e-6,15.9e-6,1.65e-6
(Alu ebenso)
Meine Frage ist nun wie ich die Eigenspannungen berechne??
Die Anfangstemp setzte ich über:
time,0.0001
IC,all,temp,450(Knoten sind vorher alle ausgewählt)
tref ist 450° (negative Dehnungen, da sich das Bauteil ja zusammenzieht)
Hier der Lösungsteil:
antype,trans
nlgeom,on
TOFFST,273 !Temperatur in Celsius
TRNOPT,full,,nodamp,,,
autots,on
solc,on
nropt,full
allsel
tref,450
time,0.0001
nsel,s,node,,all
IC,all,temp,450
allsel
solv
!Belastung
deltim,2,1,10
time,1000 !Sekunden
autots,on
kbc,1
asel,s,ext
sfa,all,,conv,20e-3,20
allsel
solve
Es folgt per etchg, tts der Wechsel auf die strukturdynamischen Elemente und die Angabe der oben beschriebenen Materialdaten.
Die Tempverteilung lese ich über:
ldread,Temp,2,,,,,rth
solve
ein.
Nur scheint hier der letzte Substep der transienten Berechnung verwendet zu werden, die gesamte Abkühlung über die Zeit wird m.M nicht berücksichtigt.
Per Schleifenprogrammierung kann ich auch nicht auf die einzelnen Substeps zurückgreifen.
Falls mir jemand Hinweise bzgl der Fehler oder evtl sogar eine Anleitung geben könnte wäre ich dankbar.
Ich komme bei diesem Thema irgendwie nicht mehr weiter(Grünes Buch und Hilfe habe ich schon durchforstet).
Gruß
Eine Antwort auf diesen Beitrag verfassen (mit Zitat/Zitat des Beitrags) IP