DKT element input data

[__Jules_]

Hello,

I am trying to use the DKT element which is provied http://projects.ce.berkeley.edu/feap/elmt07.f. I downdloaded the FEAPpv version and overwrote the elmt07.f dummy file. I used the following input file but somesthing seems to be wrong because the elmt1 subroutine is always called instead and I get the dummy output.

Code: [Select]

FEAP * * TwistedBeam_1x06_tri
  14,12,1,3,6,3
MATErial,1
USER 7
ELAStic 2.9E7 0.22 0.05 0.0
                  ! Blank termination record
COORdinates
  1  0 0 -0.55 0
  2  0 2 -0.531259204458988 -0.142350474806386
  3  0 4 -0.476313972081441 -0.275
  4  0 6 -0.388908729652601 -0.388908729652601
  5  0 8 -0.275 -0.476313972081441
  6  0 10 -0.142350474806386 -0.531259204458988
  7  0 12 0.0 -0.55
  8  0 0 0.55 0
  9  0 2 0.531259204458988 0.142350474806386
  10 0 4 0.476313972081441 0.275
  11 0 6 0.388908729652601 0.388908729652601
  12 0 8 0.275 0.476313972081441
  13 0 10 0.142350474806386 0.531259204458988
  14 0 12 0.0 0.55
                  ! Blank termination record

Code: [Select]

FEAP * * TwistedBeam_1x06_tri                                                 

    F I N I T E   E L E M E N T   A N A L Y S I S   P R O G R A M
              FEAPpv (P e r s o n a l   V e r s i o n)

             (C) Regents of the University of California
                       All Rights Reserved.

     Solution date: Thu Apr 19 10:47:00 2018

              VERSION: 4.1 Revision f   
              DATE: 20 September 2017

     Input Data Filename: ItwistedBeam_1x06_shell_tri                                                                                                     

     Number of Nodal Points  - - - - - - :       14
     Number of Elements  - - - - - - - - :       12

     Spatial Dimension of Mesh - - - - - :        3
     Degrees-of-Freedom/Node (Maximum) - :        6
     Equations/Element       (Maximum) - :        0
     Number Element Nodes    (Maximum) - :        3

     Number of Material Sets - - - - - - :        1
     Number Parameters/Set   (Program) - :      300
     Number Parameters/Set   (Users  ) - :      150

 FEAP * * TwistedBeam_1x06_tri                                                 

   M a t e r i a l    P r o p e r t i e s

     Material Number   1: Element Type: USER           : ELMT =  1
I would be happy if someone could give a hint.

[FEAP_Admin]

You should not overwrite dummy routines in FEAPpv.  The correct workflow is to build the FEAPpv library and link against it.  When linking against the library, if you provide a new version for a subroutine, the linker will use you provided routine instead of the one in the library.

That said you can see that element 7 is not being called so you need to double check that you have actually managed to compile it into your code.

[__Jules_]

I have in my solution the project of the lib so it is automatically updated. The Problem is that he is calling the subroutine for user element 1 and not for element 7 so I was thinking someting is wrong with my input file. The output also show that he somehow is not reading USER 7 from the input file.

[FEAP_Admin]

Can you add a check for isw.eq.0 to the logic block and then put a print state inside it to the effect of write(*,*) 'USER DTK Element'?  The rebuild your executable and type SHOW,ELEM at the macro prompt, to see if it is even in your code?

[FEAP_Admin]

Also, post the full output file.

[__Jules_]

I am very sure it is compiled. Anyway I added the check in elem7.f and an also in elem1.f. It is just showing the new output of elem1.f, because like I said somehow he is not reading the USER 7. The Programm things it should use USER 1. After the input file is proceeded i type SHOW,ELEM and the program crashed. I also tried this command with the executable provided by your website same happens. Here is the output:

Code: [Select]

FEAP * * TwistedBeam_1x06_tri                                                 

    F I N I T E   E L E M E N T   A N A L Y S I S   P R O G R A M
              FEAPpv (P e r s o n a l   V e r s i o n)

             (C) Regents of the University of California
                       All Rights Reserved.

     Solution date: Thu Apr 19 15:01:31 2018

              VERSION: 4.1 Revision f   
              DATE: 20 September 2017

     Input Data Filename: ItwistedBeam_1x06_shell_user7_lc2                                                                                               

     Number of Nodal Points  - - - - - - :       14
     Number of Elements  - - - - - - - - :       12

     Spatial Dimension of Mesh - - - - - :        3
     Degrees-of-Freedom/Node (Maximum) - :        6
     Equations/Element       (Maximum) - :        0
     Number Element Nodes    (Maximum) - :        3

     Number of Material Sets - - - - - - :        1
     Number Parameters/Set   (Program) - :      300
     Number Parameters/Set   (Users  ) - :      150

 FEAP * * TwistedBeam_1x06_tri                                                 

   M a t e r i a l    P r o p e r t i e s

     Material Number   1: Element Type: USER           : ELMT =  1
                               
     Degree of Freedom Assignments    Local    Global
                                     Number    Number
                                        1         1
                                        2         2
                                        3         3
                                        4         4
                                        5         5
                                        6         6

 FEAP * * TwistedBeam_1x06_tri                                                 

     Nodal Coordinates   

      node    1 Coord    2 Coord    3 Coord
         1  0.000E+00 -5.500E-01  0.000E+00
         2  2.000E+00 -5.313E-01 -1.424E-01
         3  4.000E+00 -4.763E-01 -2.750E-01
         4  6.000E+00 -3.889E-01 -3.889E-01
         5  8.000E+00 -2.750E-01 -4.763E-01
         6  1.000E+01 -1.424E-01 -5.313E-01
         7  1.200E+01  0.000E+00 -5.500E-01
         8  0.000E+00  5.500E-01  0.000E+00
         9  2.000E+00  5.313E-01  1.424E-01
        10  4.000E+00  4.763E-01  2.750E-01
        11  6.000E+00  3.889E-01  3.889E-01
        12  8.000E+00  2.750E-01  4.763E-01
        13  1.000E+01  1.424E-01  5.313E-01
        14  1.200E+01  0.000E+00  5.500E-01

 FEAP * * TwistedBeam_1x06_tri                                                 

     E l e m e n t s

   Elmt Mat Reg  Type  1 Node  2 Node  3 Node
      1   1   0 Trian       8       1       9
      2   1   0 Trian       2       9       1
      3   1   0 Trian       9       2      10
      4   1   0 Trian       3      10       2
      5   1   0 Trian      10       3      11
      6   1   0 Trian       4      11       3
      7   1   0 Trian      11       4      12
      8   1   0 Trian       5      12       4
      9   1   0 Trian      12       5      13
     10   1   0 Trian       6      13       5
     11   1   0 Trian      13       6      14
     12   1   0 Trian       7      14       6

 FEAP * * TwistedBeam_1x06_tri                                                 

  N o d a l   B. C.

      Node 1-b.c. 2-b.c. 3-b.c. 4-b.c. 5-b.c. 6-b.c.
         1      1      1      1      1      1      1
         8      1      1      1      1      1      1

 FEAP * * TwistedBeam_1x06_tri                                                 

     Nodal Forces       

      node    1 Force    2 Force    3 Force    4 Force    5 Force    6 Force

         7  0.000E+00  0.000E+00 -5.000E-01  0.000E+00  0.000E+00  0.000E+00
        14  0.000E+00  0.000E+00 -5.000E-01  0.000E+00  0.000E+00  0.000E+00

     E l e m e n t   S i z e   V a l u e s
          h-minimum =  1.1000E+00
          h-maximum =  2.2780E+00

     E q u a t i o n / P r o b l e m   S u m m a r y:

          Space dimension (ndm) =         3   Number dof (ndf) =       6
          Number of equations   =        72   Number nodes     =      14
          Average col. height   =        26   Number elements  =      12
          Number profile terms  =      1836   Number materials =       1


          Material    Element Tag   Element Type  History Terms  Element Terms
              1              1              1              0              0

 FEAP * * TwistedBeam_1x06_tri                                                 

  M a c r o   I n s t r u c t i o n s

  Macro Statement            Variable 1  Variable 2  Variable 3
       FORM residual         0.0000E+00  0.0000E+00  0.0000E+00
       TANG                  0.0000E+00  0.0000E+00  0.0000E+00
       SOLV                  0.0000E+00  0.0000E+00  0.0000E+00
       DISP ALL              0.0000E+00  0.0000E+00  0.0000E+00
       STRE ALL              0.0000E+00  0.0000E+00  0.0000E+00
       PLOT STRE             1.0000E+00  0.0000E+00  0.0000E+00
       PLOT REACtions        0.0000E+00  0.0000E+00  0.0000E+00
 *Macro   1 * form residual       v:   0.00       0.00       0.00   
                                                           t=     0.01     0.00
   Residual norm =     7.0710678E-01    1.0000000E+00
 *Macro   2 * tang                v:   0.00       0.00       0.00   
                                                           t=     0.02     0.00
                                                           t=     0.05     0.00
 *WARNING* reduced diagonal is zero in   71 equations.
 Step:      0 Iteration:     0
 Condition check: D-max 0.0000E+00; D-min 0.0000E+00; Ratio 0.0000E+00
 Maximum no. diagonal digits lost:  1
   End Triangular Decomposition                            t=     0.05     0.00
 *Macro   3 * solv                v:   0.00       0.00       0.00   
                                                           t=     0.05     0.00
   Energy convergence test
    Maximum   =    0.000000000000000E+00 Current   =    0.000000000000000E+00
    Relative  =    1.000000000000000E+00 Tolerance =    1.000000000000000E-16
 *Macro   4 * disp ALL            v:   0.00       0.00       0.00   
                                                           t=     0.05     0.00

 FEAP * * TwistedBeam_1x06_tri                                                 

  N o d a l   D i s p l a c e m e n t s     Time       0.00000E+00
                                            Prop. Ld.  1.00000E+00

   Node     1 Coord     2 Coord     3 Coord
            1 Displ     2 Displ     3 Displ     4 Displ     5 Displ     6 Displ
      1  0.0000E+00 -5.5000E-01  0.0000E+00
         0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
      2  2.0000E+00 -5.3126E-01 -1.4235E-01
         0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
      3  4.0000E+00 -4.7631E-01 -2.7500E-01
         0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
      4  6.0000E+00 -3.8891E-01 -3.8891E-01
         0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
      5  8.0000E+00 -2.7500E-01 -4.7631E-01
         0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
      6  1.0000E+01 -1.4235E-01 -5.3126E-01
         0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
      7  1.2000E+01  0.0000E+00 -5.5000E-01
         0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
      8  0.0000E+00  5.5000E-01  0.0000E+00
         0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
      9  2.0000E+00  5.3126E-01  1.4235E-01
         0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
     10  4.0000E+00  4.7631E-01  2.7500E-01
         0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
     11  6.0000E+00  3.8891E-01  3.8891E-01
         0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
     12  8.0000E+00  2.7500E-01  4.7631E-01
         0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
     13  1.0000E+01  1.4235E-01  5.3126E-01
         0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
     14  1.2000E+01  0.0000E+00  5.5000E-01
         0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00  0.0000E+00
 *Macro   5 * stre ALL            v:   0.00       0.00       0.00   
                                                           t=     0.05     0.00
 *Macro   6 * plot STRE           v:   1.00       0.00       0.00   
                                                           t=     0.07     0.00
    Minimum is   0.00E+00 Maximum is   0.00E+00
    Minimum is  -1.00E-08 Maximum is   1.00E-08
    Minimum is  -1.00E-08 Maximum is   1.00E-08
   ------ Contour Values for Plot ------
     -0.714286E-08  -0.428571E-08  -0.142857E-08   0.142857E-08   0.428571E-08
      0.714286E-08
 ** ERROR ** No plot - zero difference in values
 *Macro   7 * plot REACtions      v:   0.00       0.00       0.00   
                                                           t=     0.16     0.00
 *End of Macro Execution*                                  t=     0.19     0.00
 *Macro   1 * show elem           v:   0.00       0.00       0.00   
                                                           t=     8.87     0.00



The line in the output shows that he wants to use USER element 1 instead of 7

Code: [Select]

Material Number   1: Element Type: USER           : ELMT =  1

[Prof. R.L. Taylor]

I suggest you download a clean version of feappv and make changes to the vinput.f file in the  ./program directory to handle the problem of character length issues.

At line 99 change no:n to no:n-1 in the call to setval
At line 111 make the same change.

Hopefully this will also correct the issue of which element is called for the DKT

[__Jules_]

I did as recommended. Unfortunately this did not solved the problem. I found out that in pmatin.f in line 106 the material number instead of the user element number is given to vinput(). For now I just made a little workaround because I dont know in which variable the real user number is saved.

Code: [Select]

if(pcomp(tx(1),'user',4)) then ! Get name as 'user ma'
WRITE(tx(2), '(A)') '7'
        errck = vinput(tx(2),15,uelnum,1)
        iel   = nint(uelnum(1))
      else

[Prof. R.L. Taylor]

The data  in pmatin.f should read:

!     Input record for element type selection

301   if(ior.lt.0) write(*,2004)
      il = min(ndf+2,14)
      errck = tinput(tx(1),2,td,il)
      if(errck) go to 301

so that 16 items are read and the first two are text values

[Prof. R.L. Taylor]

there are several other problems with the pmatin.f -- replace with the attach file.

[FEAP_Admin]

I've patched the GITHUB repository (now at version 4.1i).  This version also has quick fix such that SHOW,ELEM works (assuming that you have programed output for isw.eq.0).

[__Jules_]

this is working. Great