---

Hallo,

ich kämpfe momentan mit meinen User Elements. Und zwar funktioniert die Uebergabe der Properties aus dem Inputfile an die Subroutine nicht richtig.

Inputfile:

...
*USER ELEMENT, TYPE=U1, NODES=2, PROPERTIES=3, COORDINATES=3
1,2,3
*Element, type=U1
1, 1, 2
2, 2, 3
...
*UEL PROPERTY, ELSET=USERELEMENTS
300.0, 0.3, 10.0
...

Wenn ich die Werte in der UEL ausgeben lasse via

...
WRITE(6,'(1X,A,F9.3,A,F9.3,A,F9.3,A,I3)')
    1 " :: ",PROPS(1)," ",PROPS(2)," ",PROPS(3)," ",NPROPS
...

bekomme ich folgende Ausgabe

:: 0.000    3.793    0.000  3

Die Rechnung bricht dann mit einer floating point exception ab, weil der dritte Wert eine Querschnittsfläche ist durch die dividiert werden soll...

Kann mir jemand schreiben wo der Fehler liegt? Nach der Doku dachte ich dass die Variablen über UEL PROPERTY übergeben und dann im PROPS array verfügbar sind...

Danke für die Hilfe!

Gruss
Chris

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

---

Hallo Chris,

das ist erstmal nur ne Vermutung, aber hast du auch das Elementset "USERELEMENTS" vor *UEL PROPERTY erstellt?

Ich mach das bei meinen Elementen noch im Part-Bereich der *.inp Datei.

mfg NxxN

------------------
==========
==  Dingsen  ==
==========

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

---

Hallo Nicksen,

Danke für die Antwort!

Der Auszug aus dem Inputfile ist in der Part Section. Das Elementset USERELEMENTS wird vorher erstellt mit

*Elset, elset=USERELEMENTS, internal
1, 2

Daran sollte es also eigentlich nicht liegen...

Grüsse
Chris

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

---

Hi Chris,

...okay dann ist das schon einmal ausgeschlossen.

   

Bei der UEL hast du aber das komplette Interface übernommen?!
Ich meine damit:
SUBROUTINE UEL(RHS,AMATRX,SVARS,ENERGY,NDOFEL,NRHS,NSVARS, ... usw.
und
DIMENSION RHS(MLVARX,*),AMATRX(NDOFEL,NDOFEL), ... usw.

------------------
==========
==  Dingsen  ==
==========

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

---

Hi,

den Header hab ich aus der Abaqus Doku rauskopiert:

      SUBROUTINE UEL(RHS,AMATRX,SVARS,ENERGY,NDOFEL,NRHS,NSVARS,
    1 PROPS,NPROPS,COORDS,MCRD,NNODE,U,DU,V,A,JTYPE,TIME,DTIME,
    2 KSTEP,KINC,JELEM,PARAMS,NDLOAD,JDLTYP,ADLMAG,PREDEF,NPREDF,
    3 LFLAGS,MLVARX,DDLMAG,MDLOAD,PNEWDT,JPROPS,NJPROP,PERIOD)

      INCLUDE 'ABA_PARAM.INC'     
     
      DIMENSION RHS(MLVARX,*),AMATRX(NDOFEL,NDOFEL),PROPS(*),
    1 SVARS(*),ENERGY(8),COORDS(MCRD,NNODE),U(NDOFEL),
    2 DU(MLVARX,*),V(NDOFEL),A(NDOFEL),TIME(2),PARAMS(*),
    3 JDLTYP(MDLOAD,*),ADLMAG(MDLOAD,*),DDLMAG(MDLOAD,*),
    4 PREDEF(2,NPREDF,NNODE),LFLAGS(*),JPROPS(*)

von daher sollte der komplett und richtig sein.

Gruss
Chris

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

---

Wenns dir nicht zu geheim ist, dann schick mal die inp und die subroutine.
Geht evtl schneller, wenn ich mal etwas rumspiele.

------------------
==========
==  Dingsen  ==
==========

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

---

Hallo Nicksen,

Danke für das Angebot! Das ganze Modell schick ich ungern rum, ich werd mal ein kleineres Beispiel für die Subroutine basteln...

Gruss
Chris

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

---

Das waere gut. Es reichen ja 3 Staebe, wenn man es auf zwei Dimensionen reduziert.
Hauptsache das inp hat zumindest die Angaben drin, welche im großen Modell nicht funktionieren. Bzw. solltest du das kleine Modell auch mal versuchen und wenn der Fehler dann wieder kommt, dann arbeiten wir dran.

mfg

------------------
==========
==  Dingsen  ==
==========

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

---

Hallo Nicksen,

als ich gestern Nachmittag das Minimal Beispiel zusammengebastel hab, hats natürlich funktioniert.

Verändert hab ich folgendes:

In der UEL Implizite Variablendeklaration mit IMPLICIT NONE deaktiviert.

Im INPUTFILE die Parameter als Float angegeben (300 -> 300.0)

und scheinbar noch irgendwas, weil ich den Fehler grad nicht mehr reproduzieren kann...

Danke für die Antworten!

Gruss Chris

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

---

Ach was, ich hab doch garnichts gemacht 

Viel Erfolg weiterhin!!!

und ein schoenes Wochenende

mfg

------------------
==========
==  Dingsen  ==
==========

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