[Fixed] brick20 (hexa20) mesh gives other results than tet10

[bernd]

magnitude: 87.55 mm
max z: -87.17 mm


OS: Debian GNU/Linux 10 (buster) (KDE//usr/share/xsessions/plasma)
Word size of OS: 64-bit
Word size of FreeCAD: 64-bit
Version: 0.19.19766 (Git)
Build type: Unknown
Branch: master
Hash: 2d6aa10b69f9d2c9a37ba08f33cbc5b528c2b020
Python version: 3.7.3
Qt version: 5.11.3
Coin version: 4.0.0a
OCC version: 7.3.0
Locale: German/Switzerland (de_CH)

canti_brick20.FCStd

(49.74 KiB) Downloaded 25 times

Screenshot_20200301_221305.png (174.14 KiB) Viewed 594 times

[bernd]

git commit eeb97b36 adds this example to the code and git commit 3f399736 makes a unit test ...

to run the example from Python directly use:

Code: Select all

from femexamples import ccx_cantilever_std as canti
canti.setup_cantileverhexa20faceload()

[ebrahim raeyat]

magnitude: 87.55 mm
max z: -87.17 mm

in opensees i gives this result:

-87.1642 mm!!

[bernd]

exact deformation in ccx on node 5

Code: Select all

8.13413E+00 1.73873E-03 -8.71642E+01

means -87.1642 mm too. This makes sense if the same mesh and element type is used the same (exact the same) deflection. For tetra10 element and z88 and calculix it is the same, exact the same deflection.

[ebrahim raeyat]

exact deformation in ccx on node 5

Code: Select all

8.13413E+00 1.73873E-03 -8.71642E+01

means -87.1642 mm too. This makes sense if the same mesh and element type is used the same (exact the same) deflection. For tetra10 element and z88 and calculix it is the same, exact the same deflection.

in opensees node 5:

Code: Select all

8.13413 0.00173873 -87.1642

Example8.1.py

(1.58 KiB) Downloaded 28 times

nodes.py

(7.69 KiB) Downloaded 27 times

element.py

(3.78 KiB) Downloaded 24 times

load.py

(731 Bytes) Downloaded 26 times

i could not load .tcl file, thus i renamed those to .py file! i will soon upload python file.

[bernd]

we should move all opensees related posts in the opensees topic.


EDIT: ahh you did already. OK than any disscussion in this regard will continue there ... see https://forum.freecadweb.org/viewtopic. ... 70#p373327