TYPE ORTHO not known

mgbri · September 28, 2021, 7:34am

Good morning,

I’m trying to simulate a multilayered plate with different ply orientations. But I can’t manage to obtain results. I use S8R elements. If I inspect the error file I get:

400 elements reallocated
WARNING: TYPE:ORTHO not known
WARNING: option:COMPOSITE not supported, keyword ignored
WARNING: *ORIENTATION not known. Data ignored.
WARNING: *ORIENTATION not known. Data ignored.
WARNING: *ELFILE not known. Data ignored.
WARNING: *ELPRINT not known. Data ignored.
WARNING: *ELPRINT not known. Data ignored.
WARNING: *ELPRINT not known. Data ignored.

My input file looks like this:

*MATERIAL, NAME=cfrp
    *ELASTIC, TYPE=ORTHO
        500000.,157200.,400000.,157200.,157200.,300000.,126200.,126200.,
        126200.,294.
*SHELL SECTION,COMPOSITE, ELSET=all
    0.5,,cfrp,orient
    0.5,,cfrp,ninety
    0.5,,cfrp,orient
*ORIENTATION, NAME=ninety, SYSTEM=RECTANGULAR
    0.,0.,0.,1.,0.,0.
*ORIENTATION, NAME=orient, SYSTEM=RECTANGULAR
    0.,0.,0.,0.,1.,0.
*STEP
    *STATIC, SOLVER=SPOOLES
    *EL FILE
        S
    *EL PRINT, ELSET=all
        S
    *EL PRINT, ELSET=all
        E
    *EL PRINT, ELSET=all
        ENER
    *NODE PRINT, NSET=all,OUTPUT=2D
        U
    *BOUNDARY, OP=NEW
        bc, 1,6
    *CLOAD, OP=NEW
        load, 2, 100
    *END STEP

I really can’t understand what’s wrong. This is also happening when I use TYPE=ENGINEERING CONSTANTS under the *ELASTIC card. Any idea/tips?

Best,
Massimo

Calc_em · September 28, 2021, 8:35am

Which version of CalculiX do you use ? Could you upload the whole file, with mesh and set definitions ? Maybe run tests on a single element model for now.

mgbri · September 28, 2021, 9:00am

I run the binaries contained in CalculiX Advanced Environment v0.8.0 for windows
([Release CalculiX Advanced Environment v0.8.0 · calculix/cae · GitHub]).

Here is the .imp:

*node, nset=Nall
1, 5.000000e+001, 5.000000e+001, 0.000000e+000
2, 5.000000e+001, -5.000000e+001, 0.000000e+000
3, -5.000000e+001, 5.000000e+001, 0.000000e+000
4, -5.000000e+001, -5.000000e+001, 0.000000e+000
5, -5.000000e+001, -1.666667e+001, 0.000000e+000
6, -5.000000e+001, 1.666667e+001, 0.000000e+000
7, -1.666667e+001, -5.000000e+001, 0.000000e+000
8, 1.666667e+001, -5.000000e+001, 0.000000e+000
9, 5.000000e+001, -1.666667e+001, 0.000000e+000
10, 5.000000e+001, 1.666667e+001, 0.000000e+000
11, -1.666667e+001, 5.000000e+001, 0.000000e+000
12, 1.666667e+001, 5.000000e+001, 0.000000e+000
13, -1.666667e+001, -1.666667e+001, 0.000000e+000
14, -1.666667e+001, 1.666667e+001, 0.000000e+000
15, 1.666667e+001, -1.666667e+001, 0.000000e+000
16, 1.666667e+001, 1.666667e+001, 0.000000e+000
17, -5.000000e+001, -3.333333e+001, 0.000000e+000
18, -5.000000e+001, -9.983125e-013, 0.000000e+000
19, -5.000000e+001, 3.333333e+001, 0.000000e+000
20, -3.333333e+001, -5.000000e+001, 0.000000e+000
21, -9.983125e-013, -5.000000e+001, 0.000000e+000
22, 3.333333e+001, -5.000000e+001, 0.000000e+000
23, 5.000000e+001, -3.333333e+001, 0.000000e+000
24, 5.000000e+001, -9.983125e-013, 0.000000e+000
25, 5.000000e+001, 3.333333e+001, 0.000000e+000
26, -3.333333e+001, 5.000000e+001, 0.000000e+000
27, -9.983125e-013, 5.000000e+001, 0.000000e+000
28, 3.333333e+001, 5.000000e+001, 0.000000e+000
29, -1.666667e+001, -3.333333e+001, 0.000000e+000
30, -3.333333e+001, -1.666667e+001, 0.000000e+000
31, -1.666667e+001, -4.796163e-013, 0.000000e+000
32, -3.333333e+001, 1.666667e+001, 0.000000e+000
33, -1.666667e+001, 3.333333e+001, 0.000000e+000
34, -4.867218e-013, -1.666667e+001, 0.000000e+000
35, 1.666667e+001, -3.333333e+001, 0.000000e+000
36, 3.333333e+001, -1.666667e+001, 0.000000e+000
37, 1.666667e+001, -4.760636e-013, 0.000000e+000
38, 3.333333e+001, 1.666667e+001, 0.000000e+000
39, 1.666667e+001, 3.333333e+001, 0.000000e+000
40, -4.725109e-013, 1.666667e+001, 0.000000e+000
*nset, nset=bc
7, 8, 20, 21, 22, 4, 2, 
*nset, nset=load
1, 3, 11, 12, 26, 27, 28, 
*nset, nset=all
5, 6, 17, 18, 19, 9, 10, 
23, 24, 25, 1, 7, 8, 20, 21, 
22, 13, 14, 15, 16, 29, 30, 31, 
32, 33, 34, 35, 36, 37, 38, 39, 
40, 3, 4, 2, 11, 12, 26, 27, 
28, 
*element, elset=S8R,type=S8R
13, 5, 13, 7, 4, 30, 29, 20, 17, 
14, 6, 14, 13, 5, 32, 31, 30, 18, 
15, 3, 11, 14, 6, 26, 33, 32, 19, 
16, 7, 13, 15, 8, 29, 34, 35, 21, 
17, 8, 15, 9, 2, 35, 36, 23, 22, 
18, 9, 15, 16, 10, 36, 37, 38, 24, 
19, 10, 16, 12, 1, 38, 39, 28, 25, 
20, 12, 16, 14, 11, 39, 40, 33, 27, 
21, 14, 16, 15, 13, 40, 37, 34, 31, 

*elset, elset=all
13, 14, 15, 16, 17, 18, 19, 
20, 21, 

    
*MATERIAL, NAME=cfrp
    *ELASTIC, TYPE=ORTHO
        500000.,157200.,400000.,157200.,157200.,300000.,126200.,126200.,
        126200.,294.
*SHELL SECTION,COMPOSITE, ELSET=all
    0.5,,cfrp,orient
    0.5,,cfrp,ninety
    0.5,,cfrp,orient
*ORIENTATION, NAME=ninety, SYSTEM=RECTANGULAR
    0.,0.,0.,1.,0.,0.
*ORIENTATION, NAME=orient, SYSTEM=RECTANGULAR
    0.,0.,0.,0.,1.,0.
*STEP
    *STATIC, SOLVER=SPOOLES
    *EL FILE
        S
    *EL PRINT, ELSET=all
        S
    *EL PRINT, ELSET=all
        E
    *EL PRINT, ELSET=all
        ENER
    *NODE PRINT, NSET=all,OUTPUT=2D
        U
    *BOUNDARY, OP=NEW
        bc, 1,6
    *CLOAD, OP=NEW
        load, 2, 100
    *END STEP

Nodesets “bc”, “load”, and “all” are respectively one edge of the plate (fixed), the other side of the plate (loaded to load the plate in tension membrane), and all the nodes of the mesh.

The mesh definition has been generated through a conversion of a .unv export done in Salome platform to .inp

mgbri · September 28, 2021, 9:20am

Here in addition the log of the solution:

DEBUG: Running from binaries.
INFO: Application's home directory is:
C:\calculix\cae_20200725_windows
INFO: Work directory is:
C:\Users\brivio\Desktop\plate_calculix
DEBUG: Stopping threads:
thread_1_read_cgx_stdout_1632814078

INFO: Keywords object model generated.
INFO: Loading model
C:\Users\brivio\Desktop\plate_calculix\plate.inp

DEBUG: *node, nset=Nall
DEBUG: keyword_chain: *node
DEBUG: path found: Mesh, *NODE
INFO: *NODE Nall created.

DEBUG: *nset, nset=bc
DEBUG: keyword_chain: *node, *nset
DEBUG: path found: Mesh, *NSET
INFO: *NSET bc created.

DEBUG: *nset, nset=load
DEBUG: keyword_chain: *nset
DEBUG: path found: Mesh, *NSET
INFO: *NSET load created.

DEBUG: *nset, nset=all
DEBUG: keyword_chain: *nset
DEBUG: path found: Mesh, *NSET
INFO: *NSET all created.

DEBUG: *element, elset=S8R,type=S8R
DEBUG: keyword_chain: *nset, *element
DEBUG: path found: Mesh, *ELEMENT
INFO: *ELEMENT S8R created.

DEBUG: *elset, elset=all
DEBUG: keyword_chain: *element, *elset
DEBUG: path found: Mesh, *ELSET
INFO: *ELSET all created.

DEBUG: *MATERIAL, NAME=cfrp
DEBUG: keyword_chain: *elset, *MATERIAL
DEBUG: path found: Properties, *MATERIAL
INFO: *MATERIAL cfrp created.

DEBUG: *ELASTIC, TYPE=ORTHO
DEBUG: keyword_chain: *MATERIAL, *ELASTIC
DEBUG: path found: Properties, *MATERIAL, *ELASTIC
INFO: *ELASTIC ELASTIC-1 created.

DEBUG: *SHELL SECTION,COMPOSITE, ELSET=all
DEBUG: keyword_chain: *MATERIAL, *ELASTIC, *SHELL SECTION
DEBUG: path found: Properties, Section, *SHELL SECTION
INFO: *SHELL SECTION all created.

DEBUG: *ORIENTATION, NAME=ninety, SYSTEM=RECTANGULAR
DEBUG: keyword_chain: *SHELL SECTION, *ORIENTATION
DEBUG: path found: *ORIENTATION
INFO: *ORIENTATION ninety created.

DEBUG: *ORIENTATION, NAME=orient, SYSTEM=RECTANGULAR
DEBUG: keyword_chain: *ORIENTATION
DEBUG: path found: *ORIENTATION
INFO: *ORIENTATION orient created.

DEBUG: *STEP
DEBUG: keyword_chain: *ORIENTATION, *STEP
DEBUG: path found: *STEP
INFO: *STEP STEP-1 created.

DEBUG: *STATIC, SOLVER=SPOOLES
DEBUG: keyword_chain: *STEP, *STATIC
DEBUG: path found: *STEP, Analysis type, *STATIC
INFO: *STATIC STATIC-1 created.

DEBUG: *EL FILE
DEBUG: keyword_chain: *STEP, *STATIC, *EL FILE
DEBUG: path found: *STEP, Field Output, *EL FILE
INFO: *EL FILE EL FILE-1 created.

DEBUG: *EL PRINT, ELSET=all
DEBUG: keyword_chain: *STEP, *STATIC, *EL FILE, *EL PRINT
DEBUG: path found: *STEP, Field Output, *EL PRINT
INFO: *EL PRINT all created.

DEBUG: *EL PRINT, ELSET=all
DEBUG: keyword_chain: *STEP, *STATIC, *EL FILE, *EL PRINT
DEBUG: path found: *STEP, Field Output, *EL PRINT
INFO: *EL PRINT all created.

DEBUG: *EL PRINT, ELSET=all
DEBUG: keyword_chain: *STEP, *STATIC, *EL FILE, *EL PRINT
DEBUG: path found: *STEP, Field Output, *EL PRINT
INFO: *EL PRINT all created.

DEBUG: *NODE PRINT, NSET=all,OUTPUT=3D
DEBUG: keyword_chain: *STEP, *STATIC, *EL FILE, *EL PRINT, *NODE PRINT
DEBUG: path found: *STEP, Field Output, *NODE PRINT
INFO: *NODE PRINT all created.

DEBUG: *BOUNDARY, OP=NEW
DEBUG: keyword_chain: *STEP, *STATIC, *EL FILE, *EL PRINT, *NODE PRINT, *BOUNDARY
DEBUG: path found: *STEP, Load & BC, *BOUNDARY
INFO: *BOUNDARY BOUNDARY-1 created.

DEBUG: *CLOAD, OP=NEW
DEBUG: keyword_chain: *STEP, *STATIC, *EL FILE, *EL PRINT, *NODE PRINT, *BOUNDARY, *CLOAD
DEBUG: path found: *STEP, Load & BC, *CLOAD
INFO: *CLOAD CLOAD-1 created.

DEBUG: *END STEP
DEBUG: keyword_chain: *STEP, *STATIC, *EL FILE, *EL PRINT, *NODE PRINT, *BOUNDARY, *CLOAD, *END STEP
DEBUG: path found: *STEP, *END STEP
INFO: *END STEP END STEP-1 created.
ERROR: Can't parse elements
Traceback (most recent call last):
  File "model\parsers\mesh.py", line 64, in __init__
  File "model\parsers\mesh.py", line 225, in parse_elements
ValueError: invalid literal for int() with base 10: '*elset'

INFO: Mesh parser:
40 nodes
4 nsets
1 elsets
0 surfaces
DEBUG: Killed CGX, PID=33040.
DEBUG: CGX PID=194592
DEBUG: CalculiX GraphiX WID=0x004210e4

DEBUG: Logging threads:
thread_1_read_cgx_stdout_1632820676

 The win HOME was detected:C:\Users\brivio
parameters:3 arguments:2
GL_MAX_EVAL_ORDER:8, mesh threads:1
 Try to read ccx input

******************************************************************************
*  read ccx input                                   4.4.2002 wi              *
*                                                                            *
*                                                                            *
* -only free-format is supported. fixed format might cause a prg-crash.      *
*                                                                            *
*                                                                            *
******************************************************************************
 ignore statement: *NODEPRINT

 40 nodes reallocated

 9 elements reallocated
 WARNING: TYPE:ORTHO not known
WARNING: option:COMPOSITE not supported, keyword ignored
WARNING: *ORIENTATION not known. Data ignored.
WARNING: *ORIENTATION not known. Data ignored.
WARNING: *ELFILE not known. Data ignored.
WARNING: *ELPRINT not known. Data ignored.
WARNING: *ELPRINT not known. Data ignored.
WARNING: *ELPRINT not known. Data ignored.
 create dataset:0 name= +bou
 create dataset:1 name= +clo
 create dataset:2 name= +clo
WARNING: *ENDSTEP not known. Data ignored.
 Elements: 9 Nodes:40 Datasets:3  MinElemNr: 13 MaxElemNr: 21 MinNodNr:1 MaxNodNr:40  read in 0.059000 sec
 found elements of type 1: 0
 found elements of type 2: 0
 found elements of type 3: 0
 found elements of type 4: 0
 found elements of type 5: 0
 found elements of type 6: 0
 found elements of type 7: 0
 found elements of type 8: 0
 found elements of type 9: 0
 found elements of type 10: 9
 found elements of type 11: 0
 found elements of type 12: 0
gtol calculated:1.250000e-008

INFO: C:\calculix\cae_20200725_windows\bin\ccx.exe -i C:\Users\brivio\Desktop\plate_calculix\plate

DEBUG: Logging threads:
thread_1_submit_1632820692
thread_2_read_stdout_1632820692

INFO: END

Calc_em · September 28, 2021, 2:33pm

I think that I’ve found the reason for this problem. It’s caused by orientation definition. Keep in mind that in the case of shell elements local coordinate system is created by the projection of axes defined as *Orientation and there are specific rules of this projection (described in the User’s Manual - chapter about shell elements). In certain conditions, CalculiX may not be able to do this projection. Try defining orientations in such a way that the X axis is aligned with shell normal.

mgbri · September 29, 2021, 1:15pm

Hi Calc_em, indeed this solved the problem. Thank you so much. Orientation was not defined as I thought and I was misled by the error messages.
What I actually wanted to do with orientation was the following:

*ORIENTATION, NAME=ninety, SYSTEM=RECTANGULAR
    1.,0.,0.,0.,1.,0.
*ORIENTATION, NAME=zero, SYSTEM=RECTANGULAR
    0.,1.,0.,1.,0.,0.

salah · March 8, 2023, 2:55pm

Can you help me in errors ?

salah · March 8, 2023, 2:56pm

BOUNDARY
Nfixx,1,1,0
15,2,0
1,2,0
*NSET,NSET=NALL,GENERATE
1,725
*ELSET,ELSET=EALL,GENERATE
1,658
*MATERIAL, NAME=E1
*ELASTIC,TYPE=ORTHO
17000.,15000,0.30,5000,4000,2000
*SOLID SECTION,ELSET=Cancellous,MATERIAL=E1
*MATERIAL, NAME=E2
*ELASTIC,TYPE=ORTHO
17000.,15000,0.30,5000,4000,2000
SOLID SECTION,ELSET=Cortical,MATERIAL=E2

Calc_em · March 8, 2023, 6:04pm

Add the *ORIENTATION keyword, as explained above.

salah · March 9, 2023, 6:55pm

Excuse me, i’m begining with this code. Where can i add the orientation *ORIENTATION,NAME=OR
1.,1.,0.,0.,0.,0.
for the two bone but the program gives error in abaqus and calculix

salah · March 9, 2023, 6:57pm

Can i send you the entire file?

Calc_em · March 9, 2023, 7:17pm

Sure, will be much easier to help. If the file is large, you can use some hosting website like WeTransfer, Dropbox or Google Drive and paste the link.

salah · March 9, 2023, 7:53pm

I cannot send you the whole of the file. But the principal part is :
*NSET,NSET=NALL,GENERATE
*1,725
*LSET,ELSET=EALL,GENERATE
*1,658
*MATERIAL, NAME=E1
*ELASTIC,TYPE=ORTHOTROPIC
17000.,15000,0.30,5000,4000,2000
*ORIENTATION,NAME=OR
1.,1.,0.,0.,0.,0.
*SOLID SECTION,ELSET=Cancellous,MATERIAL=E1,ORIENTATION=OR
*MATERIAL, NAME=E2
*ELASTIC,TYPE=ORTHOTROPIC
17000.,15000,0.30,5000,4000,2000
*ORIENTATION,NAME=OR
1.,1.,0.,0.,0.,0.

SOLID SECTION,ELSET=Cortical,MATERIAL=E2,ORIENTATION=OR
*STEP
*STATIC
*CLOAD
720,1,-385.
*ELPRINT,ELSET=Eall
S
*NODE FILE
U
*EL FILE
S,E
*ENDSTEP

salah · March 9, 2023, 7:56pm

we have a 2D problem with two different material (cortical bone and cancellous bone). For more detail, send me you mail to send the entire file.

salah · March 12, 2023, 11:55am

Hi, colleague Calc_em i wail your response, thank you

Calc_em · March 14, 2023, 3:23pm

What is the error message that you get ? It will be hard to help without files for testing. Maybe you could prepare some simplified ones that show the same problem.

salah · April 2, 2023, 5:16am

Thank you very much,
Please send me you mail adress to send you the data file.
Thanks

Topic		Replies	Views
Error running element level calculix input file	2	133	March 12, 2024
Orthotropic Shell/Plate Analysis	22	1821	May 25, 2022
Modelling of composites with solid sections	20	1233	May 30, 2023
Layered shell (composite) with spring element	38	1668	February 10, 2023
Buckling of composite panels	64	1580	June 25, 2024

TYPE ORTHO not known

Related topics