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
87 stars 21 forks source link

EVOL doesn't change with deformation #73

Open FEA-eng opened 10 months ago

FEA-eng commented 10 months ago

I've noticed that while EVOL is updated when elements are removed using *MODEL CHANGE, it doesn't account for deformation. In the simple example below, EVOL stays at 1000 mm^3 while it should be 512 mm^3, according to Abaqus.

*Node
1, 0.00000000E+000, 0.00000000E+000, 1.00000000E+001
2, 0.00000000E+000, 1.00000000E+001, 1.00000000E+001
3, 1.00000000E+001, 0.00000000E+000, 1.00000000E+001
4, 1.00000000E+001, 1.00000000E+001, 1.00000000E+001
5, 0.00000000E+000, 0.00000000E+000, 0.00000000E+000
6, 0.00000000E+000, 1.00000000E+001, 0.00000000E+000
7, 1.00000000E+001, 0.00000000E+000, 0.00000000E+000
8, 1.00000000E+001, 1.00000000E+001, 0.00000000E+000
*Element, Type=C3D8, Elset=Solid_part-1
1, 1, 2, 4, 3, 5, 6, 8, 7
*Nset, Nset=Internal_Selection-1_Displacement_Rotation-1
3, 4, 7, 8
*Nset, Nset=Internal_Selection-1_Displacement_Rotation-2
1, 2, 3, 4
*Nset, Nset=Internal_Selection-1_Displacement_Rotation-3
2, 4, 6, 8
*Nset, Nset=Internal_Selection-1_Displacement_Rotation-4
5, 6, 7, 8
*Nset, Nset=Internal_Selection-1_Displacement_Rotation-5
1, 2, 5, 6
*Nset, Nset=Internal_Selection-1_Displacement_Rotation-6
1, 3, 5, 7
*Elset, Elset=Internal_Selection-1_Solid_Section-1
Solid_part-1
*Elset, Elset=Internal_Selection-1_EH_Output-1
1
*Material, Name=Material-1
*Elastic
210000, 0.3
*Solid section, Elset=Internal_Selection-1_Solid_Section-1, Material=Material-1
*Step, Nlgeom
*Static
*Boundary
Internal_Selection-1_Displacement_Rotation-1, 1, 1, -1
*Boundary
Internal_Selection-1_Displacement_Rotation-2, 3, 3, -1
*Boundary
Internal_Selection-1_Displacement_Rotation-3, 2, 2, -1
*Boundary
Internal_Selection-1_Displacement_Rotation-4, 3, 3, 1
*Boundary
Internal_Selection-1_Displacement_Rotation-5, 1, 1, 1
*Boundary
Internal_Selection-1_Displacement_Rotation-6, 2, 2, 1
*El print, Elset=Internal_Selection-1_EH_Output-1, Totals=Only
EVOL
*Node file
RF, U
*El file
S, E, NOE
*End step