FEAP User Forum
FEAP => FEAPpv => Topic started by: __Jules_ on April 19, 2018, 02:30:27 AM
-
Hello,
I am trying to use the DKT element which is provied http://projects.ce.berkeley.edu/feap/elmt07.f (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.
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
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.
-
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.
-
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.
-
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?
-
Also, post the full output file.
-
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:
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
Material Number 1: Element Type: USER : ELMT = 1
-
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
-
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.
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
-
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
-
there are several other problems with the pmatin.f -- replace with the attach file.
-
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).
-
this is working. Great :)