... Are these the arrays you are looking for? Hmmmmm, I am not entirely sure. I just want the global node numbers for each 10-node terahedral element, element-by-element. My routine gives the following for elNodes (79 elements of 10 nodes each) [[ 40 19 10 20 168 106 196 171 66 107] [ 10 31 30 40 1...

@harry: is there still some problem which might be some bug in FreeCAD FEM I should have a look at? bernd Bernd, this post has evolved and certainly has nothing to do with what the title suggests. It is more a case study of the use of thermo-mechanical analysis with FC and CCX. I believe @boreilly ...

bernd wrote: ↑Tue Jul 10, 2018 11:04 am It is active by default https://github.com/FreeCAD/FreeCAD/blob ... ls.py#L268. I will deactivate it by default in my next PR.

Thanks. That makes more sense IMHO

bernd wrote: ↑Tue Jul 10, 2018 9:17 am Does it mean everything runs smooth without group meshing?

I think it does Bernd.

I will add a warning/error for the case mesh grouping is switched on, because that would crash my routine (node number > number of stress results)

Thanks again

Interestingly, if I make Solid1 steel and Solid2 concrete I get: object Name: Concrete object.References: [(<Part::PartFeature>, ('Solid2',))] CONCRETE ReferenceShape ... Type: Solid, Object name: BooleanFragments, Object label: BooleanFragments, Element name: Solid2 The limit if a node is in or ou...

Interestingly, if I make Solid1 steel and Solid2 concrete I get: object Name: Concrete object.References: [(<Part::PartFeature>, ('Solid2',))] CONCRETE ReferenceShape ... Type: Solid, Object name: BooleanFragments, Object label: BooleanFragments, Element name: Solid2 The limit if a node is in or out...

output: object: <App::MaterialObjectPython object> object Name: Concrete object.References: [(<Part::PartFeature>, ('Solid1',))] CONCRETE ReferenceShape ... Type: Solid, Object name: BooleanFragments, Object label: BooleanFragments, Element name: Solid1 <type 'list'> concrete nodes: [1L, 2L, 3L, 4L...

output: object: <App::MaterialObjectPython object> object Name: Concrete object.References: [(<Part::PartFeature>, ('Solid1',))] CONCRETE ReferenceShape ... Type: Solid, Object name: BooleanFragments, Object label: BooleanFragments, Element name: Solid1 <type 'list'> concrete nodes: [1L, 2L, 3L, 4L...

... how do I get ic, nsr etc ? Are these empty lists I need to create before? Thanks for doing this Bernd. The piece of code follows a few lines of existing code defining nsr and -yes- ic also gets initialized, as follows: if 'stress' in result_set: stress = result_set['stress'] nsr=len(stress) pri...

And here is the full code segment for identifying concrete nodes ... ic[node]==0 means undecided ic[node]==1 means concrete ic[node]==2 means not concrete for obj in FreeCAD.ActiveDocument.Objects: if obj.isDerivedFrom('App::MaterialObjectPython'): print("object: {}".format(obj)) print(&qu...