Abaqus(FEM)-Probleme bei Heat transfer Studie

Hallo zusammen,
Ich wollte in Abaqus folgendes Problem modellieren:
Ich möchte ein Stahlblech mit den Abmessungen 70 mm*70 mm* 0,5 mm erwärmen. Dies möchte ich mit dem Step HEAT TRANSFER machen. Um eine Erwärmung zu erzielen, habe ich an einer zuvor definierten Elementgruppe eine leistung in watt vorgegeben. Funktionierte für den instationären Fall auch prima. Allerdings funktioniert es nicht für steady state. Ich füge unten mal sa von mir geschriebene Input-File ein. Dazu habe ich noch ein psf skript geschrieben, um die Meshgröße, sowie stationär und instationär variieren zu können. Meine Vermutung liegt darin, dass ich den Befehl DFLUX eventuell nicht mit steady state verwenden kann bzw. vielleicht zusätzlich noch was dazufügen muss??
Ich arbeite übrigens mit Abaqus 6.9 und bin leider auch im Handbuch noch nicht fündig geworden.

Danke schon mal für alle Antworten

*Parameter
nel = 25 # number of elements
nlu = 101 # node links unten
nru = nlu + nel # node rechts unten
nlo = nlu + 100*nel # node links oben
nro = nlo + nel # node rechts oben
nhi1= 10101 # node hinten links untern
nhi2=nhi1+nel # node hinten rechts unten
nhi3=nhi1+nel*100 # node hinten links oben
nhi4=nhi3+nel # node hinten rechts oben
ellu = 101 # element links unten
elru = ellu + nel - 1 # element rechts unten
ello = 101 + 100*(nel-1) # element links oben
elro = ello + nel - 1 # element rechts oben
lat = 70.26 # latent heat
ht = ‚DELTMX = 6.0‘
X = elro - (nel/25-1)
*************************************************************************
*HEADING
STEEL HEATING 3D
*NODE
,0.,0.
,0.035,0.
,0.,0.035
,0.035,0.035
,0.,0.,0.0005
,0.035,0.,0.0005
,0.,0.035,0.0005
,0.035,0.035,0.0005
*NGEN,NSET=UNTEN
,1
,1
*NGEN,NSET=OBEN
,1
,1
*NGEN,NSET=LINKS
,100
,100
*NGEN,NSET=RECHTS
,100
,100
*NFILL,NSET=ALL
UNTEN,OBEN,100
**
*ELEMENT,TYPE=DC3D8
101,101,102,202,201,10101,10102,10202,10201
*ELSET,ELSET=X0
,100
*ELSET,ELSET=XR
,100
*ELSET,ELSET=Y0
,1
*ELSET,ELSET=YR
,1
*ELSET,ELSET=out,GENERATE
,1
*ELGEN,ELSET=ALL
101,1,1,100,100
*SOLID SECTION,ELSET=ALL,MATERIAL=STEEL
*MATERIAL,NAME=STEEL
*CONDUCTIVITY
57.,
*DENSITY
7850.,
*SPECIFIC HEAT
430.,
*LATENT HEAT
,85.0,95.0
***********************************************************
** Amplitude definition *
***********************************************************
*INITIAL CONDITIONS,TYPE=TEMPERATURE
ALL,22.0,
***********************************************************
** Output at time points *
***********************************************************
**TIME POINTS,GENERATE,NAME=tp
**0.0,5.0,0.05
***********************************************************
**RESTART,WRITE,FREQUENCY=25
*STEP,INC=100000000000
*HEAT TRANSFER,
.000001,6.0,0.001,0.01
*CONTROLS,PARAMETERS=LINE SEARCH
20
*SOLUTION TECHNIQUE,TYPE=QUASI-NEWTON,REFORM KERNEL=20
*DFLUX
out,S1,2450.
**NODE PRINT,FREQUENCY=25
**NT,
**NODE FILE,NSET=FOUT
**NT,
*OUTPUT,FIELD
*NODE OUTPUT,NSET=ALL
NT
*ELEMENT OUTPUT,ELSET=ALL
HFL,NFLUX
*END STEP

Moien

Allerdings funktioniert es nicht für steady state.

Es gibt bei der Aufgabe keinen steady state. Du hast eine Energiezufuhr, aber keine Kühlung. Die gesamte Energie geht in die Platte und die Temperatur steigt und steigt und steigt und … irgendwann steigt das System aus.

Definier eine Umgebung mit fixer Temperatur und einen Wärmeübergang zu dieser Umgebung. Dann geht das.

cu