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
75 stars 19 forks source link

Truss and shell connection overconstrained #11

Open victorkemp opened 1 year ago

victorkemp commented 1 year ago

The shell should bend like a cantilever with unconstrained rotation at the end connected to the truss elements. If that happened, nodes 49 and 51 would have opposite values for X displacement.

With the MKL solver, they have the same value, showing that rotation is incorrectly constrained. With the Spooles solver, it doesn't generate a solution.

solution

*NODE
1,1.5,0,-1
2,2.4,0,-1
3,1.8,0,-1
4,0,-0.3,0
5,2.1,-0.3,-1
6,2.4,-0.3,-1
7,1.5,-0.3,-1
8,0,0,0
9,2.1,0,-1
10,1.8,-0.3,-1
*ELEMENT,TYPE=T3D2,ELSET=Truss
1,1,8
4,7,4
*ELEMENT,TYPE=S4,ELSET=Shell
2,9,3,10,5
3,3,1,7,10
5,2,9,5,6
*MATERIAL,NAME=Material
*ELASTIC,TYPE=ISOTROPIC
200000000000,0
*SHELL SECTION,ELSET=Shell,MATERIAL=Material
0.01
*SOLID SECTION,ELSET=Truss,MATERIAL=Material
0.0001
*BOUNDARY
2,1,6,0
4,1,2,0
6,1,6,0
8,1,2,0
*STEP
*STATIC
*BOUNDARY
4,3,,2.54E-05
8,3,,2.54E-05
*NODE FILE
U
*END STEP