search

Dhondtguido / CalculiX

This repository contains the source files of CalculiX, a three-dimensional Finite Element Program (www.calculix.de).
GNU General Public License v2.0
72 stars 16 forks source link

BUG: CONTACT ELEMENTS breaks displacement on composites. #56

Open victorkemp opened 11 months ago

victorkemp commented 11 months ago

Here, the composite element should deflect but it only does when CONTACT ELEMENTS is removed from *NODE FILE. cec

*NODE
1,0.5049519538879,0.1484981536865,0.8577842116356
2,0.3799519538879,0.1484981536865,0.4721890091896
3,0.5049519538879,0.1484981536865,0.08659380674362
4,0.3799519538879,0.1484981536865,0.8577842116356
5,0.3799519538879,0.1484981536865,0.08659380674362
6,0.6299519538879,0.1484981536865,0.4721890091896
7,0.6299519538879,0.1484981536865,0.08659380674362
8,0.6299519538879,0.1484981536865,0.8577842116356
9,0,0,0
10,1,0,0
11,1,0.1,0
12,0,0.1,0
13,0,0,1
14,1,0,1
15,1,0.1,1
16,0,0.1,1
17,1,0.1,0.5
18,1,0.05,0
19,1,0,0.5
20,1,0.05,1
21,0.5,0.1,1
22,0.5,0.1,0
23,0.5,0,0
24,0.5,0,1
25,0,0.1,0.5
26,0,0.05,0
27,0,0,0.5
28,0,0.05,1
*ELEMENT,TYPE=C3D20
1,15,11,10,14,16,12,9,13,17,18,19,20,25,26,27,
28,21,22,23,24
*ELEMENT,TYPE=S8R
2,7,8,4,5,6,1,2,3
*ELSET,ELSET=Default
1
*ELSET,ELSET=Component
2
*SURFACE,NAME=pressure_faces
2,SNEG
*SURFACE,NAME=contact_faces(2)
2,SPOS
*SURFACE,NAME=contact_faces
1,S3
*ORIENTATION,NAME=O_1_L_1
0,0,1,-1,0,0
*ORIENTATION,NAME=O_1_L_2
0,0,1,-1,0,0
*MATERIAL,NAME=Material(2)
*ELASTIC,TYPE=ISOTROPIC
10000000000,0
*MATERIAL,NAME=Material_1_L_1
*ELASTIC,TYPE=ENGINEERINGCONSTANTS
100000000000,100000000000,100000000000,0,0,0,50000000000,50000000000
50000000000
*MATERIAL,NAME=Material_1_L_2
*ELASTIC,TYPE=ENGINEERINGCONSTANTS
100000000000,100000000000,100000000000,0,0,0,50000000000,50000000000
50000000000
*SOLID SECTION,ELSET=Default,MATERIAL=Material(2)
*SHELL SECTION,ELSET=Component,COMPOSITE
0.03,,Material_1_L_1,O_1_L_1
0.03,,Material_1_L_2,O_1_L_2
*BOUNDARY
1,1,6,0
4,1,6,0
8,1,6,0
13,1,3,0
14,1,3,0
15,1,3,0
16,1,3,0
20,1,3,0
21,1,3,0
24,1,3,0
28,1,3,0
*CONTACT PAIR,INTERACTION=SI_2,TYPE=SURFACE TO SURFACE
contact_faces(2),contact_faces
*SURFACE INTERACTION,NAME=SI_2
*SURFACE BEHAVIOR,PRESSURE-OVERCLOSURE=LINEAR
1000000000000
*STEP,NLGEOM=YES
*STATIC
1,1,0,0
*DSLOAD
pressure_faces,P,10000000
*NODE FILE,CONTACT ELEMENTS
U
*END STEP