simple text output

Fragen zu Solvertypen, Multicore-Rechnungen, Spannungsparametern /
Issues to solvers, multi core, stress parameter

Moderatoren: ccad, mz15, auroraIco, Lehrstuhl

Antworten
Dexter22
Newbie
Beiträge: 8
Registriert: Fr 10. Okt 2014, 20:41

simple text output

Beitrag von Dexter22 »

I've setup a simple simulation of a beam element No. 25,
1 element,
2 nodes,
Constrained at node 1 in all 6 dof
Under a force in the -ve z axis at node 2

I solve and get:

"Z88ER done
linear solver Z88R done
Solver successfull"

in z88r.log

displacements file with odd entries:


output file Z88O2.TXT: displacements, computed by Z88R V15.0B
*************
Loadcase: 1
node U(1) U(2) U(3) U(4) U(5) U(6)

1 -1.#IND000E+000 -1.#IND000E+000 -1.#IND000E+000 -1.#IND000E+000 -1.#IND000E+000 -1.#IND000E+000
2 -1.#IND000E+000 -1.#IND000E+000 -1.#IND000E+000 -1.#IND000E+000 -1.#IND000E+000 -1.#IND000E+000

Firstly, what does IND000E+000 mean?

Secondly, does anyone have any suggestions to fix this unreported error?
Benutzeravatar
SHautsch
Alumni
Beiträge: 380
Registriert: Mo 15. Apr 2013, 11:03

Re: simple text output

Beitrag von SHautsch »

Hello Dexter22,

-1.#IND000E+000 means that this value is indeterminate. Windows says "-1.#IND" while Unix says "NaN" which means "not a number".

Which solver did you use? Did you set all parameters correctly? Could you provide us with your project files, so we can have a deeper look at your problem?

Kind regards,
SHautsch
Dexter22
Newbie
Beiträge: 8
Registriert: Fr 10. Okt 2014, 20:41

Re: simple text output

Beitrag von Dexter22 »

Ok

I'm running from a command line in the win/bin folder

z88r -t -choly
then
z88r -c -choly

53.txt

Code: Alles auswählen

200000 0.3

z88elp.txt

Code: Alles auswählen

1
1 1 1.649E-005  2.278E-010 5.500E-003  2.278E-010  5.500E-003  4.556E-010  8.284E-008 0 0 1 0 0.833

z88i1.txt

Code: Alles auswählen

3 2 1 12 0
1 6 0 0 0
2 6 1 0 0
1 25
1 2
z88i2.txt

Code: Alles auswählen

2
1 1 2 0
1 2 2 0
1 3 2 0
1 4 2 0
1 5 2 0
1 6 2 0
2 3 1 -10
z88i5.txt z88int.txt - manual says that the entries in this don't matter for beam 25 structures? z88man.txt

Code: Alles auswählen

DYNAMIC START
---------------------------------------------------------------------------
Z88V14OS
---------------------------------------------------------------------------

---------------------------------------------------------------------------
GLOBAL
---------------------------------------------------------------------------

GLOBAL START
   IBFLAG          1
   IPFLAG          0
   IHFLAG          0
GLOBAL END

---------------------------------------------------------------------------
LINEAR SOLVER
---------------------------------------------------------------------------

SOLVER START
   MAXIT           10000
   EPS             1e-007
   RALPHA          0.0001
   ROMEGA          1.1
SOLVER END

---------------------------------------------------------------------------
STRESS
---------------------------------------------------------------------------

STRESS START
   KDFLAG        0
   ISFLAG        0
STRESS END

DYNAMIC END

z88mat.txt

Code: Alles auswählen

1
1 1 53.txt

z88r.dyn

Code: Alles auswählen

DYNAMIC START
---------------------------------
* Z88R.DYN build by Z88R Testmode
---------------------------------
CHOLESKY solver was selected
---------------------------------
* Language for Z88R
---------------------------------
ENGLISH
---------------------------------
* Entries for Z88R
---------------------------------
COMMON START
  MAXGS   82
  MAXIEZ  82
  MAXKOI  24
  MAXK    6
  MAXE    5
  MAXNFG  16
  MAXMAT  5
  MAXPEL  5
  MAXJNT  5
  MAXRBD  6
  MAXPR   5
  MAXGP   4
COMMON END
DYNAMIC END

z88r.log

Code: Alles auswählen

opening file Z88.DYN
Solver Cholesky, ICFLAG = 0
opening file Z88R.DYN
MAXGS   = 82
MAXKOI  = 24
MAXK    = 6
MAXE    = 5
MAXNFG  = 16
MAXMAT  = 5
MAXPEL  = 5
MAXJNT = 5
MAXESM = 3600
MAXPR = 5
MAXIEZ  = 5
file Z88R.DYN read ..seems to be o.k.
allocating dynamic memory:
memory id 1 allocated
memory id 10 allocated
memory id 11 allocated
memory id 16 allocated
memory id 30 allocated
memory id 31 allocated
memory id 32 allocated
memory id 40 allocated
memory id 41 allocated
memory id 42 allocated
memory id 43 allocated
memory id 60 allocated
memory id 61 allocated
memory id 62 allocated
memory id 63 allocated
memory id 64 allocated
memory id 65 allocated
memory id 70 allocated
memory id 71 allocated
memory id 72 allocated
memory id 73 allocated
memory id 80 allocated
memory id 81 allocated
memory id 82 allocated
memory id 83 allocated
memory id 84 allocated
memory id 85 allocated
memory id 86 allocated
memory id 87 allocated
memory id 88 allocated
memory id 89 allocated
memory id 90 allocated
memory id 91 allocated
memory id 92 allocated
memory id 93 allocated
memory id 100 allocated
memory id 101 allocated
memory id 102 allocated
memory id 103 allocated
memory id 104 allocated
memory id 105 allocated
memory id 110 allocated
memory id 111 allocated
memory id 112 allocated
memory id 113 allocated
memory id 114 allocated
memory id 115 allocated
memory id 120 allocated
memory id 121 allocated
memory id 122 allocated
memory id 123 allocated
memory id 124 allocated
memory id 125 allocated
memory id 130 allocated
memory id 131 allocated
memory id 132 allocated
memory id 133 allocated
dynamic memory totally allocated: 0 MB
leaving storage function DYN88R
start reading function RI188I
reading Z88I1.TXT
reading coordinates
reading element information
Z88I1.TXT read, leaving reading function RI188I
start Z88A
leaving Z88A
start writing function WRIA88R
writing Z88O0.TXT
vector GS needs 78 Elemente
vector KOI needs 20 Elemente
leaving writing function WRIA88R
start Z88BR

XK[1]= +0.00000E+000 , YK[1]= +0.00000E+000 , ZK[1]= +0.00000E+000 
XK[2]= +1.00000E+000 , YK[2]= +0.00000E+000 , ZK[2]= +0.00000E+000 
XCP  = +0.00000E+000 , YCP  = +1.00000E+000 , ZCP  = +0.00000E+000 
x-Achse   :A21X= +1.00000E+000 , A21Y= +0.00000E+000 , A21Z= +0.00000E+000 , RL21= +1.00000E+000 
Hilfsachse:A31X= +0.00000E+000 , A31Y= +1.00000E+000 , A31Z= +0.00000E+000 , RL31= +1.00000E+000 
y-Achse   :AYSX= +0.00000E+000 , AYSY= +0.00000E+000 , AYSZ= +1.00000E+000 , RLYS= +1.00000E+000 
z-Achse   :AZSX= +0.00000E+000 , AZSY= -1.00000E+000 , AZSZ= +0.00000E+000 , RLZS= +1.00000E+000 
VXX= +1.00000E+000 , VXY= +0.00000E+000 , VXZ= +0.00000E+000 
VYX= +0.00000E+000 , VYY= +0.00000E+000 , VYZ= +1.00000E+000 
VZX= +0.00000E+000 , VZY= -1.00000E+000 , VZZ= +0.00000E+000 
leaving Z88BR
start Z88CR
reading Z88I2.TXT
incorporating constraints pass 1
reading Z88I2.TXT
incorporating constraints pass 2
writing Z88O1.TXT
start SCAL88
leaving Z88CR
computing stresses
Z88DR done
computing nodal forces

XK[1]= +0.00000E+000 , YK[1]= +0.00000E+000 , ZK[1]= +0.00000E+000 
XK[2]= +1.00000E+000 , YK[2]= +0.00000E+000 , ZK[2]= +0.00000E+000 
XCP  = +0.00000E+000 , YCP  = +0.00000E+000 , ZCP  = +0.00000E+000 
x-Achse   :A21X= +1.00000E+000 , A21Y= +0.00000E+000 , A21Z= +0.00000E+000 , RL21= +1.00000E+000 
Hilfsachse:A31X= -1.#IND0E+000 , A31Y= -1.#IND0E+000 , A31Z= -1.#IND0E+000 , RL31= +0.00000E+000 
y-Achse   :AYSX= -1.#IND0E+000 , AYSY= -1.#IND0E+000 , AYSZ= -1.#IND0E+000 , RLYS= -1.#IND0E+000 
z-Achse   :AZSX= -1.#IND0E+000 , AZSY= -1.#IND0E+000 , AZSZ= -1.#IND0E+000 , RLZS= -1.#IND0E+000 
VXX= +1.00000E+000 , VXY= +0.00000E+000 , VXZ= +0.00000E+000 
VYX= -1.#IND0E+000 , VYY= -1.#IND0E+000 , VYZ= -1.#IND0E+000 
VZX= -1.#IND0E+000 , VZY= -1.#IND0E+000 , VZZ= -1.#IND0E+000 
Z88ER done
linear solver Z88R done
Solver successfull


I don't think I missed any files, but please let me know if I did.
mifris
Alumni
Beiträge: 153
Registriert: Mo 27. Sep 2010, 20:41

Re: simple text output

Beitrag von mifris »

Hello Dexter22,

your example has seven boundary conditions and you wrote only two in the z88i2.txt file! :-)
Therefore your example is statically undetermined.
And yes, there is no need of z88int.txt for beams, but the solver needs the information, that there is a file calls z88int.txt.

Best regards,
M. Frisch
Dexter22
Newbie
Beiträge: 8
Registriert: Fr 10. Okt 2014, 20:41

Re: simple text output

Beitrag von Dexter22 »

Ah ok, I made a silly mistake, thanks.

I've re-run the simulation but now I don't get any nodal forces, only displacements.

Do you have any more suggestions?

Best Regards,
Dexter

Z88O4.TXT

Code: Alles auswählen

output file Z88O4.TXT : nodal forces, computed by Z88R V15.0B
                        ************

starting with nodal forces computed for each element
----------------------------------------------------

element # =     1     type = beam no.25 in space
node         F(1)           F(2)           F(3)           F(4)           F(5)           F(6)
    1   -1.#IND0E+000  -1.#IND0E+000  -1.#IND0E+000  -1.#IND0E+000  -1.#IND0E+000  -1.#IND0E+000
    2   -1.#IND0E+000  -1.#IND0E+000  -1.#IND0E+000  -1.#IND0E+000  -1.#IND0E+000  -1.#IND0E+000

the nodal sums for each node
----------------------------
  node       F(1)           F(2)           F(3)           F(4)           F(5)           F(6)

    1   -1.#IND0E+000  -1.#IND0E+000  -1.#IND0E+000  -1.#IND0E+000  -1.#IND0E+000  -1.#IND0E+000
    2   -1.#IND0E+000  -1.#IND0E+000  -1.#IND0E+000  -1.#IND0E+000  -1.#IND0E+000  -1.#IND0E+000
z88r.log

Code: Alles auswählen

opening file Z88.DYN
Solver Cholesky, ICFLAG = 0
opening file Z88R.DYN
MAXGS   = 82
MAXKOI  = 24
MAXK    = 6
MAXE    = 5
MAXNFG  = 16
MAXMAT  = 5
MAXPEL  = 5
MAXJNT = 5
MAXESM = 3600
MAXPR = 5
MAXIEZ  = 5
file Z88R.DYN read ..seems to be o.k.
allocating dynamic memory:
memory id 1 allocated
memory id 10 allocated
memory id 11 allocated
memory id 16 allocated
memory id 30 allocated
memory id 31 allocated
memory id 32 allocated
memory id 40 allocated
memory id 41 allocated
memory id 42 allocated
memory id 43 allocated
memory id 60 allocated
memory id 61 allocated
memory id 62 allocated
memory id 63 allocated
memory id 64 allocated
memory id 65 allocated
memory id 70 allocated
memory id 71 allocated
memory id 72 allocated
memory id 73 allocated
memory id 80 allocated
memory id 81 allocated
memory id 82 allocated
memory id 83 allocated
memory id 84 allocated
memory id 85 allocated
memory id 86 allocated
memory id 87 allocated
memory id 88 allocated
memory id 89 allocated
memory id 90 allocated
memory id 91 allocated
memory id 92 allocated
memory id 93 allocated
memory id 100 allocated
memory id 101 allocated
memory id 102 allocated
memory id 103 allocated
memory id 104 allocated
memory id 105 allocated
memory id 110 allocated
memory id 111 allocated
memory id 112 allocated
memory id 113 allocated
memory id 114 allocated
memory id 115 allocated
memory id 120 allocated
memory id 121 allocated
memory id 122 allocated
memory id 123 allocated
memory id 124 allocated
memory id 125 allocated
memory id 130 allocated
memory id 131 allocated
memory id 132 allocated
memory id 133 allocated
dynamic memory totally allocated: 0 MB
leaving storage function DYN88R
start reading function RI188I
reading Z88I1.TXT
reading coordinates
reading element information
Z88I1.TXT read, leaving reading function RI188I
start Z88A
leaving Z88A
start writing function WRIA88R
writing Z88O0.TXT
vector GS needs 78 Elemente
vector KOI needs 20 Elemente
leaving writing function WRIA88R
start Z88BR

XK[1]= +0.00000E+000 , YK[1]= +0.00000E+000 , ZK[1]= +0.00000E+000 
XK[2]= +1.00000E+000 , YK[2]= +0.00000E+000 , ZK[2]= +0.00000E+000 
XCP  = +0.00000E+000 , YCP  = +1.00000E+000 , ZCP  = +0.00000E+000 
x-Achse   :A21X= +1.00000E+000 , A21Y= +0.00000E+000 , A21Z= +0.00000E+000 , RL21= +1.00000E+000 
Hilfsachse:A31X= +0.00000E+000 , A31Y= +1.00000E+000 , A31Z= +0.00000E+000 , RL31= +1.00000E+000 
y-Achse   :AYSX= +0.00000E+000 , AYSY= +0.00000E+000 , AYSZ= +1.00000E+000 , RLYS= +1.00000E+000 
z-Achse   :AZSX= +0.00000E+000 , AZSY= -1.00000E+000 , AZSZ= +0.00000E+000 , RLZS= +1.00000E+000 
VXX= +1.00000E+000 , VXY= +0.00000E+000 , VXZ= +0.00000E+000 
VYX= +0.00000E+000 , VYY= +0.00000E+000 , VYZ= +1.00000E+000 
VZX= +0.00000E+000 , VZY= -1.00000E+000 , VZZ= +0.00000E+000 
leaving Z88BR
start Z88CR
reading Z88I2.TXT
incorporating constraints pass 1
reading Z88I2.TXT
incorporating constraints pass 2
writing Z88O1.TXT
start SCAL88
leaving Z88CR
computing stresses
Z88DR done
computing nodal forces

XK[1]= +0.00000E+000 , YK[1]= +0.00000E+000 , ZK[1]= +0.00000E+000 
XK[2]= +1.00000E+000 , YK[2]= +0.00000E+000 , ZK[2]= +0.00000E+000 
XCP  = +0.00000E+000 , YCP  = +0.00000E+000 , ZCP  = +0.00000E+000 
x-Achse   :A21X= +1.00000E+000 , A21Y= +0.00000E+000 , A21Z= +0.00000E+000 , RL21= +1.00000E+000 
Hilfsachse:A31X= -1.#IND0E+000 , A31Y= -1.#IND0E+000 , A31Z= -1.#IND0E+000 , RL31= +0.00000E+000 
y-Achse   :AYSX= -1.#IND0E+000 , AYSY= -1.#IND0E+000 , AYSZ= -1.#IND0E+000 , RLYS= -1.#IND0E+000 
z-Achse   :AZSX= -1.#IND0E+000 , AZSY= -1.#IND0E+000 , AZSZ= -1.#IND0E+000 , RLZS= -1.#IND0E+000 
VXX= +1.00000E+000 , VXY= +0.00000E+000 , VXZ= +0.00000E+000 
VYX= -1.#IND0E+000 , VYY= -1.#IND0E+000 , VYZ= -1.#IND0E+000 
VZX= -1.#IND0E+000 , VZY= -1.#IND0E+000 , VZZ= -1.#IND0E+000 
Z88ER done
linear solver Z88R done
Solver successfull

mifris
Alumni
Beiträge: 153
Registriert: Mo 27. Sep 2010, 20:41

Re: simple text output

Beitrag von mifris »

Hello again,

hmm, that`s weird. Can you please send us your project folder?! (z88aurora-at-uni-bayreuth.de) That makes it easier for us to have a look at your example. :-)

Best regards,

M. Frisch
Dexter22
Newbie
Beiträge: 8
Registriert: Fr 10. Okt 2014, 20:41

Re: simple text output

Beitrag von Dexter22 »

Ok, I've sent the files.

Thanks for the help.

-Dexter
mifris
Alumni
Beiträge: 153
Registriert: Mo 27. Sep 2010, 20:41

Re: simple text output

Beitrag von mifris »

Thank you! We´ll check it as soon as possible! :-)

Best regards,
M. Frisch
Dexter22
Newbie
Beiträge: 8
Registriert: Fr 10. Okt 2014, 20:41

Re: simple text output

Beitrag von Dexter22 »

Thanks for your help.

like you said via email, it was the orientation control point that was wrong. I'm not sure why though:

In my case, I had a beam 25 element from (0,0,0) to (100,0,0)[on the x axis]

The point 3 (control point)

1 - Does not work at (0,1,0) [not on the x axis]
2 - Does work at (0,0,1)[not on the x axis](as expected from the theory manual.)
3 - Does not work at (1,0,0) [on the x axis] (as expected from the theory manual.)

So now my question is why does the control point not work when it is on the y axis? In the theory manual, the only requirement is that it doesn't lie on the axis of the beam element (in this case the x axis)
Antworten