---

Hallo Community,

bin sowohl in diesem Forum, als auch auf dem Gebiet Subroutinen ein Newbie und würde mich deshalb über eure Hilfe freuen.

Mittels einer Subroutine soll die Kraftverteilung durch eine Strömung auf einen Balken berechnet werden. Input ist dabei die Position einzelner Knotenpunkte auf dem Balken. Der Balken wird also in Panels unterteilt. Diese Subroutine soll dann von ABAQUS genutzt werden. die Subrotine DLOAD scheint für mein vorhaben am geeignetsten zu sein. Problem ist jetzt, ich hab keinen Schimmer wie das ganze in Fortran geschrieben wird.
Hat sowas ähnliches hier jmd schonmal gemacht bzw. könnt ihr mir Tipps, Anregungen etc. pp. geben?

Danke und Gruß

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

---

Im Handbuch gibt es mehrere Beispiele. Schau z.B. mal ins Verification Manual.

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

---

Da habe ich schon reingeschaut werde aber ehrlich gesagt noch nicht schlau daraus.

Zitat:

---

Im Handbuch gibt es mehrere Beispiele.

---

Kannst du vllt einen Link für ein Bsp posten? Sorry für dieses babysitting aber dieses Manual ist momentan noch ein echtes Labyrinth für mich.

Danke und Gruß

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

---

Im Verification Manual sind im Kapitel "4.1.3 DLOAD" 4 Beispiele verlinkt. Da drauf klicken und dich mit den Themen befassen musst du schon selbst.

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

---

komme leider in der Thematik immer noch nicht weiter. Habe jetzt versucht eine "einfache" DLOAD subroutine zu implementieren, es wird auch eine .obj-Datei generiert die ich dann in ABAQUS einbinden will, aber am Ende kommt immer.

ERROR: USER SUBROUTINE DLOAD IS MISSING

Ist das zurückzuführen auf einen Fehler im Sourcecode selbst oder kann die von mir angegebene Routine tatsächlich nicht gefunden werden?

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

---

Damit wird eindeutig die Subroutine nicht gefunden.

Entweder weil die fortran-Datei SELBST nicht gefunden wird oder die subroutine IN der fortran-Datei nicht gefunden wird.

Sieh in der log-Datei des Jobs nach. Dort sollte beim Kompilieren einiges reingeschrieben werden. Vielleicht hilft das weiter, ansonsten poste mal den Inhalt der log-Datei?

------------------
--
nix zu sehen

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

---

Hier ist die vollständige log-Datei. Ich bin absoluter Anfänger auf dem Gebiet Abaqus/Subroutinen usw. wenn also jemand essentielle Infos aus dem Dokument lesen kann außer besagtem 'ERROR', wäre ich wirklich sehr dankbar.

STEP    1    INCREMENT    1    STEP TIME    0.00   

                        S T E P      1    D Y N A M I C  A N A L Y S I S - I M P L I C I T  I N T E G R A T I O N

                                                                                         

          DYNAMIC APPLICATION PACKED SCHEME IS                TRANSIENT FIDELITY 
          HILBER-HUGHES-TAYLOR (HHT) TIME INTEGRATOR IS USED

    AUTOMATIC TIME CONTROL WITH -
          A SUGGESTED INITIAL TIME INCREMENT OF                10.0   
          AND A TOTAL TIME PERIOD OF                            10.0   
          THE MINIMUM TIME INCREMENT ALLOWED IS                1.000E-04
          THE MAXIMUM TIME INCREMENT ALLOWED IS                10.0   
          SCALE FACTOR ON TIME AVERAGE FORCE FOR HALF-
          -INCREMENT RESIDUAL COMPUTATIONS IS                1.000E+03

          HILBER-HUGHES-TAYLOR METHOD PARAMETERS:  ALPHA  IS  -5.000000E-02
                                                  BETA    IS  0.275625   
                                                  GAMMA  IS  0.550000   
          SINGULAR MASS TREATMENT TYPE IS                    ERROR
          CONSERVATIVE TIME INCREMENTATION IS USED
          STEP AMPLITUDE IS                                    STEP

CONVERGENCE TOLERANCE PARAMETERS FOR FORCE   
    CRITERION FOR RESIDUAL FORCE    FOR A NONLINEAR PROBLEM          5.000E-03
    CRITERION FOR DISP.    CORRECTION IN A NONLINEAR PROBLEM          1.000E-02
    INITIAL VALUE OF TIME AVERAGE FORCE                              1.000E-02
    AVERAGE FORCE    IS TIME AVERAGE FORCE   
    ALTERNATE CRIT. FOR RESIDUAL FORCE    FOR A NONLINEAR PROBLEM    2.000E-02
    CRITERION FOR ZERO FORCE    RELATIVE TO TIME AVRG. FORCE        1.000E-05
    CRITERION FOR RESIDUAL FORCE    WHEN THERE IS ZERO FLUX          1.000E-05
    CRITERION FOR DISP.    CORRECTION WHEN THERE IS ZERO FLUX        1.000E-03
    CRITERION FOR RESIDUAL FORCE    FOR A LINEAR INCREMENT          1.000E-08
    FIELD CONVERSION RATIO                                            1.00   
    CRITERION FOR ZERO FORCE    REL. TO TIME AVRG. MAX. FORCE        1.000E-05
    CRITERION FOR ZERO DISP.    RELATIVE TO CHARACTERISTIC LENGTH    1.000E-08

    VOLUMETRIC STRAIN COMPATIBILITY TOLERANCE FOR HYBRID SOLIDS      1.000E-05
    AXIAL STRAIN COMPATIBILITY TOLERANCE FOR HYBRID BEAMS            1.000E-05
    TRANS. SHEAR STRAIN COMPATIBILITY TOLERANCE FOR HYBRID BEAMS      1.000E-05
    SOFT CONTACT CONSTRAINT COMPATIBILITY TOLERANCE FOR P>P0          5.000E-03
    SOFT CONTACT CONSTRAINT COMPATIBILITY TOLERANCE FOR P=0.0        0.100   
    CONTACT FORCE ERROR TOLERANCE FOR CONVERT SDI=YES                1.00   
    DISPLACEMENT COMPATIBILITY TOLERANCE FOR DCOUP ELEMENTS          1.000E-05
    ROTATION COMPATIBILITY TOLERANCE FOR DCOUP ELEMENTS              1.000E-05

EQUILIBRIUM WILL BE CHECKED FOR SEVERE DISCONTINUITY ITERATIONS

TIME INCREMENTATION CONTROL PARAMETERS:
    FIRST EQUILIBRIUM ITERATION FOR CONSECUTIVE DIVERGENCE CHECK              4
    EQUILIBRIUM ITERATION AT WHICH LOG. CONVERGENCE RATE CHECK BEGINS        8
    EQUILIBRIUM ITERATION AFTER WHICH ALTERNATE RESIDUAL IS USED              9
    MAXIMUM EQUILIBRIUM ITERATIONS ALLOWED                                  16
    EQUILIBRIUM ITERATION COUNT FOR CUT-BACK IN NEXT INCREMENT              10
    MAXIMUM EQUILIB. ITERS IN TWO INCREMENTS FOR TIME INCREMENT INCREASE      4
    MAXIMUM ITERATIONS FOR SEVERE DISCONTINUITIES                            50
    MAXIMUM CUT-BACKS ALLOWED IN AN INCREMENT                                5
    MAXIMUM DISCON. ITERS IN TWO INCREMENTS FOR TIME INCREMENT INCREASE      50
    MAXIMUM CONTACT AUGMENTATIONS FOR *SURFACE BEHAVIOR,AUGMENTED LAGRANGE  50
    CONSECUTIVE INCS MEETING TIME INTEG. TOL. FOR TIME INC. INCREASE          3
    CUT-BACK FACTOR AFTER DIVERGENCE                                0.2500   
    CUT-BACK FACTOR FOR TOO SLOW CONVERGENCE                        0.5000   
    CUT-BACK FACTOR AFTER TOO MANY EQUILIBRIUM ITERATIONS            0.7500   
    CUT-BACK FACTOR AFTER TOO MANY SEVERE DISCONTINUITY ITERATIONS  0.2500   
    CUT-BACK FACTOR AFTER PROBLEMS IN ELEMENT ASSEMBLY              0.2500   
    INCREASE FACTOR AFTER TWO INCREMENTS THAT CONVERGE QUICKLY        1.500   
    MAX. TIME INCREMENT INCREASE FACTOR ALLOWED                      1.500   
    MAX. TIME INCREMENT INCREASE FACTOR ALLOWED (DYNAMICS)            1.250   
    MAX. TIME INCREMENT INCREASE FACTOR ALLOWED (DIFFUSION)          2.000   
    MINIMUM TIME INCREMENT RATIO FOR EXTRAPOLATION TO OCCUR          0.1000   
    CUT-BACK FACTOR WHEN TIME INTEGRATION ACCURACY TOL. EXCEEDED    0.8500   
    MAX. RATIO OF TIME INTEG. MEASURE TO TOL. FOR TIME INC. INCREASE 0.7500   
    INCREASE FACTOR FOR TIME INCREMENT WHEN POSSIBLE                0.8000   
    MIN. TIME INC. INCREASE FACTOR IN LINEAR TRANSIENT PROBLEMS      0.9500   
    MAX. RATIO OF TIME INCREMENT TO STABILITY LIMIT                  1.000   
    FRACTION OF STABILITY LIMIT FOR NEW TIME INCREMENT              0.9500   
    TIME INCREMENT INCREASE FACTOR BEFORE A TIME POINT                1.000   
    GLOBAL STABILIZATION CONTROL IS NOT USED

          PRINT OF INCREMENT NUMBER, TIME, ETC., EVERY    1  INCREMENTS

    THE MAXIMUM NUMBER OF INCREMENTS IN THIS STEP IS                    100

    VELOCITY-BASED PARABOLIC EXTRAPOLATION WILL BE USED

    CHARACTERISTIC ELEMENT LENGTH      6.33   

    DETAILS REGARDING ACTUAL SOLUTION WAVEFRONT REQUESTED

    DETAILED OUTPUT OF DIAGNOSTICS TO DATABASE REQUESTED

    PRINT OF INCREMENT NUMBER, TIME, ETC., TO THE MESSAGE FILE EVERY    1  INCREMENTS

    COLLECTING MODEL CONSTRAINT INFORMATION FOR OVERCONSTRAINT CHECKS

    COLLECTING STEP CONSTRAINT INFORMATION FOR OVERCONSTRAINT CHECKS

  NUMBER OF EQUATIONS =      88    NUMBER OF FLOATING PT. OPERATIONS =  1.96E+04

    CHECK POINT  START OF SOLVER   

    CHECK POINT  END OF SOLVER     

      ELAPSED USER TIME (SEC)      =  0.0000   
      ELAPSED SYSTEM TIME (SEC)    =  0.0000   
      ELAPSED TOTAL CPU TIME (SEC) =  0.0000   
      ELAPSED WALLCLOCK TIME (SEC) =          0

***ERROR: USER SUBROUTINE DLOAD IS MISSING

    ANALYSIS SUMMARY:
    TOTAL OF          0  INCREMENTS
                      0  CUTBACKS IN AUTOMATIC INCREMENTATION
                      0  ITERATIONS INCLUDING CONTACT ITERATIONS IF PRESENT
                      1  PASSES THROUGH THE EQUATION SOLVER OF WHICH
                      1  INVOLVE MATRIX DECOMPOSITION, INCLUDING
                      1  DECOMPOSITION(S) OF THE MASS MATRIX
                      1  REORDERING OF EQUATIONS TO MINIMIZE WAVEFRONT
                      0  ADDITIONAL RESIDUAL EVALUATIONS FOR LINE SEARCHES
                      0  ADDITIONAL OPERATOR EVALUATIONS FOR LINE SEARCHES
                      0  WARNING MESSAGES DURING USER INPUT PROCESSING
                      0  WARNING MESSAGES DURING ANALYSIS
                      0  ANALYSIS WARNINGS ARE NUMERICAL PROBLEM MESSAGES
                      0  ANALYSIS WARNINGS ARE NEGATIVE EIGENVALUE MESSAGES
                      1  ERROR MESSAGES

    JOB TIME SUMMARY
      USER TIME (SEC)      =  0.10000   
      SYSTEM TIME (SEC)    =  0.10000   
      TOTAL CPU TIME (SEC) =  0.20000   
      WALLCLOCK TIME (SEC) =          1

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

---

Also daraus lässt sich auch nicht mehr herauslesen als zuvor.

Dies ist aber auch die message (*.msg) Datei und nicht die log-Datei (Endung *.log).
Gibt es keine log-Datei?

------------------
--
nix zu sehen

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

---

Sorry, doch gibts:

DLOAD.for(2): error #5149: Illegal character in statement label field  [S]
Subroutine DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,COORDS,
-^
DLOAD.for(2): error #5149: Illegal character in statement label field  [u]
Subroutine DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,COORDS,
--^
DLOAD.for(2): error #5149: Illegal character in statement label field  [b]
Subroutine DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,COORDS,
---^
DLOAD.for(2): error #5149: Illegal character in statement label field  [r]
Subroutine DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,COORDS,
----^
DLOAD.for(2): error #5118: First statement in file must not be continued
Subroutine DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,COORDS,
-----^
DLOAD.for(3): error #5082: Syntax error, found END-OF-STATEMENT when expecting one of: ( % [ . = =>
    1                JLTYP,SNAME)
-----------------------------------^
DLOAD.for(12): error #5276: Unbalanced parentheses
      F=TEN*(ONE-(EXP(CONST*TIME(1)))))
--------------------------------------^
DLOAD.for(12): error #5082: Syntax error, found ')' when expecting one of: <END-OF-STATEMENT> ;
      F=TEN*(ONE-(EXP(CONST*TIME(1)))))
--------------------------------------^
DLOAD.for(14): error #6353: A RETURN statement is invalid in the main program.
      RETURN
------^
compilation aborted for DLOAD.for (code 1)
Abaqus Error: Problem during compilation - C:\Temp\DLOAD.for
Abaqus/Analysis exited with errors

Die Subroutine, die ich testweise nutzen wollte ist:

Subroutine DLOAD(F,KSTEP,KINC,TIME,NOEL,NPT,LAYER,KSPT,COORDS,
    1                JLTYP,SNAME)
C
C    EXPONENTIAL PRESSURE LOAD
C
      INCLUDE 'ABA_PARAM.INC'
C
      DIMENSION COORDS(3),TIME(2)
      CHARACTER*80 SNAME
      DATA TEN,ONE,CONST /10.,1.,-230.3/
      F=TEN*(ONE-(EXP(CONST*TIME(1)))))
      IF(NPT.EQ.1) WRITE(6,*) ' LOAD APPLIED',F, 'AT TIME=',TIME(1)
      RETURN
      END

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

---

ok ich habs in fortran nochmal neu geschrieben und jetzt spuckt er keine Fehler mehr aus, die Subroutine funktioniert auch.
Ist denn beim schreiben des obigen Beispielcodes irgendwas zu beachten? Also zB das diese vier 'C' ganz vorne stehen, oder die '1' unter 'Subroutine'. Klingt albern die Frage ich weiß, aber ich hoffe dadurch ein etwas besseres Verständnis von Fortran zu kriegen.

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

---

Ja, da gibt es definitiv einiges zu beachten, ein kleines fortran-Handbuch oder Tutorial wuerde sicher nicht schaden.

Wikipedia sagt zB folgendes:

Zitat:

---

A letter "C" in column 1 caused the entire card to be treated as a comment and ignored by the compiler. Otherwise, the card was divided into four fields.
Columns 1 to 5 were the label field: a sequence of digits here was taken as a label for the purpose of a GOTO or a FORMAT reference in a WRITE or READ statement.
Column 6 was a continuation field: a non-blank character here caused the card to be taken as a continuation of the statement on the previous card.
Columns 7 to 72 served as the statement field.
Columns 73 to 80 were ignored, so they could be used for identification information.

---

https://en.wikipedia.org/wiki/Fortran#Fixed_layout_and_punched_cards

Einige dieser Dinge sind obsolet (also _muessen_ nicht mehr beachtet werden). Es kann aber leicht sein, dass du mit einem älteren fortran-Standard kompilierst und dort sind die Spalten eben SEHR wichtig.

lg!

------------------
--
nix zu sehen

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

---

Bin leider in meinem Problem immer noch nicht weitergekommen. Kann zwar wie gesagt eine einfache zeitabhängige DLOAD Routine benutzen die auch funktioniert, aber ich würde gerne eine positionsabhängige Last einführen. Mein Balken soll am Ende durch diese Belastung anfangen zu schwingen. Leider finde ich weder ein ansatzweise ähnliches template, noch finde ich den Fortran-Code um die Positionsänderung auszulesen.

Wenn jemand jemanden kennt der jemanden kennt...ich wäre euch sehr verbunden

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

---

Wo ist genau das Problem? Die aktuelle Koordinate des Lastintegrationspunktes wird doch im Feld COORDS übergeben (NLGEOM=YES).

------------------
Pam

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

---

Das Problem ist, dass ich das in der Subroutine als Schleife schreiben will: die neuen COORDS verändern die Kraft worauf die Integrationspunkte an eine neue Position verschoben werden die wieder von der Subroutine verarbeitet wird, die Kraft verändert sich wieder and so on.

Ist für Fortgeschrittene wohl machbar aber mir bereitet es kopfschmerzen :/

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

---

Wieso Schleife? Die Subroutine wird für jeden Lastpunkt in jeder Iteration durchlaufen. Dabei bekommst du diverse Informationen (Zeit, Ort, Nummern, ...) geliefert. In Abhängigkeit davon kannst du deine Last definieren.

Wenn du möchtest, dass eine leichte Änderung der Position schon ein Update der Last erzwingen soll, brauchst du doch nur die Inkrementierung entsprechend fein einstellen. Dreh also die Max. Increment Size runter. So wird die Subroutine zwangsweise oft durchlaufen und die Daten aktualisiert.

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

---

also ich hab mich mal dran versucht aber whs hab ich momentan nur "mist-code" geschrieben:

C
C User subroutine VDLOAD
      subroutine vdload (
C Read only -
     *     nblock, ndim, stepTime, totalTime, amplitude,
     *     curCoords, velocity, dircos, jltyp, sname,
C Write only -
     *     value )
C
      include 'vaba_param.inc'
      parameter ( pi=3.14)
C
      dimension curCoords(nblock,ndim), velocity(nblock,ndim),
     *     dircos(nblock,ndim,ndim), value(nblock)
*
      character*80 sname
 
*     do k = 1, nblock
      
 
  dX=  (curCOORDS(k,1))
  dY=  (curCOORDS(k,2))
 
 
X_quadrat= dX*dX
Y_quadrat= dY*dY



R= sqrt(X_quadrat + Y_quadrat)  

      value (k)= 1000*(1/pi*((R)-))
 
  end
   
in Abaqus Standard kriege ich mit ner DLOAD funktion keine Bewegung, sondern vielmehr verformt sich mein (zylindrischer) Balken. Eigentlich wollte ich, dass die ungleichmäßige Kraftverteilung den Balken hin und her bewegt aber sie scheint gleichmäßig zu sein?

Und wenn ich wie oben die VDLOAD Routine benutze, sehe ich zwar das Spannugen in meinem Balken vorliegen, bewegen tut er sich aber ebenso wenig.

Was mach ich da nur für blöde Fehler?  

[Diese Nachricht wurde von ye43leso am 19. Mrz. 2015 editiert.]

[Diese Nachricht wurde von ye43leso am 19. Mrz. 2015 editiert.]

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

---

Wenn ich den Code richtig verstehe, ist deine Last radial wachsend (oder fallend). Das bringt doch an der Stirnseite eines 3D Balken nur axiale Lasten, womit sich kein Grund für eine Biegeschwingung ergibt.

'nur eine Vermutung. Ich kenne ja dein Modell nicht...

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

---

Unbenannt.jpg

ja da hast du schon recht, aber ich lasse die kraft eigentlich flächenmäßig überall angreifen;( Soll den Druck eines umströmenden Fluids darstellen)
Hab mal ein Bild in den Anhang gepackt vll ist der Sinn meines Modells ersichtlicher.
Ich glaube selbst nicht das der Code so stimmt aber momentan finde ich keinen gescheiten Ansatz damit ein vernünftiger Schwinger herauskommt.

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

---

Wenn das Bauteil rotationssymmetrisch ist und die Achse auch bei x=y=0 liegt, dann hebt sich der radiale Druck auf und es entsteht keine Biegung oder radial ungleichmäßige Verformung. Somit gibt es keinen Grund für eine Schwingung, oder?

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

---

ja wohl wahr. Im Umkehrschluss müsste also mit nem Offset vom Nullpunkt oder ner Einzelkrafteinleitung was passieren?

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

---

Hallo zusammen,

ich bin auf diesen Beitrag aufmerksam geworden, da ich ein ähnliches Problem habe und auch noch am Anfang stehe.
Ich möchte innerhalb einer Subroutine VDLOAD koordinatenabhängig eine Last aufbringen. Sprich sobald ein Knoten eine Bewegung erfahren hat und eine bestimmte Position erreicht, soll die Last angreifen. Dafür müsste ich also die Knotenkoordinaten "alt" und "neu" innerhalb der Berechnung für jeden nblock-Satz überprüfen, ob eine Bewegung vorliegt und dann, ob die vorab festgelegte Position für die Lasteinleitung erreicht ist.

Das Problem klingt eigtl. ziemlich simpel, ich bin mir aber schon bei Kontrolle der Positionsänderung nicht sicher, dass die richtigen Einträge miteinander verglichchen werden. Ich wäre für jeden Tipp dankbar!

Danke vorab!

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