---

Hallo liebe User,

ich simuliere gerade die Temperaturverteilung innerhalb eines Ofens und bekomme folgende Fehlermeldungen:

*** WARNING ***                        CP =    8813.583  TIME= 17:25:40
Coefficient ratio exceeds 1.0e8 - Check results.                       

*** WARNING ***                        CP =    8857.787  TIME= 17:27:58
There are 1 small equation solver pivot terms.  May be a transitional 
effect caused by the Newton-Raphson logic.                             

*** WARNING ***                        CP =    8857.857  TIME= 17:27:59
Material property KXX of material 1 of element 19061 is evaluated at a 
temperature of -277.787222, which is below the supplied temperature   
range.  Temperature range checking terminates.                         

*** ERROR ***    SUPPRESSED MESSAGE    CP =    8860.361  TIME= 17:28:26
Absolute zero or below temperature -1135.98426 on radiation           
superelement 21675 node 392.                                           

Kann mir jemand hier bitte weiter helfen?

Mit freundlichem Gruß
earn

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

---

Input für das Matrix50:

TYPE,2  !SHELL57
ASEL,S,,,85
.....
AMESH,ALL

!Elementausrichtung
ASEL,S,,,93,94,1
....
AREVERSE,ALL

/AUX12

NSLE,S,ENAME,,SHELL57
STEF,5.68E-8
EMIS,1,0.6
EMIS,2,0.8
EMIS,4,0.8
EMIS,5,0.8
VTYPE,0,20
GEOM,0
WRITE

danach die shell57-elemente löschen und se einfügen.

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

---

Hallo earn,

ich habe auch schon ähnliche Fehlermeldungen gehabt. Diese deuten auf eine fehlerhafte Strahlungsmatrix hin. Unten siehst du einen beispielhaften Auszug einer meiner Logfiles. Die Vorgehensweise entspricht dem Manual und es haut hin.

Beispiel:
...
et,1,solid70
!------------------- 1. Weitere Elementtypen fuer thermische Berechnung -------------
et,2,shell57 !--> therm. Schalenelement fuer festlegen der Strahlungsrichtung
et,3,matrix50,1 !--> Superelement, das die Matrix der Formfaktoren aufnimmt
et,4,surf152 !--> Element zur Beruecksichtigung der Konvektion
keyopt,4,4,1 !--> Element hat 4 Knoten
keyopt,4,5,1 !--> extra-Node als Senke fuer Konvektion
keyopt,4,6,0 !--> extra-Node hat Bulktemperatur
keyopt,4,8,2 !--> Konvektion mit alpha bei gemittelter Temp.diff.

!------------------- 2. Space-Node als Strahlungssenke ------------------------------
n,120000,x_sp,y_sp,z_sp !Spacenode als Strahlungssenke definiert

!------------------- 3. Aufbau der Strahlungsmatrix ----------------------------------
!------------------- 3.1 Strahlungsnormalen der Oberflaechen festlegen --------------
nsel,s,loc,y,y1 !liefert die obere Oberflaeche
nsel,a,loc,x,x1 !linke Seite des Blocks
nsel,a,loc,x,x2 !rechte Seite
nsel,a,loc,z,z1 !Vorderseite
nsel,a,loc,z,z2 !Rueckseite

type,2 !Oberflaechenknoten auf Shell57 gesetzt
mat,1
esurf !Oberflaeche beschichtet mit Shell57

!------------------- 3.2 Strahlungsmatrix mit Hilfe von /aux12 generieren -----------
/aux12
esel,s,type,,2 !alle Shell57-Elemente selektiert
nsle !alle enthaltenen Knoten
nsel,a,,,120000 !Space-Node dazu
geom !0=3D
emis,1,0.9
stef,5.67e-8 !Stefan-Boltzmann-Konstante im MKS
vtype,0 !hidden, ist allgemeiner, beruecksichtigt Verdeckungen
space,120000 !Setzt Space-Node fest
write,stra_mat !--> ins File "stra_mat.sub" geschrieben

allsel,all,all

/prep7 !wieder Preprozessor aufrufen
type,3 !Elementtyp auf matrix50 schalten
se,stra_mat !Einlesen der Superelementmatrix=181
esel,s,type,,2 !Shell57-Elemente selektieren
edele,all !loeschen dieser Elemente
allsel,all,all

!------------------- 4. Einbeziehen der Konvektion ---------------------
nsel,s,loc,y,y1 !liefert die obere Oberflaeche
nsel,a,loc,x,x1 !linke Seite des Blocks
nsel,a,loc,x,x2 !rechte Seite
nsel,a,loc,z,z1 !Vorderseite
nsel,a,loc,z,z2   !Rueckseite

type,4 !Umschalten auf surf152
mat,1 !Material

sf,all,conv,8,20 !Konvektion mit alpha=8 [W/m^2*K], 20°C RT
esurf,120000 !Oberflaeche mit surf152 beschichtet
allsel,all,all
finish
...

------------------
Gruß,
aldebaran

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

---

/AUX12 kann Strahlungsaustausch-Koeffizienten aufgrund einer Schalenelement-Anordnung bestimmen (Schalenelemente als Fläche oder Oberfläche eines Volumens). Ich habe nicht gehört, dass der SPACE NODE da irgendwie einbezogen werden müsste. Hat der einen Formfaktor? Oder hat das einen anderen Hintergrund im Beispiel?

Zu den Fehlern: Evtl. TOFF,273 im Lösungsblock vergessen?
In dem Listing von Aldebaran die SHELL57 Elemente vielleicht vor dem Einfügen des Superelementes mit SE löschen.

Sven

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