Announcement

Collapse
No announcement yet.

Abaqus UEL results

Collapse
X
  • Filter
  • Time
  • Show
Clear All
new posts

  • Abaqus UEL results

    I finished my User-Element, am in test/verify phase now. I was wondering, the results can be displayed via dummy elements, but the only thing you can see is the displacement, am I right? I also calculate the stresses in my own file, is it possible to plot them as well?

    And another question is, what for properties do you give the dummy element. I tried this and the results I could see where from the dummy element....

    Thanks for the support!

  • #2
    Hi,

    I have the same problem.

    http://www.polymerfem.com/forums/showthread.php?t=1577

    My own element is a simple flat shell, I need to plot SV field for my own elements.

    Best Regards

    Comment


    • #3
      However, for the dummy elements you need "dummy properties" as in the following example:

      ************************************************** ***** Dummy elements
      *Element, type=S4R, ELSET=E0
      1, 1, 2, 5, 4
      2, 2, 3, 6, 5
      3, 4, 5, 8, 7
      4, 5, 6, 9, 8
      5, 7, 8, 11, 10
      6, 8, 9, 12, 11
      7, 10, 11, 14, 13
      8, 11, 12, 15, 14
      ***Dummy section
      *Shell General Section, elset=E0, material=mat1
      1.0e-10,
      ***Dummy material
      *Material, name=mat1
      *Elastic
      1e-16, 0.3
      ***Real elements that are the same
      ***of dummy elements but with different id
      *USER ELEMENT, NODES=4, TYPE=U1, PROPERTIES=3, COORDINATES=3, VARIABLES=36
      1,2,3,4,5,6

      *Element, type=U1, ELSET=E1
      9, 1, 2, 5, 4
      10, 2, 3, 6, 5
      11, 4, 5, 8, 7
      12, 5, 6, 9, 8
      13, 7, 8, 11, 10
      14, 8, 9, 12, 11
      15, 10, 11, 14, 13
      16, 11, 12, 15, 14
      *UEL PROPERTY, ELSET=E1

      1e+06, 0.3, 0.1
      ************************************************** *

      then, using dummy property the corresponding stresses will be also dummy. Then you need, and me too, a way for to get plotted the stresses and the SV fields.

      Best regards

      Comment


      • #4
        Thanks ZaK,

        With field output and than SDV you can get your own SVARS as an output in the obd-file. However, when I use a dummy element, this data (from my UEL) is overwritten by the data from the dummy-element.

        ** FIELD OUTPUT
        **
        *Output, field, number interval=10, time marks=NO
        *Node Output
        RF, U, V, NT
        *Element Output, directions=YES
        HFL, LE, PE, PEEQ, S, SF, TEMP, SDV

        Maybe you know how to deal with this...

        Comment


        • #5
          Originally posted by schutter View Post
          Thanks ZaK,

          With field output and than SDV you can get your own SVARS as an output in the obd-file. However, when I use a dummy element, this data (from my UEL) is overwritten by the data from the dummy-element.
          Maybe, you can chose to get SDV for ELSET that contains your own elements and not for all.

          Originally posted by schutter View Post
          Thanks ZaK,

          ** FIELD OUTPUT
          **
          *Output, field, number interval=10, time marks=NO
          *Node Output
          RF, U, V, NT
          *Element Output, directions=YES
          HFL, LE, PE, PEEQ, S, SF, TEMP, SDV

          Maybe you know how to deal with this...
          It is right, but i don't need to plot SDV but I must use SVARS for to calculate the SF (Section force and moments). Then I need to get comunicating UEL and UGENS. In particular UGENS must use SVARS from UEL. But, I don't know yet how to do that.

          Thank you.

          Comment


          • #6
            Hi schutter,

            some news about UEL ??

            Comment


            • #7
              Hi ZaK!

              Yes, I've got the UEL working now. The only thing I need to do now is to understand and implement UVARM routine, cause I can plot the results of the SDV's but if there are more than 10 the script is becoming a real mess.

              As for your problem, if I understand your right, you need to get info from the SDV's to calculate your SF in UGENS? Is this correct?

              Comment


              • #8
                Yes, that is right. But, I would like to avoid the using of UGENS and also of dummy elements.
                Now, I can open ODB by c++ interface after analysis... But, I don't know yet how

                - to change element type from USER1 to S4 for example
                - to calculate SF field by using SDV field

                Best regards,
                zak

                Comment


                • #9
                  Dear all,

                  I am trying to output results as well from a UEL using SDV. I am trying to use the dummy element trick and I am quite confident that I understood how to do this. But when I run the simulation only my USER elements are displayed. When I tried first time I was able to see both elements overlapping, now no.

                  *Heading
                  **
                  *Node

                  ** Dummy elements
                  *Element, type=CPE4R, ELSET=E0
                  (NUMBERS)

                  **Dummy section
                  *Shell General Section, elset=E0, material=mat1
                  1.0e-10,
                  **Dummy material
                  *Material, name=mat1
                  *Elastic
                  1e-16, 0.3
                  *Depvar
                  20,
                  *User Defined Field

                  *user element, type=u1, nodes=4, coordinates=2, properties=11,
                  variables=320
                  1,2
                  1,11,12,13
                  *Element, type=u1, elset=user
                  (NUMBERS)


                  *uel property, elset=user
                  ** YMOD, Pois, sigma0, hard, DissipativelenScale(XL), energerticlenScale(smallL), thick, refstrainrate(eps0)
                  68380., 0.3, 2500., 437.34, 2., 0., 20., 1e-04,
                  **onh, eta, exponent(xn)
                  0.2, 0.1, 0.2
                  *Nset, nset=bottom, generate
                  1, 51, 1
                  *Elset, elset=bottom, generate
                  1, 50, 1
                  *Nset, nset=left, generate
                  1, 2551, 51
                  *Elset, elset=left, generate
                  1, 2451, 50
                  *Nset, nset=right, generate
                  51, 2601, 51
                  *Elset, elset=right, generate
                  50, 2500, 50
                  *Nset, nset=top, generate
                  2551, 2601, 1
                  *Elset, elset=top, generate
                  2451, 2500, 1
                  *Elset, elset=middle-elements
                  735, 736, 765, 766
                  *Nset, nset=middle-node
                  791

                  **
                  * INITIAL CONDITION, TYPE=FIELD
                  ** STEP: Step-1
                  **
                  *Step, name=Step-1, nlgeom=NO
                  *Static
                  5e-02, 1.0, 1e-10, 1e-01
                  **
                  ** BOUNDARY CONDITIONS
                  **
                  ** Name: BC-1 Type: Displacement/Rotation
                  *Boundary, TYPE=DISPLACEMENT
                  bottom, 1,2,0.0
                  top, 2, 2, 0.0
                  top, 1, 1, 10.
                  **
                  ** OUTPUT REQUESTS
                  **
                  *Restart, write, frequency=0
                  **
                  **
                  ** FIELD OUTPUT: F-Output-3
                  **
                  *Output, field
                  *Node Output
                  RF, U
                  *Element Output, elset=E0
                  MISES, S, SDV

                  *End Step


                  Comment

                  Working...
                  X