However, when I run
Code: Alles auswählen
z88r -c -siccg
Code: Alles auswählen
*** SICCG solver with sparse matrix storing ***
*** computing linear structure ***
reading solver data file Z88MAN.TXT
reading structure data file Z88I1.TXT
reading coordinates
reading element informations
reading integration order data file Z88INT.TXT
reading element parameters data file Z88ELP.TXT
reading material data file Z88MAT.TXT
>>> start Z88AI <<<
program allocated 191 MB memory
still 942 steps, pointer IEZ = 1
The z88r.log ends with
Code: Alles auswählen
start reading function RI188I
reading Z88I1.TXT
reading coordinates
reading element information
Z88I1.TXT read, leaving reading function RI188R
start Z88AI
Code: Alles auswählen
*** Cholesky solver with Jennings storing ***
*** computing linear structure ***
reading solver data file Z88MAN.TXT
reading structure data file Z88I1.TXT
reading coordinates
reading element informations
reading integration order data file Z88INT.TXT
reading element parameters data file Z88ELP.TXT
reading material data file Z88MAT.TXT
>>> start Z88A <<<
no. 182 type 7
writing Z88O0.TXT
vector GS needs 1436523 elements
vector KOI needs 1468 elements
>>> start Z88BR <<<
program allocated 92 MB memory
no. 182 type 7
>>> start Z88CR <<<
reading constraint file Z88I2.TXT
incorporating constraints pass 1
constraint no. 979 type 1incorporating constraints pass 2
writing Z88O1.TXT
start SCAL88
### diagonal element in array zero or negative..stop ###
### often caused by missing or wrong constraints ###
### recover: check constraints (underdefined ?) ###
All the elements are type7, 8-point quads, simulating 1mm-thick steel sheet.
I can view the structure in z88o, which displays the wireframe and boundary conditions correctly.
z88mat.txt:
1
1 182 51.txt
z88int.txt:
1
1 182 3 3
z88elp.txt:
1
1 182 1.000000 0 0 0 0 0 0
z88i1.txt:
2 943 182 1850 0
1 2 3.738256 -2.000000
2 2 3.083860 -3.084800
3 0 2.000000 -3.751697
...
1 7
114 142 162 128 216 217 218 219
2 7
94 149 86 66 220 221 222 223
3 7
131 117 129 169 224 225 226 227
...
z88i2.txt:
979
1 1 1 0.032369
2 1 1 0.061532
3 0 2 0.000000
3 1 2 0.000000
3 1 1 0.038771
4 0 2 0.000000
4 1 2 0.000000
4 1 1 0.022045
...
The object being simulated has lots of zero-displacement points at the bottom, and forces on all the points simulating gravity, so it shouldn't be underconstrained.
Could anyone offer suggestions as to what I might be doing wrong? I can PM the files if people want to try running it themselves, but the forum won't let me upload .txt attachments.
Thanks!
Tom