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

Invalid value for *SECTION PRINT #71

Open 3rav opened 11 months ago

3rav commented 11 months ago

For CPE4 elements, the value of total surface flux (q) is half of the correct value, but for CPE4R is OK and for 3D elements.

***********************************************************************
*NODE, NSET=Nall
       1,0.000000000000e+00,0.000000000000e+00,0.000000000000e+00
       2,2.500000000000e+02,0.000000000000e+00,0.000000000000e+00
       3,2.500000000000e+02,2.500000000000e+02,0.000000000000e+00
       4,0.000000000000e+00,2.500000000000e+02,0.000000000000e+00
       5,2.500000000000e+02,5.000000000000e+02,0.000000000000e+00
       6,0.000000000000e+00,5.000000000000e+02,0.000000000000e+00
       7,2.500000000000e+02,7.500000000000e+02,0.000000000000e+00
       8,0.000000000000e+00,7.500000000000e+02,0.000000000000e+00
       9,2.500000000000e+02,1.000000000000e+03,0.000000000000e+00
      10,0.000000000000e+00,1.000000000000e+03,0.000000000000e+00
      11,5.000000000000e+02,0.000000000000e+00,0.000000000000e+00
      12,5.000000000000e+02,2.500000000000e+02,0.000000000000e+00
      13,5.000000000000e+02,5.000000000000e+02,0.000000000000e+00
      14,5.000000000000e+02,7.500000000000e+02,0.000000000000e+00
      15,5.000000000000e+02,1.000000000000e+03,0.000000000000e+00
      16,7.500000000000e+02,0.000000000000e+00,0.000000000000e+00
      17,7.500000000000e+02,2.500000000000e+02,0.000000000000e+00
      18,7.500000000000e+02,5.000000000000e+02,0.000000000000e+00
      19,7.500000000000e+02,7.500000000000e+02,0.000000000000e+00
      20,7.500000000000e+02,1.000000000000e+03,0.000000000000e+00
      21,1.000000000000e+03,0.000000000000e+00,0.000000000000e+00
      22,1.000000000000e+03,2.500000000000e+02,0.000000000000e+00
      23,1.000000000000e+03,5.000000000000e+02,0.000000000000e+00
      24,1.000000000000e+03,7.500000000000e+02,0.000000000000e+00
      25,1.000000000000e+03,1.000000000000e+03,0.000000000000e+00
*ELEMENT, TYPE=CPE4, ELSET=Eall
     1,      1,      2,      3,      4
     2,      4,      3,      5,      6
     3,      6,      5,      7,      8
     4,      8,      7,      9,     10
     5,      2,     11,     12,      3
     6,      3,     12,     13,      5
     7,      5,     13,     14,      7
     8,      7,     14,     15,      9
     9,     11,     16,     17,     12
    10,     12,     17,     18,     13
    11,     13,     18,     19,     14
    12,     14,     19,     20,     15
    13,     16,     21,     22,     17
    14,     17,     22,     23,     18
    15,     18,     23,     24,     19
    16,     19,     24,     25,     20
***
*SURFACE, NAME=Swew
1, S1
5, S1
9, S1
13, S1
***
*NSET,NSET=Nwew 
1, 
2, 
11, 
16, 
21, 

***********************************************************************
*PHYSICAL CONSTANTS,ABSOLUTE ZERO=-273.15,STEFAN BOLTZMANN=5.669E-11
***
*MATERIAL, NAME=EPS
*CONDUCTIVITY
0.040
*SPECIFIC HEAT
1.3e9
*DENSITY
12.e-10
*SOLID SECTION, ELSET=Eall, MATERIAL=EPS

***********************************************************************
*STEP,INC=1000,NLGEOM=YES
*HEAT TRANSFER,DIRECT,STEADY STATE

***********************************************************************
*FILM
4, F3, 252.150000, 8.000000e-02
8, F3, 252.150000, 8.000000e-02
12, F3, 252.150000, 8.000000e-02
16, F3, 252.150000, 8.000000e-02

*DFLUX
1, S1, 1.000000e-02
5, S1, 1.000000e-02
9, S1, 1.000000e-02
13, S1, 1.000000e-02

***********************************************************************
*NODE FILE
NT

*NODE PRINT,NSET=Nwew, TOTALS=ONLY
RFL

***********************************************************************
*EL FILE
HFL

*SECTION PRINT,SURFACE=Swew, NAME=Cieplo
FLUX

*END STEP

Load send wew abq dflux 0.01 3D output in .dat = -1.e+4 2D output in .dat (reduced element) = -1.e+1 (-1.e+1 1000=-1.e+4) 2D output in .dat (no-reduced element) = -0.5e+1 (-0.5e+1 1000=-0.5e+4)

3D brick 1000x1000x1000 2D plane 1000x1000

oldninja commented 11 months ago

I get t he same error. I think we need to look into printoutface.f, I could not find anything related to this on sectionprints.f, so it must be on printoutface.f...