Plotting of Concrete Reinforcement Ratio

About the development of the FEM module/workbench.

Moderator: bernd

User avatar
HarryvL
Posts: 828
Joined: Sat Jan 06, 2018 7:38 pm

Re: Plotting of Concrete Reinforcement Ratio

Postby HarryvL » Mon Jul 09, 2018 10:10 am

Interestingly, if I make Solid1 steel and Solid2 concrete I get:

Code: Select all

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: 0.000000100000 in scientific: 1.0000e-07 
<type 'list'>
concrete nodes: [5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 18L, 21L, 23L, 24L, 25L, 26L, 27L, 28L, 29L, 30L, 31L, 32L, 58L, 59L, 60L, 61L, 62L, 63L, 64L, 65L, 66L, 67L, 68L, 69L, 70L, 71L, 72L, 73L, 74L, 75L, 76L, 77L, 78L, 79L, 80L, 81L, 82L, 83L, 84L, 85L, 86L, 87L, 101L, 102L, 103L, 104L, 105L, 106L, 107L, 108L, 109L, 110L, 111L, 112L, 113L]
So it seems femmesh.meshtools.get_femnodes_by_refshape(result_mesh, ref) returns 100 for Solid1 and 113 for Solid2, whereas the total mesh only has 113 nodes.
User avatar
bernd
Posts: 7051
Joined: Sun Sep 08, 2013 8:07 pm
Location: Zürich, Switzerland

Re: Plotting of Concrete Reinforcement Ratio

Postby bernd » Mon Jul 09, 2018 1:12 pm

ok, thanks for the aditional informations. I will have a look.

Bernd
User avatar
HarryvL
Posts: 828
Joined: Sat Jan 06, 2018 7:38 pm

Re: Plotting of Concrete Reinforcement Ratio

Postby HarryvL » Tue Jul 10, 2018 3:21 am

HarryvL wrote:
Mon Jul 09, 2018 10:10 am
Interestingly, if I make Solid1 steel and Solid2 concrete I get:

Code: Select all

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: 0.000000100000 in scientific: 1.0000e-07 
<type 'list'>
concrete nodes: [5L, 6L, 7L, 8L, 9L, 10L, 11L, 12L, 18L, 21L, 23L, 24L, 25L, 26L, 27L, 28L, 29L, 30L, 31L, 32L, 58L, 59L, 60L, 61L, 62L, 63L, 64L, 65L, 66L, 67L, 68L, 69L, 70L, 71L, 72L, 73L, 74L, 75L, 76L, 77L, 78L, 79L, 80L, 81L, 82L, 83L, 84L, 85L, 86L, 87L, 101L, 102L, 103L, 104L, 105L, 106L, 107L, 108L, 109L, 110L, 111L, 112L, 113L]
So it seems femmesh.meshtools.get_femnodes_by_refshape(result_mesh, ref) returns 100 for Solid1 and 113 for Solid2, whereas the total mesh only has 113 nodes.
Urrrrrrr. Actually it is 63 nodes if you count the list. It just happens so that node number 113 is part of the list. So maybe it is right after all !! ... Yup, the list for the example with Solid1=concrete and Solid2=steel also has 63 nodes and happens to end with node 100.

Sorry for dragging you into my confusion. I think the only problem was indeed the group meshing as you indicated. So, thanks for your help.

PS: what is this group meshing for? Should it perhaps be switched off by default?
User avatar
bernd
Posts: 7051
Joined: Sun Sep 08, 2013 8:07 pm
Location: Zürich, Switzerland

Re: Plotting of Concrete Reinforcement Ratio

Postby bernd » Tue Jul 10, 2018 9:17 am

HarryvL wrote:
Tue Jul 10, 2018 3:21 am
Sorry for dragging you into my confusion. I think the only problem was indeed the group meshing as you indicated. So, thanks for your help.
Does it mean everything runs smooth without group meshing? :?:


HarryvL wrote:
Tue Jul 10, 2018 3:21 am
PS: what is this group meshing for? Should it perhaps be switched off by default?
If group meshing is activated, the following happens ...
- on meshing with GMSH for every FEM document object with reference shapes a group of nodes is created
- on writing ccx input file, if there is a group with exact the same name as an object with reference shapes the group nodes are taken and used to write the input file, means no node search in the mesh. This speeds up inp file writing enormous, but seams not to work propperly in any case.

AFAIK, it should be deactivated by default. I will hava a look.
User avatar
bernd
Posts: 7051
Joined: Sun Sep 08, 2013 8:07 pm
Location: Zürich, Switzerland

Re: Plotting of Concrete Reinforcement Ratio

Postby bernd » Tue Jul 10, 2018 11:04 am

bernd wrote:
Tue Jul 10, 2018 9:17 am
HarryvL wrote:
Tue Jul 10, 2018 3:21 am
PS: what is this group meshing for? Should it perhaps be switched off by default?
AFAIK, it should be deactivated by default. I will hava a look.
It is active by default https://github.com/FreeCAD/FreeCAD/blob ... ls.py#L268. I will deactivate it by default in my next PR.
User avatar
HarryvL
Posts: 828
Joined: Sat Jan 06, 2018 7:38 pm

Re: Plotting of Concrete Reinforcement Ratio

Postby HarryvL » Tue Jul 10, 2018 1:23 pm

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
User avatar
HarryvL
Posts: 828
Joined: Sat Jan 06, 2018 7:38 pm

Re: Plotting of Concrete Reinforcement Ratio

Postby HarryvL » Tue Jul 10, 2018 1:24 pm

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
User avatar
bernd
Posts: 7051
Joined: Sun Sep 08, 2013 8:07 pm
Location: Zürich, Switzerland

Re: Plotting of Concrete Reinforcement Ratio

Postby bernd » Tue Jul 10, 2018 9:55 pm

:)
User avatar
bernd
Posts: 7051
Joined: Sun Sep 08, 2013 8:07 pm
Location: Zürich, Switzerland

Re: Plotting of Concrete Reinforcement Ratio

Postby bernd » Thu Jul 12, 2018 4:33 pm

HarryvL wrote:
Tue Jul 10, 2018 1:24 pm
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
git commit 3416a26d7