joha2 wrote: ↑Sat Jun 02, 2018 8:14 am
The code for the "heizung" example is as follows
Hi Johannes,
Nice model. I could run it and examined the results in paraview.
Nevertheless, I can reproduce my problem with accessing the inner boundaries
Please try the following code with your example and look for the result of the integral over the face 'heizung' :
Code: Select all
print(fe.assemble(1*ff.ds['aussen'](1))) # gives correct 2610000.0 mm²; area of 'outer' window and door
print(fe.assemble(1*ff.ds['innen'](1))) # gives correct 1710000.0 mm²; area of 'inner' door
print(fe.assemble(1*ff.ds['heizung'](1))) # gives 0.0 !! ; should result 2680000,00 mm²
BR,
HoWil
BTW: usually is mixing of different references a problem when creating a mesh. You have used the boolean-fragment 'raum_heizkoerper' as main Part for the gmsh-mesh 'raummesh' but than you mixed faces of 'raum_heizkoerper' and some of its child 'heizkoerper' (for 'heizung').
Tested with:
OS: Ubuntu 18.04 LTS
Word size of OS: 64-bit
Word size of FreeCAD: 64-bit
Version: 0.18.13780 (Git)
Build type: None
Branch: master
Hash: 9fb122008b99ea2d30ed3e6f7cc93a3b7717cce5
Python version: 2.7.15rc1
Qt version: 4.8.7
Coin version: 4.0.0a
OCC version: 7.2.0
Locale: English/UnitedStates (en_US)
EDIT and on
OS: Ubuntu 18.04 LTS
Version: 0.18.13886 (Git)