---

Hallo Zusammmen

Ich möchte die Verteilung einer Flächenlast auf den einzelnen Knoten in WB 14.5 untersuchen. Ich möchte z.B. wissen, wie viel Leistung pro Knoten in einer thermisch-transienten Analyse eingestellt wird, wenn ich eine Leistung pro Fläche vorgebe.

Z.B.:
- Flächenlast: 1W
- Knoten 195: 0.1w
- Knoten 196: 0.1W
- Knoten 197: 0.1W
- usw.

Mit folgenden APDL-Befehlen kann ich die gewünschten Informationen herausgeben lassen:

Code:

---

/prep7

nsel,S,,,M1                          ! Alle Knoten der Flaeche M1 auswaehlen
sf,all,hflux,10000                  ! 1W/cm2

*cfopen,vekt,data                    ! Datei  erzeugen
*cfclos                              ! Datei schliessen

*get,nanz,node,,count                ! Ausgewaehlte Knoten zaehlen
*get,nmin,node,,num,min              ! Kleinste Knotennummer bestimmen

*do,j,1,nanz,1                      ! DO-Schleife ueber alle Knoten der Flaeche M1
nsel,s,node,,nmin                    ! Knoten mit der kleinsten Nummer auswaehlen
*get,nheat,node,nmin,f,heat          ! Heat des Knotens nmin bestimmen
*cfopen,vekt,data,,append            ! Datei oeffnen zum Daten Anheften
*vwrite,nmin,nheat                  ! Variablen in Datei schreiben
(F10.0,F20.10)                      ! Datenformat definieren
*cfclos                              ! Datei schliessen
nsel,S,,,M1                          ! Alle Knoten der Flaeche Mi erneut auswaehlen
*get,nnext,node,nmin,nxth            ! Naechste Knotennummer nach nmin bestimmen
nmin=nnext                          ! Naechste Knotennummer in nmin speichern
*enddo

/solu

---

Die Werte in der Datei vekt.data zeigen allerdings, dass die Knoten KEINE Last haben:
            195.        0.0000000000
            196.        0.0000000000
            197.        0.0000000000
            usw.

Als vergleich habe ich die Simulation wiederholt und der Befehl sf,all,hflux,10000 mit dem Befehl f,all,heat,1 ersetzt. In diesem Fall wird in der Datei die richtige Last (1W pro Knoten) angezeigt:
            195.        1.0000000000
            196.        1.0000000000
            197.        1.0000000000
            usw.

In beiden Fällen erhalte ich die richtige Lösung.

Hat Jemand eine Idee, an was das liegen könnte? Wie komme ich zu den richtigen Informationen?

Vielen Dank im Voraus!

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

---

Du musst die Flächenlast zuerst mit
/SOLU
SFTRAN !funktioniert erst im Solution-Prozessor
auf Deine Knoten übertragen
Dann kannst Du die Lasten auslesen.
Beim Lösen werden die Lasten automatisch übertragen.

edit: Fehlerkorrektur

------------------
Viele Grüße         
Achim              PS: Als Anerkennung für einen Beitrag freu' nicht nur ich mich über Unities

Was ist der Unterschied zwischen Theorie und Praxis? Theorie ist, wenn man alles weiß und nichts funktioniert - Praxis ist, wenn alles funktioniert und keiner weiß warum (Autor unbekannt)

[Diese Nachricht wurde von Achim007 am 15. Mai. 2013 editiert.]

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

---

Hallo Achim007

Vielen Dank für den Tipp!

Ich habe die APDL-Befehle wie vorgeschlagen angepasst:

Code:

---

/prep7

nsel,S,,,M1                          ! Alle Knoten der Flaeche M1 auswaehlen
sf,all,hflux,10000                  ! 1W/cm2

allsel
/solu
sftran
! sbctran
nsel,S,,,M1

*cfopen,vekt,data                    ! Datei  erzeugen
*cfclos                              ! Datei schliessen

*get,nanz,node,,count                ! Ausgewaehlte Knoten zaehlen
*get,nmin,node,,num,min              ! Kleinste Knotennummer bestimmen

*do,j,1,nanz,1                      ! DO-Schleife ueber alle Knoten der Flaeche M1
nsel,s,node,,nmin                    ! Knoten mit der kleinsten Nummer auswaehlen
*get,nheat,node,nmin,f,heat          ! Heat des Knotens nmin bestimmen
*cfopen,vekt,data,,append            ! Datei oeffnen zum Daten Anheften
*vwrite,nmin,nheat                  ! Variablen in Datei schreiben
(F10.0,F20.10)                      ! Datenformat definieren
*cfclos                              ! Datei schliessen
nsel,S,,,M1                          ! Alle Knoten der Flaeche M1 erneut auswaehlen
*get,nnext,node,nmin,nxth            ! Naechste Knotennummer nach nmin bestimmen
nmin=nnext                          ! Naechste Knotennummer in nmin speichern
*enddo

---

Die erhaltenen Lasten sind allerdings immer noch Null. Keine Ahnung warum.

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

---

für die Unitis, obwohl ich Dir ja noch nicht wirklich weitergeholfen habe.
Welchen ET verwendest Du?

------------------
Viele Grüße         
Achim              PS: Als Anerkennung für einen Beitrag freu' nicht nur ich mich über Unities

Was ist der Unterschied zwischen Theorie und Praxis? Theorie ist, wenn man alles weiß und nichts funktioniert - Praxis ist, wenn alles funktioniert und keiner weiß warum (Autor unbekannt)

[Diese Nachricht wurde von Achim007 am 16. Mai. 2013 editiert.]

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

---

Die Unities sind für die konstruktive Antwort und für die Mühe, meine eher speziellen Fragen zu beantworten  .

Ich verwende SOLID90. Die Methode soll allerdings auch mit anderen ET funktionieren...

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

---

Output_Solid90.jpg

Du musst eine vorhandene Outputgröße auswählen (Siehe Anhang).

edit: habe bei der Strukturmechanik nachgeschaut, wahrscheinlich liege ich mit der Angabe doch nicht richtig.
Ich würde es mal mit einem Druck als Flächenlast versuchen und dann die Knotenkräfte ausgeben lassen
------------------
Viele Grüße         
Achim              PS: Als Anerkennung für einen Beitrag freu' nicht nur ich mich über Unities

Was ist der Unterschied zwischen Theorie und Praxis? Theorie ist, wenn man alles weiß und nichts funktioniert - Praxis ist, wenn alles funktioniert und keiner weiß warum (Autor unbekannt)

[Diese Nachricht wurde von Achim007 am 16. Mai. 2013 editiert.]

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

---

Vielen Dank für den Hinweis.
Ich werde weiter versuchen und mich ev. wieder melden.

Gruss

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

---

moin.
kann es sein, dass du einen Dreher in der Schleife hast? Mich macht stutzig, dass die Werte exakt NULL sind, numerische Lösungen geben i. allg. Werte nahe Null raus
Folgende Alternative:

Code:

---

NSEL,s,NODE,,M1                      ! M1 ist ein CM von dir
*get,nanz,node,,count                ! Ausgewaehlte Knoten zaehlen
*DIM,MY_DATA,ARRAY,NANZ,2
*DO,i_Node,1,NANZ,1
      *GET,MY_DATA(I_Node,1),NODE,,NUM,MIN       
      *GET,MY_DATA(I_Node,2),NODE,MY_DATA(I_Node,1),F,HEAT
      NSEL,U,NODE,,MY_DATA(I_Node,1)

*ENDDO
*VED,MY_DATA

---

Optional via *CFC, *VWR etc. die Daten rausschreiben. Versuch mal das Formalt (2E14.7), da du mit dem E-Format auch Werte sehr nahe Null darstellen kannst.

frohes schaffen
ral

EDIT: Korrektur

[Diese Nachricht wurde von RAL am 17. Mai. 2013 editiert.]

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